title: Emitted Radiation and Geometry

abstract: In conformal $\mathcal{N}=2$ Super Yang-Mills theory, the energy emitted by an accelerated heavy particle is computed by the one-point function of the stress tensor operator in the presence of a Wilson line. In this paper, we consider the theory on the ellipsoid and we prove a conjectured relation between the stress tensor one-point function and the first order expansion of the Wilson loop expectation value in the squashing parameter. To do this, we analyze the behavior of the Wilson loop for a small deformation of the background geometry and, at first order in the deformation, we fix the kinematics using defect CFT constraints. In the final part of the paper, we analyze the consequences of our results for the weak coupling perturbative expansion. In particular, comparing the weakly coupled matrix model with the ordinary Feynman diagram expansion, we find a natural transcendentality driven organization for the latter.


title: Evidence for reduced magnetic braking in polars from binary population models

abstract: We present the first population synthesis of synchronous magnetic cataclysmic variables, called polars, taking into account the effect of the white dwarf (WD) magnetic field on angular momentum loss. We implemented the reduced magnetic braking (MB) model proposed by Li, Wu & Wickramasinghe into the Binary Stellar Evolution (BSE) code recently calibrated for cataclysmic variable (CV) evolution. We then compared separately our predictions for polars and non-magnetic CVs with a large and homogeneous sample of observed CVs from the Sloan Digital Sky Survey. We found that the predicted orbital period distributions and space densities agree with the observations if period bouncers are excluded. For polars, we also find agreement between predicted and observed mass transfer rates, while the mass transfer rates of non-magnetic CVs with periods ${\gtrsim3}$ hr drastically disagree with those derived from observations. Our results provide strong evidence that the reduced MB model for the evolution of highly magnetized accreting WDs can explain the observed properties of polars. The remaining main issues in our understanding of CV evolution are the origin of the large number of highly magnetic WDs, the large scatter of the observed mass transfer rates for non-magnetic systems with periods ${\gtrsim3}$ hr, and the absence of period bouncers in observed samples.

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


title: Cold neutral hydrogen gas in galaxies

abstract: This review summarizes recent studies of the cold neutral hydrogen gas associated with galaxies probed via the HI 21-cm absorption line. HI 21-cm absorption against background radio-loud quasars is a powerful tool to study the neutral gas distribution and kinematics in foreground galaxies from kilo-parsec to parsec scales. At low redshifts (z<0.4), it has been used to characterize the distribution of high column density neutral gas around galaxies and study the connection of this gas with the galaxy's optical properties. The neutral gas around galaxies has been found to be patchy in distribution, with variations in optical depth observed at both kilo-parsec and parsec scales. At high redshifts (z>0.5), HI 21-cm absorption has been used to study the neutral gas in metal or Lyman-alpha absorption-selected galaxies. It has been found to be closely linked with the metal and dust content of the gas. Trends of various properties like incidence, spin temperature and velocity width of HI 21-cm absorption with redshift have been studied, which imply evolution of cold gas properties in galaxies with cosmic time. Upcoming large blind surveys of HI 21-cm absorption with next generation radio telescopes are expected to determine accurately the redshift evolution of the number density of HI 21-cm absorbers per unit redshift and hence understand what drives the global star formation rate density evolution.


title: Hic sunt dracones: Cartography of the Milky Way spiral arms and bar resonances with Gaia Data Release 2

abstract: In this paper we introduce a new method for analysing Milky Way phase-space which allows us to reveal the imprint left by the Milky Way bar and spiral arms on the stars with full phase-space data in Gaia Data Release 2. The unprecedented quality and extended spatial coverage of these data enable us to discover six prominent stellar density structures in the disc to a distance of 5 kpc from the Sun. Four of these structures correspond to the spiral arms detected previously in the gas and young stars (Scutum-Centaurus, Sagittarius, Local and Perseus). The remaining two are associated with the main resonances of the Milky Way bar where corotation is placed at around 6.2 kpc and the outer Lindblad resonance beyond the Solar radius, at around 9 kpc. For the first time we provide evidence of the imprint left by spiral arms and resonances in the stellar densities not relying on a specific tracer, through enhancing the signatures left by these asymmetries. Our method offers new avenues for studying how the stellar populations in our Galaxy are shaped.


title: Systematic Constructions of Fracton Theories

abstract: Fracton theories possess exponentially degenerate ground states, excitations with restricted mobility, and nontopological higher-form symmetries. This paper shows that such theories can be defined on arbitrary spatial lattices in three dimensions. The key element of this construction is a generalization of higher-form gauge theories to so-called $\mathfrak{F}_p$ gauge theories, in which gauge transformations of rank-$k$ fields are specified by rank-$(k - p)$ gauge parameters. The $\mathbb{Z}_2$ rank-two theory of type $\mathfrak{F}_2$, placed on a cubic lattice and coupled to scalar matter, is shown to have a topological phase exactly dual to the well-known X-cube model. Generalizations of this example yield novel fracton theories. In the continuum, the $\mathrm{U}(1)$ rank-two theory of type $\mathfrak{F}_2$ is shown to have a perturbatively gapless fracton regime that cannot be consistently interpreted as a tensor gauge theory of any kind. The compact scalar fields that naturally couple to this $\mathfrak{F}_2$ theory also show gapless fracton behavior; on a cubic lattice they have a conserved $\mathrm{U}(1)$ charge and dipole moment, but these particular charges are not necessarily conserved on more general lattices. The construction straightforwardly generalizes to $\mathfrak{F}_2$ theories of nonabelian rank-two gauge fields, giving first examples of pure nonabelian higher-rank theories.


title: The Pristine Inner Galaxy Survey (PIGS) I: Tracing the kinematics of metal-poor stars in the Galactic bulge

abstract: Our Galaxy is known to contain a central boxy/peanut-shaped bulge, yet the importance of a classical, pressure-supported component within the central part of the Milky Way is still being debated. It should be most visible at low metallicity, a regime that has not yet been studied in detail. Using metallicity-sensitive narrow-band photometry, the Pristine Inner Galaxy Survey (PIGS) has collected a large sample of metal-poor ([Fe/H] < -1.0) stars in the inner Galaxy to address this open question. We use PIGS to trace the metal-poor inner Galaxy kinematics as function of metallicity for the first time. We find that the rotational signal decreases with decreasing [Fe/H], until it becomes negligible for the most metal-poor stars. Additionally, the velocity dispersion increases with decreasing metallicity for -3.0 < [Fe/H] < -0.5, with a gradient of -44 $\pm$ 4 km$\,$s$^{-1}\,$dex$^{-1}$. These observations may signal a transition between Galactic components of different metallicities and kinematics, a different mapping onto the boxy/peanut-shaped bulge for former disk stars of different metallicities and/or the secular dynamical and gravitational influence of the bar on the pressure-supported component. Our results provide strong constraints on models that attempt to explain the properties of the inner Galaxy.


title: Capture of interstellar objects: a source of long-period comets

abstract: We simulate the passage through the Sun-Jupiter system of interstellar objects (ISOs) similar to 1I/`Oumuamua or 2I/Borisov. Capture of such objects is rare and overwhelmingly from low incoming speeds onto orbits akin to those of known long-period comets. This suggests that some of these comets could be of extra-solar origin, in particular inactive ones. Assuming ISOs follow the local stellar velocity distribution, we infer a volume capture rate of $0.051\,\mathrm{au}^3 \mathrm{yr}^{-1}$. Current estimates for orbital lifetimes and space densities then imply steady-state captured populations of $\sim10^2$ comets and $\sim10^5$ `Oumuamua-like rocks, of which 0.033% are within 6au at any time.


title: On modelling CC85 wind in an expanding local box

abstract: The Chevalier & Clegg (1985) spherical wind is widely used to model galactic outflows. Efforts are underway to understand the production of multiphase gas in such outflows. Since the important scales for the multiphase gas are very small compared to the size of the wind, it is very useful to have a local model of the cold gas interacting with a volume-filling expanding wind. In this note we explicitly derive the governing hydrodynamic equations in such a local box for the first time. We highlight that the expansion in a wind occurs only in the transverse direction, and this results in anisotropic source terms in the momentum equation. Our equations will be very useful in realistic modelling of the small scale structure of multiphase galactic outflows. [original note does not include an abstract]


title: Pristine IX: CFHT ESPaDOnS Spectroscopic Analysis of 115 Bright Metal-Poor Candidate Stars

abstract: A chemo-dynamical analysis of 115 metal-poor candidate stars selected from the narrow-band Pristine photometric survey is presented based on CFHT high-resolution ESPaDoNS spectroscopy. We have discover 28 new bright (V < 15) stars with [Fe/H]<-2.5 and 5 with [Fe/H]<-3.0 for success rates of 40% (28/70) and 19% (5/27), respectively. A detailed model atmospheres analysis is carried out for the 28 new metal-poor stars. Stellar parameters were determined from SDSS photometric colours, Gaia DR2 parallaxes, MESA/MIST stellar isochrones, and the initial Pristine survey metallicities, following a Bayesian inference method. Chemical abundances are determined for 10 elements (Na, Mg, Ca, Sc, Ti, Cr, Fe, Ni, Y, Ba). Most stars show chemical abundance patterns that are similar to the normal metal-poor stars in the Galactic halo; however, we also report the discoveries of a new r-process rich star, a new CEMP-s candidate with [Y/Ba]>0, and a metal-poor star with very low [Mg/Fe]. The kinematics and orbits for all of the highly probable metal-poor candidates are determined by combining our precision radial velocities with Gaia DR2 proper motions. Some stars show unusual kinematics for their chemistries, including planar orbits, unbound orbits, and highly elliptical orbits that plunge deeply into the Galactic bulge (Rperi < 0.5 kpc); also, eight stars have orbital energies and actions consistent with the Gaia-Enceladus accretion event. This paper contributes to our understanding of the complex chemo-dynamics of the metal-poor Galaxy, and increases the number of known bright metal-poor stars available for detailed nucleosynthetic studies.


title: Statistical localization: from strong fragmentation to strong edge modes

abstract: Certain disorder-free Hamiltonians can be non-ergodic due to a \emph{strong fragmentation} of the Hilbert space into disconnected sectors. Here, we characterize such systems by introducing the notion of `statistically localized integrals of motion' (SLIOM), whose eigenvalues label the connected components of the Hilbert space. SLIOMs are not spatially localized in the operator sense, but appear localized to sub-extensive regions when their expectation value is taken in typical states with a finite density of particles. We illustrate this general concept on several Hamiltonians, both with and without dipole conservation. Furthermore, we demonstrate that there exist perturbations which destroy these integrals of motion in the bulk of the system, while keeping them on the boundary. This results in statistically localized \emph{strong zero modes}, leading to infinitely long-lived edge magnetizations along with a thermalizing bulk, constituting the first example of such strong edge modes in a non-integrable model. We also show that in a particular example, these edge modes lead to the appearance of topological string order in a certain subset of highly excited eigenstates. Some of our suggested models can be realized in Rydberg quantum simulators.


title: A search for optical AGN variability in 35,000 low-mass galaxies with the Palomar Transient Factory

abstract: We present an analysis of the long-term optical variability for $\sim50,000$ nearby (z<0.055) galaxies from the NASA-Sloan Atlas, $35,000$ of which are low-mass ($M_{\ast}<10^{10}~M_{\odot}$). We use difference imaging of Palomar Transient Factory (PTF) R-band observations to construct light curves with typical baselines of several years. We then search for subtle variations in the nuclear light output. We determine whether detected variability is AGN-like by assessing the fit quality to a damped random walk model. We identify 424 variability-selected AGN, including 244 with stellar masses between $10^{7}$ and $10^{10}~M_{\odot}$. 75% of low-mass galaxies with AGN-like variability have narrow emission lines dominated by star formation. After controlling for nucleus magnitude, the fraction of variable AGN is constant down to $M_{\ast}=10^{9}~M_{\odot}$, suggesting no drastic decline in the BH occupation fraction down to this stellar mass regime. Combining our NASA-Sloan Atlas sample with samples of nearby galaxies with broad H$\alpha$ emission, we find no dependence of variability properties with black hole mass. However, we caution that the variable AGN fraction is strongly dependent on baseline. For baselines less than two years, the variable fraction for the full sample is 0.25%, compared to 1.0% for baselines longer than two years. Finally, comparing Stripe 82 light curves (Baldassare et al. 2018) to PTF light curves, we find populations of changing-look AGN: 8 galaxies that are variable in Stripe 82, but quiescent in PTF, and 15 galaxies where the reverse is true. Our PTF work demonstrates the promise of long-term optical variability searches in low-mass galaxies for finding AGNs missed by other selection techniques.


title: Quantum certification and benchmarking

abstract: Concomitant with the rapid development of quantum technologies, challenging demands arise concerning the certification and characterization of devices. The promises of the field can only be achieved if stringent levels of precision of components can be reached and their functioning guaranteed. This review provides a brief overview of the known characterization methods of certification, benchmarking, and tomographic recovery of quantum states and processes, as well as their applications in quantum computing, simulation, and communication.


title: Analytic solutions of neutral hyperbolic black holes with scalar hair

abstract: We find analytic solutions of hyperbolic black holes with scalar hair in anti-de Sitter (AdS) space, and they do not have spherical or planar counterparts. The system is obtained by taking a neutral limit of an Einstein-Maxwell-dilaton system whose special cases are maximal gauged supergravities, while the dilaton is kept nontrivial. There are phase transitions between these black holes and the hyperbolic Schwarzschild-AdS black hole. We discuss two AdS/CFT applications of these hyperbolic black holes. One is phase transitions of holographic Renyi entropies, and the other is phase transitions of quantum field theories in de Sitter space. In addition, we give a C-metric solution as a generalization of the hyperbolic black holes with scalar hair.


title: The universal acceleration scale from stellar feedback

abstract: It has been established for decades that rotation curves deviate from the Newtonian gravity expectation given baryons alone below a characteristic acceleration scale $g_{\dagger}\sim 10^{-8}\,\rm{cm\,s^{-2}}$, a scale promoted to a new fundamental constant in MOND. In recent years, theoretical and observational studies have shown that the star formation efficiency (SFE) of dense gas scales with surface density, SFE $\sim \Sigma/\Sigma_{\rm crit}$ with $\Sigma_{\rm crit} \sim \langle\dot{p}/m_{\ast}\rangle/(\pi\,G)\sim 1000\,\rm{M_{\odot}\,pc^{-2}}$ (where $\langle \dot{p}/m_{\ast}\rangle$ is the momentum flux output by stellar feedback per unit stellar mass in a young stellar population). We argue that the SFE, more generally, should scale with the local gravitational acceleration, i.e. that SFE $\sim g_{\rm tot}g_\mathrm{crit} \equiv (G\,M_{\rm tot}/R^{2}) / \langle\dot{p}/m_{\ast}\rangle$, where $M_{\rm tot}$ is the total gravitating mass and $g_\mathrm{crit}=\langle\dot{p}/m_{\ast}\rangle = \pi\,G\,\Sigma_{\rm crit} \approx 10^{-8}\,\rm{cm\,s^{-2}} \approx g_{\dagger}$. Hence the observed $g_\dagger$ may correspond to the characteristic acceleration scale above which stellar feedback cannot prevent efficient star formation, and baryons will eventually come to dominate. We further show how this may give rise to the observed acceleration scaling $g_{\rm obs}\sim(g_{\rm baryon}\,g_{\dagger})^{1/2}$ (where $g_{\rm baryon}$ is the acceleration due to baryons alone) and flat rotation curves. The derived characteristic acceleration $g_{\dagger}$ can be expressed in terms of fundamental constants (gravitational constant, proton mass, and Thomson cross section): $g_{\dagger}\sim 0.1\,G\,m_{p}/\sigma_{\rm T}$.


title: Revealing the Origin and Cosmic Evolution of Supermassive Black Holes

abstract: The next generation of electromagnetic and gravitational wave observatories will open unprecedented windows to the birth of the first supermassive black holes. This has the potential to reveal their origin and growth in the first billion years, as well as the signatures of their formation history in the local Universe. With this in mind, we outline three key focus areas which will shape research in the next decade and beyond: (1) What were the "seeds" of the first quasars; how did some reach a billion solar masses before z$\sim7$? (2) How does black hole growth change over cosmic time, and how did the early growth of black holes shape their host galaxies? What can we learn from intermediate mass black holes (IMBHs) and dwarf galaxies today? (3) Can we unravel the physics of black hole accretion, understanding both inflows and outflows (jets and winds) in the context of the theory of general relativity? Is it valid to use these insights to scale between stellar and supermassive BHs, i.e., is black hole accretion really scale invariant? In the following, we identify opportunities for the Canadian astronomical community to play a leading role in addressing these issues, in particular by leveraging our strong involvement in the Event Horizon Telescope, the {\it James Webb Space Telescope} (JWST), Euclid, the Maunakea Spectroscopic Explorer (MSE), the Thirty Meter Telescope (TMT), the Square Kilometer Array (SKA), the Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR), and more. We also discuss synergies with future space-based gravitational wave (LISA) and X-ray (e.g., Athena, Lynx) observatories, as well as the necessity for collaboration with the stellar and galactic evolution communities to build a complete picture of the birth of supermassive black holes, and their growth and their influence over the history of the Universe.


title: Polarised Amplitudes and Soft-Virtual Cross Sections for $b\bar b \rightarrow ZH$ at NNLO in QCD

abstract: Production of the Higgs boson, $H$ in association with a massive vector boson, $V$, i.e., the $VH$ process, plays an important role in the explorations of Higgs physics at the Large Hadron Collider, both for a precise study of Higgs' Standard Model couplings and for probing New Physics. In this publication we present the two-loop corrections in massless quantum chromodynamics (QCD) to the amplitude of the Higgs production associated with a $Z$ boson via the bottom quark-antiquark annihilation channel with a non-vanishing bottom-quark Yukawa coupling, which is a necessary ingredient of the full next-to-next-to-leading-order QCD corrections to the $VH$ process in the five-flavour scheme. The computation is performed by projecting the D-dimensional scattering amplitude directly onto an appropriate set of Lorentz structures related to the linear polarisation states of the $Z$ boson. We provide analytic expressions of the complete set of renormalised polarised amplitudes in terms of polylogarithms of maximum weight four. To give an estimation of the size of contributions from amplitudes considered in this work, we compute numerically the resulting cross sections under the soft-virtual approximation. We also take the opportunity to make a dedicated discussion regarding an interesting subtlety appearing in the conventional form factor decomposition of amplitudes involving axial currents regularised in D dimensions.


title: Brane nucleation instabilities in non-AdS/non-CFT

abstract: We speculate that the weak gravity conjecture applied to theories with holographic duals bans the existence of disordered phases at zero temperature. We test this idea by introducing a non-zero baryon chemical potential in a deformation of the $SU(N_c)\times SU(N_c)$ Klebanov-Witten gauge theory with broken supersymmetry and conformal invariance. At low temperature, a disordered phase dual to a black brane geometry is unstable for low chemical potentials and metastable for high values. In the metastable phase, states with a partial Higgsing of the gauge group are favored over the normal disordered phase. This is reflected in the properties of the effective potential for color branes in the dual geometry, where the appearance of a global minimum outside the horizon signals the onset of a brane nucleation instability. When the Higgsing involves only one of the group factors, the global minimum remains at a finite distance from the horizon, making it possible to construct holographic duals to metastable "color superconducting" states. We also consider branes dual to excitations with baryon charge, but find that the extremal geometry remains marginally stable against the emission of particles carrying baryon charge independently of the strength of the deformation.


title: A general-purpose timestep criterion for simulations with gravity

abstract: We describe a new adaptive timestep criterion for integrating gravitational motion, which uses the tidal tensor to estimate the local dynamical timescale and scales the timestep proportionally. This provides a better candidate for a truly general-purpose gravitational timestep criterion than the usual prescription derived from the gravitational acceleration, which does not respect the equivalence principle, breaks down when $\mathbf{a}=0$, and does not obey the same dimensional scaling as the true timescale of orbital motion. We implement the tidal timestep criterion in the simulation code GIZMO, and examine controlled tests of collisionless galaxy and star cluster models, as well as fully-dynamic galaxy merger and cosmological dark matter simulations. The tidal criterion estimates the dynamical time faithfully, and generally provides a more efficient timestepping scheme compared to an acceleration criterion. Specifically, the tidal criterion achieves order-of-magnitude smaller energy errors for the same number of force evaluations in potentials with inner profiles shallower than $\rho \propto r^{-1}$ (ie. where $\mathbf{a}\rightarrow 0$), such as star clusters and cored galaxies. For a given problem these advantages must be weighed against the additional overhead of computing the tidal tensor on-the-fly, but in many cases this overhead is small.


title: Quantum gravity, minimum length, and holography

abstract: The Karolyhazy uncertainty relation is the statement that if a device is used to measure a length $l$, there will be a minimum uncertainty $\delta l$ in the measurement, given by $(\delta l)^3 \sim L_P^2\; l$. This is a consequence of combining the principles of quantum mechanics and general relativity. In this note we show how this relation arises in our approach to quantum gravity, in a bottom-up fashion, from the matrix dynamics of atoms of space-time-matter. We use this relation to define a space-time-matter foam at the Planck scale, and to argue that our theory is holographic. By coarse graining over time scales larger than Planck time, one obtains the laws of quantum gravity. Quantum gravity is not a Planck scale phenomenon; rather it comes into play whenever no classical space-time background is available to describe a quantum system. Space-time and classical general relativity arise from spontaneous localisation in a highly entangled quantum gravitational system. The Karolyhazy relation continues to hold in the emergent theory. An experimental confirmation of this relation will constitute a definitive test of the quantum nature of gravity.


title: Measurements of electric quadrupole transition frequencies in $^{226}$Ra$^{+}$

abstract: We report the first driving of the $7s\ ^2S_{1/2}\rightarrow 6d\ ^2D_{3/2}$ and $7s\ ^2S_{1/2}\rightarrow 6d\ ^2D_{5/2}$ electric quadrupole (E2) transitions in Ra$^{+}$. We measure the frequencies of both E2 transitions, and two other low-lying transitions in $^{226}$Ra$^{+}$ that are important for controlling the radium ion's motional and internal states: $6d\ ^2D_{3/2}\rightarrow 7p\ ^2P_{3/2}^{0}$ and $6d\ ^2D_{5/2}\rightarrow 7p\ ^2P_{3/2}^{0}$.


title: Athermal fluctuations in disordered crystals

abstract: We analyze the fluctuations in particle positions and inter-particle forces in disordered jammed crystals in the limit of weak disorder. We demonstrate that such athermal systems are fundamentally different from their thermal counterparts, characterized by constrained fluctuations of forces perpendicular to the lattice directions. We develop a disorder perturbation expansion in polydispersity about the crystalline state, which we use to derive exact results to linear order. We show that constrained fluctuations result as a consequence of local force balance conditions, and are characterized by non-Gaussian distributions which we derive exactly. We analytically predict several properties of such systems, including the scaling of the average coordination with polydispersity and packing fraction, which we verify with numerical simulations using soft disks with one-sided harmonic interactions.


title: Universal finite-size gap scaling of the quantum Ising chain

abstract: I study the universal finite-size scaling function for the lowest gap of the quantum Ising chain with a one-parameter family of ``defect'' boundary conditions, which includes periodic, open, and antiperiodic boundary conditions as special cases. The universal behavior can be described by the Majorana fermion field theory in $1+1$ dimensions, with the mass proportional to the deviation from the critical point. Although the field theory appears to be symmetric with respect to the inversion of the mass (Kramers-Wannier duality), the actual gap is asymmetric, reflecting the spontaneous symmetry breaking in the ordered phase which leads to the two-fold ground-state degeneracy in the thermodynamic limit. The asymptotic ground-state degeneracy in the ordered phase is realized by (i) formation of a bound state at the defect (except for the periodic/antiperiodic boundary condition) and (ii) effective reversal of the fermion number parity in one of the sectors (except for the open boundary condition), resulting in a rather nontrivial crossover ``phase'' diagram in the space of the boundary condition (defect strength) and mass.


title: Fermi-Hubbard model on non-bipartite lattices: flux problem and emergent chirality

abstract: On several one-dimensional (1D) and 2D nonbipartite lattices, we study both free and Hubbard interacting lattice fermions when some magnetic fluxes are threaded or gauge fields coupled. First, we focus on finding out the optimal flux which minimizes the energy of fermions at specific fillings. For spin-$1/2$ fermions at half-filling on a ring lattice consisting of odd-numbered sites, the optimal flux turns out to be $\pm{\pi}/2$. We prove this conclusion for Hubbard interacting fermions utilizing a generalized reflection positivity technique, which can lead to further applications on 2D nonbipartite lattices such as triangular and Kagome. At half-filling the optimal flux patterns on the triangular and Kagome lattice are ascertained to be $\pm[{\pi}/2,{\pi}/2]$, $\pm[{\pi}/2,{\pi}/2,0]$, respectively (see the meaning of these notations in the main text). We also find that chirality emerges in these optimal flux states. Then, we verify these exact conclusions and further study some other fillings with the numerical exact diagonalization method. It is found that when it deviates from half-filling, Hubbard interactions can alter the optimal flux patterns on these lattices. Moreover, numerically observed emergent flux singularities driven by strong Hubbard interactions in the ground states, both in 1D and 2D, are discussed and interpreted as some kind of non-Fermi liquid feature.


title: Topological Spin Liquids: Robustness under perturbations

abstract: We study the robustness of the paradigmatic kagome Resonating Valence Bond (RVB) spin liquid and its orthogonal version, the quantum dimer model. The non-orthogonality of singlets in the RVB model and the induced finite length scale not only makes it difficult to analyze, but can also significantly affect its physics, such as how much noise resilience it exhibits. Surprisingly, we find that this is not the case: The amount of perturbations which the RVB spin liquid can tolerate is not affected by the finite correlation length, making the dimer model a viable model for studying RVB physics under perturbations. Remarkably, we find that this is a universal phenomenon protected by symmetries: First, the dominant correlations in the RVB are spinon correlations, making the state robust against doping with visons. Second, reflection symmetry stabilizes the spin liquid against doping with spinons, by forbidding mixing of the initially dominant correlations with those which lead to the breakdown of topological order.


title: Implications of the Gaia Sausage for Dark Matter Nuclear Interactions

abstract: The advent of the Gaia era has led to potentially revolutionary understanding of dark matter (DM) dynamics in our galaxy, which has important consequences for direct detection (DD) experiments. In this paper, we study how the empirical DM velocity distribution inferred from Gaia-Sausage, a dominant substructure in the solar neighborhood, affects the interpretation of DD data. We survey different classes of operators in the non-relativistic effective field theory that could arise from several relativistic benchmark models and emphasize that the Gaia velocity distribution could modify both the total number of events as well as the shape of the differential recoil spectra, the two primary observables in DD experiments. Employing the euclideanized signal method, we investigate the effects of the Gaia distribution on reconstructing DM model parameters and identifying the correct DM model given a positive signal at future DD experiments. We find that for light DM with mass ~10 GeV, the Gaia distribution poses an additional challenge for characterizing DM interactions with ordinary matter, which may be addressed by combining complementary DD experiments with different targets and lowering the detection threshold. Meanwhile, for heavy DM with mass above ~30 GeV, depending on the type of DM model, there could be a (moderate) improvement in the sensitivity at future DD experiments.


title: Heavy-light Bootstrap from Lorentzian Inversion Formula

abstract: We study heavy-light four-point function by employing Lorentzian inversion formula, where the conformal dimension of heavy operator is as large as central charge $C_T\rightarrow\infty$. We implement the Lorentzian inversion formula back and forth to reveal the universality of the lowest-twist multi-stress-tensor $T^k$ as well as large spin double-twist operators $[\mathcal{O}_H\mathcal{O}_L]_{n',J'}$. In this way, we also propose an algorithm to bootstrap the heavy-light four-point function by extracting relevant OPE coefficients and anomalous dimensions. By following the algorithm, we exhibit the explicit results in $d=4$ up to the triple-stress-tensor. Moreover, general dimensional heavy-light bootstrap up to the double-stress-tensor is also discussed, and we present an infinite series representation of the lowest-twist double-stress-tensor OPE coefficient. Exact expressions of lowest-twist double-stress-tensor OPE coefficients in $d=6,8,10$ are also obtained as further examples.


title: Asymmetric Shapley values: incorporating causal knowledge into model-agnostic explainability

abstract: Explaining AI systems is fundamental both to the development of high performing models and to the trust placed in them by their users. The Shapley framework for explainability has strength in its general applicability combined with its precise, rigorous foundation: it provides a common, model-agnostic language for AI explainability and uniquely satisfies a set of intuitive mathematical axioms. However, Shapley values are too restrictive in one significant regard: they ignore all causal structure in the data. We introduce a less restrictive framework, Asymmetric Shapley values (ASVs), which are rigorously founded on a set of axioms, applicable to any AI system, and flexible enough to incorporate any causal structure known to be respected by the data. We demonstrate that ASVs can (i) improve model explanations by incorporating causal information, (ii) provide an unambiguous test for unfair discrimination in model predictions, (iii) enable sequentially incremental explanations in time-series models, and (iv) support feature-selection studies without the need for model retraining.


title: Resolved observations at 31 GHz of spinning dust emissivity variations in $\rho$ Oph

abstract: The $\rho$ Oph molecular cloud is one of the best examples of spinning dust emission, first detected by the Cosmic Background Imager (CBI). Here we present 4.5 arcmin observations with CBI 2 that confirm 31 GHz emission from $\rho$ Oph W, the PDR exposed to B-type star HD 147889, and highlight the absence of signal from S1, the brightest IR nebula in the complex. In order to quantify an association with dust-related emission mechanisms, we calculated correlations at different angular resolutions between the 31 GHz map and proxies for the column density of IR emitters, dust radiance and optical depth templates. We found that the 31 GHz emission correlates best with the PAH column density tracers, while the correlation with the dust radiance improves when considering emission that is more extended (from the shorter baselines), suggesting that the angular resolution of the observations affects the correlation results. A proxy for the spinning dust emissivity reveals large variations within the complex, with a dynamic range of 25 at 3$\sigma$ and a variation by a factor of at least 23, at 3$\sigma$, between the peak in $\rho$ Oph W and the location of S1, which means that environmental factors are responsible for boosting spinning dust emissivities locally.


title: Structured Pruning of a BERT-based Question Answering Model

abstract: The recent trend in industry-setting Natural Language Processing (NLP) research has been to operate large %scale pretrained language models like BERT under strict computational limits. While most model compression work has focused on "distilling" a general-purpose language representation using expensive pretraining distillation, less attention has been paid to creating smaller task-specific language representations which, arguably, are more useful in an industry setting. In this paper, we investigate compressing BERT- and RoBERTa-based question answering systems by structured pruning of parameters from the underlying transformer model. We find that an inexpensive combination of task-specific structured pruning and task-specific distillation, without the expense of pretraining distillation, yields highly-performing models across a range of speed/accuracy tradeoff operating points. We start from existing full-size models trained for SQuAD 2.0 or Natural Questions and introduce gates that allow selected parts of transformers to be individually eliminated. Specifically, we investigate (1) structured pruning to reduce the number of parameters in each transformer layer, (2) applicability to both BERT- and RoBERTa-based models, (3) applicability to both SQuAD 2.0 and Natural Questions, and (4) combining structured pruning with distillation. We achieve a near-doubling of inference speed with less than a 0.5 F1-point loss in short answer accuracy on Natural Questions.


title: Compensated Phase Jump at Transition Crossing in the CERN PS

abstract: The transition energy must be crossed in the CERN PS to accelerate proton and ion beams to flat-top energy. A phase jump of the accelerating RF voltage with respect to the phase of the bunches by about 180 degrees minus twice the synchronous phase is required at the instant of transition. This phase offset is injected into the beam phase loop which locks the phase of the vector sum of the RF voltage in the cavities to the beam. The polarities of the cavity return signal and of the stable phase programme are usually flipped at transition. However, both actions are difficult to perfectly synchronize in time, causing the beam phase loop to partially lock out and relock at the new stable phase. The resulting glitch can be avoided by well-controlled phase jumps applied to both, cavity drive and return signals simultaneously. This improved implementation of transition crossing makes it virtually transparent to the beam phase loop. The new scheme has been successfully tested with proton and ion beams, and it will become fully operational in the CERN PS after the long shutdown.


title: Topological states in non-Hermitian two-dimensional Su-Schrieffer-Heeger model

abstract: A non-Hermitian topological insulator with real spectrum is interesting in the theory of non-Hermitian extension of topological systems. We find an experimentally realizable example of a two dimensional non-Hermitian topological insulator with real spectrum. We consider two-dimensional Su-Schrieffer-Heeger (SSH) model with gain and loss. We introduce non-Hermitian polarization vector to explore topological phase and show that topological edge states in the band gap exist in the system.


title: On Some Aspects of Pseudonorm Compact and Montel Operators on Locally Solid Vector Lattices

abstract: Unbounded convergences have been applied successfully to locally solid topologies on vector lattices. In the present paper, we first expose several properties of various classes of Riesz pseudonorms on vector lattices. We accomplish this by abstracting some generalities of the norm of an AM-space with strong norm unit to locally solid topologies induced by a pseudonorm. By using these classes of pseudonorms, we study compactness properties of operators (not necessarily linear) between locally solid (not necessarily Hausdorff) topologies. We study new classes of operators such as pseudonorm compact, pseudo-semicompact and pseudo-AM-compact operators as well as the classical Montel operators.


title: Free Expansion of Yukawa Gas in the Constant Plasma Background

abstract: We discuss the irreversible free expansion phenomenon for Yukawa gas and obtain the corresponding change in temperature of the gas. In our expansion scheme, the system is not allowed to exchange heat with surroundings during free expansion, therefore, present expansion scheme refers to adiabatic free expansion. Using first principle classical Molecular Dynamics (MD) simulation with reflecting boundary conditions, we show that the Yukawa gases exhibits heating effect during the process of adiabatic free expansion. We also obtain the scaling for change in temperature and establish that the change in temperature is directly proportional to the change in number density of gas. The scaling for temperature is also obtained analytically taking the mean field limit of thermodynamic model proposed by Avinash. The simulation results are consistent with the analytical results.


title: A Note on a Differential Galois Approach to Path Integrals

abstract: We point out the relevance of the Differential Galois Theory of linear differential equations for the exact semiclassical computations in path integrals in quantum mechanics. The main tool will be a necessary condition for complete integrability of classical Hamiltonian systems obtained by Ramis and myself : if a finite dimensional complex analytical Hamiltonian system is completely integrable with meromorphic first integrals, then the identity component of the Galois group of the variational equation around any integral curve must be abelian. A corollary of this result is that, for finite dimensional integrable Hamiltonian systems, the semiclassical approach is computable in closed form in the framework of the Differential Galois Theory. This explains in a very precise way the success of quantum semiclassical computations for integrable Hamiltonian systems.


title: Bayesian Temporal Factorization for Multidimensional Time Series Prediction

abstract: Large-scale and multidimensional spatiotemporal data sets are becoming ubiquitous in many real-world applications such as monitoring urban traffic and air quality. Making predictions on these time series has become a critical challenge due to not only the large-scale and high-dimensional nature but also the considerable amount of missing data. In this paper, we propose a Bayesian temporal factorization (BTF) framework for modeling multidimensional time series -- in particular spatiotemporal data -- in the presence of missing values. By integrating low-rank matrix/tensor factorization and vector autoregressive (VAR) process into a single probabilistic graphical model, this framework can characterize both global and local consistencies in large-scale time series data. The graphical model allows us to effectively perform probabilistic predictions and produce uncertainty estimates without imputing those missing values. We develop efficient Gibbs sampling algorithms for model inference and model updating for real-time prediction and test the proposed BTF framework on several real-world spatiotemporal data sets for both missing data imputation and multi-step rolling prediction tasks. The numerical experiments demonstrate the superiority of the proposed BTF approaches over existing state-of-the-art methods.


title: Assessing Spatial Information in Physical Environments

abstract: Many approaches have dealt with the hypothesis that the environment contain information, mostly focusing on how humans decode information from the environment in visual perception, navigation, and spatial decision-making. A question yet to be fully explored is how the built environment could encode forms of information in its own physical structures. This paper explores a new measure of spatial information, and applies it to twenty cities from different spatial cultures and regions of the world. Findings suggest that this methodology is able to identify similarities between cities, generating a classification scheme that opens up new questions about what we call "cultural hypothesis": the idea that spatial configurations find consistent differences between cultures and regions.


title: Thresholding Bandit Problem with Both Duels and Pulls

abstract: The Thresholding Bandit Problem (TBP) aims to find the set of arms with mean rewards greater than a given threshold. We consider a new setting of TBP, where in addition to pulling arms, one can also \emph{duel} two arms and get the arm with a greater mean. In our motivating application from crowdsourcing, dueling two arms can be more cost-effective and time-efficient than direct pulls. We refer to this problem as TBP with Dueling Choices (TBP-DC). This paper provides an algorithm called Rank-Search (RS) for solving TBP-DC by alternating between ranking and binary search. We prove theoretical guarantees for RS, and also give lower bounds to show the optimality of it. Experiments show that RS outperforms previous baseline algorithms that only use pulls or duels.


title: The theory of Besov functional calculus: developments and applications to semigroups

abstract: We extend and deepen the theory of functional calculus for semigroup generators, based on the algebra $\mathcal B$ of analytic Besov functions, which we initiated in a previous paper. In particular, we show that our construction of the calculus is optimal in several natural senses. Moreover, we clarify the structure of $\mathcal B$ and identify several important subspaces in practical terms. This leads to new spectral mapping theorems for operator semigroups and to wide generalisations of a number of basic results from semigroup theory.


title: A Novel Sensorised Insole for Sensing Feet Pressure Distributions

abstract: Wearable sensors are gaining in popularity because they enable outdoor experimental monitoring. This paper presents a cost-effective sensorised insole based on a mesh of tactile capacitive sensors. Each sensor's spatial resolution is about 4 taxels/cm^2 in order to have an accurate reconstruction of the contact pressure distribution. As a consequence, the insole provides information such as contact forces, moments, and center of pressure. To retrieve this information, a calibration technique that fuses measurements from a vacuum chamber and shoes equipped with force/torque sensors is proposed. The validation analysis shows that the best performance achieved a root mean square error (RMSE) of about 7 N for the contact forces and 2 N m for the contact moments when using the force/torque shoe data as ground truth. Thus, the insole may be an alternative to force/torque sensors for certain applications, with a considerably more cost-effective and less invasive hardware.


title: Dissipative Luttinger liquids

abstract: We investigate a one dimensional quantum fluid coupled to a dissipative bath. The quantum fluid is captured by the canonical Luttinger liquid; the bath is given by the model of Caldeira and Leggett, i.e. a tower of oscillators coupled linearly to the fluid density, $\rho$. The bath can be integrated out exactly, producing an effective interaction for the fluid that is nonlocal in time; we argue that the form corresponding to Ohmic dissipation is generic. Compared to previous works, we compute correlation functions for this minimal model without approximation, including at finite temperature $T>0$. From these and a Kubo calculation, we conclude that arbitrary dissipation destroys the perfect conductivity of the Luttinger liquid via Zeno localization, even in the absence of a spatial potential; from RG analysis of harmonic terms, we also find that the open Luttinger liquid is significantly more prone to localization by such potentials, in contrast to the usual intuition that baths make systems less localized.


title: Decay rates for the damped wave equation with finite regularity damping

abstract: Decay rates for the energy of solutions of the damped wave equation on the torus are studied. In particular, damping invariant in one direction and equal to a sum of squares of nonnegative functions with a particular number of derivatives of regularity is considered. For such damping energy decays at rate $1/t^{2/3}$. If additional regularity is assumed the decay rate improves. When such a damping is smooth the energy decays at $1/t^{4/5-\delta}$. The proof uses a positive commutator argument and relies on a pseudodifferential calculus for low regularity symbols.


title: Spectral properties of the exponential distance matrix

abstract: Given a graph $G$, the exponential distance matrix is defined entry-wise by letting the $(u,v)$-entry be $q^{\text{dist}(u,v)}$, where $\text{dist}(u,v)$ is the distance between the vertices $u$ and $v$ with the convention that if vertices are in different components, then $q^{\text{dist}(u,v)}=0$. In this paper, we will establish several properties of the characteristic polynomial (spectrum) for this matrix, give some families of graphs which are uniquely determined by their spectrum, and produce cospectral constructions.


title: Searching for the largest bound atoms in space

abstract: (abridged) Radio recombination lines (RRLs) at frequencies $\nu$ < 250 MHz trace the cold, diffuse phase of the ISM. Next generation low frequency interferometers, such as LOFAR, MWA and the future SKA, with unprecedented sensitivity, resolution, and large fractional bandwidths, are enabling the exploration of the extragalactic RRL universe. We observed the radio quasar 3C 190 (z~1.2) with the LOFAR HBA. In reducing this data for spectroscopic analysis, we have placed special emphasis on bandpass calibration. We devised cross-correlation techniques to significantly identify the presence of RRLs in a low frequency spectrum. We demonstrate the utility of this method by applying it to existing low-frequency spectra of Cassiopeia A and M 82, and to the new observations of 3C 190. RRLs have been detected in the foreground of 3C 190 at z = 1.12355 (assuming a carbon origin), owing to the first detection of RRLs outside of the local universe (first reported in Emig et al. 2019). Towards the Galactic supernova remnant Cas A, we uncover three new detections: (1) C$\epsilon$-transitions ($\Delta$n = 5) for the first time at low radio frequencies, (2) H$\alpha$-transitions at 64 MHz with a FWHM of 3.1 km/s, the most narrow and one of the lowest frequency detections of hydrogen to date, and (3) C$\alpha$ at v$_{LSR}$ = 0 km/s in the frequency range 55-78 MHz for the first time. Additionally we recover C$\alpha$, C$\beta$, C$\gamma$, and C$\delta$ from the -47 km/s and -38 km/s components. In the nearby starburst galaxy, M 82, we do not find a significant feature. Our current searches for RRLs in LOFAR observations are limited to narrow (< 100 km/s) features, owing to the relatively small number of channels available for continuum estimation. Future strategies making use of larger contiguous frequency coverage would aid calibration to deeper sensitivities and broader features.


title: The Sounds of Music : Science of Musical Scales III -- Indian Classical

abstract: In the previous articles of this series, we have discussed the development of musical scales particularly that of the heptatonic scale which forms the basis of Western classical music today. In this last article, we take a look at the basic structure of scales used in Indian classical music and how different `raga's are generated through the simple process of scale shifting.


title: Ultrafast enhancement of ferromagnetic spin exchange induced by ligand-to-metal charge transfer

abstract: Directly modifying spin exchange energies in a magnetic material with light can enable ultrafast control of its magnetic states. Current approaches rely on tuning charge hopping amplitudes that mediate exchange by optically exciting either virtual or real charge-transfer transitions (CT) between magnetic sites. Here we show that when exchange is mediated by a non-magnetic ligand, it can be substantially enhanced by optically exciting a real CT transition from the ligand to magnetic site, introducing lower order virtual hopping contributions. We demonstrate sub-picosecond enhancement in a superexchange dominated ferromagnet CrSiTe3 through this mechanism using phase-resolved coherent phonon spectroscopy. This technique can also be applied in the paramagnetic phase to disentangle light induced exchange modification from other ultrafast effects that alter the magnetization. This protocol can potentially be broadly applied to engineer thermally inaccessible spin Hamiltonians in superexchange dominated magnets.


title: Cayley Graphs on Billiard Surfaces, and Their Genus

abstract: In this article we discuss a connection between two famous constructions in mathematics: a Cayley graph of a group and a (rational) billiard surface. For each rational billiard surface, there is a natural way to draw a Cayley graph of a dihedral group on that surface. Both of these objects have the concept of "genus" attached to them. For the Cayley graph, the genus is defined to be the lowest genus amongst all surfaces that the graph can be drawn on without edge crossings. We prove that the genus of the Cayley graph associated to a billiard surface arising from a triangular billiard table is always zero or one. One reason this is interesting is that there exist triangular billiard surfaces of arbitrarily high genus , so the genus of the associated graph is usually much lower than the genus of the billiard surface.


title: SCAFFOLD: Stochastic Controlled Averaging for Federated Learning

abstract: Federated Averaging (FedAvg) has emerged as the algorithm of choice for federated learning due to its simplicity and low communication cost. However, in spite of recent research efforts, its performance is not fully understood. We obtain tight convergence rates for FedAvg and prove that it suffers from `client-drift' when the data is heterogeneous (non-iid), resulting in unstable and slow convergence. As a solution, we propose a new algorithm (SCAFFOLD) which uses control variates (variance reduction) to correct for the `client-drift' in its local updates. We prove that SCAFFOLD requires significantly fewer communication rounds and is not affected by data heterogeneity or client sampling. Further, we show that (for quadratics) SCAFFOLD can take advantage of similarity in the client's data yielding even faster convergence. The latter is the first result to quantify the usefulness of local-steps in distributed optimization.


title: Dual-path RNN: efficient long sequence modeling for time-domain single-channel speech separation

abstract: Recent studies in deep learning-based speech separation have proven the superiority of time-domain approaches to conventional time-frequency-based methods. Unlike the time-frequency domain approaches, the time-domain separation systems often receive input sequences consisting of a huge number of time steps, which introduces challenges for modeling extremely long sequences. Conventional recurrent neural networks (RNNs) are not effective for modeling such long sequences due to optimization difficulties, while one-dimensional convolutional neural networks (1-D CNNs) cannot perform utterance-level sequence modeling when its receptive field is smaller than the sequence length. In this paper, we propose dual-path recurrent neural network (DPRNN), a simple yet effective method for organizing RNN layers in a deep structure to model extremely long sequences. DPRNN splits the long sequential input into smaller chunks and applies intra- and inter-chunk operations iteratively, where the input length can be made proportional to the square root of the original sequence length in each operation. Experiments show that by replacing 1-D CNN with DPRNN and apply sample-level modeling in the time-domain audio separation network (TasNet), a new state-of-the-art performance on WSJ0-2mix is achieved with a 20 times smaller model than the previous best system.


title: Bridging Information Security and Environmental Criminology Research to Better Mitigate Cybercrime

abstract: Cybercrime is a complex phenomenon that spans both technical and human aspects. As such, two disjoint areas have been studying the problem from separate angles: the information security community and the environmental criminology one. Despite the large body of work produced by these communities in the past years, the two research efforts have largely remained disjoint, with researchers on one side not benefitting from the advancements proposed by the other. In this paper, we argue that it would be beneficial for the information security community to look at the theories and systematic frameworks developed in environmental criminology to develop better mitigations against cybercrime. To this end, we provide an overview of the research from environmental criminology and how it has been applied to cybercrime. We then survey some of the research proposed in the information security domain, drawing explicit parallels between the proposed mitigations and environmental criminology theories, and presenting some examples of new mitigations against cybercrime. Finally, we discuss the concept of cyberplaces and propose a framework in order to define them. We discuss this as a potential research direction, taking into account both fields of research, in the hope of broadening interdisciplinary efforts in cybercrime research.


title: Principled estimation of regression discontinuity designs

abstract: Regression discontinuity designs are frequently used to estimate the causal effect of election outcomes and policy interventions. In these contexts, treatment effects are typically estimated with covariates included to improve efficiency. While including covariates improves precision asymptotically, in practice, treatment effects are estimated with a small number of observations, resulting in considerable fluctuations in treatment effect magnitude and precision depending upon the covariates chosen. This practice thus incentivizes researchers to select covariates which maximize treatment effect statistical significance rather than precision. Here, I propose a principled approach for estimating RDDs which provides a means of improving precision with covariates while minimizing adverse incentives. This is accomplished by integrating the adaptive LASSO, a machine learning method, into RDD estimation using an R package developed for this purpose, adaptiveRDD. Using simulations, I show that this method significantly improves treatment effect precision, particularly when estimating treatment effects with fewer than 200 observations.


title: Particle-number distribution in large fluctuations at the tip of branching random walks

abstract: We investigate properties of the particle distribution near the tip of one-dimensional branching random walks at large times $t$, focusing on unusual realizations in which the rightmost lead particle is very far ahead of its expected position - but still within a distance smaller than the diffusion radius $\sim\sqrt{t}$. Our approach consists in a study of the generating function $G_{\Delta x}(\lambda)=\sum_n \lambda^n p_n(\Delta x)$ for the probabilities $p_n(\Delta x)$ of observing $n$ particles in an interval of given size $\Delta x$ from the lead particle to its left, fixing the position of the latter. This generating function can be expressed with the help of functions solving the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation with suitable initial conditions. In the infinite-time and large-$\Delta x$ limits, we find that the mean number of particles in the interval grows exponentially with $\Delta x$, and that the generating function obeys a nontrivial scaling law, depending on $\Delta x$ and $\lambda$ through the combined variable $[\Delta x-f(\lambda)]^{3}/\Delta x^2$, where $f(\lambda)\equiv -\ln(1-\lambda)-\ln[-\ln(1-\lambda)]$. From this property, one may conjecture that the growth of the typical particle number with the size of the interval is slower than exponential, but, surprisingly enough, only by a subleading factor at large $\Delta x$. The scaling we argue is consistent with results from a numerical integration of the FKPP equation.


title: A TESS Search for Distant Solar System Planets: A Feasibility Study

abstract: As the NASA Transiting Exoplanet Survey Satellite (TESS) fulfills its primary mission it is executing an unprecedented all-sky survey with the potential to discover distant planets in our own solar system, as well as hundreds of Transneptunian Objects (TNOs) and Centaurs. We demonstrate that shift-and-stack techniques can be used to efficiently search the Full-Frame Image (FFI) data from the TESS mission and survey the entire sky for outer Solar System objects down to $\sim22^{nd}$ magnitude.


title: Cost Sharing over Combinatorial Domains: Complement-Free Cost Functions and Beyond

abstract: We study mechanism design for combinatorial cost sharing. Imagine that multiple items or services are available to be shared among a set of interested agents. The outcome of a mechanism in this setting consists of an assignment, determining for each item the set of players who are granted service, together with respective payments. Although there are several works studying specialized versions of such problems, there has been almost no progress for general combinatorial cost sharing domains until recently \cite{DobzinskiO17}. The main goal of our work is to further understand this interplay in terms of budget balance and social cost approximation. Towards this, we provide a refinement of cross-monotonicity (trace-monotonicity) that is applicable to iterative mechanisms. The trace here refers to the order in which players become finalized. On top of this, we also provide two parameterizations of cost functions which capture the behavior of their average cost-shares. Based on our trace-monotonicity property, we design a scheme of ascending cost sharing mechanisms which is applicable to the combinatorial cost sharing setting with symmetric submodular valuations. Using our first cost function parameterization, we identify conditions under which our mechanism is weakly group-strategyproof, $O(1)$-budget-balanced and $O(H_n)$-approximate with respect to the social cost. Finally, we consider general valuation functions and exploit our second parameterization to derive a more fine-grained analysis of the Sequential Mechanism introduced by Moulin. This mechanism is budget balanced by construction, but in general only guarantees a poor social cost approximation of $n$. We identify conditions under which the mechanism achieves improved social cost approximation guarantees.


title: Decomposition of complete tripartite graphs into 5-cycles using the graph adjacency matrix

abstract: The problem of decomposing a complete tripartite graph into 5-cycles was first proposed in 1995 by Mahmoodian and Mirzakhani and since then many attempts have been made to decompose such graphs into 5-cycles. Such attempts were partially successful but parts of the problem still remain open. In this paper, we investigate the problem deeper, decompose more tripartite graphs into 5-cycles, and introduce the Graph Adjacency Matrix (GAM) method for cycle decomposition in general. GAM method transforms the cycle decomposition problem to covering squares with certain polygons. This new formulation is easier to solve and enables us to find explicit decompositions for numerous cases that were not decomposed before.


title: All of Linear Regression

abstract: Least squares linear regression is one of the oldest and widely used data analysis tools. Although the theoretical analysis of the ordinary least squares (OLS) estimator is as old, several fundamental questions are yet to be answered. Suppose regression observations $(X_1,Y_1),\ldots,(X_n,Y_n)\in\mathbb{R}^d\times\mathbb{R}$ (not necessarily independent) are available. Some of the questions we deal with are as follows: under what conditions, does the OLS estimator converge and what is the limit? What happens if the dimension is allowed to grow with $n$? What happens if the observations are dependent with dependence possibly strengthening with $n$? How to do statistical inference under these kinds of misspecification? What happens to the OLS estimator under variable selection? How to do inference under misspecification and variable selection? We answer all the questions raised above with one simple deterministic inequality which holds for any set of observations and any sample size. This implies that all our results are a finite sample (non-asymptotic) in nature. In the end, one only needs to bound certain random quantities under specific settings of interest to get concrete rates and we derive these bounds for the case of independent observations. In particular, the problem of inference after variable selection is studied, for the first time, when $d$, the number of covariates increases (almost exponentially) with sample size $n$. We provide comments on the ``right'' statistic to consider for inference under variable selection and efficient computation of quantiles.


title: Interference dislocations in condensate of indirect excitons

abstract: Phase singularities in quantum states play a significant role both in the state properties and in the transition between the states. For instance, a transition to two-dimensional superfluid state is governed by pairing of vortices and, in turn, unpaired vortices can cause dissipations for particle fluxes. Vortices and other phase defects can be revealed by characteristic features in interference patterns produced by the quantum system. We present dislocation-like phase singularities in interference patterns in a condensate of indirect excitons measured by shift-interferometry. We show that the observed dislocations in interference patterns are not associated with conventional phase defects: neither with vortices, nor with polarization vortices, nor with half-vortices, nor with skyrmions, nor with half-skyrmions. We present the origin of these new phase singularities in condensate interference patterns: the observed interference dislocations originate from converging of the condensate matter waves propagating from different sources.


title: Comparing the density of D_4 and S_4 quartic extensions of number fields

abstract: When ordered by discriminant, it is known that about 83% of quartic fields over Q have associated Galois group S_4, while the remaining 17% have Galois group D_4. We study these proportions over a general number field F. We find that asymptotically 100% of quadratic number fields have more D_4 extensions than S_4 and that the ratio between the number of D_4 and S_4 quartic extensions is biased arbitrarily in favor of D_4 extensions. Under GRH, we give a lower bound that holds for general number fields.


title: Snow Crystals

abstract: This monograph reviews our current understanding of the physical dynamics of ice crystal growth, focusing on the spontaneous formation of complex structures from water vapor (called snow crystals) as a function of temperature, supersaturation, background gas pressure, and other extrinsic parameters. Snow crystal growth is a remarkably rich and rather poorly understood phenomenon, requiring a synthesis of concepts from materials science, crystal-growth theory, statistical mechanics, diffusion-limited solidification, finite-element modeling, and molecular surface processes. Building upon recent advances in precision measurement techniques, computation modeling methods, and molecular dynamics simulations of crystalline surfaces, I believe we are moving rapidly toward the long-sought goal of developing a full physical model of snow crystal formation, using ab initio molecular dynamics simulations to create a semi-empirical characterization of the nanoscale surface attachment kinetics, and then incorporating that into a full computational model that reproduces the growth of macroscopic crystalline structures. Section 1 of this monograph deals mainly with the material properties of ice Ih in equilibrium, including thermodynamics quantities, facet surface structures, terrace step energies, and crystal twinning behaviors.


title: Optimal Paired Comparison Block Designs

abstract: For paired comparison experiments involving competing options described by two-level attributes several different methods of constructing designs having block paired observations under the main effects model are presented. These designs are compared to alternative designs available in the literature.


title: Building Information Modeling and Classification by Visual Learning At A City Scale

abstract: In this paper, we provide two case studies to demonstrate how artificial intelligence can empower civil engineering. In the first case, a machine learning-assisted framework, BRAILS, is proposed for city-scale building information modeling. Building information modeling (BIM) is an efficient way of describing buildings, which is essential to architecture, engineering, and construction. Our proposed framework employs deep learning technique to extract visual information of buildings from satellite/street view images. Further, a novel machine learning (ML)-based statistical tool, SURF, is proposed to discover the spatial patterns in building metadata. The second case focuses on the task of soft-story building classification. Soft-story buildings are a type of buildings prone to collapse during a moderate or severe earthquake. Hence, identifying and retrofitting such buildings is vital in the current earthquake preparedness efforts. For this task, we propose an automated deep learning-based procedure for identifying soft-story buildings from street view images at a regional scale. We also create a large-scale building image database and a semi-automated image labeling approach that effectively annotates new database entries. Through extensive computational experiments, we demonstrate the effectiveness of the proposed method.


title: AFP-CKSAAP: Prediction of Antifreeze Proteins Using Composition of k-Spaced Amino Acid Pairs with Deep Neural Network

abstract: Antifreeze proteins (AFPs) are the sub-set of ice binding proteins indispensable for the species living in extreme cold weather. These proteins bind to the ice crystals, hindering their growth into large ice lattice that could cause physical damage. There are variety of AFPs found in numerous organisms and due to the heterogeneous sequence characteristics, AFPs are found to demonstrate a high degree of diversity, which makes their prediction a challenging task. Herein, we propose a machine learning framework to deal with this vigorous and diverse prediction problem using the manifolding learning through composition of k-spaced amino acid pairs. We propose to use the deep neural network with skipped connection and ReLU non-linearity to learn the non-linear mapping of protein sequence descriptor and class label. The proposed antifreeze protein prediction method called AFP-CKSAAP has shown to outperform the contemporary methods, achieving excellent prediction scores on standard dataset. The main evaluater for the performance of the proposed method in this study is Youden's index whose high value is dependent on both sensitivity and specificity. In particular, AFP-CKSAAP yields a Youden's index value of 0.82 on the independent dataset, which is better than previous methods.


title: In-training Matrix Factorization for Parameter-frugal Neural Machine Translation

abstract: In this paper, we propose the use of in-training matrix factorization to reduce the model size for neural machine translation. Using in-training matrix factorization, parameter matrices may be decomposed into the products of smaller matrices, which can compress large machine translation architectures by vastly reducing the number of learnable parameters. We apply in-training matrix factorization to different layers of standard neural architectures and show that in-training factorization is capable of reducing nearly 50% of learnable parameters without any associated loss in BLEU score. Further, we find that in-training matrix factorization is especially powerful on embedding layers, providing a simple and effective method to curtail the number of parameters with minimal impact on model performance, and, at times, an increase in performance.


title: Strong submeasures and applications to non-compact dynamical systems

abstract: A strong submeasure on a compact metric space X is a sub-linear and bounded operator on the space of continuous functions on X. A strong submeasure is positive if it is non-decreasing. By Hahn-Banach theorem, a positive strong submeasure is the supremum of a non-empty collection of measures whose masses are uniformly bounded from above. There are many natural examples of continuous maps of the forms $f:U\rightarrow X$, where $X$ is a compact metric space and $U\subset X$ is an open dense subsets, where $f$ cannot extend to a reasonable function on $X$. We can mention cases such as: transcendental maps of $\mathbb{C}$, meromorphic maps on compact complex varieties, or continuous self-maps $f:U\rightarrow U$ of a dense open subset $U\subset X$ where $X$ is a compact metric space. To the maps mentioned in the previous paragraph, the use of measures is not sufficient to establish the basic properties of ergodic theory, such as the existence of invariant measures or a reasonable definition of measure theoretic entropy and topological entropy. In this paper, we show that strong submeasures can be used to completely resolve the issue and establish these basic properties. The same idea can also be applied to improve on \'etale topological entropy, previously defined by the author. In another paper we apply strong submeasures to intersection of positive closed $(1,1)$ currents on compact K\"ahler manifolds.


title: Interaction of SN~Ib 2004dk with a Previously-Expelled Envelope

abstract: The interaction between the expanding supernova (SN) ejecta with the circumstellar material (CSM) that was expelled from the progenitor prior to explosion is a long-sought phenomenon, yet observational evidence is scarce. Here we confirm a new example: SN 2004dk, originally a hydrogen-poor, helium-rich Type Ib SN that reappeared as a strong H$\alpha$-emitting point-source on narrowband H$\alpha$ images. We present follow-up optical spectroscopy that reveals the presence of a broad H$\alpha$ component with full width at half maximum of ~290 km/s in addition to the narrow H$\alpha$ +[NII] emission features from the host galaxy. Such a broad component is a clear sign of an ejecta-CSM interaction. We also present observations with the XMM-Newton Observatory, the Swift satellite, and the Chandra X-ray Observatory that span 10 days to 15 years after discovery. The detection of strong radio, X-ray, and H$\alpha$ emission years after explosion allows various constraints to be put on pre-SN mass-loss processes. We present a wind-bubble model in which the CSM is "pre-prepared" by a fast wind interacting with a slow wind. Much of the outer density profile into which the SN explodes corresponds to no steady-state mass-loss process. We estimate that the shell of compressed slow wind material was ejected ~1400 yr prior to explosion, perhaps during carbon burning, and that the SN shock had swept up about 0.04 M_sun of material. The region emitting the H$\alpha$ has a density of order $10^{-20}$ g/cc.


title: Nebula-Relay Hypothesis: Primitive Life in Nebula and Origin of Life on Earth

abstract: A modified version of panspermia theory, named Nebula-Relay hypothesis or local panspermia, is introduced to explain the origin of life on Earth. Primitive life, acting as the seeds of life on Earth, originated at pre-solar epoch through physicochemical processes and then filled in the pre-solar nebula after the death of pre-solar star. Then the history of life on the Earth can be divided into three epochs: the formation of primitive life in the pre-solar epoch; pre-solar nebula epoch; the formation of solar system and the Earth age of life. The main prediction of our model is that primitive life existed in the pre-solar nebula (even in the current nebulas) and the celestial body formed therein (i.e. solar system).


title: Demisability and survivability sensitivity to design-for-demise techniques

abstract: The paper is concerned with examining the effects that design-for-demise solutions can have not only on the demisability of components, but also on their survivability that is their capability to withstand impacts from space debris. First two models are introduced. A demisability model to predict the behaviour of spacecraft components during the atmospheric re-entry and a survivability model to assess the vulnerability of spacecraft structures against space debris impacts. Two indices that evaluate the level of demisability and survivability are also proposed. The two models are then used to study the sensitivity of the demisability and of the survivability indices as a function of typical design-for-demise options. The demisability and the survivability can in fact be influenced by the same design parameters in a competing fashion that is while the demisability is improved, the survivability is worsened and vice versa. The analysis shows how the design-for-demise solutions influence the demisability and the survivability independently. In addition, the effect that a solution has simultaneously on the two criteria is assessed. Results shows which, among the design-for-demise parameters mostly influence the demisability and the survivability. For such design parameters maps are presented, describing their influence on the demisability and survivability indices. These maps represent a useful tool to quickly assess the level of demisability and survivability that can be expected from a component, when specific design parameters are changed.


title: Plasmonic sensing with FIB-milled 2D micro-arrays of truncated gold nano-pyramids

abstract: Plasmonic platforms are a promising solution for the next generation of low-cost, integrated biomedical sensors. However, fabricating these nanostructures often requires lengthy and challenging fabrication processes. Here we exploit the peculiar features of Focused Ion Beam milling to obtain single-step patterning of a plasmonic two-dimensional (2D) array of truncated gold nano-pyramids (TNP), with gaps as small as 17 nm. We describe the formation of plasmonic bandgaps in the arrangement of crossed tapered grooves that separate the nano-pyramids, and we demonstrate refractive index sensing via dark field imaging in patterned areas of 30 $\mu$m $\times$ 30 $\mu$m.


title: Weakly maximal subgroups of branch groups

abstract: Let $G$ be a branch group acting by automorphisms on a rooted tree $T$. Stabilizers of infinite rays in $T$ are examples of weakly maximal subgroups of $G$ (subgroups that are maximal among subgroups of infinite index), but in general they are not the only examples. In this note we describe two families of weakly maximal subgroups of branch groups. We show that, for the first Grigorchuk group as well as for the torsion GGS groups, every weakly maximal subgroup belongs to one of these families. The first family is a generalization of stabilizers of rays, while the second one consists of weakly maximal subgroups with a block structure. We obtain different equivalent characterizations of these families in terms of finite generation, the existence of a trivial rigid stabilizer, the number of orbit-closures for the action on the boundary of the tree or by the means of sections.


title: The importance of magnification effects in galaxy-galaxy lensing

abstract: Magnification changes the observed number counts of galaxies on the sky. This biases the observed tangential shear profiles around galaxies, the so-called galaxy-galaxy lensing (GGL) signal, and the related excess mass profile. Correspondingly, inference of physical quantities, such as the mean mass profile of halos around galaxies, are affected by magnification effects. We use simulated shear and galaxy data of the Millennium Simulation to quantify the effect on shear and mass estimates from magnified lens and source number counts. The former are due to the large-scale matter distribution in the foreground of the lenses, the latter are caused by magnification of the source population by the matter associated with the lenses. The GGL signal is calculated from the simulations by an efficient fast-Fourier transform that can also be applied to real data. The numerical treatment is complemented by a leading-order analytical description of the magnification effects, which is shown to fit the numerical shear data well. We find the magnification effect is strongest for steep galaxy luminosity functions and high redshifts. For a lens redshift of $z_\mathrm{d}=0.83$, a limiting magnitude of $22\,\mathrm{mag}$ in the $r$-band and a source redshift of $z_\mathrm{s}=0.99$, we find that a magnification correction changes the shear profile up to $45\%$ and the mass is biased by up to $55 \%$. For medium-redshift galaxies the relative change in shear and mass is typically a few percent. As expected, the sign of the bias depends on the local slope of the lens luminosity function $\alpha_\mathrm{d}$, where the mass is biased low for $\alpha_\mathrm{d}<1$ and biased high for $\alpha_\mathrm{d}>1$. Whereas the magnification effect of sources is rarely than more $1\%$, the statistical power of future weak lensing surveys warrants correction for this effect.


title: Physics-Informed Deep Neural Network Method for Limited Observability State Estimation

abstract: The precise knowledge regarding the state of the power grid is important in order to ensure optimal and reliable grid operation. Specifically, knowing the state of the distribution grid becomes increasingly important as more renewable energy sources are connected directly into the distribution network, increasing the fluctuations of the injected power. In this paper, we consider the case when the distribution grid becomes partially observable, and the state estimation problem is under-determined. We present a new methodology that leverages a deep neural network (DNN) to estimate the grid state. The standard DNN training method is modified to explicitly incorporate the physical information of the grid topology and line/shunt admittance. We show that our method leads to a superior accuracy of the estimation when compared to the case when no physical information is provided. Finally, we compare the performance of our method to the standard state estimation approach, which is based on the weighted least squares with pseudo-measurements, and show that our method performs significantly better with respect to the estimation accuracy.


title: Addressing Troubles with Double Bubbles: Convergence and Stability at Multi-Bubble Junctions

abstract: In this report we discuss and propose a correction to a convergence and stability issue occurring in the work of Da et al.[2015], in which they proposed a numerical model to simulate soap bubbles.


title: BoTorch: A Framework for Efficient Monte-Carlo Bayesian Optimization

abstract: Bayesian optimization provides sample-efficient global optimization for a broad range of applications, including automatic machine learning, engineering, physics, and experimental design. We introduce BoTorch, a modern programming framework for Bayesian optimization that combines Monte-Carlo (MC) acquisition functions, a novel sample average approximation optimization approach, auto-differentiation, and variance reduction techniques. BoTorch's modular design facilitates flexible specification and optimization of probabilistic models written in PyTorch, simplifying implementation of new acquisition functions. Our approach is backed by novel theoretical convergence results and made practical by a distinctive algorithmic foundation that leverages fast predictive distributions, hardware acceleration, and deterministic optimization. We also propose a novel "one-shot" formulation of the Knowledge Gradient, enabled by a combination of our theoretical and software contributions. In experiments, we demonstrate the improved sample efficiency of BoTorch relative to other popular libraries.


title: Learning to Generate 6-DoF Grasp Poses with Reachability Awareness

abstract: Motivated by the stringent requirements of unstructured real-world where a plethora of unknown objects reside in arbitrary locations of the surface, we propose a voxel-based deep 3D Convolutional Neural Network (3D CNN) that generates feasible 6-DoF grasp poses in unrestricted workspace with reachability awareness. Unlike the majority of works that predict if a proposed grasp pose within the restricted workspace will be successful solely based on grasp pose stability, our approach further learns a reachability predictor that evaluates if the grasp pose is reachable or not from robot's own experience. To avoid the laborious real training data collection, we exploit the power of simulation to train our networks on a large-scale synthetic dataset. This work is an early attempt that simultaneously evaluates grasping reachability from learned knowledge while proposing feasible grasp poses with 3D CNN. Experimental results in both simulation and real-world demonstrate that our approach outperforms several other methods and achieves 82.5% grasping success rate on unknown objects.


title: Turn Skipping and the Game Coloring Number

abstract: The game coloring number gcol($G$) of a graph $G$ is a two player competitive variant of the coloring number. We introduce the preordered game coloring number to study the consequences of either player skipping any number of turns. In particular, we show that neither player can improve their performance by doing so. We use this result to show that for any induced subgraph $H \subset G$, $|G-H| = k$ implies the tight bound gcol($H$) $\leq$ gcol($G$) $-$ $2k$ and if $d_G(x)$ $\leq$ gcol($G-x$), then gcol($G-x$) $\leq$ gcol($G$) $-$ $1$.


title: Clouds in Higher Dimensions

abstract: Following some work done by Komjath and Schmerl, we extend the definition of a cloud to $\mathbb{R}^N$ for $N \geq 2$ and show that $k$ clouds cover $\mathbb{R}^2$ if and only if $k$ clouds cover $\mathbb{R}^N$. We also show that countably many clouds cover $\mathbb{R}^N$.


title: FireNet: Real-time Segmentation of Fire Perimeter from Aerial Video

abstract: In this paper, we share our approach to real-time segmentation of fire perimeter from aerial full-motion infrared video. We start by describing the problem from a humanitarian aid and disaster response perspective. Specifically, we explain the importance of the problem, how it is currently resolved, and how our machine learning approach improves it. To test our models we annotate a large-scale dataset of 400,000 frames with guidance from domain experts. Finally, we share our approach currently deployed in production with inference speed of 20 frames per second and an accuracy of 92 (F1 Score).


title: Continuous renormalization group $\beta$ function from lattice simulations

abstract: We present a real-space renormalization group transformation with continuous scale change to calculate the continuous renormalization group $\beta$ function in non-perturbative lattice simulations. Our method is motivated by the connection between Wilsonian renormalization group and the gradient flow transformation. It does not rely on the perturbative definition of the renormalized coupling and is also valid at non-perturbative fixed points. Although our method requires an additional extrapolation compared to traditional step scaling calculations, it has several advantages which compensates for this extra step even when applied in the vicinity of the perturbative fixed point. We illustrate our approach by calculating the $\beta$ function of 2-flavor QCD and show that lattice predictions from individual lattice ensembles, even without the required continuum and finite volume extrapolations, can be very close to the result of the full analysis. Thus our method provides a non-perturbative framework and intuitive understanding into the structure of strongly coupled systems, in addition to being complementary to existing lattice determinations.


title: An Efficient Numerical Method for the Generalised Kolmogorov Equation

abstract: An efficient algorithm for computing the terms appearing in the Generalised Kolmogorov Equation (GKE) written for the indefinite plane channel flow is presented. The algorithm, which features three distinct strategies for parallel computing, is designed such that CPU and memory requirements are kept to a minimum, so that high-Re wall-bounded flows can be afforded. Computational efficiency is mainly achieved by leveraging the Parseval's theorem for the two homogeneous directions available in the plane channel geometry. A speedup of 3-4 orders of magnitude, depending on the problem size, is reported in comparison to a key implementation used in the literature. Validation of the code is demonstrated by computing the residual of the GKE, and example results are presented for channel flows at Re_{\tau}=200 and Re_{\tau}=1000, where for the first time they are observed in the whole four-dimensional domain. It is shown that the space and scale properties of the scale-energy fluxes change for increasing values of the Reynolds number. Among all scale-energy fluxes, the wall-normal flux is found to show the richest behaviour for increasing streamwise scales.


title: Excitation spectrum and supersolidity of a two-leg bosonic ring ladder

abstract: We consider a system of weakly interacting bosons confined on a planar double lattice ring subjected to two artificial gauge fields. This system is known to display three phases, the Meissner phase where the flow of particles is carried at the edges of the system without transverse current, a vortex phase characterized by non-zero transverse current, and a biased-ladder phase, characterized by an imbalance of the population of the two rings. We use the Bogoliubov approximation to determine the excitation spectrum in the three phases, the dynamic structure factor and the quantum fluctuation corrections to the first-order correlation function. Our analysis reveals supersolid features as well as Josephson modes, corresponding to out-of-phase modes of the finite ring.


title: Mapping Supervised Bilingual Word Embeddings from English to low-resource languages

abstract: It is very challenging to work with low-resource languages due to the inadequate availability of data. Using a dictionary to map independently trained word embeddings into a shared vector space has proved to be very useful in learning bilingual embeddings in the past. Here we have tried to map individual embeddings of words in English and their corresponding translated words in low-resource languages like Estonian, Slovenian, Slovakian, and Hungarian. We have used a supervised learning approach. We report accuracy scores through various retrieval strategies which show that it is possible to approach challenging tasks in Natural Language Processing like machine translation for such languages, provided that we have at least some amount of proper bilingual data. We also conclude that we can follow an unsupervised learning path on monolingual text data as that is more suitable for low-resource languages.


title: The impact of catastrophic collisions and collision avoidance on a swarming behavior

abstract: Swarms of autonomous agents are useful in many applications due to their ability to accomplish tasks in a decentralized manner, making them more robust to failures. Due to the difficulty in running experiments with large numbers of hardware agents, researchers often make simplifying assumptions and remove constraints that might be present in a real swarm deployment. While simplifying away some constraints is tolerable, we feel that two in particular have been overlooked: one, that agents in a swarm take up physical space, and two, that agents might be damaged in collisions. Many existing works assume agents have negligible size or pass through each other with no added penalty. It seems possible to ignore these constraints using collision avoidance, but we show using an illustrative example that this is easier said than done. In particular, we show that collision avoidance can interfere with the intended swarming behavior and significant parameter tuning is necessary to ensure the behavior emerges as best as possible while collisions are avoided. We compare four different collision avoidance algorithms, two of which we consider to be the best decentralized collision avoidance algorithms available. Despite putting significant effort into tuning each algorithm to perform at its best, we believe our results show that further research is necessary to develop swarming behaviors that can achieve their goal while avoiding collisions with agents of non-negligible volume.


title: Nodeless superconducting gap in the candidate topological superconductor Sn$_{1-x}$In$_x$Te for x = 0.7

abstract: High-pressure synthesis techniques have allowed for the growth of Sn$_{1-x}$In$_x$Te samples beyond the ambient In-saturation limit of $x$ = 0.5 (T$_c \sim$ 4.5 K). In this study, we present measurements of the temperature dependence of the London penetration depth $\Delta\lambda(T)$ in this superconducting doped topological insulator for $x$ = 0.7, where T$_{c,onset}\approx 5$ K. The results indicate fully gapped BCS-like behavior, ruling out odd-parity $A_{2u}$ pairing; however, odd-parity $A_{1u}$ pairing is still possible. Critical field values measured below 1 K and other superconducting parameters are also presented.


title: Exporting Terrestrial Life Out of the Solar System with Gravitational Slingshots of Earthgrazing Bodies

abstract: Exporting terrestrial life out of the Solar System requires a process that both embeds microbes in boulders and ejects those boulders out of the Solar System. We explore the possibility that Earthgrazing long-period comets and interstellar objects could export life from Earth by collecting microbes from the atmosphere and receiving a gravitational slingshot effect from the Earth. We estimate the total number of exportation events over the lifetime of the Earth to be $\sim 1 - 10$ for long-period comets and $\sim 1 - 50$ for interstellar objects. If life existed above an altitude of 100 km, then the number is dramatically increased up to $\sim 10^5$ exportation events over Earth's lifetime.


abstract: High-Performance Computing (HPC) platforms enable scientific software to achieve breakthroughs in many research fields such as physics, biology, and chemistry, by employing Research Software Engineering (RSE) techniques. These include 1) novel parallelism paradigms such as Shared Memory Parallelism (with e.g. OpenMP 4.5); Distributed Memory Parallelism (with e.g. MPI 4); Hybrid Parallelism which combines them; and Heterogeneous Parallelism (for CPUs, co-processors and accelerators), 2) introducing advanced Software Engineering concepts such as Object Oriented Parallel Programming (OOPP); Parallel Unit testing; Parallel I/O Formats; Hybrid Parallel Visualization; and 3) Selecting the Best Practices in other necessary areas such as User Interface; Automatic Documentation; Version Control and Project Management. In this work we present BACKUS: Comprehensive High-Performance Research Software Engineering Approach for Simulations in Supercomputing Systems, which we found to fit best for long-lived parallel scientific codes.


title: RecSplit: Minimal Perfect Hashing via Recursive Splitting

abstract: A minimal perfect hash function bijectively maps a key set $S$ out of a universe $U$ into the first $|S|$ natural numbers. Minimal perfect hash functions are used, for example, to map irregularly-shaped keys, such as string, in a compact space so that metadata can then be simply stored in an array. While it is known that just $1.44$ bits per key are necessary to store a minimal perfect function, no published technique can go below $2$ bits per key in practice. We propose a new technique for storing minimal perfect hash functions with expected linear construction time and expected constant lookup time that makes it possible to build for the first time, for example, structures which need $1.56$ bits per key, that is, within $8.3$% of the lower bound, in less than $2$ ms per key. We show that instances of our construction are able to simultaneously beat the construction time, space usage and lookup time of the state-of-the-art data structure reaching $2$ bits per key. Moreover, we provide parameter choices giving structures which are competitive with alternative, larger-size data structures in terms of space and lookup time. The construction of our data structures can be easily parallelized or mapped on distributed computational units (e.g., within the MapReduce framework), and structures larger than the available RAM can be directly built in mass storage.


title: GEANT4 Parameter Tuning Using Professor

abstract: The Geant4 toolkit is used extensively in high energy physics to simulate the passage of particles through matter and to predict effects such as detector efficiencies and smearing. Geant4 uses many underlying models to predict particle interaction kinematics, and uncertainty in these models leads to uncertainty in high energy physics measurements. The Geant4 collaboration recently made free parameters in some models accessible through partnership with Geant4 developers. We present a study of the impact of varying parameters in three Geant4 hadronic physics models on agreement with thin target datasets and describe fits to these datasets using the Professor model tuning framework. We find that varying parameters produces substantially better agreement with some datasets, but that more degrees of freedom are required for full agreement. This work is a first step towards a common framework for propagating uncertainties in Geant4 models to high energy physics measurements, and we outline future work required to complete that goal.


title: Constructing curvelet-like bases and low-redundancy frames

abstract: We provide a detailed analysis of the obstruction (studied first by S. Durand and later by R. Yin and one of us) in the construction of multidirectional wavelet orthonormal bases corresponding to any admissible frequency partition in the framework of subband filtering with non-uniform subsampling. To contextualize our analysis, we build, in particular, multidirectional alias-free hexagonal wavelet bases and low-redundancy frames with optimal spatial decay. In addition, we show that a 2D cutting lemma can be used to subdivide the obtained wavelet systems in higher frequency rings so as to generate bases or frames that satisfy the ``parabolic scaling'' law enjoyed by curvelets and shearlets. Numerical experiments on high bit-rate image compression are conducted to illustrate the potential of the proposed systems.