title: Scattering of Dirac electrons from a skyrmion: emergence of robust skew scattering

abstract: We study electron scattering from a closed magnetic structure embedded in the top surface of a topological insulator (TI). Outside the structure there is a uniform layer of ferromagnetic insulator (FMI), leading to a positive effective mass for the Dirac electrons. The mass inside can be engineered to be negative, leading to a skyrmion structure. The geometric shape of the structure can be circular or deformed, leading to integrable or chaotic dynamics, respectively, in the classical limit. For a circular structure, the relativistic quantum scattering characteristics can be calculated analytically. For a deformed structure, we develop an efficient numerical method, the multiple multipole method, to solve the scattering wavefunctions. We find that anomalous Hall effect as characterized by strong skew scattering can arise, which is robust against structural deformation due to the resonant modes. In the short (long) wavelength regime, the resonant modes manifest themselves as confined vortices (excited edge states). The origin of the resonant states is the spin phase factor of massive Dirac electrons at the skyrmion boundary. Further, in the short wavelength regime, for a circular skyrmion, a large number of angular momentum channels contribute to the resonant modes. In this regime, in principle, classical dynamics are relevant, but we find that geometric deformations, even those as severe as leading to fully developed chaos, have little effect on the resonant modes. The vortex structure of the resonant states makes it possible to electrically ``charge'' the skyrmion, rendering feasible to manipulate its motion electrically. In the long wavelength regime, only the lowest angular momentum channels contribute to the resonant modes, making the skew scattering sharply directional. These phenomena may find applications for information storage and in Hall devices based on dynamic skyrmion.


title: Unveiling cloudy exoplanets: the influence of cloud model choices on retrieval solutions

abstract: In recent years, it has become clear that a substantial fraction of transiting exoplanets have some form of aerosol present in their atmospheres. Transit spectroscopy - mostly of hot Jupiters, but also of some smaller planets - has provided evidence for this, in the form of steep downward slopes from blue to red in the optical part of the spectrum, and muted gas absorption features throughout. Retrieval studies seeking to constrain the composition of exoplanet atmospheres must therefore account for the presence of aerosols. However, clouds and hazes are complex physical phenomena, and the transit spectra that are currently available allow us to constrain only some of their properties. Therefore, representation of aerosols in retrieval models requires that they are described by only a few parameters, and this has been done in a variety of ways within the literature. Here, I investigate a range of parameterisations for exoplanet aerosol and their effects on retrievals from transmission spectra of hot Jupiters HD 189733b and HD 209458b. I find that results qualitatively agree for the cloud/haze itself regardless of the parameterisation used, and indeed using multiple approaches provides a more holistic picture; the retrieved abundance of H2O is also very robust to assumptions about aerosols. I also find strong evidence that aerosol on HD 209458b covers less than half of the terminator region, whilst the picture is less clear for HD 189733b.


title: Massless Phases for the Villain model in $d\geq 3$

abstract: We consider the classical Villain rotator model in $\mathbb{Z}^d, d\geq 3$ at sufficiently low temperature, and prove that the truncated two-point function decays asymptotically as $|x|^{2-d}$, with an algebraic rate of convergence. We also obtain the same asymptotic decay separately for the transversal two-point functions. This quantifies the spontaneous magnetization result for the Villain model at low temperature, and rigorously establishes the Gaussian spin-wave conjecture in dimension $d\ge 3$. We believe that our method extends to finite range interactions and to other abelian spin systems and abelian gauge theory in $d\geq 3$.


title: Adiabatic theorem for closed quantum systems initialized at finite temperature

abstract: The evolution of a driven quantum system is said to be adiabatic whenever the state of the system stays close to an instantaneous eigenstate of its time-dependent Hamiltonian. The celebrated quantum adiabatic theorem ensures that such pure state adiabaticity can be maintained with arbitrary accuracy, provided one chooses a small enough driving rate. Here, we extend the notion of quantum adiabaticity to closed quantum systems initially prepared at finite temperature. In this case adiabaticity implies that the (mixed) state of the system stays close to a quasi-Gibbs state diagonal in the basis of the instantaneous eigenstates of the Hamiltonian. We prove a sufficient condition for the finite temperature adiabaticity. Remarkably, it implies that the finite temperature adiabaticity can be more robust than the pure state adiabaticity, particularly in many-body systems. We present an example of a many-body system where, in the thermodynamic limit, the finite temperature adiabaticity is maintained, while the pure state adiabaticity breaks down.

categories: quant-ph cond-mat.stat-mech math-ph math.MP


title: A deep-learning view of chemical space designed to facilitate drug discovery

abstract: Drug discovery projects entail cycles of design, synthesis, and testing that yield a series of chemically related small molecules whose properties, such as binding affinity to a given target protein, are progressively tailored to a particular drug discovery goal. The use of deep learning technologies could augment the typical practice of using human intuition in the design cycle, and thereby expedite drug discovery projects. Here we present DESMILES, a deep neural network model that advances the state of the art in machine learning approaches to molecular design. We applied DESMILES to a previously published benchmark that assesses the ability of a method to modify input molecules to inhibit the dopamine receptor D2, and DESMILES yielded a 77% lower failure rate compared to state-of-the-art models. To explain the ability of DESMILES to hone molecular properties, we visualize a layer of the DESMILES network, and further demonstrate this ability by using DESMILES to tailor the same molecules used in the D2 benchmark test to dock more potently against seven different receptors.


title: Activation Density driven Energy-Efficient Pruning in Training

abstract: Neural network pruning with suitable retraining can yield networks with considerably fewer parameters than the original with comparable degrees of accuracy. Typical pruning methods require large, fully trained networks as a starting point from which they perform a time-intensive iterative pruning and retraining procedure to regain the original accuracy. We propose a novel pruning method that prunes a network real-time during training, reducing the overall training time to achieve an efficient compressed network. We introduce an activation density based analysis to identify the optimal relative sizing or compression for each layer of the network. Our method is architecture agnostic, allowing it to be employed on a wide variety of systems. For VGG-19 and ResNet18 on CIFAR-10, CIFAR-100, and TinyImageNet, we obtain exceedingly sparse networks (up to $200 \times$ reduction in parameters and over $60 \times$ reduction in inference compute operations in the best case) with accuracy comparable to the baseline network. By reducing the network size periodically during training, we achieve total training times that are shorter than those of previously proposed pruning methods. Furthermore, training compressed networks at different epochs with our proposed method yields considerable reduction in training compute complexity ($1.6\times$ to $3.2\times$ lower) at near iso-accuracy as compared to a baseline network trained entirely from scratch.


title: Logistic Regression Regret: What's the Catch?

abstract: We address the problem of the achievable regret rates with online logistic regression. We derive lower bounds with logarithmic regret under $L_1$, $L_2$, and $L_\infty$ constraints on the parameter values. The bounds are dominated by $d/2 \log T$, where $T$ is the horizon and $d$ is the dimensionality of the parameter space. We show their achievability for $d=o(T^{1/3})$ in all these cases with Bayesian methods, that achieve them up to a $d/2 \log d$ term. Interesting different behaviors are shown for larger dimensionality. Specifically, on the negative side, if $d = \Omega(\sqrt{T})$, any algorithm is guaranteed regret of $\Omega(d \log T)$ (greater than $\Omega(\sqrt{T})$) under $L_\infty$ constraints on the parameters (and the example features). On the positive side, under $L_1$ constraints on the parameters, there exist algorithms that can achieve regret that is sub-linear in $d$ for the asymptotically larger values of $d$. For $L_2$ constraints, it is shown that for large enough $d$, the regret remains linear in $d$ but no longer logarithmic in $T$. Adapting the redundancy-capacity theorem from information theory, we demonstrate a principled methodology based on grids of parameters to derive lower bounds. Grids are also utilized to derive some upper bounds. Our results strengthen results by Kakade and Ng (2005) and Foster et al. (2018) for upper bounds for this problem, introduce novel lower bounds, and adapt a methodology that can be used to obtain such bounds for other related problems. They also give a novel characterization of the asymptotic behavior when the dimension of the parameter space is allowed to grow with $T$. They additionally establish connections to the information theory literature, demonstrating that the actual regret for logistic regression depends on the richness of the parameter class, where even within this problem, richer classes lead to greater regret.


title: Broadband X-ray analysis of 1E 1740.7-2942: constraints on spin, inclination and a tentative black hole mass

abstract: 1E 1740.7-2942 is one of the strongest hard X-ray emitters in the Galactic Centre region, believed to be a black hole in a high-mass X-ray binary system. Although extensively studied in X-rays, many aspects about the underlying nature of the system are still unknown. For example, X-ray data analyses of 1E 1740.7-2942 up to date have not yet unveiled the signature of a reflection component, whose modelling could be used to estimate parameters such as the spin of the black hole and inclination of the disc. We report here on the determination of these parameters from the analysis of the reflection component present in a public \textit{NuSTAR} observation which hasn't been subject to any previous study. We include \textit{XMM-Newton} and \textit{INTEGRAL} data to build a combined spectrum, enabling a joint analysis of both the disc and comptonisation components. Results point to a relatively high inclination disc $\gtrsim$ 50$^{\circ}$ (3 $\sigma$) and a near-maximum speed rotating black hole. The former is in agreement with a previous radio study and the latter is reported here for the first time. Lastly, we follow the methodology of recent efforts to weight black holes with only X-ray spectra and find results that suggest a black hole mass of about 5 M$_\odot$ for 1E 1740.7-2942.


title: Resonant backreaction in axion inflation

abstract: Axion inflation entails a coupling of the inflaton field to gauge fields through the Chern-Simons term. This results in a strong gauge field production during inflation, which backreacts on the inflaton equation of motion. Here we show that this strongly non-linear system generically experiences a resonant enhancement of the gauge field production, resulting in oscillatory features in the inflaton velocity as well as in the gauge field spectrum. The gauge fields source a strongly enhanced scalar power spectrum at small scales, exceeding previous estimates. For appropriate parameter choices, the collapse of these over-dense regions can lead to a large population of (light) primordial black holes with remarkable phenomenological consequences.


title: Interplay of Magnetism and Topological Superconductivity in Bilayer Kagome Metals

abstract: The binary intermetallic materials, $M_3$Sn$_2$ ($M$ = 3d transition metal) present a new class of strongly correlated systems that naturally allows for the interplay of magnetism and metallicity. Using first principles calculations we confirm that bulk Fe$_3$Sn$_2$ is a ferromagnetic metal, and show that $M$ = Ni and Cu are paramagnetic metals with non-trivial band structures. Focusing on Fe$_3$Sn$_2$ to understand the effect of enhanced correlations in an experimentally relevant atomistically thin single kagome-bilayer, our ab-initio results show that dimensional confinement naturally exposes the flatness of band structure associated with the bilayer kagome geometry in a resultant ferromagnetic Chern metal. We use a multistage minimal modeling of the magnetic bands progressively closer to the Fermi energy. This effectively captures the physics of the Chern metal with a non-zero anomalous Hall response over a material relevant parameter regime along with a possible superconducting instability of the spin-polarised band resulting in a topological superconductor.


title: Optical memories and switching dynamics of counterpropagating light states in microresonators

abstract: The Kerr nonlinearity can be a key enabler for many digital photonic circuits as it allows access to bistable states needed for all-optical memories and switches. A common technique is to use the Kerr shift to control the resonance frequency of a resonator and use it as a bistable, optically-tunable filter. However, this approach works only in a narrow power and frequency range or requires the use of an auxiliary laser. An alternative approach is to use the asymmetric bistability between counterpropagating light states resulting from the interplay between self- and cross-phase modulation, which allows light to enter a ring resonator in just one direction. Logical HIGH and LOW states can be represented and stored as the direction of circulation of light, and controlled by modulating the input power. Here we study the switching speed, operating laser frequency and power range, and contrast ratio of such a device. We reach a bitrate of 2 Mbps in our proof-of-principle device over an optical frequency range of 1 GHz and an operating power range covering more than one order of magnitude. We also calculate that integrated photonic circuits could exhibit bitrates of the order of Gbps, paving the way for the realization of robust and simple all-optical memories, switches, routers and logic gates that can operate at a single laser frequency with no additional electrical power.


title: A Multilingual View of Unsupervised Machine Translation

abstract: We present a probabilistic framework for multilingual neural machine translation that encompasses supervised and unsupervised setups, focusing on unsupervised translation. In addition to studying the vanilla case where there is only monolingual data available, we propose a novel setup where one language in the (source, target) pair is not associated with any parallel data, but there may exist auxiliary parallel data that contains the other. This auxiliary data can naturally be utilized in our probabilistic framework via a novel cross-translation loss term. Empirically, we show that our approach results in higher BLEU scores over state-of-the-art unsupervised models on the WMT'14 English-French, WMT'16 English-German, and WMT'16 English-Romanian datasets in most directions. In particular, we obtain a +1.65 BLEU advantage over the best-performing unsupervised model in the Romanian-English direction.


title: Small Data Wave Maps in Cyclic Spacetime

abstract: We study the initial value problem for the wave maps defined on the cyclic spacetime with the target Riemannian manifold that is responsive (see definition of the self coherence structure) to the parametric resonance phenomena. In particular, for arbitrary small and smooth initial data we construct blowing up solutions of the wave map if the metric of the base manifold is periodic in time.


title: $M^3$T: Multi-Modal Continuous Valence-Arousal Estimation in the Wild

abstract: This report describes a multi-modal multi-task ($M^3$T) approach underlying our submission to the valence-arousal estimation track of the Affective Behavior Analysis in-the-wild (ABAW) Challenge, held in conjunction with the IEEE International Conference on Automatic Face and Gesture Recognition (FG) 2020. In the proposed $M^3$T framework, we fuse both visual features from videos and acoustic features from the audio tracks to estimate the valence and arousal. The spatio-temporal visual features are extracted with a 3D convolutional network and a bidirectional recurrent neural network. Considering the correlations between valence / arousal, emotions, and facial actions, we also explores mechanisms to benefit from other tasks. We evaluated the $M^3$T framework on the validation set provided by ABAW and it significantly outperforms the baseline method.


title: Order Parameter and Magnetization of Antiferromagnets in Mutually Parallel Staggered and Magnetic Fields

abstract: We explore the behavior of the order parameter and the magnetization of antiferromagnetic solids subjected to mutually parallel staggered and magnetic fields. The effective field theory analysis of the partition function is taken up to the two-loop level, where the magnon-magnon interaction comes into play. These interaction effects, however, are small. A phenomenon that comes rather unexpectedly is that the finite-temperature magnetization increases with temperature when the strengths of the staggered and magnetic field are held constant.


title: Revisiting Spatial Invariance with Low-Rank Local Connectivity

abstract: Convolutional neural networks are among the most successful architectures in deep learning with this success at least partially attributable to the efficacy of spatial invariance as an inductive bias. Locally connected layers, which differ from convolutional layers only in their lack of spatial invariance, usually perform poorly in practice. However, these observations still leave open the possibility that some degree of relaxation of spatial invariance may yield a better inductive bias than either convolution or local connectivity. To test this hypothesis, we design a method to relax the spatial invariance of a network layer in a controlled manner; we create a \textit{low-rank} locally connected layer, where the filter bank applied at each position is constructed as a linear combination of basis set of filter banks with spatially varying combining weights. By varying the number of basis filter banks, we can control the degree of relaxation of spatial invariance. In experiments with small convolutional networks, we find that relaxing spatial invariance improves classification accuracy over both convolution and locally connected layers across MNIST, CIFAR-10, and CelebA datasets, thus suggesting that spatial invariance may be an overly restrictive prior.


title: Half-magnetization plateau and the origin of threefold symmetry breaking in an electrically-switchable triangular antiferromagnet

abstract: We perform high-field magnetization measurements on the triangular lattice antiferromagnet Fe$_{1/3}$NbS$_2$. We observe a plateau in the magnetization centered at approximately half the saturation magnetization over a wide range of temperature and magnetic field. From density functional theory calculations, we determine a likely set of magnetic exchange constants. Incorporating these constants into a minimal Hamiltonian model of our material, we find that the plateau and of the $Z_3$ symmetry breaking ground state both arise from interplane and intraplane antiferromagnetic interactions acting in competition. These findings are pertinent to the magneto-electric properties of Fe$_{1/3}$NbS$_2$, which allow electrical switching of antiferromagnetic textures at relatively low current densities.


title: From the Inner to Outer Milky Way: A Photometric Sample of 2.6 Million Red Clump Stars

abstract: Large pristine samples of red clump stars are highly sought after given that they are standard candles and give precise distances even at large distances. However, it is difficult to cleanly select red clumps stars because they can have the same T$_{\mathrm{eff}}$ and log $g$ as red giant branch stars. Recently, it was shown that the asteroseismic parameters, $\rm{\Delta}$P and $\rm{\Delta\nu}$, which are used to accurately select red clump stars, can be derived from spectra using the change in the surface carbon to nitrogen ratio ([C/N]) caused by mixing during the red giant branch. This change in [C/N] can also impact the spectral energy distribution. In this study, we predict the $\rm{\Delta}$P, $\rm{\Delta\nu}$, T$_{\mathrm{eff}}$ and log $g$ using 2MASS, AllWISE, \gaia, and Pan-STARRS data in order to select a clean sample of red clump stars. We achieve a contamination rate of $\sim$20\%, equivalent to what is achieved when selecting from T$_{\mathrm{eff}}$ and log $g$ derived from low resolution spectra. Finally, we present two red clump samples. One sample has a contamination rate of $\sim$ 20\% and $\sim$ 405,000 red clump stars. The other has a contamination of $\sim$ 33\% and $\sim$ 2.6 million red clump stars which includes $\sim$ 75,000 stars at distances $>$ 10 kpc. For |b|>30 degrees we find $\sim$ 15,000 stars with contamination rate of $\sim$ 9\%. The scientific potential of this catalog for studying the structure and formation history of the Galaxy is vast given that it includes millions of precise distances to stars in the inner bulge and distant halo where astrometric distances are imprecise.


title: Multilevel Acyclic Hypergraph Partitioning

abstract: A directed acyclic hypergraph is a generalized concept of a directed acyclic graph, where each hyperedge can contain an arbitrary number of tails and heads. Directed hypergraphs can be used to model data flow and execution dependencies in streaming applications. Thus, hypergraph partitioning algorithms can be used to obtain efficient parallelizations for multiprocessor architectures. However, an acyclicity constraint on the partition is necessary when mapping streaming applications to embedded multiprocessors due to resource restrictions on this type of hardware. The acyclic hypergraph partitioning problem is to partition the hypernodes of a directed acyclic hypergraph into a given number of blocks of roughly equal size such that the corresponding quotient graph is acyclic while minimizing an objective function on the partition. Here, we contribute the first n-level algorithm for the acyclic hypergraph partitioning problem. Our focus is on acyclic hypergraphs where hyperedges can have one head and arbitrary many tails. Based on this, we engineer a memetic algorithm to further reduce communication cost, as well as to improve scheduling makespan on embedded multiprocessor architectures. Experiments indicate that our algorithm outperforms previous algorithms that focus on the directed acyclic graph case which have previously been employed in the application domain. Moreover, our experiments indicate that using the directed hypergraph model for this type of application yields a significantly smaller makespan.


title: The $L^p$ Neumann problem for higher order elliptic equations

abstract: We solve the Neumann problem in the half space $\mathbb{R}^{n+1}_+$, for higher order elliptic differential equations with variable self-adjoint $t$-independent coefficients, and with boundary data in $L^p$, where $\max\bigl(1,\frac{2n}{n+2}-\varepsilon\bigr) < p < 2$. We also establish nontangential and area integral estimates on layer potentials with inputs in $L^p$ or $\dot W^{\pm1,p}$ for a similar range of~$p$, based on known bounds for $p\geq2$; in this case we may relax the requirement of self-adjointess.


title: A Translational Three-Degrees-of-Freedom Parallel Mechanism With Partial Motion Decoupling and Analytic Direct Kinematics

abstract: According to the topological design theory and method of parallel mechanism (PM) based on position and orientation characteristic (POC) equations, this paper studied a 3-DOF translational PM that has three advantages, i.e., (i) it consists of three fixed actuated prismatic joints, (ii) the PM has analytic solutions to the direct and inverse kinematic problems, and (iii) the PM is of partial motion decoupling property. Firstly, the main topological characteristics, such as the POC, degree of freedom and coupling degree were calculated for kinematic modeling. Thanks to these properties, the direct and inverse kinematic problems can be readily solved. Further, the conditions of the singular configurations of the PM were analyzed which corresponds to its partial motion decoupling property.


title: Structural diversity of disc galaxies originating in the cold gas inflow from cosmic webs

abstract: Disc galaxies show a large morphological diversity with varying contribution of three major structural components; thin discs, thick discs, and central bulges. Dominance of bulges increases with the galaxy mass (Hubble sequence) whereas thick discs are more prominent in lower mass galaxies. Because galaxies grow with the accretion of matter, this observed variety should reflect diversity in accretion history. On the basis of the prediction by the cold-flow theory for galactic gas accretion and inspired by the results of previous studies, we put a hypothesis that associates different accretion modes with different components. Namely, thin discs form as the shock-heated hot gas in high-mass halos gradually accretes to the central part, thick discs grow by the direct accretion of cold gas from cosmic webs when the halo mass is low, and finally bulges form by the inflow of cold gas through the shock-heated gas in high-redshift massive halos. We show that this simple hypothesis reproduces the mean observed variation of galaxy morphology with the galaxy mass. This scenario also predicts that thick discs are older and poorer in metals than thin discs, in agreement with the currently available observational data.


title: A polynomial algorithm for maxmin and minmax envy-free rent division on a soft budget

abstract: The current practice of envy-free rent division, lead by the fair allocation website Spliddit, is based on quasi-linear preferences. These preferences rule out agents' well documented financial constraints. To resolve this issue we consider piece-wise linear budget constrained preferences. These preferences admit differences in agents' marginal disutility of paying rent below and above a given reference, i.e., a soft budget. We construct a polynomial algorithm to calculate a maxmin utility envy-free allocation, and other related solutions, in this domain.


title: Deep Underground Neutrino Experiment (DUNE), Far Detector Technical Design Report, Volume I: Introduction to DUNE

abstract: The preponderance of matter over antimatter in the early universe, the dynamics of the supernovae that produced the heavy elements necessary for life, and whether protons eventually decay -- these mysteries at the forefront of particle physics and astrophysics are key to understanding the early evolution of our universe, its current state, and its eventual fate. The Deep Underground Neutrino Experiment (DUNE) is an international world-class experiment dedicated to addressing these questions as it searches for leptonic charge-parity symmetry violation, stands ready to capture supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. The DUNE far detector technical design report (TDR) describes the DUNE physics program and the technical designs of the single- and dual-phase DUNE liquid argon TPC far detector modules. This TDR is intended to justify the technical choices for the far detector that flow down from the high-level physics goals through requirements at all levels of the Project. Volume I contains an executive summary that introduces the DUNE science program, the far detector and the strategy for its modular designs, and the organization and management of the Project. The remainder of Volume I provides more detail on the science program that drives the choice of detector technologies and on the technologies themselves. It also introduces the designs for the DUNE near detector and the DUNE computing model, for which DUNE is planning design reports. Volume II of this TDR describes DUNE's physics program in detail. Volume III describes the technical coordination required for the far detector design, construction, installation, and integration, and its organizational structure. Volume IV describes the single-phase far detector technology. A planned Volume V will describe the dual-phase technology.


title: Interpreting the Spitzer/IRAC Colours of 7<z<9 Galaxies: Distinguishing Between Line Emission and Starlight Using ALMA

abstract: Prior to the launch of JWST, Spitzer/IRAC photometry offers the only means of studying the rest-frame optical properties of z>7 galaxies. Many such high redshift galaxies display a red [3.6] - [4.5] micron colour, often referred to as the "IRAC excess", which has conventionally been interpreted as arising from intense [OIII]+Hbeta emission within the [4.5] micron bandpass. An appealing aspect of this interpretation is similarly intense line emission seen in star-forming galaxies at lower redshift as well as the redshift-dependent behaviour of the IRAC colours beyond z~7 modelled as the various nebular lines move through the two bandpasses. In this paper we demonstrate that, given the photometric uncertainties, established stellar populations with Balmer (4000 A, rest-frame) breaks, such as those inferred at z>9 where line emission does not contaminate the IRAC bands, can equally well explain the redshift-dependent behaviour of the IRAC colours in 7<z<9 galaxies. We discuss possible ways of distinguishing between the two hypotheses using ALMA measures of [OIII] 88 micron and dust continuum fluxes. Prior to further studies with JWST, we show that the distinction is important in determining the assembly history of galaxies in the first 500 Myr.


title: Higher-order topological insulators in antiperovskites

abstract: We predict that a family of antiperovskite materials realize a higher order topological insulator phase, characterized by a previously introduced $\mathbb{Z}_4$ index. A tight binding model and a $k\cdot p$ model are used to capture the physics of the bulk, surface and hinge states of these materials. A phase diagram of the higher order and weak topological invariants is obtained for the tight binding model. The mirror Chern number is also discussed. In order to reveal the gapless hinge states in the presence of mirror Chern surface states, several ways of opening the surface gap are proposed and confirmed by calculation, including cleaving the crystal to reveal a low-symmetry surface, building a heterostructure, and applying strain. Upon opening the surface gap, we are able to study the hinge states by computing the momentum space band structure and real space distribution of mid-gap states.


title: Conformal Field Theories with Sporadic Group Symmetry

abstract: The monster sporadic group is the automorphism group of a central charge $c=24$ vertex operator algebra (VOA) or meromorphic conformal field theory (CFT). In addition to its $c=24$ stress tensor $T(z)$, this theory contains many other conformal vectors of smaller central charge; for example, it admits $48$ commuting $c=\frac12$ conformal vectors whose sum is $T(z)$. Such decompositions of the stress tensor allow one to construct new CFTs from the monster CFT in a manner analogous to the Goddard-Kent-Olive (GKO) coset method for affine Lie algebras. We use this procedure to produce evidence for the existence of a number of CFTs with sporadic symmetry groups and employ a variety of techniques, including Hecke operators, modular linear differential equations, and Rademacher sums, to compute the characters of these CFTs. Our examples include (extensions of) nine of the sporadic groups appearing as subquotients of the monster, as well as the simple groups ${}^2{E}_6(2)$ and ${F}_4(2)$ of Lie type. Many of these examples are naturally associated to McKay's $\widehat{E_8}$ correspondence, and we use the structure of Norton's monstralizer pairs more generally to organize our presentation.


title: Eigenstate Thermalization and Disorder Averaging in Gravity

abstract: Naively, a resolution of the black hole information paradox appears to involve microscopic details of a theory of quantum gravity. However, recent work has argued that a unitary Page curve can be recovered by including novel replica instantons in the gravitational path integral. Moreover, replica instantons seem to rely on disorder averaging the microscopic theory, without a definite connection to a single, underlying unitary quantum system. In this letter, we show that disorder averaging and replica instantons emerge naturally from a gravitational effective theory built out of typical microscopic states. We relate replica instantons to a moment expansion of the simple operators appearing in the Eigenstate Thermalization Hypothesis, describe Feynman rules for computing the moments, and find an elegant microcanonical description of replica instantons in terms of wormholes and Euclidean black holes.


title: Deconfined metallic quantum criticality: A $U(2)$ gauge-theoretic approach

abstract: We discuss a new class of quantum phase transitions -- Deconfined Mott Transition (DMT) -- that describe a continuous transition between a Fermi liquid metal with a generic electronic Fermi surface and an electrical insulator without Fermi surfaces of emergent neutral excitations. We construct a unified $U(2)$ gauge theory to describe a variety of metallic and insulating phases, which include Fermi liquids, fractionalized Fermi liquids (FL*), conventional insulators and quantum spin liquids, as well as the quantum phase transitions between them. Using the DMT as a basic building block, we propose a distinct quantum phase transition -- Deconfined Metal-Metal Transition (DM$^2$T) -- that describes a continuous transition between two metallic phases, accompanied by a jump in the size of their electronic Fermi surfaces (also dubbed a 'Fermi transition'). We study these new classes of deconfined metallic quantum critical points using a renormalization group framework at the leading nontrivial order in a controlled expansion, and comment on the various interesting scenarios that can emerge going beyond this leading order calculation. We also study a $U(1)\times U(1)$ gauge theory that shares a number of similarities with the $U(2)$ gauge theory and sheds important light on many phenomena related to DMT, DM$^2$T and quantum spin liquids.


title: Recurrent Neural Network Wave Functions

abstract: A core technology that has emerged from the artificial intelligence revolution is the recurrent neural network (RNN). Its unique sequence-based architecture provides a tractable likelihood estimate with stable training paradigms, a combination that has precipitated many spectacular advances in natural language processing and neural machine translation. This architecture also makes a good candidate for a variational wave function, where the RNN parameters are tuned to learn the approximate ground state of a quantum Hamiltonian. In this paper, we demonstrate the ability of RNNs to represent several many-body wave functions, optimizing the variational parameters using a stochastic approach. Among other attractive features of these variational wave functions, their autoregressive nature allows for the efficient calculation of physical estimators by providing independent samples. We demonstrate the effectiveness of RNN wave functions by calculating ground state energies, correlation functions, and entanglement entropies for several quantum spin models of interest to condensed matter physicists in one and two spatial dimensions.

categories: cond-mat.dis-nn cond-mat.str-el physics.comp-ph quant-ph


title: Deceptively cold dust in the massive starburst galaxy GN20 at $z\sim4$

abstract: We present new observations, carried out with IRAM NOEMA, of the atomic neutral carbon transitions [CI](1-0) at 492 GHz and [CI](2-1) at 809 GHz of GN20, a well-studied star-bursting galaxy at $z=4.05$. The high luminosity line ratio [CI](2-1)/[CI](1-0) implies an excitation temperature of $48^{+14}_{-9}$ K, which is significantly higher than the apparent dust temperature of $T_{\rm d}=33\pm2$ K ($\beta=1.9$) derived under the common assumption of an optically thin far-infrared dust emission, but fully consistent with $T_{\rm d}=52\pm5$ K of a general opacity model where the optical depth ($\tau$) reaches unity at a wavelength of $\lambda_0=170\pm23$ $\mu$m. Moreover, the general opacity solution returns a factor of $\sim 2\times$ lower dust mass and, hence, a lower molecular gas mass for a fixed gas-to-dust ratio, than with the optically thin dust model. The derived properties of GN20 thus provide an appealing solution to the puzzling discovery of starbursts appearing colder than main-sequence galaxies above $z>2.5$, in addition to a lower dust-to-stellar mass ratio that approaches the physical value predicted for starburst galaxies.


title: Gravitational-wave captures by intermediate-mass black holes in galactic nuclei

abstract: Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt, despite their potential role as massive seeds for quasars and sources of tidal disruption events, ultra-luminous X-ray sources, dwarf galaxy feedback, and hypervelocity stars. Gravitational wave (GW) observations can help to find and confirm the existence of IMBHs. Current and upcoming detectors, such as LIGO, Virgo, KAGRA, LISA, ET, and DECIGO promise to identify the full range from stellar-mass to supermassive black holes (SMBHs). In this paper, we address the question of whether IMBHs can produce GWs in galactic nuclei. We consider the possibility that stellar black holes (SBHs) form bound systems and later coalesce with an IMBH through gravitational captures in the dense nucleus. We show that this mechanism is efficient for IMBH masses in the range $\sim 3\times 10^3$ M$_\odot$--$2\times 10^4$ M$_\odot$. We find that the typical distributions of peak frequencies and merger timescales depend mainly on the IMBH mass. In particular, the typical peak frequency is about $0.2\,$Hz, $0.1\,$Hz, $0.09\,$Hz, and $0.05\,$Hz for $M_{\rm IMBH}=5\times 10^3$ M$_\odot$, $8\times 10^3$ M$_\odot$, $1\times 10^4$ M$_\odot$, and $2\times 10^4$ M$_\odot$, respectively. Our results show that, at design sensitivity, both DECIGO and ET should be able to detect these IMBH--SBH mergers. Furthermore, most of the mergers will appear eccentric ($e \gtrsim 0.1$), providing an indication of their dynamical origin.


title: What Drives the Redshift Evolution of Strong Emission Line Ratios?

abstract: We study the physical mechanisms that cause the offset between low-redshift and high-redshift galaxies on the [OIII]/H$\beta$ versus [NII]/H$\alpha$ ``Baldwin, Phillips & Terlevich'' (BPT) diagram using a sample of local analogues of high-redshift galaxies. These high-redshift analogue galaxies are selected from the Sloan Digital Sky Survey. Located in the same region on the BPT diagram as the ultra-violet selected galaxies at $z\sim2$, these high-redshift analogue galaxies provide an ideal local benchmark to study the offset between the local and high-redshift galaxies on the BPT diagram. We compare the nitrogen-to-oxygen ratio (N/O), the shape of the ionising radiation field, and ionisation parameters between the high-redshift analogues and a sample of local reference galaxies. The higher ionisation parameter in the high-redshift analogues is the dominant physical mechanism driving the BPT offset from low- to high-redshift, particularly at high {\nii/\ha}. Furthermore, the N/O ratio enhancement also plays a minor role to cause the BPT offset. However, the shape of the ionising radiation field is unlikely to cause the BPT offset because the high-redshift analogues have a similar hard ionising radiation field as local reference galaxies. This hard radiation field cannot be produced by the current standard stellar synthesis models. The stellar rotation and binarity may help solve the discrepancy.


title: Atmospheric Erosion by Giant Impacts onto Terrestrial Planets

abstract: We examine the mechanisms by which atmosphere can be eroded by giant impacts onto Earth-like planets with thin atmospheres, using 3D smoothed particle hydrodynamics simulations with sufficient resolution to directly model the fate of low-mass atmospheres. We present a simple scaling law to estimate the fraction lost for any impact angle and speed in this regime. In the canonical Moon-forming impact, only around 10% of the atmosphere would have been lost from the immediate effects of the collision. There is a gradual transition from removing almost none to almost all of the atmosphere for a grazing impact as it becomes more head-on or increases in speed, including complex, non-monotonic behaviour at low impact angles. In contrast, for head-on impacts, a slightly greater speed can suddenly remove much more atmosphere. Our results broadly agree with the application of 1D models of local atmosphere loss to the ground speeds measured directly from our simulations. However, previous analytical models of shock-wave propagation from an idealised point-mass impact significantly underestimate the ground speeds and hence the total erosion. The strong dependence on impact angle and the interplay of multiple non-linear and asymmetrical loss mechanisms highlight the need for 3D simulations in order to make realistic predictions.


title: The X-Ray Cavity Around Hotspot E in Cygnus A: Tunneled by a Deflected Jet

abstract: The powerful FR II radio galaxy Cygnus A exhibits primary and secondary hotspots in each lobe. A 2 Msec Chandra X-ray image of Cygnus A has revealed an approximately circular hole, with a radius of 3.9 kpc, centered on the primary hotspot in the eastern radio lobe, hotspot E. We infer the distribution of X-ray emission on our line-of-sight from an X-ray surface brightness profile of the radio lobe adjacent to the hole and use it to argue that the hole is excavated from the radio lobe. The surface brightness profile of the hole implies a depth at least 1.7 $\pm$ 0.3 times greater than its projected width, requiring a minimum depth of 13.3 $\pm$ 2.3 kpc. A similar hole observed in the 5 GHz VLA radio map reinforces the argument for a cavity lying within the lobe. We argue that the jet encounters the shock compressed intracluster medium at hotspot E, passing through one or more shocks as it is deflected back into the radio lobe. The orientation of Cygnus A allows the outflow from hotspot E to travel almost directly away from us, creating an elongated cavity, as observed. These results favor models for multiple hotspots in which an FR II jet is deflected at a primary hotspot, then travels onward to deposit the bulk of its power at a secondary hotspot, rather than the dentist drill model.


title: Localization transition on the Random Regular Graph as an unstable tricritical point in a log-normal Rosenzweig-Porter random matrix ensemble

abstract: Gaussian Rosenzweig-Porter (GRP) random matrix ensemble is the only one in which the robust multifractal phase and ergodic transition have a status of a mathematical theorem. Yet, this phase in GRP model is oversimplified: the spectrum of fractal dimensions is degenerate and the mini-band in the local spectrum is not multifractal. In this paper we suggest an extension of the GRP model by adopting a logarithmically-normal (LN) distribution of off-diagonal matrix elements. A family of such LN-RP models is parametrized by a symmetry parameter $p$ and it interpolates between the GRP at $p\rightarrow 0$ and Levy ensembles at $p\rightarrow\infty$. A special point $p=1$ is shown to be the simplest approximation to the Anderson localization model on a random regular graph.We study in detail the phase diagram of LN-RP model and show that $p=1$ is a tricritical point where the multifractal phase first collapses. This collapse is shown to be unstable with respect to the truncation of the log-normal distribution. We suggest a new criteria of stability of the non-ergodic phases and prove that the Anderson transition in LN-RP model is discontinuous at all $p>0$.


title: The role of galaxy mass on AGN emission: a view from the VANDELS survey

abstract: We present a comparative analysis of the properties of AGN emitting at radio and X-ray wavelengths. The study is performed on 907 X-ray AGN and 100 radio AGN selected on the CDFS and UDS fields and makes use of new and ancillary data available to the VANDELS collaboration. Our results indicate that the mass of the host galaxy is a fundamental quantity which determines the level of AGN activity at the various wavelengths. Indeed large stellar masses are found to be connected with AGN radio emission, as virtually all radio-active AGN reside within galaxies of M*>10^{10} Msun. Large stellar masses also seem to favour AGN activity in the X-ray, even though X-ray AGN present a mass distribution which is more spread out and with a non-negligible tail at M*<10^{9} Msun. Stellar mass alone is also observed to play a fundamental role in simultaneous radio and X-ray emission: the percentage of AGN active at both wavelengths increases from around 1% of all X-ray AGN residing within hosts of M*<10^{11} Msun to about 13% in more massive galaxies. In the case of radio-selected AGN, such a percentage moves from about 15% to about 45% (but up to 80% in the deepest fields). Neither cosmic epoch, nor radio luminosity, X-ray luminosity, Eddington ratio or star-formation rate of the hosts are found to be connected to an enhanced probability for joint radio+X-ray emission of AGN origin. Furthermore, only a loose relation is observed between X-ray and radio luminosity in those AGN which are simultaneously active at both frequencies.


title: Black hole mass function and its evolution -- the first prediction for the Einstein Telescope

abstract: The knowledge about the black hole mass function (BHMF) and its evolution would help to understand the origin of the BHs and how BH binaries formed at different stages of the history of the Universe. We demonstrate the ability of future third generation gravitational wave (GW) detector -- the Einstein Telescope (ET) to infer the slope of the BHMF and its evolution with redshift. We perform the Monte Carlo simulation of the measurements of chirp signals from binary BH systems (BBH) that could be detected by ET, including the BH masses and their luminosity distances ($d_L$). We use the mass of a primary black hole in each binary system to infer the BHMF as a power-law function with slope parameter as $\alpha$. Taking into account the bias that could be introduced by the uncertainty of measurements and by the selection effect, we carried out the numerical tests and find that only one thousand of GW events registered by ET ($\sim1\%$ amount of its yearly detection rate) could accurately infer the $\alpha$ with a precision of $\alpha\sim0.1$. Furthermore, we investigate the validity of our method to recover a scenario where $\alpha$ evolves with redshift as $\alpha(z) = \alpha_0 + \alpha_1\frac{z}{1+z}$. Taking a thousand of GW events and using $d_L$ as the redshift estimator, our tests show that one could infer the value of evolving parameter $\alpha_1$ accurately with the uncertainty level of $\sim0.5$. Our numerical tests verify the reliability of our method. The uncertainty levels of the inferred parameters can be trusted directly for the several sets of the parameter we assumed, yet shouldn't be treated as a universal level for the general case.


title: OAUNI astronomical photometry: Stellar variability of FO Aqr on 2016 low state

abstract: This work reports a photometric monitoring of the intermediate polar cataclysmic variable star FO Aqr. This is part of the ongoing OAUNI stellar variability program. Around 1200 individual measurements (or $\sim$ 9 hrs of observations) were gathered distributed in five nights. The observation epoch was coincident with the 2016 low state brightness of FO Aqr. Good quality light curves detect unambiguously the expected stellar variability of FO~Aqr. The analysis using periodograms let to determine the main two known periods, P$_{1}$~=~20.380~$\pm$~.003~min and P$_{2}$~=~11.076~$\pm$~.001~min. Our full five-nights analysis shows that P$_{1}$, the spin white dwarf period, is more prominent than P$_{2}$, the one-half of the beat period between the spin and orbital periods. Nevertheless, individual light curves show than, at least in one case, a reverse pattern is found. It suggests a dependence of the beat period with the orbital phase of the system with important variations at the same phase interval.


title: Rocking shadows in broken circumbinary discs

abstract: We use three dimensional simulations with coupled hydrodynamics and Monte Carlo radiative transfer to show that shadows cast by the inner disc in broken circumbinary discs move within a confined range of position angles on the outer disc. Over time, shadows appear to rock back and forth in azimuth as the inner disc precesses. The effect occurs because the inner disc precesses around a vector that is not the angular momentum vector of the outer disc. We relate our findings to recent observations of shadows in discs.


title: Relative Anomaly in (1+1)d Rational Conformal Field Theory

abstract: We study 't Hooft anomalies of symmetry-enriched rational conformal field theories (RCFT) in (1+1)d. Such anomalies determine whether a theory can be realized in a truly (1+1)d system with on-site symmetry, or on the edge of a (2+1)d symmetry-protected topological phase. RCFTs with the identical symmetry actions on their chiral algebras may have different 't Hooft anomalies due to additional symmetry charges on local primary operators. To compute the relative anomaly, we establish a precise correspondence between (1+1)d non-chiral RCFTs and (2+1)d doubled symmetry-enriched topological (SET) phases with a choice of symmetric gapped boundary. Based on these results we derive a general formula for the relative 't Hooft anomaly in terms of algebraic data that characterizes the SET phase and its boundary.


title: DFTpy: An efficient and object-oriented platform for orbital-free DFT simulations

abstract: In silico materials design is hampered by the computational complexity of Kohn-Sham DFT, which scales cubically with the system size. Owing to the development of new-generation kinetic energy density functionals (KEDFs), orbital-free DFT (OFDFT, a linear-scaling method) can now be successfully applied to a large class of semiconductors and such finite systems as quantum dots and metal clusters. In this work, we present DFTpy, an open source software implementing OFDFT written entirely in Python 3 and outsourcing the computationally expensive operations to third-party modules, such as NumPy and SciPy. When fast simulations are in order, DFTpy exploits the fast Fourier transforms (FFTs) from PyFFTW. New-generation, nonlocal and density-dependent-kernel KEDFs are made computationally efficient by employing linear splines and other methods for fast kernel builds. We showcase DFTpy by solving for the electronic structure of a million-atom system of aluminum metal which was computed on a single CPU. The Python 3 implementation is object-oriented, opening the door to easy implementation of new features. As an example, we present a time-dependent OFDFT implementation (hydrodynamic DFT) which we use to compute the spectra of small metal cluster recovering qualitatively the time-dependent Kohn-Sham DFT result. The Python code base allows for easy implementation of APIs. We showcase the combination of DFTpy and ASE for molecular dynamics simulations (NVT) of liquid metals. DFTpy is released under the MIT license.

categories: physics.comp-ph cond-mat.mtrl-sci physics.chem-ph quant-ph


title: Analytical model for cluster radio relics

abstract: Radio relics are vast synchrotron sources that sit on the outskirts of merging galaxy clusters. In this work we model their formation using a Press-Schechter formalism to simulate merger rates, analytical models for the intracluster medium and the shock dynamics, as well as a simple model for the cosmic-ray electrons at the merger shocks. We show that the statistical properties of the population of radio relics are strongly dependent on key physical parameters, such as the acceleration efficiency, the magnetic field strength at the relic, the geometry of the relic and the duration of the electron acceleration at merger shocks. It turns out that the flux distribution as well as the power-mass relation can constrain key parameters of the intracluster medium. With the advent of new large-area radio surveys, statistical analyses of radio relics will complement what we have learned from observations of individual objects.


title: Benchmarking $^{136}$Xe Neutrinoless $\beta\beta$ Decay Matrix Element Calculations with the $^{138}{\rm Ba}(p,t)$ Reaction

abstract: We used a high-resolution magnetic spectrograph to study neutron pair-correlated $0^+$ states in $^{136}$Ba, produced via the $^{138}{\rm Ba}(p,t)$ reaction. In conjunction with state-of-the-art shell model calculations, these data benchmark part of the dominant Gamow-Teller component of the nuclear matrix element (NME) for $^{136}$Xe neutrinoless double beta ($0\nu\beta\beta$) decay. We demonstrate for the first time an evaluation of part of a $0\nu\beta\beta$ decay NME by use of an experimental observable, presenting a new avenue of approach for more accurate calculations of $0\nu\beta\beta$ decay matrix elements.


title: Sparse PSD approximation of the PSD cone

abstract: While semidefinite programming (SDP) problems are polynomially solvable in theory, it is often difficult to solve large SDP instances in practice. One technique to address this issue is to relax the global positive-semidefiniteness (PSD) constraint and only enforce PSD-ness on smaller $k\times k$ principal submatrices --- we call this the sparse SDP relaxation. Surprisingly, it has been observed empirically that in some cases this approach appears to produce bounds that are close to the optimal objective function value of the original SDP. In this paper, we formally attempt to compare the strength of the sparse SDP relaxation vis-\`a-vis the original SDP from a theoretical perspective. In order to simplify the question, we arrive at a data independent version of it, where we compare the sizes of SDP cone and the $k$-PSD closure, which is the cone of matrices where PSD-ness is enforced on all $k\times k$ principal submatrices. In particular, we investigate the question of how far a matrix of unit Frobenius norm in the $k$-PSD closure can be from the SDP cone. We provide two incomparable upper bounds on this farthest distance as a function of $k$ and $n$. We also provide matching lower bounds, which show that the upper bounds are tight within a constant in different regimes of $k$ and $n$. Other than linear algebra techniques, we extensively use probabilistic methods to arrive at these bounds. One of the lower bounds is obtained by observing a connection between matrices in the $k$-PSD closure and matrices satisfying the restricted isometry property (RIP).


title: Desynchronization and Wave Pattern Formation in MPI-Parallel and Hybrid Memory-Bound Programs

abstract: Analytic, first-principles performance modeling of distributed-memory parallel codes is notoriously imprecise. Even for applications with extremely regular and homogeneous compute-communicate phases, simply adding communication time to computation time does often not yield a satisfactory prediction of parallel runtime due to deviations from the expected simple lockstep pattern caused by system noise, variations in communication time, and inherent load imbalance. In this paper, we highlight the specific cases of provoked and spontaneous desynchronization of memory-bound, bulk-synchronous pure MPI and hybrid MPI+OpenMP programs. Using simple microbenchmarks we observe that although desynchronization can introduce increased waiting time per process, it does not necessarily cause lower resource utilization but can lead to an increase in available bandwidth per core. In case of significant communication overhead, even natural noise can shove the system into a state of automatic overlap of communication and computation, improving the overall time to solution. The saturation point, i.e., the number of processes per memory domain required to achieve full memory bandwidth, is pivotal in the dynamics of this process and the emerging stable wave pattern. We also demonstrate how hybrid MPI-OpenMP programming can prevent desirable desynchronization by eliminating the bandwidth bottleneck among processes. A Chebyshev filter diagonalization application is used to demonstrate some of the observed effects in a realistic setting.


title: On the Number of $\tau$-Tilting Modules over Nakayama Algebras

abstract: Let $\Lambda^r_n$ be the path algebra of the linearly oriented quiver of type $\mathbb{A}$ with $n$ vertices modulo the $r$-th power of the radical, and let $\widetilde{\Lambda}^r_n$ be the path algebra of the cyclically oriented quiver of type $\widetilde{\mathbb{A}}$ with $n$ vertices modulo the $r$-th power of the radical. Adachi gave a recurrence relation for the number of $\tau$-tilting modules over $\Lambda^r_n$. In this paper, we show that the same recurrence relation also holds for the number of $\tau$-tilting modules over $\widetilde{\Lambda}^r_n$. As an application, we give a new proof for a result by Asai on recurrence formulae for the number of support $\tau$-tilting modules over $\Lambda^r_n$ and $\widetilde{\Lambda}^r_n$.


title: Learning Whole-body Motor Skills for Humanoids

abstract: This paper presents a hierarchical framework for Deep Reinforcement Learning that acquires motor skills for a variety of push recovery and balancing behaviors, i.e., ankle, hip, foot tilting, and stepping strategies. The policy is trained in a physics simulator with realistic setting of robot model and low-level impedance control that are easy to transfer the learned skills to real robots. The advantage over traditional methods is the integration of high-level planner and feedback control all in one single coherent policy network, which is generic for learning versatile balancing and recovery motions against unknown perturbations at arbitrary locations (e.g., legs, torso). Furthermore, the proposed framework allows the policy to be learned quickly by many state-of-the-art learning algorithms. By comparing our learned results to studies of preprogrammed, special-purpose controllers in the literature, self-learned skills are comparable in terms of disturbance rejection but with additional advantages of producing a wide range of adaptive, versatile and robust behaviors.


title: Mario Level Generation From Mechanics Using Scene Stitching

abstract: This paper presents a level generation method for Super Mario by stitching together pre-generated "scenes" that contain specific mechanics, using mechanic-sequences from agent playthroughs as input specifications. Given a sequence of mechanics, our system uses an FI-2Pop algorithm and a corpus of scenes to perform automated level authoring. The system outputs levels that have a similar mechanical sequence to the target mechanic sequence but with a different playthrough experience. We compare our system to a greedy method that selects scenes that maximize the target mechanics. Our system is able to maximize the number of matched mechanics while reducing emergent mechanics using the stitching process compared to the greedy approach.


title: The matrix regularization for Riemann surfaces with magnetic fluxes

abstract: We consider the matrix regularization of fields on a Riemann surface which couple to gauge fields with a nonvanishing magnetic flux. We show that such fields are described as rectangular matrices in the matrix regularization. We construct the matrix regularization explicitly for the case of the sphere and torus based on the Berezin-Toeplitz quantization, and also discuss a possible generalization to cases with higher genera. We also discuss the matrix version of the Laplacian acting on the rectangular matrices.


title: The pairing glue in cuprate superconductors from the self-energy revealed via machine learning

abstract: Recently, machine learning was applied to extract both the normal and the anomalous components of the self-energy from photoemission data at the antinodal points in Bi-based cuprate high-temperature superconductors [Y. Yamaji {\it et al.}, arXiv:1903.08060]. It was argued that both components do show prominent peaks near 50 meV, which hold information about the pairing glue, but the peaks are hidden in the actual data, which measure only the total self-energy. We analyze the self-energy within an effective fermion-boson theory. We show that soft thermal fluctuations give rise to peaks in both components of the self-energy at a frequency comparable to superconducting gap, while they cancel in the total self-energy; all irrespective of the nature of the pairing boson. However, in the quantum limit $T \to 0$ prominent peaks survive only for a very restricted subclass of pairing interactions. We argue that the way to potentially nail down the pairing boson is to determine the thermal evolution of the peaks.


title: Berry Flux Diagonalization: Application to Electric Polarization

abstract: The switching polarization of a ferroelectric is a characterization of the current that flows due to changes in polarization when the system is switched between two states. Computation of this change in polarization in crystal systems has been enabled by the modern theory of polarization, where it is expressed in terms of a change in Berry phase as the material switches. It is straightforward to compute this change of phase, but only modulo $2\pi$, requiring a branch choice from among a lattice of values separated by $2\pi$. The measured switching polarization depends on the actual path along which the material switches, which in general involves nucleation and growth of domains and is therefore quite complex. In this work, we present a physically motivated approach for predicting the experimentally measured switching polarization that involves separating the change in phase between two states into as many gauge-invariant smaller phase changes as possible. As long as the magnitudes of these smaller phase changes remain smaller than $\pi$, their sum forms a phase change which corresponds to the change one would find along any path involving minimal evolution of the atomic and electronic structure. We show that for typical ferroelectrics, including those that would have otherwise required a densely sampled path, this technique allows the switching polarization to be computed without any need for intermediate sampling between oppositely polarized states.


title: Omnidirectional flexural invisibility of multiple interacting voids in vibrating elastic plates

abstract: In elasticity, the design of a cloaking for an inclusion or a void to leave a vibrational field unperturbed by its presence, so to achieve its invisibility, is a thoroughly analyzed, but still unchallenged, mechanical problem. The 'cloaking transformation' concept, originally developed in electromagnetism and optics, is not directly applicable to elastic waves, displaying a complex vectorial nature. Consequently, all examples of elastic cloaking presented so far involve complex design and thick coating skins. These cloakings often work only for problems of unidirectional propagation, within narrow ranges of frequency, and considering only one colaked object. Here, a new method based on the concept of reinforcement, achieved via elastic stiffening and mass redistribution, is introduced to cloak multiple voids in an elastic plate. This simple technique produces invisibility of the voids to flexural waves within an extremely broad range of frequencies and thus surpassing in many aspects all existing cloaking techniques. The proposed design principle is applicable in mechanical problems ranging from the micro-scale to the scale of civil engineering. For instance, our results show how to design a perforated load-bearing building wall, vibrating during an earthquake exactly as the same wall, but unperforated, a new finding for seismic protection.


title: DropCluster: A structured dropout for convolutional networks

abstract: Dropout as a regularizer in deep neural networks has been less effective in convolutional layers than in fully connected layers. This is due to the fact that dropout drops features randomly. When features are spatially correlated as in the case of convolutional layers, information about the dropped pixels can still propagate to the next layers via neighboring pixels. In order to address this problem, more structured forms of dropout have been proposed. A drawback of these methods is that they do not adapt to the data. In this work, we introduce a novel structured regularization for convolutional layers, which we call DropCluster. Our regularizer relies on data-driven structure. It finds clusters of correlated features in convolutional layer outputs and drops the clusters randomly at each iteration. The clusters are learned and updated during model training so that they adapt both to the data and to the model weights. Our experiments on the ResNet-50 architecture demonstrate that our approach achieves better performance than DropBlock or other existing structured dropout variants. We also demonstrate the robustness of our approach when the size of training data is limited and when there is corruption in the data at test time.


title: Renofeation: A Simple Transfer Learning Method for Improved Adversarial Robustness

abstract: Fine-tuning through knowledge transfer from a pre-trained model on a large-scale dataset is a widely spread approach to effectively build models on small-scale datasets. In this work, we show that a recent adversarial attack designed for transfer learning via re-training the last linear layer can successfully deceive models trained with transfer learning via end-to-end fine-tuning. This raises security concerns for many industrial applications. In contrast, models trained with random initialization without transfer are much more robust to such attacks, although these models often exhibit much lower accuracy. To this end, we propose noisy feature distillation, a new transfer learning method that trains a network from random initialization while achieving clean-data performance competitive with fine-tuning. Code available at https://github.com/cmu-enyac/Renofeation.


title: Autonomous chaos of exciton-polariton condensates

abstract: We study the formation of chaos and strange attractors in the order parameter space of a system of two coupled, non-resonantly driven exciton-polariton condensates. The typical scenario of bifurcations experienced by the system with increasing external pumping consists of (i) formation of $\pi$-synchronized condensates at low pumping, (ii) symmetry breaking pitchfork bifurcation leading to unequal occupations of the condensates with non-trivial phase difference between them, (iii) loss of the stability of all fixed points in the system resulting in chaotic dynamics, (iv) limit cycle dynamics of the order parameter, which ends up in (v) in-phase synchronized condensates via the Hopf bifurcation from a limit cycle. The chaotic dynamics of the order parameter is evidenced by calculating the maximal Lyapunov exponent. The presence of a chaotic domain is studied as a function of polariton-polariton interaction and the Josephson coupling between the condensates. At some values of the parameters the bifurcation route is more complex and the strange attractor can coexist with the stable fixed-point lasing. We also investigate how the chaotic dynamics is reflected in the light emission spectrum from the microcavity.


title: Parallel PERM

abstract: We develop and implement a parallel flatPERM algorithm \cite{G97,PK04} with mutually interacting parallel flatPERM sequences and use it to sample self-avoiding walks in 2 and 3 dimensions. Our data show that the parallel implementation accelerates the convergence of the flatPERM algorithm. Moreover, increasing the number of interacting flatPERM sequences (rather than running longer simulations) improves the rate of convergence. This suggests that a more efficient implementation of flatPERM will be a massively parallel implementation, rather than long simulations of one, or a few parallel sequences. We also use the algorithm to estimate the growth constant of the self-avoiding walk in two and in three dimensions using simulations over 12 parallel sequences. Our best results are \[ \mu_d = \cases{ 2.6381585(1), & \hbox{if $d=2$}; \cr 4.684039(1), & \hbox{if $d=3$}. } \]


title: A Novel Evolution Strategy with Directional Gaussian Smoothing for Blackbox Optimization

abstract: We propose an improved evolution strategy (ES) using a novel nonlocal gradient operator for high-dimensional black-box optimization. Standard ES methods with $d$-dimensional Gaussian smoothing suffer from the curse of dimensionality due to the high variance of Monte Carlo (MC) based gradient estimators. To control the variance, Gaussian smoothing is usually limited in a small region, so existing ES methods lack nonlocal exploration ability required for escaping from local minima. We develop a nonlocal gradient operator with directional Gaussian smoothing (DGS) to address this challenge. The DGS conducts 1D nonlocal explorations along $d$ orthogonal directions in $\mathbb{R}^d$, each of which defines a nonlocal directional derivative as a 1D integral. We then use Gauss-Hermite quadrature, instead of MC sampling, to estimate the $d$ 1D integrals to ensure high accuracy (i.e., small variance). Our method enables effective nonlocal exploration to facilitate the global search in high-dimensional optimization. We demonstrate the superior performance of our method in three sets of examples, including benchmark functions for global optimization, and real-world science and engineering applications.


title: Bounds on the Information Divergence for Hypergeometric Distributions

abstract: The hypergeometric distributions have many important applications, but they have not had sufficient attention in information theory. Hypergeometric distributions can be approximated by binomial distributions or Poisson distributions. In this paper we present upper and lower bounds on information divergence. These bounds are important for statistical testing and a better understanding of the notion of exchange-ability.


title: Superconductivity without insulating states in twisted bilayer graphene stabilized by monolayer WSe$_2$

abstract: Magic-angle twisted bilayer graphene (TBG), with rotational misalignment close to 1.1$^\circ$, features isolated flat electronic bands that host a rich phase diagram of correlated insulating, superconducting, ferromagnetic, and topological phases. The origins of the correlated insulators and superconductivity, and the interplay between them, are particularly elusive. Both states have been previously observed only for angles within $\pm0.1^\circ$ from the magic-angle value and occur in adjacent or overlapping electron density ranges; nevertheless, it is still unclear how the two states are related. Beyond the twist angle and strain, the dependence of the TBG phase diagram on the alignment and thickness of insulating hexagonal boron nitride (hBN) used to encapsulate the graphene sheets indicates the importance of the microscopic dielectric environment. Here we show that adding an insulating tungsten-diselenide (WSe$_2$) monolayer between hBN and TBG stabilizes superconductivity at twist angles much smaller than the established magic-angle value. For the smallest angle of $\theta$ = 0.79$^\circ$, we still observe clear superconducting signatures, despite the complete absence of the correlated insulating states and vanishing gaps between the dispersive and flat bands. These observations demonstrate that, even though electron correlations may be important, superconductivity in TBG can exist even when TBG exhibits metallic behaviour across the whole range of electron density. Finite-magnetic-field measurements further reveal breaking of the four-fold spin-valley symmetry in the system, consistent with large spin-orbit coupling induced in TBG via proximity to WSe$_2$. Our results highlight the importance of symmetry breaking effects in stabilizing electronic states in TBG and open new avenues for engineering quantum phases in moir\'e systems.


title: List Decodable Subspace Recovery

abstract: Learning from data in the presence of outliers is a fundamental problem in statistics. In this work, we study robust statistics in the presence of overwhelming outliers for the fundamental problem of subspace recovery. Given a dataset where an $\alpha$ fraction (less than half) of the data is distributed uniformly in an unknown $k$ dimensional subspace in $d$ dimensions, and with no additional assumptions on the remaining data, the goal is to recover a succinct list of $O(\frac{1}{\alpha})$ subspaces one of which is nontrivially correlated with the planted subspace. We provide the first polynomial time algorithm for the 'list decodable subspace recovery' problem, and subsume it under a more general framework of list decoding over distributions that are "certifiably resilient" capturing state of the art results for list decodable mean estimation and regression.


title: Deep Underground Neutrino Experiment (DUNE), Far Detector Technical Design Report, Volume II: DUNE Physics

abstract: The preponderance of matter over antimatter in the early universe, the dynamics of the supernovae that produced the heavy elements necessary for life, and whether protons eventually decay -- these mysteries at the forefront of particle physics and astrophysics are key to understanding the early evolution of our universe, its current state, and its eventual fate. DUNE is an international world-class experiment dedicated to addressing these questions as it searches for leptonic charge-parity symmetry violation, stands ready to capture supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. The DUNE far detector technical design report (TDR) describes the DUNE physics program and the technical designs of the single- and dual-phase DUNE liquid argon TPC far detector modules. Volume II of this TDR, DUNE Physics, describes the array of identified scientific opportunities and key goals. Crucially, we also report our best current understanding of the capability of DUNE to realize these goals, along with the detailed arguments and investigations on which this understanding is based. This TDR volume documents the scientific basis underlying the conception and design of the LBNF/DUNE experimental configurations. As a result, the description of DUNE's experimental capabilities constitutes the bulk of the document. Key linkages between requirements for successful execution of the physics program and primary specifications of the experimental configurations are drawn and summarized. This document also serves a wider purpose as a statement on the scientific potential of DUNE as a central component within a global program of frontier theoretical and experimental particle physics research. Thus, the presentation also aims to serve as a resource for the particle physics community at large.


title: Experimental characterization of a section of a spherically imploding plasma liner formed by merging hypersonic plasma jets

abstract: We report experimental results on merging of hypersonic plasma jets, which is the fundamental building block for forming spherically imploding plasma liners as a potential standoff compression driver for mangeto-inertial fusion. Jets are formed and launched by contoured-gap coaxial plasma guns mounted at the six spherical chamber. First, from experiments with two and three merging jets of four different species (N, Ar, Kr, Xe), we show that (1) density spatial non-uniformities can be large (with electron-density jumps ranging from 2.9 for N to 6.6 for Xe) when shocks form upon jet merging, but smaller (density jumps <2) when shocks do not form; (2) jet impurities (20% Ti in these experiments) can increase the level of density spatial non-uniformity by increasing the collisionality of jet merging; and (3) the liner Mach number can remain high (>10), as required for plasma liners to be an effective compression driver. Second, from experiments with six and seven merging jets using Ar, we present results with improved jet-to-jet balance of <2% across jets, including (1) evidence of substantially increased balance in the jet merging and symmetry of the liner structure, and (2) potentially favorable changes in the jet-merging morphology with the addition of the seventh jet. For both experiments, we present comparison between experimental and synthetic data from three-dimensional hydrodynamic codes.


title: Bayesian Methods for the Analysis of Early-Phase Oncology Basket Trials with Information Borrowing across Cancer Types

abstract: Research in oncology has changed the focus from histological properties of tumors in a specific organ to a specific genomic aberration potentially shared by multiple cancer types. This motivates the basket trial, which assesses the efficacy of treatment simultaneously on multiple cancer types that have a common aberration. Although the assumption of homogeneous treatment effects seems reasonable given the shared aberration, in reality, the treatment effect may vary by cancer type, and potentially only a subgroup of the cancer types respond to the treatment. Various approaches have been proposed to increase the trial power by borrowing information across cancer types, which, however, tend to inflate the type I error rate. In this paper, we review some representative Bayesian information borrowing methods for the analysis of early-phase basket trials. We then propose a novel method called the Bayesian hierarchical model with a correlated prior (CBHM), which conducts more flexible borrowing across cancer types according to sample similarity. We did simulation studies to compare CBHM with independent analysis and three information borrowing approaches: the conventional Bayesian hierarchical model, the EXNEX approach and Liu's two-stage approach. Simulation results show that all information borrowing approaches substantially improve the power of independent analysis if a large proportion of the cancer types truly respond to the treatment. Our proposed CBHM approach shows an advantage over the existing information borrowing approaches, with a power similar to that of EXNEX or Liu's approach, but the potential to provide substantially better control of type I error rate.


title: Deep Underground Neutrino Experiment (DUNE), Far Detector Technical Design Report, Volume III: DUNE Far Detector Technical Coordination

abstract: The preponderance of matter over antimatter in the early universe, the dynamics of the supernovae that produced the heavy elements necessary for life, and whether protons eventually decay -- these mysteries at the forefront of particle physics and astrophysics are key to understanding the early evolution of our universe, its current state, and its eventual fate. The Deep Underground Neutrino Experiment (DUNE) is an international world-class experiment dedicated to addressing these questions as it searches for leptonic charge-parity symmetry violation, stands ready to capture supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. The DUNE far detector technical design report (TDR) describes the DUNE physics program and the technical designs of the single- and dual-phase DUNE liquid argon TPC far detector modules. Volume III of this TDR describes how the activities required to design, construct, fabricate, install, and commission the DUNE far detector modules are organized and managed. This volume details the organizational structures that will carry out and/or oversee the planned far detector activities safely, successfully, on time, and on budget. It presents overviews of the facilities, supporting infrastructure, and detectors for context, and it outlines the project-related functions and methodologies used by the DUNE technical coordination organization, focusing on the areas of integration engineering, technical reviews, quality assurance and control, and safety oversight. Because of its more advanced stage of development, functional examples presented in this volume focus primarily on the single-phase (SP) detector module.


title: Blind Source Separation for NMR Spectra with Negative Intensity

abstract: NMR spectral datasets, especially in systems with limited samples, can be difficult to interpret if they contain multiple chemical components (phases, polymorphs, molecules, crystals, glasses, etc...) and the possibility of overlapping resonances. In this paper, we benchmark several blind source separation techniques for analysis of NMR spectral datasets containing negative intensity. For benchmarking purposes, we generated a large synthetic datasbase of quadrupolar solid-state NMR-like spectra that model spin-lattice T1 relaxation or nutation tip/flip angle experiments. Our benchmarking approach focused exclusively on the ability of blind source separation techniques to reproduce the spectra of the underlying pure components. In general, we find that FastICA (Fast Independent Component Analysis), SIMPLISMA (SIMPLe-to-use-Interactive Self-modeling Mixture Analysis), and NNMF (Non-Negative Matrix Factorization) are top-performing techniques. We demonstrate that dataset normalization approaches prior to blind source separation do not considerably improve outcomes. Within the range of noise levels studied, we did not find drastic changes to the ranking of techniques. The accuracy of FastICA and SIMPLISMA degrades quickly if excess (unreal) pure components are predicted. Our results indicate poor performance of SVD (Singular Value Decomposition) methods, and we propose alternative techniques for matrix initialization. The benchmarked techniques are also applied to real solid state NMR datasets. In general, the recommendations from the synthetic datasets agree with the recommendations and results from the real data analysis. The discussion provides some additional recommendations for spectroscopists applying blind source separation to NMR datasets, and for future benchmark studies.


title: Deep Underground Neutrino Experiment (DUNE), Far Detector Technical Design Report, Volume IV: Far Detector Single-phase Technology

abstract: The preponderance of matter over antimatter in the early universe, the dynamics of the supernovae that produced the heavy elements necessary for life, and whether protons eventually decay -- these mysteries at the forefront of particle physics and astrophysics are key to understanding the early evolution of our universe, its current state, and its eventual fate. DUNE is an international world-class experiment dedicated to addressing these questions as it searches for leptonic charge-parity symmetry violation, stands ready to capture supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. Central to achieving DUNE's physics program is a far detector that combines the many tens-of-kiloton fiducial mass necessary for rare event searches with sub-centimeter spatial resolution in its ability to image those events, allowing identification of the physics signatures among the numerous backgrounds. In the single-phase liquid argon time-projection chamber (LArTPC) technology, ionization charges drift horizontally in the liquid argon under the influence of an electric field towards a vertical anode, where they are read out with fine granularity. A photon detection system supplements the TPC, directly enhancing physics capabilities for all three DUNE physics drivers and opening up prospects for further physics explorations. The DUNE far detector technical design report (TDR) describes the DUNE physics program and the technical designs of the single- and dual-phase DUNE liquid argon TPC far detector modules. Volume IV presents an overview of the basic operating principles of a single-phase LArTPC, followed by a description of the DUNE implementation. Each of the subsystems is described in detail, connecting the high-level design requirements and decisions to the overriding physics goals of DUNE.


title: Probability Distribution of Magnetic Field Strengths through the Cyclotron Lines in High-Mass X-ray Binaries

abstract: The study of variation of measured cyclotron lines is of fundamental importance to understand the physics of the accretion process in magnetized neutron star systems. We investigate the magnetic field formation, evolution and distribution for several High- Mass X-ray Binaries (HMXBs). We focus our attention on the cyclotron lines that have been detected in HMXB classes in their X-ray spectra. As has been correctly pointed out, several sources show variation in cyclotron lines, this can result due to the effect of accretion dynamics, and hence that would reflect the magnetic field characteristics. Besides, the difference in time scales of variation of accretion rate and different type of companion can be used to distinguish between magnetized neutron stars.


title: Far-forward neutrinos at the Large Hadron Collider

abstract: We present a new calculation of the energy distribution of high-energy neutrinos from the decay of charm and bottom hadrons produced at the Large Hadron Collider (LHC). In the kinematical region of very forward rapidities, heavy-flavor production and decay is a source of tau neutrinos that leads to thousands of { charged-current} tau neutrino events in a 1 m long, 1 m radius lead neutrino detector at a distance of 480 m from the interaction region. In our computation, next-to-leading order QCD radiative corrections are accounted for in the production cross-sections. Non-perturbative intrinsic-$k_T$ effects are approximated by a simple phenomenological model introducing a Gaussian $k_T$-smearing of the parton distribution functions, which might also mimic perturbative effects due to multiple initial-state soft-gluon emissions. The transition from partonic to hadronic states is described by phenomenological fragmentation functions. To study the effect of various input parameters, theoretical predictions for $D_s^\pm$ production are compared with LHCb data on double-differential cross-sections in transverse momentum and rapidity. The uncertainties related to the choice of the input parameter values, ultimately affecting the predictions of the tau neutrino event distributions, are discussed. We consider a 3+1 neutrino mixing scenario to illustrate the potential for a neutrino experiment to constrain the 3+1 parameter space using tau neutrinos and antineutrinos. We find large theoretical uncertainties in the predictions of the neutrino fluxes in the far-forward region. Untangling the effects of tau neutrino oscillations into sterile neutrinos and distinguishing a 3+1 scenario from the standard scenario with three active neutrino flavours, will be challenging due to the large theoretical uncertainties from QCD.


title: Midi-review: Status of weak scale supersymmetry after LHC Run 2 and ton-scale noble liquid WIMP searches

abstract: While LHC has discovered a very Standard Model-like Higgs boson of mass m_h~ 125 GeV, no solid signal for physics beyond the Standard Model has emerged so far at LHC or at WIMP seach experiments. For the case of weak scale supersymmetry (SUSY), LHC has found rather generally that gluinos are beyond about 2.2 TeV whilst top squark must lie beyond 1.1 TeV. These limits contradict older simplistic notions of naturalness that emerged in the 1980s-1990s, leading to the rather pessimistic view that SUSY is now excluded except for perhaps some remaining narrow corners of parameter space. Yet, this picture ignores several important developments in SUSY/string theory that emerged in the 21st century: 1. the emergence of the string theory landscape and its solution to the cosmological constant problem, 2. a more nuanced view of naturalness including the notion of "stringy naturalness", 3. the emergence of anomaly-free discrete R-symmetries and their connection to R-parity, Peccei-Quinn symmetry, the SUSY mu problem and proton decay and 4. the importance of including a solution to the strong CP problem. Rather general considerations from the string theory landscape favor large values of soft terms, subject to the vacuum selection criteria that electroweak symmetry is properly broken (no CCB minima) and the resulting magnitude of the weak scale is not too far from our measured value. Then stringy naturalness predicts a Higgs mass m_h~ 125 GeV whilst sparticle masses are typically lifted beyond present LHC bounds. In light of these refinements in theory perspective confronted by LHC and dark matter search results, we review the most likely LHC, ILC and dark matter signatures that are expected to arise from weak scale SUSY as we understand it today.


title: FiniteNet: A Fully Convolutional LSTM Network Architecture for Time-Dependent Partial Differential Equations

abstract: In this work, we present a machine learning approach for reducing the error when numerically solving time-dependent partial differential equations (PDE). We use a fully convolutional LSTM network to exploit the spatiotemporal dynamics of PDEs. The neural network serves to enhance finite-difference and finite-volume methods (FDM/FVM) that are commonly used to solve PDEs, allowing us to maintain guarantees on the order of convergence of our method. We train the network on simulation data, and show that our network can reduce error by a factor of 2 to 3 compared to the baseline algorithms. We demonstrate our method on three PDEs that each feature qualitatively different dynamics. We look at the linear advection equation, which propagates its initial conditions at a constant speed, the inviscid Burgers' equation, which develops shockwaves, and the Kuramoto-Sivashinsky (KS) equation, which is chaotic.

categories: cs.LG cs.NA math.NA physics.comp-ph physics.flu-dyn stat.ML


title: Quantum teleportation with hybrid entangled resources prepared from heralded quantum states

abstract: In this work we propose the generation of a hybrid entangled resource (HER) and its further application in a quantum teleportation scheme from an experimentally feasible point of view. The source for HER preparation is based on the four wave mixing process in a photonic crystal fiber, from which one party of its output bipartite state is used to herald a single photon or a single photon added coherent state. From the heralded state and linear optics the HER is created. In the proposed teleportation protocol Bob uses the HER to teleport qubits with different spectral properties. Bob makes a Bell measurement in the single photon basis and characterizes the scheme with its average quantum teleportation fidelity. Fidelities close to one are expected for qubits in a wide spectral range. The work also includes a discussion about the fidelity dependence on the geometrical properties of the medium through which the HER is generated. An important remark is that no spectral filtering is employed in the heralding process, which emphasizes the feasibility of this scheme without compromising photon flux.


title: Safe Wasserstein Constrained Deep Q-Learning

abstract: This paper presents a distributionally robust Q-Learning algorithm (DrQ) which leverages Wasserstein ambiguity sets to provide idealistic probabilistic out-of-sample safety guarantees during online learning. First, we follow past work by separating the constraint functions from the principal objective to create a hierarchy of machines which estimate the feasible state-action space within the constrained Markov decision process (CMDP). DrQ works within this framework by augmenting constraint costs with tightening offset variables obtained through Wasserstein distributionally robust optimization (DRO). These offset variables correspond to worst-case distributions of modeling error characterized by the TD-errors of the constraint Q-functions. This procedure allows us to safely approach the nominal constraint boundaries. Using a case study of lithium-ion battery fast charging, we explore how idealistic safety guarantees translate to generally improved safety relative to conventional methods.


title: Discrete embedded modes in the continuum in 2D lattices

abstract: We study the problem of constructing bulk and surface embedded modes (EMs) inside the quasi-continuum band of a square lattice, using a potential engineering approach \`a la Wigner and von Neumann. Building on previous results for the one-dimensional (1D) lattice, and making use of separability, we produce examples of two-dimensional envelope functions and the two-dimensional (2D) potentials that produce them. The 2D embedded mode decays like a stretched exponential, with a supporting potential that decays as a power law. The separability process can cause that a 1D impurity state (outside the 1D band) can give rise to a 2D embedded mode (inside the band). The embedded mode survives the addition of random perturbations of the potential; however, this process introduces other localized modes inside the band, and causes a general tendency towards localization of the perturbed modes.


title: Certified Robustness to Label-Flipping Attacks via Randomized Smoothing

abstract: Machine learning algorithms are known to be susceptible to data poisoning attacks, where an adversary manipulates the training data to degrade performance of the resulting classifier. In this work, we present a unifying view of randomized smoothing over arbitrary functions, and we leverage this novel characterization to propose a new strategy for building classifiers that are pointwise-certifiably robust to general data poisoning attacks. As a specific instantiation, we utilize our framework to build linear classifiers that are robust to a strong variant of label flipping, where each test example is targeted independently. In other words, for each test point, our classifier includes a certification that its prediction would be the same had some number of training labels been changed adversarially. Randomized smoothing has previously been used to guarantee---with high probability---test-time robustness to adversarial manipulation of the input to a classifier; we derive a variant which provides a deterministic, analytical bound, sidestepping the probabilistic certificates that traditionally result from the sampling subprocedure. Further, we obtain these certified bounds with minimal additional runtime complexity over standard classification and no assumptions on the train or test distributions. We generalize our results to the multi-class case, providing the first multi-class classification algorithm that is certifiably robust to label-flipping attacks.


title: Generalized metric spaces. Relations with graphs, ordered sets and automata : A survey

abstract: In this survey we present a generalization of the notion of metric space and some applications to discrete structures as graphs, ordered sets and transition systems. Results in that direction started in the middle eighties based on the impulse given by Quilliot (1983). Graphs and ordered sets were considered as kind of metric spaces, where - instead of real numbers - the values of the distance functions $d$ belong to an ordered semigroup equipped with an involution. In this frame, maps preserving graphs or posets are exactly the nonexpansive mappings (that is the maps $f$ such that $d(f(x),f(y))\leq d(x,y)$, for all $x,y$). It was observed that many known results on retractions and fixed point property for classical metric spaces (whose morphisms are the nonexpansive mappings) are also valid for these spaces. For example, the characterization of absolute retracts, by Aronszajn and Panitchpakdi (1956), the construction of the injective envelope by Isbell (1965) and the fixed point theorem of Sine and Soardi (1979) translate into the Banaschewski-Bruns theorem (1967), the MacNeille completion of a poset (1933) and the famous Tarski fixed point theorem (1955). This prompted an analysis of several classes of discrete structures from a metric point of view. In this paper, we report the results obtained over the years with a particular emphasis on the fixed point property.


title: Multiple Dirac Nodes and Symmetry Protected Dirac Nodal Line in Orthorhombic $\alpha$-RhSi

abstract: Owing to their chiral cubic structure, exotic multifold topological excitations have been predicted and recently observed in transition metal silicides like $\beta$-RhSi. Herein, we report that the topological character of RhSi is also observed in its orthorhombic $\alpha$-phase which displays multiple types of Dirac nodes very close to the Fermi level ($\varepsilon_F$) with the near absence of topologically trivial carriers. We discuss the symmetry analysis, band connectivity along high-symmetry lines using group representations, the band structure, and the nature of the Dirac points and nodal lines occurring near $\varepsilon_F$. The de Haas-van Alphen effect (dHvA) indicates a Fermi surface in agreement with the calculations. We find an elliptically-shaped nodal line very close to $\varepsilon_F$ around and near the $S$-point on the $k_y-k_z$ plane that results from the intersection of two upside-down Dirac cones. The two Dirac points of the participating Kramers degenerate bands are only 5 meV apart, hence an accessible magnetic field might induce a crossing between the spin-up partner of the upper-Dirac cone and the spin-down partner of the lower Dirac cone, possibly explaining the anomalies observed in the magnetic torque.