title: Efficient Solution of Boolean Satisfiability Problems with Digital MemComputing

abstract: Boolean satisfiability is a propositional logic problem of interest in multiple fields, e.g., physics, mathematics, and computer science. Beyond a field of research, instances of the SAT problem, as it is known, require efficient solution methods in a variety of applications. It is the decision problem of determining whether a Boolean formula has a satisfying assignment, believed to require exponentially growing time for an algorithm to solve for the worst-case instances. Yet, the efficient solution of many classes of Boolean formulae eludes even the most successful algorithms, not only for the worst-case scenarios, but also for typical-case instances. Here, we introduce a memory-assisted physical system (a digital memcomputing machine) that, when its non-linear ordinary differential equations are integrated numerically, shows evidence for polynomially-bounded scalability while solving "hard" planted-solution instances of SAT, known to require exponential time to solve in the typical case for both complete and incomplete algorithms. Furthermore, we analytically demonstrate that the physical system can efficiently solve the SAT problem in continuous time, without the need to introduce chaos or an exponentially growing energy. The efficiency of the simulations is related to the collective dynamical properties of the original physical system that persist in the numerical integration to robustly guide the solution search even in the presence of numerical errors. We anticipate our results to broaden research directions in physics-inspired computing paradigms ranging from theory to application, from simulation to hardware implementation.


title: Many-body Majorana-like zero modes without gauge symmetry breaking

abstract: Topological superconductors represent one of the key hosts of Majorana-based topological quantum computing. Typical scenarios for one-dimensional topological superconductivity assume a broken gauge symmetry associated to a superconducting state. However, no interacting one-dimensional many-body system is known to spontaneously break gauge symmetries. Here, we show that zero modes emerge in a many-body system without gauge symmetry breaking and in the absence of superconducting order. In particular, we demonstrate that Majorana zero modes of the symmetry-broken superconducting state are continuously connected to these zero-mode excitations, demonstrating that zero-bias anomalies may emerge in the absence of gauge symmetry breaking. We demonstrate that these many-body zero modes share the robustness features of the Majorana zero modes of symmetry-broken topological superconductors. We introduce a bosonization formalism to analyze these excitations and show that a ground state analogous to a topological superconducting state can be analytically found in a certain limit. Our results demonstrate that robust Majorana-like zero modes may appear in a many-body systems without gauge symmetry breaking, thus introducing a family of protected excitations with no single-particle analogs.


title: Stress Testing Method for Scenario Based Testing of Automated Driving Systems

abstract: Classical approaches and procedures for testing of automated vehicles of SAE levels 1 and 2 were based on defined scenarios with specific maneuvers, depending on the function under test. For automated driving systems (ADS) of SAE level 3+, the scenario space is infinite and calling for virtual testing and verification. However, even in simulation, the generation of safety-relevant scenarios for ADS is expensive and time-consuming. This leads to a demand for stochastic and realistic traffic simulation. Therefore, microscopic traffic flow simulation models (TFSM) are becoming a crucial part of scenario-based testing of ADS. In this paper, a co-simulation between the multi-body simulation software IPG CarMaker and the microscopic traffic flow simulation software (TFSS) PTV Vissim is used. Although the TFSS could provide realistic and stochastic behavior of the traffic participants, safety-critical scenarios (SCS) occur rarely. In order to avoid this, a novel Stress Testing Method (STM) is introduced. With this method, traffic participants are manipulated via external driver DLL interface from PTV Vissim in the vicinity of the vehicle under test in order to provoke defined critical maneuvers derived from statistical accident data on highways in Austria. These external driver models imitate human driving errors, resulting in an increase of safety-critical scenarios. As a result, the presented STM method contributes to an increase of safety-relevant scenarios for verification, testing and assessment of ADS.


title: Gelfand numbers of embeddings of Schatten classes

abstract: Let $0<p,q\leq \infty$ and denote by $\mathcal{S}_p^N$ and $\mathcal{S}_q^N$ the corresponding Schatten classes of real $N\times N$ matrices. We study the Gelfand numbers of natural identities $\mathcal{S}_p^N\hookrightarrow \mathcal{S}_q^N$ between Schatten classes and prove asymptotically sharp bounds up to constants only depending on $p$ and $q$. This extends classical results for finite-dimensional $\ell_p$ sequence spaces by E. Gluskin to the non-commutative setting and complements bounds previously obtained by B. Carl and A. Defant, A. Hinrichs and C. Michels, and J. Ch\'avez-Dom\'inguez and D. Kutzarova.


title: Assessing congressional districting in Maine and New Hampshire

abstract: We use voting precinct and election data to analyze the political geography of New Hampshire and Maine. We find that the location of dividing line between Congressional districts in both states are significantly different than what we would expect, which we argue is likely due to incumbent gerrymandering. We also discuss the limitations of classical fairness measures for plans with only two districts.


title: Symmetry-protected topological phase transitions and robust chiral order on a tunable zigzag lattice

abstract: Symmetry fractionalization, generating a large amount of symmetry-protected topological phases, provides scenarios for continuous phase transitions different from spontaneous symmetry breaking. However, it is hard to detect these symmetry-protected topological phase transitions experimentally. Motivated by the recent development of highly tunable ultracold polar molecules, we show that the setup in a zigzag optical lattice of this system provides a perfect platform to realize symmetry-protected topological phase transitions. By using infinite time-evolving block decimation, we obtain the phase diagram in a large parameter regions and find another scheme to realize the long-sought vector chiral phase, which is robust from quantum fluctuations. We discuss the existence of the chiral phase by an effective field analysis.


title: A partition-based similarity for classification distributions

abstract: Herein we define a measure of similarity between classification distributions that is both principled from the perspective of statistical pattern recognition and useful from the perspective of machine learning practitioners. In particular, we propose a novel similarity on classification distributions, dubbed task similarity, that quantifies how an optimally-transformed optimal representation for a source distribution performs when applied to inference related to a target distribution. The definition of task similarity allows for natural definitions of adversarial and orthogonal distributions. We highlight limiting properties of representations induced by (universally) consistent decision rules and demonstrate in simulation that an empirical estimate of task similarity is a function of the decision rule deployed for inference. We demonstrate that for a given target distribution, both transfer efficiency and semantic similarity of candidate source distributions correlate with empirical task similarity.


title: MARTY, a new C++ framework for automated symbolic calculations in Beyond the Standard Model physics

abstract: Theoretical calculations Beyond the Standard Model (BSM) constitute a challenge for high energy physicists, but are necessary when searching for New Physics. The predictions of a BSM scenario need to be compared with experimental data and the Standard Model values in order to identify the model that fits better what we observe in particle colliders. BSM predictions require very involved and error prone calculations of amplitudes, cross-sections and Wilson coefficients. Calculations at the one-loop level are often necessary for these quantities since some phenomenologically important processes may not occur at tree-level, such as Flavor Changing Neutral Currents (FCNC) in flavor physics that vanish at the tree-level. One-loop calculations have to be done analytically which is very time consuming and in practice rarely done for general BSM models. Here we present MARTY, a public and open-source C++ code. MARTY is the very first independent program able to calculate automatically amplitudes, squared amplitudes and Wilson coefficients at the one-loop level for general BSM models. This type of calculations requires a computer algebra system and could only be done, up to now, using Mathematica, a commercial and closed software for symbolic manipulations. MARTY does not rely on Mathematica since it re-implements its own computer algebra system also in C++, called CSL. MARTY will considerably ease BSM studies as by automating the analytical calculations the main task of the user would be the model building part. Once interfaced with other phenomenological codes in particle physics, MARTY will be incredibly efficient to make detailed predictions for general BSM models automatically.


title: On the number of binary quadratic forms having discriminant $1-4p$, $p$ prime

abstract: In this paper we obtain an asymptotic formula for the number of $\operatorname{SL}_2(\mathbb{Z})$-equivalence classes of positive definite binary quadratic forms over $\mathbb{Z}$ having bounded discriminant $\Delta = 1-4p$, with $p$ a prime. This extends work of the first named author and has an application in counting simple $(4a+1)$-knots of genus one.


title: Benchmarking magnetizabilities with recent density functionals

abstract: We have assessed the accuracy for magnetic properties of a set of 51 density functional approximations, including both recently published as well as already established functionals. The accuracy assessment considers a series of 27 small molecules and is based on comparing the predicted magnetizabilities to literature reference values calculated using coupled cluster theory with full singles and doubles and perturbative triples [CCSD(T)] employing large basis sets. The most accurate magnetizabilities, defined as the smallest mean absolute error, were obtained with the BHandHLYP functional. Three of the six studied Berkeley functionals and the three range-separated Florida functionals also yield accurate magnetizabilities. Also some older functionals like CAM-B3LYP, KT1, BHLYP (BHandH), B3LYP and PBE0 perform rather well. In contrast, unsatisfactory performance was generally obtained with Minnesota functionals, which are therefore not recommended for calculations of magnetically induced current density susceptibilities, and related magnetic properties such as magnetizabilities and nuclear magnetic shieldings. We also demonstrate that magnetizabilities can be calculated by numerical integration of the magnetizability density; we have implemented this approach as a new feature in the gauge-including magnetically induced current method (GIMIC). Magnetizabilities can be calculated from magnetically induced current density susceptibilities within this approach even when analytical approaches for magnetizabilities as the second derivative of the energy have not been implemented. The magnetizability density can also be visualized, providing additional information that is not otherwise easily accessible on the spatial origin of the magnetizabilities.


title: Collective authorship in Ukrainian science: marginal effect or new phenomenon?

abstract: One of the features of modern science is the formation of stable large collaborations of researchers working together within the projects that require the concentration of huge financial and human resources. Results of such common work are published in scientific papers by large co-authorship teams that include sometimes thousands of names. The goal of this work is to study the influence of such publications on the values of scientometric indicators calculated for individuals, research groups and science of Ukraine in general. Bibliometric data related to Ukraine, some academic institutions and selected individual researchers were collected from Scopus database and used for our study. It is demonstrated that while the relative share of publications by collective authors is comparatively small, their presence in a general pool can lead to statistically significant effects. The obtained results clearly show that traditional quantitative approaches for research assessment should be changed in order to take into account this phenomenon. Keywords: collective authorship, scientometrics, group science, Ukraine.


title: A generalized Poincar\'e-Birkhoff theorem

abstract: We prove a generalization of the classical Poincar\'e--Birkhoff theorem for Liouville domains, in arbitrary even dimensions. This is inspired by the existence of global hypersurfaces of section for the spatial case of the restricted three-body problem (as proved by the authors in arXiv:2011.10386).


title: Grothendieck's standard conjecture of Lefschetz type over finite fields

abstract: Grothendieck's standard conjecture of Lefschetz type has two main forms: the weak form $C$ and the strong form $B$. The weak form is known for varieties over finite fields as a consequence of the proof of the Weil conjectures. This suggests that the strong form of the conjecture in the same setting may be the most accessible of the standard conjectures. Here, as an advertisement for the conjecture, we explain some of its remarkable consequences.


title: On the free energy density of factor models on biregular graphs

abstract: Let $h(0),h(1),\dots,h(k)$ be a symmetric concave sequence. For a $(d,k)$-biregular factor graph $G$ and $x\in \{0,1\}^V$, we define the Hamiltonian \[H_G(x)=\sum_{f\in F} h\left(\sum_{v\in \partial f} x_v\right),\] where $V$ is the set of variable nodes, $F$ is the set of factor nodes. We prove that if $(G_n)$ is a large girth sequence of $(d,k)$-biregular factor graphs, then the free energy density of $G_n$ converges. The limiting free energy density is given by the Bethe-approximation.


title: Integrable nonunitary open quantum circuits

abstract: We explicitly construct an integrable and strongly interacting dissipative quantum circuit via a trotterization of the Hubbard model with imaginary interaction strength. To prove integrability, we build an inhomogeneous transfer matrix, from which conserved superoperator charges can be derived, in particular, the circuit's dynamical generator. After showing the trace preservation and complete positivity of local maps, we reinterpret them as the Kraus representation of the local dynamics of free fermions with single-site dephasing. The integrability of the map is broken by adding interactions to the local coherent dynamics or by removing the dephasing. In particular, even circuits built from convex combinations of local free-fermion unitaries are nonintegrable. Moreover, the construction allows us to explicitly build circuits belonging to different non-Hermitian symmetry classes, which are characterized by the behavior under transposition instead of complex conjugation. We confirm all our analytical results by using complex spacing ratios to examine the spectral statistics of the dissipative circuits.

categories: cond-mat.stat-mech cond-mat.str-el math-ph math.MP nlin.SI quant-ph


title: Tunable proximity effects and topological superconductivity in ferromagnetic hybrid nanowires

abstract: Hybrid semiconducting nanowire devices combining epitaxial superconductor and ferromagnetic insulator layers have been recently explored experimentally as an alternative platform for topological superconductivity at zero applied magnetic field. In this proof-of-principle work we show that the topological regime can be reached in actual devices depending on some geometrical constraints. To this end, we perform numerical simulations of InAs wires in which we explicitly include the superconducting Al and magnetic EuS shells, as well as the interaction with the electrostatic environment at a self-consistent mean-field level. Our calculations show that both the magnetic and the superconducting proximity effects on the nanowire can be tuned by nearby gates thanks to their ability to move the wavefunction across the wire section. We find that the topological phase is achieved in significant portions of the phase diagram only in configurations where the Al and EuS layers overlap on some wire facet, due to the rather local direct induced spin polarization and the appearance of an extra indirect exchange field through the superconductor. While of obvious relevance for the explanation of recent experiments, tunable proximity effects are of interest in the broader field of superconducting spintronics.


title: Electronic properties of InAs/EuS/Al hybrid nanowires

abstract: We study the electronic properties of InAs/EuS/Al heterostructures as explored in a recent experiment [S. Vaitiekenas \emph{et al.}, Nat. Phys. (2020)], combining both spectroscopic results and microscopic device simulations. In particular, we use angle-resolved photoemission spectroscopy to investigate the band bending at the InAs/EuS interface. The resulting band offset value serves as an essential input to subsequent microscopic device simulations, allowing us to map the electronic wave function distribution. We conclude that the magnetic proximity effects at the Al/EuS as well as the InAs/EuS interfaces are both essential to achieve topological superconductivity at zero applied magnetic field. Mapping the topological phase diagram as a function of gate voltages and proximity-induced exchange couplings, we show that the ferromagnetic hybrid nanowire with overlapping Al and EuS layers can become a topological superconductor within realistic parameter regimes, and that the topological phase can be optimized by external gating. Our work highlights the need for a combined experimental and theoretical effort for faithful device simulation.


title: Pseudo-holomorphic dynamics in the restricted three-body problem

abstract: In this article, we identify the 5-dimensional analogue of the finite energy foliations introduced by Hofer--Wysocki--Zehnder for the study of 3-dimensional Reeb flows, and show that these exist for the spatial circular restricted three-body problem (SCR3BP) whenever the planar dynamics is convex. We introduce the notion of a fiberwise-recurrent point, which may be thought of as a symplectic version of the leafwise intersections introduced by Moser, and show that they exist in abundance for a perturbative regime in the SCR3BP. We then use this foliation to induce a Reeb flow on the standard 3-sphere, via the use of pseudo-holomorphic curves, to be understood as the best approximation of the given dynamics that preserves the foliation. We discuss examples, further geometric structures, and speculate on possible applications.


title: Usefulness of adaptive strategies in asymptotic quantum channel discrimination

abstract: Adaptiveness is a key principle in information processing including statistics and machine learning. We investigate the usefulness of adaptive methods in the framework of asymptotic binary hypothesis testing, when each hypothesis represents asymptotically many independent instances of a quantum channel, and the tests are based on using the unknown channel and observing outputs. Unlike the familiar setting of quantum states as hypotheses, there is a fundamental distinction between adaptive and non-adaptive strategies with respect to the channel uses, and we introduce a number of further variants of the discrimination tasks by imposing different restrictions on the test strategies. The following results are obtained: (1) We prove that for classical-quantum channels, adaptive and non-adaptive strategies lead to the same error exponents both in the symmetric (Chernoff) and asymmetric (Hoeffding, Stein) settings. (2) The first separation between adaptive and non-adaptive symmetric hypothesis testing exponents for quantum channels, which we derive from a general lower bound on the error probability for non-adaptive strategies; the concrete example we analyze is a pair of entanglement-breaking channels. (3)We prove, in some sense generalizing the previous statement, that for general channels adaptive strategies restricted to classical feed-forward and product state channel inputs are not superior in the asymptotic limit to non-adaptive product state strategies. (4) As an application of our findings, we address the discrimination power of an arbitrary quantum channel and show that adaptive strategies with classical feedback and no quantum memory at the input do not increase the discrimination power of the channel beyond non-adaptive tensor product input strategies.


title: Fundamental Reference AGN Monitoring Experiment (FRAMEx) I: Jumping Out of the Plane with the VLBA

abstract: We present the first results from the Fundamental Reference AGN Monitoring Experiment (FRAMEx), an observational campaign dedicated to understanding the physical processes that affect the apparent positions and morphologies of AGNs. In this work, we obtained simultaneous Swift X-ray Telescope (XRT) and Very Long Baseline Array (VLBA) radio observations for a snapshot campaign of 25 local AGNs that form a volume-complete sample with hard X-ray (14-195 keV) luminosities above $10^{42}$erg s$^{-1}$, out to a distance of 40 Mpc. Despite achieving an observation depth of $\sim20$ $\mu$Jy, we find that 16 of 25 AGNs in our sample are not detected with the VLBA on milli-arcsecond (sub-parsec) scales, and the corresponding core radio luminosity upper limits are systematically below predictions from the Fundamental Plane of black hole activity. Using archival Jansky Very Large Array (VLA) radio measurements, our sample jumps back onto the Fundamental Plane, suggesting that extended radio emission is responsible for the apparent correlation between radio emission, X-ray emission, and black hole mass. We suggest that this discrepancy is likely due to extra-nuclear radio emission produced via interactions between the AGN and host environment. We compare VLBA observations of AGNs to VLA observations of nearby Galactic black holes (GBHs) and we find a mass-independent correlation between radio and X-ray luminosities of black holes of $L_\mathrm{6cm}$/$L_\mathrm{2-10 keV}$ $\sim$ 10$^{-6}$, in line with predictions for coronal emission, but allowing for the possibility of truly radio silent AGNs.


title: Quantum algorithm for nonlinear differential equations

abstract: Quantum computers are known to provide an exponential advantage over classical computers for the solution of linear differential equations in high-dimensional spaces. Here, we present a quantum algorithm for the solution of nonlinear differential equations. The quantum algorithm provides an exponential advantage over classical algorithms for solving nonlinear differential equations. Potential applications include the Navier-Stokes equation, plasma hydrodynamics, epidemiology, and more.


title: Relative Lipschitzness in Extragradient Methods and a Direct Recipe for Acceleration

abstract: We show that standard extragradient methods (i.e. mirror prox and dual extrapolation) recover optimal accelerated rates for first-order minimization of smooth convex functions. To obtain this result we provide a fine-grained characterization of the convergence rates of extragradient methods for solving monotone variational inequalities in terms of a natural condition we call relative Lipschitzness. We further generalize this framework to handle local and randomized notions of relative Lipschitzness and thereby recover rates for box-constrained $\ell_\infty$ regression based on area convexity and complexity bounds achieved by accelerated (randomized) coordinate descent for smooth convex function minimization.


title: Toward simulating Superstring/M-theory on a quantum computer

abstract: We present a novel framework for simulating matrix models on a quantum computer. Supersymmetric matrix models have natural applications to superstring/M-theory and gravitational physics, in an appropriate limit of parameters. Furthermore, for certain states in the Berenstein-Maldacena-Nastase (BMN) matrix model, several supersymmetric quantum field theories dual to superstring/M-theory can be realized on a quantum device. Our prescription consists of four steps: regularization of the Hilbert space, adiabatic state preparation, simulation of real-time dynamics, and measurements. Regularization is performed for the BMN matrix model with the introduction of energy cut-off via the truncation in the Fock space. We use the Wan-Kim algorithm for fast digital adiabatic state preparation to prepare the low-energy eigenstates of this model as well as thermofield double state. Then, we provide an explicit construction for simulating real-time dynamics utilizing techniques of block-encoding, qubitization, and quantum signal processing. Lastly, we present a set of measurements and experiments that can be carried out on a quantum computer to further our understanding of superstring/M-theory beyond analytic results.


title: Generalized Carleson perturbations of elliptic operators and applications

abstract: We extend in two directions the notion of perturbations of Carleson type for the Dirichlet problem associated to an elliptic real second-order divergence-form (possibly degenerate, not necessarily symmetric) elliptic operator. First, in addition to the classical perturbations of Carleson type, that we call additive Carleson perturbations, we introduce scalar-multiplicative and antisymmetric Carleson perturbations, which both allow non-trivial differences at the boundary. Second, we consider domains which admit an elliptic PDE in a broad sense: we count as examples the 1-sided NTA (a.k.a. uniform) domains satisfying the capacity density condition, the 1-sided chord-arc domains, the domains with low-dimensional Ahlfors-David regular boundaries, and certain domains with mixed-dimensional boundaries; thus our methods provide a unified perspective on the Carleson perturbation theory of elliptic operators. Our proofs do not introduce sawtooth domains or the extrapolation method. We also present several applications to some Dahlberg-Kenig-Pipher operators, free-boundary problems, and we provide a new characterization of $A_{\infty}$ among elliptic measures.


title: Noncoherent Multiuser Chirp Spread Spectrum: Performance with Doppler and Asynchronism

abstract: In this paper, we investigate multi user chirp spread spectrum with noncoherent detection as a continuation of our work on coherent detection in [1]. We derive the analytical bit error ratio (BER) expression for binary chirp spread spectrum (BCSS) in the presence of multiple access interference (MAI) caused by correlation with other user signals because of either asynchronism or Doppler shifts, or both, and validate with simulations. To achieve this we analyze the signal cross correlations, and compare traditional linear chirps with our recently-proposed nonlinear chirps introduced in [1] and with other nonlinear chirps from the literature. In doing so we illustrate the superior performance of our new nonlinear chirp designs in these practical conditions, for the noncoherent counterpart of [1].


title: Quantum simulation of gauge theory via orbifold lattice

abstract: We propose a new framework for simulating $\text{U}(k)$ Yang-Mills theory on a universal quantum computer. This construction uses the orbifold lattice formulation proposed by Kaplan, Katz, and Unsal, who originally applied it to supersymmetric gauge theories. Our proposed approach yields a novel perspective on quantum simulation of quantum field theories, carrying certain advantages over the usual Kogut-Susskind formulation. We discuss the application of our constructions to computing static properties and real-time dynamics of Yang-Mills theories, from glueball measurements to AdS/CFT, making use of a variety of quantum information techniques including qubitization, quantum signal processing, Jordan-Lee-Preskill bounds, and shadow tomography. The generalizations to certain supersymmetric Yang-Mills theories appear to be straightforward, providing a path towards the quantum simulation of quantum gravity via holographic duality.


title: Diffusive limits of two-parameter ordered Chinese Restaurant Process up-down chains

abstract: We construct a two-parameter family of Feller diffusions on the set of open subsets of $(0,1)$ that arise as diffusive limits of two-parameter ordered Chinese Restaurant Process up-down chains. The diffusions we construct are natural ordered analogues of Petrov's two-parameter extension of Ethier and Kurtz's infinitely-many-neutral-alleles diffusion model. Recently, there has been significant interest in ordered analogues of the diffusions Petrov constructed. Existing methods for constructing such processes have been based on pathwise methods using marked L\'evy processes and an outstanding conjecture about these processes is that they are, in fact, the diffusive limit of the ordered Chinese Restaurant Process up-down chains that we consider here. We make progress on this conjecture by showing that the diffusive limit of the ordered Chinese Restaurant Process up-down chains exists. Moreover, our methods yield a simple, explicit description of the generator of the limiting processes on a core described in terms of quasisymmetric functions.


title: Distance between reproducing kernel Hilbert spaces and geometry of finite sets in the unit ball

abstract: In this paper we study the relationships between a reproducing kernel Hilbert space, its multiplier algebra, and the geometry of the point set on which they live. We introduce a variant of the Banach-Mazur distance suited for measuring the distance between reproducing kernel Hilbert spaces, that quantifies how far two spaces are from being isometrically isomorphic as reproducing kernel Hilbert spaces. We introduce an analogous distance for multiplier algebras, that quantifies how far two algebras are from being completely isometrically isomorphic. We show that, in the setting of finite dimensional quotients of the Drury-Arveson space, two spaces are "close" to one another if and only if their multiplier algebras are "close", and that this happens if and only if the underlying point-sets are "almost congruent", meaning that one of the sets is very close to an image of the other under a biholomorphic automorphism of the unit ball. These equivalences are obtained as corollaries of quantitative estimates that we prove.


title: A geometric view on Iwasawa theory

abstract: This article extends our study of the geometry of the $p$-adic eigencurve at a point defined by a weight $1$ cuspform $f$ irregular at $p$ and having complex multiplication, and the implications in Iwasawa and in Hida theories. The novel results include the determination of the Fourier coefficients of certain non-classical $p$-adic modular forms belonging to the generalized eigenspace of $f$, in terms of $p$-adic logarithms of algebraic numbers. We also compute the "mysterious" cross-ratios of the $p$-ordinary filtrations of the Hida families containing $f$.


title: Voltage tunable quantum dot array by patterned Ge-nanowire based metal-oxide-semiconductor (MOS) devices

abstract: Semiconductor quantum dots (QDs) are being regarded as the primary unit for a wide range of advanced and emerging technologies including electronics, optoelectronics, photovoltaics and biosensing applications as well as the domain of q-bits based quantum information processing. Such QDs are suitable for several novel device applications for their unique property of confining carriers 3-dimensionally creating discrete quantum states. However, the realization of such QDs in practice exhibits serious challenge regarding their fabrication in array with desired scalability and repeatability as well as control over the quantum states at room temperature. In this context, the current work reports the fabrication of an array of highly scaled Ge-nanowire (radius ~25 nm) based vertical metal-oxide-semiconductor devices that can operate as voltage tunable quantum dots at room temperature. The electrons in such nanowire experience a geometrical confinement in the radial direction, whereas, they can be confined axially by tuning the applied bias in order to manipulate the quantum states. Such quantum confinement of electrons has been confirmed from the step-like responses in the room temperature capacitance-voltage (C-V) characteristics at relatively low frequency (200 kHz). Each of such steps has observed to encompass convolution of the quantized states occupying ~6 electronic charges. The details of such carrier confinement are analyzed in the current work by theoretically modeling the device transport properties based on non-equilibrium Green's function (NEGF) formalism.


title: On regular and random two-dimensional packing of crosses

abstract: Packing problems, even of objects with regular geometries, are in general non-trivial. For few special shapes, the features of crystalline as well as random, irregular two-dimensional (2D) packings are known. The packing of 2D crosses does not yet belong to the category of solved problems. We demonstrate in experiments with crosses of different aspect ratios (arm width to length) which packing fractions are actually achieved by random packing, and we compare them to densest regular packing structures. We determine local correlations of the orientations and positions after ensembles of randomly placed crosses were compacted in the plane until they jam. Short-range orientational order is found over 2 to 3 cross lengths. Similarly, correlations in the spatial distributions of neighbors extend over 2 to 3 crosses. Apparently, there is no simple relation between the geometries of the crosses and peaks in the spatial correlation functions. Some features of the orientational correlations, however, are intuitively evident.


title: Avoiding baryonic feedback effects on neutrino mass measurements from CMB lensing

abstract: A measurement of the sum of neutrino masses is one of the main applications of upcoming measurements of gravitational lensing of the cosmic microwave background (CMB). This measurement can be confounded by modelling uncertainties related to so-called "baryonic effects" on the clustering of matter, arising from gas dynamics, star formation, and feedback from active galactic nuclei and supernovae. In particular, a wrong assumption about the form of baryonic effects on CMB lensing can bias a neutrino mass measurement by a significant fraction of the statistical uncertainty. In this paper, we investigate three methods for mitigating this bias: (1) restricting the use of small-scale CMB lensing information when constraining neutrino mass; (2) using an external tracer to remove the low-redshift contribution to a CMB lensing map; and (3) marginalizing over a parametric model for baryonic effects on large-scale structure. We test these methods using Fisher matrix forecasts for experiments resembling the Simons Observatory and CMB-S4, using a variety of recent hydrodynamical simulations to represent the range of possible baryonic effects, and using cosmic shear measured by the Rubin Observatory's LSST as the tracer in method (2). We find that a combination of (1) and (2), or (3) on its own, will be effective in reducing the bias induced by baryonic effects on a neutrino mass measurement to a negligible level, without a significant increase in the associated statistical uncertainty.


title: NuSTAR Survey of Obscured Swift/BAT-selected Active Galactic Nuclei: II. Median High-energy Cutoff in Seyfert II Hard X-ray Spectra

abstract: Broadband X-ray spectroscopy of the X-ray emission produced in the coronae of active galactic nuclei (AGN) can provide important insights into the physical conditions very close to their central supermassive black holes. The temperature of the Comptonizing plasma that forms the corona is manifested through a high-energy cutoff that has been difficult to directly constrain even in the brightest AGN because it requires high-quality data at energies above 10 keV. In this paper we present a large collection of coronal cutoff constraints for obscured AGN based on a sample of 130 AGN selected in the hard X-ray band with Swift/BAT and observed nearly simultaneously with NuSTAR and Swift/XRT. We find that under a reasonable set of assumptions regarding partial constraints the median cutoff is well constrained to 290$\pm$20 keV, where the uncertainty is statistical and given at the 68% confidence level. We investigate the sensitivity of this result to our assumptions and find that consideration of various known systematic uncertainties robustly places the median cutoff between 240 keV and 340 keV. The central 68% of the intrinsic cutoff distribution is found to be between about 140 keV and 500 keV, with estimated uncertainties of 20 keV and 100 keV, respectively. In comparison with the literature, we find no clear evidence that the cutoffs in obscured and unobscured AGN are substantially different. Our analysis highlights the importance of carefully considering partial and potentially degenerate constraints on the coronal high-energy cutoff in AGN.


title: Cosmological simulation in tides: power spectra and halo shape responses, and shape assembly bias

abstract: The well-developed separate universe technique enables accurate calibration of the response of any observable to an isotropic long-wavelength density fluctuation. The large-scale environment also hosts tidal modes that perturb all observables anisotropically. As in the separate universe, both the long tidal and density modes can be absorbed by an effective anisotropic background, on which the interaction and evolution of the short modes change accordingly. We further develop the tidal simulation method, including proper corrections to the second order Lagrangian perturbation theory (2LPT) to generate initial conditions of the simulations. We measure the linear tidal responses of the matter power spectrum, at high redshift from our modified 2LPT, and at low redshift from the tidal simulations. Our results agree qualitatively with previous works, but exhibit quantitative differences in both cases. We also measure the linear tidal response of the halo shapes, or the shape bias, and find its universal relation with the linear halo bias, for which we provide a fitting formula. Furthermore, analogous to the assembly bias, we study the secondary dependence of the shape bias, and discover for the first time dependence on halo concentration and axis ratio. Our results provide useful insights for studies of the intrinsic alignment as source of either contamination or information. These effects need to be correctly taken into account when one uses intrinsic alignments of galaxy shapes as a precision cosmological tool.


title: Sparse PCA: Algorithms, Adversarial Perturbations and Certificates

abstract: We study efficient algorithms for Sparse PCA in standard statistical models (spiked covariance in its Wishart form). Our goal is to achieve optimal recovery guarantees while being resilient to small perturbations. Despite a long history of prior works, including explicit studies of perturbation resilience, the best known algorithmic guarantees for Sparse PCA are fragile and break down under small adversarial perturbations. We observe a basic connection between perturbation resilience and \emph{certifying algorithms} that are based on certificates of upper bounds on sparse eigenvalues of random matrices. In contrast to other techniques, such certifying algorithms, including the brute-force maximum likelihood estimator, are automatically robust against small adversarial perturbation. We use this connection to obtain the first polynomial-time algorithms for this problem that are resilient against additive adversarial perturbations by obtaining new efficient certificates for upper bounds on sparse eigenvalues of random matrices. Our algorithms are based either on basic semidefinite programming or on its low-degree sum-of-squares strengthening depending on the parameter regimes. Their guarantees either match or approach the best known guarantees of \emph{fragile} algorithms in terms of sparsity of the unknown vector, number of samples and the ambient dimension. To complement our algorithmic results, we prove rigorous lower bounds matching the gap between fragile and robust polynomial-time algorithms in a natural computational model based on low-degree polynomials (closely related to the pseudo-calibration technique for sum-of-squares lower bounds) that is known to capture the best known guarantees for related statistical estimation problems. The combination of these results provides formal evidence of an inherent price to pay to achieve robustness.


title: The density of rational points on $\mathbb{P}^1$ with three stacky points

abstract: In this paper we consider the density of rational points on the "stacky" curve $\mathcal{X}(\mathbb{P}^1;0,2;1,2;\infty,2)$ which is $\mathbb{P}^1$ with three half points, with respect to the so-called Ellenberg-Satriano-Zuerick-Brown height. In particular, we prove a conjecture of Ellenberg.


title: Impacts of the physical data model on the forward inference of initial conditions from biased tracers

abstract: We investigate the impact of each ingredient in the employed physical data model on the Bayesian forward inference of initial conditions from biased tracers at the field level. Specifically, we use dark matter halos in a given cosmological simulation volume as tracers of the underlying matter density field. We study the effect of tracer density, grid resolution, gravity model, bias model and likelihood on the inferred initial conditions. We find that the cross-correlation coefficient between true and inferred phases reacts weakly to all ingredients above, and is well predicted by the theoretical expectation derived from a Gaussian model on a broad range of scales. The bias in the amplitude of the inferred initial conditions, on the other hand, depends strongly on the bias model and the likelihood. We conclude that the bias model and likelihood hold the key to an unbiased cosmological inference. Together they must keep the systematics -- which arise from the sub-grid physics that are marginalized over -- under control in order to obtain an unbiased inference.


title: The classical double copy in curved spacetimes: Perturbative Yang-Mills from the bi-adjoint scalar

abstract: We formulate a version of the double copy for classical fields in curved spacetimes. We provide a correspondence between perturbative solutions to the bi-adjoint scalar equations and those of the Yang-Mills equations in position space. At the linear level, we show that there exists a map between these solutions for maximally symmetric spacetime backgrounds, that provides every Yang-Mills solution by the action of an appropriate differential operator on a bi-adjoint scalar solution. Given the existence of a linearized map, we show that it is possible to cast the solutions of the Yang-Mills equations at arbitrary perturbation order in terms of the corresponding bi-adjoint scalar solutions. This all-order map is reminiscent of the flat space BCJ double copy, and works for any curved spacetime where the perturbative expansion holds. We show that these results have the right flat space limit, and that the correspondence is agnostic to the choice of gauge.


title: Strongly-Interacting Ultralight Millicharged Particles

abstract: We consider the implications of an ultra-light fermionic dark matter candidate that carries baryon number. This naturally arises if dark matter has a small charge under standard model baryon number whilst having an asymmetry equal and opposite to that in the visible universe. A prototypical model is a theory of dark baryons of a non-Abelian gauge group, i.e., a dark Quantum Chromo-Dynamics (QCD). For sub-eV dark baryon masses, the inner region of dark matter halos is naturally at 'nuclear density', allowing for the formation of exotic states of matter, akin to neutron stars. The Tremaine-Gunn lower bound on the mass of fermionic dark matter, i.e., the dark baryons, is violated by the strong short-range self-interactions, cooling via emission of light dark pions, and the Cooper pairing of dark quarks that occurs at densities that are high relative to the (ultra-low) dark QCD scale. We develop the astrophysics of these STrongly-interacting Ultra-light Millicharged Particles (STUMPs) utilizing the equation of state of dense quark matter, and find halo cores consistent with observations of dwarf galaxies. These cores are prevented from core-collapse by pressure of the 'neutron star', which suggests ultra-light dark QCD as a resolution to core-cusp problem of collisionless cold dark matter. The model is distinguished from ultra-light bosonic dark matter through direct detection and collider signatures, as well as by phenomena associated with superconductivity, such as Andreev reflection and superconducting vortices.


title: A quantum optics approach to photoinduced electron transfer in cavities

abstract: We study a simple model for photoinduced electron transfer reactions for the case of many donor-acceptor pairs that are collectively and homogeneously coupled to a photon mode of a cavity. We describe both coherent and dissipative collective effects resulting from this coupling within the framework of a quantum optics Lindblad master equation. We introduce a method to derive an effective rate equation for electron transfer, by adiabatically eliminating donor and acceptor states and the cavity mode. The resulting rate equation is valid both for weak and strong coupling to the cavity mode, and describes electronic transfer through both the cavity coupled bright states and the uncoupled dark states. We derive an analytic expression for the instantaneous electron transfer rate that depends non-trivially on the time-varying number of pairs in the ground state. We find that under proper resonance conditions, and in the presence of an incoherent drive, reaction rates can be enhanced by the cavity. This enhancement persists, and can even be largest, in the weak light-matter coupling regime. We discuss how the cavity effect is relevant for realistic experiments.


title: Search for Black Hole Merger Families

abstract: The origin, environment, and evolution of stellar-mass black hole binaries are still a mystery. One of the proposed binary formation mechanisms is manifest in dynamical interactions between multiple black holes. A resulting framework of these dynamical interactions is the so-called hierarchical triple merger scenario, which happens when three black holes become gravitationally bound, causing two successive black hole mergers to occur. In such successive mergers, the black holes involved are directly related to each other, and hence this channel can be directly tested from the properties of the detected binary black hole mergers. Here we present a search for hierarchical triple mergers among events within the GWTC-1 and GWTC-2 catalogs of LIGO/Virgo, the eccentric localization of GW190521 and those found by the IAS-Princeton group. The search includes improved statistical quantification that also accounts for black hole spins. We perform our analysis for different upper bounds on the mass distribution of first generation BHs. Our results demonstrate the importance of the mass distributions' properties for constraining the hierarchical merger scenario. We present the individually significant merger pairs. The search yields interesting candidate families and hints of its future impact.


title: Shadow economy and populism-risk and uncertainty factors for establishing low-carbon economy of Balkan countries-case study for Bulgaria

abstract: The main purpose of the current publication is to formulate a scenario model for analysis of the opportunities for low-carbon economy establishment in the countries with transition economies.The model studies risk factors such as shadow economy level and populism based on the implementation and development of Balkan countries economic policy and at the same time shows future climate changes tendencies and uncertainties of climate models.A transdisciplinary approach is implemented in the study. Climate change perception and understanding about low-carbon economy are examined through the public opinion and analysis of mass-media publications.The results of the research are important in order to clarify the multicultural divergences as a factor for risk and uncertainty in the implementation process of the policy for climate change.In this way geographical aspects of risk and uncertainty, which are not only related to the economic development of the relevant countries, could be brought out.


title: A stability-driven protocol for drug response interpretable prediction (staDRIP)

abstract: Modern cancer -omics and pharmacological data hold great promise in precision cancer medicine for developing individualized patient treatments. However, high heterogeneity and noise in such data pose challenges for predicting the response of cancer cell lines to therapeutic drugs accurately. As a result, arbitrary human judgment calls are rampant throughout the predictive modeling pipeline. In this work, we develop a transparent stability-driven pipeline for drug response interpretable predictions, or staDRIP, which builds upon the PCS framework for veridical data science (Yu and Kumbier, 2020) and mitigates the impact of human judgment calls. Here we use the PCS framework for the first time in cancer research to extract proteins and genes that are important in predicting the drug responses and stable across appropriate data and model perturbations. Out of the 24 most stable proteins we identified using data from the Cancer Cell Line Encyclopedia (CCLE), 18 have been associated with the drug response or identified as a known or possible drug target in previous literature, demonstrating the utility of our stability-driven pipeline for knowledge discovery in cancer drug response prediction modeling.


title: When optimal is not the best: cost effectiveness analysis for HPV epidemic models

abstract: This paper aims to evaluate the potential cost-effectiveness of healthcare interventions against human papillomavirus (HPV). For this, we consider a two-sex epidemic model for the transmission dynamics of HPV which includes screening, vaccination of adolescent boys and girls, and vaccination of sexually active adults. We first propose public health policies using constant control parameters and develop a cost-effectiveness analysis (CEA) to identify which intervention delivers the best effectiveness for the money invested. Secondly, we consider time-dependent control parameters and formulate an optimal control problem to obtain time-dependent versions of the interventions. As in the case of constant control parameters, we perform a CEA to investigate the cost-effectiveness of the time-dependent control interventions. Our findings suggest that females' vaccination, including adolescent girls and adult women, is the most cost-effective strategy. We also compare constant against the time-dependent healthcare interventions which are optimal in the sense that they minimize the objective functional of the optimal control problem. The results indicate that time-dependent controls are not always more cost-effective than constant controls.


title: Topological order and criticality in (2+1)D monitored random quantum circuits

abstract: It has recently been discovered that random quantum circuits provide an avenue to realize rich entanglement phase diagrams, which are hidden to standard expectation values of operators. Here we study (2+1)D random circuits with random Clifford unitary gates and measurements designed to stabilize trivial area law and topologically ordered phases. With competing single qubit Pauli-Z and toric code stabilizer measurements, in addition to random Clifford unitaries, we find a phase diagram involving a tricritical point that maps to (2+1)D percolation, a possibly stable critical phase, topologically ordered, trivial, and volume law phases, and lines of critical points separating them. With Pauli-Y single qubit measurements instead, we find an anisotropic self-dual tricritical point, with dynamical exponent $z \approx 1.46$, exhibiting logarithmic violation of the area law and an anomalous exponent for the topological entanglement entropy, which thus appears distinct from any known percolation fixed point. The phase diagram also hosts a measurement-induced volume law entangled phase in the absence of unitary dynamics.


title: Black hole spin evolution in warped accretion discs

abstract: Massive black holes (BHs) inhabiting galactic nuclei can be described by two parameters only, i.e. mass and spin, that change through cosmic time in response to accretion and merger events. While most numerical simulations accurately track the BH mass, spin evolution is rarely taken into account. In this work, we implement and validate a self-consistent sub-grid model for the evolution of the BH mass and spin via gas accretion in the hydrodynamics code GIZMO. The model assumes that accretion from resolved scales does not occur instantaneously, but is mediated by a sub-grid geometrically thin $\alpha$-disc. After validating our model semi-analytically, we test it in an idealized environment consisting of a circumnuclear disc, where gas accretion onto the accretion disc is consistently determined by GIZMO. In the absence of any accretion-related feedback, the spin evolution closely traces that observed in the semi-analytical models and depends on the free parameters of our implementation, such as the initial BH spin, angular momentum of the accretion disc, and the radius at which the gas inflow circularises. In GIZMO, we also couple our model with the biconical-outflow model presented in a companion paper, wherein the feedback axis is always aligned with the BH spin. In this last case, the evolution of the central BH differs significantly from the previous cases, since the feedback process modifies the gas dynamics and its inflow rates from resolved scales. Such an interaction cannot be modeled by simple semi-analytical models and should be treated using full N-body hydrodynamical simulations.


title: Weak gravity versus de Sitter

abstract: We show that one can uncover a Dine-Seiberg problem for de Sitter critical points in supergravity theories by utilizing the magnetic weak gravity conjecture. We present a large variety of N=2 gauged supergravity models that include vector multiplets and in all cases we find that the weak gravity conjecture threatens de Sitter. A common feature in all such examples is a degenerate mass matrix for the gravitini, which we therefore deem a swampland criterion for de Sitter critical points.


title: Spectral functions across an Insulator to Superconductor Transition

abstract: In a minimal 2-band model with attractive interactions between fermions, we calculate the gap to single and two-particle excitations, the band-dependent spectral functions, the superfluid density and compressibility using quantum Monte Carlo (QMC) methods. We find Fermi and Bose insulating phases with signatures of incipient pairing evident in the single-particle spectral functions, and a superconducting state with three different spectral functions: (i) both bands show "BCS" behavior in which the minimum gap locus occurs on a closed contour on the underlying Fermi surface; (ii) both bands show "BEC" behavior in which the minimum gap occurs at a point; and (iii) band selective spectral characteristics, in which one band shows "BCS" while the other shows "BEC" behavior. At large interactions, we find a Mott phase of rung bosons in which the filling is one boson for every two sites, half the typical density constraint for Mott insulators.


title: The Efficacy of Event Isotropy as an Event Shape Observable

abstract: Event isotropy $\mathcal{I}^\text{sph}$, an event shape observable that measures the distance of a final state from a spherically symmetric state, is designed for new physics signals that are far from QCD-like. Using a new technique for producing a wide variety of signals that can range from near-spherical to jetty, we compare event isotropy to other observables. We show that thrust $T$ and the $C$ parameter (and $\lambda_\text{max}$, the largest eigenvalue of the sphericity matrix) are strongly correlated and thus redundant, to a good approximation. By contrast, event isotropy adds considerable information, often serving to break degeneracies between signals that would have almost identical $T$ and $C$ distributions. Signals with broad distributions in $T$ (or $\lambda_\text{max}$) and in $\mathcal{I}^\text{sph}$ separately often have much narrower distributions, and are more easily distinguished, in the $({\mathcal{I}^\text{sph}},\lambda_\text{max})$ plane. An intuitive, semi-analytic estimation technique clarifies why this is the case and assists with the interpretation of the distributions.


title: Stable non-supersymmetric Anti-de Sitter vacua of massive IIA supergravity

abstract: Stable non-supersymmetric Anti-de Sitter (AdS) vacua of string theory are widely believed not to exist. In this letter, we analytically compute the full bosonic Kaluza-Klein spectrum around the ${\rm G_2}$-invariant non-supersymmetric $\textrm{AdS}_4$ solution of massive IIA supergravity, and show that it is perturbatively stable. We also provide evidence that six other non-supersymmetric $\textrm{AdS}_4$ solutions of massive IIA supergravity are perturbatively stable. Since previous studies have indicated that these AdS vacua may also be non-perturbatively stable, our findings pose a challenge to the Swampland Conjecture.


title: COOL-LAMPS I. An Extraordinarily Bright Lensed Galaxy at Redshift 5.04

abstract: We report the discovery of COOL J1241+2219, a strongly-lensed galaxy at redshift $z$=5.043$\pm$0.002 with observed magnitude $z_{AB}=20.47$, lensed by a moderate-mass galaxy cluster at $z$=1.001$\pm$0.001. COOL J1241+2219 is the brightest lensed galaxy currently known at optical and near-infrared wavelengths at $z$ $\gtrsim$ 5; it is $\sim$5 times brighter than the prior record-holder lensed galaxy, and several magnitudes brighter than the brightest unlensed galaxies known at these redshifts. It was discovered as part of COOL-LAMPS, a collaboration initiated to find strongly lensed systems in recent public optical imaging data. We characterise the lensed galaxy, as well as the central galaxy of the lensing cluster using ground-based $griz$JH imaging and optical spectroscopy. We report model-based magnitudes, and derive stellar masses, dust content, metallicity and star-formation rates via stellar-population synthesis modeling. Our lens mass modeling, based on ground-based imaging, implies a median source magnification of $\sim$30, which puts the stellar mass and star formation rate (in the youngest age bin, closest to the epoch of observation) at logM$_{*}$ = $10.11^{+0.21}_{-0.26}$ and SFR = $27^{+13}_{-9}$ M$_{\odot}$/yr, respectively. We constrain a star formation history for COOL J1241+2219 consistent with constant star formation across $\sim$1 Gyr of cosmic time, and that places this galaxy on the high-mass end of the star-forming main sequence. COOL J1241+2219 is 2-4 times more luminous than a galaxy with the characteristic UV luminosity at these redshifts. The UV continuum slope $\beta$= -2.2$\pm$0.2 places this galaxy on the blue side of the observed distribution of galaxies at $z$=5, although the lack of Ly$\alpha$ emission indicates dust sufficient to suppress this emission.


title: Rise and fall of post-starburst galaxies in Magneticum Pathfinder

abstract: Post-starburst galaxies (PSBs) belong to a short-lived transition population between star-forming (SF) and quiescent galaxies. Deciphering their heavily discussed evolutionary pathways is paramount to understanding galaxy evolution. We aim to determine the dominant mechanisms governing PSB evolution in both the field and in galaxy clusters. Using the cosmological hydrodynamical simulation suite Magneticum Pathfinder, we identify 647 PSBs with $z \sim 0$ stellar mass $M_* \geq 5 \cdot 10^{10} \, \mathrm{M_{\odot}}$. We track their galactic evolution, merger history, and black hole activity over a time-span of 3.6Gyr. Additionally, we study cluster PSBs identified at different redshifts and cluster masses. Independent of environment and redshift, we find that PSBs, like SF galaxies, have frequent mergers. At z=0, 89% of PSBs have experienced mergers and 65% had at least one major merger within the last 2.5Gyr, leading to strong star formation episodes. In fact, 23% of z=0 PSBs were rejuvenated during their starburst. Following the mergers, field PSBs are generally shutdown via a strong increase in AGN feedback (power output $P_{AGN,PSB} \geq 10^{56}\,$erg/Myr). We find agreement with observations for both stellar mass functions and z = 0.9 line-of-sight phase space distributions of PSBs in galaxy clusters. Finally, we find that $z \lesssim 0.5$ cluster PSBs are predominantly infalling, especially in high mass clusters and show no signs of enhanced AGN activity. Thus, we conclude that the majority of cluster PSBs are shutdown via an environmental quenching mechanism such as ram-pressure stripping, while field PSBs are mainly quenched by AGN feedback.


title: Opacity from Loops in AdS

abstract: We investigate how quantum dynamics affects the propagation of a scalar field in Lorentzian AdS. We work in momentum space, in which the propagator admits two spectral representations (denoted "conformal" and "momentum") in addition to a closed-form one, and all have a simple split structure. Focusing on scalar bubbles, we compute the imaginary part of the self-energy $ {\rm Im} \Pi$ in the three representations, which involves the evaluation of seemingly very different objects. We explicitly prove their equivalence in any dimension, and derive some elementary and asymptotic properties of $ {\rm Im} \Pi$. Using a WKB-like approach in the timelike region, we evaluate the propagator dressed with the imaginary part of the self-energy. We find that the dressing from loops exponentially dampens the propagator when one of the endpoints is in the IR region, rendering this region opaque to propagation. This suppression may have implications for field-theoretical model-building in AdS. We argue that in the effective theory (EFT) paradigm, opacity of the IR region induced by higher dimensional operators censors the region of EFT breakdown. This confirms earlier expectations from the literature. Specializing to AdS$_5$, we determine a universal contribution to opacity from gravity.


title: On the challenges of searching for GeV-scale long-lived particles at the LHC

abstract: Many models of dark matter predict long-lived particles (LLPs) that can give rise to striking signatures at the LHC. Existing searches for displaced vertices are however tailored towards heavy LLPs. In this work we show that this bias severely affects their sensitivity to LLPs with masses at the GeV scale. To illustrate this point we consider two dark sector models with light LLPs that decay hadronically: a strongly-interacting dark sector with long-lived exotic mesons, and a Higgsed dark sector with a long-lived dark Higgs boson. We study the sensitivity of an existing ATLAS search for displaced vertices and missing energy in these two models and find that current track and vertex cuts result in very low efficiency for light LLPs. To close this gap in the current search programme we suggest two possible modifications of the vertex reconstruction and the analysis cuts. We calculate projected exclusion limits for these modifications and show that they greatly enhance the sensitivity to LLPs with low mass or short decay lengths.


title: Residues, modularity, and the Cardy limit of the 4d $\mathcal{N}=4$ superconformal index

abstract: We compute the superconformal index of the $\mathcal{N}=4$ $SU(N)$ Yang-Mills theory through a residue calculation. The method is similar in spirit to the Bethe Ansatz formalism, except that all poles are explicitly known, and we do not require specialization of any of the chemical potentials. Our expression for the index allows us to revisit the Cardy limit using modular properties of four-dimensional supersymmetric partition functions. We find that all residues contribute at leading order in the Cardy limit. In a specific region of flavour chemical potential space, close to the two unrefined points, in fact all residues contribute universally. These universal residues precisely agree with the entropy functions of the asymptotically AdS$_5$ black hole and its "twin saddle" respectively. Finally, we discuss how our formula is suited to study the implications of four-dimensional modularity for the index beyond the Cardy limit.


title: Non-isotropic feedback from accreting spinning black holes

abstract: Active galactic nuclei (AGN) are massive black holes (BHs) caught in the act of accreting gas at the centre of their host galaxies. Part of the accreting mass is converted to energy and released into the surrounding medium, in a process loosely referred to as AGN feedback. Most numerical simulations include AGN feedback as a sub-grid model, wherein energy or momentum (or both) is coupled to the nearby gas. In this work, we implement a new momentum-driven model in the hydrodynamics code GIZMO, in which accretion from large scales is mediated by a sub-grid accretion disc model, and gas particles are stochastically kicked over a bi-conical region, to mimic observed kinetic winds. The feedback cone's axis can be set either parallel to the angular momentum of the gas surrounding the BH or to the BH spin direction, which is self-consistently evolved within the accretion-disc model. Using a circumnuclear disc (CND) as a test bed, we find that (i) the conical shape of the outflow is always visible and is weakly dependent on the launching orientation and aperture, resulting in comparable mass inflows and outflows; (ii) the cone's orientation is also similar amongst our tests, and it is not always the same as the initial value, due to the interaction with the CND playing a crucial role in shaping the outflow; and (iii) the velocity of the outflow, instead, differs and strongly depends on the interplay with the CND.


title: X-ray and Radio Bursts from the Magnetar 1E1547.0-5408

abstract: We report on simultaneous radio and X-ray observations of the radio-emitting magnetar 1E1547.0-5408 on 2009 January 25 and February 3, with the 64-m Parkes radio telescope and the Chandra and XMM-Newton X-ray observatories. The magnetar was observed in a period of intense X-ray bursting activity and enhanced X-ray emission. We report here on the detection of two radio bursts from 1E1547.0-5408, reminiscent of Fast Radio Bursts (FRBs). One of the radio bursts was anticipated by ~1s (about half a rotation period of the pulsar) by a bright SGR-like X-ray burst, resulting in a F_radio/F_X ~ 10^-9. Radio pulsations were not detected during the observation showing the FRB-like radio bursts, while they were detected in the previous radio observation. We also found that the two radio bursts are neither aligned with the latter radio pulsations nor with the peak of the X-ray pulse profile (phase shift of ~0.2). Comparing the luminosity of these FRB-like bursts and those reported from SGR1935+2154, we find that the wide range in radio efficiency and/or luminosity of magnetar bursts in the Galaxy may bridge the gap between "ordinary" pulsar radio bursts and the extragalactic FRB phenomenon.


title: Spectral appearance of the planetary-surface accretion shock: Global spectra and hydrogen-line profiles and fluxes

abstract: Hydrogen-line emission from an accretion shock has recently been observed at planetary-mass objects. Our previous work predicted the shock spectrum and luminosity for a shock on the circumplanetary disc. We extend this to the planet-surface shock. We calculate the global spectral energy distribution (SED) of accreting planets by combining our model emission spectra with photospheric SEDs, and predict the line-integrated flux for several hydrogen lines, especially H alpha, but also H beta, Pa alpha, Pa beta, Pa gamma, Br alpha, and Br gamma. We apply our non-equilibrium emission model to the surface accretion shock for a wide range of accretion rates Mdot and masses M_p . Fits to formation calculations provide radii and effective temperatures. Extinction by the surrounding material is neglected, which is arguably often relevant. We find that the line luminosity increases monotonically with Mdot and M_p , depending mostly on Mdot and weakly on M_p for the relevant range of parameters. The Lyman, Balmer, and Paschen continua can exceed the photosphere. The H beta line is fainter by 0 to 1 dex than H alpha, whereas other lines are weaker (by 1 to 3 dex). Shocks on the planet or the CPD surface are distinguishable at very high spectral resolution, but the planet surface shock likely dominates if both are present. Applied to recent non-detections of H alpha, our models imply looser constraints on the Mdot of putative planets than when extrapolating fits from the stellar regime. These hydrogen-line luminosity predictions are useful for interpreting (non-)detections of accreting planets.


title: The Inverse Seesaw Family: Dirac And Majorana

abstract: After developing a general criterion for deciding which neutrino mass models belong to the category of inverse seesaw models, we apply it to obtain the Dirac analogue of the canonical Majorana inverse seesaw model. We then generalize the inverse seesaw model and obtain a class of inverse seesaw mechanisms both for Majorana and Dirac neutrinos. We further show that many of the models have double or multiple suppressions coming from tiny symmetry breaking "$\mu$-terms". These models can be tested both in colliders and with the observation of lepton flavour violating processes.


title: The Swampland at Large Number of Space-Time Dimensions

abstract: We discuss some aspects of swampland constraints - especially the swampland distance conjecture - in a large number of space-time dimensions $D$. We analyze Kaluza-Klein (KK) states at large $D$ and find that some KK spectra possess an interesting dependence on $D$. On the basis of these observations we propose a new large dimension conjecture. We apply it to KK states of compactifications to anti-de Sitter backgrounds where it predicts an upper bound on the dimension of space-time as a function of its characteristic radius. We also apply our conjecture to black hole spacetimes, whose entropies have a $D$-dependence very similar to that of the KK spectrum.


title: Entanglement order parameters and critical behavior for topological phase transitions and beyond

abstract: Topological phases are exotic quantum phases which are lacking the characterization in terms of order parameters. In this paper, we develop a unified framework based on variational iPEPS for the quantitative study of both topological and conventional phase transitions through entanglement order parameters. To this end, we employ tensor networks with suitable physical and/or entanglement symmetries encoded, and allow for order parameters detecting the behavior of any of those symmetries, both physical and entanglement ones. First, this gives rise to entanglement-based order parameters for topological phases. These topological order parameters allow to quantitatively probe topological phase transitions and to identify their universal behavior. We apply our framework to the study of the Toric Code model in different magnetic fields, which in some cases maps to the (2+1)D Ising model. We identify 3D Ising critical exponents for the entire transition, consistent with those special cases and general belief. However, we moreover find an unknown critical exponent beta=0.021. We then apply our framework of entanglement order parameters to conventional phase transitions. We construct a novel type of disorder operator (or disorder parameter), which is non-zero in the disordered phase and measures the response of the wavefunction to a symmetry twist in the entanglement. We numerically evaluate this disorder operator for the (2+1)D transverse field Ising model, where we again recover a critical exponent hitherto unknown in the model, beta=0.024, consistent with the findings for the Toric Code. This shows that entanglement order parameters can provide additional means of characterizing the universal data both at topological and conventional phase transitions, and altogether demonstrates the power of this framework to identify the universal data underlying the transition.


title: Many-body nonlocality as a resource for quantum-enhanced metrology

abstract: We demonstrate that a many-body nonlocality is a resource for ultra-precise metrology. This result is achieved by linking the sensitivity of a quantum sensor with a combination of many-body correlation functions that witness the nonlocality. We illustrate our findings with some prominent examples---a collection of spins forming an Ising chain and a gas of ultra-cold atoms in any two-mode configuration.


title: Inverse Tides in Pulsating Binary Stars

abstract: In close binary stars, the tidal excitation of pulsations typically dissipates energy, causing the system to evolve towards a circular orbit with aligned and synchronized stellar spins. However, for stars with self-excited pulsations, we demonstrate that tidal interaction with unstable pulsation modes can transfer energy in the opposite direction, forcing the spins of the stars away from synchronicity, and potentially pumping the eccentricity and spin-orbit misalignment angle. This "inverse" tidal process only occurs when the tidally forced mode amplitude is comparable to the mode's saturation amplitude, and it is thus most likely to occur in main sequence gravity mode pulsators with orbital periods of a few days. We examine the long-term evolution of inverse tidal action, finding the stellar rotation rate can potentially be driven to a very large or very small value, while maintaining a large spin-orbit misalignment angle. Several recent asteroseismic analyses of pulsating stars in close binaries have revealed extremely slow core rotation periods, which we attribute to the action of inverse tides.


title: Nonlinear evolution problems with singular coefficients in the lower order terms

abstract: We consider a Cauchy Dirichlet problem for a quasilinear second order parabolic equation with lower order term driven by a singular coefficient. We establish an existence result to such a problem and we describe the time behavior of the solution in the case of the infinite time horizon.


title: The VANDELS survey: the relation between UV continuum slope and stellar metallicity in star-forming galaxies at z~3

abstract: The estimate of stellar metallicities (Z*) of high-z galaxies are of paramount importance in order to understand the complexity of dust effects and the reciprocal interrelations among stellar mass, dust attenuation, stellar age, and metallicity. Benefiting from uniquely deep FUV spectra of >500 star-forming galaxies at redshifts 2<z<5 extracted from the VANDELS survey and stacked in bins of stellar mass (M*) and UV continuum slope (beta), we estimate their stellar metallicities Z* from stellar photospheric absorption features at 1501 and 1719 Angstrom, which are calibrated with Starburst99 models and are largely unaffected by stellar age, dust, IMF, nebular continuum or interstellar absorption. Comparing them to photometric based spectral slopes in the range 1250-1750 Angstrom, we find that the stellar metallicity increases by ~0.5 dex from beta ~ -2 to beta ~ -1 (1 < A(1600) < 3.2), and a dependence with beta holds at fixed UV absolute luminosity M(UV) and stellar mass up to 10^(9.65) Msun. As a result, the metallicity is a fundamental ingredient for properly rescaling dust corrections based on M(UV) and M*. Using the same absorption features, we analyze the mass-metallicity relation (MZR) and find it is consistent with the previous VANDELS estimation based on a global fit of the FUV spectra. Similarly, we do not find a significant evolution between z=2 and z=3.5. Finally, the slopes of our MZR and Z*-beta relation are in agreement with the predictions of well-studied semi-analytic models of galaxy formation (SAM), while some tensions with observations remain as to the absolute metallicity normalization. The relation between UV slope and stellar metallicity is fundamental for the exploitation of large volume surveys with next-generation telescopes and for the physical characterization of galaxies in the first billion years of our Universe.


title: The nature of giant clumps in high-z discs: a deep-learning comparison of simulations and observations

abstract: We use deep learning to explore the nature of observed giant clumps in high-redshift disc galaxies, based on their identification and classification in cosmological simulations. Simulated clumps are detected using the 3D gas and stellar densities in the VELA zoom-in cosmological simulation suite, with $\sim \!\! 25\ \!\rm{pc}$ maximum resolution, targeting main sequence galaxies at $1\!<\!z\!<\!3$. The clumps are classified as long-lived clumps (LLCs) or short-lived clumps (SLCs) based on their longevity in the simulations. We then train neural networks to detect and classify the simulated clumps in mock, multi-color, dusty and noisy HST-like images. The clumps are detected using an encoder-decoder convolutional neural network (CNN), and are classified according to their longevity using a vanilla CNN. Tests using the simulations show our detector and classifier to be $\sim80\%$ complete and $\sim80\%$ pure for clumps more massive than $\sim10^{7.5}\rm{M_\odot}$. When applied to observed galaxies in the CANDELS/GOODS S+N fields, we find both types of clumps to appear in similar abundances in the simulations and the observations. LLCs are, on average, more massive than SLCs by $\sim 0.5\ \rm{dex}$, and they dominate the clump population above $M_{\rm c}\gtrsim 10^{7.6}\ \rm{M_\odot}$. LLCs tend to be found closer to the galactic centre, indicating clump migration to the centre or preferential formation at smaller radii. The LLCs are found to reside in high mass galaxies, indicating better clump survivability under supernova feedback there, due to clumps being more massive in these galaxies. We find the clump properties in the simulations and the observations to agree within a factor of $\sim\! 2$.


title: Optimizing multi-rendezvous spacecraft trajectories: $\Delta V$ matrices and sequence selection

abstract: Multi-rendezvous spacecraft trajectory optimization problems are notoriously difficult to solve. For this reason, the design space is usually pruned by using heuristics and past experience. As an alternative, the current research explores some properties of $\Delta V$ matrices which provide the minimum $\Delta V$ values for a transfer between two celestial bodies for various times of departure and transfer duration values. These can assist in solving multi-rendezvous problems in an automated way. The paper focuses on the problem of, given a set of candidate objects, how to find the sequence of $N$ objects to rendezvous with that minimizes the total $\Delta V$ required. Transfers are considered as single algebraic objects corresponding to $\Delta V$ matrices, which allow intuitive concatenation via a generalized summation. Waiting times, both due to mission requirements and prospects for cheaper and faster future transfers, are also incorporated in the $\Delta V$ matrices. A transcription of the problem as a shortest path search on a graph can utilize a range of available efficient shortest path solvers. Given an efficient $\Delta V$ matrix estimator, the new paradigm proposed here is believed to offer an alternative to the pruning techniques currently used.


title: Simulation of the drawdown and its duration in L\'{e}vy models via stick-breaking Gaussian approximation

abstract: We develop a computational method for expected functionals of the drawdown and its duration in exponential L\'evy models. It is based on a novel simulation algorithm for the joint law of the state, supremum and time the supremum is attained of the Gaussian approximation of a general L\'evy process. We bound the bias for various locally Lipschitz and discontinuous payoffs arising in applications and analyse the computational complexities of the corresponding Monte Carlo and multilevel Monte Carlo estimators. Monte Carlo methods for L\'evy processes (using Gaussian approximation) have been analysed for Lipschitz payoffs, in which case the computational complexity of our algorithm is up to two orders of magnitude smaller when the jump activity is high. At the core of our approach are bounds on certain Wasserstein distances, obtained via the novel SBG coupling between a L\'evy process and its Gaussian approximation. Numerical performance, based on the implementation in the dedicated GitHub repository, exhibits a good agreement with our theoretical bounds.


title: Learning Latent Representations to Influence Multi-Agent Interaction

abstract: Seamlessly interacting with humans or robots is hard because these agents are non-stationary. They update their policy in response to the ego agent's behavior, and the ego agent must anticipate these changes to co-adapt. Inspired by humans, we recognize that robots do not need to explicitly model every low-level action another agent will make; instead, we can capture the latent strategy of other agents through high-level representations. We propose a reinforcement learning-based framework for learning latent representations of an agent's policy, where the ego agent identifies the relationship between its behavior and the other agent's future strategy. The ego agent then leverages these latent dynamics to influence the other agent, purposely guiding them towards policies suitable for co-adaptation. Across several simulated domains and a real-world air hockey game, our approach outperforms the alternatives and learns to influence the other agent.


title: Random Walk and Trapping of Interplanetary Magnetic Field Lines: Global Simulation, Magnetic Connectivity, and Implications for Solar Energetic Particles

abstract: The random walk of magnetic field lines is an important ingredient in understanding how the connectivity of the magnetic field affects the spatial transport and diffusion of charged particles. As solar energetic particles (SEPs) propagate away from near-solar sources, they interact with the fluctuating magnetic field, which modifies their distributions. We develop a formalism in which the differential equation describing the field line random walk contains both effects due to localized magnetic displacements and a non-stochastic contribution from the large-scale expansion. We use this formalism together with a global magnetohydrodynamic simulation of the inner-heliospheric solar wind, which includes a turbulence transport model, to estimate the diffusive spreading of magnetic field lines that originate in different regions of the solar atmosphere. We first use this model to quantify field line spreading at 1 au, starting from a localized solar source region, and find rms angular spreads of about 20{\deg} - 60{\deg}. In the second instance, we use the model to estimate the size of the source regions from which field lines observed at 1 au may have originated, thus quantifying the uncertainty in calculations of magnetic connectivity; the angular uncertainty is estimated to be about 20{\deg}. Finally, we estimate the filamentation distance, i.e., the heliocentric distance up to which field lines originating in magnetic islands can remain strongly trapped in filamentary structures. We emphasize the key role of slab-like fluctuations in the transition from filamentary to more diffusive transport at greater heliocentric distances.


title: Gaia-based Isochronal, Kinematic, and Spatial Analysis of the $\epsilon$ Cha Association

abstract: The precise parallax, proper motion, and photometric measurements contained in Gaia Data Release 2 (DR2) offer the opportunity to re-examine the membership and ages of Nearby Young Moving Groups (NYMGs), i.e., loose groups of stars of age $\lesssim100$ Myr in the solar vicinity. Here, we analyze the available DR2 data for members and candidate members of the $\epsilon$ Cha Association ($\epsilon$CA) which, at an estimated age of $\sim$3--5 Myr, has previously been identified as among the youngest NYMGs. The several dozen confirmed members of $\epsilon$CA include MP Mus and T Cha, two of the nearest stars of roughly solar mass that are known to host primordial protoplanetary disks, and the Herbig Ae/Be star HD 104237A. We have used Gaia DR2 data to ascertain the Galactic positions and kinematics and color-magnitude diagram positions of $\epsilon$CA members and candidates, so as to reassess their membership status and thereby refine estimates of the distance, age, multiplicity, and disk fraction of the group. Our analysis yields 36 \textit{bona fide} $\epsilon$CA members, as well as 20 provisional members, including 3 new members identified here as comoving companions to previously known $\epsilon$CA stars. We determine a mean distance to $\epsilon$CA of 101.0$\pm$4.6 pc and confirm that, at an age of $5^{+3}_{-2}$ Myr, $\epsilon$CA represents the youngest stellar group within $\sim$100 pc of Earth. We identify several new photometric binary candidates, bringing the overall multiplicity fraction (MF) of $\epsilon$CA to 40\%, intermediate between the MFs of young T associations and the field.


title: The Influence of the Buffer Size in Packet Loss for Competing Multimedia and Bursty Traffic

abstract: This work presents an analysis of the effect of the access router buffer size on packet loss rate and how it can affect the QoS of multimedia services when bursty traffic is present. VoIP traffic, real traces of viedoconferencing and videosurvellance are used in two different scenarios with medium link utilization. The study shows that the bursty nature of some applications may impair the MOS of voice calls especially when a certain number of bursts overlap. When link utilization is above 70% good values of VoIP QoS cannot be obtained.


title: doc2dial: A Goal-Oriented Document-Grounded Dialogue Dataset

abstract: We introduce doc2dial, a new dataset of goal-oriented dialogues that are grounded in the associated documents. Inspired by how the authors compose documents for guiding end users, we first construct dialogue flows based on the content elements that corresponds to higher-level relations across text sections as well as lower-level relations between discourse units within a section. Then we present these dialogue flows to crowd contributors to create conversational utterances. The dataset includes about 4800 annotated conversations with an average of 14 turns that are grounded in over 480 documents from four domains. Compared to the prior document-grounded dialogue datasets, this dataset covers a variety of dialogue scenes in information-seeking conversations. For evaluating the versatility of the dataset, we introduce multiple dialogue modeling tasks and present baseline approaches.


title: Impact of electron temperature anisotropy on the collisionless tearing mode instability in the presence of a strong guide field

abstract: We derive and analyze a dispersion relation for the growth rate of collisionless tearing modes, driven by electron inertia and accounting for equilibrium electron temperature anisotropy in a strong guide field regime. For this purpose, a new gyrofluid model is derived and subsequently simplified to make the derivation of the dispersion relation treatable analytically. The main simplifying assumptions consist in assuming cold ions, neglecting electron finite Larmor radius effects, decoupling ion gyrocenter fluctuations and considering $\beta_{\perp_e} \ll 1$, with $\beta_{\perp_e}$ indicating the ratio between the perpendicular electron thermal pressure and the magnetic pressure exerted by the guide field. This simplified version of the gyrofluid model is shown to possess a noncanonical Hamiltonian structure. The dispersion relation is obtained by applying the theory of asymptotic matching and does not predict an enhancement of the growth rate as the ratio $\Theta_e$, between perpendicular and parallel equilibrium electron temperatures, increases. This indicates a significant difference with respect to the case of absent or moderate guide field. For an equilibrium magnetic shear length of the order of the perpendicular sonic Larmor radius and at a fixed $\beta_{\perp_e}$, we obtain that the tearing mode in the strong guide field regime gets actually weakly damped, as $\Theta_e$ increases. In the isotropic limit $\Theta_e=1$, the dispersion relation reduces to a previously known formula. The analytical predictions are tested against numerical simulations showing a very good quantitative agreement. We also provide a detailed discussion of the range of validity of the derived dispersion relation and of the compatibility among the different adopted assumptions.


title: The critical number of $I_{1,t}$-free triangle-free binary matroids

abstract: A simple binary matroid, viewed as a restriction of a finite binary projective geometry $PG(n-1,2)$, is $I_{1,t}$-free if for any rank-$t$ flat of $PG(n-1,2)$, its intersection with the matroid is not a one-element set. In this paper, we show that the simple $I_{1,t}$-free and triangle-free binary matroids have bounded critical number for any $t \geq 1$.


title: Impact of the Central Frequency of Environment on Non-Markovian Dynamics in Piezoelectric Optomechanical Devices

abstract: The piezoelectric optomechanical devices supply a promising experimental platform to realize the coherent and effective control and measurement of optical circuits working in Terahertz (THz) frequencies via superconducting electron devices typically working in Radio (MHz) frequencies. However, quantum fluctuations are unavoidable when the size of mechanical oscillators enters into the nanoscale. The consequences of the noisy environment are still challenging due to the lack of analytical tools. In this paper, a semi-classical and full-quantum model of piezoelectric optomechanical systems coupled to a noisy bosonic quantum environment are introduced and solved in terms of quantum-state diffusion (QSD) trajectories in the non-Markovian regime. We show that the noisy environment, particularly the central frequency of the environment, can enhance the entanglement generation between optical cavities and LC circuits in some parameter regimes. Moreover, we observe the critical points in the coefficient functions, which can lead the different behaviors in the system. Besides, we also witness the entanglement transfers between macroscopic objects due to the memory effect of the environment. Our work can be applied in the fields of electric/ optical switches, and long-distance distribution in a large-scale quantum network.


title: An extension of the Siegel-Walfisz theorem

abstract: We extend the Siegel-Walfisz theorem to a family of integer sequences that are characterized by constraints on the size of the prime factors. Besides prime powers, this family includes smooth numbers, almost primes and practical numbers.


title: Alloying behavior of wide band gap alkaline-earth chalcogenides

abstract: Alloying is a powerful tool for tuning materials that facilitates the targeted design of desirable properties for a variety of applications. In this work, we provide a comprehensive investigation of the synthetic accessibility and electronic properties of nine alkaline-earth chalcogenide anion alloys (CaS$_{1-x}$O$_x$, CaS$_{1-x}$Se$_x$, CaS$_{1-x}$Te$_x$, SrS$_{1-x}$O$_x$, SrS$_{1-x}$Se$_x$, SrS$_{1-x}$Te$_x$, MgS$_{1- x}$O$_x$, MgS$_{1-x}$Se$_x$, and MgS$_{1-x}$Te$_x$). We show that isostructural alloying within the rock salt structure is favored for all systems except MgS$_{1-x}$Te$_x$, which is predicted to be a heterostructural alloy between the rock salt and wurtzite structures. Alloys of S and Se are shown to be readily accessible for all cations with low miscibility critical temperatures, enabling continuous tuning of electronic properties across this composition space. Alloys of S and Te have higher critical temperatures but may be accessible through non-equilibrium synthesis strategies and are predicted here to have desirable electronic properties for optoelectronics with wide band gaps and lower effective masses than alloys of S and Se. Anion alloying in MgS$_{1-x}$Te$_x$ stabilizes the wurtzite structure across a significant fraction of composition space, which may make it of particular interest as a transparent conducting material due to its lower effective masses and a higher band gap than the rock salt structure. Zero-point corrected random phase approximation (RPA) energies were computed to resolve the small polymorph energy differences of the Mg compounds and are shown to be critical for accurately describing the thermodynamic properties of the corresponding alloys.


title: Bayesian nonparametric modelling of sequential discoveries

abstract: We aim at modelling the appearance of distinct tags in a sequence of labelled objects. Common examples of this type of data include words in a corpus or distinct species in a sample. These sequential discoveries are often summarised via accumulation curves, which count the number of distinct entities observed in an increasingly large set of objects. We propose a novel Bayesian nonparametric method for species sampling modelling by directly specifying the probability of a new discovery, therefore allowing for flexible specifications. The asymptotic behavior and finite sample properties of such an approach are extensively studied. Interestingly, our enlarged class of sequential processes includes highly tractable special cases. We present a subclass of models characterized by appealing theoretical and computational properties. Moreover, due to strong connections with logistic regression models, the latter subclass can naturally account for covariates. We finally test our proposal on both synthetic and real data, with special emphasis on a large fungal biodiversity study in Finland.


title: Successful Common Envelope Ejection and Binary Neutron Star Formation in 3D Hydrodynamics

abstract: A binary neutron star merger has been observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding of the multi-spatial multi-timescale evolution for the systems that lead to compact binaries. In this paper, we report on the first successful simulations of common envelope ejection leading to binary neutron star formation in 3D hydrodynamics. We simulate the dynamical inspiral phase of the interaction between a 12$M_\odot$ red supergiant and a 1.4$M_\odot$ neutron star for different initial separations and initial conditions. For all of our simulations, we find complete envelope ejection and final orbital separations of $a_{\rm f} \approx 1.3$-$5.1 R_\odot$ depending on the simulation and criterion, leading to binary neutron stars that can merge within a Hubble time. We find $\alpha_{\rm CE}$-equivalent efficiencies of $\approx 0.1$-$2.7$ depending on the simulation and criterion, but this may be specific for these extended progenitors. We fully resolve the core of the star to $\lesssim 0.005 R_\odot$ and our 3D hydrodynamics simulations are informed by an adjusted 1D analytic energy formalism and a 2D kinematics study in order to overcome the prohibitive computational cost of simulating these systems. The framework we develop in this paper can be used to simulate a wide variety of interactions between stars, from stellar mergers to common envelope episodes leading to GW sources.


title: Steady State Analysis of Episodic Reinforcement Learning

abstract: This paper proves that the episodic learning environment of every finite-horizon decision task has a unique steady state under any behavior policy, and that the marginal distribution of the agent's input indeed converges to the steady-state distribution in essentially all episodic learning processes. This observation supports an interestingly reversed mindset against conventional wisdom: While the existence of unique steady states was often presumed in continual learning but considered less relevant in episodic learning, it turns out their existence is guaranteed for the latter. Based on this insight, the paper unifies episodic and continual RL around several important concepts that have been separately treated in these two RL formalisms. Practically, the existence of unique and approachable steady state enables a general way to collect data in episodic RL tasks, which the paper applies to policy gradient algorithms as a demonstration, based on a new steady-state policy gradient theorem. Finally, the paper also proposes and experimentally validates a perturbation method that facilitates rapid steady-state convergence in real-world RL tasks.