title: Quotient symmetry protected topological phenomena

abstract: Topological phenomena are commonly studied in phases of matter which are separated from a trivial phase by an unavoidable quantum phase transition. This can be overly restrictive, leaving out scenarios of practical relevance -- similar to the distinction between liquid water and vapor. Indeed, we show that topological phenomena can be stable over a large part of parameter space even when the bulk is strictly speaking in a trivial phase of matter. In particular, we focus on symmetry-protected topological phases which can be trivialized by extending the symmetry group. The topological Haldane phase in spin chains serves as a paradigmatic example where the $SO(3)$ symmetry is extended to $SU(2)$ by tuning away from the Mott limit. Although the Haldane phase is then adiabatically connected to a product state, we show that characteristic phenomena -- edge modes, entanglement degeneracies and bulk phase transitions -- remain parametrically stable. This stability is due to a separation of energy scales, characterized by quantized invariants which are well-defined when a subgroup of the symmetry only acts on high-energy degrees of freedom. The low-energy symmetry group is a quotient group whose emergent anomalies stabilize edge modes and unnecessary criticality, which can occur in any dimension.


title: Gravitational wave cosmology I: high frequency approximation

abstract: In this paper, we systematically study gravitational waves (GWs) produced by remote compact astrophysical sources. To describe such GWs properly, we introduce three scales, $\lambda, \; L_c$ and $L$, denoting, respectively, the typical wavelength of GWs, the scale of the cosmological perturbations, and the size of the observable universe. For GWs to be detected by the current and foreseeable detectors, the condition $\lambda \ll L_c \ll L$ holds, and such GWs can be well approximated as high-frequency GWs. In order for the backreaction of the GWs to the background to be negligible, we must assume that $\left|h_{\mu\nu}\right| \ll 1$, in addition to the condition $\epsilon \ll 1$, which are also the conditions for the linearized Einstein field equations for $h_{\mu\nu}$ to be valid, where $g_{\mu\nu} = \gamma_{\mu\nu} + \epsilon h_{\mu\nu}$, and $\gamma_{\mu\nu}$ denotes the background. To simplify the field equations, we show that the spatial, traceless, and Lorentz gauge conditions can be imposed simultaneously, even when the background is not vacuum, as long as the high-frequency GW approximation is valid. However, to develop the formulas that can be applicable to as many cases as possible, we first write down explicitly the linearized Einstein field equations by imposing only the spatial gauge. Applying the general formulas together with the geometrical optics approximation to such GWs, we find that they still move along null geodesics and its polarization bi-vector is parallel-transported, even when both the cosmological scalar and tensor perturbations are present. In addition, we also calculate the gravitational integrated Sachs-Wolfe effects, whereby the dependences of the amplitude, phase and luminosity distance of the GWs on these two kinds of perturbations are read out explicitly.

categories: astro-ph.CO astro-ph.HE gr-qc hep-ph hep-th


title: Coherent Scattering-mediated correlations between levitated nanospheres

abstract: We explore entanglement generation between multiple optically levitated nanospheres interacting with a common optical cavity via the Coherent Scattering optomechanical interaction. We derive the many-particle Hamiltonian governing the unitary evolution of the system and show that it gives rise to quantum correlations among the various partitions of the setup, following a non-Markovian dynamics of entanglement birth, death and revivals. We also consider the effects of coupling the system to external environments and show that under reasonable experimental conditions entanglement between the mechanical modes can survive even at room temperature. Its dependence upon the number of nanoparticles, their initial temperature and coupling strength is studied. A numerical toolbox to simulate the closed and open dynamics of Gaussian optomechanical states and their informational measures is developed.


title: Thermal buckling transition of crystalline membranes in a field

abstract: Two dimensional crystalline membranes in isotropic embedding space exhibit a flat phase with anomalous elasticity, relevant e.g., for graphene. Here we study their thermal fluctuations in the absence of exact rotational invariance in the embedding space. An example is provided by a membrane in an orientational field, tuned to a critical buckling point by application of in-plane stresses. Through a detailed analysis, we show that the transition is in a new universality class. The self-consistent screening method predicts a second order transition, with modified anomalous elasticity exponents at criticality, while the RG suggests a weakly first order transition.


title: Collider constraints on axion-like particles

abstract: The current status and future prospects of searches for axion-like particles (ALPs) at colliders, mostly focused on the CERN LHC, are summarized. Constraints on ALPs with masses above a few GeV that couple to photons, as well as to Z or Higgs bosons, have been set at the LHC through searches for new $a\to\gamma\gamma$ resonances in di-, tri-, and four-photon final states. Inclusive and exclusive diphotons in proton-proton and lead-lead collisions, pp, PbPb $\to a \to \gamma\gamma (+X)$, as well as exotic Z and Higgs boson decays, pp $\to \mathrm{Z},\mathrm{H}\to a\gamma \to 3\gamma$ and pp $\to \mathrm{H}\to aa \to 4\gamma$, have been analyzed. Exclusive searches in PbPb collisions provide the best exclusion limits for ALP masses $m_a\approx 5-$100 GeV, whereas the other channels are the most competitive ones over $m_a\approx 100$ GeV$-$2.6 TeV. Integrated ALP production cross sections up to $\sim$100 nb are excluded at 95% confidence level, corresponding to constraints on axion-photon couplings down to $g_{a\gamma}\approx$ 0.05 TeV$^{-1}$, over broad mass ranges. Factors of 10$-$100 improvements in these limits are expected at the LHC approaching $g_{a\gamma}\approx 10^{-3}$ TeV$^{-1}$ over $m_a\approx 1$ GeV$-$5 TeV in the next decade.


title: Good and Proper: Self-similarity of N-body Simulations with Proper Force Softening

abstract: Analysis of self-similarity in scale-free $N$-body simulations reveals the spatial and temporal scales for which statistics measured in cosmological simulations are converged to the physical continuum limit. We examine how the range of scales in which the two-point correlation function is converged depends on the force softening length and whether it is held constant in comoving or proper coordinates. We find that a proper softening that reaches roughly 1/30th of the inter-particle spacing by the end of the simulation resolves the same spatial and temporal scales as a comoving softening of the same length while using a third fewer time steps, for a range of scale factors typical to $\Lambda$CDM simulations. We additionally infer an inherent resolution limit, set by the particle mass and scaling as $a^{-1/2}$, beyond which reducing the softening does not improve the resolution. We postulate a mapping of these results with spectral index $n=-2$ to $\Lambda$CDM simulations.


title: Statistical strong lensing. I. Constraints on the inner structure of galaxies from samples of a thousand lenses

abstract: Context: The number of known strong gravitational lenses is expected to grow substantially in the next few years. The statistical combination of large samples of lenses has the potential of providing strong constraints on the inner structure of galaxies. Aims: We investigate to what extent we can calibrate stellar mass measurements and constrain the average dark matter density profile of galaxies by statistically combining strong lensing data from thousands of lenses. Methods: We generate mock samples of axisymmetric lenses. We assume that, for each lens, we have measurements of two image positions of a strongly lensed background source, as well as magnification information from full surface brightness modelling, and a stellar population synthesis-based estimate of the lens stellar mass. We then fit models describing the distribution of the stellar population synthesis mismatch parameter $\alpha_{sps}$ (the ratio between the true stellar mass and the stellar population synthesis-based estimate) and dark matter density profile of the population of lenses to an ensemble of 1000 mock lenses. Results: The average $\alpha_{sps}$, projected dark matter mass and dark matter density slope can be obtained with great precision and accuracy, compared with current constraints. A flexible model and the knowledge of the lens detection efficiency as a function of image configuration are required in order to avoid a biased inference. Conclusions: Statistical strong lensing inferences from upcoming surveys have the potential to calibrate stellar mass measurements and to constrain the inner dark matter density profile of massive galaxies.


title: Flips in symmetric separated set-systems

abstract: For a positive integer $n$, a collection $S$ of subsets of $[n]=\{1,\ldots,n\}$ is called symmetric if $X\in S$ implies $X^\ast\in S$, where $X^\ast:=\{i\in [n]\colon n-i+1\notin X\}$ (the involution $\ast$ was introduced by Karpman). Leclerc and Zelevinsky showed that the set of maximal strongly (resp. weakly) separated collections in $2^{[n]}$ is connected via flips, or mutations, ``in the presence of six (resp. four) witnesses''. We give a symmetric analog of those results, by showing that each maximal symmetric strongly (weakly) separated collection in $2^{[n]}$ can be obtained from any other one by a series of special symmetric local transformations, so-called symmetric flips. Also we establish the connectedness via symmetric flips for the class of maximal symmetric $r$-separated collections in $2^{[n]}$ when $n,r$ are even (where sets $A,B\subseteq [n]$ are called $r$-separated if there are no elements $i_0<i_1< \cdots <i_{r+1}$ in $[n]$ which alternate in $A\setminus B$ and $B\setminus A$). This is related to a symmetric version of higher Bruhat orders. These results are obtained as consequences of our study of related geometric objects: symmetric rhombus and combined tilings and symmetric cubillages.


title: Adaptive Doubly Robust Estimator from Non-stationary Logging Policy under a Convergence of Average Probability

abstract: Adaptive experiments, including efficient average treatment effect estimation and multi-armed bandit algorithms, have garnered attention in various applications, such as social experiments, clinical trials, and online advertisement optimization. This paper considers estimating the mean outcome of an action from samples obtained in adaptive experiments. In causal inference, the mean outcome of an action has a crucial role, and the estimation is an essential task, where the average treatment effect estimation and off-policy value estimation are its variants. In adaptive experiments, the probability of choosing an action (logging policy) is allowed to be sequentially updated based on past observations. Due to this logging policy depending on the past observations, the samples are often not independent and identically distributed (i.i.d.), making developing an asymptotically normal estimator difficult. A typical approach for this problem is to assume that the logging policy converges in a time-invariant function. However, this assumption is restrictive in various applications, such as when the logging policy fluctuates or becomes zero at some periods. To mitigate this limitation, we propose another assumption that the average logging policy converges to a time-invariant function and show the doubly robust (DR) estimator's asymptotic normality. Under the assumption, the logging policy itself can fluctuate or be zero for some actions. We also show the empirical properties by simulations.


title: EEG-based Texture Roughness Classification in Active Tactile Exploration with Invariant Representation Learning Networks

abstract: During daily activities, humans use their hands to grasp surrounding objects and perceive sensory information which are also employed for perceptual and motor goals. Multiple cortical brain regions are known to be responsible for sensory recognition, perception and motor execution during sensorimotor processing. While various research studies particularly focus on the domain of human sensorimotor control, the relation and processing between motor execution and sensory processing is not yet fully understood. Main goal of our work is to discriminate textured surfaces varying in their roughness levels during active tactile exploration using simultaneously recorded electroencephalogram (EEG) data, while minimizing the variance of distinct motor exploration movement patterns. We perform an experimental study with eight healthy participants who were instructed to use the tip of their dominant hand index finger while rubbing or tapping three different textured surfaces with varying levels of roughness. We use an adversarial invariant representation learning neural network architecture that performs EEG-based classification of different textured surfaces, while simultaneously minimizing the discriminability of motor movement conditions (i.e., rub or tap). Results show that the proposed approach can discriminate between three different textured surfaces with accuracies up to 70%, while suppressing movement related variability from learned representations.


title: Photon-photon polaritons in chi(2) microresonators

abstract: We consider a high-Q microresonator with $\chi^{(2)}$ nonlinearity under conditions when the coupling rates between the sidebands around the pump and second harmonic exceed the damping rates, implying the strong coupling regime (SC). Using the dressed-resonator approach we demonstrate that this regime leads to the dominance of the Hermitian part of the operator driving the side-band dynamics over its non-Hermitian part responsible for the parametric gain. This has allowed us to introduce and apply the cross-area concept of the polariton quasi-particles and define their effective masses in the context of $\chi^{(2)}$ ring-microresonators. We further use polaritons to predict the modified spectral response of the resonator to a weak probe field, and to reveal splitting of the bare-resonator resonances, avoided crossings, and Rabi dynamics. Polariton basis also allows deriving a discrete sequence of the parametric thresholds for the generation of sidebands of different orders.


title: K2-138 g: Spitzer Spots a Sixth Planet for the Citizen Science System

abstract: $K2$ greatly extended $Kepler$'s ability to find new planets, but it was typically limited to identifying transiting planets with orbital periods below 40 days. While analyzing $K2$ data through the Exoplanet Explorers project, citizen scientists helped discover one super-Earth and four sub-Neptune sized planets in the relatively bright ($V=12.21$, $K=10.3$) K2-138 system, all which orbit near 3:2 mean motion resonances. The $K2$ light curve showed two additional transit events consistent with a sixth planet. Using $Spitzer$ photometry, we validate the sixth planet's orbital period of $41.966\pm0.006$ days and measure a radius of $3.44^{+0.32}_{-0.31}\,R_{\oplus}$, solidifying K2-138 as the $K2$ system with the most currently known planets. There is a sizeable gap between the outer two planets, since the fifth planet in the system, K2-138 f, orbits at 12.76 days. We explore the possibility of additional non-transiting planets in the gap between f and g. Due to the relative brightness of the K2-138 host star, and the near resonance of the inner planets, K2-138 could be a key benchmark system for both radial velocity and transit timing variation mass measurements, and indeed radial velocity masses for the inner four planets have already been obtained. With its five sub-Neptunes and one super-Earth, the K2-138 system provides a unique test bed for comparative atmospheric studies of warm to temperate planets of similar size, dynamical studies of near resonant planets, and models of planet formation and migration.


title: Tuning the flat bands of the kagome metal CoSn with Fe, In, or Ni doping

abstract: CoSn is a Pauli paramagnet with relatively flat d-bands centered about 100 meV below the Fermi energy Ef. Single crystals of CoSn lightly doped with Fe, In, or Ni are investigated using x-ray and neutron scattering, magnetic susceptibility and magnetization, ac susceptibility, specific heat and resistivity measurements. Within the rigid band approximation, hole doping with a few percent of Fe or In should move the flat bands closer to Ef, whereas electron doping with Ni should move the flat bands further away from Ef. We provide evidence that this indeed occurs. Fe and In doping drive CoSn toward magnetism, while Ni doping suppresses CoSn's already weak magnetic response. The resulting ground state is different for Fe versus In doping. For Fe-doped crystals, Co1-xFexSn, with 0.02 < x < 0.27, the magnetic and specific heat data are consistent with the formation of a spin glass, with a glass transition temperature, Tg, ranging from 1 K for x=0.02 to 10 K for x= 0.27. Powder and single crystal neutron diffraction found no evidence of long-range magnetic order below Tg with x = 0.17. For In-doped crystals, CoSn1-yIny, both the magnetic susceptibility and the Sommerfeld coefficient, gamma, increase substantially relative to pure CoSn, but with no clear indication of a magnetic transition for 0.05 < y < 0.2. CoSn crystals doped with Ni (Co0.93Ni0.07Sn) have a significantly smaller magnetic susceptibility and gamma than pure CoSn, consistent with the flat bands further from Ef.


title: The Dependence of the Type Ia Supernova Host Bias on Observation or Fitting Technique

abstract: More luminous Type Ia supernovae (SNe Ia) prefer less massive hosts and regions of higher star formation. This correlation is inverted during width-color-luminosity light curve standardization resulting in step-like biases of distance measurements with respect to host properties. Using the PISCO supernova host sample and SDSS, GALEX, and 2MASS photometry, we compare host stellar mass and specific star formation rate (sSFR) from different observation methods, including local vs. global, and fitting techniques to measure their impact on the host step biases. Mass step measurements for all our mass samples are consistent within a 1$\sigma$ significance from -0.03$\pm$0.02 mag to -0.04$\pm$0.02 mag. Including or excluding UV information had no effect on measured mass step size or location. Specific SFR (sSFR) step sizes are more significant than mass step measurements and varied from $0.05\pm0.03$ mag (H$\alpha$) and $0.06\pm0.02$ mag (UV) for a 51 host sample. The sSFR step location is influenced by mass sample used to normalize star formation and by sSFR tracer choice. The step size is reduced to 0.04$\pm$0.03 mag when using all available 73 hosts with H$\alpha$ measurements. This 73 PISCO host subsample overall lacked a clear step signal, but here we are searching for whether different choices of mass or sSFR estimation can create a step signal. We find no evidence that different observation or fitting techniques choice can create a distance measurement step in either mass or sSFR.


title: Conceptual 12M: Pushing Web-Scale Image-Text Pre-Training To Recognize Long-Tail Visual Concepts

abstract: The availability of large-scale image captioning and visual question answering datasets has contributed significantly to recent successes in vision-and-language pre-training. However, these datasets are often collected with overrestrictive requirements inherited from their original target tasks (e.g., image caption generation), which limit the resulting dataset scale and diversity. We take a step further in pushing the limits of vision-and-language pre-training data by relaxing the data collection pipeline used in Conceptual Captions 3M (CC3M) [Sharma et al. 2018] and introduce the Conceptual 12M (CC12M), a dataset with 12 million image-text pairs specifically meant to be used for vision-and-language pre-training. We perform an analysis of this dataset and benchmark its effectiveness against CC3M on multiple downstream tasks with an emphasis on long-tail visual recognition. Our results clearly illustrate the benefit of scaling up pre-training data for vision-and-language tasks, as indicated by the new state-of-the-art results on both the nocaps and Conceptual Captions benchmarks.


title: Bulk geometry in gauge/gravity duality and color degrees of freedom

abstract: U($N$) supersymmetric Yang-Mills theory naturally appears as the low-energy effective theory of a system of $N$ D-branes and open strings between them. Transverse spatial directions emerge from scalar fields, which are $N\times N$ matrices with color indices; roughly speaking, the eigenvalues are the locations of D-branes. In the past, it was argued that this simple 'emergent space' picture cannot be used in the context of gauge/gravity duality, because the ground-state wave function delocalizes at large $N$, leading to a conflict with the locality in the bulk geometry. In this paper we show that this conventional wisdom is not correct: the ground-state wave function does not delocalize, and there is no conflict with the locality of the bulk geometry. This conclusion is obtained by clarifying the meaning of the 'diagonalization of a matrix' in Yang-Mills theory, which is not as obvious as one might think. This observation opens up the prospect of characterizing the bulk geometry via the color degrees of freedom in Yang-Mills theory, all the way down to the center of the bulk.


title: Automated Detection of Equine Facial Action Units

abstract: The recently developed Equine Facial Action Coding System (EquiFACS) provides a precise and exhaustive, but laborious, manual labelling method of facial action units of the horse. To automate parts of this process, we propose a Deep Learning-based method to detect EquiFACS units automatically from images. We use a cascade framework; we firstly train several object detectors to detect the predefined Region-of-Interest (ROI), and secondly apply binary classifiers for each action unit in related regions. We experiment with both regular CNNs and a more tailored model transferred from human facial action unit recognition. Promising initial results are presented for nine action units in the eye and lower face regions. Code for the project is publicly available.


title: The *-Edge-Reinforced Random Walk

abstract: We define a linearly reinforced process called the *-Edge-Reinforced Random Walk (*-ERRW ) which can be seen as a Yaglom reversible, hence non-reversible, extension of the Edge-Reinforced Random Walk (ERRW) introduced by Coppersmith and Diaconis in 1986. This family of processes also generalizes the r-dependent ERRW introduced by Bacallado (2009). Under some assumptions on the initial weights, the *-ERRW is partially exchangeable in the sense of Diaconis and Freedman (1980), and thus it is a random walk in a random environment. The main result of the paper gives the explicit expression of the mixing law, hence extending the "magic formula" of Coppersmith and Diaconis from the case of mixtures of reversible Markov chains to the case of mixtures of Yaglom reversible Markov chains.


title: Scanning the Cycle: Timing-based Authentication on PLCs

abstract: Programmable Logic Controllers (PLCs) are a core component of an Industrial Control System (ICS). However, if a PLC is compromised or the commands sent across a network from the PLCs are spoofed, consequences could be catastrophic. In this work, a novel technique to authenticate PLCs is proposed that aims at raising the bar against powerful attackers while being compatible with real-time systems. The proposed technique captures timing information for each controller in a non-invasive manner. It is argued that Scan Cycle is a unique feature of a PLC that can be approximated passively by observing network traffic. An attacker that spoofs commands issued by the PLCs would deviate from such fingerprints. To detect replay attacks a PLC Watermarking technique is proposed. PLC Watermarking models the relationship between the scan cycle and the control logic by modeling the input/output as a function of request/response messages of a PLC. The proposed technique is validated on an operational water treatment plant (SWaT) and smart grid (EPIC) testbed. Results from experiments indicate that PLCs can be distinguished based on their scan cycle timing characteristics.


title: Dark matter candidates in the NMSSM with RH neutrino superfields

abstract: R-parity conserving supersymmetric models with right-handed (RH) neutrinos are very appealing since they could naturally explain neutrino physics and also provide a good dark matter (DM) candidate such as the lightest supersymmetric particle (LSP). In this work we consider the next-to-minimal supersymmetric standard model (NMSSM) plus RH neutrino superfields, with effective Majorana masses dynamically generated at the electroweak scale (EW). We perform a scan of the relevant parameter space and study both possible DM candidates: RH sneutrino and neutralino. Especially for the case of RH sneutrino DM we analyse the intimate relation between both candidates to obtain the correct amount of relic density. Besides the well-known resonances, annihilations through scalar quartic couplings and coannihilation mechanisms with all kind of neutralinos, are crucial. Finally, we present the impact of current and future direct and indirect detection experiments on both DM candidates.


title: RFI Mitigation for One-bit UWB Radar Systems

abstract: Radio frequency interference (RFI) mitigation is critical to the proper operation of ultra-wideband (UWB) radar systems since RFI can severely degrade the radar imaging capability and target detection performance. In this paper, we address the RFI mitigation problem for one-bit UWB radar systems. A one-bit UWB system obtains its signed measurements via a low-cost and high rate sampling scheme, referred to as the Continuous Time Binary Value (CTBV) technology. This sampling strategy compares the signal to a known threshold varying with slow-time and therefore can be used to achieve a rather high sampling rate and quantization resolution with rather simple and affordable hardware. This paper establishes a proper data model for the RFI sources and proposes a novel RFI mitigation method for the one-bit UWB radar system that uses the CTBV sampling technique. Specifically, we first model the RFI sources as a sum of sinusoids with frequencies fixed during the coherent processing interval (CPI) and we exploit the sparsity of the RFI spectrum. We extend a majorization-minimization based 1bRELAX algorithm, referred to as 1bMMRELAX, to estimate the RFI source parameters from the signed measurements obtained by using the CTBV sampling strategy. We also devise a new fast frequency initialization method based on the Alternating Direction Method of Multipliers (ADMM) methodology for the extended 1bMMRELAX algorithm to significantly improve its computational efficiency. Moreover, an ADMM-based sparse method is introduced to recover the desired radar echoes using the estimated RFI parameters. Both simulated and experimental results are presented to demonstrate that our proposed algorithm outperforms the existing digital integration method, especially for severe RFI cases.


title: The *-Vertex-Reinforced Jump Process I : exchangeability and mixing measure

abstract: We introduce a non-reversible generalization of the Vertex-Reinforced Jump Process (VRJP), which we call *-Vertex-Reinforced Jump Process (*-VRJP). It can be seen as the continuous-time counterpart of the *-Edge-Reinforced Random Walk (*-ERRW) (see Bacallado (2009), Bacallado, Sabot and Tarr\`es (2020)), which is itself a non-reversible, and in fact Yaglom reversible, generalization of the original ERRW introduced by Coppersmith and Diaconis in 1986. In contrast to the classical VRJP, the *-VRJP is not exchangeable after time change, which leads to several difficulties and new phenomena. However, we show that with some appropriate randomization of the initial local time, it becomes partially exchangeable after time change. We provide a representation of the "randomized" *-VRJP as mixture of Yaglom reversible Markov jump processes with an explicit mixing measure, and we prove that the non-randomized *-VRJP can be written as a mixture of conditioned Markov processes.


title: Tunable quantum confinement of neutral excitons using electric fields and exciton-charge interactions

abstract: Quantum confinement is the discretization of energy when motion of particles is restricted to length scales smaller than their de Broglie wavelength. The experimental realization of this effect has had wide ranging impact in diverse fields of physics and facilitated the development of new technologies. In semiconductor physics, quantum confinement of optically excited quasiparticles, such as excitons or trions, is typically achieved by modulation of material properties - an approach crucially limited by the lack of insitu tunability and scalability of confining potentials. Achieving fully tunable quantum confinement of optical excitations has therefore been an outstanding goal in quantum photonics. Here, we demonstrate electrically controlled quantum confinement of neutral excitons in a gate-defined monolayer p-i-n diode. A combination of dc Stark shift induced by large in-plane fields and a previously unknown confining mechanism based on repulsive interaction between excitons and free charges ensures tight exciton confinement in the narrow neutral region. Quantization of exciton motion manifests in multiple discrete, spectrally narrow, voltage-dependent optical resonances that emerge below the free exciton resonance. Our measurements reveal several unique physical features of these quantum confined excitons, including an in-plane dipolar character, one-dimensional center-of-mass confinement, and strikingly enhanced exciton size in the presence of magnetic fields. Our method provides an experimental route towards creating scalable arrays of identical single photon sources, which will constitute building blocks of strongly correlated photonic systems.


title: BEDS: Bagging ensemble deep segmentation for nucleus segmentation with testing stage stain augmentation

abstract: Reducing outcome variance is an essential task in deep learning based medical image analysis. Bootstrap aggregating, also known as bagging, is a canonical ensemble algorithm for aggregating weak learners to become a strong learner. Random forest is one of the most powerful machine learning algorithms before deep learning era, whose superior performance is driven by fitting bagged decision trees (weak learners). Inspired by the random forest technique, we propose a simple bagging ensemble deep segmentation (BEDs) method to train multiple U-Nets with partial training data to segment dense nuclei on pathological images. The contributions of this study are three-fold: (1) developing a self-ensemble learning framework for nucleus segmentation; (2) aggregating testing stage augmentation with self-ensemble learning; and (3) elucidating the idea that self-ensemble and testing stage stain augmentation are complementary strategies for a superior segmentation performance. Implementation Detail: https://github.com/xingli1102/BEDs.


title: Generalization in Quantum Machine Learning: a Quantum Information Perspective

abstract: Quantum classification and hypothesis testing are two tightly related subjects, the main difference being that the former is data driven: how to assign to quantum states $\rho(x)$ the corresponding class $c$ (or hypothesis) is learnt from examples during training, where $x$ can be either tunable experimental parameters or classical data "embedded" into quantum states. Does the model generalize? This is the main question in any data-driven strategy, namely the ability to predict the correct class even of previously unseen states. Here we establish a link between quantum machine learning classification and quantum hypothesis testing (state and channel discrimination) and then show that the accuracy and generalization capability of quantum classifiers depend on the (R\'enyi) mutual informations $I(C{:}Q)$ and $I_2(X{:}Q)$ between the quantum state space $Q$ and the classical parameter space $X$ or class space $C$. Based on the above characterization, we then show how different properties of $Q$ affect classification accuracy and generalization, such as the dimension of the Hilbert space, the amount of noise, and the amount of neglected information from $X$ via, e.g., pooling layers. Moreover, we introduce a quantum version of the Information Bottleneck principle that allows us to explore the various tradeoffs between accuracy and generalization. Finally, in order to check our theoretical predictions, we study the classification of the quantum phases of an Ising spin chain, and we propose the Variational Quantum Information Bottleneck (VQIB) method to optimize quantum embeddings of classical data to favor generalization.


title: SIS Epidemic Model: Birth-and-Death Markov Chain Approach

abstract: We are interested in describing the infected size of the SIS Epidemic model using Birth-Death Markov process. The Susceptible-Infected-Susceptible (SIS) model is defined within a population of constant size $M$; the size is kept constant by replacing each death with a newborn healthy individual. The life span of each individual in the population is modelled by an exponential distribution with parameter $\alpha$; and the disease spreads within the population is modelled by a Poisson process with a rate $\lambda_{I}$. $\lambda_{I}=\beta I(1-\frac{I}{M}) $ is similar to the instantaneous rate in the logistic population growth model. The analysis is focused on the disease outbreak, where the reproduction number $R=\frac{\beta}{\alpha} $ is greater than one. As methodology, we use both numerical and analytical approaches. The analysis relies on the stationary distribution for Birth and Death Markov process. The numerical approach creates sample path simulations into order show the infected size dynamics, and the relationship between infected size and $R$. As $M$ becomes large, some stable statistical characteristics of the infected size distribution can be deduced. And the infected size is shown analytically to follow a normal distribution with mean $(1-\frac{1}{R}) M$ and Variance $\frac{M}{R} $.


title: Using Distance Correlation for Efficient Bayesian Optimization

abstract: We propose a novel approach for Bayesian optimization, called $\textsf{GP-DC}$, which combines Gaussian processes with distance correlation. It balances exploration and exploitation automatically, and requires no manual parameter tuning. We evaluate $\textsf{GP-DC}$ on a number of benchmark functions and observe that it outperforms state-of-the-art methods such as $\textsf{GP-UCB}$ and max-value entropy search, as well as the classical expected improvement heuristic. We also apply $\textsf{GP-DC}$ to optimize sequential integral observations with a variable integration range and verify its empirical efficiency on both synthetic and real-world datasets.


title: Big Data meets Causal Survey Research: Understanding Nonresponse in the Recruitment of a Mixed-mode Online Panel

abstract: Survey scientists increasingly face the problem of high-dimensionality in their research as digitization makes it much easier to construct high-dimensional (or "big") data sets through tools such as online surveys and mobile applications. Machine learning methods are able to handle such data, and they have been successfully applied to solve \emph{predictive} problems. However, in many situations, survey statisticians want to learn about \emph{causal} relationships to draw conclusions and be able to transfer the findings of one survey to another. Standard machine learning methods provide biased estimates of such relationships. We introduce into survey statistics the double machine learning approach, which gives approximately unbiased estimators of causal parameters, and show how it can be used to analyze survey nonresponse in a high-dimensional panel setting.


title: Integer colorings with no rainbow 3-term arithmetic progression

abstract: In this paper, we study the rainbow Erd\H{o}s-Rothschild problem with respect to 3-term arithmetic progressions. We obtain the asymptotic number of $r$-colorings of $[n]$ without rainbow 3-term arithmetic progressions, and we show that the typical colorings with this property are 2-colorings. We also prove that $[n]$ attains the maximum number of rainbow 3-term arithmetic progression-free $r$-colorings among all subsets of $[n]$. Moreover, the exact number of rainbow 3-term arithmetic progression-free $r$-colorings of $\mathbb{Z}_p$ is obtained, where $p$ is any prime and $\mathbb{Z}_p$ is the cyclic group of order $p$.


title: Higher algebra of $A_\infty$ and $\Omega B As$-algebras in Morse theory II

abstract: This paper introduces the notion of $n$-morphisms between two $A_\infty$-algebras, such that 0-morphisms correspond to standard $A_\infty$-morphisms and 1-morphisms correspond to $A_\infty$-homotopies between $A_\infty$-morphisms. The set of higher morphisms between two $A_\infty$-algebras then defines a simplicial set which has the property of being an algebraic $\infty$-category. The operadic structure of $n-A_\infty$-morphisms is also encoded by new families of polytopes, which we call the $n$-multiplihedra and which generalize the standard multiplihedra. These are constructed from the standard simplices and multiplihedra by lifting the Alexander-Whitney map to the level of simplices. Rich combinatorics arise in this context, as conveniently described in terms of overlapping partitions. Shifting from the $A_\infty$ to the $\Omega B As$ framework, we define the analogous notion of $n$-morphisms between $\Omega B As$-algebras, which are again encoded by the $n$-multiplihedra, endowed with a refined cell decomposition by stable gauged ribbon tree type. We then realize this higher algebra of $A_\infty$ and $\Omega B As$-algebras in Morse theory. Given two Morse functions $f$ and $g$, we construct $n-\Omega B As$-morphisms between their respective Morse cochain complexes endowed with their $\Omega B As$-algebra structures, by counting perturbed Morse gradient trees associated to an admissible simplex of perturbation data. We moreover show that the simplicial set consisting of higher morphisms defined by a count of perturbed Morse gradient trees is a contractible Kan complex.


title: One-shot action recognition in challenging therapy scenarios

abstract: One-shot action recognition aims to recognize new action categories from a single reference example, typically referred to as the anchor example. This work presents a novel approach for one-shot action recognition in the wild that computes motion representations robust to variable kinematic conditions. One-shot action recognition is then performed by evaluating anchor and target motion representations. We also develop a set of complementary steps that boost the action recognition performance in the most challenging scenarios. Our approach is evaluated on the public NTU-120 one-shot action recognition benchmark, outperforming previous action recognition models. Besides, we evaluate our framework on a real use-case of therapy with autistic people. These recordings are particularly challenging due to high-level artifacts from the patient motion. Our results provide not only quantitative but also online qualitative measures, essential for the patient evaluation and monitoring during the actual therapy.


title: Fractional charge and fractional statistics in the quantum Hall effects

abstract: Quasiparticles with fractional charge and fractional statistics are key features of the fractional quantum Hall effect. We discuss in detail the definitions of fractional charge and statistics and the ways in which these properties may be observed. In addition to theoretical foundations, we review the present status of the experiments in the area. We also discuss the notions of non-Abelian statistics and attempts to find experimental evidence for the existence of non-Abelian quasiparticles in certain quantum Hall systems.


title: Local Monodromy of 1-Dimensional p-Divisible Groups

abstract: Let $G$ be a $p$-divisible group over a complete discrete valuation ring $R$ of characteristic $p$. The generic fiber of $G$ determines a Galois representation $\rho$. The image of $\rho$ admits a ramification filtration and a Lie filtration. We relate these filtrations in the case $G$ is one dimensional, giving an equicharacteristic version of Sen's theorem in this setting. This result generalizes results of Gross and Chai. Additionally, we prove that the representation associated to the \'etale part of $G$ is irreducible.


title: Mobile Computational Photography: A Tour

abstract: The first mobile camera phone was sold only 20 years ago, when taking pictures with one's phone was an oddity, and sharing pictures online was unheard of. Today, the smartphone is more camera than phone. How did this happen? This transformation was enabled by advances in computational photography -the science and engineering of making great images from small form factor, mobile cameras. Modern algorithmic and computing advances, including machine learning, have changed the rules of photography, bringing to it new modes of capture, post-processing, storage, and sharing. In this paper, we give a brief history of mobile computational photography and describe some of the key technological components, including burst photography, noise reduction, and super-resolution. At each step, we may draw naive parallels to the human visual system.


title: Towards AIOps in Edge Computing Environments

abstract: Edge computing was introduced as a technical enabler for the demanding requirements of new network technologies like 5G. It aims to overcome challenges related to centralized cloud computing environments by distributing computational resources to the edge of the network towards the customers. The complexity of the emerging infrastructures increases significantly, together with the ramifications of outages on critical use cases such as self-driving cars or health care. Artificial Intelligence for IT Operations (AIOps) aims to support human operators in managing complex infrastructures by using machine learning methods. This paper describes the system design of an AIOps platform which is applicable in heterogeneous, distributed environments. The overhead of a high-frequency monitoring solution on edge devices is evaluated and performance experiments regarding the applicability of three anomaly detection algorithms on edge devices are conducted. The results show, that it is feasible to collect metrics with a high frequency and simultaneously run specific anomaly detection algorithms directly on edge devices with a reasonable overhead on the resource utilization.


title: Impartial selection with prior information

abstract: We study the problem of {\em impartial selection}, a topic that lies at the intersection of computational social choice and mechanism design. The goal is to select the most popular individual among a set of community members. The input can be modeled as a directed graph, where each node represents an individual, and a directed edge indicates nomination or approval of a community member to another. An {\em impartial mechanism} is robust to potential selfish behavior of the individuals and provides appropriate incentives to voters to report their true preferences by ensuring that the chance of a node to become a winner does not depend on its outgoing edges. The goal is to design impartial mechanisms that select a node with an in-degree that is as close as possible to the highest in-degree. We measure the efficiency of such a mechanism by the difference of these in-degrees, known as its {\em additive} approximation. In particular, we study the extent to which prior information on voters' preferences could be useful in the design of efficient deterministic impartial selection mechanisms with good additive approximation guarantees. We consider three models of prior information, which we call the {\em opinion poll}, the {\em a prior popularity}, and the {\em uniform} model. We analyze the performance of a natural selection mechanism that we call {\em approval voting with default} (AVD) and show that it achieves a $O(\sqrt{n\ln{n}})$ additive guarantee for opinion poll and a $O(\ln^2n)$ for a priori popularity inputs, where $n$ is the number of individuals. We consider this polylogarithmic bound as our main technical contribution. We complement this last result by showing that our analysis is close to tight, showing an $\Omega(\ln{n})$ lower bound. This holds in the uniform model, which is the simplest among the three models.


title: Domain Impression: A Source Data Free Domain Adaptation Method

abstract: Unsupervised Domain adaptation methods solve the adaptation problem for an unlabeled target set, assuming that the source dataset is available with all labels. However, the availability of actual source samples is not always possible in practical cases. It could be due to memory constraints, privacy concerns, and challenges in sharing data. This practical scenario creates a bottleneck in the domain adaptation problem. This paper addresses this challenging scenario by proposing a domain adaptation technique that does not need any source data. Instead of the source data, we are only provided with a classifier that is trained on the source data. Our proposed approach is based on a generative framework, where the trained classifier is used for generating samples from the source classes. We learn the joint distribution of data by using the energy-based modeling of the trained classifier. At the same time, a new classifier is also adapted for the target domain. We perform various ablation analysis under different experimental setups and demonstrate that the proposed approach achieves better results than the baseline models in this extremely novel scenario.


title: L\'evy Processes, Generalized Moments and Uniform Integrability

abstract: We give new proofs of certain equivalent conditions for the existence of generalized moments of a L\'evy process $(X_t)_{t\geq 0}$; in particular, the existence of a generalized $g$-moment is equivalent to the uniform integrability of $(g(X_t))_{t\in [0,1]}$. As a consequence, certain functions of a L\'evy process which are integrable and local martingales are already true martingales. Our methods extend to moments of stochastically continuous additive processes, and we give new, short proofs for the characterization of lattice distributions and the transience of L\'evy processes.


title: An Efficient Diagnosis Algorithm for Inconsistent Constraint Sets

abstract: Constraint sets can become inconsistent in different contexts. For example, during a configuration session the set of customer requirements can become inconsistent with the configuration knowledge base. Another example is the engineering phase of a configuration knowledge base where the underlying constraints can become inconsistent with a set of test cases. In such situations we are in the need of techniques that support the identification of minimal sets of faulty constraints that have to be deleted in order to restore consistency. In this paper we introduce a divide-and-conquer based diagnosis algorithm (FastDiag) which identifies minimal sets of faulty constraints in an over-constrained problem. This algorithm is specifically applicable in scenarios where the efficient identification of leading (preferred) diagnoses is crucial. We compare the performance of FastDiag with the conflict-directed calculation of hitting sets and present an in-depth performance analysis that shows the advantages of our approach.


title: Difference of convex algorithms for bilevel programs with applications in hyperparameter selection

abstract: In this paper, we present difference of convex algorithms for solving bilevel programs in which the upper level objective functions are difference of convex functions, and the lower level programs are fully convex. This nontrivial class of bilevel programs provides a powerful modelling framework for dealing with applications arising from hyperparameter selection in machine learning. Thanks to the full convexity of the lower level program, the value function of the lower level program turns out to be convex and hence the bilevel program can be reformulated as a difference of convex bilevel program. We propose two algorithms for solving the reformulated difference of convex program and show their convergence under very mild assumptions. Finally we conduct numerical experiments to a bilevel model of support vector machine classification.


title: Nonlinear reconstruction of features in the primordial power spectrum from large-scale structure

abstract: Potential features in the primordial power spectrum have been searched for in galaxy surveys in recent years since these features can assist in understanding the nature of inflation. The null detection to date suggests that any such features should be fairly weak, and next-generation galaxy surveys, with their unprecedented sizes and precisions, are in a position to place stronger constraints than before. However, even if such primordial features once existed in the early Universe, they would have been significantly damped in the nonlinear regime at low redshift due to structure formation, which makes them difficult to be directly detected in real observations. A potential way to tackle this challenge for probing the features is to undo the cosmological evolution, i.e., using reconstruction to obtain an approximate linear density field. By employing a set of N-body simulations, here we show that a recently-proposed nonlinear reconstruction algorithm can effectively retrieve damped oscillatory features from halo catalogues and improve the accuracy of the measurement of feature parameters (assuming that such primordial features do exist). We do a Fisher analysis to forecast how nonlinear reconstruction affects the constraining power, and find that it can lead to significantly more robust constraints on the feature amplitude for a DESI-like survey. Comparing nonlinear reconstruction with other ways of improving constraints, such as increasing the survey volume and range of scales, this shows that it is possible to achieve what the latter do, but at a lower cost.


title: Divide-and-Conquer MCMC for Multivariate Binary Data

abstract: The analysis of large scale medical claims data has the potential to improve quality of care by generating insights which can be used to create tailored medical programs. In particular, the multivariate probit model can be used to investigate the correlation between multiple binary responses of interest in such data, e.g. the presence of multiple chronic conditions. Bayesian modeling is well suited to such analyses because of the automatic uncertainty quantification provided by the posterior distribution. A complicating factor is that large medical claims datasets often do not fit in memory, which renders the estimation of the posterior using traditional Markov Chain Monte Carlo (MCMC) methods computationally infeasible. To address this challenge, we extend existing divide-and-conquer MCMC algorithms to the multivariate probit model, demonstrating, via simulation, that they should be preferred over mean-field variational inference when the estimation of the latent correlation structure between binary responses is of primary interest. We apply this algorithm to a large database of de-identified Medicare Advantage claims from a single large US health insurance provider, where we find medically meaningful groupings of common chronic conditions and asses the impact of the urban-rural health gap by identifying underutilized provider specialties in rural areas.


title: BORE: Bayesian Optimization by Density-Ratio Estimation

abstract: Bayesian optimization (BO) is among the most effective and widely-used blackbox optimization methods. BO proposes solutions according to an explore-exploit trade-off criterion encoded in an acquisition function, many of which are computed from the posterior predictive of a probabilistic surrogate model. Prevalent among these is the expected improvement (EI) function. The need to ensure analytical tractability of the predictive often poses limitations that can hinder the efficiency and applicability of BO. In this paper, we cast the computation of EI as a binary classification problem, building on the link between class-probability estimation and density-ratio estimation, and the lesser-known link between density-ratios and EI. By circumventing the tractability constraints, this reformulation provides numerous advantages, not least in terms of expressiveness, versatility, and scalability.


title: Observational evidence of fractality in the large-scale distribution of galaxies

abstract: Using a sample of 133 991 galaxies distributed in the sky region $100^{\circ} <\alpha<270^{\circ}$ and $7^{\circ}<\delta<65^{\circ}$, extracted from the SDSS NASA/AMES Value Added Galaxy Catalog (AMES-VAGC), we estimate the fractal dimension using two different methods. First, using an algorithm to estimate the correlation dimension. The second method, in a novel approach, creates a graph from the data and estimates the graph dimension purely from connectivity information. In both methods we found a dimension $D\approx 2$ in scales below 20 Mpc, which agrees with previous works. This result shows the non-homogeneity of galaxies distribution at certain scales.


title: Power Minimization in Vehicular Cloud Architecture

abstract: Modern vehicles equipped with on-board units (OBU) are playing an essential role in the smart city revolution. The vehicular processing resources, however, are not used to their fullest potential. The concept of vehicular clouds is proposed to exploit the underutilized vehicular resources to supplement cloud computing services to relieve the burden on cloud data centers and improve quality of service. In this paper we introduce a vehicular cloud architecture supported by fixed edge computing nodes and the central cloud. A mixed integer linear programming (MLP) model is developed to optimize the allocation of the computing demands in the distributed architecture while minimizing power consumption. The results show power savings as high as 84% over processing in the conventional cloud. A heuristic with performance approaching that of the MILP model is developed to allocate computing demands in real time.


title: Improving Hierarchical Adversarial Robustness of Deep Neural Networks

abstract: Do all adversarial examples have the same consequences? An autonomous driving system misclassifying a pedestrian as a car may induce a far more dangerous -- and even potentially lethal -- behavior than, for instance, a car as a bus. In order to better tackle this important problematic, we introduce the concept of hierarchical adversarial robustness. Given a dataset whose classes can be grouped into coarse-level labels, we define hierarchical adversarial examples as the ones leading to a misclassification at the coarse level. To improve the resistance of neural networks to hierarchical attacks, we introduce a hierarchical adversarially robust (HAR) network design that decomposes a single classification task into one coarse and multiple fine classification tasks, before being specifically trained by adversarial defense techniques. As an alternative to an end-to-end learning approach, we show that HAR significantly improves the robustness of the network against $\ell_2$ and $\ell_{\infty}$ bounded hierarchical attacks on the CIFAR-10 and CIFAR-100 dataset.


title: A Visibility Roadmap Sampling Approach for a Multi-Robot Visibility-Based Pursuit-Evasion Problem

abstract: Given a two-dimensional polygonal space, the multi-robot visibility-based pursuit-evasion problem tasks several pursuer robots with the goal of establishing visibility with an arbitrarily fast evader. The best known complete algorithm for this problem takes time doubly exponential in the number of robots. However, sampling-based techniques have shown promise in generating feasible solutions in these scenarios. One of the primary drawbacks to employing existing sampling-based methods is that existing algorithms have long execution times and high failure rates for complex environments. This paper addresses that limitation by proposing a new algorithm that takes an environment as its input and returns a joint motion strategy which ensures that the evader is captured by one of the pursuers. Starting with a single pursuer, we sequentially construct Sample-Generated Pursuit-Evasion Graphs to create such a joint motion strategy. This sequential graph structure ensures that our algorithm will always terminate with a solution, regardless of the complexity of the environment. We describe an implementation of this algorithm and present quantitative results that show significant improvement in comparison to the existing algorithm.


title: Coordinated Receding-Horizon Control of Battery Electric Vehicle Speed and Gearshift Using Relaxed Mixed Integer Nonlinear Programming

abstract: In this paper, we propose an approach to coordinated receding-horizon control of vehicle speed and transmission gearshift for automated battery electric vehicles (BEVs) to achieve improved energy efficiency. The introduction of multi-speed transmissions in BEVs creates an opportunity to manipulate the operating point of electric motors under given vehicle speed and acceleration command, thus providing the potential to further improve the energy efficiency. However, co-optimization of vehicle speed and transmission gearshift leads to a mixed integer nonlinear program (MINLP), solving which can be computationally very challenging. In this paper, we propose a novel continuous relaxation technique to treat such MINLPs that makes it possible to compute solutions with conventional nonlinear programming solvers. After analyzing its theoretical properties, we use it to solve the optimization problem involved in coordinated receding-horizon control of BEV speed and gearshift. Through simulation studies, we show that co-optimizing vehicle speed and transmission gearshift can achieve considerably greater energy efficiency than optimizing them sequentially, and the proposed relaxation technique can reduce the online computational cost to a level that is comparable to the time available for real-time implementation.


title: Graded local cohomology modules with respect to the linked ideals

abstract: Let $R=\oplus_{n\in \N_0}R_n$ be a standard graded ring, $M$ be a finitely generated graded $R$-module and $R_+:=\oplus_{n\in \N}R_n$ denotes the irrelevant ideal of $R$. In this paper, considering the new concept of linkage of ideals over a module, we study the graded components $H^i_{\fa}(M)_n$ when $\fa$ is an h-linked ideal over $M$. More precisely, we show that $H^i_{\fa}(M)$ is tame in each of the following cases: \begin{itemize} \item [(i)] $i=f_{\fa}^{R_+}(M)$, the first integer $i$ for which $R_+\nsubseteq \sqrt{0:H^i_{\fa}(M)}$; \item [(ii)] $i=\cd(R_+,M)$, the last integer $i$ for which $H^{i}_{R_+}(M)\neq 0$, and $\fa=\fb+R_+$ where $\fb$ is an h-linked ideal with $R_+$ over $M$. \end{itemize} Also, among other things, we describe the components $H^i_{\fa}(M)_n$ where $\fa$ is radically h-$M$-licci with respect to $R_+$ of length 2.


title: Relativistic quantum fluid with boost invariance

abstract: We study a relativistic fluid with longitudinal boost invariance in a quantum-statistical framework as an example of a solvable non-equilibrium problem. For the free quantum field, we calculate the exact form of the expectation values of the stress-energy tensor and the entropy current. For the stress-energy tensor, we find that a finite value can be obtained only by subtracting the vacuum of the density operator at some fixed proper time \tau_0. As a consequence, the stress-energy tensor acquires non-trivial quantum corrections to the classical free-streaming form.


title: Designing Approximately Optimal Search on Matching Platforms

abstract: We study the design of a decentralized two-sided matching market in which agents' search is guided by the platform. There are finitely many agent types, each with (potentially random) preferences drawn from known type-specific distributions. Equipped with knowledge of these distributions, the platform guides the search process by determining the meeting rate between each pair of types from the two sides. Focusing on symmetric pairwise preferences in a continuum model, we first characterize the unique stationary equilibrium that arises given a feasible set of meeting rates. We then introduce the platform's optimal directed search problem, which involves optimizing meeting rates to maximize equilibrium social welfare. We first show that incentive issues arising from congestion and cannibalization make the design problem fairly intricate. Nonetheless, we develop an efficiently computable search design whose corresponding equilibrium achieves at least 1/4 the social welfare of the optimal design. In fact, our construction always recovers at least 1/4 the first-best social welfare, where agents' incentives are disregarded. Our directed search design is simple and easy-to-implement, as its corresponding bipartite graph consists of disjoint stars. Furthermore, our design implies the platform can substantially limit choice and yet induce an equilibrium with an approximately optimal welfare. Finally, we show that approximation is likely the best we can hope for by establishing that the problem of designing optimal directed search is NP-hard to even approximate beyond a certain constant factor.


title: Growth optimization of TaN for superconducting spintronics

abstract: We have optimized the growth of superconducting TaN thin films on \ch{SiO2} substrates via dc magnetron sputtering and extract a maximum superconducting transition temperature of $T_{\mathrm{c}}=5$ K as well as a maximum critical field $\mu_0H_{\mathrm{c2}}=(13.8\pm0.1)$ T. To investigate the impact of spin-orbit interaction in superconductor/ferromagnet heterostructures, we then analyze the magnetization dynamics of both normal state and superconducting TaN/\ch{Ni80Fe20}(Permalloy, Py)-bilayers as a function of temperature using broadband ferromagnetic resonance (bbFMR) spectroscopy. The phase sensitive detection of the microwave transmission signal is used to quantitatively extract the inverse current-induced torques of the bilayers. The results are compared to our previous study on NbN/Py-bilayers. In the normal state of TaN, we detect a positive damping-like current-induced torque $\sigma_{\mathrm{d}}$ from the inverse spin Hall effect (iSHE) and a small field-like torque $\sigma_{\mathrm{f}}$ attributed to the inverse Rashba-Edelstein effect (iREE) at the TaN/Py-interface. In the superconducting state of TaN, we detect a negative $\sigma_{\mathrm{d}}$ attributed to the quasiparticle mediated inverse spin Hall effect (QMiSHE) and the unexpected manifestation of a large positive field-like $\sigma_{\mathrm{f}}$ of unknown origin matching our previous results for NbN/Py-bilayers.


title: Theory and simulation of shock waves: Entropyproduction and energy conversion

abstract: We have considered a shock wave as a surface of discontinuity and computed the entropy production using non-equilibrium thermodynamics for surfaces. The results from this method, which we call the "Gibbs excess method" (GEM), were compared with results from three alternative methods, all based on the entropy balance in the shock front region, but with different assumptions about local equilibrium. Non-equilibrium molecular dynamics (NEMD) simulations were used to simulate a thermal blast in a one-component gas consisting of particles interacting with the Lennard-Jones/spline potential. This provided data for the theoretical analysis. Two cases were studied, a weak shock with Mach number $M \approx 2$ and a strong shock with $M \approx 6$ and with a Prandtl number of the gas $Pr \approx 1.4$ in both cases. The four theoretical methods gave consistent results for the time-dependent surface excess entropy production for both Mach numbers. The internal energy was found to deviate only slightly from equilibrium values in the shock front. The pressure profile was found to be consistent with the Navier-Stokes equations. The entropy production in the weak and strong shocks were approximately proportional to the square of the Mach number and decayed with time at approximately the same relative rate. In both cases, some 97 \% of the total entropy production in the gas occurred in the shock wave. The GEM showed that most of the shock's kinetic energy was converted reversibly into enthalpy and entropy, and a small amount was dissipated as produced entropy. The shock waves traveled at almost constant speed and we found that the overpressure determined from NEMD simulations agreed well with the Rankine-Hugoniot conditions for steady-state shocks.


title: Linear Nearest Neighbor Flocks with All Distinct Agents

abstract: This paper analyzes the global dynamics of 1-dimensional agent arrays with nearest neighbor linear couplings. The equations of motion are second order linear ODEs with constant coeffcients. The novel part of this research is that the couplings are different for each distinct agent. We allow the forces to depend on the positions and velocity (damping terms) but the magnitudes of both the position and velocity couplings are different for each agent. We, also, do not assume that the forces are "Newtonian" (i.e. the force due to A on B equals the minus the force of B on A) as this assumption does not apply to certain situations, such as traffic modeling. For example, driver A reacting to driver B does not imply the opposite reaction in driver B. There are no known analytical means to solve these systems, even though they are linear, and so relatively little is known about them. This paper is a generalization of previous work that computed the global dynamics of 1-dimensional sequences of identical agents [3] assuming periodic boundary conditions. In this paper, we push that method further, similar to [2], and use an extended periodic boundary condition to to gain quantitative insights to the systems under consideration. We find that we can approximate the global dynamics of such a system by carefully analyzing the low-frequency behavior of the system with (generalized) periodic boundary conditions.


title: Monitoring behavioural responses during pandemic via reconstructed contact matrices from online and representative surveys

abstract: The unprecedented behavioural responses of societies have been evidently shaping the COVID-19 pandemic, yet it is a significant challenge to accurately monitor the continuously changing social mixing patterns in real-time. Contact matrices, usually stratified by age, summarise interaction motifs efficiently, but their collection relies on conventional representative survey techniques, which are expensive and slow to obtain. Here we report a data collection effort involving over $2.3\%$ of the Hungarian population to simultaneously record contact matrices through a longitudinal online and sequence of representative phone surveys. To correct non-representative biases characterising the online data, by using census data and the representative samples we develop a reconstruction method to provide a scalable, cheap, and flexible way to dynamically obtain closer-to-representative contact matrices. Our results demonstrate the potential of combined online-offline data collections to understand the changing behavioural responses determining the future evolution of the outbreak, and inform epidemic models with crucial data.


title: Sample-path large deviations for stochastic evolutions driven by the square of a Gaussian process

abstract: Recently, a number of physical models has emerged described by a random process with increments given by a quadratic form of a fast Gaussian process. We find that the rate function which describes sample-path large deviations for such a process can be computed from the large domain size asymptotic of a certain Fredholm determinant. The latter can be evaluated analytically using a theorem of Widom which generalizes the celebrated Szeg\H{o}-Kac formula to the multi-dimensional case. This provides a large class of random dynamical systems with time scale separation for which an explicit sample-path large deviation functional can be found. Inspired by problems in hydrodynamics and atmosphere dynamics, we construct a simple example with a single slow degree of freedom driven by the square of a fast multi-variate Gaussian process and analyse its large deviations functional using our general results. Even though the noiseless limit of this example has a single fixed point, the corresponding large deviations effective potential has multiple fixed points. In other words, it is the addition of noise that leads to metastability. We use the explicit answers for the rate function to construct instanton trajectories connecting the metastable states.


title: Estimate Three-Phase Distribution Line Parameters With Physics-Informed Graphical Learning Method

abstract: Accurate estimates of network parameters are essential for modeling, monitoring, and control in power distribution systems. In this paper, we develop a physics-informed graphical learning algorithm to estimate network parameters of three-phase power distribution systems. Our proposed algorithm uses only readily available smart meter data to estimate the three-phase series resistance and reactance of the primary distribution line segments. We first develop a parametric physics-based model to replace the black-box deep neural networks in the conventional graphical neural network (GNN). Then we derive the gradient of the loss function with respect to the network parameters and use stochastic gradient descent (SGD) to estimate the physical parameters. Prior knowledge of network parameters is also considered to further improve the accuracy of estimation. Comprehensive numerical study results show that our proposed algorithm yields high accuracy and outperforms existing methods.


title: Deep Learning Approaches for Forecasting Strawberry Yields and Prices Using Satellite Images and Station-Based Soil Parameters

abstract: Computational tools for forecasting yields and prices for fresh produce have been based on traditional machine learning approaches or time series modelling. We propose here an alternate approach based on deep learning algorithms for forecasting strawberry yields and prices in Santa Barbara county, California. Building the proposed forecasting model comprises three stages: first, the station-based ensemble model (ATT-CNN-LSTM-SeriesNet_Ens) with its compound deep learning components, SeriesNet with Gated Recurrent Unit (GRU) and Convolutional Neural Network LSTM with Attention layer (Att-CNN-LSTM), are trained and tested using the station-based soil temperature and moisture data of SantaBarbara as input and the corresponding strawberry yields or prices as output. Secondly, the remote sensing ensemble model (SIM_CNN-LSTM_Ens), which is an ensemble model of Convolutional NeuralNetwork LSTM (CNN-LSTM) models, is trained and tested using satellite images of the same county as input mapped to the same yields and prices as output. These two ensembles forecast strawberry yields and prices with minimal forecasting errors and highest model correlation for five weeks ahead forecasts.Finally, the forecasts of these two models are ensembled to have a final forecasted value for yields and prices by introducing a voting ensemble. Based on an aggregated performance measure (AGM), it is found that this voting ensemble not only enhances the forecasting performance by 5% compared to its best performing component model but also outperforms the Deep Learning (DL) ensemble model found in literature by 33% for forecasting yields and 21% for forecasting prices


title: Maximal operators and Fourier restriction on the moment curve

abstract: We bound certain r-maximal restriction operators on the moment curve.


title: Optimizing Large-Scale Hyperparameters via Automated Learning Algorithm

abstract: Modern machine learning algorithms usually involve tuning multiple (from one to thousands) hyperparameters which play a pivotal role in terms of model generalizability. Black-box optimization and gradient-based algorithms are two dominant approaches to hyperparameter optimization while they have totally distinct advantages. How to design a new hyperparameter optimization technique inheriting all benefits from both approaches is still an open problem. To address this challenging problem, in this paper, we propose a new hyperparameter optimization method with zeroth-order hyper-gradients (HOZOG). Specifically, we first exactly formulate hyperparameter optimization as an A-based constrained optimization problem, where A is a black-box optimization algorithm (such as deep neural network). Then, we use the average zeroth-order hyper-gradients to update hyperparameters. We provide the feasibility analysis of using HOZOG to achieve hyperparameter optimization. Finally, the experimental results on three representative hyperparameter (the size is from 1 to 1250) optimization tasks demonstrate the benefits of HOZOG in terms of simplicity, scalability, flexibility, effectiveness and efficiency compared with the state-of-the-art hyperparameter optimization methods.


title: Robust Model Predictive Path Integral Control: Analysis and Performance Guarantees

abstract: In this paper we propose a novel decision making architecture for Robust Model Predictive Path Integral control (RMPPI) and investigate its performance guarantees and applicability to off-road navigation. Key building blocks of the proposed architecture are an augmented state space representation of the system consisting of nominal and actual dynamics, a placeholder for different types of tracking controllers, a safety logic for nominal state propagation, and an importance sampling scheme that takes into account the capabilities of the underlying tracking control. Using these ingredients, we derive a bound on the free energy growth of the dynamical system which is a function of task constraint satisfaction level, the performance of the underlying tracking controller, and the sampling error of the stochastic optimization used within RMPPI. To validate the bound on free energy growth, we perform experiments in simulation using two types of tracking controllers, namely the iterative Linear Quadratic Gaussian and Contraction-Metric based control. We further demonstrate the applicability of RMPPI in real hardware using the GT AutoRally vehicle. Our experiments demonstrate that RMPPI outperforms MPPI and Tube-MPPI by alleviating issues of the aforementioned model predictive controllers related to either lack of robustness or excessive conservatism. RMPPI provides the best of the two worlds in terms of agility and robustness to disturbances.


title: Kubo Combinatorics for Turbulence Scaling Laws

abstract: We present an extension to Kolmogorov's refined similarity hypothesis for universal fully developed turbulence. The extension is applied within Z. She and E. Leveque's multifractal model of inertial range scaling and its generalizations. Our modification rectifies an apparent gap between the implicit continuum of length scales in Obukhov's conception of a turbulent energy cascade, and scaling law models derived from Kolmogorov's refined similarity hypothesis that lack infinite divisibility. The development has relevance to universal fully developed turbulence, a state we describe explicitly in terms of the coupling between velocity fluctuations and averaged energy dissipation at all orders. This description is unique and leads to a reparametrization of the She-Leveque model that preserves its original forecasts and is infinitely divisible.


title: Joint Continuous and Discrete Model Selection via Submodularity

abstract: In model selection problems for machine learning, the desire for a well-performing model with meaningful structure is typically expressed through a regularized optimization problem. In many scenarios, however, the meaningful structure is specified in some discrete space, leading to difficult nonconvex optimization problems. In this paper, we connect the model selection problem with structure-promoting regularizers to submodular function minimization with continuous and discrete arguments. In particular, we leverage the theory of submodular functions to identify a class of these problems that can be solved exactly and efficiently with an agnostic combination of discrete and continuous optimization routines. We show how simple continuous or discrete constraints can also be handled for certain problem classes and extend these ideas to a robust optimization framework. We also show how some problems outside of this class can be embedded within the class, further extending the class of problems our framework can accommodate. Finally, we numerically validate our theoretical results with several proof-of-concept examples with synthetic and real-world data, comparing against state-of-the-art algorithms.


title: Proactive DP: A Multple Target Optimization Framework for DP-SGD

abstract: We introduce a multiple target optimization framework for DP-SGD referred to as pro-active DP. In contrast to traditional DP accountants, which are used to track the expenditure of privacy budgets, the pro-active DP scheme allows one to a-priori select parameters of DP-SGD based on a fixed privacy budget (in terms of $\epsilon$ and $\delta$) in such a way to optimize the anticipated utility (test accuracy) the most. To achieve this objective, we first propose significant improvements to the moment account method, presenting a closed-form $(\epsilon,\delta)$-DP guarantee that connects all parameters in the DP-SGD setup. We show that DP-SGD is $(\epsilon<0.5,\delta=1/N)$-DP if $\sigma=\sqrt{2(\epsilon +\ln(1/\delta))/\epsilon}$ with $T$ at least $\approx 2k^2/\epsilon$ and $(2/e)^2k^2-1/2\geq \ln(N)$, where $T$ is the total number of rounds, and $K=kN$ is the total number of gradient computations where $k$ measures $K$ in number of epochs of size $N$ of the local data set. We prove that our expression is close to tight in that if $T$ is more than a constant factor $\approx 4$ smaller than the lower bound $\approx 2k^2/\epsilon$, then the $(\epsilon,\delta)$-DP guarantee is violated. The above DP guarantee can be enhanced in thatDP-SGD is $(\epsilon, \delta)$-DP if $\sigma = \sqrt{2(\epsilon+\ln(1/\delta))/\epsilon}$ with $T$ at least $\approx 2k^2/\epsilon$ together with two additional, less intuitive, conditions that allow larger $\epsilon\geq 0.5$. Our DP theory allows us to create a utility graph and DP calculator. These tools link privacy and utility objectives and search for optimal experiment setups, efficiently taking into account both accuracy and privacy objectives, as well as implementation goals. We furnish a comprehensive implementation flow of our proactive DP, with rigorous experiments to showcase the proof-of-concept.


title: On the Convergence of Stochastic Gradient Descent with Bandwidth-based Step Size

abstract: We investigate the stochastic gradient descent (SGD) method where the step size lies within a banded region instead of being given by a fixed formula. The optimal convergence rate under mild conditions and large initial step size is proved. Our analysis provides comparable theoretical error bounds for SGD associated with a variety of step sizes. In addition, the convergence rates for some existing step size strategies, e.g., triangular policy and cosine-wave, can be revealed by our analytical framework under the boundary constraints. The bandwidth-based step size provides efficient and flexible step size selection in optimization. We also propose a $1/t$ up-down policy and give several non-monotonic step sizes. Numerical experiments demonstrate the efficiency and significant potential of the bandwidth-based step-size in many applications.


title: Consistent Lock-free Parallel Stochastic Gradient Descent for Fast and Stable Convergence

abstract: Stochastic gradient descent (SGD) is an essential element in Machine Learning (ML) algorithms. Asynchronous parallel shared-memory SGD (AsyncSGD), including synchronization-free algorithms, e.g. HOGWILD!, have received interest in certain contexts, due to reduced overhead compared to synchronous parallelization. Despite that they induce staleness and inconsistency, they have shown speedup for problems satisfying smooth, strongly convex targets, and gradient sparsity. Recent works take important steps towards understanding the potential of parallel SGD for problems not conforming to these strong assumptions, in particular for deep learning (DL). There is however a gap in current literature in understanding when AsyncSGD algorithms are useful in practice, and in particular how mechanisms for synchronization and consistency play a role. We focus on the impact of consistency-preserving non-blocking synchronization in SGD convergence, and in sensitivity to hyper-parameter tuning. We propose Leashed-SGD, an extensible algorithmic framework of consistency-preserving implementations of AsyncSGD, employing lock-free synchronization, effectively balancing throughput and latency. We argue analytically about the dynamics of the algorithms, memory consumption, the threads' progress over time, and the expected contention. We provide a comprehensive empirical evaluation, validating the analytical claims, benchmarking the proposed Leashed-SGD framework, and comparing to baselines for training multilayer perceptrons (MLP) and convolutional neural networks (CNN). We observe the crucial impact of contention, staleness and consistency and show how Leashed-SGD provides significant improvements in stability as well as wall-clock time to convergence (from 20-80% up to 4x improvements) compared to the standard lock-based AsyncSGD algorithm and HOGWILD!, while reducing the overall memory footprint.


title: The Galactic HII Region Luminosity Function at Radio and Infrared Wavelengths

abstract: The Galactic HII region luminosity function (LF) is an important metric for understanding global star formation properties of the Milky Way, but only a few studies have been done and all use relatively small numbers of HII regions. We use a sample of 797 first Galactic quadrant HII regions compiled from the WISE Catalog of Galactic HII Regions to examine the form of the LF at multiple infrared and radio wavelengths. Our sample is statistically complete for all regions powered by single stars of type O9.5V and earlier. We fit the LF at each wavelength with single and double power laws. Averaging the results from all wavelengths, the mean of the best-fit single power law index is $\langle\alpha\rangle=-1.75\,\pm\,0.01$. The mean best-fit double power law indices are $\langle\alpha_1\rangle=-1.40\,\pm\,0.03$ and $\langle\alpha_2\rangle=-2.33\,\pm\,0.04$. We conclude that neither a single nor a double power law is strongly favored over the other. The LFs show some variation when we separate the HII region sample into subsets by heliocentric distance, physical size, Galactocentric radius, and location relative to the spiral arms, but blending individual HII regions into larger complexes does not change the value of the power law indices of the best-fit LF models. The consistency of the power law indices across multiple wavelengths suggests that the LF is independent of wavelength. This implies that infrared and radio tracers can be employed in place of H$\alpha$.


title: On intrinsic negative curves

abstract: Let $\mathbb K$ be an algebraically closed field of characteristic $0$. A curve of $(\mathbb K^*)^2$ arising from a Laurent polynomial in two variables is {\em intrinsic negative} if its tropical compactification has negative self-intersection. The aim of this note is to start a systematic study of these curves and to relate them with the problem of computing Seshadri constants of toric surfaces.


title: Resolution analysis of inverting the generalized $N$-dimensional Radon transform in $\mathbb R^n$ from discrete data

abstract: Let $\mathcal R$ denote the generalized Radon transform (GRT), which integrates over a family of $N$-dimensional smooth submanifolds $\mathcal S_{\tilde y}\subset\mathcal U$, $1\le N\le n-1$, where an open set $\mathcal U\subset\mathbb R^n$ is the image domain. The submanifolds are parametrized by points $\tilde y\subset\tilde{\mathcal V}$, where an open set $\tilde{\mathcal V}\subset\mathbb R^n$ is the data domain. The continuous data are $g={\mathcal R} f$, and the reconstruction is $\check f=\mathcal R^*\mathcal B g$. Here $\mathcal R^*$ is a weighted adjoint of $\mathcal R$, and $\mathcal B$ is a pseudo-differential operator. We assume that $f$ is a conormal distribution, $\text{supp}(f)\subset\mathcal U$, and its singular support is a smooth hypersurface $\mathcal S\subset\mathcal U$. Discrete data consists of the values of $g$ on a lattice $\tilde y^j$ with the step size $O(\epsilon)$. Let $\check f_\epsilon=\mathcal R^*\mathcal B g_\epsilon$ denote the reconstruction obtained by applying the inversion formula to an interpolated discrete data $g_\epsilon(\tilde y)$. Pick a generic pair $(x_0,\tilde y_0)$, where $x_0\in\mathcal S$, and $\mathcal S_{\tilde y_0}$ is tangent to $\mathcal S$ at $x_0$. The main result of the paper is the computation of the limit $$ f_0(\check x):=\lim_{\epsilon\to0}\epsilon^\kappa \check f_\epsilon(x_0+\epsilon\check x). $$ Here $\kappa\ge 0$ is selected based on the strength of the reconstructed singularity, and $\check x$ is confined to a bounded set. The limiting function $f_0(\check x)$, which we call the discrete transition behavior, allows computing the resolution of reconstruction.


title: Coordinated tractions increase the size of a collectively moving pack in a cell monolayer

abstract: Cells in an epithelial monolayer coordinate motion with their neighbors giving rise to collectively moving packs of sizes spanning multiple cell diameters. The physical mechanism controlling the pack size, however, remains unclear. A potential mechanism comes from assuming that cell-substrate traction forces persist over some time scale: with large enough persistence time, collective cell packs emerge. To test this hypothesis, we measured the velocity and net traction of each cell. The data showed that in addition to having some temporal persistence, tractions were spatially correlated, suggesting that cells coordinate with their neighbors to apply tractions in the same direction. Chemical inhibitors and activators of actomyosin contraction were used to determine effects of altering the traction persistence and alignment. Numerical simulations based on the self-propelled Voronoi model, augmented to include both traction persistence and alignment and calibrated against the experimental data, matched the experimentally measured pack size. The model identified that if there were no alignment of traction between neighboring cells, the size of the collective pack would be substantially smaller than observed in the experiments. Hence, combining experiments and a simple mechanical model, this study confirms the long-standing assumption of traction persistence and adds the notion of traction alignment between neighbors. Together, persistence and alignment are two factors controlling the size of a collectively moving cell pack.


title: Surface Densities Prewet a Near-Critical Membrane

abstract: Recent work has highlighted roles for thermodynamic phase behavior in diverse cellular processes. Proteins and nucleic acids can phase separate into three-dimensional liquid droplets in the cytoplasm and nucleus and the plasma membrane of animal cells appears tuned close to a two-dimensional liquid-liquid critical point. In some examples, cytoplasmic proteins aggregate at plasma membrane domains, forming structures such as the post-synaptic density and diverse signaling clusters. Here we examine the physics of these surface densities, employing minimal simulations of co-acervating polymers coupled to an Ising membrane surface in conjunction with a complementary Landau theory. We argue that these surface densities are a novel phase reminiscent of pre-wetting, in which a molecularly thin three-dimensional liquid forms on a usually solid surface. However, in surface densities the solid surface is replaced by a membrane with an independent propensity to phase separate. We show that proximity to criticality in the membrane dramatically increases the parameter regime in which a pre-wetting-like transition occurs, leading to a broad region where coexisting surface phases can form even when a bulk phase is unstable. Our simulations naturally exhibit three surface phase coexistence even though both the membrane and the polymer bulk can only display two phase coexistence on their own. We argue that the physics of these surface densities enables diverse functions seen in Eukaryotic cells.


title: On robustly convergent and efficient iterative methods for anisotropic radiative transfer

abstract: This paper considers the iterative solution of linear systems arising from discretization of the anisotropic radiative transfer equation with discontinuous elements on the sphere. In order to achieve robust convergence behavior in the discretization parameters and the physical parameters we develop preconditioned Richardson iterations in Hilbert spaces. We prove convergence of the resulting scheme. The preconditioner is constructed in two steps. The first step borrows ideas from matrix splittings and ensures mesh independence. The second step uses a subspace correction technique to reduce the influence of the optical parameters. The correction spaces are build from low-order spherical harmonics approximations generalizing well-known diffusion approximations. We discuss in detail the efficient implementation and application of the discrete operators. In particular, for the considered discontinuous spherical elements, the scattering operator becomes dense and we show that $\mathcal{H}$- or $\mathcal{H}^2$-matrix compression can be applied in a black-box fashion to obtain almost linear or linear complexity when applying the corresponding approximations. The effectiveness of the proposed method is shown in numerical examples.


title: Spacewalker: Rapid UI Design Exploration Using Lightweight Markup Enhancement and Crowd Genetic Programming

abstract: User interface design is a complex task that involves designers examining a wide range of options. We present Spacewalker, a tool that allows designers to rapidly search a large design space for an optimal web UI with integrated support. Designers first annotate each attribute they want to explore in a typical HTML page, using a simple markup extension we designed. Spacewalker then parses the annotated HTML specification, and intelligently generates and distributes various configurations of the web UI to crowd workers for evaluation. We enhanced a genetic algorithm to accommodate crowd worker responses from pairwise comparison of UI designs, which is crucial for obtaining reliable feedback. Based on our experiments, Spacewalker allows designers to effectively search a large design space of a UI, using the language they are familiar with, and improve their design rapidly at a minimal cost.


title: Scaling limits of external multi-particle DLA on the plane and the supercooled Stefan problem

abstract: We consider (a variant of) the external multi-particle diffusion-limited aggregation (MDLA) process of Rosenstock and Marquardt on the plane. Based on the recent findings of [11], [10] in one space dimension it is natural to conjecture that the scaling limit of the growing aggregate in such a model is given by the growing solid phase in a suitable "probabilistic" formulation of the single-phase supercooled Stefan problem for the heat equation. To address this conjecture, we extend the probabilistic formulation from [10] to multiple space dimensions. We then show that the equation that characterizes the growth rate of the solid phase in the supercooled Stefan problem is satisfied by the scaling limit of the external MDLA process with an inequality, which can be strict in general. In the course of the proof, we establish two additional results interesting in their own right: (i) the stability of a "crossing property" of planar Brownian motion and (ii) a rigorous connection between the probabilistic solutions to the supercooled Stefan problem and its classical and weak solutions.


title: Reaching Consensus for Asynchronous Distributed Key Generation

abstract: We give a protocol for Asynchronous Distributed Key Generation (A-DKG) that is optimally resilient (can withstand $f<\frac{n}{3}$ faulty parties), has a constant expected number of rounds, has $\tilde{O}(n^3)$ expected communication complexity, and assumes only the existence of a PKI. Prior to our work, the best A-DKG protocols required $\Omega(n)$ expected number of rounds, and $\Omega(n^4)$ expected communication. Our A-DKG protocol relies on several building blocks that are of independent interest. We define and design a Proposal Election (PE) protocol that allows parties to retrospectively agree on a valid proposal after enough proposals have been sent from different parties. With constant probability the elected proposal was proposed by a non-faulty party. In building our PE protocol, we design a Verifiable Gather protocol which allows parties to communicate which proposals they have and have not seen in a verifiable manner. The final building block to our A-DKG is a Validated Asynchronous Byzantine Agreement (VABA) protocol. We use our PE protocol to construct a VABA protocol that does not require leaders or an asynchronous DKG setup. Our VABA protocol can be used more generally when it is not possible to use threshold signatures.


title: Modeling Extremes with d-max-decreasing Neural Networks

abstract: We propose a novel neural network architecture that enables non-parametric calibration and generation of multivariate extreme value distributions (MEVs). MEVs arise from Extreme Value Theory (EVT) as the necessary class of models when extrapolating a distributional fit over large spatial and temporal scales based on data observed in intermediate scales. In turn, EVT dictates that $d$-max-decreasing, a stronger form of convexity, is an essential shape constraint in the characterization of MEVs. As far as we know, our proposed architecture provides the first class of non-parametric estimators for MEVs that preserve these essential shape constraints. We show that our architecture approximates the dependence structure encoded by MEVs at parametric rate. Moreover, we present a new method for sampling high-dimensional MEVs using a generative model. We demonstrate our methodology on a wide range of experimental settings, ranging from environmental sciences to financial mathematics and verify that the structural properties of MEVs are retained compared to existing methods.


title: Astrometric Membership Tests for the Zinn-Newell-Gibson "UV-Bright" Stars in Galactic Globular Clusters

abstract: In 1972, Zinn, Newell, & Gibson (ZNG) published a list of 156 candidate "UV-bright" stars they had found in 27 Galactic globular clusters (GCs), based on photographs in the U and V bands. UV-bright stars lie above the horizontal branch (HB) and blueward of the asymptotic-giant branch (AGB) and red-giant branch in the clusters' color-magnitude diagrams. They are in rapid evolutionary phases--if they are members and not unrelated bright foreground stars. The ZNG list has inspired numerous follow-up studies, aimed at understanding late stages of stellar evolution. However, the ZNG candidates were presented only in finding charts, and celestial coordinates were not given. Using my own collection of CCD frames in u and V, I have identified all of the ZNG objects, and have assembled their coordinates, parallaxes, and proper motions from the recent Gaia Early Data Release 3 (EDR3). Based on the Gaia astrometry, I have determined which objects are probable cluster members (45% of the sample). For the members, using photometry from EDR3, I have assigned the stars to various evolutionary stages, including luminous post-AGB stars, and stars above the HB. I point out several ZNG stars of special interest that have still, to my knowledge, never been studied in detail. This study is an adjunct to a forthcoming survey of the Galactic GCs in the uBVI photometric system, designed for detection of low-gravity stars with large Balmer discontinuities.


title: Balanced Operators and Operator Domains

abstract: We shall say that a densely defined closed operator $T$ on a Hilbert space is balanced if $\cD(T)=\cD(T^*)$. Balanced operators are described in terms of their phase operators abnd their moduli. Examples of balanced operators are developed. A characterization of the domain equality $\cD(A)=\cD(B)$ for positive self-adjoint operators $A$ and $B$ with bounded inverses is given in terms of their spectral measures.


title: Effect of bulk viscosity in cosmic acceleration

abstract: In this paper, we have investigated the late time cosmic acceleration issue in the context of $f(R,T)$ gravity. The matter field is considered to be that of viscous fluid. The model has been framed as a mathematical formalism and the effect of viscous fluid on the cosmic expansion has been shown. The equation of state parameter indicates the quintessence behaviour of the Universe at late time. The theoretical results obtained here shows its alignment with the cosmological observations result.