title: Disconnecting the Dots: Re-Examining the Nature of Stellar "Strings" in the Milky Way

abstract: Recent analyses of Gaia data have resulted in the identification of new stellar structures, including a new class of extended stellar filaments called stellar "strings", first proposed by Kounkel and Covey 2019. We explore the spatial, kinematic, and chemical composition of strings to demonstrate that these newfound structures are largely inconsistent with being physical objects whose members share a common origin. Examining the 3D spatial distribution of string members, we find that the spatial dispersion around the claimed string spine does not improve in the latest Gaia DR3 data release -- despite tangible gains in the signal-to-noise of the parallax measurements -- counter to expectations of a bona fide structure. Using the radial velocity dispersion of the strings (averaging $\rm\sigma_{V_r} = 16\;km\;s^{-1}$) to estimate their virial masses, we find that all strings are gravitationally unbound. Given the finding that the strings are dispersing, the reported stellar ages of the strings are typically $120\times$ larger than their measured dispersal times. Finally, we validate prior work that stellar strings are more chemically homogeneous than their local field stars, but show it is possible to obtain the same signatures of chemical homogeneity by drawing random samples of stars from spatially, temporally, and kinematically unrelated open clusters. Our results show that while some strings may be composed of real sub-structures, there is no consistent evidence for larger string-like connections over the sample. These results underline the need for caution in over-interpreting the significance of these strings and their role in understanding the star formation history of the Milky Way.


title: Implications of a Temperature Dependent IMF II: An Updated View of the Star-Forming Main Sequence

abstract: The stellar initial mass function (IMF) is predicted to depend upon the temperature of gas in star-forming molecular clouds. The introduction of an additional parameter, $T_{IMF}$ , into photometric template fitting, allows galaxies to be fit with a range of IMFs. Three surprising new features appear: (1) most star-forming galaxies are best fit with a bottom-lighter IMF than the Milky Way; (2) most star-forming galaxies at fixed redshift are fit with a very similar IMF; and (3) the most massive star-forming galaxies at fixed redshift instead exhibit a less bottom-light IMF, similar to that measured in quiescent galaxies. Additionally, since stellar masses and star formation rates both depend on the IMF, these results slightly modify the resulting relationship, while yielding similar qualitative characteristics to previous studies.


title: Simultaneous X-ray and optical spectroscopy of V404 Cygni supports the multi-phase nature of X-ray binary accretion disc winds

abstract: Observational signatures of accretion disc winds have been found in a significant number of low-mass X-ray binaries at either X-ray or optical wavelengths. The 2015 outburst of the black hole transient V404 Cygni provided a unique opportunity for studying both types of outflows in the same system. We used contemporaneous X-ray (Chandra Observatory) and optical (Gran Telescopio Canarias, GTC) spectroscopy, in addition to hard X-ray light curves (INTEGRAL). We show that the kinetic properties of the wind, as derived from P-Cyg profiles detected in the optical range at low hard X-ray fluxes and in a number of X-ray transitions during luminous flares, are remarkably similar. Furthermore, strictly simultaneous data taken at intermediate hard X-ray fluxes show consistent emission line properties between the optical and the X-ray emission lines, which most likely arise in the same accretion disc wind. We discuss several scenarios to explain the properties of the wind, favouring the presence of a dynamic, multi-phase outflow during the entire outburst of the system. This study, together with the growing number of wind detections with fairly similar characteristic velocities at different wavelengths, suggest that wind-type X-ray binary outflows might be predominantly multi-phase in nature.


title: Finding lensed direct-collapse black holes and supermassive primordial stars

abstract: Direct-collapse black holes (DCBHs) may be the seeds of the first quasars, over 200 of which have now been detected at $z > 6$ . The James Webb Space Telescope (JWST) could detect DCBHs in the near infrared (NIR) at $z \lesssim 20$ and probe the evolution of primordial quasars at their earliest stages, but only in narrow fields that may not capture many of them. Wide-field NIR surveys by Euclid and the Nancy Grace Roman Survey Telescope (RST) would enclose far greater numbers of DCBHs but only directly detect them at $z \lesssim 6 - 8$ because of their lower sensitivities. However, their large survey areas will cover thousands of galaxy clusters and massive galaxies that could gravitationally lense flux from DCBHs, boosting them above current Euclid and RST detection limits and revealing more of them than could otherwise be detected. Here, we estimate the minimum number density of strongly lensed DCBHs and supermassive primordial stars required for detection in surveys by Euclid, RST and JWST at $z \lesssim 20$. We find that for reasonable estimates of host halo numbers RST, Euclid, and JWST could potentially find hundreds of strongly-lensed DCBHs at $z = 7 - 20$. RST would detect the most objects at $z \lesssim 10$ and JWST would find the most at higher redshifts. Lensed supermassive primordial stars could potentially also be found, but in fewer numbers because of their short lifetimes.


title: Predictive wavefront control on Keck II adaptive optics bench: on-sky coronagraphic results

abstract: The behavior of an adaptive optics (AO) system for ground-based high contrast imaging (HCI) dictates the achievable contrast of the instrument. In conditions where the coherence time of the atmosphere is short compared to the speed of the AO system, the servo-lag error can become the dominant error term of the AO system. While the AO system measures the wavefront error and subsequently applies a correction (typically taking a total of one or a few milliseconds), the atmospheric turbulence above the telescope has changed resulting in the servo-lag error. In addition to reducing the Strehl ratio, the servo-lag error causes a build-up of speckles along the direction of the dominant wind vector in the coronagraphic image, severely limiting the contrast at small angular separations. One strategy to mitigate this problem is to predict the evolution of the turbulence over the delay time. Our predictive wavefront control algorithm minimizes, in a mean square sense, the wavefront error over the delay and has been implemented on the Keck II AO bench. In this paper, we report on the latest results of our algorithm and discuss updates to the algorithm itself. We explore how to tune various filter parameters based on both daytime laboratory tests and on-sky tests. We show a reduction in the residual-mean-square wavefront error for the predictor compared to the leaky integrator (the standard controller for Keck) implemented on Keck for three separate nights. Finally, we present contrast improvements for daytime and on-sky tests for the first time. Using the L-band vortex coronagraph for Keck's NIRC2 instrument, we find a contrast gain of up to 2 at a separation of 3 lambda/D and up to 3 for larger separations (3-7 lambda/D).


title: The Great Dimming of Betelgeuse seen by the Himawari-8 meteorological satellite

abstract: Betelgeuse, one of the most studied red supergiant stars, dimmed in the optical by ~1.2 mag between late 2019 and early 2020, reaching an historical minimum called "the Great Dimming." Thanks to enormous observational effort to date, two hypotheses remain that can explain the Dimming: a decrease in the effective temperature and an enhancement of the extinction caused by newly produced circumstellar dust. However, the lack of multi-wavelength monitoring observations, especially in the mid infrared where emission from circumstellar dust can be detected, has prevented us from closely examining these hypotheses. Here we present 4.5-year, 16-band photometry of Betelgeuse between 2017-2021 in the 0.45-13.5 micron wavelength range making use of images taken by the Himawari-8 geostationary meteorological satellite. By examining the optical and near-infrared light curves, we show that both a decreased effective temperature and increased dust extinction may have contributed by almost the same amount to the Great Dimming. Moreover, using the mid-infrared light curves, we find that the enhanced circumstellar extinction actually contributed to the Dimming. Thus, the Dimming event of Betelgeuse provides us an opportunity to examine the mechanism responsible for the mass loss of red supergiants, which affects the fate of massive stars as supernovae.


title: Scalar field damping at high temperatures

abstract: The motion of a scalar field that interacts with a hot plasma, like the inflaton during reheating, is damped, which is a dissipative process. At high temperatures the damping can be described by a local term in the effective equation of motion. The damping coefficient is sensitive to multiple scattering. In the loop expansion its computation would require an all-order resummation. Instead we solve an effective Boltzmann equation, similarly to the computation of transport coefficients. For an interaction with another scalar field we obtain a simple relation between the damping coefficient and the bulk viscosity, so that one can make use of known results for the latter. The numerical prefactor of the damping coefficient turns out to be rather large, of order $ 10 ^ 4 $.


title: The Covariance of Photometric and Spectroscopic Two-Point Statistics: Implications for Cosmological Parameter Inference

abstract: To combine information from measurements of the redshift-space power spectrum from spectroscopic data with angular weak lensing, galaxy clustering and galaxy-galaxy lensing power spectra from photometric surveys (i.e. the $3 \times 2$ point statistics), we must account for the covariance between the two probes. Currently any covariance between the two types of measurements is neglected as existing photometric and spectroscopic surveys largely probe different cosmological volumes. This will cease to be the case as data arrives from Stage-IV surveys. In this paper we derive an analytic expression for the covariance between photometric 2D angular power spectra and the 3D redshift-space power spectrum for Gaussian fields under the plane-parallel approximation. We find that the two probes are covariant on large radial scales, but because the information content of these modes is extremely low due to sample variance, we forecast that it is safe to neglect this covariance when performing cosmological parameter inference.


title: Signatures of reionization feedback in the near-infrared background

abstract: The reionization of the intergalactic medium at redshifts $z\gtrsim 6$ is expected to have a lasting impact on galaxies residing in low-mass dark matter halos. Unable to accrete or retain gas photo-heated to temperatures $T \gtrsim 10^4$ K, the star formation histories of faint galaxies in the early Universe are expected to decline as they exhaust their gas supply, resulting in a turn-over in the galaxy luminosity function (LF) and a potential solution to the missing satellites problem in the local group. Unfortunately, there are several challenges to constraining `reionization feedback' empirically, most notably that galaxies in low-mass halos are intrinsically faint, and that there are other physical mechanisms capable of inducing a turn-over in the LF. In this work, we investigate a new signature of reionization feedback that is in principle distinct from other processes: as faint galaxies are quenched by reionization, their stellar populations passively age and grow redder while the brighter galaxies nearby continue to form stars at an increasing rate and so remain relatively blue. We find that this contrast, between quenched and un-quenched galaxies induces a scale and colour-dependent signature in the present-day near-infrared background comparable to the expected sensitivity of NASA's upcoming SPHEREx mission. Whereas models with pure mass suppression largely affect the signal at wavelengths $\lesssim 2 \mu\rm{m}$, $\sim 5$%-level differences in the background persist out to $\sim 5 \mu\rm{m}$ for reionization feedback models. Finally, the power spectra of intensity ratio maps exhibit larger variations, and may thus be a promising target for future analyses.


title: Phase-transition-like behavior in information retrieval of a quantum scrambled random circuit system

abstract: Information in a chaotic quantum system will scramble across the system, preventing any local measurement from reconstructing it. The scrambling dynamics is key to understanding a wide range of quantum many-body systems. Here we use Holevo information to quantify the scrambling dynamics, which shows a phase-transition-like behavior. When applying long random Clifford circuits to a large system, no information can be recovered from a subsystem of less than half the system size. When exceeding half the system size, the amount of stored information grows by two bits of classical information per qubit until saturation through another sharp unanalytical change. We also study critical behavior near the transition points. Finally, we use coherent information to quantify the scrambling of quantum information in the system, which shows similar phase-transition-like behavior.


title: Echoes from backreacting exotic compact objects

abstract: The possible detection of echoes in late gravitational-wave signals is the most promising way to test horizonless alternatives to general relativistic black holes, and probe the physics of these hypothetical ultra-compact objects. While there is currently no evidence for the presence of such signatures, better accuracy is expected with the growing wealth of data from gravitational waves observatories. So far, several searches for these specific signals have been performed considering equidistant intervals between consecutive echoes, i.e. quasi-periodic wave-forms, and ignoring possible backreaction effects of the incoming waves. Here we study scalar perturbations in exotic compact object scenarios that account for possible backreaction phenomena. In particular, we find that if one considers the increase of the central object mass due to the partial absorption of the energy carried by the perturbation, the echo signal can be quite different and non-periodic. Apart from this simple scenario, we also consider the case in which, in order to preserve its compactness above the black hole limit, the compact object absorption shuts down in a finite amount of time or leads to an expansion. In both these cases we find interesting new features that should be taken into account in future searches.


title: Lamellar fluctuations melt ferroelectricity

abstract: We consider a standard Ginzburg-Landau model of a ferroelectric whose electrical polarization is coupled to gradients of elastic strain. At the harmonic level, such flexoelectric interaction is known to hybridize acoustic and optic phonon modes and lead to phases with modulated lattice structures that precede the symmetry broken state for sufficiently large couplings. Here, we use the we use the self-consistent phonon approximation to calculate the effects of thermal and quantum polarization fluctuations on the bare modes to show that such long-range modulated order is unstable at all temperatures. We discuss the implications for the nearly ferroelectric SrTiO$_3$ and KTaO$_3$, and propose that these systems are melted versions of an underlying modulated state which is dominated by non-zero momentum thermal fluctuations except at the very lowest temperatures.


title: Infrared Spectral Energy Distribution and Variability of Active Galactic Nuclei: Clues to the Structure of Circumnuclear Material

abstract: The active galactic nucleus (AGN) phenomena results from a supermassive black hole accreting its surrounding gaseous and dusty material. The infrared (IR) regime provides most of the information to characterize the dusty structures that bridge from the galaxy to the black hole, providing clues to the black hole growth and host galaxy evolution. Over the past several decades, with the commissioning of various ground, airborne and space IR observing facilities, our interpretations of the AGN circumnuclear structures have advanced significantly through improved understanding of how their dust emission changes as a function of wavelength and how the heating of the dusty structures responds to variations of the energy released from the central engine. In this review, we summarize the current observational knowledge of the AGN IR broad-band spectral energy distributions (SEDs) and the IR time variability behavior covering large ranges of AGN luminosity and redshift, and discuss some first-order insights into the obscuring structures and host galaxy IR properties that can be obtained by integrating the relevant observations into a coherent picture.


title: Momentum Stiefel Optimizer, with Applications to Suitably-Orthogonal Attention, and Optimal Transport

abstract: The problem of optimization on Stiefel manifold, i.e., minimizing functions of (not necessarily square) matrices that satisfy orthogonality constraints, has been extensively studied. Yet, a new approach is proposed based on, for the first time, an interplay between thoughtfully designed continuous and discrete dynamics. It leads to a gradient-based optimizer with intrinsically added momentum. This method exactly preserves the manifold structure but does not require additional operation to keep momentum in the changing (co)tangent space, and thus has low computational cost and pleasant accuracy. Its generalization to adaptive learning rates is also demonstrated. Notable performances are observed in practical tasks. For instance, we found that placing orthogonal constraints on attention heads of trained-from-scratch Vision Transformer [Dosovitskiy et al. 2022] could markedly improve its performance, when our optimizer is used, and it is better that each head is made orthogonal within itself but not necessarily to other heads. This optimizer also makes the useful notion of Projection Robust Wasserstein Distance [Paty & Cuturi 2019; Lin et al. 2020] for high-dim. optimal transport even more effective.

categories: cs.LG cs.NA math.DS math.NA math.OC stat.ML


title: Private and Byzantine-Proof Cooperative Decision-Making

abstract: The cooperative bandit problem is a multi-agent decision problem involving a group of agents that interact simultaneously with a multi-armed bandit, while communicating over a network with delays. The central idea in this problem is to design algorithms that can efficiently leverage communication to obtain improvements over acting in isolation. In this paper, we investigate the stochastic bandit problem under two settings - (a) when the agents wish to make their communication private with respect to the action sequence, and (b) when the agents can be byzantine, i.e., they provide (stochastically) incorrect information. For both these problem settings, we provide upper-confidence bound algorithms that obtain optimal regret while being (a) differentially-private and (b) tolerant to byzantine agents. Our decentralized algorithms require no information about the network of connectivity between agents, making them scalable to large dynamic systems. We test our algorithms on a competitive benchmark of random graphs and demonstrate their superior performance with respect to existing robust algorithms. We hope that our work serves as an important step towards creating distributed decision-making systems that maintain privacy.


title: Kinetic theory of {\it tilted} Dirac cone materials

abstract: We formulate the Boltzmann kinetic equations for interacting tilted Dirac fermions in two space dimensions characterized by a tilt parameter $0\le\zeta<1$. Solving the linearized Boltzmann equation, we find that the broadening of the Drude pole is enhanced by $\kappa(\zeta)\times(1-\zeta^2)^{-1/2}$, where the $\kappa$ is interaction-induced enhancement factor. The intensity of the Drude pole is also anisotropically enhanced by $(1-\zeta^2)^{-1}$. The ubiquitous "redshift" factors $(1-\zeta^2)^{1/2}$ can be regarded as a manifestation of an underlying spacetime structure in such solids. The additional broadening $\kappa$ indicates that interaction effects are more pronounced for electrons in a $\zeta$-deformed Minkowski spacetime of tilted Dirac fermions.


title: Orthotropic metamaterials with freely tailorable elastic constants

abstract: An orthotropic metamaterial is composed of elements arrayed periodically in space. The element includes two cuboid structures. The first structure is the basic structure of the element, and the second structure is the transformation of the first structure of the element. The first structure of the element is a cuboid structure composed of 24 bars connected by 8 nodes, and the second structure of the element is a cuboid structure composed of 36 bars connected by 14 nodes. This metamaterial has 6 independent elastic constants, so there is a large degree of freedom in material design. Using a simple universal design method, a metamaterial with tailored elastic constants can be designed. Therefore, it has great application value in the fields of mechanical metamaterials, elastic wave metamaterials, acoustic metamaterials, and seismic metamaterials, and has also laid the foundation for realizing the dream of controlling elastic waves, acoustic waves and vibrations.


title: Multiscale Voxel Based Decoding For Enhanced Natural Image Reconstruction From Brain Activity

abstract: Reconstructing perceived images from human brain activity monitored by functional magnetic resonance imaging (fMRI) is hard, especially for natural images. Existing methods often result in blurry and unintelligible reconstructions with low fidelity. In this study, we present a novel approach for enhanced image reconstruction, in which existing methods for object decoding and image reconstruction are merged together. This is achieved by conditioning the reconstructed image to its decoded image category using a class-conditional generative adversarial network and neural style transfer. The results indicate that our approach improves the semantic similarity of the reconstructed images and can be used as a general framework for enhanced image reconstruction.


title: Temperature Screening and Cross-Field Impurity Accumulation from a Thermodynamic Perspective

abstract: In a variety of different systems, high-Z ion species show a marked tendency to accumulate in regions of high plasma density. It has previously been suggested that the apparent universality of this behavior could be explained thermodynamically, in terms of the maximum-entropy state attainable when the system must obey an ambipolarity condition. However, the previous treatment did not allow for the possibility of temperature gradients. Here, tools from non-equilibrium (Onsager) thermodynamics are used to show that ambipolarity continues to play a key role in producing this behavior in the presence of temperature gradients, and to recover well-known temperature screening effects that appear in these cases.


title: Real-time equation-of-motion CC cumulant and CC Green's function simulations of photoemission spectra of water and water dimer

abstract: Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses results obtained with the real-time equation-of-motion CC cumulant approach (RT-EOM-CC), and CC Green's function (CCGF) approaches in applications to the water and water dimer molecules. We compare the ionization potentials obtained with these methods for the valence region with the results obtained with the CCSD(T) formulation as a difference of energies for N and N-1 electron systems. All methods show good agreement with each other. They also agree well with experiment, with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.


title: Effect of matrix sparsity and quantum noise on quantum random walk linear solvers

abstract: We study the effects of quantum noise in hybrid quantum-classical solver for sparse systems of linear equations using quantum random walks, applied to stoquastic Hamiltonian matrices. In an ideal noiseless quantum computer, sparse matrices achieve solution vectors with lower relative error than dense matrices. However, we find quantum noise reverses this effect, with overall error increasing as sparsity increases. We identify invalid quantum random walks as the cause of this increased error and propose a revised linear solver algorithm which improves accuracy by mitigating these invalid walks.


title: Direction and Trajectory Tracking Control for Nonholonomic Spherical Robot by Combining Sliding Mode Controller and Model Prediction Controller

abstract: Spherical robot is a nonlinear, nonholonomic and unstable system which increases the difficulty of the direction and trajectory tracking problem. In this study, we propose a new direction controller HTSMC, an instruction planning controller MPC, and a trajectory tracking framework MHH. The HTSMC is designed by integrating a fast terminal algorithm, a hierarchical method, the motion features of a spherical robot, and its dynamics. In addition, the new direction controller has an excellent control effect with a quick response speed and strong stability. MPC can obtain optimal commands that are then transmitted to the velocity and direction controller. Since the two torque controllers in MHH are all Lyapunov-based sliding mode controllers, the MHH framework may achieve optimal control performance while assuring stability. Finally, the two controllers eliminate the requirement for MPC's stability and dynamic constraints. Finally, hardware experiments demonstrate the efficacy of the HTSMC, MPC, and MHH.


title: Who is we? Disambiguating the referents of first person plural pronouns in parliamentary debates

abstract: This paper investigates the use of first person plural pronouns as a rhetorical device in political speeches. We present an annotation schema for disambiguating pronoun references and use our schema to create an annotated corpus of debates from the German Bundestag. We then use our corpus to learn to automatically resolve pronoun referents in parliamentary debates. We explore the use of data augmentation with weak supervision to further expand our corpus and report preliminary results.


title: Characterizing the Variable X-ray and UV-Optical Flux Behavior of Blazars

abstract: The variability of blazars in the X-ray and optical regions informs both the physics of their emitting region and places demands on the observer if a program requires that the object be bright or faint. The extensive simultaneous X-ray and optical observation by Swift provides the best insight into the variable nature of these objects. This program uses \textit{Swift} data for 19 X-ray-bright blazars, generally at $z > 0.1$, to determine their variability properties. The analysis is based on structure functions and provides insight into the nature of the variability and how it depends on time, luminosity, and redshift. We also consider strategies for observing blazars at or above average brightness, given a time delay between planning an observation and obtaining the data. This is critical to observations with orbiting X-ray telescopes, current or future. The variability in the soft X-ray band is typically three to eight times larger than at UV-optical wavelengths, at fixed time differences (i.e., 30 or 100 days). There is almost no difference in the amplitude of variation (X-ray or UV-optical) as a function of redshift (time delay of 30 days) and a modest positive correlation with luminosity. In the X-ray band, blazars that become brighter than normal typically remain bright for at least 2-3 months, although with significant flickering. One can avoid observing objects that are significantly below the average X-ray flux by scheduling the observation when the $F_X > 0.9 F_{X,avg}$, which requires monitoring observations near the time of the scheduling activity.


title: 370 New Eclipsing Binary Candidates from TESS Sectors 1-26

abstract: We present 370 candidate eclipsing binaries (EBs), identified from ~510,000 short cadence TESS light curves. Our statistical criteria identify 5,105 light curves with features consistent with eclipses (~1% of the initial sample). After visual confirmation of the light curves, we have a final sample of 2,288 EB candidates. Among these, we find 370 sources that were not included in the catalog recently published by Prsa et al. We publish our full sample of 370 new EB candidates, and statistical features used for their identification, reported per observation sector.


title: The Fay relations satisfied by the elliptic associator

abstract: Let $A_\tau$ denote the elliptic associator constructed by Enriquez, a power series in two non-commutative variables $a,b$ defined as an iterated integral of the Kronecker function $F_\tau$. We study a family of {\it Fay relations} satisfied by $A_\tau$, derived from the original Fay relation satisfied by the $F_\tau$. The Fay relations of $A_\tau$ were studied by Broedel, Matthes and Schlotterer, and determined up to non-explicit correction terms that arise from the necessity of regularizing the non-convergent integral. Here we study a reduced version $\bar{A}_\tau$ mod $2\pi i$. We recall a different construction of $\bar{A}_\tau$ in three steps, due to Matthes, Lochak and the author: first one defines the reduced {\it elliptic generating series} $\bar{E}_\tau$ which comes from the reduced Drinfeld associator $\overline{\Phi}_{KZ}$ and whose coefficients generate the same ring $\bar{R}$ as those of $\bar{A}_\tau$; then one defines $\Psi$ to be the automorphism of the free associative ring $\bar{R}\langle\langle a,b\rangle\rangle$ defined by $\Psi(a)=\bar{E}_\tau$ and $\Psi([a,b])=[a,b]$; finally one shows that the reduced elliptic associator $\bar{A}_\tau$ is equal to $\Psi\bigl({{ad(b)}\over{e^{ad(b)}-1}}(a)\bigr)$. Using this construction and mould theory and working with Lie-like versions of the elliptic generating series and associator, we prove the following results: (1) a mould satisfies the Fay relations if and only if a closely related mould satisfies the "swap circ-neutrality" relations defining the elliptic Kashiwara-Vergne Lie algebra $krv_{ell}$, (2) the reduced elliptic generating series satisfies a family of Fay relations with extremely simple correction terms coming directly from those of the Drinfeld associator, and (3) the correction terms for the Fay relations satisfied by the reduced elliptic associator can be deduced explicitly from these.


title: The Effect of Increased Access to IVF on Women's Careers

abstract: Motherhood is the main contributor to gender gaps in the labor market. IVF is a method of assisted reproduction that can delay fertility, which results in decreased motherhood income penalty. In this research, I estimate the effects of expanded access to in vitro fertilization (IVF) arising from state insurance mandates. I use a difference-in-differences model to estimate the effect of increased IVF accessibility for delaying childbirth and decreasing the motherhood income penalty. Using the fertility supplement dataset from the Current Population Survey (CPS), I estimate how outcomes change in states when they implement their mandates compared to how outcomes change in states that are not changing their policies. The results indicate that IVF mandates increase the probability of motherhood by 38 by 3.1 percentage points (p<0.01). However, the results provide no evidence that IVF insurance mandates impact women's earnings.


title: Solutions in Nonlinear Electrodynamics and their double copy regular black holes

abstract: We study solutions in non-linear electrodynamics (NED) and establish several general results. We show, that the $SO(2)$ electric-magnetic duality symmetry is restrictive enough to allow for reconstruction of the NED Lagrangian from the spherically-symmetric electrostatic (Coulomb-like) solution -- although there are infinitely many different NED theories admitting a given solution, there exists a unique $SO(2)$ invariant one among them. We introduce a general algorithm for constructing new $SO(2)$ invariant NED theories in the conventional approach, where only a few examples are available. We also show how to derive the Lagrangian of the $SO(2)$ invariant theory admitting a given electrostatic solution. We further show on a simple example that some NED theories may require sources (particles) of finite (non-zero) size. Such a non-zero size source not only regularizes the infinite energy of the point charge but also satisfies the condition of regularity, that the electric field is zero at the origin. The latter condition was identified earlier as necessary and sufficient for the NED solution to generate a regular black hole via so-called double copy construction and is also satisfied by solitons. We propose a large class of solitonic NED solutions that give rise to regular black holes via double copy construction and contain solutions of Maxwell and Born-Infeld as different limits. This class of NED solutions acquires two new properties in the limit where the corresponding regular black hole's asymptotics becomes Minkowski: it gives rise to regular higher-spin black holes via generalization of double copy -- ``higher-copy'' construction, and for very short distances changes the sign of the force becoming repulsive/attractive for opposite/similar signs of charges.


title: Introducing k4.0s: a Model for Mixed-Criticality Container Orchestration in Industry 4.0 (extended)

abstract: Time predictable edge cloud is seen as the answer for many arising needs in Industry 4.0 environments, since it is able to provide flexible, modular, and reconfigurable services with low latency and reduced costs. Orchestration systems are becoming the core component of clouds since they take decisions on the placement and lifecycle of software components. Current solutions start introducing real-time containers support for time predictability; however, these approaches lack of determinism as well as support for workloads requiring multiple levels of assurance/criticality. In this paper, we present k4.0s, an orchestration model for real-time and mixed-criticality environments, which includes timeliness, criticality and network requirements. The model leverages new abstractions for both node and jobs, e.g., node assurance, and requires novel monitoring strategies. We sketch an implementation of the proposal based on Kubernetes, and present an experimentation motivating the need for node assurance levels and adequate monitoring.


title: Optimizing Objective Functions from Trained ReLU Neural Networks via Sampling

abstract: This paper introduces scalable, sampling-based algorithms that optimize trained neural networks with ReLU activations. We first propose an iterative algorithm that takes advantage of the piecewise linear structure of ReLU neural networks and reduces the initial mixed-integer optimization problem (MIP) into multiple easy-to-solve linear optimization problems (LPs) through sampling. Subsequently, we extend this approach by searching around the neighborhood of the LP solution computed at each iteration. This scheme allows us to devise a second, enhanced algorithm that reduces the initial MIP problem into smaller, easier-to-solve MIPs. We analytically show the convergence of the methods and we provide a sample complexity guarantee. We also validate the performance of our algorithms by comparing them against state-of-the-art MIP-based methods. Finally, we show computationally how the sampling algorithms can be used effectively to warm-start MIP-based methods.


title: Intrinsically harmonic forms and characterization of flat circle bundles

abstract: We establish a criterion for the flatness of a principal circle bundle in terms of the intrinsically harmonic form problem. It states that the flatness is equivalent to the intrinsic harmonicity of a certain natural associated form.


title: Sensing Eating Events in Context: A Smartphone-Only Approach

abstract: While the task of automatically detecting eating events has been examined in prior work using various wearable devices, the use of smartphones as standalone devices to infer eating events remains an open issue. This paper proposes a framework that infers eating vs. non-eating events from passive smartphone sensing and evaluates it on a dataset of 58 college students. First, we show that time of the day and features from modalities such as screen usage, accelerometer, app usage, and location are indicative of eating and non-eating events. Then, we show that eating events can be inferred with an AUROC (area under the receiver operating characteristics curve) of 0.65 using subject-independent machine learning models, which can be further improved up to 0.81 for subject-dependent and 0.81 for hybrid models using personalization techniques. Moreover, we show that users have different behavioral and contextual routines around eating episodes requiring specific feature groups to train fully personalized models. These findings are of potential value for future mobile food diary apps that are context-aware by enabling scalable sensing-based eating studies using only smartphones; detecting under-reported eating events, thus increasing data quality in self report-based studies; providing functionality to track food consumption and generate reminders for on-time collection of food diaries; and supporting mobile interventions towards healthy eating practices.


title: Constrained Langevin Algorithms with L-mixing External Random Variables

abstract: Langevin algorithms are gradient descent methods augmented with additive noise, and are widely used in Markov Chain Monte Carlo (MCMC) sampling, optimization, and machine learning. In recent years, the non-asymptotic analysis of Langevin algorithms for non-convex learning has been extensively explored. For constrained problems with non-convex losses over a compact convex domain with IID data variables, the projected Langevin algorithm achieves a deviation of $O(T^{-1/4} (\log T)^{1/2})$ from its target distribution [27] in $1$-Wasserstein distance. In this paper, we obtain a deviation of $O(T^{-1/2} \log T)$ in $1$-Wasserstein distance for non-convex losses with $L$-mixing data variables and polyhedral constraints (which are not necessarily bounded). This improves on the previous bound for constrained problems and matches the best-known bound for unconstrained problems.


title: Active metasurfaces: lighting the path to commercial success

abstract: Active optical metasurfaces are rapidly emerging as a major frontier in photonics research, development, and commercialization. They promise compact, light-weight, and energy-efficient reconfigurable optical systems with unprecedented performance and functions that can be dynamically defined on-demand. Compared to their passive counterparts, the reconfiguration capacity of active metasurfaces also set additional challenges in scalable design, manufacturing, and control toward their practical deployment. This perspective aims to review the state-of-the-art of active metasurface technologies and their applications while highlighting key research advances essential to enabling their transition from laboratory curiosity to commercial reality.


title: On Sums of hyper-Kloosterman Sums

abstract: A formula of Kuznetsov allows one to interpret a smooth sum of Kloosterman sums as a sum over the spectrum of $GL(2)$ automorphic forms. In this paper, we construct a similar formula for the first hyper-Kloosterman sums using $GL(3)$ automorphic forms, resolving a long-standing problem of Bump, Friedberg and Goldfeld. Along the way, we develop what are apparently new bounds for the order derivatives of the classical $J$-Bessel function, and we conclude with a discussion of the original method of Bump, Friedberg and Goldfeld.


title: Unsupervised learning of features and object boundaries from local prediction

abstract: A visual system has to learn both which features to extract from images and how to group locations into (proto-)objects. Those two aspects are usually dealt with separately, although predictability is discussed as a cue for both. To incorporate features and boundaries into the same model, we model a layer of feature maps with a pairwise Markov random field model in which each factor is paired with an additional binary variable, which switches the factor on or off. Using one of two contrastive learning objectives, we can learn both the features and the parameters of the Markov random field factors from images without further supervision signals. The features learned by shallow neural networks based on this loss are local averages, opponent colors, and Gabor-like stripe patterns. Furthermore, we can infer connectivity between locations by inferring the switch variables. Contours inferred from this connectivity perform quite well on the Berkeley segmentation database (BSDS500) without any training on contours. Thus, computing predictions across space aids both segmentation and feature learning, and models trained to optimize these predictions show similarities to the human visual system. We speculate that retinotopic visual cortex might implement such predictions over space through lateral connections.


title: FadMan: Federated Anomaly Detection across Multiple Attributed Networks

abstract: Anomaly subgraph detection has been widely used in various applications, ranging from cyber attack in computer networks to malicious activities in social networks. Despite an increasing need for federated anomaly detection across multiple attributed networks, only a limited number of approaches are available for this problem. Federated anomaly detection faces two major challenges. One is that isolated data in most industries are restricted share with others for data privacy and security. The other is most of the centralized approaches training based on data integration. The main idea of federated anomaly detection is aligning private anomalies from local data owners on the public anomalies from the attributed network in the server through public anomalies to federate local anomalies. In each private attributed network, the detected anomaly subgraph is aligned with an anomaly subgraph in the public attributed network. The significant public anomaly subgraphs are selected for federated private anomalies while preventing local private data leakage. The proposed algorithm FadMan is a vertical federated learning framework for public node aligned with many private nodes of different features, and is validated on two tasks correlated anomaly detection on multiple attributed networks and anomaly detection on an attributeless network using five real-world datasets. In the first scenario, FadMan outperforms competitive methods by at least 12% accuracy at 10% noise level. In the second scenario, by analyzing the distribution of abnormal nodes, we find that the nodes of traffic anomalies are associated with the event of postgraduate entrance examination on the same day.


title: Exposing Gravitational Waves below the Quantum Shot Noise

abstract: The sensitivities of ground-based gravitational-wave (GW) detectors are limited by quantum shot noise at a few hundred Hertz and above. Nonetheless, one can use a quantum-correlation technique proposed by Martynov, et al. [Phys. Rev. A 95, 043831 (2017)] to remove the expectation value of the shot noise, thereby exposing underlying classical signals in the cross spectrum formed by cross-correlating the two outputs in a GW interferometer's anti-symmetric port. We explore here the prospects and analyze the sensitivity of using quantum correlation to detect astrophysical GW signals. Conceptually, this technique is similar to the correlation of two different GW detectors as it utilizes the fact that a GW signal will be correlated in the two outputs but the shot noise will be uncorrelated. Quantum correlation also has its unique advantages as it requires only a single interferometer to make a detection. Therefore, quantum correlation could increase the duty cycle, enhance the search efficiency, and enable the detection of highly polarized signals. In particular, we show that quantum correlation could be especially useful for detecting post-merger remnants of binary neutron stars with both short ($< 1\,{\rm s}$) and intermediate ($\sim 10-10^4\,{\rm s}$) durations and setting upper limits on continuous emissions from unknown pulsars.


title: Generalized Reductions: Making any Hierarchical Clustering Fair and Balanced with Low Cost

abstract: Clustering is a fundamental building block of modern statistical analysis pipelines. Fair clustering has seen much attention from the machine learning community in recent years. We are some of the first to study fairness in the context of hierarchical clustering, after the results of Ahmadian et al. from NeurIPS in 2020. We evaluate our results using Dasgupta's cost function, perhaps one of the most prevalent theoretical metrics for hierarchical clustering evaluation. Our work vastly improves the previous $O(n^{5/6}poly\log(n))$ fair approximation for cost to a near polylogarithmic $O(n^\delta poly\log(n))$ fair approximation for any constant $\delta\in(0,1)$. This result establishes a cost-fairness tradeoff and extends to broader fairness constraints than the previous work. We also show how to alter existing hierarchical clusterings to guarantee fairness and cluster balance across any level in the hierarchy.


title: AGILE Observations of GRB 220101A: A "New Year's Burst" with an Exceptionally Huge Energy Release

abstract: We report the AGILE observations of GRB 220101A, which took place at the beginning of 1st January 2022 and was recognized as one of the most energetic gamma-ray bursts (GRBs) ever detected since their discovery. The AGILE satellite acquired interesting data concerning the prompt phase of this burst, providing an overall temporal and spectral description of the event in a wide energy range, from tens of keV to tens of MeV. Dividing the prompt emission into three main intervals, we notice an interesting spectral evolution, featuring a notable hardening of the spectrum in the central part of the burst. The average fluxes encountered in the different time intervals are relatively moderate, with respect to those of other remarkable bursts, and the overall fluence exhibits a quite ordinary value among the GRBs detected by MCAL. However, GRB 220101A is the second farthest event detected by AGILE, and the burst with the highest isotropic equivalent energy of the whole MCAL GRB sample, releasing E_iso=2.54x10^54 erg and exhibiting an isotropic luminosity of L_iso=2.34x10^52 erg/s (both in the 400 keV - 10 MeV energy range). We also analyzed the first 10^6 s of the afterglow phase, using the publicly available Swift XRT data, carrying out a theoretical analysis of the afterglow, based on the forward shock model. We notice that GRB 220101A is with high probability surrounded with a wind-like density medium, and that the energy carried by the initial shock shall be a fraction of the total E_iso, presumably near 50%.


title: Energy dynamics, heat production and heat-work conversion with qubits: towards the development of quantum machines

abstract: We present an overview of recent advances in the study of energy dynamics and mechanisms for energy conversion in qubit systems with special focus on realizations in superconducting quantum circuits. We briefly introduce the relevant theoretical framework to analyze heat generation, energy transport and energy conversion in these systems with and without time-dependent driving considering the effect of equilibrium and non-equilibrium environments. We analyze specific problems and mechanisms under current investigation in the context of qubit systems. These include the problem of energy dissipation and possible routes for its control, energy pumping between driving sources and heat pumping between reservoirs, implementation of thermal machines and mechanisms for energy storage. We highlight the underlying fundamental phenomena related to geometrical and topological properties, as well as many-body correlations. We also present an overview of recent experimental activity in this field.


title: Collisional S-Matrix for the Vibrational Dynamics of H+H2 by Quantum Computing

abstract: An algorithm and a system of quantum circuits is developed and applied to compute accurately the S matrix for the transitions between vibrational states of H2 for collisions with H. The algorithm was applied to 100 eV laboratory collision energy at a quantum circuit simulator. The effects of the discretized dissociative continuum to the transition cross sections are carefully studied and accuracy and convergence of the results with the chosen parameters of the algorithm and the collision system are verified by comparison with a solution of the time-dependent Schrodinger equation using the classical algorithm as well as comparison with a few results from the literature.


title: Robust Phi-Divergence MDPs

abstract: In recent years, robust Markov decision processes (MDPs) have emerged as a prominent modeling framework for dynamic decision problems affected by uncertainty. In contrast to classical MDPs, which only account for stochasticity by modeling the dynamics through a stochastic process with a known transition kernel, robust MDPs additionally account for ambiguity by optimizing in view of the most adverse transition kernel from a prescribed ambiguity set. In this paper, we develop a novel solution framework for robust MDPs with s-rectangular ambiguity sets that decomposes the problem into a sequence of robust Bellman updates and simplex projections. Exploiting the rich structure present in the simplex projections corresponding to phi-divergence ambiguity sets, we show that the associated s-rectangular robust MDPs can be solved substantially faster than with state-of-the-art commercial solvers as well as a recent first-order solution scheme, thus rendering them attractive alternatives to classical MDPs in practical applications.


title: Dynamic bonding influenced by the proximity of adatoms to one-atom high step edges

abstract: Low-temperature scanning tunneling microscopy is used here to study dynamic bonding of gold atoms on surfaces under low coordination conditions. In the experiments, using an atomically-sharp gold tip, a gold adatom is deposited onto a gold surface with atomic precision either on the first hollow site near a step edge, or far away from it. Classical molecular dynamics simulations at 4.2 K and density functional theory calculations serve to elucidate the difference in the bonding behavior between these two different placements, while also providing information on the crystalline classification of the STM tips based on their experimental performance.


title: Multimodal Masked Autoencoders Learn Transferable Representations

abstract: Building scalable models to learn from diverse, multimodal data remains an open challenge. For vision-language data, the dominant approaches are based on contrastive learning objectives that train a separate encoder for each modality. While effective, contrastive learning approaches introduce sampling bias depending on the data augmentations used, which can degrade performance on downstream tasks. Moreover, these methods are limited to paired image-text data, and cannot leverage widely-available unpaired data. In this paper, we investigate whether a large multimodal model trained purely via masked token prediction, without using modality-specific encoders or contrastive learning, can learn transferable representations for downstream tasks. We propose a simple and scalable network architecture, the Multimodal Masked Autoencoder (M3AE), which learns a unified encoder for both vision and language data via masked token prediction. We provide an empirical study of M3AE trained on a large-scale image-text dataset, and find that M3AE is able to learn generalizable representations that transfer well to downstream tasks. Surprisingly, we find that M3AE benefits from a higher text mask ratio (50-90%), in contrast to BERT whose standard masking ratio is 15%, due to the joint training of two data modalities. We also provide qualitative analysis showing that the learned representation incorporates meaningful information from both image and language. Lastly, we demonstrate the scalability of M3AE with larger model size and training time, and its flexibility to train on both paired image-text data as well as unpaired data.


title: ALMA: Hierarchical Learning for Composite Multi-Agent Tasks

abstract: Despite significant progress on multi-agent reinforcement learning (MARL) in recent years, coordination in complex domains remains a challenge. Work in MARL often focuses on solving tasks where agents interact with all other agents and entities in the environment; however, we observe that real-world tasks are often composed of several isolated instances of local agent interactions (subtasks), and each agent can meaningfully focus on one subtask to the exclusion of all else in the environment. In these composite tasks, successful policies can often be decomposed into two levels of decision-making: agents are allocated to specific subtasks and each agent acts productively towards their assigned subtask alone. This decomposed decision making provides a strong structural inductive bias, significantly reduces agent observation spaces, and encourages subtask-specific policies to be reused and composed during training, as opposed to treating each new composition of subtasks as unique. We introduce ALMA, a general learning method for taking advantage of these structured tasks. ALMA simultaneously learns a high-level subtask allocation policy and low-level agent policies. We demonstrate that ALMA learns sophisticated coordination behavior in a number of challenging environments, outperforming strong baselines. ALMA's modularity also enables it to better generalize to new environment configurations. Finally, we find that while ALMA can integrate separately trained allocation and action policies, the best performance is obtained only by training all components jointly. Our code is available at https://github.com/shariqiqbal2810/ALMA


title: Network Digital Twin: Context, Enabling Technologies and Opportunities

abstract: The proliferation of emergent network applications (e.g., telesurgery, metaverse) is increasing the difficulty of managing modern communication networks. These applications entail stringent network requirements (e.g., ultra-low deterministic latency), which hinders network operators to manage their resources efficiently. In this article, we introduce the network digital twin (NDT), a renovated concept of classical network modeling tools whose goal is to build accurate data-driven network models that can operate in real-time. We describe the general architecture of the NDT and argue that modern machine learning (ML) technologies enable building some of its core components. Then, we present a case study that leverages a ML-based NDT for network performance evaluation and apply it to routing optimization in a QoS-aware use case. Lastly, we describe some key open challenges and research opportunities yet to be explored to achieve effective deployment of NDTs in real-world networks.


title: A comparison of Fourier and POD mode decomposition methods for high-speed Hall thruster video

abstract: Hall thrusters are susceptible to large-amplitude plasma oscillations that impact thruster performance and lifetime and are also difficult to model. High-speed cameras are a popular tool to study these dynamics due to their spatial resolution and are a popular, nonintrusive complement to in-situ probes. High-speed video of thruster oscillations can be isolated (decomposed) into coherent structures (modes) with algorithms that help us better understand the evolution and interactions of each. This work provides an introduction, comparison, and step-by-step tutorial on established Fourier and newer Proper Orthogonal Decomposition (POD) algorithms as applied to high-speed video of the unshielded H6 6-kW laboratory model Hall thruster. From this dataset, both sets of algorithms identify and characterize $m=0$ and $m>0$ modes in the discharge channel and cathode regions of the thruster plume, as well as mode hopping between the $m=3$ and $m=4$ rotating spokes in the channel. The Fourier methods are ideal for characterizing linear modal structures and also provide intuitive dispersion relationships. By contrast, the POD method tailors a basis set using energy minimization techniques that better captures the nonlinear nature of these structures and with a simpler implementation. Together, the Fourier and POD methods provide a more complete toolkit for studying Hall thruster plasma instabilities and mode dynamics. Specifically, we recommend first applying POD first to quickly identify the nature and location of global dynamics and then using Fourier methods to isolate dispersion plots and other wave-based physics.


title: Targeted Adaptive Design

abstract: Modern advanced manufacturing and advanced materials design often require searches of relatively high-dimensional process control parameter spaces for settings that result in optimal structure, property, and performance parameters. The mapping from the former to the latter must be determined from noisy experiments or from expensive simulations. We abstract this problem to a mathematical framework in which an unknown function from a control space to a design space must be ascertained by means of expensive noisy measurements, which locate optimal control settings generating desired design features within specified tolerances, with quantified uncertainty. We describe targeted adaptive design (TAD), a new algorithm that performs this sampling task efficiently. TAD creates a Gaussian process surrogate model of the unknown mapping at each iterative stage, proposing a new batch of control settings to sample experimentally and optimizing the updated log-predictive likelihood of the target design. TAD either stops upon locating a solution with uncertainties that fit inside the tolerance box or uses a measure of expected future information to determine that the search space has been exhausted with no solution. TAD thus embodies the exploration-exploitation tension in a manner that recalls, but is essentially different from, Bayesian optimization and optimal experimental design.


title: StarGraph: Knowledge Representation Learning based on Incomplete Two-hop Subgraph

abstract: Conventional representation learning algorithms for knowledge graphs (KG) map each entity to a unique embedding vector, ignoring the rich information contained in the neighborhood. We propose a method named StarGraph, which gives a novel way to utilize the neighborhood information for large-scale knowledge graphs to obtain entity representations. An incomplete two-hop neighborhood subgraph for each target node is at first generated, then processed by a modified self-attention network to obtain the entity representation, which is used to replace the entity embedding in conventional methods. We achieved SOTA performance on ogbl-wikikg2 and got competitive results on fb15k-237. The experimental results proves that StarGraph is efficient in parameters, and the improvement made on ogbl-wikikg2 demonstrates its great effectiveness of representation learning on large-scale knowledge graphs. The code is now available at \url{https://github.com/hzli-ucas/StarGraph}.


title: MIP-GNN: A Data-Driven Framework for Guiding Combinatorial Solvers

abstract: Mixed-integer programming (MIP) technology offers a generic way of formulating and solving combinatorial optimization problems. While generally reliable, state-of-the-art MIP solvers base many crucial decisions on hand-crafted heuristics, largely ignoring common patterns within a given instance distribution of the problem of interest. Here, we propose MIP-GNN, a general framework for enhancing such solvers with data-driven insights. By encoding the variable-constraint interactions of a given mixed-integer linear program (MILP) as a bipartite graph, we leverage state-of-the-art graph neural network architectures to predict variable biases, i.e., component-wise averages of (near) optimal solutions, indicating how likely a variable will be set to 0 or 1 in (near) optimal solutions of binary MILPs. In turn, the predicted biases stemming from a single, once-trained model are used to guide the solver, replacing heuristic components. We integrate MIP-GNN into a state-of-the-art MIP solver, applying it to tasks such as node selection and warm-starting, showing significant improvements compared to the default setting of the solver on two classes of challenging binary MILPs.


title: KL-Entropy-Regularized RL with a Generative Model is Minimax Optimal

abstract: In this work, we consider and analyze the sample complexity of model-free reinforcement learning with a generative model. Particularly, we analyze mirror descent value iteration (MDVI) by Geist et al. (2019) and Vieillard et al. (2020a), which uses the Kullback-Leibler divergence and entropy regularization in its value and policy updates. Our analysis shows that it is nearly minimax-optimal for finding an $\varepsilon$-optimal policy when $\varepsilon$ is sufficiently small. This is the first theoretical result that demonstrates that a simple model-free algorithm without variance-reduction can be nearly minimax-optimal under the considered setting.


title: Exemplar Free Class Agnostic Counting

abstract: We tackle the task of Class Agnostic Counting, which aims to count objects in a novel object category at test time without any access to labeled training data for that category. All previous class agnostic counting methods cannot work in a fully automated setting, and require computationally expensive test time adaptation. To address these challenges, we propose a visual counter which operates in a fully automated setting and does not require any test time adaptation. Our proposed approach first identifies exemplars from repeating objects in an image, and then counts the repeating objects. We propose a novel region proposal network for identifying the exemplars. After identifying the exemplars, we obtain the corresponding count by using a density estimation based Visual Counter. We evaluate our proposed approach on FSC-147 dataset, and show that it achieves superior performance compared to the existing approaches.


title: When Is Recoverable Consensus Harder Than Consensus?

abstract: We study the ability of different shared object types to solve recoverable consensus using non-volatile shared memory in a system with crashes and recoveries. In particular, we compare the difficulty of solving recoverable consensus to the difficulty of solving the standard wait-free consensus problem in a system with halting failures. We focus on the model where individual processes may crash and recover and the large class of object types that are equipped with a read operation. We characterize the readable object types that can solve recoverable consensus among a given number of processes. Using this characterization, we show that the number of processes that can solve consensus using a readable type can be larger than the number of processes that can solve recoverable consensus using that type, but only slightly larger.


title: Bias in the representative volume element method: periodize the ensemble instead of its realizations

abstract: We study the Representative Volume Element (RVE) method, which is a method to approximately infer the effective behavior $a_{\text{hom}}$ of a stationary random medium. The latter is described by a coefficient field $a(x)$ generated from a given ensemble $\langle\cdot\rangle$ and the corresponding linear elliptic operator $-\nabla\cdot a\nabla$. In line with the theory of homogenization, the method proceeds by computing $d = 3$ correctors (d denoting the space dimension).To be numerically tractable, this computation has to be done on a finite domain: the so-called "representative" volume element, i. e. a large box with, say, periodic boundary conditions. The main message of this article is: Periodize the ensemble instead of its realizations. By this we mean that it is better to sample from a suitably periodized ensemble than to periodically extend the restriction of a realization $a(x)$ from the whole-space ensemble $\langle\cdot\rangle$. We make this point by investigating the bias (or systematic error), i. e. the difference between $a_{\text{hom}}$ and the expected value of the RVE method, in terms of its scaling w. r. t. the lateral size $L$ of the box. In case of periodizing $a(x)$, we heuristically argue that this error is generically $O(L^{-1})$. In case of a suitable periodization of $\langle\cdot\rangle$, we rigorously show that it is $O(L^{-d})$. In fact, we give a characterization of the leading-order error term for both strategies, and argue that even in the isotropic case it is generically non-degenerate. We carry out the rigorous analysis in the convenient setting of ensembles $\langle\cdot\rangle$ of Gaussian type with integrable covariance, which allow for a straightforward periodization and which make the Price theorem and the Malliavin calculus available for optimal stochastic estimates of correctors.


title: Taylor microscale and effective Reynolds number near the Sun from PSP

abstract: The Taylor microscale is a fundamental length scale in turbulent fluids, representing the end of fluid properties and onset of dissipative processes. The Taylor microscale can also be used to evaluate the Reynolds number in classical turbulence theory. Although the solar wind is weakly collisional, it approximately behaves as a magnetohydrodynamic (MHD) fluid at scales larger than the kinetic scale. As a result, classical fluid turbulence theory and formalisms are often used to study turbulence in the MHD range. Therefore, a Taylor microscale can be used to estimate an effective Reynolds number in the solar wind. NASA's Parker Solar Probe (PSP) has reached progressively closer to the Sun than any other spacecraft before. The collected data have revealed many new findings in the near-Sun solar wind. Here, we use the PSP data to estimate the Taylor microscale and effective Reynolds number near the Sun. We find that the Taylor microscale and Reynolds number are small compared to the corresponding near-Earth values, indicating a solar wind that has been less processed by turbulence, with very small-scale dissipative processes near the Sun.


title: NuSTAR Observations of AGN with Low Observed X-ray to [OIII] Luminosity Ratios: Heavily Obscured AGN or Turned-Off AGN?

abstract: Type 2 active galactic nuclei (AGN) show signatures of accretion onto a supermassive black hole through strong, high-ionization, narrow emission lines extended on scales of 100s to 1000s of parsecs, but they lack the broad emission lines from close in to the black hole that characterize type 1 AGN. The lack of broad emission could indicate obscuration of the innermost nuclear regions, or could indicate that the black hole is no longer strongly accreting. Since high-energy X-rays can penetrate thick obscuring columns, they have the power to distinguish these two scenarios. We present high-energy NuSTAR observations of 9 Seyfert 2 AGN from the IRAS 12 micron survey, supplemented with low-energy X-ray observations from Chandra, XMM-Newton, and Swift. The galaxies were selected to have anomalously low observed 2-10 keV luminosities compared to their [O III] optical luminosities, a traditional diagnostic of heavily obscured AGN, reaching into the Compton-thick regime for the highest hydrogen column densities ($N_{\rm H} > 1.5 \times 10^{24}\, {\rm cm}^{-2}$). Based on updated [O III] luminosities and intrinsic X-ray luminosities based on physical modeling of the hard X-ray spectra, we find that one galaxy was misclassified as type 2 (NGC 5005) and most of the remaining AGN are obscured, including three confirmed as Compton-thick (IC 3639, NGC 1386, and NGC 3982). One galaxy, NGC 3627, appears to be recently deactivated. Compared to the original sample the 9 AGN were selected from, this is a rate of approximately 1%. We also find a new X-ray changing-look AGN in NGC 6890.


title: Diffusion-LM Improves Controllable Text Generation

abstract: Controlling the behavior of language models (LMs) without re-training is a major open problem in natural language generation. While recent works have demonstrated successes on controlling simple sentence attributes (e.g., sentiment), there has been little progress on complex, fine-grained controls (e.g., syntactic structure). To address this challenge, we develop a new non-autoregressive language model based on continuous diffusions that we call Diffusion-LM. Building upon the recent successes of diffusion models in continuous domains, Diffusion-LM iteratively denoises a sequence of Gaussian vectors into word vectors, yielding a sequence of intermediate latent variables. The continuous, hierarchical nature of these intermediate variables enables a simple gradient-based algorithm to perform complex, controllable generation tasks. We demonstrate successful control of Diffusion-LM for six challenging fine-grained control tasks, significantly outperforming prior work.


title: On the Hunt for the Origins of the Orphan--Chenab Stream: Detailed Element Abundances with APOGEE and Gaia

abstract: Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan--Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the parent of the OC stream, we use astrometry from the early third data release of ESA's Gaia mission (Gaia EDR3) and radial velocity information from the SDSS-IV APOGEE survey to find up to 13 stars that are likely members of the OC stream. We use the APOGEE survey to study the chemical nature (for up to 13 stars) of the OC stream in the $\alpha$ (O, Mg, Ca, Si, Ti, S), odd-Z (Al, K, V), Fe-peak (Fe, Ni, Mn, Co, Cr) and neutron capture (Ce) elemental groups. We find that the stars that make up the OC stream are not consistent with a mono-metallic population and have a median metallicity of --1.92~dex with a dispersion of 0.28 dex. Our results also indicate that the $\alpha$-elements are depleted compared to the known Milky Way populations and that its [Mg/Al] abundance ratio is not consistent with second generation stars from globular clusters. The detailed chemical pattern of these stars indicates that the OC stream progenitor is very likely to be a dwarf spheroidal galaxy with a mass of ~10$^6$ M$_\odot$.


title: Controllable Text Generation with Neurally-Decomposed Oracle

abstract: We propose a general and efficient framework to control auto-regressive generation models with NeurAlly-Decomposed Oracle (NADO). Given a pre-trained base language model and a sequence-level boolean oracle function, we propose to decompose the oracle function into token-level guidance to steer the base model in text generation. Specifically, the token-level guidance is approximated by a neural model trained with examples sampled from the base model, demanding no additional auxiliary labeled data. Based on posterior regularization, we present the closed-form optimal solution to incorporate the token-level guidance into the base model for controllable generation. We further provide a theoretical analysis of how the approximation quality of NADO affects the controllable generation results. Experiments conducted on two applications: (1) text generation with lexical constraints and (2) machine translation with formality control demonstrate that our framework efficiently guides the base model towards the given oracle while maintaining high generation quality.


title: Selective Inference for Sparse Multitask Regression with Applications in Neuroimaging

abstract: Multi-task learning is frequently used to model a set of related response variables from the same set of features, improving predictive performance and modeling accuracy relative to methods that handle each response variable separately. Despite the potential of multi-task learning to yield more powerful inference than single-task alternatives, prior work in this area has largely omitted uncertainty quantification. Our focus in this paper is a common multi-task problem in neuroimaging, where the goal is to understand the relationship between multiple cognitive task scores (or other subject-level assessments) and brain connectome data collected from imaging. We propose a framework for selective inference to address this problem, with the flexibility to: (i) jointly identify the relevant covariates for each task through a sparsity-inducing penalty, and (ii) conduct valid inference in a model based on the estimated sparsity structure. Our framework offers a new conditional procedure for inference, based on a refinement of the selection event that yields a tractable selection-adjusted likelihood. This gives an approximate system of estimating equations for maximum likelihood inference, solvable via a single convex optimization problem, and enables us to efficiently form confidence intervals with approximately the correct coverage. Applied to both simulated data and data from the Adolescent Brain Cognitive Development (ABCD) study, our selective inference methods yield tighter confidence intervals than commonly used alternatives, such as data splitting. We also demonstrate through simulations that multi-task learning with selective inference can more accurately recover true signals than single-task methods.


title: Data-driven quantification of model-form uncertainty in Reynolds-averaged simulations of wind farms

abstract: Computational fluid dynamics using the Reynolds-averaged Navier-Stokes (RANS) remains the most cost-effective approach to study wake flows and power losses in wind farms. The underlying assumptions associated with turbulence closures are one of the biggest sources of errors and uncertainties in the model predictions. This work aims to quantify model-form uncertainties in RANS simulations of wind farms at high Reynolds numbers under neutrally stratified conditions by perturbing the Reynolds stress tensor through a data-driven machine-learning technique. To this end, a two-step feature-selection method is applied to determine key features of the model. Then, the extreme gradient boosting algorithm is validated and employed to predict the perturbation amount and direction of the modeled Reynolds stress toward the limiting states of turbulence on the barycentric map. This procedure leads to a more accurate representation of the Reynolds stress anisotropy. The data-driven model is trained on high-fidelity data obtained from large-eddy simulation of a specific wind farm, and it is tested on two other (unseen) wind farms with distinct layouts to analyze its performance in cases with different turbine spacing and partial wake. The results indicate that, unlike the data-free approach in which a uniform and constant perturbation amount is applied to the entire computational domain, the proposed framework yields an optimal estimation of the uncertainty bounds for the RANS-predicted quantities of interest, including the wake velocity, turbulence intensity, and power losses in wind farms.


title: Electrostatic correlation induced ion condensation and charge inversion in multivalent electrolytes

abstract: The study of the electrical double layer lies at the heart of colloidal and interfacial science. The standard mean-field Poisson-Boltzmann (PB) theory is incapable of modeling many phenomena originated from ion correlation. An important example is charge inversion or overcharging of electrical double layers in multivalent electrolyte solutions. Existing theories aiming to include correlations cannot capture the nonmonotonic dependence of charge inversion on salt concentration because they have not systematically accounted for the inhomogeneous nature of correlations from surface to the bulk and the excluded volume effect of ions and solvent molecules. In this work, we modify the Gaussian renormalized fluctuation theory by including the excluded volume effect to study ion condensation and charge inversion. A boundary layer approach is developed to accurately model the giant difference in ion correlations between the condensed layer near the surface and the diffuse layer outside. The theory is used to study charge inversion in multivalent electrolytes and their mixtures. We predict a surface charge induced formation of a three-dimensional condensed layer, which is necessary but not sufficient for charge inversion. The value of the effective surface potential is found to depend non-monotonically on the bulk salt concentration. Our results also show a non-monotonic reduction in charge inversion in monovalent and multivalent electrolyte mixtures. Our work is the first to qualitatively reproduce experimental and simulation observations and explains the underlying physics.


title: A short note on inf-sup conditions for the Taylor-Hood family $Q_k$-$Q_{k-1}$

abstract: We discuss two types of discrete inf-sup conditions for the Taylor-Hood family $Q_k$-$Q_{k-1}$ for all $k\in \mathbb{N}$ with $k\ge 2$ in 2D and 3D. While in 2D all results hold for a general class of hexahedral meshes, the results in 3D are restricted to meshes of parallelepipeds. The analysis is based on an element-wise technique as opposed to the widely used macroelement technique. This leads to inf-sup conditions on each element of the subdivision as well as to inf-sup conditions on the whole computational domain.


title: Will Bilevel Optimizers Benefit from Loops

abstract: Bilevel optimization has arisen as a powerful tool for solving a variety of machine learning problems. Two current popular bilevel optimizers AID-BiO and ITD-BiO naturally involve solving one or two sub-problems, and consequently, whether we solve these problems with loops (that take many iterations) or without loops (that take only a few iterations) can significantly affect the overall computational efficiency. Existing studies in the literature cover only some of those implementation choices, and the complexity bounds available are not refined enough to enable rigorous comparison among different implementations. In this paper, we first establish unified convergence analysis for both AID-BiO and ITD-BiO that are applicable to all implementation choices of loops. We then specialize our results to characterize the computational complexity for all implementations, which enable an explicit comparison among them. Our result indicates that for AID-BiO, the loop for estimating the optimal point of the inner function is beneficial for overall efficiency, although it causes higher complexity for each update step, and the loop for approximating the outer-level Hessian-inverse-vector product reduces the gradient complexity. For ITD-BiO, the two loops always coexist, and our convergence upper and lower bounds show that such loops are necessary to guarantee a vanishing convergence error, whereas the no-loop scheme suffers from an unavoidable non-vanishing convergence error. Our numerical experiments further corroborate our theoretical results.


title: Characterizing and mitigating coherent errors in a trapped ion quantum processor using hidden inverses

abstract: Quantum computing testbeds exhibit high-fidelity quantum control over small collections of qubits, enabling performance of precise, repeatable operations followed by measurements. Currently, these noisy intermediate-scale devices can support a sufficient number of sequential operations prior to decoherence such that near term algorithms can be performed with proximate accuracy (like chemical accuracy for quantum chemistry). While the results of these algorithms are imperfect, these imperfections can help bootstrap quantum computer testbed development. Demonstrations of these algorithms over the past few years, coupled with the idea that imperfect algorithm performance can be caused by several dominant noise sources in the quantum processor, which can be measured and calibrated during algorithm execution or in post-processing, has led to the use of noise mitigation to improve computational results. Conversely, benchmark algorithms coupled with noise mitigation can help diagnose the nature of noise, whether systematic or purely random. Here, we outline the use of coherent noise mitigation techniques as a characterization tool in trapped-ion testbeds. We perform model-fitting of the noisy data to determine the noise source based on realistic noise models and demonstrate that systematic noise amplification coupled with error mitigation schemes provides useful data for noise model deduction. Further, in order to connect lower level noise model details with application specific performance of near term algorithms, we experimentally construct the loss landscape of a variational algorithm under various injected noise sources coupled with error mitigation techniques. This type of connection enables application-aware hardware codesign, in which the most important noise sources in specific applications, like quantum chemistry, become foci of improvement in subsequent hardware generations.


title: Fast and Light-Weight Answer Text Retrieval in Dialogue Systems

abstract: Dialogue systems can benefit from being able to search through a corpus of text to find information relevant to user requests, especially when encountering a request for which no manually curated response is available. The state-of-the-art technology for neural dense retrieval or re-ranking involves deep learning models with hundreds of millions of parameters. However, it is difficult and expensive to get such models to operate at an industrial scale, especially for cloud services that often need to support a big number of individually customized dialogue systems, each with its own text corpus. We report our work on enabling advanced neural dense retrieval systems to operate effectively at scale on relatively inexpensive hardware. We compare with leading alternative industrial solutions and show that we can provide a solution that is effective, fast, and cost-efficient.


title: Dynamic Phases and Reentrant Hall Effect for Vortices and Skyrmions on Periodic Pinning Arrays

abstract: We consider the motion of superconducting vortices and skyrmions on a square substrate near the first commensurate matching field. Slightly above commensuration, a series of dynamic phases appear including interstitial flow, and there is a transition from fluid flow to soliton flow that generates negative differential conductivity. Slightly below commensuration, vacancy depinning occurs. The dynamic phase transitions produce features in the velocity-force curves, differential mobility, and velocity fluctuations. When a Magnus force is also present, as in certain superconducting vortex or skyrmion systems, there is an expansion of the fluid state, and at lower drives there is a finite Hall angle in the fluid phase but a vanishing Hall angle in the soliton phase, giving rise to a reentrant Hall effect. We also find a regime where the Hall motion of the particles exhibits the same dynamic phases, including soliton motion at a finite angle that produces negative differential conductivity in the Hall response but not in the longitudinal response. Our results suggest that vortices driven over periodic pinning may be an ideal system for determining if a vortex Hall effect is occurring and would also be relevant for skyrmions at smaller Magnus forces.


title: Sparse Conditional Hidden Markov Model for Weakly Supervised Named Entity Recognition

abstract: Weakly supervised named entity recognition methods train label models to aggregate the token annotations of multiple noisy labeling functions (LFs) without seeing any manually annotated labels. To work well, the label model needs to contextually identify and emphasize well-performed LFs while down-weighting the under-performers. However, evaluating the LFs is challenging due to the lack of ground truths. To address this issue, we propose the sparse conditional hidden Markov model (Sparse-CHMM). Instead of predicting the entire emission matrix as other HMM-based methods, Sparse-CHMM focuses on estimating its diagonal elements, which are considered as the reliability scores of the LFs. The sparse scores are then expanded to the full-fledged emission matrix with pre-defined expansion functions. We also augment the emission with weighted XOR scores, which track the probabilities of an LF observing incorrect entities. Sparse-CHMM is optimized through unsupervised learning with a three-stage training pipeline that reduces the training difficulty and prevents the model from falling into local optima. Compared with the baselines in the Wrench benchmark, Sparse-CHMM achieves a 3.01 average F1 score improvement on five comprehensive datasets. Experiments show that each component of Sparse-CHMM is effective, and the estimated LF reliabilities strongly correlate with true LF F1 scores.


title: Learning to Find Proofs and Theorems by Learning to Refine Search Strategies: The Case of Loop Invariant Synthesis

abstract: We propose a new approach to automated theorem proving where an AlphaZero-style agent is self-training to refine a generic high-level expert strategy expressed as a nondeterministic program. An analogous teacher agent is self-training to generate tasks of suitable relevance and difficulty for the learner. This allows leveraging minimal amounts of domain knowledge to tackle problems for which training data is unavailable or hard to synthesize. As a specific illustration, we consider loop invariant synthesis for imperative programs and use neural networks to refine both the teacher and solver strategies.


title: Semi-supervised Semantics-guided Adversarial Training for Trajectory Prediction

abstract: Predicting the trajectories of surrounding objects is a critical task for self-driving vehicles and many other autonomous systems. Recent works demonstrate that adversarial attacks on trajectory prediction, where small crafted perturbations are introduced to history trajectories, may significantly mislead the prediction of future trajectories and induce unsafe planning. However, few works have addressed enhancing the robustness of this important safety-critical task.In this paper, we present a novel adversarial training method for trajectory prediction. Compared with typical adversarial training on image tasks, our work is challenged by more random input with rich context and a lack of class labels. To address these challenges, we propose a method based on a semi-supervised adversarial autoencoder, which models disentangled semantic features with domain knowledge and provides additional latent labels for the adversarial training. Extensive experiments with different types of attacks demonstrate that our Semisupervised Semantics-guided Adversarial Training (SSAT) method can effectively mitigate the impact of adversarial attacks by up to 73% and outperform other popular defense methods. In addition, experiments show that our method can significantly improve the system's robust generalization to unseen patterns of attacks. We believe that such semantics-guided architecture and advancement on robust generalization is an important step for developing robust prediction models and enabling safe decision-making.


title: Adhesive Transfer operates during Galling

abstract: In order to reduce cobalt within the primary circuit of pressurised water reactors (PWRs), wear-resistant steels are being researched and developed. In particular interest is the understanding of galling mechanisms, an adhesive wear mechanism which is particularly prevalent in PWR valves. Here we show that large shear stresses and adhesive transfer occur during galling by exploiting the 2 wt per cent manganese difference between 304L and 316L stainless steels, even at relatively low compressive stresses of 50 MPa. Through these findings, the galling mechanisms of stainless steels can be better understood, which may help with the development of galling resistant stainless steels.


title: Competitive Gradient Optimization

abstract: We study the problem of convergence to a stationary point in zero-sum games. We propose competitive gradient optimization (CGO ), a gradient-based method that incorporates the interactions between the two players in zero-sum games for optimization updates. We provide continuous-time analysis of CGO and its convergence properties while showing that in the continuous limit, CGO predecessors degenerate to their gradient descent ascent (GDA) variants. We provide a rate of convergence to stationary points and further propose a generalized class of $\alpha$-coherent function for which we provide convergence analysis. We show that for strictly $\alpha$-coherent functions, our algorithm convergences to a saddle point. Moreover, we propose optimistic CGO (OCGO), an optimistic variant, for which we show convergence rate to saddle points in $\alpha$-coherent class of functions.


title: Absorbing the Arrow of Electromagnetic Radiation

abstract: We argue that the asymmetry between diverging and converging electromagnetic waves is just one of many asymmetries in observed phenomena that can be explained by a past hypothesis and statistical postulate (together assigning probabilities to different states of matter and field in the early universe). The arrow of electromagnetic radiation is thus absorbed into a broader account of temporal asymmetries in nature. We give an accessible introduction to the problem of explaining the arrow of radiation and compare our preferred strategy for explaining the arrow to three alternatives: (i) modifying the laws of electromagnetism by adding a radiation condition requiring that electromagnetic fields always be attributable to past sources, (ii) removing electromagnetic fields and having particles interact directly with one another through retarded action-at-a-distance, (iii) adopting the Wheeler-Feynman approach and having particles interact directly through half-retarded half-advanced action-at-a-distance. In addition to the asymmetry between diverging and converging waves, we also consider the related asymmetry of radiation reaction.


title: Two-Leg Deep Space Relay Architectures: Performance, Challenges, and Perspectives

abstract: In this paper, architectures for interplanetary communications that feature the use of a data relay are investigated. In the considered "two-leg" architecture, a spacecraft orbiting the Earth, or in orbit at a Lagrange point, receives data from a deep space probe (leg-1) and relays them towards ground (leg-2). Different wireless technologies for the interplanetary link, namely, radio frequencies above the Ka band and optical frequencies, are considered. Moreover, the cases of transparent and regenerative relaying as well as different different orbital configurations are addressed, offering a thorough analysis of such systems from different viewpoints. Results show that, under certain constraints in terms of pointing accuracy and onboard antenna size, the adoption of a two-leg architecture can achieve the data rates supported by direct space-to-Earth link configurations with remarkably smaller ground station antennas.


title: Freezing Sets for Arbitrary Digital Dimension

abstract: We show how to obtain freezing sets for digital images in $\mathbb{Z}^n$ for arbitrary $n$, using the $c_1$ and $c_n$ adjacencies.


title: FedControl: When Control Theory Meets Federated Learning

abstract: To date, the most popular federated learning algorithms use coordinate-wise averaging of the model parameters. We depart from this approach by differentiating client contributions according to the performance of local learning and its evolution. The technique is inspired from control theory and its classification performance is evaluated extensively in IID framework and compared with FedAvg.


title: Provably Sample-Efficient RL with Side Information about Latent Dynamics

abstract: We study reinforcement learning (RL) in settings where observations are high-dimensional, but where an RL agent has access to abstract knowledge about the structure of the state space, as is the case, for example, when a robot is tasked to go to a specific room in a building using observations from its own camera, while having access to the floor plan. We formalize this setting as transfer reinforcement learning from an abstract simulator, which we assume is deterministic (such as a simple model of moving around the floor plan), but which is only required to capture the target domain's latent-state dynamics approximately up to unknown (bounded) perturbations (to account for environment stochasticity). Crucially, we assume no prior knowledge about the structure of observations in the target domain except that they can be used to identify the latent states (but the decoding map is unknown). Under these assumptions, we present an algorithm, called TASID, that learns a robust policy in the target domain, with sample complexity that is polynomial in the horizon, and independent of the number of states, which is not possible without access to some prior knowledge. In synthetic experiments, we verify various properties of our algorithm and show that it empirically outperforms transfer RL algorithms that require access to "full simulators" (i.e., those that also simulate observations).


title: The branching number of intermediate growth trees

abstract: We introduce an "intermediate branching number"(IBN) which captures the branching of intermediate growth trees, similar in spirit to the well-studied branching number of exponential growth trees. We show that the IBN is the critical threshold for several random processes on trees, and analyze the IBN on some examples of interest. Our main result is an algorithm to find spherically symmetric trees with large IBN inside some permutation wreath products. We demonstrate the usefulness of these trees to the study of intermediate growth groups by using them to get the first tight bounds for the firefighter problem on some inetrmediate growth groups.


title: Global classical solutions to the compressible micropolar viscous fluids with large oscillations and vacuum

abstract: In this paper, we consider the three dimensional Cauchy problem of the compressible micropolar viscous flows, we prove the existence of unique global classical solution for smooth initial data with small initial energy but possibly large oscillations, the initial density may allowed to contain vacuum states. Furthermore, the large-time behavior of the solution is obtained.


title: Deterministic Langevin Monte Carlo with Normalizing Flows for Bayesian Inference

abstract: We propose a general purpose Bayesian inference algorithm for expensive likelihoods, replacing the stochastic term in the Langevin equation with a deterministic density gradient term. The particle density is evaluated from the current particle positions using a Normalizing Flow (NF), which is differentiable and has good generalization properties in high dimensions. We take advantage of NF preconditioning and NF based Metropolis-Hastings updates for a faster convergence. We show on various examples that the method is competitive against state of the art sampling methods.

categories: stat.ML cond-mat.stat-mech cs.LG physics.data-an stat.CO


title: Exact constraints and appropriate norms in machine learned exchange-correlation functionals

abstract: Machine learning techniques have received growing attention as an alternative strategy for developing general-purpose density functional approximations, augmenting the historically successful approach of human designed functionals derived to obey mathematical constraints known for the exact exchange-correlation functional. More recently efforts have been made to reconcile the two techniques, integrating machine learning and exact-constraint satisfaction. We continue this integrated approach, designing a deep neural network that exploits the exact constraint and appropriate norm philosophy to deorbitalize the strongly constrained and appropriately normed SCAN functional. The deep neural network is trained to replicate the SCAN functional from only electron density and local derivative information, avoiding use of the orbital dependent kinetic energy density. The performance and transferability of the machine learned functional are demonstrated for molecular and periodic systems.