title: Designing Inherently Interpretable Machine Learning Models

abstract: Interpretable machine learning (IML) becomes increasingly important in highly regulated industry sectors related to the health and safety or fundamental rights of human beings. In general, the inherently IML models should be adopted because of their transparency and explainability, while black-box models with model-agnostic explainability can be more difficult to defend under regulatory scrutiny. For assessing inherent interpretability of a machine learning model, we propose a qualitative template based on feature effects and model architecture constraints. It provides the design principles for high-performance IML model development, with examples given by reviewing our recent works on ExNN, GAMI-Net, SIMTree, and the Aletheia toolkit for local linear interpretability of deep ReLU networks. We further demonstrate how to design an interpretable ReLU DNN model with evaluation of conceptual soundness for a real case study of predicting credit default in home lending. We hope that this work will provide a practical guide of developing inherently IML models in high risk applications in banking industry, as well as other sectors.


title: UnProjection: Leveraging Inverse-Projections for Visual Analytics of High-Dimensional Data

abstract: Projection techniques are often used to visualize high-dimensional data, allowing users to better understand the overall structure of multi-dimensional spaces on a 2D screen. Although many such methods exist, comparably little work has been done on generalizable methods of inverse-projection -- the process of mapping the projected points, or more generally, the projection space back to the original high-dimensional space. In this paper we present NNInv, a deep learning technique with the ability to approximate the inverse of any projection or mapping. NNInv learns to reconstruct high-dimensional data from any arbitrary point on a 2D projection space, giving users the ability to interact with the learned high-dimensional representation in a visual analytics system. We provide an analysis of the parameter space of NNInv, and offer guidance in selecting these parameters. We extend validation of the effectiveness of NNInv through a series of quantitative and qualitative analyses. We then demonstrate the method's utility by applying it to three visualization tasks: interactive instance interpolation, classifier agreement, and gradient visualization.


title: Eigen's Paradox and the Quasispecies Model in a Non-Archimedean Framework

abstract: In this article we present a new p-adic generalization of the Eigen-Schuster model where the genomes (sequences) are represented by words written in the alphabet 0,1,,p-1, where p is a prime number, with a time variable length. The time evolution of the concentration of a sequence is controlled by a p-adic evolution equation. The long term behavior of the concentration of a sequence depends on a fitness function f, a mutation measure Q, and an initial concentration distribution. The new model provides essentially two types of asymptotic scenarios for evolution. If the complexity of sequences grows at the right pace, then in the long term the survival is assured. This agrees with the fact that larger genome size improves the replication fidelity. In other case, the sequences cannot copy themselves with sufficiently fidelity, and in the long term they will not survive. Eigen's paradox is one, among the infinitely many, possible scenarios of the evolution in the long term. The mathematical formulation of this fact requires solving the Cauchy problem for the p-adic Eigen-Schuster model in a rigorous mathematical way. This requires imposing restrictions on the fitness function and on the mutation measure, among other conditions. The study of the mentioned initial value problem requires techniques of p-adic wavelets and heat kernels developed in the last 35 years.


title: Superdense beaming of axion dark matter in the vicinity of the light cylinder of pulsars

abstract: In this article we treat the non-adiabatic photon-to-axion resonant conversion of curvature radiation, synchrotron emission and inverse Compton scattering dominating the spectral density function of pulsars. First, we introduce emission models and benchmark observational data.vWe adopt a state-of-the-art density profile that relieves tension with the quantum electrodynamics vacuum polarization effect in highly magnetic stars, leading to efficient mixing. Then, we estimate the dark matter flux induced by photon-axion oscillation across the light cylinder of the neutron star. We find that pulsars might produce axion overdensities many orders of magnitude over the occupation number of dark matter in the Galactic halo within a broad parameter space. We point out possible new methods for axion detection derived from these results and other future lines of work.

categories: astro-ph.HE astro-ph.CO astro-ph.SR hep-ph hep-th


title: Demonstration of entanglement and coherence in GHZ-like state when exposed to classical environments with power-law noise

abstract: Entanglement and coherence protection are investigated using the dynamical map of three non-interacting qubits that are initially prepared as maximally entangled GHZ-like states coupled to external fields in solid-state and superconducting materials. Thermal fluctuations and resistance in these materials produce power-law (PL) noise, which we assume controls external fields in three different configurations with single or multiple noise sources. The genuine response of isolated environments to entanglement and coherence retention is analyzed. We also briefly discuss the initial purity and relative efficiency of the GHZ-like states. Unlike the multiple PL noise sources, exposure, the GHZ-class state remains partially entangled and coherent for an indefinite time when subject to single noise source. However, long-term non-local correlation and coherence are still feasible under multiple noise sources. Due to the lack of back-action of the environments, the conversion of the free state into the resource GHZ-class states is not allowed. The parameter optimization, in addition to the noise phase, regulates disorders, noise effects, and memory properties in the current environments. Unlike the bipartite states and tripartite W-type states, the GHZ-like state has been shown to have positive traits of preserving entanglement and coherence against PL noise.


title: Sunyaev Zel'dovich effect in galaxy clusters cavities: thermal or non-thermal origin?

abstract: Several galaxy clusters host X-ray cavities, often filled with relativistic electrons emitting in the radio band. In the cluster MS 0735.6+7421 the cavities have been detected through the Sunyaev Zel'dovich (SZ) effect, but it has not been possible to determine if this effect is thermal (produced by a very high temperature gas filling the cavity) or non-thermal (produced by the relativistic electrons that produce the diffuse radio emission detected in the cavity). In this paper we discuss the role of the density of the high temperature gas inside the cavities in determining whether the dominant SZ effect is the thermal or the non-thermal one, and how it can be possible to distinguish between the two possibilities, discussing the role of observations at higher energy bands.


title: TOI-2257 b: A highly eccentric long-period sub-Neptune transiting a nearby M dwarf

abstract: Thanks to the relative ease of finding and characterizing small planets around M dwarf stars, these objects have become cornerstones in the field of exoplanet studies. The current paucity of planets in long-period orbits around M dwarfs make such objects particularly compelling as they provide clues about the formation and evolution of these systems. In this study, we present the discovery of TOI-2257 b (TIC 198485881), a long-period (35 d) sub-Neptune orbiting an M3 star at 57.8pc. Its transit depth is about 0.4%, large enough to be detected with medium-size, ground-based telescopes. The long transit duration suggests the planet is in a highly eccentric orbit ($e \sim 0.5$), which would make it the most eccentric planet that is known to be transiting an M-dwarf star. We combined TESS and ground-based data obtained with the 1.0-m SAINT-EX, 0.60-m TRAPPIST-North and 1.2-m FLWO telescopes to find a planetary size of 2.2 $R_{\oplus}$ and an orbital period of 35.19 days. In addition, we make use of archival data, high-resolution imaging, and vetting packages to support our planetary interpretation. With its long period and high eccentricity, TOI-2257 b falls in a novel slice of parameter space. Despite the planet's low equilibrium temperature ($\sim$ 256 K), its host star's small size ($R_* = 0.311 \pm{0.015}$) and relative infrared brightness (K$_{mag}$ = 10.7) make it a suitable candidate for atmospheric exploration via transmission spectroscopy.


title: Spiking Generative Adversarial Networks With a Neural Network Discriminator: Local Training, Bayesian Models, and Continual Meta-Learning

abstract: Neuromorphic data carries information in spatio-temporal patterns encoded by spikes. Accordingly, a central problem in neuromorphic computing is training spiking neural networks (SNNs) to reproduce spatio-temporal spiking patterns in response to given spiking stimuli. Most existing approaches model the input-output behavior of an SNN in a deterministic fashion by assigning each input to a specific desired output spiking sequence. In contrast, in order to fully leverage the time-encoding capacity of spikes, this work proposes to train SNNs so as to match distributions of spiking signals rather than individual spiking signals. To this end, the paper introduces a novel hybrid architecture comprising a conditional generator, implemented via an SNN, and a discriminator, implemented by a conventional artificial neural network (ANN). The role of the ANN is to provide feedback during training to the SNN within an adversarial iterative learning strategy that follows the principle of generative adversarial network (GANs). In order to better capture multi-modal spatio-temporal distribution, the proposed approach -- termed SpikeGAN -- is further extended to support Bayesian learning of the generator's weight. Finally, settings with time-varying statistics are addressed by proposing an online meta-learning variant of SpikeGAN. Experiments bring insights into the merits of the proposed approach as compared to existing solutions based on (static) belief networks and maximum likelihood (or empirical risk minimization).


title: Rovibrational-Specific QCT and Master Equation Study on $\text{N}_2(\text{X}^1\Sigma_g^+)$+$\text{O}({}^3\text{P})$ and $\text{NO}(\text{X}^2\Pi)$+$\text{N}({}^4\text{S})$ Systems in High-Energy Collisions

abstract: This work presents a detailed investigation of the energy transfer and dissociation mechanisms in $\text{N}_2(\text{X}^1\Sigma_g^+)$+$\text{O}({}^3\text{P})$ and $\text{NO}(\text{X}^2\Pi)$+$\text{N}({}^4\text{S})$ systems using rovibrational-specific quasi-classical trajectory (QCT) and master equation analyses. The complete set of state-to-state kinetic data, obtained via QCT, allows for an in-depth investigation of the Zel'dovich mechanism leading to the formation of $\text{NO}$ molecules at microscopic and macroscopic scales. The master equation analysis demonstrates that the low-lying vibrational states of $\text{N}_2$ and $\text{NO}$ have dominant contributions to the $\text{NO}$ formation and the corresponding extinction of $\text{N}_2$ through the exchange process. For the considered temperature range, it is found that while nearly 50% of the dissociation processes for $\text{N}_2$ and $\text{NO}$ occurs in the molecular quasi-steady-state (QSS) regime, the amount of the Zel'dovich reaction is zero. Using the QSS approximation to model the Zel'dovich mechanism leads to an overestimation of $\text{NO}$ production by more than a factor of 4 in the high-temperature range. The breakdown of this well-known approximation has profound consequences for the approaches that heavily rely on the validity of QSS assumption in hypersonic applications. The investigation of the rovibrational state population dynamics reveals substantial similarity among different chemical systems for the energy transfer and the dissociation processes, providing promising physical foundations for the use of reduced-order strategies to other chemical systems without significant loss of accuracy.


title: Black Hole Images as Tests of General Relativity: Effects of Spacetime Geometry

abstract: The images of supermassive black holes surrounded by optically-thin, radiatively-inefficient accretion flows, like those observed with the Event Horizon Telescope, are characterized by a bright ring of emission surrounding the black-hole shadow. In the Kerr spacetime this bright ring, when narrow, closely traces the boundary of the shadow and can, with appropriate calibration, serve as its proxy. The present paper expands the validity of this statement by considering two particular spacetime geometries: a solution to the field equations of a modified gravity theory and another that parametrically deviates from Kerr but recovers the Kerr spacetime when its deviation parameters vanish. A covariant, axisymmetric analytic model of the accretion flow based on conservation laws and spanning a broad range of plasma conditions is utilized to calculate synthetic non-Kerr black-hole images, which are then analysed and characterized. We find that in all spacetimes: (i) it is the gravitationally-lensed unstable photon orbit that plays the critical role in establishing the diameter of the rings observed in black-hole images, not the event horizon or the innermost stable circular orbit, (ii) bright rings in these images scale in size with, and encompass, the boundaries of the black-hole shadows, even when deviating significantly from Kerr, and (iii) uncertainties in the physical properties of the accreting plasma introduce subdominant corrections to the relation between the diameter of the image and the diameter of the black-hole shadow. These results provide important new theoretical justification for using black-hole images to probe and test the spacetimes of supermassive black holes.


title: The Influence of 10 Unique Chemical Elements in Shaping the Distribution of Kepler Planets

abstract: The chemical abundances of planet-hosting stars offer a glimpse into the composition of planet-forming environments. To further understand this connection, we make the first ever measurement of the correlation between planet occurrence and chemical abundances for ten different elements (C, Mg, Al, Si, S, K, Ca, Mn, Fe, and Ni). Leveraging data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) and Gaia to derive precise stellar parameters ($\sigma_{R_\star}\approx2.3\%$, $\sigma_{M_\star}\approx4.5\%$) for a sample of 1,018 Kepler Objects of Interest, we construct a sample of well-vetted Kepler planets with precisely measured radii ($\sigma_{R_p}\approx3.4\%$). After controlling for biases in the Kepler detection pipeline and the selection function of the APOGEE survey, we characterize the relationship between planet occurrence and chemical abundance as the number density of nuclei of each element in a star's photosphere raised to a power, $\beta$. $\beta$ varies by planet type, but is consistent within our uncertainties across all ten elements. For hot planets ($P$ = 1-10 days), an enhancement in any element of 0.1 dex corresponds to an increased occurrence of $\approx$20% for Super-Earths ($R_p=1-1.9R_\oplus$) and $\approx$60% for Sub-Neptunes ($R_p=1.9-4R_\oplus$). Trends are weaker for warm ($P$ = 10-100 days) planets of all sizes and for all elements, with the potential exception of Sub-Saturns ($R_p=4-8R_\oplus$). Finally, we conclude this work with a caution to interpreting trends between planet occurrence and stellar age due to degeneracies caused by Galactic chemical evolution and make predictions for planet occurrence rates in nearby open clusters to facilitate demographics studies of young planetary systems.


title: Meta-Learning to Improve Pre-Training

abstract: Pre-training (PT) followed by fine-tuning (FT) is an effective method for training neural networks, and has led to significant performance improvements in many domains. PT can incorporate various design choices such as task and data reweighting strategies, augmentation policies, and noise models, all of which can significantly impact the quality of representations learned. The hyperparameters introduced by these strategies therefore must be tuned appropriately. However, setting the values of these hyperparameters is challenging. Most existing methods either struggle to scale to high dimensions, are too slow and memory-intensive, or cannot be directly applied to the two-stage PT and FT learning process. In this work, we propose an efficient, gradient-based algorithm to meta-learn PT hyperparameters. We formalize the PT hyperparameter optimization problem and propose a novel method to obtain PT hyperparameter gradients by combining implicit differentiation and backpropagation through unrolled optimization. We demonstrate that our method improves predictive performance on two real-world domains. First, we optimize high-dimensional task weighting hyperparameters for multitask pre-training on protein-protein interaction graphs and improve AUROC by up to 3.9%. Second, we optimize a data augmentation neural network for self-supervised PT with SimCLR on electrocardiography data and improve AUROC by up to 1.9%.


title: Wannier Diagram and Brown-Zak Fermions of Graphene on Hexagonal Boron-Nitride

abstract: The moir\'e potential of graphene on hexagonal boron nitride (hBN) generates a supercell sufficiently large as to thread a full magnetic flux quantum $\Phi_0$ for experimentally accessible magnetic field strengths. Close to rational fractions of $\Phi_0$, $p/q \cdot\Phi_0$, magnetotranslation invariance is restored giving rise to Brown-Zak fermions featuring the same dispersion relation as in the absence of the field. Employing a highly efficient numerical approach we have performed the first realistic simulation of the magnetoconductance for a 250 nm wide graphene ribbon on hexagonal boron nitride using a full ab-initio derived parametrization including strain. The resulting Hofstadter butterfly is analyzed in terms of a novel Wannier diagram for Landau spectra of Dirac particles that includes the lifting of the spin and valley degeneracy by the magnetic field and the moir\'e potential. This complex diagram can account for many experimentally observed features on a single-particle level, such as spin and valley degeneracy lifting and a non-periodicidy in $\Phi_0$.


title: CW Operation of the European XFEL: SC-Gun Injector Optimization, S2E Calculations and SASE Performance

abstract: This note presents optimization results for an injector with a super-conducting gun for the cw operation of the European XFEL, the corresponding start-to-end simulation calculations for the beam transport to the undulators and SASE calculations of the X-ray intensities for the achievable photon energies. The optimizations for an injector with a super-conducting gun are done mainly for an already experimentally shown accelerating gradient of 40 MV/m. Different working points with respect to transverse emittance and bunch length are chosen to see whether the lower transverse emittance of the longer bunches survives the necessary stronger bunch compression to achieve comparable peak currents in the undulators. Also, results for possible further improvement of the gun gradient to 50 MV/m and a modified transverse laser profile (truncated Gaussian) are shown. The S2E calculations are taking into account all relevant collective effects, including a thorough treatment of the impact of the so called Micro-Bunching Instability. For that purpose, tracking with the real number of particles is done, with a resolution in space and time fine enough to calculate instability growth model-free.


title: Factorization Algebras for Bulk-Boundary Systems

abstract: Costello and Gwilliam have given both 1) a general definition of perturbative quantum gauge theory on a manifold M and 2) a construction of a factorization algebra of quantum observables assigned to every quantum gauge theory. In this dissertation, we extend these constructions to a certain general class of field theories on manifolds with boundary.


title: Universal Path Gain Laws for Common Wireless Communication Environments

abstract: Simple and accurate expressions for path gain are derived from electromagnetic fundamentals in a wide variety of common environments, including Line-of-Sight (LOS) and Non-Line-of-Sight (NLOS) indoor urban canyons, urban/rural macro, outdoor-indoor and suburban streets with vegetation. Penetration into a scattering region, sometimes aided by guiding, is the "universal" phenomenon shared by the diverse morphologies. Root Mean Square (RMS) errors against extensive measurements are under 5 dB, better than 3GPP models by 1-12 dB RMS, depending on environment. In urban canyons the models have 4.7 dB RMS error, as compared to 7.9 dB from linear fit to data and 13.9/17.2 dB from LOS/NLOS 3GPP models. The theoretical path gains depend on distance as a power law with exponents from a small set {1.5, 2, 2.5, 4}, specific to each morphology. This provides a theoretical justification for widely used power law empirical models. Only coarse environmental data is needed as parameters: street width, building height, vegetation depth, wall material and antenna heights.


title: Truly Low-Space Element Distinctness and Subset Sum via Pseudorandom Hash Functions

abstract: We consider low-space algorithms for the classic Element Distinctness problem: given an array of $n$ input integers with $O(\log n)$ bit-length, decide whether or not all elements are pairwise distinct. Beame, Clifford, and Machmouchi [FOCS 2013] gave an $\tilde O(n^{1.5})$-time randomized algorithm for Element Distinctness using only $O(\log n)$ bits of working space. However, their algorithm assumes a random oracle (in particular, read-only random access to polynomially many random bits), and it was asked as an open question whether this assumption can be removed. In this paper, we positively answer this question by giving an $\tilde O(n^{1.5})$-time randomized algorithm using $O(\log ^3 n\log \log n)$ bits of space, with one-way access to random bits. As a corollary, we also obtain a $\operatorname{\mathrm{poly}}(n)$-space $O^*(2^{0.86n})$-time randomized algorithm for the Subset Sum problem, removing the random oracles required in the algorithm of Bansal, Garg, Nederlof, and Vyas [STOC 2017]. The main technique underlying our results is a pseudorandom hash family based on iterative restrictions, which can fool the cycle-finding procedure in the algorithms of Beame et al. and Bansal et al.


title: Increasing Liquid State Machine Performance with Edge-of-Chaos Dynamics Organized by Astrocyte-modulated Plasticity

abstract: The liquid state machine (LSM) combines low training complexity and biological plausibility, which has made it an attractive machine learning framework for edge and neuromorphic computing paradigms. Originally proposed as a model of brain computation, the LSM tunes its internal weights without backpropagation of gradients, which results in lower performance compared to multi-layer neural networks. Recent findings in neuroscience suggest that astrocytes, a long-neglected non-neuronal brain cell, modulate synaptic plasticity and brain dynamics, tuning brain networks to the vicinity of the computationally optimal critical phase transition between order and chaos. Inspired by this disruptive understanding of how brain networks self-tune, we propose the neuron-astrocyte liquid state machine (NALSM) that addresses under-performance through self-organized near-critical dynamics. Similar to its biological counterpart, the astrocyte model integrates neuronal activity and provides global feedback to spike-timing-dependent plasticity (STDP), which self-organizes NALSM dynamics around a critical branching factor that is associated with the edge-of-chaos. We demonstrate that NALSM achieves state-of-the-art accuracy versus comparable LSM methods, without the need for data-specific hand-tuning. With a top accuracy of 97.61% on MNIST, 97.51% on N-MNIST, and 85.84% on Fashion-MNIST, NALSM achieved comparable performance to current fully-connected multi-layer spiking neural networks trained via backpropagation. Our findings suggest that the further development of brain-inspired machine learning methods has the potential to reach the performance of deep learning, with the added benefits of supporting robust and energy-efficient neuromorphic computing on the edge.


title: Two neural-network-based methods for solving obstacle problems

abstract: Two neural-network-based numerical schemes are proposed to solve the classical obstacle problems. The schemes are based on the universal approximation property of neural networks, and the cost functions are taken as the energy minimization of the obstacle problems. We rigorously prove the convergence of the two schemes and derive the convergence rates with the number of neurons $N$. In the simulations, we use two example problems (1-D & 2-D) to verify the convergence rate of the methods and the quality of the results.


title: AIRCC-Clim: a user-friendly tool for generating regional probabilistic climate change scenarios and risk measures

abstract: Complex physical models are the most advanced tools available for producing realistic simulations of the climate system. However, such levels of realism imply high computational cost and restrictions on their use for policymaking and risk assessment. Two central characteristics of climate change are uncertainty and that it is a dynamic problem in which international actions can significantly alter climate projections and information needs, including partial and full compliance of global climate goals. Here we present AIRCC-Clim, a simple climate model emulator that produces regional probabilistic climate change projections of monthly and annual temperature and precipitation, as well as risk measures, based both on standard and user-defined emissions scenarios for six greenhouse gases. AIRCC-Clim emulates 37 atmosphere-ocean coupled general circulation models with low computational and technical requirements for the user. This standalone, user-friendly software is designed for a variety of applications including impact assessments, climate policy evaluation and integrated assessment modelling.


title: Modelling COVID-19 Pandemic Dynamics Using Transparent, Interpretable, Parsimonious and Simulatable (TIPS) Machine Learning Models: A Case Study from Systems Thinking and System Identification Perspectives

abstract: Since the outbreak of COVID-19, an astronomical number of publications on the pandemic dynamics appeared in the literature, of which many use the susceptible infected removed (SIR) and susceptible exposed infected removed (SEIR) models, or their variants, to simulate and study the spread of the coronavirus. SIR and SEIR are continuous-time models which are a class of initial value problems (IVPs) of ordinary differential equations (ODEs). Discrete-time models such as regression and machine learning have also been applied to analyze COVID-19 pandemic data (e.g. predicting infection cases), but most of these methods use simplified models involving a small number of input variables pre-selected based on a priori knowledge, or use very complicated models (e.g. deep learning), purely focusing on certain prediction purposes and paying little attention to the model interpretability. There have been relatively fewer studies focusing on the investigations of the inherent time-lagged or time-delayed relationships e.g. between the reproduction number (R number), infection cases, and deaths, analyzing the pandemic spread from a systems thinking and dynamic perspective. The present study, for the first time, proposes using systems engineering and system identification approach to build transparent, interpretable, parsimonious and simulatable (TIPS) dynamic machine learning models, establishing links between the R number, the infection cases and deaths caused by COVID-19. The TIPS models are developed based on the well-known NARMAX (Nonlinear AutoRegressive Moving Average with eXogenous inputs) model, which can help better understand the COVID-19 pandemic dynamics. A case study on the UK COVID-19 data is carried out, and new findings are detailed. The proposed method and the associated new findings are useful for better understanding the spread dynamics of the COVID-19 pandemic.


title: A Harder-Narasimhan stratification of the $B_{dR}^+$-Grassmannian

abstract: We establish a Harder-Narasimhan formalism for modifications of $G$-bundles on the Fargues-Fontaine curve. The semi-stable stratum of the associated stratification of the $B_{dR}^+$-Grassmannian coincides with the variant of the weakly admissible locus defined by Viehmann, and its classical points agree with those of the basic Newton stratum. When restricted to minuscule affine Schubert cells, the stratification corresponds to the Harder-Narasimhan stratification of Dat, Orlik and Rapoport. We also study basic geometric properties of the strata, and the relation to the Hodge-Newton decomposition.


title: Nonsingular Cosmology from an Interacting Vacuum

abstract: We examine the dynamics of FLRW cosmologies in which the vacuum interacts with a perfect fluid through an energy exchange, focusing on the exploration of nonsingular configurations, including cyclic and bouncing models. We consider two specific choices for the energy transfer. In the first case, the energy transfer is proportional to a linear combination of the vacuum and fluid energy densities which makes the conservation equations exactly integrable. The resulting Friedmann equation can be interpreted as an energy constraint equation with an effective potential for the scale factor that may include an infinite barrier forcing a bounce at small values of the scale factor, as well as a potential well allowing for cycling solutions. In the second case, the energy transfer is a nonlinear combination of the vacuum and fluid energy densities. Nonetheless even in this case the dynamics can be partially integrated, leading to a first integral, reducing the number of degrees of freedom. We show that also in this nonlinear case bouncing and cycling cosmologies may arise. In both cases the structure of the resulting phase space allows for nonsingular orbits with an early accelerated phase around a single bounce, connected via a decelerated matter-dominated era to a late-time accelerated phase dominated by an effective cosmological constant.


title: Quantitative unique continuation for Robin boundary value problems on $C^{1,1}$ domains

abstract: In this paper, we prove two unique continuation results for second order elliptic equations with Robin boundary conditions on $C^{1,1}$ domains. The first one is a sharp vanishing order estimate of Robin problems with Lipschitz coefficients and differentiable, sign-changing potentials. This generalizes the result for the "Robin eigenfunctions" in [26], which deals with the case with constant potentials. The second result is a unique continuation result from the boundary -- any non-trivial solution cannot vanish at infinite order from the boundary or vanish on an open subset on the boundary. Such result generalizes the one in [1] for the Laplace equation on $C^{1,1}$ domains with zero Neumann boundary conditions.


title: Regularization for Shuffled Data Problems via Exponential Family Priors on the Permutation Group

abstract: In the analysis of data sets consisting of (X, Y)-pairs, a tacit assumption is that each pair corresponds to the same observation unit. If, however, such pairs are obtained via record linkage of two files, this assumption can be violated as a result of mismatch error rooting, for example, in the lack of reliable identifiers in the two files. Recently, there has been a surge of interest in this setting under the term "Shuffled data" in which the underlying correct pairing of (X, Y)-pairs is represented via an unknown index permutation. Explicit modeling of the permutation tends to be associated with substantial overfitting, prompting the need for suitable methods of regularization. In this paper, we propose a flexible exponential family prior on the permutation group for this purpose that can be used to integrate various structures such as sparse and locally constrained shuffling. This prior turns out to be conjugate for canonical shuffled data problems in which the likelihood conditional on a fixed permutation can be expressed as product over the corresponding (X,Y)-pairs. Inference is based on the EM algorithm in which the intractable E-step is approximated by the Fisher-Yates algorithm. The M-step is shown to admit a significant reduction from $n^2$ to $n$ terms if the likelihood of (X,Y)-pairs has exponential family form as in the case of generalized linear models. Comparisons on synthetic and real data show that the proposed approach compares favorably to competing methods.


title: Nearly Optimal Algorithms for Level Set Estimation

abstract: The level set estimation problem seeks to find all points in a domain ${\cal X}$ where the value of an unknown function $f:{\cal X}\rightarrow \mathbb{R}$ exceeds a threshold $\alpha$. The estimation is based on noisy function evaluations that may be acquired at sequentially and adaptively chosen locations in ${\cal X}$. The threshold value $\alpha$ can either be \emph{explicit} and provided a priori, or \emph{implicit} and defined relative to the optimal function value, i.e. $\alpha = (1-\epsilon)f(x_\ast)$ for a given $\epsilon > 0$ where $f(x_\ast)$ is the maximal function value and is unknown. In this work we provide a new approach to the level set estimation problem by relating it to recent adaptive experimental design methods for linear bandits in the Reproducing Kernel Hilbert Space (RKHS) setting. We assume that $f$ can be approximated by a function in the RKHS up to an unknown misspecification and provide novel algorithms for both the implicit and explicit cases in this setting with strong theoretical guarantees. Moreover, in the linear (kernel) setting, we show that our bounds are nearly optimal, namely, our upper bounds match existing lower bounds for threshold linear bandits. To our knowledge this work provides the first instance-dependent, non-asymptotic upper bounds on sample complexity of level-set estimation that match information theoretic lower bounds.


title: Harvesting BAT-GUANO with NITRATES (Non-Imaging Transient Reconstruction And TEmporal Search): Detecting and localizing the faintest GRBs with a likelihood framework

abstract: The detection of the gravitational wave counterpart GRB 170817A, underluminous compared to the cosmological GRB population by a factor of 10,000, motivates significant effort in detecting and localizing a dim, nearby, and slightly off-axis population of short GRBs. Swift/BAT is the most sensitive GRB detector in operation, and the only one that regularly localizes GRBs to arcminute precision, critical to rapid followup studies. However, the utility of BAT in targeted sub-threshold searches had been historically curtailed by the unavailability of the necessary raw data for analysis. The new availability of time-tagged event (TTE) data from the GUANO system (arXiv:2005.01751), motivates renewed focus on developing sensitive targeted search analysis techniques to maximally exploit these data. While computationally cheap, we show that the typical coded-mask deconvolution imaging is limited in its sensitivity due to several factors. We formalize a maximum likelihood framework for the analysis of BAT data wherein signals are forward modelled through the full instrument response, and -- coupled with the development of new response models -- demonstrate its superior sensitivity to typical imaging via archival comparisons, injection campaigns, and, after implementing as a targeted search, a large number of low-latency GRB discoveries and confirmed arcminute localizations to date. We also demonstrate independent localization of some out-of-FOV GRBs for the first time. NITRATES's increased sensitivity boosts the discovery rate of GRB 170817A-like events in BAT by a factor of at least $3-4$x, along with enabling joint analyses and searches with other GRB, GW, neutrino, and FRB instruments. We provide public access to the response functions and search pipeline code.


title: The Curtis-Wellington spectral sequence through cohomology

abstract: We study stable homotopy through unstable methods applied to its representing infinite loop space, as pioneered by Curtis and Wellington. Using cohomology instead of homology, we find a width filtration whose subquotients are simple quotients of Dickson algebras. We make initial calculations and determine towers in the resulting width spectral sequence. We also make calculations related to the image of $J$ and conjecture that it is captured exactly by the lowest filtration in the width spectral sequence.


title: Multiple generalized cluster structures on $D(\text{GL}_n)$

abstract: We produce a large class of generalized cluster structures on the Drinfeld double of $\text{GL}_n$ that are compatible with Poisson brackets given by Belavin-Drinfeld classification. The resulting construction is compatible with the previous results on cluster structures on $\text{GL}_n$.


title: A Temporal Scale to Track the Spectral Transitions in Low-Mass X-ray Binaries

abstract: The results of a temporal analysis of observations for a sample of nine low-mass X-ray binaries (LMXBs) are presented. Of these sources, five host a neutron star (NS) primary (4U1608-52, Aql X-1, 4U1705-44, GX17+2 and Cyg X-2), and four host a black hole (BH) (GX339-4, XTE J1859+226, H1743-322 and MAXI J1659-152). The NS group includes three Atolls and two Z-type sources. We utilized archival Proportional Counter Array (PCA)/RXTE data to construct high-resolution lightcurves. A wavelet transform of the lightcurves is deployed to extract a minimal time scale (MTS) associated with the spectral state of the sources. The MTS, together with the fractional root-mean-square (RMS) and hardness ratios, is used to construct RMS-MTS and hardness-MTS diagrams that enable a direct comparison of the evolution of spectral transitions in the target sources. Observations with high fractional RMS and high hardness cluster in a broad region occupied jointly by BH and NS sources. For low fractional RMS observations the Atolls exhibit large MTS whereas Z-type sources exhibit small MTS. This new feature raises the possibility of discriminating between these two types of sources. Moreover, in the hardness-MTS plane, BH sources are the sole occupiers of the low-hardness and small-MTS domain thus potentially signaling a unique property for distinguishing BH and NS hosts in LMXBs.


title: Data-Driven System Identification of 6-DoF Ship Motion in Waves with Neural Networks

abstract: Critical evaluation and understanding of ship responses in the ocean is important for not only the design and engineering of future platforms but also the operation and safety of those that are currently deployed. Simulations or experiments are typically performed in nominal sea conditions during ship design or prior to deployment and the results may not be reflective of the instantaneous state of the vessel and the ocean environment while deployed. Short-term temporal predictions of ship responses given the current wave environment and ship state would enable enhanced decision-making onboard for both manned and unmanned vessels. However, the current state-of-the-art in numerical hydrodynamic simulation tools are too computationally expensive to be employed for real-time ship motion forecasting and the computationally efficient tools are too low fidelity to provide accurate responses. A methodology is developed with long short-term memory (LSTM) neural networks to represent the motions of a free running David Taylor Model Basin (DTMB) 5415 destroyer operating at 20 knots in Sea State 7 stern-quartering irregular seas. Case studies are performed for both course-keeping and turning circle scenarios. An estimate of the vessel's encounter frame is made with the trajectories observed in the training dataset. Wave elevation time histories are given by artificial wave probes that travel with the estimated encounter frame and serve as input into the neural network, while the output is the 6-DOF temporal ship motion response. Overall, the neural network is able to predict the temporal response of the ship due to unseen waves accurately, which makes this methodology suitable for system identification and real-time ship motion forecasting. The methodology, the dependence of model accuracy on wave probe and training data quantity and the estimated encounter frame are all detailed.


title: Emergence and structure of decentralised trade networks around dark web marketplaces

abstract: Dark web marketplaces (DWMs) are online platforms that facilitate illicit trade among millions of users generating billions of dollars in annual revenue. Recently, two interview-based studies have suggested that DWMs may also promote the emergence of direct user-to-user (U2U) trading relationships. Here, we quantify the scale of, and thoroughly investigate, U2U trading around DWMs by analysing 31 million Bitcoin transactions among users of 40 DWMs between June 2011 and Jan 2021. We find that half of the DWM users trade through U2U pairs generating a total trading volume greater than DWMs themselves. We then show that hundreds of thousands of DWM users form stable trading pairs that are persistent over time. Users in stable pairs are typically the ones with the largest trading volume on DWMs. Then, we show that new U2U pairs often form while both users are active on the same DWM, suggesting the marketplace may serve as a catalyst for new direct trading relationships. Finally, we reveal that stable U2U pairs tend to survive DWM closures and that they were not affected by COVID-19, indicating that their trading activity is resilient to external shocks. Our work unveils sophisticated patterns of trade emerging in the dark web and highlights the importance of investigating user behaviour beyond the immediate buyer-seller network on a single marketplace.


title: Hadronic structure on the light-front II: QCD strings, Wilson lines and potentials

abstract: This is the second paper on hadronic wave functions in the light front formulation. We first consider only confinement effects, using the classical Nambu-Gotto string with massive end points in the light cone gauge. We derive the light-front Hamiltonian and show how to solve it, using an expansion in suitable basis functions. We next discuss the correlators of Wilson lines on the light front, leading to an effective Hamiltonian that includes spin effects. At the end, we consider a separate problem of instanton-induced interactions, at the origin of the pion as a massless Goldstone mode on the light front.


title: A Time and Cost Effective Method for Entropic Coefficient Determination of a Large Commercial Battery Cell

abstract: The entropic coefficient of a lithium-ion battery cell is used to calculate the reversible heat of a battery during operation, which is a non-negligible part for the battery thermal modelling. The contribution of this article is to propose a novel method to establish the entropic coefficient profile of a 26 Ah commercial pouch cell, and compare the results with those obtained from the traditional potentiometric and calorimetric methods, and all are found to be in a good agreement. The originality of this work is to use a method which consists of supplying a square pulse current wave form at a certain frequency and thus the resulting heat variation could be successfully linked to the input current using Fourier analysis. The current magnitude used were 1 C and 1.5 C, which are representative for the normal operation current in an electrified vehicle application. The method proposed is found to be cost efficient with a short experiment time and simple experiment setup. In fact, it can be used to characterize cells that are already mounted in a pack without the access to a climate chamber or calorimeter.


title: A Framework for Real-World Multi-Robot Systems Running Decentralized GNN-Based Policies

abstract: GNNs are a paradigm-shifting neural architecture to facilitate the learning of complex multi-agent behaviors. Recent work has demonstrated remarkable performance in tasks such as flocking, multi-agent path planning and cooperative coverage. However, the policies derived through GNN-based learning schemes have not yet been deployed to the real-world on physical multi-robot systems. In this work, we present the design of a system that allows for fully decentralized execution of GNN-based policies. We create a framework based on ROS2 and elaborate its details in this paper. We demonstrate our framework on a case-study that requires tight coordination between robots, and present first-of-a-kind results that show successful real-world deployment of GNN-based policies on a decentralized multi-robot system relying on Adhoc communication. A video demonstration of this case-study, as well as the accompanying source code repository, can be found online. https://www.youtube.com/watch?v=COh-WLn4iO4 https://github.com/proroklab/ros2_multi_agent_passage https://github.com/proroklab/rl_multi_agent_passage


title: Location inference on social media data for agile monitoring of public health crises: An application to opioid use and abuse during the Covid-19 pandemic

abstract: The Covid-19 pandemic has intersected with the opioid epidemic to create a unique public health crisis, with the health and economic consequences of the virus and associated lockdowns compounding pre-existing social and economic stressors associated with rising opioid and heroin use and abuse. In order to better understand these interlocking crises, we use social media data to extract qualitative and quantitative insights on the experiences of opioid users during the Covid-19 pandemic. In particular, we use an unsupervised learning approach to create a rich geolocated data source for public health surveillance and analysis. To do this we first infer the location of 26,000 Reddit users that participate in opiate-related sub-communities (subreddits) by combining named entity recognition, geocoding, density-based clustering, and heuristic methods. Our strategy achieves 63 percent accuracy at state-level location inference on a manually-annotated reference dataset. We then leverage the geospatial nature of our user cohort to answer policy-relevant questions about the impact of varying state-level policy approaches that balance economic versus health concerns during Covid-19. We find that state government strategies that prioritized economic reopening over curtailing the spread of the virus created a markedly different environment and outcomes for opioid users. Our results demonstrate that geospatial social media data can be used for agile monitoring of complex public health crises.


title: Quantum signatures in nonlinear gravitational waves

abstract: The effective quantum field theory description of gravity, despite its non-renormalizability, allows for predictions beyond classical general relativity. As we enter the age of gravitational wave astronomy, an important and timely question is whether measurable quantum predictions that depart from classical gravity, analogous to quantum optics effects which cannot be explained by classical electrodynamics, can be found. In this work, we investigate quantum signatures in gravitational waves using tools from quantum optics. Squeezed-coherent gravitational waves, which can exhibit sub-Poissonian graviton statistics, can enhance or suppress the signal measured by an interferometer, a characteristic effect of quantum squeezing. Moreover, we show that Gaussian gravitational wave quantum states can be reconstructed from measurements over an ensemble of optical fields interacting with a single copy of the gravitational wave, thus opening the possibility of detecting quantum features of gravity beyond classical general relativity.


title: Game of Life on Graphs

abstract: We consider a specific graph dynamical system inspired by the famous Conway's Game of Life in this work. We study the properties of the dynamical system on different graphs and introduce a new efficient heuristic for graph isomorphism testing. We use the evolution of our system to extract features from a graph in a deterministic way and observe that the extracted features are unique and the distance induced by that features satisfy triangle inequality for all connected graphs with up to ten vertices.


title: Analysis and Numerical Solution of a Modular Convex Nash Equilibrium Problem

abstract: We investigate a modular convex Nash equilibrium problem involving nonsmooth functions acting on linear mixtures of strategies, as well as smooth coupling functions. An asynchronous block-iterative decomposition method is proposed to solve it.


title: Improving the Cook et al. Proximity Bound Given Integral Valued Constraints

abstract: Consider a linear program of the form $\max\;c^{\top}x:Ax\leq b$, where $A$ is an $m\times n$ integral matrix. In 1986 Cook, Gerards, Schrijver, and Tardos proved that, given an optimal solution $x^{*}$, if an optimal integral solution $z^{*}$ exists, then it may be chosen such that $\left\Vert x^{*}-z^{*}\right\Vert _{\infty}<n\Delta$, where $\Delta$ is the largest magnitude of any subdeterminant of $A$. Since then an open question has been to improve this bound, assuming that $b$ is integral valued too. In this manuscript we show that $n\Delta$ can be replaced with $\frac{n}{2}\cdot\Delta$ whenever $n\geq2$. We also show that, in certain circumstances, the factor $n$ can be removed entirely.


title: Optimal bailout strategies resulting from the drift controlled supercooled Stefan problem

abstract: We consider the problem faced by a central bank which bails out distressed financial institutions that pose systemic risk to the banking sector. In a structural default model with mutual obligations, the central agent seeks to inject a minimum amount of cash in order to limit defaults to a given proportion of entities. We prove that the value of the central agent's control problem converges as the number of defaultable institutions goes to infinity, and that it satisfies a drift controlled version of the supercooled Stefan problem. We compute optimal strategies in feedback form by solving numerically a regularized version of the corresponding mean field control problem using a policy gradient method. Our simulations show that the central agent's optimal strategy is to subsidise banks whose equity values lie in a non-trivial time-dependent region.


title: Towards the 5/6-Density Conjecture of Pinwheel Scheduling

abstract: Pinwheel Scheduling aims to find a perpetual schedule for unit-length tasks on a single machine subject to given maximal time spans (a.k.a. frequencies) between any two consecutive executions of the same task. The density of a Pinwheel Scheduling instance is the sum of the inverses of these task frequencies; the 5/6-Conjecture (Chan and Chin, 1993) states that any Pinwheel Scheduling instance with density at most 5/6 is schedulable. We formalize the notion of Pareto surfaces for Pinwheel Scheduling and exploit novel structural insights to engineer an efficient algorithm for computing them. This allows us to (1) confirm the 5/6-Conjecture for all Pinwheel Scheduling instances with at most 12 tasks and (2) to prove that a given list of only 23 schedules solves all schedulable Pinwheel Scheduling instances with at most 5 tasks.


title: PatchGame: Learning to Signal Mid-level Patches in Referential Games

abstract: We study a referential game (a type of signaling game) where two agents communicate with each other via a discrete bottleneck to achieve a common goal. In our referential game, the goal of the speaker is to compose a message or a symbolic representation of "important" image patches, while the task for the listener is to match the speaker's message to a different view of the same image. We show that it is indeed possible for the two agents to develop a communication protocol without explicit or implicit supervision. We further investigate the developed protocol and show the applications in speeding up recent Vision Transformers by using only important patches, and as pre-training for downstream recognition tasks (e.g., classification). Code available at https://github.com/kampta/PatchGame.


title: A Recommendation System to Enhance Midwives' Capacities in Low-Income Countries

abstract: Maternal and child mortality is a public health problem that disproportionately affects low- and middle-income countries. Every day, 800 women and 6,700 newborns die from complications related to pregnancy or childbirth. And for every maternal death, about 20 women suffer serious birth injuries. However, nearly all of these deaths and negative health outcomes are preventable. Midwives are key to revert this situation, and thus it is essential to strengthen their capacities and the quality of their education. This is the aim of the Safe Delivery App, a digital job aid and learning tool to enhance the knowledge, confidence and skills of health practitioners. Here, we use the behavioral logs of the App to implement a recommendation system that presents each midwife with suitable contents to continue gaining expertise. We focus on predicting the click-through rate, the probability that a given user will click on a recommended content. We evaluate four deep learning models and show that all of them produce highly accurate predictions.


title: Flatfield Ultrafast Imaging with Single-Shot Non-Synchronous Array Photography

abstract: We present a method for acquiring a sequence of time-resolved images in a single shot, called Single-Shot Non-Synchronous Array Photography (SNAP). In SNAP, a pulsed laser beam is split by a diffractive optical element into an array of angled beamlets whose illumination fronts remain perpendicular to the optical axis. Different time delays are imparted to each beamlet by an echelon, enabling them to probe a scene in rapid succession. The beamlets are then imaged onto different regions of a camera by a lenslet array. Because the illumination fronts remain flat (head-on) independently of beamlet angle, the temporal resolution of SNAP is fundamentally limited only by the laser pulse duration, akin to a "global shutter" in conventional imaging. We demonstrate SNAP by capturing the evolution of a laser induced plasma filament for 20 frames at an average rate of 4.2 Tfps.


title: Wide binaries from the H3 survey: the thick disk and halo have similar wide binary fractions

abstract: Due to the different environments in the Milky Way's disk and halo, comparing wide binaries in the disk and halo is key to understanding wide binary formation and evolution. By using Gaia Early Data Release 3, we search for resolved wide binary companions in the H3 survey, a spectroscopic survey that has compiled $\sim$150,000 spectra for thick-disk and halo stars to date. We identify 800 high-confidence (a contamination rate of 4%) wide binaries and two resolved triples, with binary separations mostly between $10^3$-$10^5$ AU and a lowest [Fe/H] of $-2.7$. Based on their Galactic kinematics, 33 of them are halo wide binaries, and most of those are associated with the accreted Gaia-Sausage-Enceladus galaxy. The wide binary fraction in the thick disk decreases toward the low metallicity end, consistent with the previous findings for the thin disk. Our key finding is that the halo wide binary fraction is consistent with the thick-disk stars at a fixed [Fe/H]. There is no significant dependence of the wide binary fraction on the $\alpha$-captured abundance. Therefore, the wide binary fraction is mainly determined by the iron abundance, not their disk or halo origin nor the $\alpha$-captured abundance. Our results suggest that the formation environments play a major role for the wide binary fraction, instead of other processes like radial migration that only apply to disk stars.


title: The eccentricity distribution of wide binaries and their individual measurements

abstract: Eccentricity of wide binaries is difficult to measure due to their long orbital periods. With Gaia's high-precision astrometric measurements, eccentricity of a wide binary can be constrained by the angle between the separation vector and the relative velocity vector (the $v$-$r$ angle). In this paper, by using the $v$-$r$ angles of wide binaries in Gaia Early Data Release 3, we develop a Bayesian approach to measure the eccentricity distribution as a function of binary separations. Furthermore, we infer the eccentricities of individual wide binaries and make them publicly available. Our results show that the eccentricity distribution of wide binaries at $10^2$ AU is close to uniform and becomes superthermal at $>10^{3}$ AU, suggesting two formation mechanisms dominating at different separation regimes. The close binary formation, most likely disk fragmentation, results in a uniform eccentricity distribution at $<10^{2}$ AU. The wide binary formation that leads to highly eccentric wide binaries at $>10^{3}$ AU may be turbulent fragmentation and/or the dynamical unfolding of compact triples. With Gaia, measuring eccentricities is now possible for a large number of wide binaries, opening a new window to understanding binary formation and evolution.


title: M5-branes Probing Flux Backgrounds

abstract: We analyze the global symmetries and anomalies of 4d $\mathcal{N} = 1$ field theories that arise from a stack of $N$ M5-branes probing a class of flux backgrounds. These backgrounds consist of a resolved $\mathbb{C}^2 / \mathbb{Z}_k$ singularity fibered over a smooth Riemann surface of genus $g \geq 2$, supported by a non-trivial $G_4$-flux configuration labeled by a collection of $2(k-1)$ flux quanta, $\{N_i\}$. For $k=2$, this setup defines a non-trivial superconformal field theory (SCFT) in the IR, which is holographically dual to an explicit $AdS_5$ solution first described by Gauntlett, Martelli, Sparks, and Waldram. The generalization to $k \geq 3$ is hard to tackle directly within holography. Instead, in this paper we lay the groundwork for a systematic analysis of such a generalization by adopting anomaly inflow methods to identify continuous and discrete global symmetries of the 4d field theories. We also compute the 't Hooft anomalies for continuous symmetries at leading order in the limit of large $N$, $N_i$.


title: Neutron-Mirror Neutron oscillations in Matter

abstract: The possibility that a neutron can be transformed to a hidden sector particle remains intriguingly open. Proposed theoretical models conjecture that the hidden sector can be represented by a mirror sector, and the neutron n can oscillate into its sterile mirror twin n', exactly or nearly degenerate in mass with n. Oscillations n - n' can take place in vacuum and in the environment of the regular matter and the magnetic field where only neutron will be subject of interaction with the environment. We describe the propagation of the oscillating n - n' system as a particle of the cold neutron beam passing through the dense absorbing materials in connection with the possible regeneration type of experiments where the effect of n -> n' -> n transformation can be observed.


title: Ultrahigh-sensitivity optical power monitor for Si photonic circuits

abstract: A phototransistor is a promising candidate as an optical power monitor in Si photonic circuits since the internal gain of photocurrent enables high sensitivity. However, state-of-the-art waveguide-coupled phototransistors suffer from a responsivity of lower than $10^3$ A/W, which is insufficient for detecting very low power light. Here, we present a waveguide-coupled phototransistor consisting of an InGaAs ultrathin channel on a Si waveguide working as a gate electrode to increase the responsivity. The Si waveguide gate underneath the InGaAs ultrathin channel enables the effective control of transistor current without optical absorption by the gate metal. As a result, our phototransistor achieved the highest responsivity of approximately $10^6$ A/W among the waveguide-coupled phototransistors, allowing us to detect light of 621 fW propagating in the Si waveguide. The high responsivity and the reasonable response time of approximately 100 $\mu$s make our phototransistor promising as an effective optical power monitor in Si photonics circuits.


title: Effect of He Self-organized pattern plasma-activated media with different conductivity on cancer cells

abstract: The self-organized pattern (SOP) phenomenon is prevalent in plasma, while knowledge about SOP discharge affecting reactive species generated plasma-activated media (PAM) for cancer therapy is poorly documented. The aim of this study focused on the effect of SOP discharge modes on reactive oxygen and nitrogen species (ROS, RNS) in He SOP plasma-activated media with different conductivity (saline solution and deionized (DI) water), and employed them to breast cancer MDA-MB-231 and pancreatic BxPC-3 cancer cells. Optical emission spectrum and Fluorimetric analysis were used to identify and quantify ROS and RNS generated in He SOP plasma-activated saline solution and DI water. Furthermore, He SOP plasma discharge modes are capable of efficiently controlling the ROS and RNS concentration in the plasma-activated saline solution and DI water, which contribute to the cytotoxic effect. On the other hand, stainless steel and copper were used as a lower electrode to compare their effect on cell viability. Taken together, our findings provide insight into potential mechanisms involved in cell death after treatment with He SOP plasma-activated media.


title: Experimental and computational studies on the reactivity of methanimine radical cation (H$_2$CNH$^{+\cdot}$) and its isomer aminomethylene (HCNH$_2^{+\cdot}$) with C$_2$H$_2$

abstract: Experimental and theoretical studies are presented on the reactivity of the radical cation isomers methanimine and aminomethylene with ethyne. Selective isomer generation is performed via dissociative photoionization of suitable neutral precursors and via direct photoionization of methanimine. Reactive cross sections and product branching ratios are measured as a function of photon and collision energies. Results are discussed in light of ab initio calculations of reaction mechanisms. The major channels, for both isomers, are due to H atom elimination from covalently bound adducts to give [C3NH4]+. Theoretical calculations show that while for the reaction of aminomethylene with acetylene any of the three lowest energy [C3NH4]+ isomers can form via barrierless and exothermic pathways, for the methanimine reagent the only barrierless pathway is the one leading to the production of protonated vinyl cyanide (CH2CHCNH+), a prototypical branched nitrile species that has been proposed as a likely intermediate in star forming regions and in the atmosphere of Titan. The astrochemical implications of the results are briefly addressed.

categories: physics.chem-ph astro-ph.EP astro-ph.GA astro-ph.SR


title: X-PIPS-M Data Summary

abstract: A compilation of descriptive statistics of responses to the X-PIPS-M suite of surveys for students, student instructors, and faculty instructors of introductory mathematics courses. Data were obtained from Precalculus, Calculus 1, and Calculus 2 courses at 12 universities during the 2017-18 and 2018-19 academic years, as part of the Progress through Calculus research project (NSF DUE IUSE #1430540).


title: A Helioscope for Gravitationally Bound Millicharged Particles

abstract: Particles may be emitted efficiently from the solar interior if they are sufficiently light and weakly coupled to the solar plasma. In a narrow region of phase space, they are emitted with velocities smaller than the escape velocity of the solar system, thereby populating a gravitationally bound density that can accumulate over the solar lifetime, referred to as a "solar basin." Detection strategies that can succeed in spite of (or even be enhanced by) the low particle velocities are therefore poised to explore new regions of parameter space when taking this solar population into account. Here we identify "direct deflection" as a powerful method to detect such a population of millicharged particles. This approach involves distorting the local flow of gravitationally bound millicharges with an oscillating electromagnetic field and measuring these distortions with a resonant LC circuit. Since it is easier to distort the flow of slowly moving particles, the signal is parametrically enhanced by the small solar escape velocity near Earth. The proposed setup can probe couplings an order of magnitude smaller than other methods for millicharge masses ranging from 100 meV to 100 eV and can operate concurrently as a search for sub-GeV millicharged dark matter. The signal power scales as the millicharge coupling to the eighth power, meaning that even with conservative assumptions, direct deflection could begin to explore new regions of parameter space. We also highlight novel features of millicharge solar basins, including those associated with the phase space distribution and the possibility for the occupation number to vastly exceed that of a thermal distribution.


title: Demonstration of MHz frequency domain multiplexing readout of 37 transition edge sensors for high-resolution x-ray imaging spectrometers

abstract: We report on the development and demonstration of a MHz frequency domain multiplexing (FDM) technology to read out arrays of cryogenic transition edge sensor (TES) X-ray microcalorimeters. In our FDM scheme, TESs are AC-biased at different resonant frequencies in the low MHz range through an array of high-$Q$ LC resonators. The current signals of all TESs are summed at superconducting quantum interference devices (SQUIDs). We have demonstrated multiplexing for a readout of 31 pixels using room temperature electronics, high-$Q$ LC filters and TES arrays developed at SRON, and SQUID arrays from VTT. We repeated this on a second setup with 37 pixels. The summed X-ray spectral resolutions $@$ 5.9 keV are $\Delta E_{\rm 31 pix ~MUX}=2.14\pm0.03$ eV and $\Delta E_{\rm 37 pix ~MUX}=2.23\pm0.03$ eV. The demonstrated results are comparable with other multiplexing approaches. There is potential to further improve the spectral resolution and to increase the number of multiplexed TESs, and to open up applications for TES X-ray microcalorimeters.


title: Turbulent Disk Viscosity and the Bifurcation of Planet Formation Histories

abstract: ALMA observations of dust ring/gap structures in a minority but growing sample of protoplanetary disks can be explained by the presence of planets at large disk radii - yet the origins of these planets remains debated. We perform planet formation simulations using a semi-analytic model of the HL Tau disk to follow the growth and migration of hundreds of planetary embryos initially distributed throughout the disk, assuming either a high or low turbulent $\alpha$ viscosity. We have discovered that there is a bifurcation in the migration history of forming planets as a consequence of varying the disk viscosity. In our high viscosity disks, inward migration prevails and yields compact planetary systems, tempered only by planet trapping at the water iceline around 5 au. In our lower viscosity models however, low mass planets can migrate outward to twice their initial orbital radii, driven by a radially extended region of strong outward-directed corotation torques located near the heat transition (where radiative heating of the disk by the star is comparable to viscous heating) - before eventually migrating inwards. We derive analytic expressions for the planet mass at which the corotation torque dominates, and find that this "corotation mass" scales as $M_{\rm p, corot} \sim \alpha^{2/3}$. If disk winds dominate the corotation torque, the corotation mass scales linearly with wind strength. We propose that the observed bifurcation in disk demographics into a majority of compact dust disks and a minority of extended ring/gap systems is a consequence of a distribution of viscosity across the disk population.


title: Conformal Bootstrap Deformations

abstract: We explore the space of extremal functionals in the conformal bootstrap. By recasting the bootstrap problem as a set of non-linear equations parameterized by the CFT data, we find an efficient algorithm for converging to the extremal solution corresponding to the boundary of allowed regions in the parameter space of CFTs. Furthermore, by deforming these solutions, we demonstrate that certain solutions corresponding to known theories are continuously connected. Employing these methods, we will explore the space of non-unitary CFTs in the context of modular as well as correlation function bootstrap. In two dimensions, we show that the extremal solution corresponding to the Ising model is connected to that of the Yang-Lee minimal model. By deforming this solution to three dimensions, we provide evidence that the CFT data obtained in this way is compatible with the $\epsilon$-expansion for a non-unitary theory.


title: Search for an Axion-Like Particle in B Meson Decays

abstract: Axion-like particles (ALPs) are predicted in many extensions of the Standard Model, and their masses can naturally be well below the electroweak scale. In the presence of couplings to electroweak bosons, these particles could be emitted in flavor-changing $B$ meson decays. We report herein a search for an ALP, $a$, in the reaction $B^\pm\rightarrow K^\pm a$, $a\rightarrow\gamma\gamma$ using data collected by the BABAR experiment at SLAC. No significant signal is observed, and 90% confidence level upper limits on the ALP coupling to electroweak bosons are derived as a function of ALP mass, improving current constraints by several orders of magnitude in the range $0.175\,\,\mathrm{GeV} < m_a < 4.78\,\,\mathrm{GeV}$.


title: Blandford-Znajek jets in galaxy formation simulations: exploring the diversity of outflows produced by spin-driven AGN jets in Seyfert galaxies

abstract: Recent observations of Seyfert galaxies indicate that low power, misaligned jets can undergo significant interaction with the gas in the galactic disc and may be able to drive large-scale, multiphase outflows. We apply our novel sub-grid model for Blandford-Znajek jets to simulations of the central regions of Seyferts, in which a black hole is embedded in a dense, sub-kpc circumnuclear disc (CND) and surrounded by a dilute circumgalactic medium (CGM). We find that the variability of the accretion flow is highly sensitive both to the jet power and to the CND thermodynamics and, ultimately, is determined by the complex interplay between jet-driven outflows and backflows. Even at moderate Eddington ratios, AGN jets are able to significantly alter the thermodynamics and kinematics of CNDs and entrain up to 10% of their mass in the outflow. Mass outflow rates and kinetic powers of the warm outflowing component are in agreement with recent observations for black holes with similar bolometric luminosities, with outflow velocities that are able to reach 500 km/s. Depending on their power and direction, jets are able to drive a wide variety of large-scale outflows, ranging from light, hot and collimated structures to highly mass-loaded, multiphase, bipolar winds. This diversity of jet-driven outflows highlights the importance of applying physically motivated models of AGN feedback to realistic galaxy formation contexts. Such simulations will play a crucial role in accurately interpreting the wealth of data that next generation facilities such as JWST, SKA and Athena will provide.


title: Supermassive black holes in cosmological simulations II: the AGN population and predictions for upcoming X-ray missions

abstract: In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modeling of feedback from active galactic nuclei (AGN). The amount of energy released by AGN feedback is proportional to the mass that has been accreted onto the BHs, but the exact sub-grid modeling of AGN feedback differs in all simulations. Whilst modern simulations reliably produce populations of quiescent massive galaxies at z<2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN population of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z=0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGN of L_{x, 2-10 keV}~10^43-10^44 erg/s. AGN feedback in some simulations prevents the existence of any bright AGN with L_{x, 2-10 keV}>=10^45 erg/s (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z<=2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous (L_{x, 2-10 keV}~10^43-10^45 erg/s) and fainter (L_{x, 2-10 keV}~10^42-10^43 erg/s) AGN, and at both low and high redshift. These differences can aid us in improving future BH and galaxy subgrid modeling in simulations. Upcoming X-ray missions (e.g., Athena, AXIS, and LynX) will bring faint AGN to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGN in future surveys spans at least one order of magnitude across the simulations, at any redshift.


title: Dynamical Mass of the Young Substellar Companion HD 984 B

abstract: Model-independent masses of substellar companions are critical tools to validate models of planet and brown dwarf cooling, test their input physics, and determine the formation and evolution of these objects. In this work, we measure the dynamical mass and orbit of the young substellar companion HD 984 B. We obtained new high-contrast imaging of the HD 984 system with Keck/NIRC2 which expands the baseline of relative astrometry from 3 to 8 years. We also present new radial velocities of the host star with the Habitable-Zone Planet Finder spectrograph at the Hobby-Eberly Telescope. Furthermore, HD 984 exhibits a significant proper motion difference between Hipparcos and Gaia EDR3. Our joint orbit fit of the relative astrometry, proper motions, and radial velocities yields a dynamical mass of $61 \pm 4$ $\mathrm{M_{Jup}}$ for HD 984 B, placing the companion firmly in the brown dwarf regime. The new fit also reveals a higher eccentricity for the companion ($e = 0.76 \pm 0.05$) compared to previous orbit fits. Given the broad age constraint for HD 984, this mass is consistent with predictions from evolutionary models. HD 984 B's dynamical mass places it among a small but growing list of giant planet and brown dwarf companions with direct mass measurements.


title: HI and H$_2$ gas evolution over cosmic times: ColdSIM

abstract: We present first results of the evolution of cold cosmic gas obtained through a set of state-of-the-art numerical simulations (ColdSIM). We model time-dependent atomic and molecular non-equilibrium chemistry coupled to HI and H$_2$ self-shielding, various UV backgrounds as suggested by the recent literature, H$_2$ dust grain catalysis, photoelectric heating, cosmic-ray heating, as well as hydrodynamics, star formation and feedback effects. By means of such non-equilibrium calculations we are finally able to reproduce the latest HI and H$_2$ observational data. The neutral-gas mass density parameter results around $\Omega_{\rm neutral} \!\sim\! 10^{-3}$ and increases from lower to higher redshift ($z$). The molecular-gas mass density parameter shows peak values of $ \Omega_{\rm H_2} \! \sim \! 10^{-4}$, while expected H$_2$ fractions can be as high as 50% of the cold gas mass at $ z\!\sim$4-8, in line with the latest high-$z$ measurements. Both observed HI and H$_2$ trends are well reproduced by our non-equilibrium H$_2$-based star formation modelling. H$_2$ depletion times remain below the Hubble time and comparable to the dynamical time at all epochs. These findings suggest that, besides HI, non-equilibrium H$_2$ analyses are key probes for assessing the cold gas and the role of UV background radiation. Abridged.


title: Giant proximity exchange and flat Chern band in 2D magnet-semiconductor heterostructures

abstract: Van der Waals (vdW) heterostructures formed by two-dimensional magnets and semiconductors have provided a fertile ground for fundamental science and for spintronics. We present first-principles calculations finding a proximity exchange splitting of 14 meV equivalent to an effective Zeeman field of 120 T in the vdW magnet-semiconductor heterostructure MoS$_2$/CrBr$_3$, leading to a 2D spin-polarized half-metal with carrier densities ranging up to $10^{13}$ cm$^{-2}$. We consequently explore the effect of large exchange coupling on the electronic bandstructure when the magnetic layer hosts chiral spin textures such as skyrmions. A flat Chern band is found at a "magic" value of magnetization $\overline{m} \sim 0.2$ for Schr\"odinger electrons, and it generally occurs for Dirac electrons. The magnetic proximity induced anomalous Hall effect enables transport-based detection of chiral spin textures, and flat Chern bands provides an avenue for engineering various strongly correlated states.


title: Non-circular flows in HIghMass galaxies in a test of the late accretion hypothesis

abstract: We present H-alpha velocity maps for the HIghMass galaxies UGC 7899, UGC 8475, UGC 9037 and UGC 9334, obtained with the SITELLE Imaging Fourier Transform Spectrometer on the Canada-France-Hawaii Telescope, to search for kinematic signatures of late gas accretion to explain their large atomic gas reservoirs. The maps for UGC 7899, UGC 9037, and UGC 9334 are amenable to disk wide radial flow searches with the DiskFit algorithm, and those for UGC 7899 and UGC 9037 are also amenable to inner-disk kinematic analyses. We find no evidence for outer disk radial flows down to Vr ~ 20 km/s in UGC 9037 and UGC 9334, but hints of such flows in UGC 7899. Conversely, we find clear signatures of inner (r ~ 5 kpc) noncircularities in UGC 7899 and UGC 9037 that can be modelled as either bisymmetric (which could be produced by a bar) or radial flows. Comparing these models to the structure implied by photometric disk-bulge-bar decompositions, we favour inner radial flows in UGC 7899 and an inner bar in UGC 9037. With hints of outer disk radial flows and an outer disk warp, UGC 7899 is the best candidate for late accretion among the galaxies examined, but additional modelling is required to disentangle potential degeneracies between these signatures in H I and H-alpha velocity maps. Our search provides only weak = constraints on hot-mode accretion models that could explain the unusually high H I content of HIghMass galaxies.


title: Identification of single spectral lines in large spectroscopic surveys using UMLAUT: an Unsupervised Machine Learning Algorithm based on Unbiased Topology

abstract: The identification of an emission line is unambiguous when multiple spectral features are clearly visible in the same spectrum. However, in many cases, only one line is detected, making it difficult to correctly determine the redshift. We developed a freely available unsupervised machine-learning algorithm based on unbiased topology (UMLAUT) that can be used in a very wide variety of contexts, including the identification of single emission lines. To this purpose, the algorithm combines different sources of information, such as the apparent magnitude, size and color of the emitting source, and the equivalent width and wavelength of the detected line. In each specific case, the algorithm automatically identifies the most relevant ones (i.e., those able to minimize the dispersion associated with the output parameter). The outputs can be easily integrated into different algorithms, allowing us to combine supervised and unsupervised techniques and increasing the overall accuracy. We tested our software on WISP (WFC3 IR Spectroscopic Parallel) survey data. WISP represents one of the closest existing analogs to the near-IR spectroscopic surveys that are going to be performed by the future Euclid and Roman missions. These missions will investigate the large-scale structure of the universe by surveying a large portion of the extragalactic sky in near-IR slitless spectroscopy, detecting a relevant fraction of single emission lines. In our tests, UMLAUT correctly identifies real lines in 83.2% of the cases. The accuracy is slightly higher (84.4%) when combining our unsupervised approach with a supervised approach we previously developed.


title: The transversity parton distribution function of the nucleon using the pseudo-distribution approach

abstract: We present a determination of the non-singlet transversity parton distribution function (PDF) of the nucleon, normalized with respect to the tensor charge at $\mu^2=2$ GeV$^2$ from lattice quantum chromodynamics. We apply the pseudo-distribution approach, using a gauge ensemble with a lattice spacing of 0.094 fm and the light quark mass tuned to a pion mass of 358 MeV. We extract the transversity PDF from the analysis of the short-distance behavior of the Ioffe-time pseudo-distribution using the leading-twist next-to-leading order (NLO) matching coefficients calculated for transversity. We reconstruct the $x$-dependence of the transversity PDF through an expansion in a basis of Jacobi polynomials in order to reduce the PDF ansatz dependence. Within the limitations imposed by a heavier-than-physical pion mass and a fixed lattice spacing, we present a comparison of our estimate for the valence transversity PDF with the recent global fit results based on single transverse spin asymmetry. We find the intrinsic nucleon sea to be isospin symmetric with respect to transversity.


title: Neutron-capture elements record the ordered chemical evolution of the disc over time

abstract: An ensemble of chemical abundances probing different nucleosynthetic channels can be leveraged to build a comprehensive understanding of the chemical and structural evolution of the Galaxy. Using GALAH DR3 data, we seek to trace the enrichment by the supernovae Ia, supernovae II, asymptotic giant branch stars, and neutron-star mergers and/or collapsars nucleosynthetic sources by studying the [Fe/H], [$\alpha$/Fe], [Ba/Fe], and [Eu/Fe] chemical compositions of $\sim$50,000 red giant stars, respectively. Employing small [Fe/H]-[$\alpha$/Fe] cells, which serve as an effective reference-frame of supernovae contributions, we characterise the abundance-age profiles for [Ba/Fe] and [Eu/Fe]. Our results disclose that these age-abundance relations vary across the [Fe/H]-[$\alpha$/Fe] plane. Within cells, we find negative age-[Ba/Fe] relations and flat age-[Eu/Fe] relations. Across cells, we see the slope of the age-[Ba/Fe] relations evolve smoothly and the [Eu/Fe] relations vary in amplitude. We subsequently model our empirical findings in a theoretical setting using the flexible Chempy Galactic chemical evolution (GCE) code, using the mean [Fe/H], [Mg/Fe], [Ba/Fe], and age values for stellar populations binned in [Fe/H], [Mg/Fe], and age space. We find that within a one-zone framework, an ensemble of GCE model parameters vary to explain the data. Using present day orbits from \textit{Gaia} EDR3 measurements we infer that the GCE model parameters, which set the observed chemical abundance distributions, vary systematically across mean orbital radii. Under our modelling assumptions, the observed chemical abundances are consistent with a small gradient in the high mass end of the initial mass function (IMF) across the disc, where the IMF is more top heavy towards the inner disc and more bottom heavy in the outer disc.


title: An efficient method for fitting radiation-mediated shocks to gamma-ray burst data: The Kompaneets RMS approximation

abstract: Shocks that occur below a gamma-ray burst (GRB) jet photosphere are mediated by radiation. Such radiation-mediated shocks (RMSs) could be responsible for shaping the prompt GRB emission. Although well studied theoretically, RMS models have not yet been fitted to data due to the computational cost of simulating RMSs from first principles. Here, we bridge the gap between theory and observations by developing an approximate method capable of accurately reproducing radiation spectra from mildly relativistic (in the shock frame) or slower RMSs, called the Kompaneets RMS approximation (KRA). The approximation is based on the similarities between thermal Comptonization of radiation and the bulk Comptonization that occurs inside an RMS. We validate the method by comparing simulated KRA radiation spectra to first-principle radiation-hydrodynamics simulations, finding excellent agreement both inside the RMS and in the RMS downstream. The KRA is then applied to a shock scenario inside a GRB jet, allowing for fast and efficient fitting to GRB data. We illustrate the capabilities of the developed method by performing a fit to a non-thermal spectrum in GRB 150314A. The fit allows us to uncover the physical properties of the RMS responsible for the prompt emission, such as the shock speed and the upstream plasma temperature.


title: The galaxy power spectrum on the lightcone: deep, wide-angle redshift surveys and the turnover scale

abstract: We derive expressions for the survey-window convolved galaxy power spectrum in real space for a full sky and deep redshift survey, but taking into account the geometrical lightcone effect. We investigate the impact of using the standard mean redshift approximation as a function of survey depth, and show that this assumption can lead to both an overall amplitude suppression and scale-dependent error when compared to the `true' spectrum. However, we also show that by using a carefully chosen `effective fixed-time', one can find a range of scales where the approximation to the full model is highly accurate, but only on a more restricted set of scales. We validate the theory by constructing dark matter and galaxy lightcone mock surveys from a large $N$-body simulation with a high cadence of snapshots. We do this by solving the light cone equation exactly for every particle, where the particle worldlines are obtained in a piecewise fashion with cubic interpolation between neighbouring snapshots. We find excellent agreement between our measurements and the theory ($\sim \pm 5\%$) over scales $(0.004 h{\rm Mpc}^{-1} \leq k \leq 0.54 h {\rm Mpc}^{-1})$ and for a variety of magnitude limits. Finally, we look to see how accurately we can measure the turnover scale of the galaxy power spectrum $k_0$. Using the lightcone mocks we show that one can detect the turnover scale with a probability $P \geq 95\%$ in an all-sky catalogue limited to an apparent magnitude $m_{\rm lim}\sim 21$. We also show that the detection significance would remain high for surveys with $m_{\rm lim}\sim22$ and $20\%$ sky coverage.


title: Phase diagrams of SO($N$) Majorana-Hubbard models: Dimerization, internal symmetry breaking, and fluctuation-induced first-order transitions

abstract: We study the zero-temperature phase diagrams of Majorana-Hubbard models with SO($N$) symmetry on two-dimensional honeycomb and $\pi$-flux square lattices, using mean-field and renormalization group approaches. The models can be understood as real counterparts of the SU($N$) Hubbard-Heisenberg models, and may be realized in Abrikosov vortex phases of topological superconductors, or in fractionalized phases of strongly-frustrated spin-orbital magnets. In the weakly-interacting limit, the models feature stable and fully symmetric Majorana semimetal phases. Increasing the interaction strength beyond a finite threshold for large $N$, we find a direct transition towards dimerized phases, which can be understood as staggered valence bond solid orders, in which part of the lattice symmetry is spontaneously broken and the Majorana fermions acquire a mass gap. For small to intermediate $N$, on the other hand, phases with spontaneously broken SO($N$) symmetry, which can be understood as generalized N\'eel antiferromagnets, may be stabilized. These antiferromagnetic phases feature fully gapped fermion spectra for even $N$, but gapless Majorana modes for odd $N$. While the transitions between Majorana semimetal and dimerized phases are strongly first order, the transitions between Majorana semimetal and antiferromagnetic phases are continuous for small $N \leq 3$ and weakly first order for intermediate $N \geq 4$. The weakly-first-order nature of the latter transitions arises from fixed-point annihilation in the corresponding effective field theory, which contains a real symmetric tensorial order parameter coupled to the gapless Majorana degrees of freedom, realizing interesting examples of fluctuation-induced first-order transitions.


title: Spatial regionalization based on optimal information compression

abstract: Regionalization, spatially contiguous clustering, provides a means to reduce the effect of noise in sampled data and identify homogeneous areas for policy development among many other applications. Existing regionalization methods require user input such as the number of regions or a similarity measure between regions, which does not allow for the extraction of the natural regions defined solely by the data itself. Here we view the problem of regionalization as one of data compression and develop an efficient, parameter-free regionalization algorithm based on the minimum description length principle. We demonstrate that our method is capable of recovering planted spatial clusters in noisy synthetic data, and that it can meaningfully coarse-grain real demographic data. Using our description length formulation, we find that spatial ethnoracial data in U.S. metropolitan areas has become less compressible over the period from 1980 to 2010, reflecting the rising complexity of urban segregation patterns in these metros.


title: On the Origin of Pulsar and Magnetar Magnetic Fields

abstract: In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue that these expectations should apply to the convective proto-neutron stars at the centers of core-collapse supernovae. We analyze a suite of three-dimensional simulations of core collapse, featuring a realistic equation of state and full neutrino transport, in this context. All our progenitor models, ranging from 9 solar masses to 25 solar masses, including one with initial rotation, have sufficiently vigorous proto-neutron-star convection to generate dipole fields of order ~10^15 gauss, if the modified Rossby number resides in the critical range. Thus, the magnetar/radio pulsar dichotomy may arise naturally in part from the distribution of core rotation rates in massive stars.


title: A New Constraint on the Nuclear Equation of State from Statistical Distributions of Compact Remnants of Supernovae

abstract: Understanding how matter behaves at the highest densities and temperatures is a major open problem in both nuclear physics and relativistic astrophysics. This physics is often encapsulated in the so-called high-temperature nuclear equation of state, which influences compact binary mergers, core-collapse supernovae, and many more phenomena. One such case is the type (either black hole or neutron star) and mass of the remnant of the core collapse of a massive star. For each of six candidate equations of state, we use a very large suite of spherically symmetric supernova models to generate a suite of synthetic populations of such remnants. We then compare these synthetic populations to the observed remnant population. We thus provide a novel constraint on the high-temperature nuclear equation of state and describe which EOS candidates are more or less favored by this metric.


title: Orbital Dynamics and the Evolution of Planetary Habitability in the AU Mic System

abstract: The diversity of planetary systems that have been discovered are revealing the plethora of possible architectures, providing insights into planet formation and evolution. They also increase our understanding of system parameters that may affect planetary habitability, and how such conditions are influenced by initial conditions. The AU~Mic system is unique among known planetary systems in that it is a nearby, young, multi-planet transiting system. Such a young and well characterized system provides an opportunity to study orbital dynamical and habitability studies for planets in the very early stages of their evolution. Here, we calculate the evolution of the Habitable Zone of the system through time, including the pre-main sequence phase that the system currently resides in. We discuss the planetary atmospheric processes occurring for an Earth-mass planet during this transitionary period, and provide calculations of the climate state convergence age for both volatile rich and poor initial conditions. We present results of an orbital dynamical analysis of the AU~Mic system that demonstrate the rapid eccentricity evolution of the known planets, and show that terrestrial planets within the Habitable Zone of the system can retain long-term stability. Finally, we discuss follow-up observation prospects, detectability of possible Habitable Zone planets, and how the AU Mic system may be used as a template for studies of planetary habitability evolution.


title: Learning and Predicting from Dynamic Models for COVID-19 Patient Monitoring

abstract: COVID-19 has challenged health systems to learn how to learn. This paper describes the context, methods and challenges for learning to improve COVID-19 care at one academic health center. Challenges to learning include: (1) choosing a right clinical target; (2) designing methods for accurate predictions by borrowing strength from prior patients' experiences; (3) communicating the methodology to clinicians so they understand and trust it; (4) communicating the predictions to the patient at the moment of clinical decision; and (5) continuously evaluating and revising the methods so they adapt to changing patients and clinical demands. To illustrate these challenges, this paper contrasts two statistical modeling approaches - prospective longitudinal models in common use and retrospective analogues complementary in the COVID-19 context - for predicting future biomarker trajectories and major clinical events. The methods are applied to and validated on a cohort of 1,678 patients who were hospitalized with COVID-19 during the early months of the pandemic. We emphasize graphical tools to promote physician learning and inform clinical decision making.


title: Classification of the factors of functioning of the individual human thermoregulation system and prediction of its state (local thermophysical aspects)

abstract: A classification of the tasks of human thermoregulation is proposed, which allows combining various scientific studies into one system, which takes into account environmental factors, biological feedback in the body and different properties of the constituent parts of the human thermoregulation system. The proposed multifactorial model of the human thermoregulation system is based on combining the classification of tasks represented by a group of main factors (temperature, time and the characteristic size of the system under consideration), with the classification of tasks in the studied sciences. The article discusses a number of the most significant scientific works for each area of research. Identified and systematized the most important problems for study in each of them from the point of view of heat and mass transfer and energy conversion at different levels of the description of biosystems in the human body.


title: Hunting for open clusters in Gaia EDR3: $628$ new open clusters found with OCfinder

abstract: The improvements in the precision of the published data in \textit{Gaia} EDR3 with respect to \textit{Gaia} DR2, particularly for parallaxes and proper motions, offer the opportunity to increase the number of known open clusters in the Milky Way by detecting farther and fainter objects that have so far go unnoticed. Our aim is to keep completing the open cluster census in the Milky Way with the detection of new stellar groups in the Galactic disc. We use \textit{Gaia} EDR3 up to magnitude $G = 18$ mag, increasing in one unit the magnitude limit and therefore the search volume explored in our previous studies. We use the \texttt{OCfinder} method to search for new open clusters in \textit{Gaia} EDR3 using a Big Data environment. As a first step, \texttt{OCfinder} identifies stellar statistical overdensities in the five dimensional astrometric space (position, parallax and proper motions) using the \texttt{DBSCAN} clustering algorithm. Then, these overdensities are classified into random statistical overdensities or real physical open clusters using a deep artificial neural network trained on well-characterised $G$, $G_{\rm BP} - G_{\rm RP}$ colour-magnitude diagrams. We report the discovery of $664$ new open clusters within the Galactic disc, most of them located beyond $1$ kpc from the Sun. From the estimation of ages, distances and line-of-sight extinctions of these open clusters, we see that young clusters align following the Galactic spiral arms while older ones are dispersed in the Galactic disc. Furthermore, we find that most open clusters are located at low Galactic altitudes with the exception of a few groups older than $1$ Gyr. We show the success of the \texttt{OCfinder} method leading to the discovery of a total of $1\,310$ open clusters (joining the discoveries here with the previous ones based on \textit{Gaia} DR2)[abridged]