title: Combinatorial Howe duality of symplectic type

abstract: We give a new combinatorial interpretation of Howe dual pairs of the form $(\g,{\rm Sp}_{2\ell})$, where $\g$ is a Lie (super)algebra of classical type. This is done by establishing a symplectic analogue of the RSK algorithm associated to this pair, in a uniform way which does not depend on $\g$. We introduce an analogue of jeu de taquin sliding for spinor model of irreducible characters of a Lie superalgebra $\g$ to define $P$-tableau and show that the associated $Q$-tableau is given by a symplectic tableau due to King.


title: Phase stability of Au-Li binary systems studied using neural network potential

abstract: The miscibility of Au and Li exhibits a potential application as an adhesion layer and electrode material in secondary batteries. Here, to explore alloying properties, we constructed a neural network potential (NNP) of Au-Li binary systems based on density functional theory (DFT) calculations. To accelerate construction of NNPs, we proposed an efficient and inexpensive method of structural dataset generation. The predictions by the constructed NNP on lattice parameters and phonon properties agree well with those obtained by DFT calculations. We also investigated the mixing energy of Au$_{1-x}$Li$_{x}$ with fine composition grids, showing excellent agreement with DFT verifications. We found the existence of various compositions with structures on and slightly above the convex hull, which can explain the lack of consensus on the Au-Li stable phases in previous studies. Moreover, we newly found Au$_{0.469}$Li$_{0.531}$ as a stable phase, which has never been reported elsewhere. Finally, we examined the alloying process starting from the phase separated structure to the complete mixing phase. We found that when multiple adjacent Au atoms dissolved into Li, the alloying of the entire Au/Li interface started from the dissolved region. This paper demonstrates the applicability of NNPs toward miscible phases and provides the understanding of the alloying mechanism.


title: Experimental Analysis of Legendre Decomposition in Machine Learning

abstract: In this technical report, we analyze Legendre decomposition for non-negative tensor in theory and application. In theory, the properties of dual parameters and dually flat manifold in Legendre decomposition are reviewed, and the process of tensor projection and parameter updating is analyzed. In application, a series of verification experiments and clustering experiments with parameters on submanifold were carried out, hoping to find an effective lower dimensional representation of the input tensor. The experimental results show that the parameters on submanifold have no ability to be directly used as low-rank representations. Combined with analysis, we connect Legendre decomposition with neural networks and low-rank representation applications, and put forward some promising prospects.


title: Inter-Image Communication for Weakly Supervised Localization

abstract: Weakly supervised localization aims at finding target object regions using only image-level supervision. However, localization maps extracted from classification networks are often not accurate due to the lack of fine pixel-level supervision. In this paper, we propose to leverage pixel-level similarities across different objects for learning more accurate object locations in a complementary way. Particularly, two kinds of constraints are proposed to prompt the consistency of object features within the same categories. The first constraint is to learn the stochastic feature consistency among discriminative pixels that are randomly sampled from different images within a batch. The discriminative information embedded in one image can be leveraged to benefit its counterpart with inter-image communication. The second constraint is to learn the global consistency of object features throughout the entire dataset. We learn a feature center for each category and realize the global feature consistency by forcing the object features to approach class-specific centers. The global centers are actively updated with the training process. The two constraints can benefit each other to learn consistent pixel-level features within the same categories, and finally improve the quality of localization maps. We conduct extensive experiments on two popular benchmarks, i.e., ILSVRC and CUB-200-2011. Our method achieves the Top-1 localization error rate of 45.17% on the ILSVRC validation set, surpassing the current state-of-the-art method by a large margin. The code is available at https://github.com/xiaomengyc/I2C.


title: Magnetar giant flare origin for GRB 200415A inferred from a new scaling relation and its implications

abstract: Soft gamma-ray repeaters (SGRs) are mainly a Galactic population and originate from neutron stars with intense ($B\simeq 10^{15}{\rm \ G}$) magnetic fields ('magnetars'). Occasionally, a giant flare occurs with enormous intensity, displaying a short hard spike, followed by a weaker, oscillatory phase which exhibits the rotational period of the neutron star. If the magnetar giant flares occur in nearby galaxies, they would appear as cosmic short-hard gamma-ray bursts (GRBs) without detecting the weak oscillatory phase. Recently, a short-hard GRB named GRB 200415A was detected, with a position coincident with the Sculptor Galaxy (NGC 253), rasing the question whether it is a classic short GRB or a magentar giant flare. Here we show that magnetar giant flares follow a scaling relation between the spectral peak energy and the isotropic energy in $1\,{\rm keV} - 10\,{\rm MeV}$, i.e., $E_{\rm p}\propto E_{\rm iso}^{1/4}$, and locate in a distinct region of the $E_{\rm p}-E_{\rm iso}$ plane from that of classic short GRBs. The relation can be well understood in the model that giant flares arise from the photosphere emissions of relativistically expanding fireball. GRB 200415A, together with two other candidate giant flares (GRB 051103 and GRB 070201) follow this relation, which strongly favor the giant flare origin of these GRBs. The GeV emission detected by \emph{Fermi}/LAT from GRB 200415A at $18-285 \,$s can also be explained in the giant flare scenario. The total energy in the GeV emission implies a baryon load of $\sim 10^{23}{\rm g}$ in the giant flare fireball of GRB 200415A.


title: Determination of the strong vertices of doubly heavy baryons with pseudoscalar mesons in QCD

abstract: The strong coupling constant of doubly heavy baryons with light pseudoscalar mesons $ \pi $ and $ K $ are computed within the light cone sum rules. We take into account two-particle and three-particle distribution amplitudes of the said pseudoscalar mesons. We compare our result with the one existing in the literature.


title: Multiparametric Cardiac 18F-FDG PET in Humans: Kinetic Model Selection and Identifiability Analysis

abstract: Cardiac 18F-FDG PET has been used in clinics to assess myocardial glucose metabolism. Its ability for imaging myocardial glucose transport, however, has rarely been exploited in clinics. Using the dynamic FDG-PET scans of ten patients with coronary artery disease, we investigate in this paper appropriate dynamic scan and kinetic modeling protocols for efficient quantification of myocardial glucose transport. Three kinetic models and the effect of scan duration were evaluated by using statistical fit quality, assessing the impact on kinetic quantification, and analyzing the practical identifiability. The results show that the kinetic model selection depends on the scan duration. The reversible two-tissue model was needed for a one-hour dynamic scan. The irreversible two-tissue model was optimal for a scan duration of around 10-15 minutes. If the scan duration was shortened to 2-3 minutes, a one-tissue model was the most appropriate. For global quantification of myocardial glucose transport, we demonstrated that an early dynamic scan with a duration of 10-15 minutes and irreversible kinetic modeling was comparable to the full one-hour scan with reversible kinetic modeling. Myocardial glucose transport quantification provides an additional physiological parameter on top of the existing assessment of glucose metabolism and has the potential to enable single tracer multiparametric imaging in the myocardium.


title: The view from a Wigner bubble

abstract: In a recent no-go theorem [Bong et al, Nature Physics (2020)], we proved that the predictions of unitary quantum mechanics for an extended Wigner's friend scenario are incompatible with any theory satisfying three metaphysical assumptions, the conjunction of which we call "Local Friendliness": Absoluteness of Observed Events, Locality and No-Superdeterminism. In this paper (based on an invited talk for the QBism jubilee at the 2019 Vaxjo conference) I discuss the implications of this theorem for QBism, as seen from the perspective of experimental metaphysics. I argue that the key distinction between QBism and realist interpretations of quantum mechanics is best understood in terms of their adherence to different theories of truth: the pragmatist versus the correspondence theories. I argue that a productive pathway to resolve the measurement problem within a pragmatist view involves taking seriously the perspective of quantum betting agents, even those in what I call a "Wigner bubble". The notion of reality afforded by QBism, I propose, will correspond to the invariant elements of any theory that has pragmatic value to all rational agents -- that is, the elements that are invariant upon changes of agent perspectives. The classical notion of `event' is not among those invariants, even when those events are observed by some agent. Neither are quantum states. Nevertheless, I argue that far from solipsism, a personalist view of quantum states is an expression of its precise opposite: Copernicanism.


title: FATNN: Fast and Accurate Ternary Neural Networks

abstract: Ternary Neural Networks (TNNs) have received much attention due to being potentially orders of magnitude faster in inference, as well as more power efficient, than full-precision counterparts. However, 2 bits are required to encode the ternary representation with only 3 quantization levels leveraged. As a result, conventional TNNs have similar memory consumption and speed compared with the standard 2-bit models, but have worse representational capability. Moreover, there is still a significant gap in accuracy between TNNs and full-precision networks, hampering their deployment to real applications. To tackle these two challenges, in this work, we first show that, under some mild constraints, computational complexity of the ternary inner product can be reduced by a factor of 2. Second, to mitigate the performance gap, we elaborately design an implementation-dependent ternary quantization algorithm. The proposed framework is termed Fast and Accurate Ternary Neural Networks (FATNN). Experiments on image classification demonstrate that our FATNN surpasses the state-of-the-arts by a significant margin in accuracy. More importantly, speedup evaluation compared with various precisions is analyzed on several platforms, which serves as a strong benchmark for further research.


title: On Uniformly Sampling Traces of a Transition System (Extended Version)

abstract: A key problem in constrained random verification (CRV) concerns generation of input stimuli that result in good coverage of the system's runs in targeted corners of its behavior space. Existing CRV solutions however provide no formal guarantees on the distribution of the system's runs. In this paper, we take a first step towards solving this problem. We present an algorithm based on Algebraic Decision Diagrams for sampling bounded traces (i.e. sequences of states) of a sequential circuit with provable uniformity (or bias) guarantees, while satisfying given constraints. We have implemented our algorithm in a tool called TraceSampler. Extensive experiments show that TraceSampler outperforms alternative approaches that provide similar uniformity guarantees.


title: Sampling Based Approximate Skyline Calculation on Big Data

abstract: The existing algorithms for processing skyline queries cannot adapt to big data. This paper proposes two approximate skyline algorithms based on sampling. The first algorithm obtains a fixed size sample and computes the approximate skyline on the sample. The error of the first algorithm is relatively small in most cases, and is almost independent of the input relation size. The second algorithm returns an $(\epsilon,\delta)$-approximation for the exact skyline. The size of sample required by the second algorithm can be regarded as a constant relative to the input relation size, so is the running time. Experiments verify the error analysis of the first algorithm and show that the second algorithm is much faster than the existing skyline algorithms.


title: A matrix concentration inequality for products

abstract: We present a non-asymptotic concentration inequality for the random matrix product \begin{equation}\label{eq:Zn} Z_n = \left(I_d-\alpha X_n\right)\left(I_d-\alpha X_{n-1}\right)\cdots \left(I_d-\alpha X_1\right), \end{equation} where $\left\{X_k \right\}_{k=1}^{+\infty}$ is a sequence of bounded independent random positive semidefinite matrices with common expectation $\mathbb{E}\left[X_k\right]=\Sigma$. Under these assumptions, we show that, for small enough positive $\alpha$, $Z_n$ satisfies the concentration inequality \begin{equation}\label{eq:CTbound} \mathbb{P}\left(\left\Vert Z_n-\mathbb{E}\left[Z_n\right]\right\Vert \geq t\right) \leq 2d^2\cdot\exp\left(\frac{-t^2}{\alpha \sigma^2} \right) \quad \text{for all } t\geq 0, \end{equation} where $\sigma^2$ denotes a variance parameter.


title: Local Temperature Scaling for Probability Calibration

abstract: For semantic segmentation, label probabilities are often uncalibrated as they are typically only the by-product of a segmentation task. Intersection over Union (IoU) and Dice score are often used as criteria for segmentation success, while metrics related to label probabilities are not often explored. However, probability calibration approaches have been studied, which match probability outputs with experimentally observed errors. These approaches mainly focus on classification tasks, but not on semantic segmentation. Thus, we propose a learning-based calibration method that focuses on multi-label semantic segmentation. Specifically, we adopt a convolutional neural network to predict local temperature values for probability calibration. One advantage of our approach is that it does not change prediction accuracy, hence allowing for calibration as a post-processing step. Experiments on the COCO, CamVid, and LPBA40 datasets demonstrate improved calibration performance for a range of different metrics. We also demonstrate the good performance of our method for multi-atlas brain segmentation from magnetic resonance images.


title: Settling SETH vs. Approximate Sparse Directed Unweighted Diameter (up to (NU)NSETH)

abstract: We prove several tight results on the fine-grained complexity of approximating the diameter of a graph. First, we prove that, for any $\varepsilon>0$, assuming the Strong Exponential Time Hypothesis (SETH), there are no near-linear time $2-\varepsilon$-approximation algorithms for the Diameter of a sparse directed graph, even in unweighted graphs. This result shows that a simple near-linear time 2-approximation algorithm for Diameter is optimal under SETH, answering a question from a survey of Rubinstein and Vassilevska-Williams (SIGACT '19) for the case of directed graphs. In the same survey, Rubinstein and Vassilevska-Williams also asked if it is possible to show that there are no $2-\varepsilon$ approximation algorithms for Diameter in a directed graph in $O(n^{1.499})$ time. We show that, assuming a hypothesis called NSETH, one cannot use a deterministic SETH-based reduction to rule out the existence of such algorithms. Extending the techniques in these two results, we characterize whether a $2-\varepsilon$ approximation algorithm running in time $O(n^{1+\delta})$ for the Diameter of a sparse directed unweighted graph can be ruled out by a deterministic SETH-based reduction for every $\delta\in(0,1)$ and essentially every $\varepsilon\in(0,1)$, assuming NSETH. This settles the SETH-hardness of approximating the diameter of sparse directed unweighted graphs for deterministic reductions, up to NSETH. We make the same characterization for randomized SETH-based reductions, assuming another hypothesis called NUNSETH. We prove additional hardness and non-reducibility results for undirected graphs.


title: Mitigating the Effect of Nanoscale Porosity on Thermoelectric Power Factor of Si

abstract: The addition of porosity to thermoelectric materials can significantly increase the figure of merit, ZT, by reducing the thermal conductivity. Unfortunately, porosity is also detrimental to the thermoelectric power factor in the numerator of the figure of merit ZT. In this manuscript we derive strategies to recoup electrical performance in nanoporous Si by fine tuning the carrier concentration and through judicious design of the pore size and shape so as to provide energy selective electron filtering. In this study, we considered phosphorus doped silicon containing discrete pores that are either spheres, cylinders, cubes, or triangular prisms. The effects from these pores are compared with those from extended pores with circular, square and triangular cross sectional shape, and infinite length perpendicular to the electrical current. A semiclassical Boltzmann transport equation is used to model Si thermoelectric power factor. This model reveals three key results: The largest enhancement in Seebeck coefficient occurs with cubic pores. The fractional improvement is about 15% at low carrier concentration ($< 10^{20}\ \mathrm{1/cm^3}$) up to 60% at high carrier population with characteristic length around $\sim 1\ \mathrm{nm}$. To obtain the best energy filtering effect at room temperature, nanoporous Si needs to be doped to higher carrier concentration than is optimal for bulk Si. Finally, in $n$-type Si thermoelectrics the electron filtering effect that can be generated with nanoscale porosity is significantly lower than the ideal filtering effect; nevertheless, the enhancement in the Seebeck coefficient that can be obtained is large enough to offset the reduction in electrical conductivity caused by porosity.


title: Symmetry driven topological phases in XAgBi (X=Ba,Sr): An Ab-initio hybrid functional calculations

abstract: Density functional theory (DFT) approaches have been ubiquitously used to predict topological order and non-trivial band crossings in real materials, like Dirac, Weyl semimetals and so on. However, use of less accurate exchange-correlation functional often yields false prediction of non-trivial band order leading to misguide the experimental judgment about such materials. Using relatively more accurate hybrid functional exchange-correlation, we explore a set of (already) experimentally synthesized materials (crystallizing in space group P6_3/mmc) Our calculations based on more accurate functional helps to correct various previous predictions for this material class. Based on point group symmetry analysis and ab-initio calculations, we systematically show how lattice symmetry breaking via alloy engineering manifests different fermionic behavior, namely Dirac, triple point and Weyl in a single material. Out of various compounds, XAgBi (X=Ba,Sr) turn out to be two ideal candidates, in which the topological nodal point lie very close to the Fermi level, within minimal/no extra Fermi pocket. We further studied the surface states and Fermi arc topology on the surface of Dirac, triple point and Weyl semimetallic phases of BaAgBi. We firmly believe that, while the crystal symmetry is essential to protect the band crossings, the use of accurate exchange correlation functional in any DFT calculation is an important necessity for correct prediction of band order which can be trusted and explored in future experiments.


title: A Bayesian Approach to Spherical Factor Analysis for Binary Data

abstract: Factor models are widely used across diverse areas of application for purposes that include dimensionality reduction, covariance estimation, and feature engineering. Traditional factor models can be seen as an instance of linear embedding methods that project multivariate observations onto a lower dimensional Euclidean latent space. This paper discusses a new class of geometric embedding models for multivariate binary data in which the embedding space correspond to a spherical manifold, with potentially unknown dimension. The resulting models include traditional factor models as a special case, but provide additional flexibility. Furthermore, unlike other techniques for geometric embedding, the models are easy to interpret, and the uncertainty associated with the latent features can be properly quantified. These advantages are illustrated using both simulation studies and real data on voting records from the U.S. Senate.


title: Facial Expression Retargeting from Human to Avatar Made Easy

abstract: Facial expression retargeting from humans to virtual characters is a useful technique in computer graphics and animation. Traditional methods use markers or blendshapes to construct a mapping between the human and avatar faces. However, these approaches require a tedious 3D modeling process, and the performance relies on the modelers' experience. In this paper, we propose a brand-new solution to this cross-domain expression transfer problem via nonlinear expression embedding and expression domain translation. We first build low-dimensional latent spaces for the human and avatar facial expressions with variational autoencoder. Then we construct correspondences between the two latent spaces guided by geometric and perceptual constraints. Specifically, we design geometric correspondences to reflect geometric matching and utilize a triplet data structure to express users' perceptual preference of avatar expressions. A user-friendly method is proposed to automatically generate triplets for a system allowing users to easily and efficiently annotate the correspondences. Using both geometric and perceptual correspondences, we trained a network for expression domain translation from human to avatar. Extensive experimental results and user studies demonstrate that even nonprofessional users can apply our method to generate high-quality facial expression retargeting results with less time and effort.


title: Baryon and Lepton Number Violation from Gravitational Waves

abstract: We describe a unique gravitational wave signature for a class of models with a vast hierarchy between the symmetry breaking scales. The unusual shape of the signal is a result of the overlapping contributions to the stochastic gravitational wave background from cosmic strings produced at a high scale and a cosmological phase transition at a low scale. We apply this idea to a simple model with gauged baryon and lepton number, in which the high-scale breaking of lepton number is motivated by the seesaw mechanism for the neutrinos, whereas the low scale of baryon number breaking is required by the observed dark matter relic density. The novel signature can be searched for in upcoming gravitational wave experiments.


title: Dynamically Constrained Motion Planning Networks for Non-Holonomic Robots

abstract: Reliable real-time planning for robots is essential in today's rapidly expanding automated ecosystem. In such environments, traditional methods that plan by relaxing constraints become unreliable or slow-down for kinematically constrained robots. This paper describes the algorithm Dynamic Motion Planning Networks (Dynamic MPNet), an extension to Motion Planning Networks, for non-holonomic robots that address the challenge of real-time motion planning using a neural planning approach. We propose modifications to the training and planning networks that make it possible for real-time planning while improving the data efficiency of training and trained models' generalizability. We evaluate our model in simulation for planning tasks for a non-holonomic robot. We also demonstrate experimental results for an indoor navigation task using a Dubins car.


title: Statistical thermodynamics of the diced lattice

abstract: In this work we analyze the statistical thermodynamics of Diced lattice carriers employing a Greens function formulation to examine the grand potential, Helmholtz free energy, the grand and ordinary partition functions and entropy. This facilitates the calculation of the specific heat, and all evaluations are carried out for both the degenerate and nondegenerate statistical regimes.


title: The Stellar Mass Assembly of Low Redshift, Massive, Central Galaxies in SDSS and the TNG300 simulation

abstract: The stellar mass assembly of galaxies can be affected by both secular and environmental processes. In this study, for the first time, we investigate the stellar mass assembly of $\sim90,000$ low redshift, central galaxies selected from SDSS group catalogues (M$_{\rm Stellar}\gtrsim10^{9.5}$M$_{\odot}$, M$_{\rm Halo}\gtrsim10^{12}$M$_{\odot}$) as a function of both stellar and halo mass. We use estimates of the times at which 10, 50 and 90 per cent of the stellar mass was assembled from photometric spectral energy distribution fitting, allowing a more complete investigation than single stellar ages alone. We consider trends in both stellar and halo mass simultaneously, finding dependencies of all assembly times on both. We find that galaxies with higher stellar masses (at constant halo mass) have on average older lookback times, similar to previous studies of galaxy assembly. We also find that galaxies at higher halo mass (at constant stellar mass) have younger lookback times, possibly due to a larger reservoir of gas for star formation. An exception to this is a sub sample with high stellar-to-halo mass ratios, which are likely massive, field spirals. We compare these observed trends to those predicted by the TNG300 simulation, finding good agreement overall as a function of either stellar or halo mass. However, some differences in the assembly times (of up to $\sim 3$ Gyr) appear when considering both stellar and halo mass simultaneously, noticeably at intermediate stellar masses (M$_{\rm Stellar} \sim 10^{11}$ M$_{\odot}$). These discrepancies are possibly linked to the quenched fraction of galaxies and the kinetic mode AGN feedback implemented in TNG300.


title: The Star-Formation Properties of the Observed and Simulated AGN Universe: BAT vs EAGLE

abstract: In this paper we present data from 72 low redshift, hard X-ray selected AGN taken from the {\it Swift}-BAT 58 month catalogue. We utilise spectral energy distribution fitting to the optical to IR photometry in order to estimate host galaxy properties. We compare this observational sample to a volume and flux matched sample of AGN from the EAGLE hydrodynamical simulations in order to verify how accurately the simulations can reproduce observed AGN host galaxy properties. After correcting for the known +0.2 dex offset in the SFRs between EAGLE and previous observations, we find agreement in the SFR and X-ray luminosity distributions; however we find that the stellar masses in EAGLE are $0.2 - 0.4$ dex greater than the observational sample, which consequently leads to lower sSFRs. We compare these results to our previous study at high redshift, finding agreement in both the observations and simulations, whereby the widths of sSFR distributions are similar ($\sim0.4-0.6$ dex) and the median of the SFR distributions lie below the star forming main sequence by $\sim0.3-0.5$ dex across all samples. We also use EAGLE to select a sample of AGN host galaxies at high and low redshift and follow their characteristic evolution from $z=8$ to $z=0$. We find similar behaviour between these two samples, whereby star formation is quenched when the black hole goes through its phase of most rapid growth. Utilising EAGLE we find that 23\% of AGN selected at $z\sim0$ are also AGN at high redshift, and that their host galaxies are among the most massive objects in the simulation. Overall we find EAGLE reproduces the observations well, with some minor inconsistencies ($\sim$ 0.2 dex in stellar masses and $\sim$ 0.4 dex in sSFRs).


title: Ultra-peripheral nuclear collisions as a source of heavy quarkonia

abstract: Heavy quarkonium production in ultra-peripheral nuclear collisions (UPC) is described within the QCD dipole formalism. Realistic quarkonium wave functions in the $Q\bar Q$ rest frame are calculated by solving the Schr\"odinger equation with a subsequent Lorentz boost to high energy. We rely on several realistic $Q\bar Q$ potentials, which allow to describe well the quarkonium masses and decay widths, as well as data on diffractive electroproduction of quarkonia on protons. Nuclear effects are calculated with the phenomenological dipole cross sections fitted to DIS data. The higher twist quark shadowing related to the lowest $Q\bar Q$ Fock component of the photon, as well as the leading twist gluon shadowing, related to higher components containing gluons, are included. The results for coherent and incoherent photoproduction of charmonia and bottomonia in UPC of heavy nuclei are in good accord with available data from the LHC. They can also be verified in future experiments at electron-ion colliders.


title: A Longitudinal Method for Simultaneous Whole-Brain and Lesion Segmentation in Multiple Sclerosis

abstract: In this paper we propose a novel method for the segmentation of longitudinal brain MRI scans of patients suffering from Multiple Sclerosis. The method builds upon an existing cross-sectional method for simultaneous whole-brain and lesion segmentation, introducing subject-specific latent variables to encourage temporal consistency between longitudinal scans. It is very generally applicable, as it does not make any prior assumptions on the scanner, the MRI protocol, or the number and timing of longitudinal follow-up scans. Preliminary experiments on three longitudinal datasets indicate that the proposed method produces more reliable segmentations and detects disease effects better than the cross-sectional method it is based upon.


title: Mean-field spin-oscillation dynamics beyond the single-mode approximation for a harmonically trapped spin-1 Bose-Einstein condensate

abstract: Compared to single-component Bose-Einstein condensates, spinor Bose-Einstein condensates display much richer dynamics. In addition to density oscillations, spinor Bose-Einstein condensates exhibit intriguing spin dynamics that is associated with population transfer between different hyperfine components. This work analyzes the validity of the widely employed single-mode approximation when describing the spin dynamics in response to a quench of the system Hamiltonian. The single-mode approximation assumes that the different hyperfine states all share the same time-independent spatial mode. This implies that the resulting spin Hamiltonian only depends on the spin interaction strength and not on the density interaction strength. Taking the spinor sodium Bose-Einstein condensate in the $f=1$ hyperfine manifold as an example and working within the mean-field theory framework, it is found numerically that the single-mode approximation misses, in some parameter regimes, intricate details of the spin and spatial dynamics. We develop a physical picture that explains the observed phenomenon. Moreover, using that the population oscillations described by the single-mode approximation enter into the effective potential felt by the mean-field spinor, we derive a semi-quantitative condition for when dynamical mean-field induced corrections to the single-mode approximation are relevant. Our mean-field results have implications for a variety of published and planned experimental studies.


title: Turbulence dictates the fate of virus-containing droplets in violent expiratory events

abstract: Violent expiratory events, such as coughing and sneezing, are highly nontrivial examples of a two-phase mixture of liquid droplets dispersed into an unsteady turbulent airflow. Understanding the physical mechanisms determining the dispersion and evaporation process of respiratory droplets has recently become a priority given the global emergency caused by the SARS-CoV-2 infection. By means of high-resolution direct numerical simulations (DNS) of the expiratory airflow and a comprehensive Lagrangian model for the droplet dynamics, we identify the key role of turbulence on the fate of exhaled droplets. Due to the considerable spread in the initial droplet size, we show that the droplet evaporation time is controlled by the combined effect of turbulence and droplet inertia. This mechanism is clearly highlighted when comparing the DNS results with those obtained using coarse-grained descriptions that are employed in the majority of the current state-of-the-art investigations, resulting in errors up to $100\%$ when the turbulent fluctuations are filtered or completely averaged out.


title: The effective theory of nuclear scattering for a WIMP of arbitrary spin

abstract: We introduce a systematic approach to characterize the most general non-relativistic WIMP-nucleus interaction allowed by Galilean invariance for a WIMP of arbitrary spin $j_\chi$ in the approximation of one-nucleon currents. Five nucleon currents arise from the nonrelativistic limit of the free nucleon Dirac bilinears. Our procedure consists in (1) organizing the WIMP currents according to the rank of the $2 j_\chi+1$ irreducible operator products of up to $2 j_\chi$ WIMP spin vectors, and (2) coupling each of the WIMP currents to each of the five nucleon currents. The transferred momentum $q$ appears to a power fixed by rotational invariance. For a WIMP of spin $j_\chi$ we find a basis of 4+20$j_\chi$ independent operators that exhaust all the possible operators that drive elastic WIMP-nucleus scattering in the approximation of one-nucleon currents. By comparing our operator basis, which is complete, to the operators already introduced in the literature we show that some of the latter for $j_\chi=1$ were not independent and some were missing. We provide explicit formulas for the squared scattering amplitudes in terms of the nuclear response functions, which are available in the literature for most of the targets used in WIMP direct detection experiments.


title: Multiply positive functions, critical exponent phenomena, and the Jain-Karlin-Schoenberg kernel

abstract: We prove the converse to a result of Karlin [Trans. AMS 1964], and also strengthen his result and two results of Schoenberg [Ann. of Math. 1955]. One of the latter results concerns zeros of Laplace transforms of multiply positive functions. The other results study which powers $\alpha$ of two specific kernels are totally non-negative of order $p\geq 2$ (denoted TN$_p$); both authors showed this happens for $\alpha\geq p-2$, and Schoenberg proved that it does not for $\alpha<p-2$. We show more strongly that for every $p \times p$ submatrix of either kernel, up to a 'shift', its $\alpha$th power is totally positive of order $p$ (TP$_p$) for every $\alpha > p-2$, and is not TN$_p$ for every non-integer $\alpha\in(0,p-2)$. In particular, these results reveal 'critical exponent' phenomena in total positivity. We also prove the converse to a 1968 result of Karlin, revealing yet another critical exponent phenomenon - for Laplace transforms of all Polya Frequency (PF) functions. We further classify the powers preserving all TN$_p$ Hankel kernels on intervals, and isolate individual kernels encoding these powers. We then transfer results on preservers by Polya-Szego (1925), Loewner/Horn (1969), and Khare-Tao (in press), from positive matrices to Hankel TN$_p$ kernels. Another application constructs individual matrices encoding the Loewner convex powers. This complements Jain's results (2020) for Loewner positivity, which we strengthen to total positivity. Remarkably, these (strengthened) results of Jain, those of Schoenberg and Karlin, the latter's converse, and the above Hankel kernels all arise from a single symmetric rank-two kernel and its powers: $\max(1+xy,0)$. We also provide a novel characterization of PF functions and sequences of order $p\geq 3$, following Schoenberg's 1951 result for $p=2$. We correct a small gap in his paper, in the classification of discontinuous PF functions.


title: The Language Interpretability Tool: Extensible, Interactive Visualizations and Analysis for NLP Models

abstract: We present the Language Interpretability Tool (LIT), an open-source platform for visualization and understanding of NLP models. We focus on core questions about model behavior: Why did my model make this prediction? When does it perform poorly? What happens under a controlled change in the input? LIT integrates local explanations, aggregate analysis, and counterfactual generation into a streamlined, browser-based interface to enable rapid exploration and error analysis. We include case studies for a diverse set of workflows, including exploring counterfactuals for sentiment analysis, measuring gender bias in coreference systems, and exploring local behavior in text generation. LIT supports a wide range of models--including classification, seq2seq, and structured prediction--and is highly extensible through a declarative, framework-agnostic API. LIT is under active development, with code and full documentation available at https://github.com/pair-code/lit.


title: Invariant learning based multi-stage identification for Lithium-ion battery performance degradation

abstract: By informing accurate performance (e.g., capacity), health state management plays a significant role in safeguarding battery and its powered system. While most current approaches are primarily based on data-driven methods, lacking in-depth analysis of battery performance degradation mechanism may discount their performances. To fill in the research gap about data-driven battery performance degradation analysis, an invariant learning based method is proposed to investigate whether the battery performance degradation follows a fixed behavior. First, to unfold the hidden dynamics of cycling battery data, measurements are reconstructed in phase subspace. Next, a novel multi-stage division strategy is put forward to judge the existent of multiple degradation behaviors. Then the whole aging procedure is sequentially divided into several segments, among which cycling data with consistent degradation speed are assigned in the same stage. Simulations on a well-know benchmark verify the efficacy of the proposed multi-stages identification strategy. The proposed method not only enables insights into degradation mechanism from data perspective, but also will be helpful to related topics, such as stage of health.


title: Leveraging Automated Mixed-Low-Precision Quantization for tiny edge microcontrollers

abstract: The severe on-chip memory limitations are currently preventing the deployment of the most accurate Deep Neural Network (DNN) models on tiny MicroController Units (MCUs), even if leveraging an effective 8-bit quantization scheme. To tackle this issue, in this paper we present an automated mixed-precision quantization flow based on the HAQ framework but tailored for the memory and computational characteristics of MCU devices. Specifically, a Reinforcement Learning agent searches for the best uniform quantization levels, among 2, 4, 8 bits, of individual weight and activation tensors, under the tight constraints on RAM and FLASH embedded memory sizes. We conduct an experimental analysis on MobileNetV1, MobileNetV2 and MNasNet models for Imagenet classification. Concerning the quantization policy search, the RL agent selects quantization policies that maximize the memory utilization. Given an MCU-class memory bound of 2MB for weight-only quantization, the compressed models produced by the mixed-precision engine result as accurate as the state-of-the-art solutions quantized with a non-uniform function, which is not tailored for CPUs featuring integer-only arithmetic. This denotes the viability of uniform quantization, required for MCU deployments, for deep weights compression. When also limiting the activation memory budget to 512kB, the best MobileNetV1 model scores up to 68.4% on Imagenet thanks to the found quantization policy, resulting to be 4% more accurate than the other 8-bit networks fitting the same memory constraints.


title: Prediction on Properties of Rare-earth 2-17-X Magnets Ce2Fe17-xCoxCN : A Combined Machine-learning and Ab-initio Study

abstract: We employ a combination of machine learning and first-principles calculations to predict magnetic properties of rare-earth lean magnets. For this purpose, based on training set constructed out of experimental data, the machine is trained to make predictions on magnetic transition temperature (Tc), largeness of saturation magnetization ({\mu}0Ms), and nature of the magnetocrystalline anisotropy (Ku). Subsequently, the quantitative values of {\mu}0Ms and Ku of the yet-to-be synthesized compounds, screened by machine learning, are calculated by first-principles density functional theory. The applicability of the proposed technique of combined machine learning and first-principles calculations is demonstrated on 2-17-X magnets, Ce2Fe17-xCoxCN. Further to this study, we explore stability of the proposed compounds by calculating vacancy formation energy of small atom interstitials (N/C). Our study indicates a number of compounds in the proposed family, offers the possibility to become solution of cheap, and efficient permanent magnet.


title: Insights into the structural, electronic and optical properties of X$_2$MgZ$_4$($X=$ Sc, Y; $Z=$ S, Se) spinel compounds: Materials for the future optoelectronic applications

abstract: Direct bandgap bulk materials are very important for the optical applications. It is therefore important to predict new materials with the desired properties. In the present work, density functional theory is applied to study different physical properties of X$_2$MgZ$_4$($X=$ Sc, Y; $Z=$ S, Se) spinel compounds. Generalized gradient approximation is used to analyze the structural and elastic parameters while modified Becke Johnson exchange potential is applied to calculate electric band profiles and optical properties. All the studied compounds are stable in the cubic structure. Also the energy bandgap is of direct nature. Therefore these compounds can find useful applications in the optoelectrics devices. Optical properties of the compounds are studied in terms of dielectric function, refractive index, extinction coefficient, optical conductivity and reflectivity. The transport parameters like electrical conductivity, Seebeck coefficient, and thermal conductivity are also evaluated.


title: On Higher order Poincar\'e Inequalities with radial derivatives and Hardy improvements on the hyperbolic space

abstract: In this paper we prove higher order Poincar\'e inequalities involving radial derivatives namely, \begin{equation*} \int_{\mathbb{H}^{N}} |\nabla_{r,\mathbb{H}^{N}}^{k} u|^2 \, {\rm d}v_{\mathbb{H}^{N}} \geq \bigg(\frac{N-1}{2}\bigg)^{2(k-l)} \int_{\mathbb{H}^{N}} |\nabla_{r,\mathbb{H}^{N}}^{l} u|^2 \, {\rm d}v_{\mathbb{H}^{N}} \ \ \text{ for all } u\in H^k(\mathbb{H}^{N}), \end{equation*} where underlying space is $N$-dimensional hyperbolic space $\mathbb{H}^{N}$, $0\leq l<k$ are integers and the constant $\big(\frac{N-1}{2}\big)^{2(k-l)}$ is sharp. Furthermore we improve the above inequalities by adding Hardy-type remainder terms and the sharpness of some constants is also discussed.


title: Comments on Mathematical Modeling of Current Source Matrix Converter with Venturini and SVM

abstract: In this paper, authors want to comment on a recently published article describing the Mathematical Modeling of Current Source Matrix Converter (CSMC) with two modulation strategies, namely: Venturini and Space Vector Modulation (SVM). Reported flaws are broadly classified into two (2) categories; namely: (1) Technical Flaws and (2) Rubric and Grammar Flaws. This paper not only reports these flaws but also provides suggestive rectifications for correction and improvement of the published article.


title: Open Set Recognition with Conditional Probabilistic Generative Models

abstract: Deep neural networks have made breakthroughs in a wide range of visual understanding tasks. A typical challenge that hinders their real-world applications is that unknown samples may be fed into the system during the testing phase, but traditional deep neural networks will wrongly recognize these unknown samples as one of the known classes. Open set recognition (OSR) is a potential solution to overcome this problem, where the open set classifier should have the flexibility to reject unknown samples and meanwhile maintain high classification accuracy in known classes. Probabilistic generative models, such as Variational Autoencoders (VAE) and Adversarial Autoencoders (AAE), are popular methods to detect unknowns, but they cannot provide discriminative representations for known classification. In this paper, we propose a novel framework, called Conditional Probabilistic Generative Models (CPGM), for open set recognition. The core insight of our work is to add discriminative information into the probabilistic generative models, such that the proposed models can not only detect unknown samples but also classify known classes by forcing different latent features to approximate conditional Gaussian distributions. We discuss many model variants and provide comprehensive experiments to study their characteristics. Experiment results on multiple benchmark datasets reveal that the proposed method significantly outperforms the baselines and achieves new state-of-the-art performance.


title: Power and spectral characterization of photonic integrated circuit based axicon like lens

abstract: We demonstrate an on-chip Silicon-on-Insulator (SOI) axicon etched using a low resolution (200 nm feature size, 250 nm gap) deep-ultraviolet lithographic fabrication. The axicon consists of circular gratings with seven stages of 1x2 multimode interferometers. We present a technique to apodize the gratings azimuthally by breaking up the circles into arcs which successfully increased the penetration depth in the gratings from $\approx$5 $\mu$m to $\approx$55 $\mu$m. We characterize the device's performance by coupling 1300$\pm$50 nm swept source laser in to the chip from the axicon, and measuring the out-coupled light from a grating coupler. Further, we also present the implementation of balanced homodyne detection method for the spectral characterization of the device and show that the position of the output lobe of the axicon does not change significantly with wavelength.


title: Learning to Reason in Round-based Games: Multi-task Sequence Generation for Purchasing Decision Making in First-person Shooters

abstract: Sequential reasoning is a complex human ability, with extensive previous research focusing on gaming AI in a single continuous game, round-based decision makings extending to a sequence of games remain less explored. Counter-Strike: Global Offensive (CS:GO), as a round-based game with abundant expert demonstrations, provides an excellent environment for multi-player round-based sequential reasoning. In this work, we propose a Sequence Reasoner with Round Attribute Encoder and Multi-Task Decoder to interpret the strategies behind the round-based purchasing decisions. We adopt few-shot learning to sample multiple rounds in a match, and modified model agnostic meta-learning algorithm Reptile for the meta-learning loop. We formulate each round as a multi-task sequence generation problem. Our state representations combine action encoder, team encoder, player features, round attribute encoder, and economy encoders to help our agent learn to reason under this specific multi-player round-based scenario. A complete ablation study and comparison with the greedy approach certify the effectiveness of our model. Our research will open doors for interpretable AI for understanding episodic and long-term purchasing strategies beyond the gaming community.


title: Object Detection for Graphical User Interface: Old Fashioned or Deep Learning or a Combination?

abstract: Detecting Graphical User Interface (GUI) elements in GUI images is a domain-specific object detection task. It supports many software engineering tasks, such as GUI animation and testing, GUI search and code generation. Existing studies for GUI element detection directly borrow the mature methods from computer vision (CV) domain, including old fashioned ones that rely on traditional image processing features (e.g., canny edge, contours), and deep learning models that learn to detect from large-scale GUI data. Unfortunately, these CV methods are not originally designed with the awareness of the unique characteristics of GUIs and GUI elements and the high localization accuracy of the GUI element detection task. We conduct the first large-scale empirical study of seven representative GUI element detection methods on over 50k GUI images to understand the capabilities, limitations and effective designs of these methods. This study not only sheds the light on the technical challenges to be addressed but also informs the design of new GUI element detection methods. We accordingly design a new GUI-specific old-fashioned method for non-text GUI element detection which adopts a novel top-down coarse-to-fine strategy, and incorporate it with the mature deep learning model for GUI text detection.Our evaluation on 25,000 GUI images shows that our method significantly advances the start-of-the-art performance in GUI element detection.


title: An Inter- and Intra-Band Loss for Pansharpening Convolutional Neural Networks

abstract: Pansharpening aims to fuse panchromatic and multispectral images from the satellite to generate images with both high spatial and spectral resolution. With the successful applications of deep learning in the computer vision field, a lot of scholars have proposed many convolutional neural networks (CNNs) to solve the pansharpening task. These pansharpening networks focused on various distinctive structures of CNNs, and most of them are trained by L2 loss between fused images and simulated desired multispectral images. However, L2 loss is designed to directly minimize the difference of spectral information of each band, which does not consider the inter-band relations in the training process. In this letter, we propose a novel inter- and intra-band (IIB) loss to overcome the drawback of original L2 loss. Our proposed IIB loss can effectively preserve both inter- and intra-band relations and can be directly applied to different pansharpening CNNs.


title: Schatten class positive Toeplitz operators on Bergman spaces of the Siegel upper half-space

abstract: We characterize Schatten class membership of positive Toeplitz operators defined on the Bergman spaces over the Siegel upper half-space in terms of averaging functions and Berezin transforms in the range of $0<p<\infty$.


title: Fully $S$-coidempotent modules

abstract: Let R be be a commutative ring with identity and S be a multiplicatively closed subset of R. In this article we introduce the concepts of S-coidempotent submodules and fully S-coidempotent R-modules as generalizations of coidempotent submodules and fully coidempotent R-modules. We explore some basic properties of these classes of R-modules.


title: Quantization dimension and stability for infinite self-similar measures with respect to geometric mean error

abstract: Let $\mu$ be a Borel probability measure associated with an iterated function system consisting of a countably infinite number of contracting similarities and an infinite probability vector. In this paper, we study the quantization dimension of the measure $\mu$ with respect to the geometric mean error. The quantization for infinite systems is different from the well-known finite case investigated by Graf and Luschgy. That is, many tools which are used in the finite setting, for example, existence of finite maximal antichains, fail in the infinite case. We prove that the quantization dimension of the measure $\mu$ is equal to its Hausdorff dimension which extends a well-known result of Graf and Luschgy for the finite case to an infinite setting. In the last section, we discuss the stability of quantization dimension for infinite systems.


title: Molecular dynamics simulation of crack growth in mono-crystal nickel with voids and inclusions

abstract: In this study, the crack propagation of the pre-cracked mono-crystal nickel with the voids and inclusions has been investigated by molecular dynamics simulations. Different sizes of voids, inclusions and materials of inclusions are used to fully study the effect of the voids and inclusions during the crack propagation process. The dislocations evolution, stress distribution and crack length are analyzed as the associated mechanical properties. The results indicate that the voids and inclusions can change the path of crack propagation of the pre-cracked mono-crystal nickel. Moreover, the results show that the voids and inclusions can lead a better resistance to plastic deformation of the mono-crystal and the inclusions can make the system more difficult to fracture.


title: Quantum coherence, quantum Fisher information and teleportation in the Ising-Heisenberg spin chain model of a heterotrimetallic Fe-Mn-Cu coordination polymer with magnetic impurity

abstract: The effect of non-uniform magnetic fields on the Ising-Heisenberg chain of a heterotrimetallic coordination compound $\mathrm{Fe-Mn-Cu}$, modeling a magnetic impurity on one dimer is studied. This impurity is configured by imposing non-uniform magnetic fields on each sites of $j-$th interstitial ionic dimer $\mathrm{Mn}^{2+}-\mathrm{Cu}^{2+}$ of the chain model. The quantum coherence and pairwise entanglement between spin-1/2 magnetic $\mathrm{Mn}^{2+}-\mathrm{Cu}^{2+}$ ion dimers clearly depend on the site which the impurity is located. It is demonstrated that when the magnetic impurity is considered for one magnetic dimer, by altering Ising nodal exchange interaction and Heisenberg anisotropy parameter of the impurity dimer, the entanglement can be enhanced to the maximal value 1 at a special fixed magnetic field. Moreover, we find that the quantum Fisher information of the model with magnetic impurity behaves considerably different from the original model. Besides of the quantum resources like concurrence, we prove that the quantum Fisher information can be used as new quantum tool for estimating the quantum phase transition in the model under consideration. On the other hand, thermal teleportation can be significantly optimized by adjusting the magnetic impurity, and a strong increase in the average fidelity is observed. Finally, the magnetic impurity can be manipulated to locally control the thermal entanglement, coherence, quantum Fisher information and teleportation unlike the original model where it is done globally.


title: Neumann-Rosochatius system for (m,n) string in $AdS_3 \times S^3$ with mixed flux

abstract: $SL(2,\mathbb{Z})$ invariant action for probe $(m,n)$ string in $AdS_3\times S^3\times T^4$ with mixed three-form fluxes can be described by an integrable deformation of an one-dimensional Neumann-Rosochatius (NR) system. We present the deformed features of the integrable model and study general class of rotating and pulsating solutions by solving the integrable equations of motion. For the rotating string, the explicit solutions can be expressed in terms of elliptic functions. We make use of the integrals of motion to find out the scaling relation among conserved charges for the particular case of constant radii solutions. Then we study the closed $(m,n)$ string pulsating in $R_t\times S^3$. We find the string profile and calculate the total energy of such pulsating string in terms of oscillation number $(\cal{N})$ and angular momentum $(\cal{J})$.


title: The Italian bondage and reinforcement numbers of digraphs

abstract: An \textit{Italian dominating function} on a digraph $D$ with vertex set $V(D)$ is defined as a function $f : V(D) \rightarrow \{0, 1, 2\}$ such that every vertex $v \in V(D)$ with $f(v) = 0$ has at least two in-neighbors assigned $1$ under $f$ or one in-neighbor $w$ with $f(w) = 2$. The \textit{weight} of an Italian dominating function $f$ is the value $\omega(f) = f(V(D)) = \sum_{u \in V(D)} f(u)$. The \textit{Italian domination number} of a digraph $D$, denoted by $\gamma_I(D)$, is the minimum taken over the weights of all Italian dominating functions on $D$. The \textit{Italian bondage number} of a digraph $D$, denoted by $b_I(D)$, is the minimum number of arcs of $A(D)$ whose removal in $D$ results in a digraph $D'$ with $\gamma_I(D') > \gamma_I(D)$. The \textit{Italian reinforcement number} of a digraph $D$, denoted by $r_I(D)$, is the minimum number of extra arcs whose addition to $D$ results in a digraph $D'$ with $\gamma_I(D') < \gamma_I(D)$. In this paper, we initiate the study of Italian bondage and reinforcement numbers in digraphs and present some bounds for $b_I(D)$ and $r_I(D)$. We also determine the Italian bondage and reinforcement numbers of some classes of digraphs.


title: Coded Elastic Computing on Machines with Heterogeneous Storage and Computation Speed

abstract: We study the optimal design of heterogeneous Coded Elastic Computing (CEC) where machines have varying computation speeds and storage. CEC introduced by Yang et al. in 2018 is a framework that mitigates the impact of elastic events, where machines can join and leave at arbitrary times. In CEC, data is distributed among machines using a Maximum Distance Separable (MDS) code such that subsets of machines can perform the desired computations. However, state-of-the-art CEC designs only operate on homogeneous networks where machines have the same speeds and storage. This may not be practical. In this work, based on an MDS storage assignment, we develop a novel computation assignment approach for heterogeneous CEC networks to minimize the overall computation time. We first consider the scenario where machines have heterogeneous computing speeds but same storage and then the scenario where both heterogeneities are present. We propose a novel combinatorial optimization formulation and solve it exactly by decomposing it into a convex optimization problem for finding the optimal computation load and a "filling problem" for finding the exact computation assignment. A low-complexity "filling algorithm" is adapted and can be completed within a number of iterations equals at most the number of available machines.


title: Debris cloud of India Anti-Satellite Test to Microsat-R Satellite

abstract: Understanding the motion of debris cloud produced by the anti-satellite test can help us to know the danger of these tests. This study presents the orbit status of 57 fragments observed by the CelesTrak and presented in the NORAD Two-Line Element Sets of India Anti-Satellite Test. There are 10 of these observed fragments have altitudes of the apogee larger than 1000.0km, the maximum one is 1725.7km. We also numerical calculated the number of debris, the results show that the number of debris with the diameter larger than 0.2m is 14, the number of debris with the diameter larger than 0.01m is 6587, and the number of debris with the diameter larger than 0.001m is 7.22e+5. The results of the secondary collision of the debris will produced more fragments in the space. The life time of the fragments depends on the initial orbit parameters and the sizes of the debris.


title: Detectability of ultra-compact X-ray binaries as LISA sources

abstract: Ultra-compact X-ray binaries (UCXBs) are low-mass X-ray binaries with ultra-short orbital periods (usually less than 1 hour) and hydrogen-poor donor stars, which are proposed to be the potential LISA sources. In this work, we firstly employ the MESA code to examine the parameter space of the progenitors of UCXBs that LISA will detect. Our simulations indicate that the initial binaries with a neutron star and a $0.4-3.3~M_{\odot}$ companion star in an orbit of initial orbital period smaller than the bifurcation period could evolve into UCXBs, some of which will emit gravitational wave signals that can be detectable by the LISA. However, the initial orbital periods of the binaries that will evolve into UCXB-LISA sources in a distance of 10 kpc are located in a very narrow range, i.e. the formation of these LISA source requires an extremely fine-tuning of initial parameter. According to the characteristic stains and the derived maximum detectable distances, four sources among eight UCXBs with the observed distances are expected to be detected by the LISA. Based on the parameter space given by the detailed binary evolution models and the rapid binary star evolution code, the birthrate of UCXBs appearing as LISA sources in the Galaxy is estimated to be $(2-2.6)\times10^{\rm -6}~\rm yr^{-1}$. Considering the contribution of UCXBs in the globular clusters, the number of UCXB-LISA sources can reach $240 - 320$. Although the formation condition is severe, the detectability of UCXBs by the LISA is still optimistic and significant because they provide an opportunity to pursue full multi-messenger investigations.


title: Topological edge states of a graphene zigzag nanoribbon with spontaneous edge magnetism

abstract: The topological phases of graphene with spin-orbit coupling, an exchange field, and a staggered-sublattice potential determine the properties of the edge states of the zigzag nanoribbon. In the presence of the Hubbard interaction, the spontaneous magnetization at the zigzag terminations induces sizable magnetic moments at the lattice sites in the bulk region. Thus, the exchange field and staggered-sublattice potential in the bulk region are effectively changed, which in turn change the topological phase. Within a certain parameter regime, quasi-stable excited states of the zigzag nanoribbon exist, which have a different magnetism configuration at the zigzag terminations from the ground state. The quasi-stable excited states could effectively suppress the finite size effect of the topological edge states. The investigation of the topological edge states in the presence of interaction helps the engineering of spintronic nanodevices based on realistic materials.


title: Further Unifying the Landscape of Cell Probe Lower Bounds

abstract: In a landmark paper, P\v{a}tra\c{s}cu demonstrated how a single lower bound for the static data structure problem of reachability in the butterfly graph, could be used to derive a wealth of new and previous lower bounds via reductions. These lower bounds are tight for numerous static data structure problems. Moreover, he also showed that reachability in the butterfly graph reduces to dynamic marked ancestor, a classic problem used to prove lower bounds for dynamic data structures. Unfortunately, P\v{a}tra\c{s}cu's reduction to marked ancestor loses a $\lg \lg n$ factor and therefore falls short of fully recovering all the previous dynamic data structure lower bounds that follow from marked ancestor. In this paper, we revisit P\v{a}tra\c{s}cu's work and give a new lossless reduction to dynamic marked ancestor, thereby establishing reachability in the butterfly graph as a single seed problem from which a range of tight static and dynamic data structure lower bounds follow.


title: 1000 Days to First Light: Construction of the Perth-Lowell Telescope Facility 1968-71

abstract: Negotiations began 1n 1968 for a telescope facility at Perth Observatory for NASA's International Planetary Patrol Network. 1,000 days later the telescope saw first light. The facility bears no resemblance to other observatories. Inside a dome, the telescope sits on a 42 ft tall concrete pier with a wrap-around staircase and concrete legs. Surrounding forest is similar in height to the dome, the design of which is counter intuitive. This study investigated why, at the risk of compromising performance, there was a departure from standard design, and to to identify drivers for the decision making. Observatory visitors learn of a government architect, Tadeusz Andrzejaczek who made whimsical, successive increases to the height of the structure. Though designed in collaboration with Acting Government Astronomer, Bertrand Harris, it is improbable that a public servant architect would have such influence over a scientific installation. Vibration amelioration was met by designing massive strength and rigidity into the structure. Thermal expansion and wind stresses were reduced using features such as shade fins and protective walls, and ground thermal disturbance was addressed by simply making it tall. Seeing measurements were not a significant design consideration. The facility exists with its current floor height because of successive approvals for modification. The initial design was by Harris and requests for redesigns came from him but in close negotiation the Andrzejaczek who desired a structure of futuristic shape and proportions. Harris's designs were influenced by his personal English background and the Old Perth Observatory where he worked as an astronomer. Andrzejaczek's design was influenced by an observatory in his birth city, his alignment with contemporary designers and his artistic flair.


title: Tail risk forecasting using Bayesian realized EGARCH models

abstract: This paper develops a Bayesian framework for the realized exponential generalized autoregressive conditional heteroskedasticity (realized EGARCH) model, which can incorporate multiple realized volatility measures for the modelling of a return series. The realized EGARCH model is extended by adopting a standardized Student-t and a standardized skewed Student-t distribution for the return equation. Different types of realized measures, such as sub-sampled realized variance, sub-sampled realized range, and realized kernel, are considered in the paper. The Bayesian Markov chain Monte Carlo (MCMC) estimation employs the robust adaptive Metropolis algorithm (RAM) in the burn in period and the standard random walk Metropolis in the sample period. The Bayesian estimators show more favourable results than maximum likelihood estimators in a simulation study. We test the proposed models with several indices to forecast one-step-ahead Value at Risk (VaR) and Expected Shortfall (ES) over a period of 1000 days. Rigorous tail risk forecast evaluations show that the realized EGARCH models employing the standardized skewed Student-t distribution and incorporating sub-sampled realized range are favored, compared to a range of models.


title: Unidirectional Oriented Water Wire in Short Nanotube

abstract: The orientation of water molecules is the key factor for the fast transport of water in small nanotubes. It has been accepted that the bidirectional water burst in short nanotubes can be transformed into unidirectional transport when the orientation of water molecules is maintained in long nanotubes under the external field. In this work, based on molecular dynamics simulations and first-principles calculations, we showed without external field, it only needs 21 water molecules to maintain the unidirectional single file water intrinsically in carbon nanotube at seconds. Detailed analysis indicates that the surprising result comes from the step by step process for the flip of water chain, which is different with the perceived concerted mechanism. Considering the thickness of cell membrane (normally 5-10 nm) is larger than the length threshold of the unidirectional water wire, this study suggests it may not need the external field to maintain the unidirectional flow in the water channel at the macroscopic timescale.

categories: physics.comp-ph physics.atm-clus physics.chem-ph quant-ph


title: ASAP-Net: Attention and Structure Aware Point Cloud Sequence Segmentation

abstract: Recent works of point clouds show that mulit-frame spatio-temporal modeling outperforms single-frame versions by utilizing cross-frame information. In this paper, we further improve spatio-temporal point cloud feature learning with a flexible module called ASAP considering both attention and structure information across frames, which we find as two important factors for successful segmentation in dynamic point clouds. Firstly, our ASAP module contains a novel attentive temporal embedding layer to fuse the relatively informative local features across frames in a recurrent fashion. Secondly, an efficient spatio-temporal correlation method is proposed to exploit more local structure for embedding, meanwhile enforcing temporal consistency and reducing computation complexity. Finally, we show the generalization ability of the proposed ASAP module with different backbone networks for point cloud sequence segmentation. Our ASAP-Net (backbone plus ASAP module) outperforms baselines and previous methods on both Synthia and SemanticKITTI datasets (+3.4 to +15.2 mIoU points with different backbones). Code is availabe at https://github.com/intrepidChw/ASAP-Net


title: Identification of MISO systems in Minimal Realization Form

abstract: The paper is concerned with identifying transfer functions of individual input channels in minimal realization form of a Multi-Input Single Output (MISO) from the input-output data corrupted by the error in all the variables. Such a framework is commonly referred to as error-in-variables (EIV). A common approach in the existing methods for identification of MISO systems is to estimate a non-minimal order transfer function under a subset of simplistic assumptions like homoskedastic error variances, known order, and delay. In this work, we deal with the challenging problem of identifying order, delay in each input of minimal realization form separately while estimating the transfer functions. We also estimate the heteroskedastic noise variances in each of the multiple inputs and output variables. An automated approach for the identification of MISO systems of minimal realization form in the EIV framework is proposed. Numerical case studies are presented to illustrate the efficacy of the proposed algorithm in identifying the transfer function along with the order, delay, and noise variances.


title: Polaritonic cylinders as multifunctional metamaterials: Single scattering and effective medium description

abstract: Polaritonic materials, owing to a strong phonon-polariton resonance in the THz and far-infrared parts of the electromagnetic spectrum, offer both high-index dielectric and metallic response in this regime. This complex response makes them suitable candidates for the design of metamaterial-related phenomena and applications. Here we show that one type of polaritonic-material-based structures that are particularly suitable for the achievement of a wide range of metamaterial properties are systems of polaritonic rods. To study the interplay between the material and the structural resonances in such systems we employ as model systems rods of LiF and SiC and we calculate first the scattering properties of a single rod, identifying and discussing the behavior of the different resonances for different rod diameters. To analyze the response of ensembles of polaritonic rods we employ an effective medium approach based on the Coherent Potential Approximation (CPA), which is shown to be superior to the simple Maxwell-Garnett approximation for polaritonic and high-index dielectric metamaterials. Calculating and analyzing the CPA effective parameters, we found that our systems exhibit a large variety of interesting metamaterial properties, including hyperbolic dispersion, epsilon-near-zero and negative refractive index response. This rich response, achievable in almost any system of polaritonic rods, is highly engineerable by properly selecting the radius and the filling ratio of the rods, making polaritonic rod systems an ideal platform for demonstration of multifunctional metamaterials.


title: Stability of EHI and regularity of MMD spaces

abstract: We show that elliptic Harnack inequality is stable under form-bounded perturbations for strongly local Dirichlet forms on complete locally compact separable metric spaces that satisfy metric doubling property (or equivalently, relative ball connectedness property).


title: Wishart and random density matrices: Analytical results for the mean-square Hilbert-Schmidt distance

abstract: Hilbert-Schmidt distance is one of the prominent distance measures in quantum information theory which finds applications in diverse problems, such as construction of entanglement witnesses, quantum algorithms in machine learning, and quantum state tomography. In this work, we calculate exact and compact results for the mean square Hilbert-Schmidt distance between a random density matrix and a fixed density matrix, and also between two random density matrices. In the course of derivation, we also obtain corresponding exact results for the distance between a Wishart matrix and a fixed Hermitian matrix, and two Wishart matrices. We verify all our analytical results using Monte Carlo simulations. Finally, we apply our results to investigate the Hilbert-Schmidt distance between reduced density matrices generated using coupled kicked tops.

categories: quant-ph math-ph math.MP physics.data-an stat.AP


title: "Conserved charges" of the Bondi-Metzner-Sachs algebra in the Brans-Dicke theory

abstract: The asymptotic symmetries in Brans-Dicke theory are analyzed using Penrose's conformal completion method, which is independent of the coordinate system used. These symmetries indeed include supertranslations and Lorentz transformations for an asymptotically flat spacetime. With the Wald-Zoupas formalism, "conserved charges" and fluxes of the Bondi-Metzner-Sachs algebra are computed. The scalar degree of freedom contributes only to the Lorentz boost charge, although it plays a role in various fluxes. The flux-balance laws are further used to constrain displacement memory, spin memory and center-of-mass memory effects.


title: Elliptic Quantum Curves of Class $\mathcal{S}_k$

abstract: Quantum curves arise from Seiberg-Witten curves associated to 4d $\mathcal{N}=2$ gauge theories by promoting coordinates to non-commutative operators. In this way the algebraic equation of the curve is interpreted as an operator equation where a Hamiltonian acts on a wave-function with zero eigenvalue. We find that this structure generalises when one considers torus-compactified 6d $\mathcal{N}=(1,0)$ SCFTs. The corresponding quantum curves are elliptic in nature and hence the associated eigenvectors/eigenvalues can be expressed in terms of Jacobi forms. In this paper we focus on the class of 6d SCFTs arising from M5 branes transverse to a $\mathbb{C}^2/\mathbb{Z}_k$ singularity. In the limit where the compactified 2-torus has zero size, the corresponding 4d $\mathcal{N}=2$ theories are known as class $\mathcal{S}_k$. We explicitly show that the eigenvectors associated to the quantum curve are expectation values of codimension 2 surface operators, while the corresponding eigenvalues are codimension 4 Wilson surface expectation values.


title: HOSE-Net: Higher Order Structure Embedded Network for Scene Graph Generation

abstract: Scene graph generation aims to produce structured representations for images, which requires to understand the relations between objects. Due to the continuous nature of deep neural networks, the prediction of scene graphs is divided into object detection and relation classification. However, the independent relation classes cannot separate the visual features well. Although some methods organize the visual features into graph structures and use message passing to learn contextual information, they still suffer from drastic intra-class variations and unbalanced data distributions. One important factor is that they learn an unstructured output space that ignores the inherent structures of scene graphs. Accordingly, in this paper, we propose a Higher Order Structure Embedded Network (HOSE-Net) to mitigate this issue. First, we propose a novel structure-aware embedding-to-classifier(SEC) module to incorporate both local and global structural information of relationships into the output space. Specifically, a set of context embeddings are learned via local graph based message passing and then mapped to a global structure based classification space. Second, since learning too many context-specific classification subspaces can suffer from data sparsity issues, we propose a hierarchical semantic aggregation(HSA) module to reduces the number of subspaces by introducing higher order structural information. HSA is also a fast and flexible tool to automatically search a semantic object hierarchy based on relational knowledge graphs. Extensive experiments show that the proposed HOSE-Net achieves the state-of-the-art performance on two popular benchmarks of Visual Genome and VRD.


title: Towards Geometry Guided Neural Relighting with Flash Photography

abstract: Previous image based relighting methods require capturing multiple images to acquire high frequency lighting effect under different lighting conditions, which needs nontrivial effort and may be unrealistic in certain practical use scenarios. While such approaches rely entirely on cleverly sampling the color images under different lighting conditions, little has been done to utilize geometric information that crucially influences the high-frequency features in the images, such as glossy highlight and cast shadow. We therefore propose a framework for image relighting from a single flash photograph with its corresponding depth map using deep learning. By incorporating the depth map, our approach is able to extrapolate realistic high-frequency effects under novel lighting via geometry guided image decomposition from the flashlight image, and predict the cast shadow map from the shadow-encoding transformed depth map. Moreover, the single-image based setup greatly simplifies the data capture process. We experimentally validate the advantage of our geometry guided approach over state-of-the-art image-based approaches in intrinsic image decomposition and image relighting, and also demonstrate our performance on real mobile phone photo examples.


title: Balanced Depth Completion between Dense Depth Inference and Sparse Range Measurements via KISS-GP

abstract: Estimating a dense and accurate depth map is the key requirement for autonomous driving and robotics. Recent advances in deep learning have allowed depth estimation in full resolution from a single image. Despite this impressive result, many deep-learning-based monocular depth estimation (MDE) algorithms have failed to keep their accuracy yielding a meter-level estimation error. In many robotics applications, accurate but sparse measurements are readily available from Light Detection and Ranging (LiDAR). Although they are highly accurate, the sparsity limits full resolution depth map reconstruction. Targeting the problem of dense and accurate depth map recovery, this paper introduces the fusion of these two modalities as a depth completion (DC) problem by dividing the role of depth inference and depth regression. Utilizing the state-of-the-art MDE and our Gaussian process (GP) based depth-regression method, we propose a general solution that can flexibly work with various MDE modules by enhancing its depth with sparse range measurements. To overcome the major limitation of GP, we adopt Kernel Interpolation for Scalable Structured (KISS)-GP and mitigate the computational complexity from O(N^3) to O(N). Our experiments demonstrate that the accuracy and robustness of our method outperform state-of-the-art unsupervised methods for sparse and biased measurements.


title: Prevalence, Contents and Automatic Detection of KL-SATD

abstract: When developers use different keywords such as TODO and FIXME in source code comments to describe self-admitted technical debt (SATD), we refer it as Keyword-Labeled SATD (KL-SATD). We study KL-SATD from 33 software repositories with 13,588 KL-SATD comments. We find that the median percentage of KL-SATD comments among all comments is only 1,52%. We find that KL-SATD comment contents include words expressing code changes and uncertainty, such as remove, fix, maybe and probably. This makes them different compared to other comments. KL-SATD comment contents are similar to manually labeled SATD comments of prior work. Our machine learning classifier using logistic Lasso regression has good performance in detecting KL-SATD comments (AUC-ROC 0.88). Finally, we demonstrate that using machine learning we can identify comments that are currently missing but which should have a SATD keyword in them. Automating SATD identification of comments that lack SATD keywords can save time and effort by replacing manual identification of comments. Using KL-SATD offers a potential to bootstrap a complete SATD detector.


title: Comment on "Solar Wind and Heavy Ion Properties of Interplanetary Coronal Mass Ejections" by M.J.Owens

abstract: A recent paper by Owens (2018) presents a statistical analysis of the properties of interplanetary coronal mass ejections (ICMEs) and variations in their compositions and ion charge-state signatures using data from the Advanced Composition Explorer (ACE) spacecraft. However, the article contains several serious flaws, which we will discuss here.


title: Overcoming the numerical sign problem in the Wigner dynamics via adaptive particle annihilation

abstract: The infamous numerical sign problem poses a fundamental obstacle to particle-based stochastic Wigner simulations in high dimensional phase space. Although the existing particle annihilation via uniform mesh significantly alleviates the sign problem when dimensionality D $\le$ 4, the mesh size grows dramatically when D $\ge$ 6 due to the curse of dimensionality and consequently makes the annihilation very inefficient. In this paper, we propose an adaptive particle annihilation algorithm, termed Sequential-clustering Particle Annihilation via Discrepancy Estimation (SPADE), to overcome the sign problem. SPADE follows a divide-and-conquer strategy: Adaptive clustering of particles via controlling their number-theoretic discrepancies and independent random matching in each group, and it may learn the minimal amount of particles that can accurately capture the non-classicality of the Wigner function. Combining SPADE with the variance reduction technique based on the stationary phase approximation, we attempt to simulate the proton-electron couplings in 6-D and 12-D phase space. A thorough performance benchmark of SPADE is provided with the reference solutions in 6-D phase space produced by a characteristic-spectral-mixed scheme under a $73^3 \times 80^3$ uniform grid, which fully explores the limit of grid-based deterministic Wigner solvers.


title: On some invariants of cubic fourfolds

abstract: For a general cubic fourfold $X \subset \mathbb{P}^5$, we compute the Hodge numbers of the locus $S \subset F$ of lines of second type. We also give an upper bound of 6 for the degree of irrationality of the Fano scheme of lines of any smooth cubic hypersurface.


title: Implications on Feature Detection when using the Benefit-Cost Ratio

abstract: In many practical machine learning applications, there are two objectives: one is to maximize predictive accuracy and the other is to minimize costs of the resulting model. These costs of individual features may be financial costs, but can also refer to other aspects, like for example evaluation time. Feature selection addresses both objectives, as it reduces the number of features and can improve the generalization ability of the model. If costs differ between features, the feature selection needs to trade-off the individual benefit and cost of each feature. A popular trade-off choice is the ratio of both, the BCR (benefit-cost ratio). In this paper we analyze implications of using this measure with special focus to the ability to distinguish relevant features from noise. We perform a simulation study for different cost and data settings and obtain detection rates of relevant features and empirical distributions of the trade-off ratio. Our simulation study exposed a clear impact of the cost setting on the detection rate. In situations with large cost differences and small effect sizes, the BCR missed relevant features and preferred cheap noise features. We conclude that a trade-off between predictive performance and costs without a controlling hyperparameter can easily overemphasize very cheap noise features. While the simple benefit-cost ratio offers an easy solution to incorporate costs, it is important to be aware of its risks. Avoiding costs close to 0, rescaling large cost differences, or using a hyperparameter trade-off are ways to counteract the adverse effects exposed in this paper.


title: The network footprint of replication in popular DBMSs

abstract: Database replication is an important component of reliable, disaster tolerant and highly available distributed systems. However, data replication also causes communication and processing overhead. Quantification of these overheads is crucial in choosing a suitable DBMS form several available options and capacity planning. In this paper, we present results from a comparative empirical analysis of replication activities of three commonly used DBMSs - MySQL, PostgreSQL and Cassandra under text as well as image traffic. In our experiments, the total traffic with two replicas (which is the norm) was as much as $300$\% higher than the total traffic with no replica. Furthermore, activation of the compression option for replication traffic, built into MySQL, reduced the total network traffic by as much as $20$\%. We also found that average CPU utilization and memory utilization were not impacted by the number of replicas or the dataset.


title: On the Galois Theory of Generalized Laguerre Polynomials and Trimmed Exponential

abstract: Inspired by the work of Schur on the Taylor series of the exponential and Laguerre polynomials, we study the Galois theory of trimmed exponentials $f_{n,n+k}=\sum_{i=0}^{k} \frac{x^{i}}{(n+i)!}$ and of the generalized Laguerre polynomials $L^{(n)}_k$ of degree $k$. We show that if $n$ is chosen uniformly from $\{1,\ldots, x\}$, then, asymptotically almost surely, for all $k\leq x^{o(1)}$ the Galois groups of $f_{n,n+k}$ and of $L_{k}^{(n)}$ are the full symmetric group $S_k$.


title: The Mathematical Foundations of Physical Systems Modeling Languages

abstract: Modern modeling languages for general physical systems, such as Modelica, Amesim, or Simscape, rely on Differential Algebraic Equations (DAE), i.e., constraints of the form f(dot{x},x,u)=0. This drastically facilitates modeling from first principles of the physics and the reuse of models. In this paper we develop the mathematical theory needed to establish the development of compilers and tools for DAE based physical modeling languages on solid mathematical bases. Unlike Ordinary Differential Equations, DAE exhibit subtle issues because of the notion of differentiation index and related latent equations -- ODE are DAE of index zero for which no latent equation needs to be considered. Prior to generating execution code and calling solvers, the compilation of such languages requires a nontrivial \emph{structural analysis} step that reduces the differentiation index to a level acceptable by DAE solvers. The models supported by tools of the Modelica class involve multiple modes with mode-dependent DAE based dynamics and state-dependent mode switching. Multimode DAE are much more difficult than DAE. The main difficulty is the handling of the events of mode change. Unfortunately, the large literature devoted to the mathematical analysis of DAEs does not cover the multimode case, typically saying nothing about mode changes. This lack of foundations causes numerous difficulties to the existing modeling tools. Some models are well handled, others are not, with no clear boundary between the two classes. In this paper, we develop a comprehensive mathematical approach supporting compilation and code generation for this class of languages. Its core is the structural analysis of multimode DAE systems. As a byproduct of this structural analysis, we propose well sound criteria for accepting or rejecting models. For our mathematical development, we rely on nonstandard analysis, which allows us to cast hybrid systems dynamics to discrete time dynamics with infinitesimal step size, thus providing a uniform framework for handling both continuous dynamics and mode change events.


title: Noncritical Weighted Hardy's Inequalities with compact perturbations

abstract: Let $\Omega$ be a bounded domain of $\mathbb{R}^N$ whose boundary is a $\mathbb{C}^2$ compact manifolds. In the present paper we shall study a variational problem relating the weighted Hardy inequalities with sharp missing terms. As weights we adopted powers of the distance function $\delta(x)$ to the boundary $\partial\Omega$.


title: Caching Placement and Resource Allocation for Cache-Enabling UAV NOMA Networks

abstract: This article investigates the cache-enabling unmanned aerial vehicle (UAV) cellular networks with massive access capability supported by non-orthogonal multiple access (NOMA). The delivery of a large volume of multimedia contents for ground users is assisted by a mobile UAV base station, which caches some popular contents for wireless backhaul link traffic offloading. In cache-enabling UAV NOMA networks, the caching placement of content caching phase and radio resource allocation of content delivery phase are crucial for network performance. To cope with the dynamic UAV locations and content requests in practical scenarios, we formulate the long-term caching placement and resource allocation optimization problem for content delivery delay minimization as a Markov decision process (MDP). The UAV acts as an agent to take actions for caching placement and resource allocation, which includes the user scheduling of content requests and the power allocation of NOMA users. In order to tackle the MDP, we propose a Q-learning based caching placement and resource allocation algorithm, where the UAV learns and selects action with \emph{soft ${\varepsilon}$-greedy} strategy to search for the optimal match between actions and states. Since the action-state table size of Q-learning grows with the number of states in the dynamic networks, we propose a function approximation based algorithm with combination of stochastic gradient descent and deep neural networks, which is suitable for large-scale networks. Finally, the numerical results show that the proposed algorithms provide considerable performance compared to benchmark algorithms, and obtain a trade-off between network performance and calculation complexity.


title: High Rate Hybrid MnO2@CNT Fabric Anode for Li-Ion Batteries: Properties and Lithium Storage Mechanism by In-Situ Synchrotron X-Ray Scattering

abstract: High-performance anodes for rechargeable Li-ion battery are produced by nanostructuring of the transition metal oxides on a conductive support. Here, we demonstrate a hybrid material of MnO2 directly grown onto fabrics of carbon nanotube fibres, which exhibits notable specific capacity over 1100 and 500 mAh/g at a discharge current density of 25 mA/g and 5 A/g, respectively, with coulombic efficiency of 97.5 %. Combined with 97 % capacity retention after 1500 cycles at a current density of 5 A/g, both capacity and stability are significantly above literature data. Detailed investigations involving electrochemical and in situ synchrotron X-ray scattering study reveal that during galvanostatic cycling, MnO2 undergoes an irreversible phase transition to LiMnO2, which stores lithium through an intercalation process, followed by conversion mechanism and pseudocapacitive processes. This mechanism is further confirmed by Raman spectroscopy and X-ray photoelectron spectroscopy. The fraction of pseudocapacitive charge storage ranges from 27% to 83%, for current densities from 25 mA/g to 5 A/g. Firm attachment of the active material to the built-in current collector makes the electrodes flexible and mechanically robust, and ensures that the low charge transfer resistance and the high electrode surface area remain after irreversible phase transition of the active material and extensive cycling.


title: MAYONNAISE: a morphological components analysis pipeline for circumstellar disks and exoplanets imaging in the near infrared

abstract: Imaging circumstellar disks in the near-infrared provides unprecedented information about the formation and evolution of planetary systems. However, current post-processing techniques for high-contrast imaging using ground-based telescopes have a limited sensitivity to extended signals and their morphology is often plagued with strong morphological distortions. Moreover, it is challenging to disentangle planetary signals from the disk when the two components are close or intertwined. We propose a pipeline that is capable of detecting a wide variety of disks and preserving their shapes and flux distributions. By construction, our approach separates planets from disks. After analyzing the distortions induced by the current angular differential imaging (ADI) post-processing techniques, we establish a direct model of the different components constituting a temporal sequence of high-contrast images. In an inverse problem framework, we jointly estimate the starlight residuals and the potential extended sources and point sources hidden in the images, using low-complexity priors for each signal. To verify and estimate the performance of our approach, we tested it on VLT/SPHERE-IRDIS data, in which we injected synthetic disks and planets. We also applied our approach on observations containing real disks. Our technique makes it possible to detect disks from ADI datasets of a contrast above $3\times10^{-6}$ with respect to the host star. As no specific shape of the disks is assumed, we are capable of extracting a wide diversity of disks, including face-on disks. The intensity distribution of the detected disk is accurately preserved and point sources are distinguished, even close to the disk.


title: Fast and Eager k-Medoids Clustering: O(k) Runtime Improvement of the PAM, CLARA, and CLARANS Algorithms

abstract: Clustering non-Euclidean data is difficult, and one of the most used algorithms besides hierarchical clustering is the popular algorithm Partitioning Around Medoids (PAM), also simply referred to as k-medoids clustering. In Euclidean geometry the mean-as used in k-means-is a good estimator for the cluster center, but this does not exist for arbitrary dissimilarities. PAM uses the medoid instead, the object with the smallest dissimilarity to all others in the cluster. This notion of centrality can be used with any (dis-)similarity, and thus is of high relevance to many domains and applications. A key issue with PAM is its high run time cost. We propose modifications to the PAM algorithm that achieve an O(k)-fold speedup in the second ("SWAP") phase of the algorithm, but will still find the same results as the original PAM algorithm. If we relax the choice of swaps performed (while retaining comparable quality), we can further accelerate the algorithm by eagerly performing additional swaps in each iteration. With the substantially faster SWAP, we can now explore faster initialization strategies, because (i) the classic ("BUILD") initialization now becomes the bottleneck, and (ii) our swap is fast enough to compensate for worse starting conditions. We also show how the CLARA and CLARANS algorithms benefit from the proposed modifications. While we do not study the parallelization of our approach in this work, it can easily be combined with earlier approaches to use PAM and CLARA on big data (some of which use PAM as a subroutine, hence can immediately benefit from these improvements), where the performance with high k becomes increasingly important. In experiments on real data with k=100,200, we observed a 458x respectively 1191x speedup compared to the original PAM SWAP algorithm, making PAM applicable to larger data sets, and in particular to higher k.


title: PyMGRIT: A Python Package for the parallel-in-time method MGRIT

abstract: In this paper, we introduce the Python framework PyMGRIT, which implements the multigrid-reduction-in-time (MGRIT) algorithm for solving the (non-)linear systems arising from the discretization of time-dependent problems. The MGRIT algorithm is a reduction-based iterative method that allows parallel-in-time simulations, i. e., calculating multiple time steps simultaneously in a simulation, by using a time-grid hierarchy. The PyMGRIT framework features many different variants of the MGRIT algorithm, ranging from different multigrid cycle types and relaxation schemes, as well as various coarsening strategies, including time-only and space-time coarsening, to using different time integrators on different levels in the multigrid hierachy. PyMGRIT allows serial runs for prototyping and testing of new approaches, as well as parallel runs using the Message Passing Interface (MPI). Here, we describe the implementation of the MGRIT algorithm in PyMGRIT and present the usage from both user and developer point of views. Three examples illustrate different aspects of the package, including pure time parallelism as well as space-time parallelism by coupling PyMGRIT with PETSc or Firedrake, which enable spatial parallelism through MPI.