title: Turbulent dissipation, CH$^+$ abundance, H$_2$ line luminosities, and polarization in the cold neutral medium

abstract: In the cold neutral medium, high out-of-equilibrium temperatures are created by intermittent dissipation processes, including shocks, viscous heating, and ambipolar diffusion. The high-temperature excursions are thought to explain the enhanced abundance of CH$^{+}$ observed along diffuse molecular sight-lines. Intermittent high temperatures should also have an impact on H$_2$ line luminosities. We carry out simulations of MHD turbulence in molecular clouds including heating and cooling, and post-process them to study H$_2$ line emission and hot-gas chemistry, particularly the formation of CH$^+$. We explore multiple magnetic field strengths and equations of state. We use a new H$_2$ cooling function for $n_{\rm H} \leq 10^5\,{\rm cm}^{-3}$, $T\leq 5000\,{\rm K}$, and variable H$_2$ fraction. We make two important simplifying assumptions: (i) the ${\rm H}_2/{\rm H}$ fraction is fixed everywhere, and (ii) we exclude from our analysis regions where the ion-neutral drift velocity is calculated to be greater than 5 km/s. Our models produce H$_2$ emission lines in accord with many observations, although extra excitation mechanisms are required in some clouds. For realistic r.m.s. magnetic field strengths ($\approx 10$ $\mu$G) and velocity dispersions, we reproduce observed CH$^+$ abundances. These findings contrast with those of Valdivia et al. (2017). Comparison of predicted dust polarization with observations by {\it Planck} suggests that the mean field $\gtrsim 5 \mu$G, so that the turbulence is sub-Alfv\'enic. We recommend future work treating ions and neutrals as separate fluids to more accurately capture the effects of ambipolar diffusion on CH$^+$ abundance.


title: Powderday: Dust Radiative Transfer for Galaxy Simulations

abstract: We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.


title: Leptoquarks in Oblique Corrections and Higgs Signal Strength: Status and Prospects

abstract: Leptoquarks (LQs) are predicted within Grand Unified Theories and are well motivated by the current flavor anomalies. In this article we investigate the impact of scalar LQs on Higgs decays and oblique corrections as complementary observables in the search for them. Taking into account all five LQ representations under the Standard Model gauge group and including the most general mixing among them, we calculate the effects in $h\to\gamma\gamma$, $h\to gg$, $h \to Z \gamma$ and the Peskin-Takeuchi parameters $S$, $T$ and $U$. We find that these observables depend on the same Lagrangian parameters, leading to interesting correlations among them. While the current experimental bounds only yield weak constraints on the model, these correlations can be used to distinguish different LQ representations at future colliders (ILC, CLIC, FCC-ee and FCC-hh), whose discovery potential we are going to discuss.


title: The Large Magellanic Cloud stellar content with SMASH: I. Assessing the stability of the Magellanic spiral arms

abstract: The Large Magellanic Cloud (LMC) is the closest and most studied example of an irregular galaxy. Among its principal defining morphological features, its off-centred bar and single spiral arm stand out, defining a whole family of galaxies known as the Magellanic spirals (Sm). These structures are thought to be triggered by tidal interactions and possibly maintained via gas accretion. However, it is still unknown whether they are long-lived stable structures. In this work, by combining photometry that reaches down to the oldest main sequence turn-off in the colour-magnitude diagrams (CMD, up to a distance of $\sim$4.4 kpc from the LMC centre) from the SMASH survey and CMD fitting techniques, we find compelling evidence supporting the long-term stability of the LMC spiral arm, dating the origin of this structure to more than 2~Gyr ago. The evidence suggests that the close encounter between the LMC and the Small Magellanic Cloud (SMC) that produced the gaseous Magellanic Stream and its Leading Arm (LA) also triggered the formation of the LMC's spiral arm. Given the mass difference between the Clouds and the notable consequences of this interaction, we can speculate that this should have been one of their closest encounters. These results set important constraints on the timing of LMC-SMC collisions, as well as on the physics behind star formation induced by tidal encounters.


title: Predictions of the Nancy Grace Roman Space Telescope Galactic Exoplanet Survey II: Free-Floating Planet Detection Rates

abstract: The Nancy Grace Roman Space Telescope (Roman) will perform a Galactic Exoplanet Survey (RGES) to discover bound exoplanets with semi-major axes greater than 1 au using gravitational microlensing. Roman will even be sensitive to planetary mass objects that are not gravitationally bound to any host star. Such free-floating planetary mass objects (FFPs) will be detected as isolated microlensing events with timescales shorter than a few days. A measurement of the abundance and mass function of FFPs is a powerful diagnostic of the formation and evolution of planetary systems, as well as the physics of the formation of isolated objects via direct collapse. We show that Roman will be sensitive to FFP lenses that have masses from that of Mars ($0.1 M_\oplus$) to gas giants ($M\gtrsim100M_\oplus$) as isolated lensing events with timescales from a few hours to several tens of days, respectively. We investigate the impact of the detection criteria on the survey, especially in the presence of finite-source effects for low-mass lenses. The number of detections will depend on the abundance of such FFPs as a function of mass, which is at present poorly constrained. Assuming that FFPs follow the fiducial mass function of cold, bound planets adapted from Cassan et al. (2012), we estimate that Roman will detect $\sim250$ FFPs with masses down to that of Mars (including $\sim 60$ with masses $\le M_\oplus$). We also predict that Roman will improve the upper limits on FFP populations by at least an order of magnitude compared to currently-existing constraints.

categories: astro-ph.EP astro-ph.GA astro-ph.IM astro-ph.SR


title: ZTF20aajnksq (AT2020blt): A Fast Optical Transient at $z \approx 2.9$ With No Detected Gamma-Ray Burst Counterpart

abstract: We present ZTF20aajnksq (AT2020blt), a fast-fading ($\Delta r=2.4$ mag in $\Delta t=1.3$ days) red ($g-r\approx0.6$ mag) and luminous ($M_{1626}=-25.9$) optical transient at $z=2.9$ discovered by the Zwicky Transient Facility (ZTF). AT2020blt shares several features in common with afterglows to long-duration gamma-ray bursts (GRBs): (1) an optical light curve well-described by a broken power-law with a break at $t_\mathrm{j}=1$ day (observer-frame); (2) a luminous $(L_X = 10^{46}$ $\mathrm{erg}$ $\mathrm{s}^{-1})$ X-ray counterpart; and (3) luminous ($L_\nu = 4 \times 10^{31}$ $\mathrm{erg}$ $\mathrm{sec}^{-1}$ $\mathrm{Hz}^{-1}$ at 10 GHz) radio emission. However, no GRB was detected in the 0.74d between the last ZTF non-detection ($r > 20.64$) and the first ZTF detection ($r = 19.57$), with an upper limit on the isotropic-equivalent gamma-ray energy release of $E_{\gamma,\mathrm{iso}} < 7 \times 10^{52}$ erg. AT2020blt is thus the third afterglow-like transient discovered without a detected GRB counterpart (after PTF11agg and ZTF19abvizsw) and the second (after ZTF19abvizsw) with a redshift measurement. We conclude that the properties of AT2020blt are consistent with a classical (initial Lorentz factor $\Gamma_0 \gtrsim 100$) on-axis GRB that was missed by high-energy satellites. Furthermore, by estimating the rate of transients with light curves similar to that of AT2020blt in ZTF high-cadence data, we agree with previous results that there is no evidence for an afterglow-like phenomenon that is significantly more common than classical GRBs. We conclude by discussing the status and future of fast-transient searches in wide-field high-cadence optical surveys.


title: Warping Wormholes with Dust: a Metric Construction of the Python's Lunch

abstract: We show how wormholes in three spacetime dimensions can be customizably warped using pressureless matter. In particular, we exhibit a large new class of solutions in (2+1)-dimensional general relativity with energy-momentum tensor describing a negative cosmological constant and positive-energy dust. From this class of solutions, we construct wormhole geometries and study their geometric and holographic properties, including Ryu-Takayanagi surfaces, entanglement wedge cross sections, mutual information, and outer entropy. Finally, we construct a Python's Lunch geometry: a wormhole in asymptotically anti-de Sitter space with a local maximum in size near its middle.


title: Broken into Pieces: ATLAS and Aliqa Uma as One Single Stream

abstract: We present the first spectroscopic measurements of the ATLAS and Aliqa Uma streams from the Southern Stellar Stream Spectroscopic Survey ($S^5$), in combination with the photometric data from the Dark Energy Survey and astrometric data from $Gaia$. From the coherence of spectroscopic members in radial velocity and proper motion, we find out that these two systems are extremely likely to be one stream with discontinuity in morphology and density on the sky (the "kink" feature). We refer to this entire stream as the ATLAS-Aliqa Uma stream, or the AAU stream. We perform a comprehensive exploration of the effect of baryonic substructures and find that only an encounter with the Sagittarius dwarf $\sim 0.5$ Gyr ago can create a feature similar to the observed "kink". In addition, we also identify two gaps in the ATLAS component associated with the broadening in the stream width (the "broadening" feature). These gaps have likely been created by small mass perturbers, such as dark matter halos, as the AAU stream is the most distant cold stream known with severe variations in both the stream surface density and the stream track on the sky. With the stream track, stream distance and kinematic information, we determine the orbit of the AAU stream and find that it has been affected by the Large Magellanic Cloud, resulting in a misalignment between the proper motion and stream track. Together with the Orphan-Chenab Stream, AAU is the second stream pair that has been found to be a single stream separated into two segments by external perturbation.


title: Evidence for supernova feedback sustaining gas turbulence in nearby star-forming galaxies

abstract: It is well known that gas in galaxy discs is highly turbulent, but there is much debate on which mechanism can energetically maintain this turbulence. Among the possible candidates, supernova (SN) explosions are likely the primary drivers but doubts remain on whether they can be sufficient in regions of moderate star formation activity, in particular in the outer parts of discs. In this paper, we measure the SN efficiency $\eta$, namely the fraction of the total SN energy needed to sustain turbulence in galaxies, and verify that SNe can indeed be the sole driving mechanism. The key novelty of our approach is that we take into account the increased turbulence dissipation timescale associated to the flaring in outer regions of gaseous discs. We analyse the distribution and kinematics of HI and CO in 10 nearby star-forming galaxies to obtain the radial profiles of the kinetic energy per unit area, for both the atomic gas and the molecular gas. We use a theoretical model to reproduce the observed energy with the sum of turbulent energy from SNe, as inferred from the observed star formation rate (SFR) surface density, and the gas thermal energy. We find that the observed kinetic energy is remarkably well reproduced by our model across the whole extent of the galactic discs, assuming $\eta$ constant with the galactocentric radius. Taking into account the uncertainties on the SFR surface density and on the atomic gas phase, we obtain that the median SN efficiencies for our sample of galaxies are $\langle \eta_\mathrm{atom} \rangle=0.015_{-0.008}^{+0.018}$ for the atomic gas and $\langle \eta_\mathrm{mol} \rangle = 0.003_{-0.002}^{+0.006}$ for the molecular gas. We conclude that SNe alone can sustain gas turbulence in nearby galaxies with only few percent of their energy and that there is essentially no need for any further source of energy.


title: Simulations Predicting the Ability of Multi-Color Simultaneous Photometry to Distinguish TESS Candidate Exoplanets from False Positives

abstract: The Transiting Exoplanet Survey Satellite (TESS) is currently concluding its 2-year primary science mission searching 85% of the sky for transiting exoplanets. TESS has already discovered well over one thousand TESS objects of interest (TOIs), but these candidate exoplanets must be distinguished from astrophysical false positives using other instruments or techniques. The 3-band Multi-color Simultaneous Camera for Studying Atmospheres of Transiting Planets (MuSCAT), as well as the 4-band MuSCAT2, can be used to validate TESS discoveries. Transits of exoplanets are achromatic when observed in multiple bandpasses, while transit depths for false positives often vary with wavelength. We created software tools to simulate MuSCAT/MuSCAT2 TESS follow-up observations and reveal which planet candidates can be efficiently distinguished from blended eclipsing binary (BEB) false positives using these two instruments, and which must be validated using other techniques. We applied our software code to the Barclay et al. (2018) predicted TESS discoveries, as well as to TOIs downloaded from the ExoFOP-TESS website. We estimate that MuSCAT (MuSCAT2 values in parentheses) will be able to use its multi-color capabilities to distinguish BEB false positives for $\sim$17% ($\sim$18%) of all TESS discoveries, and $\sim$13% ($\sim$15%) of $R_{\rm pl} < 4R_\oplus$ discoveries. Our TOI analysis shows that MuSCAT (MuSCAT2) can distinguish BEB false positives for $\sim$55% ($\sim$52%) of TOIs with transit depths greater than 0.001, for $\sim$64% ($\sim$61%) of TOIs with transit depths greater than 0.002, and for $\sim$70% ($\sim$68%) of TOIs with transit depth greater than 0.003. Our work shows that MuSCAT and MuSCAT2 can validate hundreds of $R_{\rm pl} < 4R_\oplus$ candidate exoplanets, thus supporting the TESS mission in achieving its Level 1 Science Requirement.


title: Critical Fermi surfaces in generic dimensions arising from transverse gauge field interactions

abstract: We study critical Fermi surfaces in generic dimensions arising from coupling finite-density fermions with transverse gauge fields, by applying the dimensional regularization scheme developed previously [Phys. Rev. B 92, 035141 (2015)]. We consider the cases of $U(1)$ and $U(1)\times U(1)$ transverse gauge couplings, and extract the nature of the renormalization group (RG) flow fixed points as well as the critical scalings. Our analysis allows us to treat a critical Fermi surface of a generic dimension $m$ perturbatively in an expansion parameter $\epsilon =\left (2-m \right ) /\left (m+1 \right).$ One of our key results is that although the two-loop corrections do not alter the existence of an RG flow fixed line for certain $U(1)\times U(1)$ theories, which was identified earlier for $m=1$ at one-loop order, the third-order diagrams do. However, this fixed line feature is also obtained for $m>1$, where the answer is one-loop exact due to UV/IR mixing.


title: Non-linearly realized discrete symmetries

abstract: While non-linear realizations of continuous symmetries feature derivative interactions and have no potential, non-linear realizations of discrete symmetries feature non-derivative interactions and have a highly suppressed potential. These Goldstone bosons of discrete symmetries have a non-zero potential, but the potential generated from quantum corrections is inherently very highly suppressed. We explore various discrete symmetries and to what extent the potential is suppressed for each of them.


title: Exploring Chemical Homogeneity in Dwarf Galaxies: A VLT-MUSE study of JKB 18

abstract: Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$\sigma$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.


title: Tracking the Footprints of Spin Fluctuations: A MultiMethod, MultiMessenger Study of the Two-Dimensional Hubbard Model

abstract: The Hubbard model represents the fundamental model for interacting quantum systems and electronic correlations. Using the two-dimensional half-filled Hubbard model at weak coupling as a testing ground, we perform a comparative study of a comprehensive set of state of the art quantum many-body methods. Upon cooling into its insulating antiferromagnetic ground-state, the model hosts a rich sequence of distinct physical regimes with crossovers between a high-temperature incoherent regime, an intermediate temperature metallic regime and a low-temperature insulating regime with a pseudogap created by antiferromagnetic fluctuations. We assess the ability of each method to properly address these physical regimes and crossovers through the computation of several observables probing both quasiparticle properties and magnetic correlations, with two numerically exact methods (diagrammatic and determinantal quantum Monte Carlo) serving as a benchmark. By combining computational results and analytical insights, we elucidate the nature and role of spin fluctuations in each of these regimes. Based on this analysis, we explain how quasiparticles can coexist with increasingly long-range antiferromagnetic correlations, and why dynamical mean-field theory is found to provide a remarkably accurate approximation of local quantities in the metallic regime. We also critically discuss whether imaginary time methods are able to capture the non-Fermi liquid singularities of this fully nested system.


title: Is HR 6819 a triple system containing a black hole? -- An alternative explanation

abstract: HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components.


title: Populating the upper black hole mass gap through stellar collisions in young star clusters

abstract: Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae. However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called "upper-mass gap" would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through collisions of young massive stars in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as $\sim20\%$ of all BH progenitors undergo one or more collisions prior to stellar collapse and up to $\sim1\%$ of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the "multiple-generation" merger channel described in other analyses. This has clear relevance for the formation of very massive BH binaries as recently detected by LIGO/Virgo in GW190521. Finally, we describe how stellar collisions in clusters may provide a unique pathway to pair-instability supernovae and briefly discuss the expected rate of these events and other electromagnetic transients.


title: Survival of the fractional Josephson effect in time-reversal-invariant topological superconductors

abstract: Time-reversal-invariant topological superconductor (TRITOPS) wires host Majorana Kramers pairs that have been predicted to mediate a fractional Josephson effect with $4\pi$ periodicity in the superconducting phase difference. We explore the TRITOPS fractional Josephson effect in the presence of time-dependent `local mixing' perturbations that instantaneously preserve time-reversal symmetry. Specifically, we show that just as such couplings render braiding of Majorana Kramers pairs non-universal, the Josephson current becomes either aperiodic or $2\pi$-periodic (depending on conditions that we quantify) unless the phase difference is swept sufficiently quickly. We further analyze topological superconductors with $\mathcal{T}^2 = +1$ time-reversal symmetry and reveal a rich interplay between interactions and local mixing that can be experimentally probed in nanowire arrays.


title: Dark matter capture in celestial objects: light mediators, self-interactions, and complementarity with direct detection

abstract: We generalize the formalism for DM capture in celestial bodies to account for arbitrary mediator mass, and update the existing and projected astrophysical constraints on DM-nucleon scattering cross section from observations of neutron stars. We show that the astrophysical constraints on the DM-nucleon interaction strength, that were thought to be the most stringent, drastically weaken for light mediators and can be completely voided. For asymmetric DM, existing astrophysical constraints are completely washed out for mediators lighter than 5 MeV, and for annihilating DM the projected constraints are washed out for mediators lighter than 0.25 MeV. Related terrestrial direct detection bounds also weaken, but in a complementary fashion; they supersede the astrophysical capture bounds for small or large DM mass, respectively for asymmetric or annihilating DM. Repulsive self-interactions of DM have an insignificant impact on the total capture rate, but a significant impact on the black hole formation criterion. This further weakens the constraints on DM-nucleon interaction strength for asymmetric self-repelling DM, whereas constraints remain unaltered for annihilating self-repelling DM. We use the correct Hawking evaporation rate of the newly formed black hole, that was approximated as a blackbody in previous studies, and show that, despite a more extensive alleviation of collapse as a result, the observation of a neutron star collapse can probe a wide range of DM self-interaction strengths.


title: Black Rubber and the Non-linear Elastic Response of Scale Invariant Solids

abstract: We discuss the nonlinear elastic response in scale invariant solids. Following previous work, we split the analysis into two basic options: according to whether scale invariance (SI) is a manifest or a spontaneously broken symmetry. In the latter case, one can employ effective field theory methods, whereas in the former we use holographic methods. We focus on a simple class of holographic models that exhibit elastic behaviour, and obtain their nonlinear stress-strain curves as well as an estimate of the elasticity bounds - the maximum possible deformation in the elastic (reversible) regime. The bounds differ substantially in the manifest or spontaneously broken SI cases, even when the same stress-strain curve is assumed in both cases. Additionally, the hyper-elastic subset of models (that allow for large deformations) is found to have stress-strain curves akin to natural rubber. The holographic instances in this category, which we dub black rubber, display richer stress-strain curves - with two different power-law regimes at different magnitudes of the strain.


title: The relative orientation between the magnetic field and gradients of surface brightness within thin velocity slices of 12CO and 13CO emission from the Taurus molecular cloud

abstract: We examine the role of the interstellar magnetic field to modulate the orientation of turbulent flows within the Taurus molecular cloud using spatial gradients of thin velocity slices of 12CO and 13CO antenna temperatures. Our analysis accounts for the random errors of the gradients that arise from the thermal noise of the spectra. The orientations of the vectors normal to the antenna temperature gradient vectors are compared to the magnetic field orientations that are calculated from Planck 353~GHz polarization data. These relative orientations are parameterized with the projected Rayleigh statistic and mean resultant vector. For 12CO, 28 percent and 39 percent of the cloud area exhibit strongly parallel or strongly perpendicular relative orientations respectively. For the lower opacity 13CO emission, strongly parallel and strongly perpendicular orientations are found in 7 percent and 43 percent of the cloud area respectively. For both isotopologues, strongly parallel or perpendicular alignments are restricted to localized regions with low levels of turbulence. If the relative orientations serve as an observational proxy to the Alfvenic Mach number then our results imply local variations of the Alfvenic Mach number throughout the cloud.


title: From magnetoelectric response to optical activity

abstract: We apply a microscopic theory of polarization and magnetization to crystalline insulators at zero temperature and consider the orbital electronic contribution of the linear response to spatially varying, time-dependent electromagnetic fields. The charge and current density expectation values generally depend on both the microscopic polarization and magnetization fields, and on the microscopic free charge and current densities. But contributions from the latter vanish in linear response for the class of insulators we consider. Thus we need only consider the former, which can be decomposed into "site" polarization and magnetization fields, from which "site multipole moments" can be constructed. Macroscopic polarization and magnetization fields follow, and we identify the relevant contributions to them; for electromagnetic fields varying little over a lattice constant these are the electric and magnetic dipole moments per unit volume, and the electric quadrupole moment per unit volume. A description of optical activity and related magneto-optical phenomena follows from the response of these macroscopic quantities to the electromagnetic field and, while in this paper we work within the independent particle and frozen-ion approximations, both optical rotary dispersion and circular dichroism can be described with this strategy. Earlier expressions describing the magnetoelectric effect are recovered as the zero frequency limit of our more general equations. Since our site quantities are introduced with the use of Wannier functions, the site multipole moments and their macroscopic analogs are generally gauge dependent. However, the resulting macroscopic charge and current densities, together with the optical effects to which they lead, are gauge invariant, as would be physically expected.


title: KiDS+GAMA: The weak lensing calibrated stellar-to-halo mass relation of central and satellite galaxies

abstract: We simultaneously present constraints on the stellar-to-halo mass relation for central and satellite galaxies through a weak lensing analysis of spectroscopically classified galaxies. Using overlapping data from the fourth data release of the Kilo-Degree Survey (KiDS), and the Galaxy And Mass Assembly survey (GAMA), we find that satellite galaxies are hosted by halo masses that are $0.53 \pm 0.39$ dex (68\% confidence, $3\sigma$ detection) smaller than those of central galaxies of the same stellar mass (for a stellar mass of $\log(M_{\star}/M_{\odot}) = 10.6$). This is consistent with galaxy formation models, whereby infalling satellite galaxies are preferentially stripped of their dark matter. We find consistent results with similar uncertainties when comparing constraints from a standard azimuthally averaged galaxy-galaxy lensing analysis and a two-dimensional likelihood analysis of the full shear field. As the latter approach is somewhat biased due to the lens incompleteness and as it does not provide any improvement to the precision when applied to actual data, we conclude that stacked tangential shear measurements are best-suited for studies of the galaxy-halo connection.


title: An adapted filter function for density split statistics in weak lensing

abstract: Context. The density split statistics in weak gravitational lensing analyses probes the correlation between regions of different (foreground) galaxy number densities and their weak lensing signal, measured by the shape distortion of background galaxies. Aims. In this paper, we reconsider density split statistics, by constructing a new angular filter function that is adapted to the expected relation between galaxy number density and shear pattern, in a way that the filter weighting the galaxy number density is matched to the filter that is used to quantify the shear signal. Methods. We use the results of numerical ray-tracing simulations, specifically through Millennium Simulation supplemented by a galaxy distribution based on a semi-analytic model, to construct a matched pair of adapted filter functions for the galaxy density and the tangential shear signal. We compare the performance of our new filter to the previously used top-hat filter, applying both to a different and independent set of numerical simulations (SLICS, cosmo-SLICS). Results. We show that the adapted filter yields a better correlation between the total matter and the galaxy distribution. Furthermore, the adapted filter provides a larger signal-to-noise ratio to constrain the bias between the total matter and the galaxy distribution, and we show that it is, in general, a more sensitive discriminator between different cosmologies, with the exception of cosmologies with very large $\sigma_8$ values. All analyses lead to the conclusion that our adapted filter should be favored in future density split statistic works.


title: Speck of Chaos

abstract: It has been shown that, despite being local, a perturbation applied to a single site of the one-dimensional XXZ model is enough to bring this interacting integrable spin-1/2 system to the chaotic regime. Here, we show that this is not unique to this model, but happens also to the Ising model in a transverse field and to the spin-1 Lai-Sutherland chain. The larger the system is, the smaller the amplitude of the local perturbation for the onset of chaos. We focus on two indicators of chaos, the correlation hole, which is a dynamical tool, and the distribution of off-diagonal elements of local observables, which is used in the eigenstate thermalization hypothesis. Both methods avoid spectrum unfolding and can detect chaos even when the eigenvalues are not separated by symmetry sectors.


title: Constraining Non-Relativistic RG Flows with Holography

abstract: We examine non-relativistic holographic RG flows by working with Einstein-Maxwell-scalar theories which support geometries that break Lorentz invariance at some energy scale. We adopt the superpotential formalism, which helps us characterize the radial flow in this setup and bring to light a number of generic features. In particular, we identify several quantities that behave monotonically under RG flow. As an example, we show that the index of refraction is generically monotonic. We also construct a combination of the superpotentials that flows monotonically in Einstein-scalar theories supporting non-relativistic solutions, and which reduces to the known c-function in the relativistic limit. Interestingly, such quantity also exhibits monotonicity in a variety of black hole solutions to the full Einstein-Maxwell-scalar theory, hinting at a deeper structure. Finally, we comment on the breakdown of such monotonicity conditions and on the relation to a candidate c-function obtained previously from entanglement entropy.


title: Let Effective Field Theory of Inflation flow: stochastic generation of models with red/blue tensor tilt

abstract: We extend the method of Flow Equations to the Effective Field Theory framework of inflation, in order to investigate the observable predictions of a very broad class of inflationary models. Focusing our attention on the gravitational-wave sector, we derive a general expression for the consistency relation for effective models and provide a numerical implementation which allows to study how the generated models populate the $(r,nt)$ plane. We analyse $5 \times 10^{4}$ realizations of inflationary scenarios that respect the Null-Energy Condition ($\epsilon > 0$) and $5 \times 10^{4}$ realizations that violate it ($\epsilon > 0$). In both cases, 90% of the viable models are below the most recent upper bound on the tensor-to-scalar ratio from Planck and BICEP2/Keck Array BK15 data: $r_{0.002} < 0.056$ at 95 % CL. We find that general EFT inflationary models with $\epsilon > 0$ are typically characterized by $n_{t} < 0$, whereas the vast majority of NEC-violating models presents a blue-tilted spectrum ($n_{t} > 0$). Since a blue tensor spectral index implies more power on small scales, this result is of considerable interest in view of a possible direct detection of the primordial gravitational-wave background.


title: AI Feynman 2.0: Pareto-optimal symbolic regression exploiting graph modularity

abstract: We present an improved method for symbolic regression that seeks to fit data to formulas that are Pareto-optimal, in the sense of having the best accuracy for a given complexity. It improves on the previous state-of-the-art by typically being orders of magnitude more robust toward noise and bad data, and also by discovering many formulas that stumped previous methods. We develop a method for discovering generalized symmetries (arbitrary modularity in the computational graph of a formula) from gradient properties of a neural network fit. We use normalizing flows to generalize our symbolic regression method to probability distributions from which we only have samples, and employ statistical hypothesis testing to accelerate robust brute-force search.

categories: cs.LG cs.AI cs.IT math.IT physics.comp-ph stat.ML


title: Quiver Mutations, Seiberg Duality and Machine Learning

abstract: We initiate the study of applications of machine learning to Seiberg duality, focusing on the case of quiver gauge theories, a problem also of interest in mathematics in the context of cluster algebras. Within the general theme of Seiberg duality, we define and explore a variety of interesting questions, broadly divided into the binary determination of whether a pair of theories picked from a series of duality classes are dual to each other, as well as the multi-class determination of the duality class to which a given theory belongs. We study how the performance of machine learning depends on several variables, including number of classes and mutation type (finite or infinite). In addition, we evaluate the relative advantages of Naive Bayes classifiers versus Convolutional Neural Networks. Finally, we also investigate how the results are affected by the inclusion of additional data, such as ranks of gauge/flavor groups and certain variables motivated by the existence of underlying Diophantine equations. In all questions considered, high accuracy and confidence can be achieved.


title: Spin-orbit alignment of the $\beta$ Pictoris planetary system

abstract: A crucial diagnostic that can tell us about processes involved in the formation and dynamical evolutionof planetary systems is the angle between the rotation axis of a star and a planet's orbital angular momentum vector ("spin-orbit" alignment or "obliquity"). Here we present the first spin-orbit alignment measurement for a wide-separation exoplanetary system, namely on the directly-imaged planet $\beta$ Pictoris b. We use VLTI/GRAVITY spectro-interferometry with an astrometric accuracy of 1 $\mu$as (microarcsecond) in the Br$\gamma$ photospheric absorption line to measure the photocenter displacement associated with the stellar rotation. Taking inclination constraints from astroseismology into account, we constrain the 3-dimensional orientation of the stellar spin axis and find that $\beta$ Pic b orbits its host star on a prograde orbit. The angular momentum vectors of the stellar photosphere, the planet, and the outer debris disk are well-aligned with mutual inclinations $<3\pm5^{\circ}$, which indicates that $\beta$ Pic b formed in a system without significant primordial misalignments. Our results demonstrate the potential of infrared interferometry to measure the spin-orbit alignment for wide-separation planetary systems, probing a highly complementary regime to the parameter space accessible with the Rossiter-McLaughlin effect. If the low obliquity is confirmed by measurements on a larger sample of wide-separation planets, it would lend support to theories that explain the obliquity in Hot Jupiter systems with dynamical scattering and the Kozai-Lidov mechanism.


title: UK COVID-19 Lockdown: What are the impacts on air pollution

abstract: A country-wide `lock-down' imposed on the 23red March 2020 in the UK had a significant impact on the UK's anthropogenic movements. The closure of work-places and restrictions imposed on visiting friends and family has radically reduced the amount of traffic on the roads. In this short communication, we use data from UK air-quality sensors to quantify air pollution trends pre- and post-lock-down. While we detect large falls in nitrogen dioxide at levels not seen over the last decade, trends in other pollutants were mixed especially when compared to historical data. It suggests that the implication that lock-down was beneficial for the environment was not so obvious.


title: Understanding the electromagnetic response of graphene/metallic nanostructures hybrids of different dimensionality

abstract: Plasmonic excitations such as surface-plasmon-polaritons (SPPs) and graphene-plasmons (GPs), carry large momenta and are thus able to confine electromagnetic fields to small dimensions. This property makes them ideal platforms for subwavelength optical control and manipulation at the nanoscale. The momenta of these plasmons are even further increased if a scheme of metal-insulator-metal and graphene-insulator-metal are used for SPPs and GPs, respectively. However, with such large momenta, their far-field excitation becomes challenging. In this work, we consider hybrids of graphene and metallic nanostructures and study the physical mechanisms behind the interaction of far-field light with the supported high momenta plasmon modes. While there are some similarities in the properties of GPs and SPPs, since both are of the plasmon-polariton type, their physical properties are also distinctly different. For GPs we find two different physical mechanism related to either GPs confined to isolated cavities, or large area collective grating couplers. Strikingly, we find that although the two systems are conceptually different, under specific conditions they can behave similarly. By applying the same study to SPPs, we find a different physical behavior, which fundamentally stems from the different dispersion relations of SPPs as compared to GPs. Furthermore, these hybrids produce large field enhancements that can also be electrically tuned and modulated making them the ideal candidates for a variety of plasmonic devices.


title: Construction of potential functions associated with a given energy spectrum -- An inverse problem. II

abstract: We continue our solution of the inverse problem started by the first author in [Int. J. Mod. Phys. A 35, xxxx (2020), in production]. Additional potential functions for exactly solvable problems that correspond to the same energy spectrum formula but for different energy polynomials and bases are found. In this work, we obtain a class of potential functions associated with the Wilson polynomial and "Jacobi basis".


title: Classification of critically fixed anti-rational maps

abstract: We show that there is a one-to-one correspondence between conjugacy classes of critically fixed anti-rational maps and equivalence classes of certain plane graphs. We furthermore prove that critically fixed anti-rational maps are combinatorially equivalent to "topological Schottky maps" associated to these graphs, given in each face by a topological reflection in its boundary. As a corollary, we obtain a similar classification of critically fixed anti-polynomials by certain plane trees. One of the main technical tools is an analogue of Thurston's characterization of rational maps in the orientation-reversing case.


title: Combinatorics of antiprism triangulations

abstract: The antiprism triangulation provides a natural way to subdivide a simplicial complex $\Delta$, similar to barycentric subdivision, which appeared independently in combinatorial algebraic topology and computer science. It can be defined as the simplicial complex of chains of multi-pointed faces of $\Delta$, from a combinatorial point of view, and by successively applying the antiprism construction, or balanced stellar subdivisions, on the faces of $\Delta$, from a geometric point of view. This paper studies enumerative invariants associated to this triangulation, such as the transformation of the $h$-vector of $\Delta$ under antiprism triangulation, and algebraic properties of its Stanley--Reisner ring. Among other results, it is shown that the $h$-polynomial of the antiprism triangulation of a simplex is real-rooted and that the antiprism triangulation of $\Delta$ has the almost strong Lefschetz property over ${\mathbb R}$ for every shellable complex $\Delta$. Several related open problems are discussed.


title: Quantitative non-divergence and lower bounds for points with algebraic coordinates near manifolds

abstract: Point counting estimates are a key stepping stone to various results in metric Diophantine approximation. In this paper we use the quantitative non-divergence estimates originally developed by Kleinbock and Margulis to improve lower bounds by Bernik, G\"{o}tze et al. for the number of points with algebraic conjugate coordinates close to a given manifold. In the process, we also improve on a Khinchin-Groshev-type theorem for a problem of constrained approximation by polynomials.


title: Characterization of space-momentum entangled photon with a time resolving CMOS SPAD array

abstract: Single photon avalanche diode arrays can provide both the spatial and temporal information of each detected photon. We present here the characterization of entangled light with a sensor specifically designed for quantum imaging applications. The sensors is time-tagging each detection events at the pixel level with sub-nanosecond accuracy, within frames of 50 ns. The spatial correlations between any number of detections in a defined temporal window can thus be directly extracted from the data. We show the ability of the sensor to characterize space-momentum entangled photon pairs emitted by spontaneous parametric downconversion. Their entanglement is demonstrated by violating an EPR-type inequality.


title: Bootstrapping Complete The Look at Pinterest

abstract: Putting together an ideal outfit is a process that involves creativity and style intuition. This makes it a particularly difficult task to automate. Existing styling products generally involve human specialists and a highly curated set of fashion items. In this paper, we will describe how we bootstrapped the Complete The Look (CTL) system at Pinterest. This is a technology that aims to learn the subjective task of "style compatibility" in order to recommend complementary items that complete an outfit. In particular, we want to show recommendations from other categories that are compatible with an item of interest. For example, what are some heels that go well with this cocktail dress? We will introduce our outfit dataset of over 1 million outfits and 4 million objects, a subset of which we will make available to the research community, and describe the pipeline used to obtain and refresh this dataset. Furthermore, we will describe how we evaluate this subjective task and compare model performance across multiple training methods. Lastly, we will share our lessons going from experimentation to working prototype, and how to mitigate failure modes in the production environment. Our work represents one of the first examples of an industrial-scale solution for compatibility-based fashion recommendation.


title: Reflection in Game-Based Learning: A Survey of Programming Games

abstract: Reflection is a critical aspect of the learning process. However, educational games tend to focus on supporting learning concepts rather than supporting reflection. While reflection occurs in educational games, the educational game design and research community can benefit from more knowledge of how to facilitate player reflection through game design. In this paper, we examine educational programming games and analyze how reflection is currently supported. We find that current approaches prioritize accuracy over the individual learning process and often only support reflection post-gameplay. Our analysis identifies common reflective features, and we develop a set of open areas for future work. We discuss these promising directions towards engaging the community in developing more mechanics for reflection in educational games.


title: Inferences About the Magnetic Field Structure of a CME with Both In Situ and Faraday Rotation Constraints

abstract: On 2012 August 2, two CMEs (CME-1 and CME-2) erupted from the west limb of the Sun as viewed from Earth, and were observed in images from the white light coronagraphs on the SOHO and STEREO spacecraft. These events were also observed by the Very Large Array (VLA), which was monitoring the Sun at radio wavelengths, allowing time-dependent Faraday rotation observations to be made of both events. We use the white-light imaging and radio data to model the 3-D field geometry of both CMEs, assuming a magnetic flux rope geometry. For CME-2, we also consider 1 au in situ field measurements in the analysis, as this CME hits STEREO-A on August~6, making this the first CME with observational constraints from stereoscopic coronal imaging, radio Faraday rotation, and in situ plasma measurements combined. The imaging and in situ observations of CME-2 provide two clear predictions for the radio data; namely that VLA should observe positive rotation measures (RMs) when the radio line of sight first encounters the CME, and that the sign should reverse to negative within a couple hours. The initial positive RMs are in fact observed. The expected sign reversal is not, but the VLA data unfortunately end too soon to be sure of the significance of this discrepancy. We interpret an RM increase prior to the expected occultation time of the CME as a signature of a sheath region of deflected field ahead of the CME itself.


title: Tunable anomalous Hall transport in bulk and two-dimensional 1$T$-CrTe$_{2}$: A first-principles study

abstract: Layered materials with robust magnetic ordering have been attracting significant research interest. In recent experiments, a new layered material 1$T$-CrTe$_{2}$ has been synthesized and exhibits ferromagnetism above the room temperature. Here, based on first-principles calculations, we investigate the electronic, magnetic, and transport properties of 1$T$-CrTe$_{2}$, both in the bulk and in the two-dimensional (2D) limit. We show that 1$T$-CrTe$_{2}$ can be stable in the monolayer form, and has a low exfoliation energy. The monolayer structure is an intrinsic ferromagnetic metal, which maintains a high Curie temperature above the room temperature. Particularly, we reveal interesting features in the anomalous Hall transport. We show that in the ground state, both bulk and monolayer 1$T$-CrTe$_{2}$ possess vanishing anomalous Hall effect, because the magnetization preserves one vertical mirror symmetry. The anomalous Hall conductivity can be made sizable by tuning the magnetization direction or by uniaxial strains that break the mirror symmetry. The room-temperature 2D ferromagnetism and the tunable anomalous Hall effect make the material a promising platform for nanoscale device applications.


title: Null sets and combinatorial covering properties

abstract: A subset of the Cantor cube is null-additive if its algebraic sum with any null set is null. We construct a set of cardinality continuum such that: all continuous images of the set into the Cantor cube are null-additive, it contains a homeomorphic copy of a set that is not null-additive, and it has the property $\gamma$, a strong combinatorial covering property. We also construct a nontrivial subset of the Cantor cube with the property $\gamma$ that is not null additive. Set-theoretic assumptions used in our constructions are far milder than used earlier by Galvin--Miller and Bartoszy\'nski--Rec{\l}aw, to obtain sets with analogous properties. We also consider products of Sierpi\'nski sets in the context of combinatorial covering properties.


title: On the Manhattan pinball problem

abstract: We consider the periodic Manhattan lattice with alternating orientations going north-south and east-west. Place obstructions on vertices independently with probability $0<p<1$. A particle is moving on the edges with unit speed following the orientation of the lattice and it will turn only when encountering an obstruction. The problem is that for which value of $p$ is the trajectory of the particle closed almost surely. We prove this for $p>\frac{1}{2}-\varepsilon$ with some $\varepsilon>0$.


title: A Gaussian particle distribution for branching Brownian motion with an inhomogeneous branching rate

abstract: Motivated by the goal of understanding the evolution of populations undergoing selection, we consider branching Brownian motion in which particles independently move according to one-dimensional Brownian motion with drift, each particle may either split into two or die, and the difference between the birth and death rates is a linear function of the position of the particle. We show that, under certain assumptions, after a sufficiently long time, the empirical distribution of the positions of the particles is approximately Gaussian. This provides mathematically rigorous justification for results in the biology literature indicating that the distribution of the fitness levels of individuals in a population over time evolves like a Gaussian traveling wave.


title: Higher order analysis on the existence of periodic solutions in continuous differential equations via degree theory

abstract: Recently, the higher order averaging method for studying periodic solutions of both Lipschitz differential equations and discontinuous piecewise smooth differential equations was developed in terms of Brouwer degree theory. Between the Lipschitz and the discontinuous piecewise smooth differential equations, there is a huge class of differential equations lacking in a higher order analysis on the existence of periodic solutions, namely the class of continuous non-Lipschitz differential equations. In this paper, based on the coincidence degree theory for nonlinear operator equations, we perform a higher order analysis of continuous (non-Lipschitz) perturbed differential equations and derive sufficient conditions for the existence of periodic solutions for such systems. We apply our results to study continuous (non-Lipschitz) higher order perturbations of a harmonic oscillator.


title: Weighted QMIX: Expanding Monotonic Value Function Factorisation for Deep Multi-Agent Reinforcement Learning

abstract: QMIX is a popular $Q$-learning algorithm for cooperative MARL in the centralised training and decentralised execution paradigm. In order to enable easy decentralisation, QMIX restricts the joint action $Q$-values it can represent to be a monotonic mixing of each agent's utilities. However, this restriction prevents it from representing value functions in which an agent's ordering over its actions can depend on other agents' actions. To analyse this representational limitation, we first formalise the objective QMIX optimises, which allows us to view QMIX as an operator that first computes the $Q$-learning targets and then projects them into the space representable by QMIX. This projection returns a representable $Q$-value that minimises the unweighted squared error across all joint actions. We show in particular that this projection can fail to recover the optimal policy even with access to $Q^*$, which primarily stems from the equal weighting placed on each joint action. We rectify this by introducing a weighting into the projection, in order to place more importance on the better joint actions. We propose two weighting schemes and prove that they recover the correct maximal action for any joint action $Q$-values, and therefore for $Q^*$ as well. Based on our analysis and results in the tabular setting, we introduce two scalable versions of our algorithm, Centrally-Weighted (CW) QMIX and Optimistically-Weighted (OW) QMIX and demonstrate improved performance on both predator-prey and challenging multi-agent StarCraft benchmark tasks.


title: Predictive Complexity Priors

abstract: Specifying a Bayesian prior is notoriously difficult for complex models such as neural networks. Reasoning about parameters is made challenging by the high-dimensionality and over-parameterization of the space. Priors that seem benign and uninformative can have unintuitive and detrimental effects on a model's predictions. For this reason, we propose predictive complexity priors: a functional prior that is defined by comparing the model's predictions to those of a reference model. Although originally defined on the model outputs, we transfer the prior to the model parameters via a change of variables. The traditional Bayesian workflow can then proceed as usual. We apply our predictive complexity prior to high-dimensional regression, reasoning over neural network depth, and sharing of statistical strength for few-shot learning.


title: DS6, Deformation-aware Semi-supervised Learning: Application to Small Vessel Segmentation with Noisy Training Data

abstract: Blood vessels of the brain provide the human brain with the required nutrients and oxygen. As a vulnerable part of the cerebral blood supply, pathology of small vessels can cause serious problems such as Cerebral Small Vessel Diseases (CSVD). It has also been shown that CSVD is related to neurodegeneration, such as Alzheimer's disease. With the advancement of 7 Tesla MRI systems, higher spatial image resolution can be achieved, enabling the depiction of very small vessels in the brain. Non-Deep Learning-based approaches for vessel segmentation, e.g., Frangi's vessel enhancement with subsequent thresholding, are capable of segmenting medium to large vessels but often fail to segment small vessels. The sensitivity of these methods to small vessels can be increased by extensive parameter tuning or by manual corrections, albeit making them time-consuming, laborious, and not feasible for larger datasets. This paper proposes a deep learning architecture to automatically segment small vessels in 7 Tesla 3D Time-of-Flight (ToF) Magnetic Resonance Angiography (MRA) data. The algorithm was trained and evaluated on a small imperfect semi-automatically segmented dataset of only 11 subjects; using six for training, two for validation, and three for testing. The deep learning model based on U-Net Multi-Scale Supervision was trained using the training subset and was made equivariant to elastic deformations in a self-supervised manner using deformation-aware learning to improve the generalisation performance. The proposed technique was evaluated quantitatively and qualitatively against the test set and achieved a Dice score of 80.44 $\pm$ 0.83. Furthermore, the result of the proposed method was compared against a selected manually segmented region (62.07 resultant Dice) and has shown a considerable improvement (18.98\%) with deformation-aware learning.


title: Supervision Accelerates Pre-training in Contrastive Semi-Supervised Learning of Visual Representations

abstract: We investigate a strategy for improving the efficiency of contrastive learning of visual representations by leveraging a small amount of supervised information during pre-training. We propose a semi-supervised loss, SuNCEt, based on noise-contrastive estimation and neighbourhood component analysis, that aims to distinguish examples of different classes in addition to the self-supervised instance-wise pretext tasks. On ImageNet, we find that SuNCEt can be used to match the semi-supervised learning accuracy of previous contrastive approaches while using less than half the amount of pre-training and compute. Our main insight is that leveraging even a small amount of labeled data during pre-training, and not only during fine-tuning, provides an important signal that can significantly accelerate contrastive learning of visual representations. Our code is available online at github.com/facebookresearch/suncet.


title: Karp's patching algorithm on dense digraphs

abstract: We consider the following question. We are given a dense digraph $D$ with minimum in- and out-degree at least $\alpha n$, where $\alpha>1/2$ is a constant. The edges of $D$ are given edge costs $C(e),e\in E(D)$, where $C(e)$ is an independent copy of the uniform $[0,1]$ random variable $U$. Let $C(i,j),i,j\in[n]$ be the associated $n\times n$ cost matrix where $C(i,j)=\infty$ if $(i,j)\notin E(D)$. We show that w.h.p. the patching algorithm of Karp finds a tour for the asymmetric traveling salesperson problem that is asymptotically equal to that of the associated assignment problem. Karp's algorithm runs in polynomial time.


title: Pervasive Communications Technologies For Managing Pandemics

abstract: Pandemics always have had serious consequences unless they were effectively contained. Recent experiences with COVID-19 show that by using a smart and swift approach to deal with pandemics, avoids overwhelming of healthcare systems, and reduces the loss of precious life. This paper is about using smart technologies such as Mobile Edge Clouds (MEC), Internet of Things (IoT), and Artificial Intelligence (AI), as an approach to effectively manage pandemics. IoT provides pervasive connectivity among various devices and can be used for collecting information such as location and symptoms of potentially infected individuals. MECs provide cloud services on the edge, integrating IoT infrastructure and execution of sophisticated AI algorithms in the Cloud. In this paper, we develop a prototype to demonstrate the convergence of pervasive technologies to support research in managing pandemics. Low-cost Single Board Computers (SBC) based clusters are integrated within MEC to support remote medical teams in the field. The prototype implements a lightweight Docker container orchestrated by Kubernetes eco-system which is deployed on the clusters. The prototype successfully demonstrates that mobile medical facilities can utilize the proposed solution to collect information and execute AI algorithms while on the go. Finally, we present a discussion on the role of converging pervasive technologies on managing pandemics.


title: Optical detection of the 1.1-day variability at the white dwarf GD 394 with TESS

abstract: Recent discoveries have demonstrated that planetary systems routinely survive the post-main sequence evolution of their host stars, leaving the resulting white dwarf with a rich circumsteller environment. Among the most intriguing of such hosts is the hot white dwarf GD 394, exhibiting a unique $1.150\pm0.003$ d flux variation detected in Extreme Ultraviolet Explorer (EUVE) observations in the mid 1990s. The variation has eluded a satisfactory explanation, but hypotheses include channeled accretion producing a dark spot of metals, occultation by a gas cloud from an evaporating planet, or heating from a flux tube produced by an orbiting iron-cored planetesimal. We present observations obtained with the Transiting Exoplanet Survey Satellite (TESS) of GD 394. The space-based optical photometry demonstrates a $0.12\pm0.01$ % flux variation with a period of $1.146\pm 0.001$ d, consistent with the EUVE period and the first re-detection of the flux variation outside of the extreme ultraviolet. We describe the analysis of the TESS light curve and measurement of the optical variation, and discuss the implications of our results for the various physical explanations put forward for the variability of GD 394.


title: Semantic Linking Maps for Active Visual Object Search

abstract: We aim for mobile robots to function in a variety of common human environments. Such robots need to be able to reason about the locations of previously unseen target objects. Landmark objects can help this reasoning by narrowing down the search space significantly. More specifically, we can exploit background knowledge about common spatial relations between landmark and target objects. For example, seeing a table and knowing that cups can often be found on tables aids the discovery of a cup. Such correlations can be expressed as distributions over possible pairing relationships of objects. In this paper, we propose an active visual object search strategy method through our introduction of the Semantic Linking Maps (SLiM) model. SLiM simultaneously maintains the belief over a target object's location as well as landmark objects' locations, while accounting for probabilistic inter-object spatial relations. Based on SLiM, we describe a hybrid search strategy that selects the next best view pose for searching for the target object based on the maintained belief. We demonstrate the efficiency of our SLiM-based search strategy through comparative experiments in simulated environments. We further demonstrate the real-world applicability of SLiM-based search in scenarios with a Fetch mobile manipulation robot.


title: Data-Driven Game Development: Ethical Considerations

abstract: In recent years, the games industry has made a major move towards data-driven development, using data analytics and player modeling to inform design decisions. Data-driven techniques are beneficial as they allow for the study of player behavior at scale, making them very applicable to modern digital game development. However, with this move towards data driven decision-making comes a number of ethical concerns. Previous work in player modeling as well as work in the fields of AI and machine learning have demonstrated several ways in which algorithmic decision-making can be flawed due to data or algorithmic bias or lack of data from specific groups. Further, black box algorithms create a trust problem due to lack of interpretability and transparency of the results or models developed based on the data, requiring blind faith in the results. In this position paper, we discuss several factors affecting the use of game data in the development cycle. In addition to issues raised by previous work, we also raise issues with algorithms marginalizing certain player groups and flaws in the resulting models due to their inability to reason about situational factors affecting players' decisions. Further, we outline some work that seeks to address these problems and identify some open problems concerning ethics and game data science.


title: Environmental Dependence of Self-Regulating Black-hole Feedback in Massive Galaxies

abstract: In the universe's most massive galaxies, kinetic feedback from a central supermassive black hole appears to limit star formation. Abundant circumstantial evidence suggests that accumulation of cold gas near the central black hole strongly boosts the feedback output, keeping the ambient medium in a state marginally unstable to condensation and formation of cold gas clouds. However, the ability of that mechanism to self-regulate may depend on numerous environmental factors, including the depth of the potential well and the pressure of the surrounding circumgalactic medium (CGM). Here we present a suite of numerical simulations that explores the dependence of cold-fuelled bipolar kinetic feedback on those environmental factors. Halo mass in this simulation suite ranges from $2 \times 10^{12} \, M_\odot$ to $8 \times 10^{14} \, M_\odot$. We include the spatially extended mass and energy input from the massive galaxy's old stellar population, which is capable of sweeping gas out of the galaxy and away from the central black hole if the confining CGM pressure is sufficiently low. Our simulations show that this feedback mechanism is tightly self-regulating in a massive galaxy with a deep central potential and low CGM pressure, permitting only small amounts of multiphase gas to accumulate and allowing almost no star formation. In a massive galaxy of similar mass but a shallower central potential and greater CGM pressure the same feedback mechanism is more episodic, producing extended multiphase gas and occasionally allowing small rates of star formation ($\sim 0.1 \, M_\odot \, {\rm yr}^{-1}$). At the low-mass end of the explored range the mechanism becomes implausibly explosive, perhaps because the ambient gas initially has no angular momentum, which would have reduced the amount of condensed gas capable of fueling feedback.


title: Reparameterized Variational Divergence Minimization for Stable Imitation

abstract: While recent state-of-the-art results for adversarial imitation-learning algorithms are encouraging, recent works exploring the imitation learning from observation (ILO) setting, where trajectories \textit{only} contain expert observations, have not been met with the same success. Inspired by recent investigations of $f$-divergence manipulation for the standard imitation learning setting(Ke et al., 2019; Ghasemipour et al., 2019), we here examine the extent to which variations in the choice of probabilistic divergence may yield more performant ILO algorithms. We unfortunately find that $f$-divergence minimization through reinforcement learning is susceptible to numerical instabilities. We contribute a reparameterization trick for adversarial imitation learning to alleviate the optimization challenges of the promising $f$-divergence minimization framework. Empirically, we demonstrate that our design choices allow for ILO algorithms that outperform baseline approaches and more closely match expert performance in low-dimensional continuous-control tasks.


title: Learning to infer in recurrent biological networks

abstract: A popular theory of perceptual processing holds that the brain learns both a generative model of the world and a paired recognition model using variational Bayesian inference. Most hypotheses of how the brain might learn these models assume that neurons in a population are conditionally independent given their common inputs. This simplification is likely not compatible with the type of local recurrence observed in the brain. Seeking an alternative that is compatible with complex inter-dependencies yet consistent with known biology, we argue here that the cortex may learn with an adversarial algorithm. Many observable symptoms of this approach would resemble known neural phenomena, including wake/sleep cycles and oscillations that vary in magnitude with surprise, and we describe how further predictions could be tested. We illustrate the idea on recurrent neural networks trained to model image and video datasets. This framework for learning brings variational inference closer to neuroscience and yields multiple testable hypotheses.


title: Overgroups of regular unipotent elements in simple algebraic groups

abstract: We investigate positive-dimensional closed reductive subgroups of almost simple algebraic groups containing a regular unipotent element. Our main result states that such subgroups do not lie inside proper parabolic subgroups unless possibly when their connected component is a torus. This extends the earlier result of Testerman and Zalesski treating connected reductive subgroups.


title: From spirals to lenticulars: evidence from the rotation curves and mass models of three early-type galaxies

abstract: Rotation curves have traditionally been difficult to trace for early-type galaxies (ETGs) because they often lack a high-density disk of cold gas as in late-type galaxies (LTGs). We derive rotation curves for three lenticular galaxies from the ATLAS3D survey, combining CO data in the inner parts with deep HI data in the outer regions, extending out to 10-20 effective radii. We also use Spitzer photometry at 3.6 um to decompose the rotation curves into the contributions of baryons and dark matter (DM). We find that (1) the rotation-curve shapes of these ETGs are similar to those of LTGs of similar mass and surface brightness; (2) the dynamically$-$inferred stellar mass-to-light ratios are small for quiescent ETGs but similar to those of star-forming LTGs; (3) the DM halos follow the same scaling relations with galaxy luminosity as those of LTGs; (4) one galaxy (NGC 3626) is poorly fitted by cuspy DM profiles, suggesting that DM cores may exist in high-mass galaxies too. Our results indicate that these lenticular galaxies have recently transitioned from LTGs to ETGs without altering their DM halo structure (e.g., via a major merger) and could be faded spirals. We also confirm that ETGs follow the same radial acceleration relation as LTGs, reinforcing the notion that this is a universal law for all galaxy types.


title: FLAMBE: Structural Complexity and Representation Learning of Low Rank MDPs

abstract: In order to deal with the curse of dimensionality in reinforcement learning (RL), it is common practice to make parametric assumptions where values or policies are functions of some low dimensional feature space. This work focuses on the representation learning question: how can we learn such features? Under the assumption that the underlying (unknown) dynamics correspond to a low rank transition matrix, we show how the representation learning question is related to a particular non-linear matrix decomposition problem. Structurally, we make precise connections between these low rank MDPs and latent variable models, showing how they significantly generalize prior formulations for representation learning in RL. Algorithmically, we develop FLAMBE, which engages in exploration and representation learning for provably efficient RL in low rank transition models.


title: Automatically Learning Compact Quality-aware Surrogates for Optimization Problems

abstract: Solving optimization problems with unknown parameters often requires learning a predictive model to predict the values of the unknown parameters and then solving the problem using these values. Recent work has shown that including the optimization problem as a layer in the model training pipeline results in predictions of the unobserved parameters that lead to higher decision quality. Unfortunately, this process comes at a large computational cost because the optimization problem must be solved and differentiated through in each training iteration; furthermore, it may also sometimes fail to improve solution quality due to non-smoothness issues that arise when training through a complex optimization layer. To address these shortcomings, we learn a low-dimensional surrogate model of a large optimization problem by representing the feasible space in terms of meta-variables, each of which is a linear combination of the original variables. By training a low-dimensional surrogate model end-to-end, and jointly with the predictive model, we achieve: i) a large reduction in training and inference time; and ii) improved performance by focusing attention on the more important variables in the optimization and learning in a smoother space. Empirically, we demonstrate these improvements on a non-convex adversary modeling task, a submodular recommendation task and a convex portfolio optimization task.


title: Inequalities from Lorentz-Finsler norms

abstract: We show that Lorentz-Finsler geometry offers a powerful tool in obtaining inequalities. With this aim, we first point out that a series of famous inequalities such as: the (weighted) arithmetic-geometric mean inequality, Acz\'el's, Popoviciu's and Bellman's inequalities, are all particular cases of a reverse Cauchy-Schwarz, respectively, of a reverse triangle inequality holding in Lorentz-Finsler geometry. Then, we use the same method to prove some completely new inequalities, including two refinements of Acz\'el's inequality.


title: Fast, Lifetime-Preserving Readout for High-Coherence Quantum Annealers

abstract: We demonstrate, for the first time, that a quantum flux parametron (QFP) is capable of acting as both isolator and amplifier in the readout circuit of a capacitively shunted flux qubit (CSFQ). By treating the QFP like a tunable coupler and biasing it such that the coupling is off, we show that $T_1$ of the CSFQ is not impacted by Purcell loss from its low-Q readout resonator ($Q_e = 760$) despite being detuned by only $40$ MHz. When annealed, the QFP amplifies the qubit's persistent current signal such that it generates a flux qubit-state-dependent frequency shift of $85$ MHz in the readout resonator, which is over $9$ times its linewidth. The device is shown to read out a flux qubit in the persistent current basis with fidelities surpassing $98.6\%$ with only $80$ ns integration, and reaches fidelities of $99.6\%$ when integrated for $1$ $\mu$s. This combination of speed and isolation is critical to the readout of high-coherence quantum annealers.


title: Kaczmarz-type inner-iteration preconditioned flexible GMRES methods for consistent linear systems

abstract: We propose using greedy and randomized Kaczmarz inner-iterations as preconditioners for the right-preconditioned flexible GMRES method to solve consistent linear systems, with a parameter tuning strategy for adjusting the number of inner iterations and the relaxation parameter. We also present theoretical justifications of the right-preconditioned flexible GMRES for solving consistent linear systems. Numerical experiments on overdetermined and underdetermined linear systems show that the proposed method is superior to the GMRES method preconditioned by NE-SOR inner iterations in terms of total CPU time.


title: Central Limit Theorem for Symmetric Exchangeable Random Variables

abstract: A central limit theorem for arrays of symmetric row-wise exchangeable random variables is presented. The result is valid for finite and infinite extendable and non-extendable sequences. Unlike most reported versions of the central limit theorem valid only for partial sums of finite non-extendable sequences, ours applies to the entire sum and mimics the classical theorem in this sense.


title: A New One-Point Residual-Feedback Oracle For Black-Box Learning and Control

abstract: Zeroth-order optimization (ZO) algorithms have been recently used to solve black-box or simulation-based learning and control problems, where the gradient of the objective function cannot be easily computed but can be approximated using the objective function values. Many existing ZO algorithms adopt two-point feedback schemes due to their fast convergence rate compared to one-point feedback schemes. However, two-point schemes require two evaluations of the objective function at each iteration, which can be impractical in applications where the data are not all available a priori, e.g., in online optimization. In this paper, we propose a novel one-point feedback scheme that queries the function value once at each iteration and estimates the gradient using the residual between two consecutive points. When optimizing a deterministic Lipschitz function, we show that the query complexity of ZO with the proposed one-point residual feedback matches that of ZO with the existing two-point schemes. Moreover, the query complexity of the proposed algorithm can be improved when the objective function has Lipschitz gradient. Then, for stochastic bandit optimization problems where only noisy objective function values are given, we show that ZO with one-point residual feedback achieves the same convergence rate as that of two-point scheme with uncontrollable data samples. We demonstrate the effectiveness of the proposed one-point residual feedback via extensive numerical experiments.


title: Pure Point Diffraction and Mean, Besicovitch and Weyl Almost Periodicity

abstract: We show that a translation bounded measure has pure point diffraction if and only if it is mean almost periodic. We then go on and show that a translation bounded measure solves what we call the phase problem if and only if it is Besicovitch almost periodic. Finally, we show that a translation bounded measure solves the phase problem independent of the underlying van Hove sequence if and only if it is Weyl almost periodic. These results solve fundamental issues in the theory of pure point diffraction and answer questions of Lagarias.


title: Cooperative Multi-Agent Reinforcement Learning with Partial Observations

abstract: In this paper, we propose a distributed zeroth-order policy optimization method for Multi-Agent Reinforcement Learning (MARL). Existing MARL algorithms often assume that every agent can observe the states and actions of all the other agents in the network. This can be impractical in large-scale problems, where sharing the state and action information with multi-hop neighbors may incur significant communication overhead. The advantage of the proposed zeroth-order policy optimization method is that it allows the agents to compute the local policy gradients needed to update their local policy functions using local estimates of the global accumulated rewards that depend on partial state and action information only and can be obtained using consensus. Specifically, to calculate the local policy gradients, we develop a new distributed zeroth-order policy gradient estimator that relies on one-point residual-feedback which, compared to existing zeroth-order estimators that also rely on one-point feedback, significantly reduces the variance of the policy gradient estimates improving, in this way, the learning performance. We show that the proposed distributed zeroth-order policy optimization method with constant stepsize converges to the neighborhood of a policy that is a stationary point of the global objective function. The size of this neighborhood depends on the agents' learning rates, the exploration parameters, and the number of consensus steps used to calculate the local estimates of the global accumulated rewards. Moreover, we provide numerical experiments that demonstrate that our new zeroth-order policy gradient estimator is more sample-efficient compared to other existing one-point estimators.


title: "And then they died": Using Action Sequences for Data Driven,Context Aware Gameplay Analysis

abstract: Many successful games rely heavily on data analytics to understand players and inform design. Popular methodologies focus on machine learning and statistical analysis of aggregated data. While effective in extracting information regarding player action, much of the context regarding when and how those actions occurred is lost. Qualitative methods allow researchers to examine context and derive meaningful explanations about the goals and motivations behind player behavior, but are difficult to scale. In this paper, we build on previous work by combining two existing methodologies: Interactive Behavior Analytics (IBA) and sequence analysis (SA), in order to create a novel, mixed methods, human-in-the-loop data analysis methodology that uses behavioral labels and visualizations to allow analysts to examine player behavior in a way that is context sensitive, scalable, and generalizable. We present the methodology along with a case study demonstrating how it can be used to analyze behavioral patterns of teamwork in the popular multiplayer game Defense of the Ancients 2 (DotA 2).


title: Simple mechanism that breaks the Hall effect linearity at low temperatures

abstract: Hall resistance $R_{\textrm{xy}}$ is commonly suggested to be linear-in-magnetic-field $B$, provided the field is small. We argue here that at low temperatures this linearity is broken due to weak localization/antilocalization phenomena in inhomogeneous systems, while in a uniform medium these effects do not affect the linear field dependence of $R_{\textrm{xy}}(B)$. We calculate the Hall resistance for different two-component media using a mean-field approach and show that this non-linearity is experimentally observable.


title: Pure point spectrum for dynamical systems and mean almost periodicity

abstract: We consider metrizable ergodic topological dynamical systems over locally compact, $\sigma$-compact abelian groups. We study pure point spectrum via suitable notions of almost periodicity for the points of the dynamical system. More specifically, we characterize pure point spectrum via mean almost periodicity of generic points. We then go on and show how Besicovitch almost periodic points determine both eigenfunctions and the measure in this case. After this, we characterize those systems arising from Weyl almost periodic points and use this to characterize weak and Bohr almost periodic systems. Finally, we consider applications to aperiodic order.


title: Tilted Halved Hexagons: Hexagons, Semi-hexagons, and Halved Hexagons Under One Roof

abstract: We investigate a new family of regions that is the universal generalization of three well-known region families in the field of enumeration of tilings: the quasi-regular hexagons, the semi-hexagons, and the halved hexagons. We prove a simple product formula for the number of tilings of these new regions. Our main result also yields the enumerations of two special classes of plane partitions with restricted parts.


title: When Painlev\'e-Gullstrand coordinates fail

abstract: Painlev\'e-Gullstrand coordinates, a very useful tool in spherical horizon thermodynamics, fail in anti-de Sitter space and in the inner region of Reissner-Nordstr\"om. We predict this breakdown to occur in any region containing negative Misner-Sharp-Hernandez quasilocal mass because of repulsive gravity stopping the motion of PG observers in radial free fall with zero initial velocity. PG coordinates break down also in the static Einstein universe for completely different reasons. The more general Martel-Poisson family of charts, which normally has PG coordinates as a limit, is reported for static cosmologies (de Sitter, anti-de Sitter and the static Einstein universe).


title: Stability of Internal States in Recurrent Neural Networks Trained on Regular Languages

abstract: We provide an empirical study of the stability of recurrent neural networks trained to recognize regular languages. When a small amount of noise is introduced into the activation function, the neurons in the recurrent layer tend to saturate in order to compensate the variability. In this saturated regime, analysis of the network activation shows a set of clusters that resemble discrete states in a finite state machine. We show that transitions between these states in response to input symbols are deterministic and stable. The networks display a stable behavior for arbitrarily long strings, and when random perturbations are applied to any of the states, they are able to recover and their evolution converges to the original clusters. This observation reinforces the interpretation of the networks as finite automata, with neurons or groups of neurons coding specific and meaningful input patterns.


title: Matrix Completion with Quantified Uncertainty through Low Rank Gaussian Copula

abstract: Modern large scale datasets are often plagued with missing entries. For tabular data with missing values, a flurry of imputation algorithms solve for a complete matrix which minimizes some penalized reconstruction error. However, almost none of them can estimate the uncertainty of its imputations. This paper proposes a probabilistic and scalable framework for missing value imputation with quantified uncertainty. Our model, the Low Rank Gaussian Copula, augments a standard probabilistic model, Probabilistic Principal Component Analysis, with marginal transformations for each column that allow the model to better match the distribution of the data. It naturally handles Boolean, ordinal, and real-valued observations and quantifies the uncertainty in each imputation. The time required to fit the model scales linearly with the number of rows and the number of columns in the dataset. Empirical results show the method yields state-of-the-art imputation accuracy across a wide range of data types, including those with high rank. Our uncertainty measure predicts imputation error well: entries with lower uncertainty do have lower imputation error (on average). Moreover, for real-valued data, the resulting confidence intervals are well-calibrated.


title: A spectrally accurate method for the dielectric obstacle scattering problem and applications to the inverse problem

abstract: We analyze the inverse problem to reconstruct the shape of a three dimensional homogeneous dielectric obstacle from the knowledge of noisy far field data. The forward problem is solved by a system of second kind boundary integral equations. For the numerical solution of these coupled integral equations we propose a fast spectral algorithm by transporting these equations onto the unit sphere. We review the differentiability properties of the boundary to far field operator and give a characterization of the adjoint operator of the first Fr\'echet derivative. Using these results we discuss the implementation of the iteratively regularized Gauss-Newton method for the numerical solution of the inverse problem and give numerical results for star-shaped obstacles.


title: A Methodology for Assessing the Environmental Effects Induced by ICT Services. Part I: Single Services

abstract: Information and communication technologies (ICT) are increasingly seen as key enablers for climate change mitigation measures. They can make existing products and activities more efficient or substitute them altogether. Consequently, different initiatives have started to estimate the environmental effects of ICT services. Such assessments, however, lack scientific rigor and often rely on crude assumptions and methods, leading to inaccurate or even misleading results. The few methodological attempts that exist do not address several crucial aspects, and are thus insufficient to foster good as-sessment practice. Starting from such a high level standard from the European Telecommunication Standardisation Institute (ETSI) and the International Telecommunication Union (ITU), this article identifies the shortcomings of existing methodologies and proposes solutions. It addresses several aspects for the assessment of single ICT services: the goal and scope definition (analyzing differences between ICT substitution and optimization, the time perspective of the assessment, the challenge of a hypothetical baseline for the situation without the ICT solution, and the differences between modelling and case studies) as well as the often ignored influence of rebound effects and the difficult extrapolation from case studies to larger populations.


title: A Global Non-Hydrostatic Atmospheric Model with a Mass and Energy Conserving Vertically-Implicit-Correction (VIC) Scheme

abstract: Global non-hydrostatic atmospheric models are becoming increasingly important for studying the climates of planets and exoplanets. However, such models suffer from computational difficulties due to the large aspect ratio between the horizontal and vertical directions. To overcome this problem, we developed a global model using a vertically-implicit-correction (VIC) scheme in which the integration time step is no longer limited by the propagation of acoustic waves in the vertical. We proved that our model, based on the $\rm Athena^{++}$ framework and its extension for planetary atmospheres - SNAP (Simulating Non-hydrostatic Atmosphere on Planets), rigorously conserves mass and energy in finite volume simulations. We found that traditional numerical stabilizers such as hyper-viscosity and divergence damping are not needed when using the VIC scheme, which greatly simplifies the numerical implementation and improves stability. We present simulation results ranging from 1D linear waves to 3D global circulations with and without the VIC scheme. These tests demonstrate that our formulation correctly tracks local turbulent motions, produces Kelvin-Helmholtz instability, and generates a super-rotating jet on hot Jupiters. Employing this VIC scheme improves the computational efficiency of global simulations by more than two orders of magnitude compared to an explicit model and facilitates the capability of simulating a wide range of planetary atmospheres both regionally and globally.


title: Amortized Causal Discovery: Learning to Infer Causal Graphs from Time-Series Data

abstract: On time-series data, most causal discovery methods fit a new model whenever they encounter samples from a new underlying causal graph. However, these samples often share relevant information which is lost when following this approach. Specifically, different samples may share the dynamics which describe the effects of their causal relations. We propose Amortized Causal Discovery, a novel framework that leverages such shared dynamics to learn to infer causal relations from time-series data. This enables us to train a single, amortized model that infers causal relations across samples with different underlying causal graphs, and thus leverages the shared dynamics information. We demonstrate experimentally that this approach, implemented as a variational model, leads to significant improvements in causal discovery performance, and show how it can be extended to perform well under added noise and hidden confounding.


title: Schr\"odinger evolution of the Hawking state

abstract: A Schr\"odinger-picture description of the evolving quantum state of Hawking radiation is given, based on an ADM decomposition using time slicings that smoothly cross the horizon. This treatment avoids requiring a role for trans-planckian modes, which can be viewed as artifacts of Hawking's original calculation, and also supports arguments that radiation from black holes is produced in a "quantum atmosphere" with thickness comparable to the horizon size, rather than microscopically far from it. Particularly explicit formulas are given for the two-dimensional analog of the Schwarzschild geometry. This analysis is expected to generalize to other black holes, and to cosmology. The resulting quantum evolution also provides important background for investigating corrections to the Hawking process, as are necessary for restoring unitary evolution of black hole decay.


title: Bounded confidence dynamics and graph control: enforcing consensus

abstract: A generic feature of bounded confidence type models is the formation of clusters of agents. We propose and study a variant of bounded confidence dynamics with the goal of inducing unconditional convergence to a consensus. The defining feature of these dynamics, which we name the No one left behind dynamics, is the introduction of a local control on the agents which preserves the connectivity of the interaction network. We rigorously demonstrate that these dynamics result in unconditional convergence to a consensus. The qualitative nature of our argument prevents us quantifying how fast a consensus emerges, however we present numerical evidence that sharp convergence rates would be challenging to obtain for such dynamics. Finally, we propose a relaxed version of the control. The dynamics that result maintain many of the qualitative features of the bounded confidence dynamics yet ultimately still converge to a consensus as the control still maintains connectivity of the interaction network.


title: An Integer Linear Programming Framework for Mining Constraints from Data

abstract: Structured output prediction problems (e.g., sequential tagging, hierarchical multi-class classification) often involve constraints over the output label space. These constraints interact with the learned models to filter infeasible solutions and facilitate in building an accountable system. However, although constraints are useful, they are often based on hand-crafted rules. This raises a question -- \emph{can we mine constraints and rules from data based on a learning algorithm?} In this paper, we present a general framework for mining constraints from data. In particular, we consider the inference in structured output prediction as an integer linear programming (ILP) problem. Then, given the coefficients of the objective function and the corresponding solution, we mine the underlying constraints by estimating the outer and inner polytopes of the feasible set. We verify the proposed constraint mining algorithm in various synthetic and real-world applications and demonstrate that the proposed approach successfully identifies the feasible set at scale. In particular, we show that our approach can learn to solve 9x9 Sudoku puzzles and minimal spanning tree problems from examples without providing the underlying rules. Our algorithm can also integrate with a neural network model to learn the hierarchical label structure of a multi-label classification task. Besides, we provide a theoretical analysis about the tightness of the polytopes and the reliability of the mined constraints.