title: Homological mirror symmetry for invertible polynomials in two variables

abstract: In this paper, we give a proof of homological mirror symmetry for two variable invertible polynomials, where the symmetry group on the $B$-side is taken to be maximal. The proof involves an explicit gluing construction of the Milnor fibres, and as an application, we prove derived equivalences between certain nodal stacky curves, some of whose irreducible components have non-trivial generic stabiliser.


title: Reconfigurable Parallel Architecture of High Speed Round Robin Arbiter

abstract: With a view to managing the increasing traffic in computer networks, round robin arbiter has been proposed to work with packet switching system to have increased speed in providing access and scheduling. Round robin arbiter is a doorway to a particular bus based on request along with equal priority and gives turns to devices connected to it in a cyclic order. Considering the rapid growth in computer networking and the emergence of computer automation which will need much more access to the existing limited resources, this paper emphasizes on designing a reconfigurable round robin arbiter over FPGA which takes parallel requests and processes them with high efficiency and less delay than existing designs. Proposed round robin arbiter encounters with 4 to 12 devices. Results show that with 200% increment in the number of connected devices, only 2.69% increment has been found in the delay. With less delay, proposed round robin arbiter exhibits high speed performance with higher traffic, which is a new feature in comparison with the existing designs.


title: The dominance of big teams in China's scientific output

abstract: Modern science is dominated by scientific productions from teams. A recent finding shows that teams with both large and small sizes are essential in research, prompting us to analyze the extent to which a country's scientific work is carried out by big/small teams. Here, using over 26 million publications from Web of Science, we find that China's research output is more dominated by big teams than the rest of the world, which is particularly the case in fields of natural science. Despite the global trend that more papers are done by big teams, China's drop in small team output is much steeper. As teams in China shift from small to large size, the team diversity that is essential for innovative works does not increase as much as that in other countries. Using the national average as the baseline, we find that the National Natural Science Foundation of China (NSFC) supports fewer small team works than the National Science Foundation of U.S. (NSF) does, implying that big teams are more preferred by grant agencies in China. Our finding provides new insights into the concern of originality and innovation in China, which urges a need to balance small and big teams.


title: Learning from Suspected Target: Bootstrapping Performance for Breast Cancer Detection in Mammography

abstract: Deep learning object detection algorithm has been widely used in medical image analysis. Currently all the object detection tasks are based on the data annotated with object classes and their bounding boxes. On the other hand, medical images such as mammography usually contain normal regions or objects that are similar to the lesion region, and may be misclassified in the testing stage if they are not taken care of. In this paper, we address such problem by introducing a novel top likelihood loss together with a new sampling procedure to select and train the suspected target regions, as well as proposing a similarity loss to further identify suspected targets from targets. Mean average precision (mAP) according to the predicted targets and specificity, sensitivity, accuracy, AUC values according to classification of patients are adopted for performance comparisons. We firstly test our proposed method on a private dense mammogram dataset. Results show that our proposed method greatly reduce the false positive rate and the specificity is increased by 0.25 on detecting mass type cancer. It is worth mention that dense breast typically has a higher risk for developing breast cancers and also are harder for cancer detection in diagnosis, and our method outperforms a reported result from performance of radiologists. Our method is also validated on the public Digital Database for Screening Mammography (DDSM) dataset, brings significant improvement on mass type cancer detection and outperforms the most state-of-the-art work.


title: Adaptive Millimeter-Wave Communications Exploiting Mobility and Blockage Dynamics

abstract: Mobility may degrade the performance of next-generation vehicular networks operating at the millimeter-wave spectrum: frequent loss of alignment and blockages require repeated beam training and handover, thus incurring huge overhead. In this paper, an adaptive and joint design of beam training, data transmission and handover is proposed, that exploits the mobility process of mobile users and the dynamics of blockages to optimally trade-off throughput and power consumption. At each time slot, the serving base station decides to perform either beam training, data communication, or handover when blockage is detected. The problem is cast as a partially observable Markov decision process, and solved via an approximate dynamic programming algorithm based on PERSEUS [2]. Numerical results show that the PERSEUS-based policy performs near-optimally, and achieves a 55% gain in spectral efficiency compared to a baseline scheme with periodic beam training. Inspired by its structure, an adaptive heuristic policy is proposed with low computational complexity and small performance degradation.


title: Unsupervised Domain Adaptation for Mammogram Image Classification: A Promising Tool for Model Generalization

abstract: Generalization is one of the key challenges in the clinical validation and application of deep learning models to medical images. Studies have shown that such models trained on publicly available datasets often do not work well on real-world clinical data due to the differences in patient population and image device configurations. Also, manually annotating clinical images is expensive. In this work, we propose an unsupervised domain adaptation (UDA) method using Cycle-GAN to improve the generalization ability of the model without using any additional manual annotations.


title: Combinatorial Nullstellensatz and DP-coloring of Graphs

abstract: We initiate the study of applying the Combinatorial Nullstellensatz to the DP-coloring of graphs even though, as is well-known, the Alon-Tarsi theorem does not apply to DP-coloring. We define the notion of good covers of prime order which allows us to apply the Combinatorial Nullstellensatz to DP-coloring. We apply these tools to DP-coloring of the cones of certain bipartite graphs and uniquely 3-colorable graphs. We also extend a result of Akbari, Mirrokni, and Sadjad (2006) on unique list colorability to the context of DP-coloring. We establish a sufficient algebraic condition for a graph $G$ to satisfy $\chi_{DP}(G) \leq 3$, and we completely determine the DP-chromatic number of squares of all cycles.


title: Warwick Electron Microscopy Datasets

abstract: Large, carefully partitioned datasets are essential to train neural networks and standardize performance benchmarks. As a result, we have set up new repositories to make our electron microscopy datasets available to the wider community. There are three main datasets containing 19769 scanning transmission electron micrographs, 17266 transmission electron micrographs, and 98340 simulated exit wavefunctions, and multiple variants of each dataset for different applications. To visualize image datasets, we trained variational autoencoders to encode data as 64-dimensional multivariate normal distributions, which we cluster in two dimensions by t-distributed stochastic neighbor embedding. In addition, we have improved dataset visualization with variational autoencoders by introducing encoding normalization and regularization, adding an image gradient loss, and extending t-distributed stochastic neighbor embedding to account for encoded standard deviations. Our datasets, source code, pretrained models, and interactive visualizations are openly available at https://github.com/Jeffrey-Ede/datasets.


title: A POD-Galerkin reduced order model of a turbulent convective buoyant flow of sodium over a backward-facing step

abstract: A Finite-Volume based POD-Galerkin reduced order modeling strategy for steady-state Reynolds averaged Navier--Stokes (RANS) simulation is extended for low-Prandtl number flow. The reduced order model is based on a full order model for which the effects of buoyancy on the flow and heat transfer are characterized by varying the Richardson number. The Reynolds stresses are computed with a linear eddy viscosity model. A single gradient diffusion hypothesis, together with a local correlation for the evaluation of the turbulent Prandtl number, is used to model the turbulent heat fluxes. The contribution of the eddy viscosity and turbulent thermal diffusivity fields are considered in the reduced order model with an interpolation based data-driven method. The reduced order model is tested for buoyancy-aided turbulent liquid sodium flow over a vertical backward-facing step with a uniform heat flux applied on the wall downstream of the step. The wall heat flux is incorporated with a Neumann boundary condition in both the full order model and the reduced order model. The velocity and temperature profiles predicted with the reduced order model for the same and new Richardson numbers inside the range of parameter values are in good agreement with the RANS simulations. Also, the local Stanton number and skin friction distribution at the heated wall are qualitatively well captured. Finally, the reduced order simulations, performed on a single core, are about $10^5$ times faster than the RANS simulations that are performed on eight cores.


title: A Framework for Interdomain and Multioutput Gaussian Processes

abstract: One obstacle to the use of Gaussian processes (GPs) in large-scale problems, and as a component in deep learning system, is the need for bespoke derivations and implementations for small variations in the model or inference. In order to improve the utility of GPs we need a modular system that allows rapid implementation and testing, as seen in the neural network community. We present a mathematical and software framework for scalable approximate inference in GPs, which combines interdomain approximations and multiple outputs. Our framework, implemented in GPflow, provides a unified interface for many existing multioutput models, as well as more recent convolutional structures. This simplifies the creation of deep models with GPs, and we hope that this work will encourage more interest in this approach.


title: International Large Detector: Interim Design Report

abstract: The ILD detector is proposed for an electron-positron collider with collision centre-of-mass energies from 90~\GeV~to about 1~\TeV. It has been developed over the last 10 years by an international team of scientists with the goal to design and eventually propose a fully integrated detector, primarily for the International Linear Collider, ILC. In this report the fundamental ideas and concepts behind the ILD detector are discussed and the technologies needed for the realisation of the detector are reviewed. The document starts with a short review of the science goals of the ILC, and how the goals can be achieved today with the detector technologies at hand. After a discussion of the ILC and the environment in which the experiment will take place, the detector is described in more detail, including the status of the development of the technologies foreseen for each subdetector. The integration of the different sub-systems into an integrated detector is discussed, as is the interface between the detector and the collider. This is followed by a concise summary of the benchmarking which has been performed in order to find an optimal balance between performance and cost. To the end the costing methodology used by ILD is presented, and an updated cost estimate for the detector is presented. The report closes with a summary of the current status and of planned future actions.


title: Inferring the location of reflecting surfaces exploiting loudspeaker directivity

abstract: Accurate sound field reproduction in rooms is often limited by the lack of knowledge of the room characteristics. Information about the room shape or nearby reflecting boundaries can, in principle, be used to improve the accuracy of the reproduction. In this paper, we propose a method to infer the location of nearby reflecting boundaries from measurements on a microphone array. As opposed to traditional methods, we explicitly exploit the loudspeaker directivity model (beyond omnidirectional radiation) and the microphone array geometry. This approach does not require noiseless timing information of the echoes as input, nor a tailored loudspeaker-wall-microphone measurement step. Simulations show the proposed model outperforms current methods that disregard directivity in reverberant environments.


title: Spectroscopic Signatures of Electron-Phonon Coupling in Silicon-Vacancy Centers in Diamond

abstract: Vacancy centers in diamond have proven to be a viable solid-state platform for quantum coherent opto-electronic applications. Among the variety of vacancy centers, silicon-vacancy (SiV) centers have recently attracted much attention as an inversion-symmetric system that is less susceptible to electron-phonon interactions. Nevertheless, phonon-mediated processes still degrade the coherent properties of SiV centers, however characterizing their electron-phonon coupling is extremely challenging due to their weak spectroscopic signatures and remains an open experimental problem. In this paper we theoretically investigate signatures of electron-phonon coupling in simulated linear and nonlinear spectra of SiV centers. We demonstrate how even extremely weak electron-phonon interactions, such as in SiV centers, may be completely characterized via nonlinear spectroscopic techniques and even resolved between different fine-structure transitions.


title: Kraken reveals itself -- the merger history of the Milky Way reconstructed with the E-MOSAICS simulations

abstract: Globular clusters (GCs) formed when the Milky Way experienced a phase of rapid assembly. We use the wealth of information contained in the Galactic GC population to quantify the properties of the satellite galaxies from which the Milky Way assembled. To achieve this, we train an artificial neural network on the E-MOSAICS cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies. The network uses the ages, metallicities, and orbital properties of GCs that formed in the same progenitor galaxies to predict the stellar masses and accretion redshifts of these progenitors. We apply the network to Galactic GCs associated with five progenitors: {\it Gaia}-Enceladus, the Helmi streams, Sequoia, Sagittarius, and the recently discovered, `low-energy' GCs, which provide an excellent match to the predicted properties of the enigmatic galaxy `Kraken'. The five galaxies cover a narrow stellar mass range [$M_\star=(0.6{-}4.6)\times10^8~{\rm M}_\odot$], but have widely different accretion redshifts ($z_{\rm acc}=0.57{-}2.65$). All accretion events represent minor mergers, but Kraken likely represents the most major merger ever experienced by the Milky Way, with stellar and virial mass ratios of $r_{M_\star}=1$:$31^{+34}_{-16}$ and $r_{M_{\rm h}}=1$:$7^{+4}_{-2}$, respectively. The progenitors match the $z=0$ relation between GC number and halo virial mass, but have elevated specific frequencies, suggesting an evolution with redshift. Even though these progenitors likely were the Milky Way's most massive accretion events, they contributed a total mass of only $\log{(M_{\rm \star,tot}/{\rm M}_\odot)}=9.0\pm0.1$, similar to the stellar halo. This implies that the Milky Way grew its stellar mass mostly by in-situ star formation. We conclude by organising these accretion events into the most detailed reconstruction to date of the Milky Way's merger tree.


title: The expansion of stripped-envelope stars: consequences for supernovae and gravitational-wave progenitors

abstract: Massive binaries that merge as compact objects are the progenitors of gravitational-wave sources. Most of these binaries experience one or more phases of mass transfer, during which one of the stars loses part or all of its outer envelope and becomes a stripped-envelope star. The evolution of the size of these stripped stars is crucial in determining whether they experience further interactions and their final fate. We present new calculations of stripped-envelope stars based on binary evolution models computed with MESA. We use these to investigate their radius evolution as a function of mass and metallicity. We further discuss their pre-supernova observable characteristics and potential consequences of their evolution on the properties of supernovae from stripped stars. At high metallicity we find that practically all of the hydrogen-rich envelope is removed, in agreement with earlier findings. Only progenitors with initial masses below 10\Msun expand to large radii (up to 100\Rsun), while more massive progenitors stay compact. At low metallicity, a substantial amount of hydrogen remains and the progenitors can, in principle, expand to giant sizes (> 400\Rsun), for all masses we consider. This implies that they can fill their Roche lobe anew. We show that the prescriptions commonly used in population synthesis models underestimate the stellar radii by up to two orders of magnitude. We expect that this has consequences for the predictions for gravitational-wave sources from double neutron star mergers, in particular for their metallicity dependence.


title: Octagon at finite coupling

abstract: We study a special class of four-point correlation functions of infinitely heavy half-BPS operators in planar N=4 SYM which admit factorization into a product of two octagon form factors. We demonstrate that these functions satisfy a system of nonlinear integro-differential equations which are powerful enough to fully determine their dependence on the 't Hooft coupling and two cross ratios. At weak coupling, solution to these equations yields a known series representation of the octagon in terms of ladder integrals. At strong coupling, we develop a systematic expansion of the octagon in the inverse powers of the coupling constant and calculate accompanying expansion coefficients analytically. We examine the strong coupling expansion of the correlation function in various kinematical regions and observe a perfect agreement both with the expected asymptotic behavior dictated by the OPE and with results of numerical evaluation. We find that, surprisingly enough, the strong coupling expansion is Borel summable. Applying the Borel-Pade summation method, we show that the strong coupling expansion correctly describes the correlation function over a wide region of the 't Hooft coupling.


title: A Semi-Automated Computational Approach for Infrared Dark Cloud Localization: A Catalog of Infrared Dark Clouds

abstract: The field of computer vision has greatly matured in the past decade, and many of the methods and techniques can be useful for astronomical applications. One example is in searching large imaging surveys for objects of interest, especially when it is difficult to specify the characteristics of the objects being searched for. We have developed a method using contour finding and convolution neural networks (CNNs) to search for Infrared Dark Clouds (IRDCs) in the Spitzer Galactic plane survey data. IRDCs can vary in size, shape, orientation, and optical depth, and are often located near regions with complex emission from molecular clouds and star formation, which can make the IRDCs difficult to reliably identify. False positives can occur in regions where emission is absent, rather than from a foreground IRDC. The contour finding algorithm we implemented found most closed figures in the mosaic and we developed rules to filter out some of the false positive before allowing the CNNs to analyze them. The method was applied to the Spitzer data in the Galactic plane surveys, and we have constructed a catalog of IRDCs which includes additional parts of the Galactic plane that were not included in earlier surveys.


title: Pre-Common-Envelope Mass Loss from Coalescing Binary Systems

abstract: Binary systems undergoing unstable Roche Lobe overflow spill gas into their circumbinary environment as their orbits decay toward coalescence. In this paper, we use a suite of hydrodynamic models of coalescing binaries involving an extended donor and a more compact accretor. We focus on the period of unstable Roche Lobe overflow that ends as the accretor plunges within the envelope of the donor at the onset of a common envelope phase. During this stage, mass is removed from the donor and flung into the circumbinary environment. Across a wide range of binary mass ratios, we find that the mass expelled as the separation decreases from the Roche limit to the donor's original radius is of the order of 25% of the accretor's mass. We study the kinematics of this ejecta and its dependencies on binary properties and find that it assembles into a toroidal circumbinary distribution. These circumbinary tori have approximately constant specific angular momentum due to momentum transport by spiral shocks launched from the orbiting binary. We show that an analytic model with these torus properties captures many of the main features of the azimuthally-averaged profiles of our hydrodynamic simulations. Our results, in particular the simple relationship between accretor mass and expelled mass and its spatial distribution, may be useful in interpreting stellar coalescence transients like luminous red novae, and in initializing hydrodynamic simulations of the subsequent common envelope phase.


title: A new method for deriving composition of S-type asteroids from noisy and incomplete near-infrared spectra

abstract: The surface composition of S-type asteroids can be determined using band parameters extracted from their near-infrared (NIR) spectra (0.7-2.50 $\mu$m) along with spectral calibrations derived from laboratory samples. In the past, these empirical equations have been obtained by combining NIR spectra of meteorite samples with information about their composition and mineral abundance. For these equations to give accurate results, the characteristics of the laboratory spectra they are derived from should be similar to those of asteroid spectral data (i.e., similar signal-to-noise ratio (S/N) and wavelength range). Here we present new spectral calibrations that can be used to determine the mineral composition of ordinary chondrite-like S-type asteroids. Contrary to previous work, the S/N of the ordinary chondrite spectra used in this study has been decreased to recreate the S/N typically observed among asteroid spectra, allowing us to obtain more realistic results. In addition, the new equations have been derived for five wavelength ranges encompassed between 0.7 and 2.50 $\mu$m, making it possible to determine the composition of asteroids with incomplete data. The new spectral calibrations were tested using band parameters measured from the NIR spectrum of asteroid (25143) Itokawa, and comparing the results with laboratory measurements of the returned samples. We found that the spectrally derived olivine and pyroxene chemistry, which are given by the molar contents of fayalite (Fa) and ferrosilite (Fs), are in excellent agreement with the mean values measured from the samples (Fa$_{28.6\pm1.1}$ and Fs$_{23.1\pm2.2}$), with a maximum difference of 0.6 mol\% for Fa and 1.4 mol\% for Fs.


title: Towards a general parametrization of the warm dark matter halo mass function

abstract: Studies of flux anomalies statistics and perturbations in stellar streams have the potential to constrain models of warm dark matter (WDM), including sterile neutrinos. Producing these constraints requires a parametrization of the WDM mass function relative to that of the cold dark matter (CDM) equivalent. We use five WDM models with half-mode masses, $M_\mathrm{hm}=[1.3,35]\times10^{8}$~$M_{\odot}$, spread across simulations of the Local Group, lensing ellipticals and the $z=2$ universe, to generate such a parametrization: we fit parameters to a functional form for the WDM-to-CDM halo mass function ratio, $n_\mathrm{WDM}(M_{X})/n_\mathrm{CDM}(M_{X})$, of ($1+(\alpha M_\mathrm{hm}/M_{X})^{\beta})^{\gamma}$. For $M_{X}\equiv$ virial mass of central halos we obtain $\alpha=2.3$, $\beta=0.8$, and $\gamma=-1.0$, and this fit is steeper than the extended Press-Schechter formalism predicts. For $M_{X}\equiv$ mass of subhalos we instead obtain $\alpha=4.2$, $\beta=2.5$ and $\gamma=-0.2$; in both mass definitions the scatter is $\sim20$~per~cent. The second fit typically underestimates the relative abundance of $z=2$ WDM subhaloes at the tens of per cent level. We caution that robust constraints will require bespoke simulations and a careful definition of halo mass, particularly for subhalos of mass $<10^{8}M_{\odot}$.


title: A free-floating or wide-orbit planet in the microlensing event OGLE-2019-BLG-0551

abstract: High-cadence observations of the Galactic bulge by the microlensing surveys led to the discovery of a handful of extremely short-timescale microlensing events that can be attributed to free-floating or wide-orbit planets. Here, we report the discovery of another strong free-floating planet candidate, which was found from the analysis of the gravitational microlensing event OGLE-2019-BLG-0551. The light curve of the event is characterized by a very short duration (<3 d) and a very small amplitude (< 0.1 mag). From modeling of the light curve, we find that the Einstein timescale, tE = 0.381 +/- 0.017 d, is much shorter, and the angular Einstein radius, thetaE = 4.35 +/- 0.34 uas, is much smaller than those of typical lensing events produced by stellar-mass lenses (tE ~ 20 d, thetaE ~ 0.3 mas), indicating that the lens is very likely to be a planetary-mass object. We conduct an extensive search for possible signatures of a companion star in the light curve of the event, finding no significant evidence for the putative host star. For the first time, we also demonstrate that the angular Einstein radius of the lens does not depend on blending in the low-magnification events with strong finite source effects.

categories: astro-ph.EP astro-ph.GA astro-ph.IM astro-ph.SR


title: Generating metal-polluting debris in white dwarf planetary systems from small-impact crater ejecta

abstract: Metal pollution in white dwarf photospheres originates from the accretion of some combination of planets, moons, asteroids, comets, boulders, pebbles and dust. When large bodies reside in dynamically stagnant locations -- unable themselves to pollute nor even closely approach the white dwarf -- then smaller reservoirs of impact debris may become a complementary or the primary source of metal pollutants. Here, we take a first step towards exploring this possibility by computing limits on the recoil mass that escapes the gravitational pull of the target object following a single impact onto an atmosphere-less surface. By considering vertical impacts only with the full-chain analytical prescription from Kurosawa & Takada (2019), we provide lower bounds for the ejected mass for basalt, granite, iron and water-rich target objects across the radii range 10^{0-3} km. Our use of the full-chain prescription as opposed to physical experiments or hydrocode simulations allows us to quickly sample a wide range of parameter space appropriate to white dwarf planetary systems. Our numerical results could be used in future studies to constrain freshly-generated small debris reservoirs around white dwarfs given a particular planetary system architecture, bombardment history, and impact geometries.


title: First extraction of the $\Lambda$ polarising fragmentation function from Belle $e^+e^-$ data

abstract: We present a thorough phenomenological analysis of the experimental data from Belle Collaboration for the transverse $\Lambda$ and $\bar\Lambda$ polarisation, measured in $e^+e^-$ annihilation processes, for the case of inclusive (plus a jet) and associated production with a light charged hadron. This allows for the first ever extraction of the quark polarising fragmentation function for a $\Lambda$ hyperon, a transverse momentum dependent distribution giving the probability that an unpolarised quark fragments into a transversely polarised spin-1/2 hadron.


title: Turbulent mixing of r-process elements in the Milky Way

abstract: We study turbulent gas diffusion affects on $r$-process abundances in Milky Way stars, by a combination of an analytical approach and a Monte Carlo simulation. Higher $r$-process event rates and faster diffusion, lead to more efficient mixing corresponding to a reduced scatter of $r$-process abundances and causing $r$-process enriched stars to start appearing at lower metallicities. We use three independent observations to constrain the model parameters: (i) the scatter of radioactively stable $r$-process element abundances, (ii) the largest $r$-process enrichment values observed in any solar neighborhood stars and (iii) the isotope abundance ratios of different radioactive $r$-process elements ($^{244}$Pu/$^{238}$U and $^{247}$Cm/$^{238}$U) at the early solar system as compared to their formation. Our results indicate that the Galactic $r$-process rate and the diffusion coefficient are respectively $r<4\times 10^{-5}\mbox{ yr}^{-1}, D>0.1 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ ($r<4\times 10^{-6}\mbox{ yr}^{-1}, D>0.5 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ for collapsars or similarly prolific $r$-process sources) with allowed values satisfying an approximate anti-correlation such that $D\approx r^{-2/3}$, implying that the time between two $r$-process events that enrich the same location in the Galaxy, is $\tau_{\rm mix}\approx 100-200\mbox{ Myr}$. This suggests that a fraction of $\sim 0.8$ ($\sim 0.5$) of the observed $^{247}$Cm ($^{244}$Pu) abundance is dominated by one $r$-process event in the early solar system. Radioactively stable element abundances are dominated by contributions from $\sim 10$ different events in the early solar system. For metal poor stars (with [Fe/H]$\lesssim -2$), their $r$-process abundances are dominated by either a single or several events, depending on the star formation history.


title: Improved Simulation of Quantum Circuits by Fewer Gaussian Eliminations

abstract: We show that the cost of strong simulation of quantum circuits using $t$ $T$ gate magic states exhibits non-trivial reductions on its upper bound for $t=1$, $t=2$, $t=3$, and $t=6$ with odd-prime-qudits. This agrees with previous numerical bounds found for qubits. We define simulation cost by the number of terms that require Gaussian elimination of a $t \times t$ matrix and so capture the cost of simulation methods that proceed by computing stabilizer inner products or evaluating quadratic Gauss sums. Prior numerical searchs for qubits were unable to converge beyond $t=7$. We effectively increase the space searched for these non-trivial reductions by $>10^{10^4}$ and extend the bounds to $t=14$ for qutrits. This is accomplished by using the Wigner-Weyl-Moyal formalism to algebraically find bounds instead of relying on numerics. We find a new reduction in the upper bound from the $12$-qutrit magic state of ${3^{\sim 0.469t}}$, which improves on the bound obtained from the $6$-qutrit magic state of ${3^{\sim 0.482t}}$.


title: When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

abstract: In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rotation. We study the relative contributions of galactic rotation and local noncircular motions to the circulation of gas, $\Gamma$, a macroscopic measure of local rotation, defined as the line integral of the velocity field around a closed path. We measure the circulation distribution as a function of spatial scale in a set of simulated disk galaxies and we model the velocity field as the sum of galactic rotation and a Gaussian random field. The random field is parameterized by a broken power law in Fourier space, with a break at the scale $\lambda_c$. We define the spatial scale $\lambda_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $\Gamma$. For our simulated galaxies, the gas dynamics at the scale of molecular clouds is usually dominated by noncircular motions, but in the center of galactic disks galactic rotation is still relevant. Our model shows that the transfer of rotation from large scales breaks at the scale $\lambda_c$ and this transition is necessary to reproduce the circulation distribution. We find that $\lambda_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.


title: The Strength of the Dynamical Spiral Perturbation in the Galactic Disk

abstract: The mean Galactocentric radial velocities $\langle v_{R}\rangle(R,\varphi)$ of luminous red giant stars within the mid-plane of the Milky Way reveal a spiral signature, which could plausibly reflect the response to a non-axisymmetric perturbation of the gravitational potential in the Galactic disk. We apply a simple steady-state toy model of a logarithmic spiral to interpret these observations, and find a good qualitative and quantitative match. Presuming that the amplitude of the gravitational potential perturbation is proportionate to that in the disk's surface mass density, we estimate the surface mass density amplitude to be $\Sigma_{\rm max} (R_{\odot})\approx 5.5\,\rm M_{\odot}\,pc^{-2}$ at the solar radius when choosing a fixed pattern speed of $\Omega_{\mathrm p}=12\,\rm km\,s^{-1}\,kpc^{-1}$. Combined with the local disk density, this implies a surface mass density contrast between the arm and inter-arm regions of approximately $\pm 10\%$ at the solar radius, with an increases towards larger radii. Our model constrains the pitch angle of the dynamical spiral arms to be approximately $12^{\circ}$.


title: Maximum velocity quantum circuits

abstract: We consider the long-time limit of out-of-time-order correlators (OTOCs) in two classes of quantum lattice models with time evolution governed by local unitary quantum circuits and maximal butterfly velocity $v_{B} = 1$. Using a transfer matrix approach, we present analytic results for the long-time value of the OTOC on and inside the light cone. First, we consider `dual-unitary' circuits with various levels of ergodicity, including the integrable and non-integrable kicked Ising model, where we show exponential decay away from the light cone and relate both the decay rate and the long-time value to those of the correlation functions. Second, we consider a class of kicked XY models similar to the integrable kicked Ising model, again satisfying $v_{B}=1$, highlighting that maximal butterfly velocity is not exclusive to dual-unitary circuits.

categories: quant-ph cond-mat.stat-mech nlin.CD nlin.SI


title: $\Delta(27)$ framework for cobimaximal neutrino mixing models

abstract: We propose a simple framework based on $\Delta(27)$ that leads to the successful cobimaximal lepton mixing ansatz, thus providing a predictive explanation for leptonic mixing observables. We explore first the effective neutrino mass operators, then present a specific model realization based on type I seesaw, and also propose a model with radiative 1-loop seesaw which features viable dark matter candidates.


title: 11D Supergravity on $AdS_4 \times S^7$ versus $AdS_7 \times S^4$

abstract: The maximally supersymmetric Freund-Rubin vacua for eleven dimensional supergravity, namely $AdS_4 \times S^7$ and $AdS_7 \times S^4$, admit an analytic continuation to $S^4 \times S^7$. From the full harmonic expansions on $S^4 \times S^7$, it is shown that by analytical continuation to either $AdS_4$, or to $AdS_7$, the detailed structure of the Kaluza-Klein spectrum can be obtained for both vacua in a unified manner. The results are shown to be related by a simple rule which interchanges the spacetime and internal space representations. We also obtain the linearized field equations for the singletons and doubletons but they can be gauged away by fixing certain Stuckelberg shift symmetries inherited from the Kaluza-Klein reduction.


title: A pair of TESS planets spanning the radius valley around the nearby mid-M dwarf LTT 3780

abstract: We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbital periods of $P_b=0.77$ days, $P_c=12.25$ days and sizes $r_{p,b}=1.33\pm 0.07$ R$_{\oplus}$, $r_{p,c}=2.30\pm 0.16$ R$_{\oplus}$, the two planets span the radius valley in period-radius space around low mass stars thus making the system a laboratory to test competing theories of the emergence of the radius valley in that stellar mass regime. By combining 63 precise radial-velocity measurements from HARPS and HARPS-N, we measure planet masses of $m_{p,b}=2.62^{+0.48}_{-0.46}$ M$_{\oplus}$ and $m_{p,c}=8.6^{+1.6}_{-1.3}$ M$_{\oplus}$, which indicates that LTT 3780b has a bulk composition consistent with being Earth-like, while LTT 3780c likely hosts an extended H/He envelope. We show that the recovered planetary masses are consistent with predictions from both photoevaporation and from core-powered mass loss models. The brightness and small size of LTT 3780, along with the measured planetary parameters, render LTT 3780b and c as accessible targets for atmospheric characterization of planets within the same planetary system and spanning the radius valley.


title: UV line driven disc wind as the origin of ultrafast outflows in AGN

abstract: UltraFast Outflows (UFO) are observed in some active galactic nuclei (AGN), with blueshifted and highly ionised Fe-K absorption features. AGN typically have a UV bright accretion flow, so UV line driving is an obvious candidate for launching these winds. However this mechanism requires material with UV opacity, in apparent conflict with the observed high ionisation state of the wind. In this paper we synthesise the X-ray energy spectra resulting from different lines of sight through a state of the art radiation hydrodynamics UV line driven disc wind simulation. We demonstrate that there are some lines of sight which only intercept highly ionised and fast outflowing material. The cooler material required for the UV line driving acceleration is out of the line of sight, close to the disc, shielded from the X-rays by a failed wind. We fit these simulated wind spectra to data from the archetypal UFO source PG 1211+143 and show that they broadly reproduce the depth and velocity of the iron absorption lines seen. This directly demonstrates that UV line driving is a viable mechanism to launch even the fastest UFOs. We simulate microcalorimeter observations of this wind and show that their high energy resolution can resolve the detailed structure in the wind and recover the wind energetics when combined with models which correctly estimate the line formation radius of the wind. New data from microcalorimeters will pave the way for physical predictions of AGN wind feedback in cosmological simulations.


title: Mean-field entanglement transitions in random tree tensor networks

abstract: Entanglement phase transitions in quantum chaotic systems subject to projective measurements and in random tensor networks have emerged as a new class of critical points separating phases with different entanglement scaling. We propose a mean-field theory of such transitions by studying the entanglement properties of random tree tensor networks. As a function of bond dimension, we find a phase transition separating area-law from logarithmic scaling of the entanglement entropy. Using a mapping onto a replica statistical mechanics model defined on a Cayley tree and the cavity method, we analyze the scaling properties of such transitions. Our approach provides a tractable, mean-field-like example of an entanglement transition. We verify our predictions numerically by computing directly the entanglement of random tree tensor network states.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.str-el quant-ph


title: Spin-caloritronic transport in hexagonal graphene nanoflakes

abstract: We investigate the spin-dependent thermoelectric effect of graphene flakes with magnetic edges in the ballistic regime. Employing static, respectively, dynamic mean-field theory we first show that magnetism appears at the zigzag edges for a window of Coulomb interactions that increases significantly with increasing flake size. We then use the Landauer formalism in the framework of the non-equilibrium Green's function method to calculate the spin and charge currents in magnetic hexagonal graphene flakes by varying the temperature of the junction for different flake sizes. While in non-magnetic gated graphene the temperature gradient drives a charge current, we observe a significant spin current for hexagonal graphene flakes with magnetic zigzag edges. Specifically, we show that in the "meta" configuration of a hexagonal flake subject to weak Coulomb interactions, a pure spin current can be driven just by a temperature gradient in a temperature range that is promising for device applications. Bigger flakes are found to yield a bigger window of Coulomb interactions where such spin currents are induced by the magnetic zigzag edges, and larger values of the current.


title: The CARMENES search for exoplanets around M dwarfs. Two planets on the opposite sides of the radius gap transiting the nearby M dwarf LTT 3780

abstract: We present the discovery and characterisation of two transiting planets observed by the Transiting Exoplanet Survey Satellite (TESS) orbiting the nearby (d ~ 22 pc), bright (J ~ 9 mag) M3.5 dwarf LTT 3780 (TOI-732). We confirm both planets and their association with LTT 3780 via ground-based photometry and determine their masses using precise radial velocities measured with the CARMENES spectrograph. Precise stellar parameters determined from CARMENES high resolution spectra confirm that LTT 3780 is a mid-M dwarf with an effective temperature of T_eff = 3360 +\- 51 K, a surface gravity of log(g) = 4.81 +/- 0.04 (cgs), and an iron abundance of [Fe/H] = 0.09 +/- 0.16 dex, with an inferred mass of M_star = 0.379 +/- 0.016 M_sun and a radius of R_star = 0.382 +/- 0.012 R_sun. The ultra-short-period planet LTT 3780 b (P_b = 0.77 d) with a radius of 1.35^{+0.06}_{-0.06} R_earth, a mass of 2.34^{+0.24}_{-0.23} M_earth, and a bulk density of 5.24^{+0.94}_{-0.81} g cm^{-3} joins the population of Earth-size planets with rocky, terrestrial composition. The outer planet, LTT 3780 c, with an orbital period of 12.25 d, radius of 2.42^{+0.10}_{-0.10} R_earth, mass of 6.29^{+0.63}_{-0.61} M_earth, and mean density of 2.45^{+0.44}_{-0.37} g cm^{-3} belongs to the population of dense sub-Neptunes. With the two planets located on opposite sides of the radius gap, this planetary system is an excellent target for testing planetary formation, evolution and atmospheric models. In particular, LTT 3780 c is an ideal object for atmospheric studies with the James Webb Space Telescope.


title: Resonant capture in quadruple stellar systems

abstract: Some quadruple star systems in the hierarchical 2+2 configuration exhibit orbit-orbit resonances between the two compact binaries. We show that the most important resonances occur at period ratios of 1:1, 3:2 and 2:1. We describe the conditions required for capture and show that they can be satisfied at the 3:2 and 2:1 resonances in binaries that migrate significantly in semimajor axis after circularization, probably through magnetic braking or gravitational radiation.


title: Efficient variational contraction of two-dimensional tensor networks with a non-trivial unit cell

abstract: Tensor network states provide an efficient class of states that faithfully capture strongly correlated quantum models and systems in classical statistical mechanics. While tensor networks can now be seen as becoming standard tools in the description of such complex many-body systems, close to optimal variational principles based on such states are less obvious to come by. In this work, we generalize a recently proposed variational uniform matrix product state algorithm for capturing one-dimensional quantum lattices in the thermodynamic limit, to the study of regular two-dimensional tensor networks with a non-trivial unit cell. A key property of the algorithm is a computational effort that scales linearly rather than exponentially in the size of the unit cell. We demonstrate the performance of our approach on the computation of the classical partition functions of the antiferromagnetic Ising model and interacting dimers on the square lattice, as well as of a quantum doped resonating valence bond state.


title: Calculation of Optical Response Functions of Dilute-N GaPAsN Lattice-matched to Si

abstract: Dilute-N GaPAsN alloys have great potential for optoelectronics lattice-matched to Si. However, there is a lack of systematic calculation of the optical response of these alloys. The present paper uses the sp^3d^5s^*s_N tight-binding model to calculate the fullband electronic structure of dilute-N GaPAsN, and then calculate the optical response functions considering direct transitions within the electric dipole approximation. Good agreement is obtained for the dielectric function in comparison to available optical data for dilute nitrides. To achieve this, the sp^3d^5s^* parameters for GaP and GaAs are optimized for their optical properties in comparison to published data, which are then used as the basis for the sp^3d^5s^*s_N parameters for dilute-N GaPN and GaAsN. The calculated absorption between the valence band and the newly formed lowest conduction band of the dilute nitrides increases as the N fraction increases, in agreement with experiments, mainly due to the net increase in their coupling in the entire Brillouin zone, supported by the calculated momentum matrix element in the present work.


title: Higher-order topological band structures

abstract: The interplay of topology and symmetry in a material's band structure may result in various patterns of topological states of different dimensionality on the boundary of a crystal. The protection of these "higher-order" boundary states comes from topology, with constraints imposed by symmetry. We review the bulk-boundary correspondence of topological crystalline band structures, which relates the topology of the bulk band structure to the pattern of the boundary states. Furthermore, recent advances in the K-theoretic classification of topological crystalline band structures are discussed.


title: Cubic Dirac operators and the strange Freudenthal-de Vries formula for colour Lie algebras

abstract: The aim of this paper is to define cubic Dirac operators for colour Lie algebras. We give a necessary and sufficient condition to construct a colour Lie algebra from an $\epsilon$-orthogonal representation of an $\epsilon$-quadratic colour Lie algebra. This is used to prove a strange Freudenthal-de Vries formula for basic colour Lie algebras as well as a Parthasarathy formula for cubic Dirac operators of colour Lie algebras. We calculate the cohomology induced by this Dirac operator, analogously to the algebraic Vogan conjecture proved by Huang and Pand\v{z}i\'c.


title: Central extensions and bounded cohomology

abstract: It was shown by Gersten that a central extension of a finitely generated group is quasi-isometrically trivial provided that its Euler class is bounded. We say that a finitely generated group $G$ satisfies Property QITB (quasi-isometrically trivial implies bounded) if the Euler class of any quasi-isometrically trivial central extension of $G$ is bounded. We exhibit a finitely generated group $G$ which does not satisfy Property QITB. This answers a question by Neumann and Reeves, and provides partial answers to related questions by Wienhard and Blank. We also prove that Property QITB holds for a large class of groups, including amenable groups, right-angled Artin groups, relatively hyperbolic groups with amenable peripheral subgroups, and 3-manifold groups. Finally, we show that Property QITB holds for every finitely presented group if a conjecture by Gromov on bounded primitives of differential forms holds as well.


title: Momentum-dependent flow fluctuations as a hydrodynamic response to initial geometry

abstract: We propose a redefinition of the principal component analysis (PCA) of anisotropic flow that makes it more directly connected to fluctuations of the initial geometry of the system. Then, using state-of-the-art hydrodynamic simulations, we make an explicit connection between flow fluctuations and a cumulant expansion of the initial transverse geometry. In particular, we show that the second principal component of elliptic flow is generated by higher-order cumulants, and therefore probes smaller length scales of the initial state. With this information, it will be possible to put new constraints on properties of the early-time dynamics of a heavy-ion collision, including small-scale structure, as well as properties of the quark-gluon plasma.


title: A ringed pole-on outflow from DO Tauri revealed by ALMA

abstract: We present new ALMA Band 6 observations including the CO(2-1) line and 1.3 mm continuum emission from the surroundings of the young stellar object DO Tauri. The ALMA CO molecular data show three different series of rings at different radial velocities. These rings have radii around 220 au and 800 au. We make individual fits to the rings and note that their centers are aligned with DO Tauri and its optical high-velocity jet. In addition, we notice that the velocity of these structures increases with the separation from the young star. We discuss the data under the hypothesis that the rings represent velocity cuts through three outflowing shells that are possibly driven by a wide-angle wind, dragging the environment material along a direction close to the line of sight (i=19{\deg}). We estimate the dynamical ages, the mass, the momentum and the energy of each individual outflow shell and those of the whole outflow. The results are in agreement with those found in outflows from Class II sources. We make a rough estimate for the size of the jet/wind launching region, which needs to be of <15 au. We report the physical characteristics of DO Tauri's disk continuum emission (almost face-on and with a projected major axis in the north-south direction) and its velocity gradient orientation (north-south), indicative of disk rotation for a 1-2 Msun central star. Finally we show an HST [SII] image of the optical jet and report a measurement of its orientation in the plane of the sky.


title: Decision-Making for Automated Vehicles Using a Hierarchical Behavior-Based Arbitration Scheme

abstract: Behavior planning and decision-making are some of the biggest challenges for highly automated systems. A fully automated vehicle (AV) is confronted with numerous tactical and strategical choices. Most state-of-the-art AV platforms implement tactical and strategical behavior generation using finite state machines. However, these usually result in poor explainability, maintainability and scalability. Research in robotics has raised many architectures to mitigate these problems, most interestingly behavior-based systems and hybrid derivatives. Inspired by these approaches, we propose a hierarchical behavior-based architecture for tactical and strategical behavior generation in automated driving. It is a generalizing and scalable decision-making framework, utilizing modular behavior blocks to compose more complex behaviors in a bottom-up approach. The system is capable of combining a variety of scenario- and methodology-specific solutions, like POMDPs, RRT* or learning-based behavior, into one understandable and traceable architecture. We extend the hierarchical behavior-based arbitration concept to address scenarios where multiple behavior options are applicable but have no clear priority against each other. Then, we formulate the behavior generation stack for automated driving in urban and highway environments, incorporating parking and emergency behaviors as well. Finally, we illustrate our design in an explanatory evaluation.


title: Online Agnostic Boosting via Regret Minimization

abstract: Boosting is a widely used machine learning approach based on the idea of aggregating weak learning rules. While in statistical learning numerous boosting methods exist both in the realizable and agnostic settings, in online learning they exist only in the realizable case. In this work we provide the first agnostic online boosting algorithm; that is, given a weak learner with only marginally-better-than-trivial regret guarantees, our algorithm boosts it to a strong learner with sublinear regret. Our algorithm is based on an abstract (and simple) reduction to online convex optimization, which efficiently converts an arbitrary online convex optimizer to an online booster. Moreover, this reduction extends to the statistical as well as the online realizable settings, thus unifying the 4 cases of statistical/online and agnostic/realizable boosting.


title: Wheel graph strategy for PEV localization of networks

abstract: Investigation of eigenvector localization properties of complex networks is not only important for gaining insight into fundamental network problems such as network centrality measure, spectral partitioning, development of approximation algorithms, but also is crucial for understanding many real-world phenomena such as disease spreading, criticality in brain network dynamics. For a network, an eigenvector is said to be localized when most of its components take value near to zero, with a few components taking very high values. In this article, we devise a methodology to construct a principal eigenvector (PEV) localized network from a given input network. The methodology relies on adding a small component having a wheel graph to the given input network. By extensive numerical simulation and an analytical formulation based on the largest eigenvalue of the input network, we compute the size of the wheel graph required to localize the PEV of the combined network. Using the susceptible-infected-susceptible model, we demonstrate the success of this method for various models and real-world networks consider as input networks. We show that on such PEV localized networks, the disease gets localized within a small region of the network structure before the outbreaks. The study is relevant in controlling spreading processes on complex systems represented by networks.


title: Non-Hermitian topology of one-dimensional spin-torque oscillator arrays

abstract: Magnetic systems have been extensively studied both from a fundamental physics perspective and as building blocks for a variety of applications. Their topological properties, in particular those of excitations, remain relatively unexplored due to their inherently dissipative nature. The recent introduction of non-Hermitian topological classifications opens up new opportunities for engineering topological phases in dissipative systems. Here, we propose a magnonic realization of a non-Hermitian topological system. A crucial ingredient of our proposal is the injection of spin current into the magnetic system, which alters and can even change the sign of terms describing dissipation. We show that the magnetic dynamics of an array of spin-torque oscillators can be mapped onto a non-Hermitian Su-Schrieffer-Heeger model exhibiting topologically protected edge states. Using exact diagonalization of the linearized dynamics and numerical solutions of the non-linear equations of motion, we find that a topological magnonic phase can be accessed by tuning the spin current injected into the array. In the topologically nontrivial regime, a single spin-torque oscillator on the edge of the array is driven into auto-oscillation and emits a microwave signal, while the bulk oscillators remain inactive. Our findings have practical utility for memory devices and spintronics neural networks relying on spin-torque oscillators as constituent units.


title: Which Pull Requests Get Accepted and Why? A study of popular NPM Packages

abstract: Background: Pull Request (PR) Integrators often face challenges in terms of multiple concurrent PRs, so the ability to gauge which of the PRs will get accepted can help them balance their workload. PR creators would benefit from knowing if certain characteristics of their PRs may increase the chances of acceptance. Aim: We modeled the probability that a PR will be accepted within a month after creation using a Random Forest model utilizing 50 predictors representing properties of the author, PR, and the project to which PR is submitted. Method: 483,988 PRs from 4218 popular NPM packages were analysed and we selected a subset of 14 predictors sufficient for a tuned Random Forest model to reach high accuracy. Result: An AUC-ROC value of 0.95 was achieved predicting PR acceptance. The model excluding PR properties that change after submission gave an AUC-ROC value of 0.89. We tested the utility of our model in practical scenarios by training it with historical data for the NPM package \textit{bootstrap} and predicting if the PRs submitted in future will be accepted. This gave us an AUC-ROC value of 0.94 with all 14 predictors, and 0.77 excluding PR properties that change after its creation. Conclusion: PR integrators can use our model for a highly accurate assessment of the quality of the open PRs and PR creators may benefit from the model by understanding which characteristics of their PRs may be undesirable from the integrators' perspective. The model can be implemented as a tool, which we plan to do as a future work.


title: Efficient algorithms for the Potts model on small-set expanders

abstract: An emerging trend in approximate counting is to show that certain `low-temperature' problems are easy on typical instances, despite worst-case hardness results. For the class of regular graphs one usually shows that expansion can be exploited algorithmically, and since random regular graphs are good expanders with high probability the problem is typically tractable. Inspired by approaches used in subexponential-time algorithms for Unique Games, we develop an approximation algorithm for the partition function of the ferromagnetic Potts model on graphs with a small-set expansion condition. In such graphs it may not suffice to explore the state space of the model close to ground states, and a novel feature of our method is to efficiently find a larger set of `pseudo-ground states' such that it is enough to explore the model around each pseudo-ground state.


title: Towards Probability-based Safety Verification of Systems with Components from Machine Learning

abstract: Machine learning (ML) has recently created many new success stories. Hence, there is a strong motivation to use ML technology in software-intensive systems, including safety-critical systems. This raises the issue of safety verification of MLbased systems, which is currently thought to be infeasible or, at least, very hard. We think that it requires taking into account specific properties of ML technology such as: (i) Most ML approaches are inductive, which is both their power and their source of error. (ii) Neural networks (NN) resulting from deep learning are at the current state of the art not transparent. Consequently, there will always be errors remaining and, at least for deep NNs (DNNs), verification of their internal structure is extremely hard. In general, safety engineering cannot provide full guarantees that no harm will ever occur. That is why probabilities are used, e.g., for specifying a risk or a Tolerable Hazard Rate (THR). In this vision paper, we propose verification based on probabilities of errors both estimated by controlled experiments and output by the inductively learned classifier itself. Generalization error bounds may propagate to the probabilities of a hazard, which must not exceed a THR. As a result, the quantitatively determined bound on the probability of a classification error of an ML component in a safety-critical system contributes in a well-defined way to the latter's overall safety verification.


title: Real-World Human-Robot Collaborative Reinforcement Learning

abstract: The intuitive collaboration of humans and intelligent robots (embodied AI) in the real-world is an essential objective for many desirable applications of robotics. Whilst there is much research regarding explicit communication, we focus on how humans and robots interact implicitly, on motor adaptation level. We present a real-world setup of a human-robot collaborative maze game, designed to be non-trivial and only solvable through collaboration, by limiting the actions to rotations of two orthogonal axes, and assigning each axes to one player. This results in neither the human nor the agent being able to solve the game on their own. We use deep reinforcement learning for the control of the robotic agent, and achieve results within 30 minutes of real-world play, without any type of pre-training. We then use this setup to perform systematic experiments on human/agent behaviour and adaptation when co-learning a policy for the collaborative game. We present results on how co-policy learning occurs over time between the human and the robotic agent resulting in each participant's agent serving as a representation of how they would play the game. This allows us to relate a person's success when playing with different agents than their own, by comparing the policy of the agent with that of their own agent.


title: Reinforcement co-Learning of Deep and Spiking Neural Networks for Energy-Efficient Mapless Navigation with Neuromorphic Hardware

abstract: Energy-efficient mapless navigation is crucial for mobile robots as they explore unknown environments with limited on-board resources. Although the recent deep reinforcement learning (DRL) approaches have been successfully applied to navigation, their high energy consumption limits their use in several robotic applications. Here, we propose a neuromorphic approach that combines the energy-efficiency of spiking neural networks with the optimality of DRL and benchmark it in learning control policies for mapless navigation. Our hybrid framework, spiking deep deterministic policy gradient (SDDPG), consists of a spiking actor network (SAN) and a deep critic network, where the two networks were trained jointly using gradient descent. The co-learning enabled synergistic information exchange between the two networks, allowing them to overcome each other's limitations through a shared representation learning. To evaluate our approach, we deployed the trained SAN on Intel's Loihi neuromorphic processor. When validated on simulated and real-world complex environments, our method on Loihi consumed 75 times less energy per inference as compared to DDPG on Jetson TX2, and also exhibited a higher rate of successful navigation to the goal, which ranged from 1% to 4.2% and depended on the forward-propagation timestep size. These results reinforce our ongoing efforts to design brain-inspired algorithms for controlling autonomous robots with neuromorphic hardware.


title: Quantum Correlations Dynamics In Two Coupled Semiconductor InAs Quantum Dots

abstract: We investigate the dynamics of quantum discord and concurrence between two excitonic qubits placed inside two coupled semiconductor quantum dots independently interacting with dephasing reservoirs. We explore their behavior against the dimensionless time and the temperature in both Markovian and non-Markovian environments. Moreover, we analyze the external electric field effects and the F\"orster interaction effects on these correlations. We show that, although the quantum correlations amount is strongly influenced by the variation of the electric field and the F\"orster interaction, their non-Markovian behavior is still preserved under the variation of these two parameters. Furthermore, we show that for large values of temperature and dimensionless time, unlike concurrence which vanishes, nonzero discord can still be observed.


title: Time-domain sparsity promoting least-squares reverse time migration with source estimation

abstract: Least-squares reverse time migration is well-known for its capability to generate artifact-free true-amplitude subsurface images through fitting observed data in the least-squares sense. However, when applied to realistic imaging problems, this approach is faced with issues related to overfitting and excessive computational costs induced by many wave-equation solves. The fact that the source function is unknown complicates this situation even further. Motivated by recent results in stochastic optimization and transform-domain sparsity-promotion, we demonstrate that the computational costs of inversion can be reduced significantly while avoiding imaging artifacts and restoring amplitudes. While powerful, these new approaches do require accurate information on the source-time function, which is often lacking. Without this information, the imaging quality deteriorates rapidly. We address this issue by presenting an approach where the source-time function is estimated on the fly through a technique known as variable projection. Aside from introducing negligible computational overhead, the proposed method is shown to perform well on imaging problems with noisy data and problems that involve complex settings such as salt. In either case, the presented method produces high resolution high-amplitude fidelity images including an estimates for the source-time function. In addition, due to its use of stochastic optimization, we arrive at these images at roughly one to two times the cost of conventional reverse time migration involving all data.


title: A note on a Bonnet-Myers type diameter bound for graphs with positive entropic Ricci curvature

abstract: An equivalent definition of entropic Ricci curvature on discrete spaces was given in terms of the global gradient estimate. With a particular choice of the density function $\rho$, we obtain a localized gradient estimate, which in turns allow us to derive a Bonnet-Myers type diameter bound for graphs with positive entropic Ricci curvature. However, the case of the hypercubes indicates that the bound may be not optimal (where $\theta$ is chosen to be logarithmic mean by default). If $\theta$ is arithmetic mean, the Bakry-\'Emery criterion can be recovered and the diameter bound is optimal as it can be attained by the hypercubes.


title: A note on Kawaguchi-Silverman conjecture

abstract: We collect some results on endomorphisms on projective varieties related with the Kawaguchi-Silverman conjecture. We discuss certain condition on automorphism groups of projective varieties and positivity conditions on leading real eigendivisors of self-morphisms. We prove Kawaguchi-Silverman conjecture for endomorphisms on projective bundles on a smooth Fano variety of Picard number one. In the last section, we discuss endomorphisms and augmented base loci of their eigendivisors.


title: Multichannel Singing Voice Separation by Deep Neural Network Informed DOA Constrained CNMF

abstract: This work addresses the problem of multichannel source separation combining two powerful approaches, multichannel spectral factorization with recent monophonic deep-learning (DL) based spectrum inference. Individual source spectra at different channels are estimated with a Masker-Denoiser Twin Network (MaD TwinNet), able to model long-term temporal patterns of a musical piece. The monophonic source spectrograms are used within a spatial covariance mixing model based on Complex Non-Negative Matrix Factorization (CNMF) that predicts the spatial characteristics of each source. The proposed framework is evaluated on the task of singing voice separation with a large multichannel dataset. Experimental results show that our joint DL+CNMF method outperforms both the individual monophonic DL-based separation and the multichannel CNMF baseline methods.


title: Understanding Contexts Inside Robot and Human Manipulation Tasks through a Vision-Language Model and Ontology System in a Video Stream

abstract: Manipulation tasks in daily life, such as pouring water, unfold intentionally under specialized manipulation contexts. Being able to process contextual knowledge in these Activities of Daily Living (ADLs) over time can help us understand manipulation intentions, which are essential for an intelligent robot to transition smoothly between various manipulation actions. In this paper, to model the intended concepts of manipulation, we present a vision dataset under a strictly constrained knowledge domain for both robot and human manipulations, where manipulation concepts and relations are stored by an ontology system in a taxonomic manner. Furthermore, we propose a scheme to generate a combination of visual attentions and an evolving knowledge graph filled with commonsense knowledge. Our scheme works with real-world camera streams and fuses an attention-based Vision-Language model with the ontology system. The experimental results demonstrate that the proposed scheme can successfully represent the evolution of an intended object manipulation procedure for both robots and humans. The proposed scheme allows the robot to mimic human-like intentional behaviors by watching real-time videos. We aim to develop this scheme further for real-world robot intelligence in Human-Robot Interaction.


title: Dopant levels in large nanocrystals using stochastic optimally tuned range-separated hybrid density functional theory

abstract: We apply a stochastic version of an optimally tuned range-separated hybrid functional to provide insight on the electronic properties of P- and B- doped Si nanocrystals of experimentally relevant sizes. We show that we can use the range-separation parameter for undoped systems to calculate accurate results for dopant activation energies. We apply this strategy for tuning functionals to study doped nanocrystals up to 2.5 nm in diameter at the hybrid functional level. In this confinement regime, the P- and B- dopants have large activation energies and have strongly localized states that lie deep within the energy gaps. Structural relaxation plays a greater role for B-substituted dopants and contributes to the increase in activation energy when the B dopant is near the nanocrystal surface.


title: Large mass minimizers for isoperimetric problems with integrable nonlocal potentials

abstract: This paper is concerned with volume-constrained minimization problems derived from Gamow's liquid drop model for the atomic nucleus, involving the competition of a perimeter term and repulsive nonlocal potentials. We consider a large class of potentials, given by general radial nonnegative kernels which are integrable on $\mathbb{R}^n$, such as Bessel potentials, and study the behavior of the problem for large masses (i.e., volumes). Contrarily to the small mass case, where the nonlocal term becomes negligible compared to the perimeter, here the nonlocal term explodes compared to it. However, using the integrability of those kernels, we rewrite the problem as the minimization of the difference between the classical perimeter and a nonlocal perimeter, which converges to a multiple of the classical perimeter as the mass goes to infinity. Renormalizing to a fixed volume, we show that, if the first moment of the kernels is smaller than an explicit threshold, the problem admits minimizers of arbitrarily large mass, which contrasts with the usual case of Riesz potentials. In addition, we prove that, any sequence of minimizers converges to the ball as the mass goes to infinity. Finally, we study the stability of the ball, and show that our threshold on the first moment of the kernels is sharp in the sense that large balls go from stable to unstable. A direct consequence of the instability of large balls above this threshold is that there exist nontrivial compactly supported kernels for which the problems admit minimizers which are not balls, that is, symmetry breaking occurs.


title: Discrimination and estimation of incoherent sources under misalignment

abstract: Spatially resolving two incoherent point sources whose separation is well below the diffraction limit dictated by classical optics has recently been shown possible using techniques that decompose the incoming radiation into orthogonal transverse modes. Such a demultiplexing procedure, however, must be perfectly calibrated to the transverse profile of the incoming light as any misalignment of the modes effectively restores the diffraction limit for small source separations. We study by how much can one mitigate such an effect at the level of measurement which, after being imperfectly demultiplexed due to inevitable misalignment, may still be partially corrected by linearly transforming the relevant dominating transverse modes. We consider two complementary tasks: the estimation of the separation between the two sources and the discrimination between one and two incoherent point sources. We show that, although one cannot fully restore super-resolving powers even when the value of the misalignment is perfectly known its negative impact on the ultimate sensitivity can be significantly reduced. In the case of estimation we analytically determine the exact relation between the minimal resolvable separation as a function of misalignment whereas for discrimination we analytically determine the relation between misalignment and the probability of error, as well as numerically determine how the latter scales in the limit of long interrogation times.


title: Dust dynamics and vertical settling in gravitoturbulent protoplanetary discs

abstract: Gravitational instability (GI) controls the dynamics of young massive protoplanetary discs. Apart from facilitating gas accretion on to the central protostar, it must also impact on the process of planet formation: directly through fragmentation, and indirectly through the turbulent concentration of small solids. To understand the latter process, it is essential to determine the dust dynamics in such a turbulent flow. For that purpose, we conduct a series of 3D shearing box simulations of coupled gas and dust, including the gas's self-gravity and scanning a range of Stokes numbers, from 0.001 to ~0.2. First, we show that the vertical settling of dust in the midplane is significantly impeded by gravitoturbulence, with the dust scale-height roughly 0.6 times the gas scale height for centimetre grains. This is a result of the strong vertical diffusion issuing from (a) small-scale inertial-wave turbulence feeding off the GI spiral waves and (b) the larger-scale vertical circulations that naturally accompany the spirals. Second, we show that at R=50 AU concentration events involving sub-metre particles and yielding order 1 dust to gas ratios are rare and last for less than an orbit. Moreover, dust concentration is less efficient in 3D than in 2D simulations. We conclude that GI is not especially prone to the turbulent accumulation of dust grains. Finally, the large dust scale-height measured in simulations could be, in the future, compared with that of edge-on discs seen by ALMA, thus aiding detection and characterisation of GI in real systems.


title: Long-term Spatial Modeling for Characteristics of Extreme Heat Events

abstract: There is increasing evidence that global warming manifests itself in more frequent warm days and that heat waves will become more frequent. Presently, a formal definition of a heat wave is not agreed upon in the literature. To avoid this debate, we consider extreme heat events, which, at a given location, are well-defined as a run of consecutive days above an associated local threshold. Characteristics of EHEs are of primary interest, such as incidence and duration, as well as the magnitude of the average exceedance and maximum exceedance above the threshold during the EHE. Using approximately 60-year time series of daily maximum temperature data collected at 18 locations in a given region, we propose a spatio-temporal model to study the characteristics of EHEs over time. The model enables prediction of the behavior of EHE characteristics at unobserved locations within the region. Specifically, our approach employs a two-state space-time model for EHEs with local thresholds where one state defines above threshold daily maximum temperatures and the other below threshold temperatures. We show that our model is able to recover the EHE characteristics of interest and outperforms a corresponding autoregressive model that ignores thresholds based on out-of-sample prediction.


title: A Robust Functional EM Algorithm for Incomplete Panel Count Data

abstract: Panel count data describes aggregated counts of recurrent events observed at discrete time points. To understand dynamics of health behaviors, the field of quantitative behavioral research has evolved to increasingly rely upon panel count data collected via multiple self reports, for example, about frequencies of smoking using in-the-moment surveys on mobile devices. However, missing reports are common and present a major barrier to downstream statistical learning. As a first step, under a missing completely at random assumption (MCAR), we propose a simple yet widely applicable functional EM algorithm to estimate the counting process mean function, which is of central interest to behavioral scientists. The proposed approach wraps several popular panel count inference methods, seamlessly deals with incomplete counts and is robust to misspecification of the Poisson process assumption. Theoretical analysis of the proposed algorithm provides finite-sample guarantees by expanding parametric EM theory to our general non-parametric setting. We illustrate the utility of the proposed algorithm through numerical experiments and an analysis of smoking cessation data. We also discuss useful extensions to address deviations from the MCAR assumption and covariate effects.


title: Central moments of the free energy of the O'Connell-Yor polymer

abstract: Sepp\"al\"ainen and Valk\'o showed in \cite{SV} that for a suitable choice of parameters, the variance growth of the free energy of the stationary O'Connell-Yor polymer is governed by the exponent $2/3$, characteristic of models in the KPZ universality class. We develop exact formulas based on Gaussian integration by parts to relate the cumulants of the free energy, $\log Z_{n,t}^\theta$, to expectations of products of quenched cumulants of the time of the first jump from the boundary into the system, $s_0$. We then use these formulas to obtain estimates for the $k$-th central moment of $\log Z_{n,t}^\theta$ as well as the $k$-th annealed moment of $s_0$ for $k> 2$, with nearly optimal exponents $(1/3)k+\epsilon$ and $(2/3)k+\epsilon$, respectively.


title: A Semi-supervised Graph Attentive Network for Financial Fraud Detection

abstract: With the rapid growth of financial services, fraud detection has been a very important problem to guarantee a healthy environment for both users and providers. Conventional solutions for fraud detection mainly use some rule-based methods or distract some features manually to perform prediction. However, in financial services, users have rich interactions and they themselves always show multifaceted information. These data form a large multiview network, which is not fully exploited by conventional methods. Additionally, among the network, only very few of the users are labelled, which also poses a great challenge for only utilizing labeled data to achieve a satisfied performance on fraud detection. To address the problem, we expand the labeled data through their social relations to get the unlabeled data and propose a semi-supervised attentive graph neural network, namedSemiGNN to utilize the multi-view labeled and unlabeled data for fraud detection. Moreover, we propose a hierarchical attention mechanism to better correlate different neighbors and different views. Simultaneously, the attention mechanism can make the model interpretable and tell what are the important factors for the fraud and why the users are predicted as fraud. Experimentally, we conduct the prediction task on the users of Alipay, one of the largest third-party online and offline cashless payment platform serving more than 4 hundreds of million users in China. By utilizing the social relations and the user attributes, our method can achieve a better accuracy compared with the state-of-the-art methods on two tasks. Moreover, the interpretable results also give interesting intuitions regarding the tasks.


title: Volume Above Distance Below

abstract: Given a pair of metric tensors $g_1 \ge g_0$ on a Riemannian manifold, $M$, it is well known that $\operatorname{Vol}_1(M) \ge \operatorname{Vol}_0(M)$. Furthermore one has rigidity: the volumes are equal if and only if the metric tensors are the same $g_1=g_0$. Here we prove that if $g_j \ge g_0$ and $\operatorname{Vol}_1(M)\to \operatorname{Vol}_0(M)$ then $(M,g_j)$ converge to $(M,g_0)$ in the volume preserving intrinsic flat sense. Well known examples demonstrate that one need not obtain smooth, $C^0$, Lipschitz, or even Gromov-Hausdorff convergence in this setting. Our theorem may also be applied as a tool towards proving other open conjectures concerning the geometric stability of a variety of rigidity theorems in Riemannian geometry. To complete our proof, we provide a novel way of estimating the intrinsic flat distance between Riemannian manifolds which is interesting in its own right.


title: Automatic continuity for homeomorphism groups of noncompact manifolds

abstract: We extend the proof of automatic continuity for homeomorphism groups of manifolds to non-compact manifolds and manifolds with marked points and their mapping class groups. Specifically, we show that, for any manifold $M$ homeomorphic to the interior of a compact manifold, and a set $X \subset M$ homeomorphic to the union of a Cantor set and finite set, the relative homeomorphism group $\mathrm{Homeo}(M, X)$ and the mapping class group $\mathrm{Homeo}(M, X)/\mathrm{Homeo}_0(M,X)$ have the property that any homomorphism from such a group to any separable topological group is necessarily continuous.


title: LiDARNet: A Boundary-Aware Domain Adaptation Model for Point Cloud Semantic Segmentation

abstract: We present a boundary-aware domain adaptation model for LiDAR scan full-scene semantic segmentation (LiDARNet). Our model can extract both the domain private features and the domain shared features with a two-branch structure. We embedded Gated-SCNN into the segmentor component of LiDARNet to learn boundary information while learning to predict full-scene semantic segmentation labels. Moreover, we further reduce the domain gap by inducing the model to learn a mapping between two domains using the domain shared and private features. Additionally, we introduce a new dataset (SemanticUSL\footnote{The access address of SemanticUSL:\url{https://unmannedlab.github.io/research/SemanticUSL}}) for domain adaptation for LiDAR point cloud semantic segmentation. The dataset has the same data format and ontology as SemanticKITTI. We conducted experiments on real-world datasets SemanticKITTI, SemanticPOSS, and SemanticUSL, which have differences in channel distributions, reflectivity distributions, diversity of scenes, and sensors setup. Using our approach, we can get a single projection-based LiDAR full-scene semantic segmentation model working on both domains. Our model can keep almost the same performance on the source domain after adaptation and get an 8\%-22\% mIoU performance increase in the target domain.


title: Fast and efficient deterministic quantum state transfer between two remote mechanical resonators

abstract: The main challenge in deterministic quantum state transfer in long-distance quantum communications is the transmission losses in the communication channel. To overcome this limitation, here we use the adiabatic theorem and find a lossless evolution path between two remote mechanical modes. By adiabatic variation of the effective coupling strengths between the two nodes and the intermediate optical channel modes, we engineer a transmission path for the quantum state transfer that is decoupled from the decaying fiber modes. Using our proposed method we show that one obtains a quantum state transfer with high efficiency. Furthermore, to bypass the slow nature of the adiabatic process and its sensitivity to the mechanical damping and noise as well as the strength of the driving pulses, we develop the shortcut to adiabatic passage protocol for our proposed quantum state transfer. Our results show that the shortcut to adiabaticity provides an efficient and fast quantum state transfer even for small values of the coupling strength. We show that the performance of our protocol for long-distance quantum communications remains efficient for transferring the quantum states between two remote mechanical resonators being hundred meters apart.


title: Deep Survival Machines: Fully Parametric Survival Regression and Representation Learning for Censored Data with Competing Risks

abstract: We describe a new approach to estimating relative risks in time-to-event prediction problems with censored data in a fully parametric manner. Our approach does not require making strong assumptions of constant proportional hazard of the underlying survival distribution, as required by the Cox-proportional hazard model. By jointly learning deep nonlinear representations of the input covariates, we demonstrate the benefits of our approach when used to estimate survival risks through extensive experimentation on multiple real world datasets with different levels of censoring. We further demonstrate advantages of our model in the competing risks scenario. To the best of our knowledge, this is the first work involving fully parametric estimation of survival times with competing risks in the presence of censoring.


title: Shortcut to geostrophy in wave-driven rotating turbulence: the quartetic instability

abstract: We report on laboratory experiments of wave-driven rotating turbulence. A set of wavemakers produces inertial-wave beams that interact nonlinearly in the central region of a water tank mounted on a rotating platform. The forcing thus injects energy into inertial waves only. For moderate forcing amplitude, part of the energy of the forced inertial waves is transferred to subharmonic waves, through a standard triadic resonance instability. This first step is broadly in line with the theory of weak turbulence. Surprisingly however, stronger forcing does not lead to an inertial-wave turbulence regime. Instead, most of the kinetic energy condenses into a vertically invariant geostrophic flow, even though the latter is unforced. We show that resonant quartets of inertial waves can trigger an instability -- the ``quartetic instability'' -- that leads to such spontaneous emergence of geostrophy. In the present experiment, this instability sets in as a secondary instability of the classical triadic instability.


title: A Study of the Fundamental Performance Characteristics of GPUs and CPUs for Database Analytics (Extended Version)

abstract: There has been significant amount of excitement and recent work on GPU-based database systems. Previous work has claimed that these systems can perform orders of magnitude better than CPU-based database systems on analytical workloads such as those found in decision support and business intelligence applications. A hardware expert would view these claims with suspicion. Given the general notion that database operators are memory-bandwidth bound, one would expect the maximum gain to be roughly equal to the ratio of the memory bandwidth of GPU to that of CPU. In this paper, we adopt a model-based approach to understand when and why the performance gains of running queries on GPUs vs on CPUs vary from the bandwidth ratio (which is roughly 16x on modern hardware). We propose Crystal, a library of parallel routines that can be combined together to run full SQL queries on a GPU with minimal materialization overhead. We implement individual query operators to show that while the speedups for selection, projection, and sorts are near the bandwidth ratio, joins achieve less speedup due to differences in hardware capabilities. Interestingly, we show on a popular analytical workload that full query performance gain from running on GPU exceeds the bandwidth ratio despite individual operators having speedup less than bandwidth ratio, as a result of limitations of vectorizing chained operators on CPUs, resulting in a 25x speedup for GPUs over CPUs on the benchmark.


title: Spatially resolved velocity structure in jets of DF Tau and UY Aur A

abstract: Young stars accrete mass and angular momentum from their circumstellar disks. Some of them also drive outflows, which can be distinguished in optical forbidden emission lines (FELs). We analyze a sample of binary T Tauri stars observed with long-slit spectroscopy by the Hubble Space Telescope (HST), searching for spatially resolved outflows. We detect resolved [O I] emission in two cases out of twenty one. In DF Tau we resolve high and medium velocity outflows in a jet and counterjet out to 60 au. The outflows are accelerated within the inner 12 au and retain a constant speed thereafter. In UY Aur, we detect a blue- and a red-shifted outflow from UY Aur A, as well as a blue-shifted jet from UY Aur B. All of these features have been seen in [Fe II] with data taken ten years apart indicating that the underlying outflow pattern is stable on these timescales.


title: A Correlation Method for Deriving UA Intra-Molecular Potentials from AA Molecular Dynamics Simulations: Application to Alkanes

abstract: United-Atom (UA) force fields are usually constructed using a Boltzmann-inverse method based on distributions obtained from Monte-Carlo simulations. A new method of constructing UA force fields from All-Atom (AA) molecular dynamics simulations is proposed. In this method, one determines time correlations between oscillators of the same type: between a CC bond and the connected CH bonds; between a valence angle CCC and all valence angles sharing at least one hydrogen with it; between a dihedral angle CCCC and the dihedral angles sharing the two central carbons with it. In the case of no correlation between lengths or angles, energies of the oscillators are independent. In this case, the AA and UA components of the force fields are identical. When a correlation is total, the UA potential of an object is the sum of the energies of the oscillators coupled to it. Partial correlations are also possible. Several kinematic matching tests of the AA molecule with the corresponding UA molecule are considered: mean periods of vibration of the oscillators; energies associated with presence probabilities; mean end-to-end lengths of the molecules. The method is, as an example, applied to aliphatic molecules, simulated with the OPLS-AA force field.


title: RandomNet: Towards Fully Automatic Neural Architecture Design for Multimodal Learning

abstract: Almost all neural architecture search methods are evaluated in terms of performance (i.e. test accuracy) of the model structures that it finds. Should it be the only metric for a good autoML approach? To examine aspects beyond performance, we propose a set of criteria aimed at evaluating the core of autoML problem: the amount of human intervention required to deploy these methods into real world scenarios. Based on our proposed evaluation checklist, we study the effectiveness of a random search strategy for fully automated multimodal neural architecture search. Compared to traditional methods that rely on manually crafted feature extractors, our method selects each modality from a large search space with minimal human supervision. We show that our proposed random search strategy performs close to the state of the art on the AV-MNIST dataset while meeting the desirable characteristics for a fully automated design process.


title: Assessing Software Defection Prediction Performance: Why Using the Matthews Correlation Coefficient Matters

abstract: Context: There is considerable diversity in the range and design of computational experiments to assess classifiers for software defect prediction. This is particularly so, regarding the choice of classifier performance metrics. Unfortunately some widely used metrics are known to be biased, in particular F1. Objective: We want to understand the extent to which the widespread use of the F1 renders empirical results in software defect prediction unreliable. Method: We searched for defect prediction studies that report both F1 and the Matthews correlation coefficient (MCC). This enabled us to determine the proportion of results that are consistent between both metrics and the proportion that change. Results: Our systematic review identifies 8 studies comprising 4017 pairwise results. Of these results, the direction of the comparison changes in 23% of the cases when the unbiased MCC metric is employed. Conclusion: We find compelling reasons why the choice of classification performance metric matters, specifically the biased and misleading F1 metric should be deprecated.


title: A comparison of computed and experimental neutron diffraction intensity at large momentum for MnO and NiO

abstract: Magnetic neutron scattering measures spin-spin correlations giving information about the long-range spin order as well as the shape of the spin density in magnetic materials. Similarly, detailed first principles calculations directly compute the spin density in materials. In this work, the authors carefully compare experimentally measured magnetic neutron intensities to three levels of theory: density functional theory in two approximations, and fixed-node diffusion Monte Carlo. While each theory performs similarly for the simple antiferromagnet MnO, there are significant differences between density functional theory and diffusion Monte Carlo in NiO. For both materials, fixed-node diffusion Monte Carlo shows the lowest RMS error with respect to experiment for the form factor and the intensity. Through connection of the intensities and form factors to the real-space spin densities, it is shown that diffusion Monte Carlo spin density becomes more diffuse by spreading spin away from the bond directions. This benchmark is of importance when considering the efficacy of an ab initio calculation in capturing both the long-range magnetic correlations and the microscopic spin details of magnetic systems.


title: Variational inference formulation for a model-free simulation of a dynamical system with unknown parameters by a recurrent neural network

abstract: We propose a recurrent neural network for a "model-free" simulation of a dynamical system with unknown parameters without prior knowledge. The deep learning model aims to jointly learn the nonlinear time marching operator and the effects of the unknown parameters from a time series dataset. We assume that the time series data set consists of an ensemble of trajectories for a range of the parameters. The learning task is formulated as a statistical inference problem by considering the unknown parameters as random variables. A latent variable is introduced to model the effects of the unknown parameters, and a variational inference method is employed to simultaneously train probabilistic models for the time marching operator and an approximate posterior distribution for the latent variable. Unlike the classical variational inference, where a factorized distribution is used to approximate the posterior, we employ a feedforward neural network supplemented by an encoder recurrent neural network to develop a more flexible probabilistic model. The approximate posterior distribution makes an inference on a trajectory to identify the effects of the unknown parameters. The time marching operator is approximated by a recurrent neural network, which takes a latent state sampled from the approximate posterior distribution as one of the input variables, to compute the time evolution of the probability distribution conditioned on the latent variable. In the numerical experiments, it is shown that the proposed variational inference model makes a more accurate simulation compared to the standard recurrent neural networks. It is found that the proposed deep learning model is capable of correctly identifying the dimensions of the random parameters and learning a representation of complex time series data.


title: Higgs modes in two-dimensional spatially-indirect exciton condensates

abstract: Higgs modes in condensed matter physics have drawn attention because of analogies to the Higgs bosons of particle physics. Here we use a microscopic time-dependent mean-field theory to study the collective mode spectra of two-dimensional spatially indirect exciton (electron-hole pair) condensates, focusing on the Higgs modes, i.e., those that have a large weight in electron-hole pair amplitude response functions. We find that in the low exciton density (BEC) limit, the dominant Higgs modes of spatially indirect exciton condensates correspond to adding electron-hole pairs that are orthogonal to the condensed pair state. We comment on the previously studied Higgs-like collective excitations of superconductors in light of this finding.