title: Calibrated Adversarial Refinement for Stochastic Semantic Segmentation

abstract: In semantic segmentation tasks, input images can often have more than one plausible interpretation, thus allowing for multiple valid labels. To capture such ambiguities, recent work has explored the use of probabilistic networks that can learn a distribution over predictions. However, these do not necessarily represent the empirical distribution accurately. In this work, we present a strategy for learning a calibrated predictive distribution over semantic maps, where the probability associated with each prediction reflects its ground truth correctness likelihood. To this end, we propose a novel two-stage, cascaded approach for calibrated adversarial refinement: (i) a standard segmentation network is trained with categorical cross entropy to predict a pixelwise probability distribution over semantic classes and (ii) an adversarially trained stochastic network is used to model the inter-pixel correlations to refine the output of the first network into coherent samples. Importantly, to calibrate the refinement network and prevent mode collapse, the expectation of the samples in the second stage is matched to the probabilities predicted in the first. We demonstrate the versatility and robustness of the approach by achieving state-of-the-art results on the multigrader LIDC dataset and on a modified Cityscapes dataset with injected ambiguities. In addition, we show that the core design can be adapted to other tasks requiring learning a calibrated predictive distribution by experimenting on a toy regression dataset. We provide an open source implementation of our method at https://github.com/EliasKassapis/CARSSS.


title: An $\epsilon$-regularity criterion and estimates of the regular set for Navier-Stokes flows in terms of initial data

abstract: We prove an $\epsilon$-regularity criterion for the 3D Navier-Stokes equations in terms of initial data. It shows that if a scaled local $L^2$ norm of initial data is sufficiently small around the origin, a suitable weak solution is regular in a set enclosed by a paraboloid started from the origin. The result is applied to the estimate of the regular set for local energy solutions with initial data in weighted $L^2$ spaces. We also apply this result to studying energy concentration near a possible blow-up time and regularity of forward discretely self-similar solutions.


title: Increasing Subsequences and Kronecker Coefficients

abstract: It has been conjectured by W. Chen that the distribution of the length of the longest increasing subsequence in a uniformly random permutation is log-concave. We propose a stronger version of this conjecture which involves the Kronecker coefficients of the symmetric group.


title: Note on arXiv:2005.05301, 'Preparation of the Neutrino-4 experiment on search for sterile neutrino and the obtained results of measurements'

abstract: We comment on the claimed observation [arXiv:arXiv:2005.05301] of sterile neutrino oscillations by the Neutrino-4 collaboration. Such a claim, which requires the existence of a new fundamental particle, demands a level of rigor commensurate with its impact. The burden lies with the Neutrino-4 collaboration to provide the information necessary to prove the validity of their claim to the community. In this note, we describe aspects of both the data and analysis method that might lead to an oscillation signature arising from a null experiment and describe additional information needed from the Neutrino-4 collaboration to support the oscillation claim. Additionally, as opposed to the assertion made by the Neutrino-4 collaboration, we also show that the method of 'coherent summation' using the $L/E$ parameter produces similar results to the methods used by the PROSPECT and the STEREO collaborations.


title: Contact model for elastically anisotropic bodies and efficient implementation into the discrete-element method

abstract: We introduce a contact law for the normal force generated between two contacting, elastically anisotropic bodies of arbitrary geometry. The only requirement is that their surfaces be smooth and frictionless. This anisotropic contact law is obtained from a simplification of the exact solution to the continuum elasticity problem and takes the familiar form of Hertz' contact law, with the only difference being the orientation-dependence of the material-specific contact modulus. The contact law is remarkably accurate when compared with the exact solution, for a wide range of materials and surface geometries. We describe a computationally efficient implementation of the contact law into a discrete element method code, taking advantage of the precomputation of the contact modulus over all possible orientations. Finally, we showcase two application examples based on real materials where elastic anisotropy of the particles induces noticeable effects on macroscopic behavior. Notably, the second example demonstrates the ability to engineer tunable vibrational band gaps in a one-dimensional granular crystal by mere rotation of the constituent spherical particles.


title: Near Conformal Perturbation Theory in SYK Type Models

abstract: We present a systematic procedure to extract the dynamics of the low energy soft mode in SYK type models with a single energy scale $J$ and emergent reparametrization symmetry in the IR. This is given in the framework of the perturbation theory scheme of arXiv:1608.07567 based on specific (off-shell) breaking of conformal invariance in the UV, adjusted to yield the exact large-$N$ saddle point. While this breaking formally vanishes on-shell, it has a non-trivial effect on correlation functions and the effective action. In particular, it leads to the Schwarzian action with a specific coupling to bi-local matter. The method is applied to the evaluation of $O(1)$ corrections to the correlation function of bi-locals. As a byproduct we confirm precise agreement with the explicit, symmetry breaking procedure. We provide a verification in the large $q$ limit (Liouville theory), where the correlators can be calculated exactly at all length scales. In this case, our scheme illuminates how the enhanced $O(J)$ and the subleading $O(1)$ contributions originate from the Schwarzian dynamics of the soft mode and its interaction with $h=2$ (bi-local) matter.


title: Thickening of the diagonal and interleaving distance

abstract: Given a topological space $X$, a thickening kernel is a monoidal presheaf on $(\mathbb{R}_{\geq0},+)$ with values in the monoidal category of derived kernels on $X$. A bi-thickening kernel is defined on $(\mathbb{R},+)$. To such a thickening kernel, one naturally associates an interleaving distance on the derived category of sheaves on $X$. We prove that a thickening kernel exists and is unique as soon as it is defined on an interval containing $0$, allowing us to construct (bi-)thickenings in two different situations. First, when $X$ is a ``good'' metric space, starting with small usual thickenings of the diagonal. The associated interleaving distance satisfies the stability property and Lipschitz kernels give rise to Lipschitz maps. Second, by using [GKS12], when $X$ is a manifold and one is given a non-positive Hamiltonian isotopy on the cotangent bundle. In case $X$ is a complete Riemannian manifold having a strictly positive convexity radius, we prove that it is a good metric space and that the two bi-thickening kernels of the diagonal, one associated with the distance, the other with the geodesic flow, coincide.


title: Entanglement of pseudo-Hermitian random states

abstract: In a recent paper (arXiv:1905.07348v1), Dyson scheme to deal with density matrix of non-Hermitian Hamiltonians has been used to investigate the entanglement of states of a PT-symmetric bosonic system. They found that von Neumann entropy can show a different behavior in the broken and unbroken regime. We show that their results can be recast in terms of an abstract model of pseudo-Hermitian random matrices. It is found however that, although the formalism is practically the same, the entanglement is not of Fock states but of Bell states.


title: Magnify Your Population: Statistical Downscaling to Augment the Spatial Resolution of Socioeconomic Census Data

abstract: Fine resolution estimates of demographic and socioeconomic attributes are crucial for planning and policy development. While several efforts have been made to produce fine-scale gridded population estimates, socioeconomic features are typically not available at scales finer than Census units, which may hide local heterogeneity and disparity. In this paper we present a new statistical downscaling approach to derive fine-scale estimates of key socioeconomic attributes. The method leverages demographic and geographical extensive covariates available at multiple scales and additional Census covariates only available at coarse resolution, which are included in the model hierarchically within a "forward learning" approach. For each selected socioeconomic variable, a Random Forest model is trained on the source Census units and then used to generate fine-scale gridded predictions, which are then adjusted to ensure the best possible consistency with the coarser Census data. As a case study, we apply this method to Census data in the United States, downscaling the selected socioeconomic variables available at the block group level, to a grid of ~300 spatial resolution. The accuracy of the method is assessed at both spatial scales, first computing a pseudo cross-validation coefficient of determination for the predictions at the block group level and then, for extensive variables only, also for the (unadjusted) predicted counts summed by block group. Based on these scores and on the inspection of the downscaled maps, we conclude that our method is able to provide accurate, smoother, and more detailed socioeconomic estimates than the available Census data.


title: Learning dynamics for improving control of overactuated flying systems

abstract: Overactuated omnidirectional flying vehicles are capable of generating force and torque in any direction, which is important for applications such as contact-based industrial inspection. This comes at the price of an increase in model complexity. These vehicles usually have non-negligible, repetitive dynamics that are hard to model, such as the aerodynamic interference between the propellers. This makes it difficult for high-performance trajectory tracking using a model-based controller. This paper presents an approach that combines a data-driven and a first-principle model for the system actuation and uses it to improve the controller. In a first step, the first-principle model errors are learned offline using a Gaussian Process (GP) regressor. At runtime, the first-principle model and the GP regressor are used jointly to obtain control commands. This is formulated as an optimization problem, which avoids ambiguous solutions present in a standard inverse model in overactuated systems, by only using forward models. The approach is validated using a tilt-arm overactuated omnidirectional flying vehicle performing attitude trajectory tracking. The results show that with our proposed method, the attitude trajectory error is reduced by 32% on average as compared to a nominal PID controller.


title: Inferring Causal Networks of Dynamical Systems through Transient Dynamics and Perturbation

abstract: Inferring causal relations from time series measurements is an ill-posed mathematical problem, where typically an infinite number of potential solutions can reproduce the given data. We explore in depth a strategy to disambiguate between possible underlying causal networks by perturbing the network, where the actuations are either targeted or applied at random. The resulting transient dynamics provide the critical information necessary to infer causality. Two methods are shown to provide accurate causal reconstructions: Granger causality (GC) with perturbations, and our proposed perturbation cascade inference (PCI). Perturbed GC is capable of inferring smaller networks under low coupling strength regimes. Our proposed PCI method demonstrated consistently strong performance in inferring causal relations for small (2-5 node) and large (10-20 node) networks, with both linear and nonlinear dynamics. Thus the ability to apply a large and diverse set of perturbations/actuations to the network is critical for successfully and accurately determining causal relations and disambiguating between various viable networks.


title: Logical Neural Networks

abstract: We propose a novel framework seamlessly providing key properties of both neural nets (learning) and symbolic logic (knowledge and reasoning). Every neuron has a meaning as a component of a formula in a weighted real-valued logic, yielding a highly intepretable disentangled representation. Inference is omnidirectional rather than focused on predefined target variables, and corresponds to logical reasoning, including classical first-order logic theorem proving as a special case. The model is end-to-end differentiable, and learning minimizes a novel loss function capturing logical contradiction, yielding resilience to inconsistent knowledge. It also enables the open-world assumption by maintaining bounds on truth values which can have probabilistic semantics, yielding resilience to incomplete knowledge.


title: Lithium in the context of Multiple Populations in Globular Clusters

abstract: Multiple Populations represent the standard for Globular Clusters (GC): a fraction (10-50%) of their stars have the same elemental abundances of halo stars of similar metallicity (first generation, or 1G), but the other stars (second generation, 2G) are characterised by patterns of light elements abundances which resemble those typical of gas processed by proton-capture reactions at high temperature. Consequently, we should naively expect that Lithium is destroyed in the 2G stars, but instead it is generally observed, at abundances only slightly depleted with respect to the 1G stars. After discussing the models for the formation of multiple populations, I examine the role of dilution with pristine gas and the possible role of the Asymptotic Giant Branch (AGB) scenario in accounting for the Lithium patterns in GC stars. Super-AGB and AGB yields of Lithium, produced by the Cameron Fowler mechanism in the `Hot Bottom Burning' convective envelopes, may help to explain the peculiar high Lithium in a few extreme 2G stars. On the other hand, modeling the abundances in mild 2G stars depends explicitly on whether the initial Li in the gas forming the 1G stars and the diluting gas of the 2G ones is that predicted by the Big Bang nucleosynthesis or the ~3 times smaller value observed at the surface of halo dwarfs.


title: Fractional statistics in anyon collisions

abstract: Two-dimensional systems can host exotic particles called anyons whose quantum statistics are neither bosonic nor fermionic. For example, the elementary excitations of the fractional quantum Hall effect at filling factor $\nu=1/m$ (where m is an odd integer) have been predicted to obey abelian fractional statistics, with a phase $\varphi$ associated with the exchange of two particles equal to $\pi/m$. However, despite numerous experimental attempts, clear signatures of fractional statistics remain elusive. Here we experimentally demonstrate abelian fractional statistics at filling factor $\nu=1/3$ by measuring the current correlations resulting from the collision between anyons at a beam-splitter. By analyzing their dependence on the anyon current impinging on the splitter and comparing with recent theoretical models, we extract $\varphi=\pi/3$, in agreement with predictions.


title: The principles of adaptation in organisms and machines II: Thermodynamics of the Bayesian brain

abstract: This article reviews how organisms learn and recognize the world through the dynamics of neural networks from the perspective of Bayesian inference, and introduces a view on how such dynamics is described by the laws for the entropy of neural activity, a paradigm that we call thermodynamics of the Bayesian brain. The Bayesian brain hypothesis sees the stimulus-evoked activity of neurons as an act of constructing the Bayesian posterior distribution based on the generative model of the external world that an organism possesses. A closer look at the stimulus-evoked activity at early sensory cortices reveals that feedforward connections initially mediate the stimulus-response, which is later modulated by input from recurrent connections. Importantly, not the initial response, but the delayed modulation expresses animals' cognitive states such as awareness and attention regarding the stimulus. Using a simple generative model made of a spiking neural population, we reproduce the stimulus-evoked dynamics with the delayed feedback modulation as the process of the Bayesian inference that integrates the stimulus evidence and a prior knowledge with time-delay. We then introduce a thermodynamic view on this process based on the laws for the entropy of neural activity. This view elucidates that the process of the Bayesian inference works as the recently-proposed information-theoretic engine (neural engine, an analogue of a heat engine in thermodynamics), which allows us to quantify the perceptual capacity expressed in the delayed modulation in terms of entropy.


title: Gravitational Production of Hidden Photon Dark Matter in Light of the XENON1T Excess

abstract: Recently, the XENON1T experiment has reported an excess in the electronic recoil events. The excess is consistent with the interpretation of absorption of 3keV bosonic dark matter, for example, hidden photon dark matter with kinetic mixing of the order of $10^{-15}$. We point out that the minimally gravitational production provides a viable mechanism for obtaining a correct relic hidden photon abundance. We present parameter dependence of the hidden photon dark matter abundance on the inflationary scale $H_{\rm inf}$ and also the reheating temperature $T_{\rm R}$. We show that the inflationary Hubble scale and reheating temperature are both bounded from below, $H_{\rm inf}\gtrsim 7\times 10^{11}$GeV, $T_{\rm R}\gtrsim 10^2$GeV. In particular, the high-scale inflation is consistent with 3keV hidden photon dark matter.


title: Environment Shaping in Reinforcement Learning using State Abstraction

abstract: One of the central challenges faced by a reinforcement learning (RL) agent is to effectively learn a (near-)optimal policy in environments with large state spaces having sparse and noisy feedback signals. In real-world applications, an expert with additional domain knowledge can help in speeding up the learning process via \emph{shaping the environment}, i.e., making the environment more learner-friendly. A popular paradigm in literature is \emph{potential-based reward shaping}, where the environment's reward function is augmented with additional local rewards using a potential function. However, the applicability of potential-based reward shaping is limited in settings where (i) the state space is very large, and it is challenging to compute an appropriate potential function, (ii) the feedback signals are noisy, and even with shaped rewards the agent could be trapped in local optima, and (iii) changing the rewards alone is not sufficient, and effective shaping requires changing the dynamics. We address these limitations of potential-based shaping methods and propose a novel framework of \emph{environment shaping using state abstraction}. Our key idea is to compress the environment's large state space with noisy signals to an abstracted space, and to use this abstraction in creating smoother and more effective feedback signals for the agent. We study the theoretical underpinnings of our abstraction-based environment shaping, and show that the agent's policy learnt in the shaped environment preserves near-optimal behavior in the original environment.


title: Collider Constraints on a Dark Matter Interpretation of the XENON1T Excess

abstract: In light of the excess in the low-energy electron recoil events reported by XENON1T, many new physics scenarios have been proposed as a possible origin of the excess. One possible explanation is that the excess is a result of a fast moving dark matter (DM), with velocity $v\sim0.05-0.20$ and mass between 1 MeV and 10 GeV, scattering off an electron. Assuming the fast moving DM-electron interaction is mediated by a vector particle, we derive collider constraints on the said DM-electron interaction. The bounds on DM-electron coupling is then used to constrain possible production mechanisms of the fast moving DM. We find that the preferred mass of the vector mediator is relatively light ($\lesssim$ 1 GeV) and the coupling of the vector to the electron is much smaller than the coupling to the fast moving DM.


title: Discrete correlations of order 2 of generalised Rudin-Shapiro sequences: a combinatorial approach

abstract: We introduce a family of block-additive automatic sequences, that are obtained by allocating a weight to each couple of digits, and defining the $n$th term of the sequence as being the total weight of the integer $n$ written in base $k$. Under an additional difference condition on the weight function, these sequences can be interpreted as generalised Rudin-Shapiro sequences, and we prove that they have the same correlations of order 2 as sequences of symbols chosen uniformly and independently at random. The speed of convergence is very fast and is independent of the prime factor decomposition of $k$. This extends recent work of Tahay. The proof relies on direct observations about base-$k$ representations of integers and combinatorial considerations. We also provide extensions of our results to higher-dimensional block-additive sequences.


title: MANTRA: A Machine Learning reference lightcurve dataset for astronomical transient event recognition

abstract: We introduce MANTRA, an annotated dataset of 4869 transient and 71207 non-transient object lightcurves built from the Catalina Real Time Transient Survey. We provide public access to this dataset as a plain text file to facilitate standardized quantitative comparison of astronomical transient event recognition algorithms. Some of the classes included in the dataset are: supernovae, cataclysmic variables, active galactic nuclei, high proper motion stars, blazars and flares. As an example of the tasks that can be performed on the dataset we experiment with multiple data pre-processing methods, feature selection techniques and popular machine learning algorithms (Support Vector Machines, Random Forests and Neural Networks). We assess quantitative performance in two classification tasks: binary (transient/non-transient) and eight-class classification. The best performing algorithm in both tasks is the Random Forest Classifier. It achieves an F1-score of 96.25% in the binary classification and 52.79% in the eight-class classification. For the eight-class classification, non-transients ( 96.83% ) is the class with the highest F1-score, while the lowest corresponds to high-proper-motion stars ( 16.79% ); for supernovae it achieves a value of 54.57% , close to the average across classes. The next release of MANTRA includes images and benchmarks with deep learning models.


title: Joint Object Detection and Multi-Object Tracking with Graph Neural Networks

abstract: Object detection and data association are critical components in multi-object tracking (MOT) systems. Despite the fact that the two components are dependent on each other, prior works often design detection and data association modules separately which are trained with separate objectives. As a result, one cannot back-propagate the gradients and optimize the entire MOT system, which leads to sub-optimal performance. To address this issue, recent works simultaneously optimize detection and data association modules under a joint MOT framework, which has shown improved performance in both modules. In this work, we propose a new instance of joint MOT approach based on Graph Neural Networks (GNNs). The key idea is that GNNs can model relations between variable-sized objects in both the spatial and temporal domains, which is essential for learning discriminative features for detection and data association. Through extensive experiments on the MOT15/16/17/20 datasets, we demonstrate the effectiveness of our GNN-based joint MOT approach and show state-of-the-art performance for both detection and MOT tasks. Our code is available at: https://github.com/yongxinw/GSDT


title: Provably Efficient Reinforcement Learning for Discounted MDPs with Feature Mapping

abstract: Modern tasks in reinforcement learning have large state and action spaces. To deal with them efficiently, one often uses predefined feature mapping to represent states and actions in a low-dimensional space. In this paper, we study reinforcement learning for discounted Markov Decision Processes (MDPs), where the transition kernel can be parameterized as a linear function of certain feature mapping. We propose a novel algorithm that makes use of the feature mapping and obtains a $\tilde O(d\sqrt{T}/(1-\gamma)^2)$ regret, where $d$ is the dimension of the feature space, $T$ is the time horizon and $\gamma$ is the discount factor of the MDP. To the best of our knowledge, this is the first polynomial regret bound without accessing the generative model or making strong assumptions such as ergodicity of the MDP. By constructing a special class of MDPs, we also show that for any algorithms, the regret is lower bounded by $\Omega(d\sqrt{T}/(1-\gamma)^{1.5})$. Our upper and lower bound results together suggest that the proposed reinforcement learning algorithm is near-optimal up to a $(1-\gamma)^{-0.5}$ factor.


title: Steiner's Hat: a Constant-Area Deltoid Associated with the Ellipse

abstract: The Negative Pedal Curve (NPC) of the Ellipse with respect to a boundary point M is a 3-cusp closed-curve which is the affine image of the Steiner Deltoid. Over all M the family has invariant area and displays an array of interesting properties.


title: Diffusions interacting through a random matrix: universality via stochastic Taylor expansion

abstract: Consider $(X_{i}(t))$ solving a system of $N$ stochastic differential equations interacting through a random matrix $\mathbf J = (J_{ij})$ with independent (not necessarily identically distributed) random coefficients. We show that the trajectories of averaged observables of $(X_i(t))$, initialized from some $\mu$ independent of $\mathbf J$, are universal, i.e., only depend on the choice of the distribution $\mathbf{J}$ through its first and second moments (assuming e.g., sub-exponential tails). We take a general combinatorial approach to proving universality for dynamical systems with random coefficients, combining a stochastic Taylor expansion with a moment matching-type argument. Concrete settings for which our results imply universality include aging in the spherical SK spin glass, and Langevin dynamics and gradient flows for symmetric and asymmetric Hopfield networks.


title: Lattice dynamics in the conformational environment of multilayered hexagonal boron nitride (h-BN) conveys to peculiar infrared optical responses

abstract: Stacking mismatches in hexagonal boron nitride (h-BN) nanostructures affect their photonic, mechanical, and thermal properties. To access information about the stacked configuration of layered ensembles, highly sophisticated techniques like X-ray photoemission spectroscopy or electron microscopy are necessary. Here, instead, by taking advantage of the geometrical and chemical nature of h-BN, we show how simple structural models, based on shortened interplanar distances, can produce effective charge densities. Accounting these in the non-analytical part of the lattice dynamical description makes it possible to access information about the composition of differently stacked variants in experimental samples characterized by infrared spectroscopy. The results are obtained by density functional theory and confirmed by various functionals and pseudopotential approximations. Even though the method is shown on h-BN, the conclusions are more general and show how effective dielectric models can be considered as valuable theoretical pathways for the vibrational structure of any layered material.

categories: cond-mat.mtrl-sci cond-mat.mes-hall physics.app-ph physics.chem-ph physics.comp-ph


title: Experience Replay with Likelihood-free Importance Weights

abstract: The use of past experiences to accelerate temporal difference (TD) learning of value functions, or experience replay, is a key component in deep reinforcement learning. Prioritization or reweighting of important experiences has shown to improve performance of TD learning algorithms.In this work, we propose to reweight experiences based on their likelihood under the stationary distribution of the current policy. Using the corresponding reweighted TD objective, we implicitly encourage small approximation errors on the value function over frequently encountered states. We use a likelihood-free density ratio estimator over the replay buffer to assign the prioritization weights. We apply the proposed approach empirically on two competitive methods, Soft Actor Critic (SAC) and Twin Delayed Deep Deterministic policy gradient (TD3) -- over a suite of OpenAI gym tasks and achieve superior sample complexity compared to other baseline approaches.


title: Variational Orthogonal Features

abstract: Sparse stochastic variational inference allows Gaussian process models to be applied to large datasets. The per iteration computational cost of inference with this method is $\mathcal{O}(\tilde{N}M^2+M^3),$ where $\tilde{N}$ is the number of points in a minibatch and $M$ is the number of `inducing features', which determine the expressiveness of the variational family. Several recent works have shown that for certain priors, features can be defined that remove the $\mathcal{O}(M^3)$ cost of computing a minibatch estimate of an evidence lower bound (ELBO). This represents a significant computational savings when $M\gg \tilde{N}$. We present a construction of features for any stationary prior kernel that allow for computation of an unbiased estimator to the ELBO using $T$ Monte Carlo samples in $\mathcal{O}(\tilde{N}T+M^2T)$ and in $\mathcal{O}(\tilde{N}T+MT)$ with an additional approximation. We analyze the impact of this additional approximation on inference quality.


title: ObjectNav Revisited: On Evaluation of Embodied Agents Navigating to Objects

abstract: We revisit the problem of Object-Goal Navigation (ObjectNav). In its simplest form, ObjectNav is defined as the task of navigating to an object, specified by its label, in an unexplored environment. In particular, the agent is initialized at a random location and pose in an environment and asked to find an instance of an object category, e.g., find a chair, by navigating to it. As the community begins to show increased interest in semantic goal specification for navigation tasks, a number of different often-inconsistent interpretations of this task are emerging. This document summarizes the consensus recommendations of this working group on ObjectNav. In particular, we make recommendations on subtle but important details of evaluation criteria (for measuring success when navigating towards a target object), the agent's embodiment parameters, and the characteristics of the environments within which the task is carried out. Finally, we provide a detailed description of the instantiation of these recommendations in challenges organized at the Embodied AI workshop at CVPR 2020 http://embodied-ai.org .


title: Multiple populations in massive star clusters under the magnifying glass of photometry: Theory and tools

abstract: The existence of star-to-star light-element abundance variations in massive Galactic and extragalactic star clusters has fairly recently superseded the traditional paradigm of individual clusters hosting stars with the same age, and uniform chemical composition. Several scenarios have been put forward to explain the origin of this multiple stellar population phenomenon, but so far all have failed to reproduce the whole range of key observations. Complementary to high-resolution spectroscopy, which has first revealed and characterized chemically the presence of multiple populations in Galactic globular clusters, photometry has been instrumental in investigating this phenomenon in much larger samples of stars --adding a number of crucial observational constraints and correlations with global cluster properties-- and in the discovery and characterization of multiple populations also in Magellanic Clouds' intermediate age clusters. The purpose of this review is to present the theoretical underpinning and application of the photometric techniques devised to identify and study multiple populations in resolved star clusters. These methods have played and continue to play a crucial role in advancing our knowledge of the cluster multiple population phenomenon, and promise to extend the scope of these investigations to resolved clusters even beyond the Local Group, with the launch of the James Webb Space Telescope.


title: Deep Reinforcement Learning Control for Radar Detection and Tracking in Congested Spectral Environments

abstract: In this paper, dynamic non-cooperative coexistence between a cognitive pulsed radar and a nearby communications system is addressed by applying nonlinear value function approximation via deep reinforcement learning (Deep RL) to develop a policy for optimal radar performance. The radar learns to vary the bandwidth and center frequency of its linear frequency modulated (LFM) waveforms to mitigate mutual interference with other systems and improve target detection performance while also maintaining sufficient utilization of the available frequency bands required for a fine range resolution. We demonstrate that our approach, based on the Deep Q-Learning (DQL) algorithm, enhances important radar metrics, including SINR and bandwidth utilization, more effectively than policy iteration or sense-and-avoid (SAA) approaches in a variety of realistic coexistence environments. We also extend the DQL-based approach to incorporate Double Q-learning and a recurrent neural network to form a Double Deep Recurrent Q-Network (DDRQN). We demonstrate the DDRQN results in favorable performance and stability compared to DQL and policy iteration. Finally, we demonstrate the practicality of our proposed approach through a discussion of experiments performed on a software defined radar (SDRadar) prototype system. Our experimental results indicate that the proposed Deep RL approach significantly improves radar detection performance in congested spectral environments when compared to policy iteration and SAA.


title: Partial regularity of a nematic liquid crystal model with kinematic transport effects

abstract: In this paper, we will establish the global existence of a suitable weak solution to the Erickson--Leslie system modeling hydrodynamics of nematic liquid crystal flows with kinematic transports for molecules of various shapes in ${\mathbb{R}^3}$, which is smooth away from a closed set of (parabolic) Hausdorff dimension at most $\frac{15}{7}$.


title: Forward trijet production in p-p and p-Pb collisions at LHC

abstract: We calculate various azimuthal angle distributions for three jets produced in the forward rapidity region with transverse momenta $p_T>20\,\mathrm{GeV}$ in proton-proton (p-p) and proton-lead (p-Pb) collisions at center of mass energy $5.02\,\,\mathrm{TeV}$. We use the multi-parton extension of the so-called small-$x$ Improved Transverse Momentum Dependent factorization (ITMD). We study effects related to change from the standard $k_T$-factorization to ITMD factorization as well as changes as one goes from p-p collision to p-Pb. We observe rather large differences in the distribution when we change the factorization approach, which allows to both improve the small-$x$ TMD gluon distributions as well as validate and improve the factorization approach. We also see significant depletion of the nuclear modification ratio, indicating a possibility of searches for saturation effects using trijet final states in a more exclusive way than for dijets.


title: Boundary Regularized Building Footprint Extraction From Satellite Images Using Deep Neural Network

abstract: In recent years, an ever-increasing number of remote satellites are orbiting the Earth which streams vast amount of visual data to support a wide range of civil, public and military applications. One of the key information obtained from satellite imagery is to produce and update spatial maps of built environment due to its wide coverage with high resolution data. However, reconstructing spatial maps from satellite imagery is not a trivial vision task as it requires reconstructing a scene or object with high-level representation such as primitives. For the last decade, significant advancement in object detection and representation using visual data has been achieved, but the primitive-based object representation still remains as a challenging vision task. Thus, a high-quality spatial map is mainly produced through complex labour-intensive processes. In this paper, we propose a novel deep neural network, which enables to jointly detect building instance and regularize noisy building boundary shapes from a single satellite imagery. The proposed deep learning method consists of a two-stage object detection network to produce region of interest (RoI) features and a building boundary extraction network using graph models to learn geometric information of the polygon shapes. Extensive experiments show that our model can accomplish multi-tasks of object localization, recognition, semantic labelling and geometric shape extraction simultaneously. In terms of building extraction accuracy, computation efficiency and boundary regularization performance, our model outperforms the state-of-the-art baseline models.


title: Inference with Artificial Neural Networks on Analog Neuromorphic Hardware

abstract: The neuromorphic BrainScaleS-2 ASIC comprises mixed-signal neurons and synapse circuits as well as two versatile digital microprocessors. Primarily designed to emulate spiking neural networks, the system can also operate in a vector-matrix multiplication and accumulation mode for artificial neural networks. Analog multiplication is carried out in the synapse circuits, while the results are accumulated on the neurons' membrane capacitors. Designed as an analog, in-memory computing device, it promises high energy efficiency. Fixed-pattern noise and trial-to-trial variations, however, require the implemented networks to cope with a certain level of perturbations. Further limitations are imposed by the digital resolution of the input values (5 bit), matrix weights (6 bit) and resulting neuron activations (8 bit). In this paper, we discuss BrainScaleS-2 as an analog inference accelerator and present calibration as well as optimization strategies, highlighting the advantages of training with hardware in the loop. Among other benchmarks, we classify the MNIST handwritten digits dataset using a two-dimensional convolution and two dense layers. We reach 98.0% test accuracy, closely matching the performance of the same network evaluated in software.


title: Does Matter Matter? Using the mass distribution to distinguish neutron stars and black holes

abstract: Gravitational-wave detectors have opened a new window through which we can observe black holes (BHs) and neutron stars (NSs). Analyzing the 11 detections from LIGO/Virgo's first gravitational-wave catalog, GWTC-1, we investigate whether the power-law fit to the BH mass spectrum can also accommodate the binary neutron star (BNS) event GW170817, or whether we require an additional feature, such as a mass gap, in between the NS and BH populations. We find that with respect to the power-law fit to binary black hole (BBH) masses, GW170817 is an outlier at the 0.13\% level, suggesting a distinction between NS and BH masses. A single power-law fit across the entire mass range is in mild tension with: (a) the detection of one source in the BNS mass range ($\sim 1$--$2.5 \,M_\odot$), (b) the absence of detections in the "mass-gap" range ($\sim 2.5$--$5 \,M_\odot$), and (c) the detection of 10 sources in the BBH mass range ($\gtrsim 5 \,M_\odot$). Instead, the data favor models with a feature between NS and BH masses, including a mass gap (Bayes factor of 4.6) and a break in the power law, with a steeper slope at NS masses compared to BH masses (91\% credibility). We estimate the merger rates of compact binaries based on our fit to the global mass distribution, finding $\mathcal{R}_\mathrm{BNS} = 871^{+3015}_{-805} \ \mathrm{Gpc}^{-3} \ \mathrm{yr}^{-1}$ and $\mathcal{R}_\mathrm{BBH} = 47.5^{+57.9}_{-28.8} \ \mathrm{Gpc}^{-3} \ \mathrm{yr}^{-1}$. We conclude that, even in the absence of any prior knowledge of the difference between NSs and BHs, the gravitational-wave data alone already suggest two distinct populations of compact objects.


title: $\alpha$-event Characterization and Rejection in Point-Contact HPGe Detectors

abstract: P-type point contact (PPC) HPGe detectors are a leading technology for rare event searches due to their excellent energy resolution, low thresholds, and multi-site event rejection capabilities. We have characterized a PPC detector's response to $\alpha$ particles incident on the sensitive passivated and p+ surfaces, a previously poorly-understood source of background. The detector studied is identical to those in the MAJORANA DEMONSTRATOR experiment, a search for neutrinoless double-beta decay ($0\nu\beta\beta$) in $^{76}$Ge. $\alpha$ decays on most of the passivated surface exhibit significant energy loss due to charge trapping, with waveforms exhibiting a delayed charge recovery (DCR) signature caused by the slow collection of a fraction of the trapped charge. The DCR is found to be complementary to existing methods of $\alpha$ identification, reliably identifying $\alpha$ background events on the passivated surface of the detector. We demonstrate effective rejection of all surface $\alpha$ events (to within statistical uncertainty) with a loss of only 0.2% of bulk events by combining the DCR discriminator with previously-used methods. The DCR discriminator has been used to reduce the background rate in the $0\nu\beta\beta$ region of interest window by an order of magnitude in the MAJORANA DEMONSTRATOR, and will be used in the upcoming LEGEND-200 experiment.


title: New Perspectives on the so-called Fermi's Golden Rule in Quantum Mechanics including Adiabatic Following

abstract: A novel and readily understandable derivation of the Golden Rule of time dependent perturbation theory is presented. The derivation is based on adiabatic turning on of the perturbation as used, for instance, in some formal developments of scattering theory. Energy conservation is expressed in terms of an intuitively and physically appealing Lorentzian line shape rather than the artificial, oscillatory sin(x)=x type line shape that appears in conventional derivations. The conditions for the derivation's validity are compactly and conveniently expressed in the frequency/energy domain rather than in the usual time domain. The derivation also highlights how, along with energy conservation, the transition rate faithfully and instantaneously follows the variations in the square of the perturbing potential as one may expect in the adiabatic limit. In the first instance, the adiabatic turning on is achieved, as usual, by a single exponential time variation. But we demonstrate that the instantaneous following of the square of the perturbing potential by the transition rate is more general and that one can derive the Golden Rule for a general slowly varying time dependent perturbation. This allows one to derive generalisations of the simple decay law, originally derived in the classic paper by Weisskopf and Wigner; a tutorial exposition of the essence of this classic work is provided. The Oppenheimer method for applying the Golden rule to problems, such as the electric field ionisation of atoms, in which the perturbing potential can also create the final states, is reviewed. No use of an energy gap condition is needed to derive our results on adiabatic behaviour in contrast to the origianl derivations of the adiabatic theorem in quantum mechanics.


title: Numerical aspects of integration in semi-closed option pricing formulas for stochastic volatility jump diffusion models

abstract: In mathematical finance, a process of calibrating stochastic volatility (SV) option pricing models to real market data involves a numerical calculation of integrals that depend on several model parameters. This optimization task consists of large number of integral evaluations with high precision and low computational time requirements. However, for some model parameters, many numerical quadrature algorithms fail to meet these requirements. We can observe an enormous increase in function evaluations, serious precision problems and a significant increase of computational time. In this paper we numerically analyse these problems and show that they are especially caused by inaccurately evaluated integrands. We propose a fast regime switching algorithm that tells if it is sufficient to evaluate the integrand in standard double arithmetic or if a higher precision arithmetic has to be used. We compare and recommend numerical quadratures for typical SV models and different parameter values, especially for problematic cases.


title: On the Global Optimality of Model-Agnostic Meta-Learning

abstract: Model-agnostic meta-learning (MAML) formulates meta-learning as a bilevel optimization problem, where the inner level solves each subtask based on a shared prior, while the outer level searches for the optimal shared prior by optimizing its aggregated performance over all the subtasks. Despite its empirical success, MAML remains less understood in theory, especially in terms of its global optimality, due to the nonconvexity of the meta-objective (the outer-level objective). To bridge such a gap between theory and practice, we characterize the optimality gap of the stationary points attained by MAML for both reinforcement learning and supervised learning, where the inner-level and outer-level problems are solved via first-order optimization methods. In particular, our characterization connects the optimality gap of such stationary points with (i) the functional geometry of inner-level objectives and (ii) the representation power of function approximators, including linear models and neural networks. To the best of our knowledge, our analysis establishes the global optimality of MAML with nonconvex meta-objectives for the first time.


title: Exothermic Dark Matter for XENON1T Excess

abstract: Motivated by the recent excess in the electron recoil from XENON1T experiment, we consider the possibility of exothermic dark matter, which is composed of two states with mass splitting. The heavier state down-scatters off the electron into the lighter state, making an appropriate recoil energy required for the Xenon excess even for the standard Maxwellian velocity distribution of dark matter. Accordingly, we determine the mass difference between two component states of dark matter to the peak electron recoil energy at about $2.5\,{\rm keV}$ up to the detector resolution, accounting for the recoil events over $E_R=2-3\,{\rm keV}$, which are most significant. We include the effects of the phase-space enhancement and the atomic excitation factor to calculate the required scattering cross section for the Xenon excess. We discuss the implications of dark matter interactions in the effective theory for exothermic dark matter and a massive $Z'$ mediator and provide microscopic models realizing the required dark matter and electron couplings to $Z'$.


title: Scintillation can explain the spectral structure of the bright radio burst from SGR 1935+2154

abstract: The discovery of a fast radio burst (FRB) associated with a magnetar in the Milky Way by the Canadian Hydrogen Intensity Mapping Experiment FRB collaboration (CHIME/FRB) and the Survey for Transient Astronomical Radio Emission 2 (STARE2) has provided an unprecedented opportunity to refine FRB emission models. The burst discovered by CHIME/FRB shows two components with different spectra. We explore interstellar scintillation as the origin for this variation in spectral structure. Modeling a weak scattering screen in the supernova remnant associated with the magnetar, we find that a superluminal apparent transverse velocity of the emission region of $>9.5\,c$ is needed to explain the spectral variation. Alternatively, the two components could have originated from independent emission regions spaced by $>8.3\times10^4\,$km. These scenarios may arise in "far-away" models where the emission originates from well beyond the magnetosphere of the magnetar (for example through a synchrotron-maser mechanism set up by an ultra-relativistic radiative shock), but not in "close-in" models of emission from within the magnetosphere. If further radio observations of the magnetar confirm scintillation as the source for the observed variation in spectral structure, this scattering model thus constrains the location of the emission region.


title: Luminescence Enhancement in One-Dimensional Mie-Resonant Arrays

abstract: In this paper, we demonstrate the infrared photoluminescence emission from Ge(Si) quantum dots enhanced with collective Mie modes of silicon nanopillars. We show that the excitation of band edge dipolar modes of a linear nanopillar array results in strong reshaping of the photoluminescence spectra. Among other collective modes, the magnetic dipolar mode with the polarization along the array axis contributes the most to the emission spectrum exhibiting an experimentally measured Q-factor of around 500 for an array of 11 pillars. The results belong to the first experimental evidence of light emission enhancement of quantum emitters applying collective Mie resonances and therefore represent an important contribution to the new field of active all-dielectric meta-optics.


title: Two-fluid simulations of the magnetic field evolution in neutron star cores in the weak-coupling regime

abstract: In a previous paper, we reported simulations of the evolution of the magnetic field in neutron star cores through ambipolar diffusion, taking the neutrons as a motionless uniform background. However, in real neutron stars, neutrons are free to move, and a strong composition gradient leads to stable stratification (stability against convective motions) both of which might impact on the time-scales of evolution. Here we address these issues by providing the first long-term two-fluid simulations of the evolution of an axially symmetric magnetic field in a neutron star core composed of neutrons, protons, and electrons with density and composition gradients. Again, we find that the magnetic field evolves towards barotropic "Grad-Shafranov equilibria", in which the magnetic force is balanced by the degeneracy pressure gradient and gravitational force of the charged particles. However, the evolution is found to be faster than in the case of motionless neutrons, as the movement of charged particles (which are coupled to the magnetic field, but are also limited by the collisional drag forces exerted by neutrons) is less constrained, since neutrons are now allowed to move. The possible impact of non-axisymmetric instabilities on these equilibria, as well as beta decays, proton superconductivity, and neutron superfluidity, are left for future work.


title: Network-based brain computer interfaces: principles and applications

abstract: Brain-computer interfaces (BCIs) make possible to interact with the external environment by decoding the mental intention of individuals. BCIs can therefore be used to address basic neuroscience questions but also to unlock a variety of applications from exoskeleton control to neurofeedback (NFB) rehabilitation. In general, BCI usability critically depends on the ability to comprehensively characterize brain functioning and correctly identify the user s mental state. To this end, much of the efforts have focused on improving the classification algorithms taking into account localized brain activities as input features. Despite considerable improvement BCI performance is still unstable and, as a matter of fact, current features represent oversimplified descriptors of brain functioning. In the last decade, growing evidence has shown that the brain works as a networked system composed of multiple specialized and spatially distributed areas that dynamically integrate information. While more complex, looking at how remote brain regions functionally interact represents a grounded alternative to better describe brain functioning. Thanks to recent advances in network science, i.e. a modern field that draws on graph theory, statistical mechanics, data mining and inferential modelling, scientists have now powerful means to characterize complex brain networks derived from neuroimaging data. Notably, summary features can be extracted from these networks to quantitatively measure specific organizational properties across a variety of topological scales. In this topical review, we aim to provide the state-of-the-art supporting the development of a network theoretic approach as a promising tool for understanding BCIs and improve usability.


title: Efficient Spatially Adaptive Convolution and Correlation

abstract: Fast methods for convolution and correlation underlie a variety of applications in computer vision and graphics, including efficient filtering, analysis, and simulation. However, standard convolution and correlation are inherently limited to fixed filters: spatial adaptation is impossible without sacrificing efficient computation. In early work, Freeman and Adelson have shown how steerable filters can address this limitation, providing a way for rotating the filter as it is passed over the signal. In this work, we provide a general, representation-theoretic, framework that allows for spatially varying linear transformations to be applied to the filter. This framework allows for efficient implementation of extended convolution and correlation for transformation groups such as rotation (in 2D and 3D) and scale, and provides a new interpretation for previous methods including steerable filters and the generalized Hough transform. We present applications to pattern matching, image feature description, vector field visualization, and adaptive image filtering.


title: Expert-Supervised Reinforcement Learning for Offline Policy Learning and Evaluation

abstract: Offline Reinforcement Learning (RL) is a promising approach for learning optimal policies in environments where direct exploration is expensive or unfeasible. However, the adoption of such policies in practice is often challenging, as they are hard to interpret within the application context, and lack measures of uncertainty for the learned policy value and its decisions. To overcome these issues, we propose an Expert-Supervised RL (ESRL) framework which uses uncertainty quantification for offline policy learning. In particular, we have three contributions: 1) the method can learn safe and optimal policies through hypothesis testing, 2) ESRL allows for different levels of risk averse implementations tailored to the application context, and finally, 3) we propose a way to interpret ESRL's policy at every state through posterior distributions, and use this framework to compute off-policy value function posteriors. We provide theoretical guarantees for our estimators and regret bounds consistent with Posterior Sampling for RL (PSRL). Sample efficiency of ESRL is independent of the chosen risk aversion threshold and quality of the behavior policy.


title: Facing the Hard Problems in FGVC

abstract: In fine-grained visual categorization (FGVC), there is a near-singular focus in pursuit of attaining state-of-the-art (SOTA) accuracy. This work carefully analyzes the performance of recent SOTA methods, quantitatively, but more importantly, qualitatively. We show that these models universally struggle with certain "hard" images, while also making complementary mistakes. We underscore the importance of such analysis, and demonstrate that combining complementary models can improve accuracy on the popular CUB-200 dataset by over 5%. In addition to detailed analysis and characterization of the errors made by these SOTA methods, we provide a clear set of recommended directions for future FGVC researchers.


title: FUGIN hot core survey. I. Survey method and initial results for $l = 10^\circ-20^\circ$

abstract: We have developed a method to make a spectral-line-based survey of hot cores, which represent an important stage of high-mass star formation, and applied the method to the data of the FUGIN (FOREST Unbiased Galactic plane Imaging survey with the Nobeyama 45-m telescope) survey. First, we select hot core candidates by searching the FUGIN data for the weak hot core tracer lines (HNCO and CH$_3$CN) by stacking, and then we conduct follow-up pointed observations on these candidates in C$^{34}$S, SO, OCS, HC$_3$N, HNCO, CH$_3$CN, and CH$_3$OH $J=2-1$ and $J=8-7$ lines to confirm and characterize them. We applied this method to the $l = 10^\circ-20^\circ$ portion of the FUGIN data and identified 22 "Hot Cores" (compact sources with more than two significant detection of the hot core tracer lines, i.e., SO, OCS, HC$_3$N, HNCO, CH$_3$CN, or CH$_3$OH $J=8-7$ lines) and 14 "Dense Clumps" (sources with more than two significant detection of C$^{34}$S, CH$_3$OH $J=2-1$, or the hot core tracer lines). The identified Hot Cores are found associated with signposts of high-mass star formation such as ATLASGAL clumps, WISE HII regions, and Class II methanol masers. For those associated with ATLASGAL clumps, their bolometric luminosity to clump mass ratios are consistent with the star formation stages centered at the hot core phase. The catalog of FUGIN Hot Cores provides a useful starting point for further statistical studies and detailed observations of high-mass star forming regions.


title: Adversarial Robustness of Deep Sensor Fusion Models

abstract: We experimentally study the robustness of deep camera-LiDAR fusion architectures for 2D object detection in autonomous driving. First, we find that the fusion model is usually both more accurate, and more robust against single-source attacks than single-sensor deep neural networks. Furthermore, we show that without adversarial training, early fusion is more robust than late fusion, whereas the two perform similarly after adversarial training. However, we note that single-channel adversarial training of deep fusion is often detrimental even to robustness. Moreover, we observe cross-channel externalities, where single-channel adversarial training reduces robustness to attacks on the other channel. Additionally, we observe that the choice of adversarial model in adversarial training is critical: using attacks restricted to cars' bounding boxes is more effective in adversarial training and exhibits less significant cross-channel externalities. Finally, we find that joint-channel adversarial training helps mitigate many of the issues above, but does not significantly boost adversarial robustness.


title: Uniqueness and stability of an inverse problem for a semi-linear wave equation

abstract: We consider the recovery of a potential associated with a semi-linear wave equation on $\mathbb{R}^{n+1}$, $n\geq 1$. We show a H\"older stability estimate for the recovery of an unknown potential $a$ of the wave equation $\square u +a u^m=0$ from its Dirichlet-to-Neumann map. We show that an unknown potential $a(x,t)$, supported in $\Omega\times[t_1,t_2]$, of the wave equation $\square u +a u^m=0$ can be recovered in a H\"older stable way from the map $u|_{\partial \Omega\times [0,T]}\mapsto \langle\psi,\partial_\nu u|_{\partial \Omega\times [0,T]}\rangle_{L^2(\partial \Omega\times [0,T])}$. This data is equivalent to the inner product of the Dirichlet-to-Neumann map with a measurement function $\psi$. We also prove similar stability result for the recovery of $a$ when there is noise added to the boundary data. The method we use is constructive and it is based on the higher order linearization. As a consequence, we also get a uniqueness result. We also give a detailed presentation of the forward problem for the equation $\square u +a u^m=0$.


title: Instant 3D Object Tracking with Applications in Augmented Reality

abstract: Tracking object poses in 3D is a crucial building block for Augmented Reality applications. We propose an instant motion tracking system that tracks an object's pose in space (represented by its 3D bounding box) in real-time on mobile devices. Our system does not require any prior sensory calibration or initialization to function. We employ a deep neural network to detect objects and estimate their initial 3D pose. Then the estimated pose is tracked using a robust planar tracker. Our tracker is capable of performing relative-scale 9-DoF tracking in real-time on mobile devices. By combining use of CPU and GPU efficiently, we achieve 26-FPS+ performance on mobile devices.


title: Cation Order Control of Correlations in Double Perovskite Sr$_2$VNbO$_6$

abstract: Double perovskites extend the design space for new materials, and they often host phenomena that don't exist in their parent perovskite compounds. Here, we present a detailed first principles study of the correlated double perovskite Sr$_2$VNbO$_6$, where inter-cationic charge transfer and strength of electronic correlations depend strongly on the cation order. By using Density Functional Theory + Embedded Dynamical Mean Field Theory, we show that this compound has a completely different electronic structure than either of its parent compounds despite V and Nb being from the same group in the periodic table. We explain how the electronic correlations' effect on the crystal structural parameters determines on which side of the Hund's metal-Mott insulator transition the material is. Our results demonstrate the emergence of Hund's metallic behavior in a double perovskite that has $d^1$ parents, and underlines the importance of electronic correlation effects on the crystal structure.


title: Abelian fixed point free endomorphisms and the Yang-Baxter equation

abstract: We obtain a simple family of solutions to the set-theoretic Yang-Baxter equation, one which depends only on considering special endomorphisms of a finite group. We show how such an endomorphism gives rise to two non-degenerate solutions to the Yang-Baxter equation, solutions which are inverse to each other. We give concrete examples using dihedral, alternating, symmetric, and metacyclic groups.


title: Analytical Survival Analysis of the Ornstein-Uhlenbeck Process

abstract: We use asymptotic methods from the theory of differential equations to obtain an analytical expression for the survival probability of an Ornstein-Uhlenbeck process with a potential defined over a broad domain. We form a uniformly continuous analytical solution covering the entire domain by asymptotically matching approximate solutions in an interior region, centered around the origin, to those in boundary layers, near the lateral boundaries of the domain. The analytic solution agrees extremely well with the numerical solution and takes into account the non-negligible leakage of probability that occurs at short times when the stochastic process begins close to one of the boundaries. Given the range of applications of Ornstein-Uhlenbeck processes, the analytic solution is of broad relevance across many fields of natural and engineering science.


title: Spectral Bias and Task-Model Alignment Explain Generalization in Kernel Regression and Infinitely Wide Neural Networks

abstract: Generalization beyond a training dataset is a main goal of machine learning, but theoretical understanding of generalization remains an open problem for many models. The need for a new theory is exacerbated by recent observations in deep neural networks where overparameterization leads to better performance, contradicting the conventional wisdom from classical statistics. In this paper, we investigate generalization error for kernel regression, which, besides being a popular machine learning method, also includes infinitely overparameterized neural networks trained with gradient descent. We use techniques from statistical mechanics to derive an analytical expression for generalization error applicable to any kernel or data distribution. We present applications of our theory to real and synthetic datasets, and for many kernels including those that arise from training deep neural networks in the infinite-width limit. We elucidate an inductive bias of kernel regression to explain data with "simple functions", which are identified by solving a kernel eigenfunction problem on the data distribution. This notion of simplicity allows us to characterize whether a kernel is compatible with a learning task, facilitating good generalization performance from a small number of training examples. We show that more data may impair generalization when noisy or not expressible by the kernel, leading to non-monotonic learning curves with possibly many peaks. To further understand these phenomena, we turn to the broad class of rotation invariant kernels, which is relevant to training deep neural networks in the infinite-width limit, and present a detailed mathematical analysis of them when data is drawn from a spherically symmetric distribution and the number of input dimensions is large.


title: On weak separation property for self-affine Jordan arcs

abstract: We consider self-affine arcs in $\mathbb R^2$ and prove that violation of "inner" weak separation property for such arcs implies that the arc is a parabolic segment. Therefore, if a self-affine Jordan arc is not a parabolic segment, then it is the attractor of some multizipper.


title: Combining Neural Language Models for WordSense Induction

abstract: Word sense induction (WSI) is the problem of grouping occurrences of an ambiguous word according to the expressed sense of this word. Recently a new approach to this task was proposed, which generates possible substitutes for the ambiguous word in a particular context using neural language models, and then clusters sparse bag-of-words vectors built from these substitutes. In this work, we apply this approach to the Russian language and improve it in two ways. First, we propose methods of combining left and right contexts, resulting in better substitutes generated. Second, instead of fixed number of clusters for all ambiguous words we propose a technique for selecting individual number of clusters for each word. Our approach established new state-of-the-art level, improving current best results of WSI for the Russian language on two RUSSE 2018 datasets by a large margin.


title: A stabilized finite element method for inverse problems subject to the convection-diffusion equation. II: convection-dominated regime

abstract: We consider the numerical approximation of the ill-posed data assimilation problem for stationary convection-diffusion equations and extend our previous analysis in [Numer. Math. 144, 451--477, 2020] to the convection-dominated regime. Slightly adjusting the stabilized finite element method proposed for dominant diffusion, we draw upon a local error analysis to obtain quasi-optimal convergence along the characteristics of the convective field through the data set. The weight function multiplying the discrete solution is taken to be Lipschitz and a corresponding super approximation result (discrete commutator property) is proven. The effect of data perturbations is included in the analysis and we conclude the paper with some numerical experiments.


title: Simple and Effective VAE Training with Calibrated Decoders

abstract: Variational autoencoders (VAEs) provide an effective and simple method for modeling complex distributions. However, training VAEs often requires considerable hyperparameter tuning to determine the optimal amount of information retained by the latent variable. We study the impact of calibrated decoders, which learn the uncertainty of the decoding distribution and can determine this amount of information automatically, on the VAE performance. While many methods for learning calibrated decoders have been proposed, many of the recent papers that employ VAEs rely on heuristic hyperparameters and ad-hoc modifications instead. We perform the first comprehensive comparative analysis of calibrated decoder and provide recommendations for simple and effective VAE training. Our analysis covers a range of image and video datasets and several single-image and sequential VAE models. We further propose a simple but novel modification to the commonly used Gaussian decoder, which computes the prediction variance analytically. We observe empirically that using heuristic modifications is not necessary with our method. Project website is at https://orybkin.github.io/sigma-vae/


title: Machine learning active-nematic hydrodynamics

abstract: Hydrodynamic theories effectively describe many-body systems out of equilibrium in terms of a few macroscopic parameters. However, such hydrodynamic parameters are difficult to derive from microscopics. Seldom is this challenge more apparent than in active matter where the energy cascade mechanisms responsible for autonomous large-scale dynamics are poorly understood. Here, we use active nematics to demonstrate that neural networks can extract the spatio-temporal variation of hydrodynamic parameters directly from experiments. Our algorithms analyze microtubule-kinesin and actin-myosin experiments as computer vision problems. Unlike existing methods, neural networks can determine how multiple parameters such as activity and elastic constants vary with ATP and motor concentration. In addition, we can forecast the evolution of these chaotic many-body systems solely from image-sequences of their past by combining autoencoder and recurrent networks with residual architecture. Our study paves the way for artificial-intelligence characterization and control of coupled chaotic fields in diverse physical and biological systems even when no knowledge of the underlying dynamics exists.

categories: cond-mat.soft cond-mat.mtrl-sci nlin.CD physics.flu-dyn


title: Asymptotic commutativity of quantized spaces: the case of $\mathbb{CP}^{p,q}$

abstract: We present a procedure for quantizing complex projective spaces $\mathbb{CP}^{p,q}$, $q\ge 1$, as well as construct relevant star products on these spaces. The quantization is made unique with the demand that it preserves the full isometry algebra of the metric. Although the isometry algebra, namely $su(p+1,q)$, is preserved by the quantization, the Killing vectors generating these isometries pick up quantum corrections. The quantization procedure is an extension of one applied recently to Euclidean $AdS_2$, where it was found that all quantum corrections to the Killing vectors vanish in the asymptotic limit, in addition to the result that the star product trivializes to pointwise product in the limit. In other words, the space is asymptotically anti-de Sitter making it a possible candidate for the $AdS/CFT$ correspondence principle. In this article, we find indications that the results for quantized Euclidean $AdS_2$ can be extended to quantized $\mathbb{CP}^{p,q}$, i.e., noncommutativity is restricted to a limited neighborhood of some origin, and these quantum spaces approach $\mathbb{CP}^{p,q}$ in the asymptotic limit.


title: Long-Horizon Visual Planning with Goal-Conditioned Hierarchical Predictors

abstract: The ability to predict and plan into the future is fundamental for agents acting in the world. To reach a faraway goal, we predict trajectories at multiple timescales, first devising a coarse plan towards the goal and then gradually filling in details. In contrast, current learning approaches for visual prediction and planning fail on long-horizon tasks as they generate predictions (1) without considering goal information, and (2) at the finest temporal resolution, one step at a time. In this work we propose a framework for visual prediction and planning that is able to overcome both of these limitations. First, we formulate the problem of predicting towards a goal and propose the corresponding class of latent space goal-conditioned predictors (GCPs). GCPs significantly improve planning efficiency by constraining the search space to only those trajectories that reach the goal. Further, we show how GCPs can be naturally formulated as hierarchical models that, given two observations, predict an observation between them, and by recursively subdividing each part of the trajectory generate complete sequences. This divide-and-conquer strategy is effective at long-term prediction, and enables us to design an effective hierarchical planning algorithm that optimizes trajectories in a coarse-to-fine manner. We show that by using both goal-conditioning and hierarchical prediction, GCPs enable us to solve visual planning tasks with much longer horizon than previously possible.


title: Linear cycles of consecutive lengths

abstract: A well-known result of Verstra\"ete \cite{V00} shows that for each integer $k\geq 2$ every graph $G$ with average degree at least $8k$ contains cycles of $k$ consecutive even lengths, the shortest of which is at most twice the radius of $G$. We establish two extensions of Verstra\"ete's result for linear cycles in linear $r$-uniform hypergraphs. We show that for any fixed integers $r\geq 3,k\geq 2$, there exist constants $c_1=c_1(r)$ and $c_2=c_2(r,k)$, such that every linear $r$-uniform hypergraph $G$ with average degree $d(G)\geq c_1 k$ contains linear cycles of $k$ consecutive even lengths, the shortest of which is at most $2\lceil \frac{ \log n}{\log (d(G)/k)-c_2}\rceil$. In particular, as an immediate corollary, we retrieve the current best known upper bound on the linear Tur\'an number of $C^r_{2k}$ with improved coefficients. Furthermore, we show that for any fixed integers $r\geq 3,k\geq 2$, there exist constants $c_3=c_3(r)$ and $c_4=c_4(r)$ such that every $n$-vertex linear $r$-uniform graph with average degree $d(G)\geq c_3k$, contains linear cycles of $k$ consecutive lengths, the shortest of which has length at most $6\lceil \frac{\log n}{\log (d(G)/k)-c_4} \rceil +6$. Both the degree condition and the shortest length among the cycles guaranteed are best possible up to a constant factor.


title: Correspondence principle for many-body scars in ultracold Rydberg atoms

abstract: The theory of quantum scarring -- a remarkable violation of quantum unique ergodicity -- rests on two complementary pillars: the existence of unstable classical periodic orbits and the so-called quasimodes, i.e., the non-ergodic states that strongly overlap with a small number of the system's eigenstates. Recently, interest in quantum scars has been revived in a many-body setting of Rydberg atom chains. While previous theoretical works have identified periodic orbits for such systems using time-dependent variational principle (TDVP), the link between periodic orbits and quasimodes has been missing. Here we provide a conceptually simple analytic construction of quasimodes for the non-integrable Rydberg atom model, and prove that they arise from a "requantisation" of previously established periodic orbits when quantum fluctuations are restored to all orders. Our results shed light on the TDVP classical system simultaneously playing the role of both the mean-field approximation and the system's classical limit, thus allowing us to firm up the analogy between the eigenstate scarring in the Rydberg atom chains and the single-particle quantum systems.

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


title: Feature Expansive Reward Learning: Rethinking Human Input

abstract: When a person is not satisfied with how a robot performs a task, they can intervene to correct it. Reward learning methods enable the robot to adapt its reward function online based on such human input, but they rely on handcrafted features. When the correction cannot be explained by these features, recent work in deep Inverse Reinforcement Learning (IRL) suggests that the robot could ask for task demonstrations and recover a reward defined over the raw state space. Our insight is that rather than implicitly learning about the missing feature(s) from demonstrations, the robot should instead ask for data that explicitly teaches it about what it is missing. We introduce a new type of human input in which the person guides the robot from states where the feature being taught is highly expressed to states where it is not. We propose an algorithm for learning the feature from the raw state space and integrating it into the reward function. By focusing the human input on the missing feature, our method decreases sample complexity and improves generalization of the learned reward over the above deep IRL baseline. We show this in experiments with a physical 7DOF robot manipulator, as well as in a user study conducted in a simulated environment.


title: What Happens when Separate and Unequal School Districts Merge?

abstract: We study the welfare effects of school district consolidation, i.e. the integration of disjoint school districts into a centralised clearinghouse. We show theoretically that, in the worst-case scenario, district consolidation may unambiguously reduce students' welfare, even if the student-optimal stable matching is consistently chosen. However, on average all students experience expected welfare gains from district consolidation, particularly those who belong to smaller and over-demanded districts. Using data from the Hungarian secondary school assignment mechanism, we compute the actual welfare gains from district consolidation in Budapest and compare these to our theoretical predictions. We empirically document substantial welfare gains from district consolidation for students, equivalent to attending a school five kilometres closer to the students' home addresses. As an important building block of our empirical strategy, we describe a method to consistently estimate students' preferences over schools and vice versa that does not fully assume that students report their preferences truthfully in the student-proposing deferred acceptance algorithm.


title: Precision measurement of the cosmic-ray electron and positron fluxes as a function of time and energy with the Alpha Magnetic Spectrometer on the International Space Station

abstract: This thesis presents an analysis of the cosmic-ray electron and positron flux using the AMS-02 detector on the International Space Station as a function of time and energy. The time-averaged flux is integrated over 6.5 years of AMS-02 science data and provides the electron and positron flux with unprecedented accuracy, covering the energy range from 0.5 GeV to 1 TeV. In total 28.39 million events were identified as electrons and 1.95 million as positrons. For each of the 88 Bartels rotation periods (27 days), within the 6.5 years, an individual electron and positron flux is derived spanning the energy range from 1 - 50 GeV.


title: Meta Transfer Learning for Emotion Recognition

abstract: Deep learning has been widely adopted in automatic emotion recognition and has lead to significant progress in the field. However, due to insufficient annotated emotion datasets, pre-trained models are limited in their generalization capability and thus lead to poor performance on novel test sets. To mitigate this challenge, transfer learning performing fine-tuning on pre-trained models has been applied. However, the fine-tuned knowledge may overwrite and/or discard important knowledge learned from pre-trained models. In this paper, we address this issue by proposing a PathNet-based transfer learning method that is able to transfer emotional knowledge learned from one visual/audio emotion domain to another visual/audio emotion domain, and transfer the emotional knowledge learned from multiple audio emotion domains into one another to improve overall emotion recognition accuracy. To show the robustness of our proposed system, various sets of experiments for facial expression recognition and speech emotion recognition task on three emotion datasets: SAVEE, EMODB, and eNTERFACE have been carried out. The experimental results indicate that our proposed system is capable of improving the performance of emotion recognition, making its performance substantially superior to the recent proposed fine-tuning/pre-trained models based transfer learning methods.


title: Automated Detection of COVID-19 from CT Scans Using Convolutional Neural Networks

abstract: COVID-19 is an infectious disease that causes respiratory problems similar to those caused by SARS-CoV (2003). Currently, swab samples are being used for its diagnosis. The most common testing method used is the RT-PCR method, which has high specificity but variable sensitivity. AI-based detection has the capability to overcome this drawback. In this paper, we propose a prospective method wherein we use chest CT scans to diagnose the patients for COVID-19 pneumonia. We use a set of open-source images, available as individual CT slices, and full CT scans from a private Indian Hospital to train our model. We build a 2D segmentation model using the U-Net architecture, which gives the output by marking out the region of infection. Our model achieves a sensitivity of 96.428% (95% CI: 88%-100%) and a specificity of 88.39% (95% CI: 82%-94%). Additionally, we derive a logic for converting our slice-level predictions to scan-level, which helps us reduce the false positives.


title: Influence of the theoretical model eye on the numerical evaluation of fractal intraocular lenses

abstract: In this work we present the numerical evaluation of a new design of fractal intraocular lens studied through a ray-tracing program. To determine the monochromatic and polychromatic performance of these lenses in different theoretical model eyes the Modulation Transfer Function (MTF) and the area above the MTF (AMTF) have been used. These merit functions show the same behavior for different values of asphericity (Q), independently from the theoretical model eye, even though there are differences due to the Spherical Aberration (SA) considered in each model.


title: A new trifocal corneal inlay for presbyopia

abstract: Corneal inlays (CIs) are the most recent surgical procedure for the treatment of presbyopia in patients who want complete independence from the use of glasses or contact lenses. Although refractive surgery in presbyopic patients is mostly performed in combination with cataract surgery, when the implantation of an intraocular lens is not necessary, the option of CIs has the advantage of being minimally invasive. Current designs of CIs are, either: small aperture devices, or refractive devices, however, both methods do not have good performance simultaneously at intermediate and near distances in eyes that are unable to accommodate. In the present study, we propose the first design of a trifocal CI, allowing good vision, at the same time, at far, intermediate and near vision in presbyopic eyes. We first demonstrate the good performance of the new inlay in comparison with a commercially available CI by using optical design software. We next confirm experimentally the image forming capabilities of our proposal employing an adaptive optics based optical simulator. This new design also has a number of parameters that can be varied to make personalized trifocal CI, opening up a new avenue for the treatment of presbyopia.


title: Traffic congestion anomaly detection and prediction using deep learning

abstract: Congestion prediction represents a major priority for traffic management centres around the world to ensure timely incident response handling. The increasing amounts of generated traffic data have been used to train machine learning predictors for traffic, however, this is a challenging task due to inter-dependencies of traffic flow both in time and space. Recently, deep learning techniques have shown significant prediction improvements over traditional models, however, open questions remain around their applicability, accuracy and parameter tuning. This paper brings two contributions in terms of: 1) applying an outlier detection an anomaly adjustment method based on incoming and historical data streams, and 2) proposing an advanced deep learning framework for simultaneously predicting the traffic flow, speed and occupancy on a large number of monitoring stations along a highly circulated motorway in Sydney, Australia, including exit and entry loop count stations, and over varying training and prediction time horizons. The spatial and temporal features extracted from the 36.34 million data points are used in various deep learning architectures that exploit their spatial structure (convolutional neuronal networks), their temporal dynamics (recurrent neuronal networks), or both through a hybrid spatio-temporal modelling (CNN-LSTM). We show that our deep learning models consistently outperform traditional methods, and we conduct a comparative analysis of the optimal time horizon of historical data required to predict traffic flow at different time points in the future. Lastly, we prove that the anomaly adjustment method brings significant improvements to using deep learning in both time and space.


title: Oscillation inequalities on real and ergodic $H^1$ spaces

abstract: Let $(x_n)$ be a sequence and $\rho\geq 1$. For a fixed sequences $n_1<n_2<n_3<\dots$, and $M$ define the oscillation operators $$\mathcal{O}_\rho (x_n)=\left(\sum_{k=1}^\infty\sup_{\substack{n_k\leq m< n_{k+1}\\m\in M}}\left|x_m-x_{n_k}\right|^\rho\right)^{1/\rho}.$$ Let $(X,\mathscr{B} ,\mu , \tau)$ be a dynamical system with $(X,\mathscr{B} ,\mu )$ a probability space and $\tau$ a measurable, invertible, measure preserving point transformation from $X$ to itself.\\ Suppose that the sequences $(n_k)$ and $M$ are lacunary. Then we prove the following results for $\rho\geq 2$: (i) Define $\phi_n(x)=\frac{1}{n}\chi_{[0,n]}(x)$ on $\mathbb{R}$. Then there exists a constant $C>0$ such that $\|\mathcal{O}_\rho (\phi_n\ast f)\|_{L^1(\mathbb{R})}\leq C\|f\|_{H^1(\mathbb{R})}$ for all $f\in H^1(\mathbb{R})$. (ii) Let $A_nf(x)=\frac{1}{n}\sum_{k=1}^nf(\tau^kx)$ be the usual ergodic averages in ergodic theory. Then $\|\mathcal{O}_\rho (A_nf)\|_{L^1(X)}\leq C\|f\|_{H^1(X)}$ for all $f\in H^1(X)$. (iii) If $[f(x)\log (x)]^+$ is integrable, then $\mathcal{O}_\rho (A_nf)$ is integrable.


title: Electromagnetic potentials and radiation of the Proca field created by a massless charge

abstract: We study a theory less explored in the literature, that of a massive abelian gauge boson interacting with a massless fermion. In this framework we calculate the electromagnetic potentials created by a point like massless charge. We determine also the radiation associated with a change in the momentum of the charge in the low energy regime. The results show both similarities and discrepancies with respect to quantum electrodynamics.


title: Surface cluster algebra expansion formulae via loop graphs

abstract: In 2011 Musiker, Schiffler and Williams obtained expansion formulae for cluster algebras from orientable surfaces. For singly and doubly notched arcs these formulae required the notion of $\gamma$-symmetric perfect matchings and $\gamma$-compatible pairs of $\gamma$-symmetric perfect matchings, respectively. We simplify and unify these approaches by considering good matchings of loop graphs.


title: Interpreting LIGO/Virgo "Mass-Gap" events as lensed Neutron Star-Black Hole binaries

abstract: Gravitational lensing allows the detection of binary black holes (BBH) at cosmological distances with chirp masses that appear to be enhanced by $1+z$ in the range $1<z<4$, in good agreement with the reported BBH masses. We propose this effect also accounts for the puzzling "mass gap" events (MG) newly reported by LIGO/Virgo, as distant, lensed NSBH events with $1<z<4$. The fitted mass of the neutron star member becomes $(1+z)\times 1.4M_\odot$, and is therefore misclassified as a low mass black hole. In this way, we derive a redshift of $z\simeq 3.5$ and $z\simeq 1.0$ for two newly reported "mass asymmetric" events GW190412 \& GW190814, by interpreting them as lensed NSBH events, comprising a stellar mass black hole and neutron star. Over the past year an additional 31 BBH events and 5 MG events have been reported with high probability ($>95\%$), from which we infer a factor $\simeq 5$ higher intrinsic rate of NSBH events than BBH events, reflecting a higher proportion of neutron stars formed by early star formation. We predict a distinctive locus for lensed NSBH events in the observed binary mass plane, spanning $1<z<4$ with a narrow mass ratio, $q \simeq 0.2$, that can be readily tested when the waveform data are unlocked. All such events may show disrupted NS emission and are worthy of prompt follow-up as the high lensing magnification means EM detections are not prohibitive despite the high redshifts that we predict. Such lensed NSBH events provide an exciting prospect of directly charting the history of coalescing binaries via the cosmological redshift of their waveforms, determined relative to the characteristic mass of the neutron star member.


title: Dynamical chemistry: non-equilibrium effective actions for reactive fluids

abstract: We present two approaches for describing chemical reactions taking place in fluid phase. The first method mirrors the usual derivation of the hydrodynamic equations of motion by relating conserved---or to account for chemical reactions, non-conserved---currents to local-equilibrium parameters. The second method involves a higher-brow approach in which we attack the same problem from the perspective of non-equilibrium effective field theory (EFT). Non-equilibrium effective actions are defined using the in-in formalism on the Schwinger-Keldysh contour and are therefore capable of describing thermal fluctuations and dissipation as well as quantum effects. The non-equilibrium EFT approach is especially powerful as all terms in the action are fully specified by the symmetries of the system; in particular the second law of thermodynamics does not need to be included by hand, but is instead derived from the action itself. We find that the equations of motion generated by both methods agree, but the EFT approach yields certain advantages. To demonstrate some of these advantages we construct a quadratic action that is valid to very small distance scales---much smaller than the scales at which ordinary hydrodynamic theories break down. Such an action captures the full thermodynamic and quantum behavior of reactions and diffusion at quadratic order. Finally, taking the low-frequency and low-wavenumber limit, we reproduce the linearized version of the well-known reaction-diffusion equations as a final coherence check.


title: Quenching and morphological evolution due to circumgalactic gas expulsion in a simulated galaxy with a controlled assembly history

abstract: We examine the influence of dark matter halo assembly on the evolution of a simulated $\sim L^\star$ galaxy. Starting from a zoom-in simulation of a star-forming galaxy evolved with the EAGLE galaxy formation model, we use the genetic modification technique to create a pair of complementary assembly histories: one in which the halo assembles later than in the unmodified case, and one in which it assembles earlier. Delayed assembly leads to the galaxy exhibiting a greater present-day star formation rate than its unmodified counterpart, whilst in the accelerated case the galaxy quenches at $z\simeq 1$, and becomes spheroidal. We simulate each assembly history nine times, adopting different seeds for the random number generator used by EAGLE's stochastic subgrid implementations of star formation and feedback. The systematic changes driven by differences in assembly history are significantly stronger than the random scatter induced by this stochasticity. The sensitivity of $\sim L^\star$ galaxy evolution to dark matter halo assembly follows from the close coupling of the growth histories of the central black hole (BH) and the halo, such that earlier assembly fosters the formation of a more massive BH, and more efficient expulsion of circumgalactic gas. In response to this expulsion, the circumgalactic medium reconfigures at a lower density, extending its cooling time and thus inhibiting the replenishment of the interstellar medium. Our results indicate that halo assembly history significantly influences the evolution of $\sim L^\star$ central galaxies, and that the expulsion of circumgalactic gas is a crucial step in quenching them.


title: Certified variational quantum algorithms for eigenstate preparation

abstract: Solutions to many-body problem instances often involve an intractable number of degrees of freedom and admit no known approximations in general form. In practice, representing quantum-mechanical states of a given Hamiltonian using available numerical methods, in particular those based on variational Monte Carlo simulations, become exponentially more challenging with increasing system size. Recently quantum algorithms implemented as variational models have been proposed to accelerate such simulations. The variational ansatz states are characterized by a polynomial number of parameters devised in a way to minimize the expectation value of a given Hamiltonian, which is emulated by local measurements. In this study, we develop a means to certify the termination of variational algorithms. We demonstrate our approach by applying it to three models: the transverse field Ising model, the model of one-dimensional spinless fermions with competing interactions, and the Schwinger model of quantum electrodynamics. By means of comparison, we observe that our approach shows better performance near critical points in these models. We hence take a further step to improve the applicability and to certify the results of variational quantum simulators.

categories: quant-ph cond-mat.dis-nn cond-mat.str-el cs.LG


title: The Contour Method: a new approach to finding modes of non-adiabatic stellar pulsations

abstract: The contour method is a new approach to calculating the non-adiabatic pulsation frequencies of stars. These frequencies can be found by solving for the complex roots of a characteristic equation constructed from the linear non-adiabatic stellar pulsation equations. A complex-root solver requires an initial trial frequency for each non adiabatic root. A standard method for obtaining initial trial frequencies is to use a star's adiabatic pulsation frequencies, but this method can fail to converge to non-adiabatic roots, especially as the growth and/or damping rate of the pulsations becomes large. The contour method provides an alternative way for obtaining initial trial frequencies that robustly converges to non-adiabatic roots, even for stellar models with extremely non-adiabatic pulsations and thus large growth/damping rates. We describe the contour method implemented in the GYRE stellar pulsation code and use it to calculate the non-adiabatic pulsation frequencies of $10\,\rm{M_{\odot}}$ and $20\,\rm{M_{\odot}}$ $\beta$ Cephei star models, and of a $0.9\,\rm{M_{\odot}}$ extreme helium star model.


title: Dark Radiation from Inflationary Fluctuations

abstract: Light new vector bosons can be produced gravitationally through quantum fluctuations during inflation; if these particles are feebly coupled and cosmologically metastable, they can account for the observed dark matter abundance. However, in minimal anomaly free $U(1)$ extensions to the Standard Model, these vectors generically decay to neutrinos if at least one neutrino mass eigenstate is sufficiently light. If these decays occur between neutrino decoupling and CMB freeze out, the resulting radiation energy density can contribute to $\Delta N_{\rm eff}$ at levels that can ameliorate the Hubble tension and be discovered with future CMB and relic neutrino detection experiments. Since the additional neutrinos are produced from vector decays after BBN, this scenario predicts $\Delta N_{\rm eff} > 0$ at recombination, but $\Delta N_{\rm eff} = 0$ during BBN. Furthermore, due to a fortuitous cancellation, the contribution to $\Delta N_{\rm eff}$ is approximately mass independent.


title: Rapid Elimination of Small Dust Grains in Molecular Clouds

abstract: We argue that impact velocities between dust grains with sizes less than $\sim 0.1$ $\mu m$ in molecular cloud cores are dominated by drift arising from ambipolar diffusion. This effect is due to the size dependence of the dust coupling to the magnetic field and the neutral gas. Assuming perfect sticking in collisions up to $\approx 50$ m/s, we show that this effect causes rapid depletion of small grains - consistent with starlight extinction and IR/microwave emission measurements, both in the core center ($n \sim 10^{6}$ cm$^{-3}$) and envelope ($n \sim 10^{4}$ cm$^{-3}$). The upper end of the size distribution does not change significantly if only velocities arising from this effect are considered. We consider the impact of an evolved dust size distribution on the gas temperature, and argue that if the depletion of small dust grains occurs as would be expected from our model, then the cosmic ray ionization rate must be well below $10^{-16}$ s$^{-1}$ at a number density of $10^{5}$ cm$^{-3}$.