title: Coupling a mobile hole to an antiferromagnetic spin background: Transient dynamics of a magnetic polaron

abstract: Understanding the interplay between charge and spin and its effects on transport is a ubiquitous challenge in quantum many-body systems. In the Fermi-Hubbard model, this interplay is thought to give rise to magnetic polarons, whose dynamics may explain emergent properties of quantum materials such as high-temperature superconductivity. In this work, we use a cold-atom quantum simulator to directly observe the formation dynamics and subsequent spreading of individual magnetic polarons. Measuring the density- and spin-resolved evolution of a single hole in a 2D Hubbard insulator with short-range antiferromagnetic correlations reveals fast initial delocalization and a dressing of the spin background, indicating polaron formation. At long times, we find that dynamics are slowed down by the spin exchange time, and they are compatible with a polaronic model with strong density and spin coupling. Our work enables the study of out-of-equilibrium emergent phenomena in the Fermi-Hubbard model, one dopant at a time.


title: A DSO Framework for Comprehensive Market Participation of DER Aggregators

abstract: In this paper, a distribution system operator (DSO) framework is proposed to optimally coordinate distributed energy resources (DER) aggregators' comprehensive participation in retail energy market as well as wholesale energy and regulation markets. Various types of DER aggregators, including energy storage aggregators (ESAGs), dispatchable distributed generation aggregators (DDGAGs), electric vehicles charging stations (EVCSs), and demand response aggregators (DRAGs), are modeled in the proposed DSO framework. Distribution network constraints are considered by using a linearized power flow. The problem is modeled using mixed-integer linear programming (MILP) which can be solved by commercial solvers. Case studies are performed to analyze the interactions between DER aggregators and wholesale/retail electricity markets.


title: Corona Games: Masks, Social Distancing and Mechanism Design

abstract: Pandemic response is a complex affair. Most governments employ a set of quasi-standard measures to fight COVID-19 including wearing masks, social distancing, virus testing and contact tracing. We argue that some non-trivial factors behind the varying effectiveness of these measures are selfish decision-making and the differing national implementations of the response mechanism. In this paper, through simple games, we show the effect of individual incentives on the decisions made with respect to wearing masks and social distancing, and how these may result in a sub-optimal outcome. We also demonstrate the responsibility of national authorities in designing these games properly regarding the chosen policies and their influence on the preferred outcome. We promote a mechanism design approach: it is in the best interest of every government to carefully balance social good and response costs when implementing their respective pandemic response mechanism.


title: Towards Early Diagnosis of Epilepsy from EEG Data

abstract: Epilepsy is one of the most common neurological disorders, affecting about 1% of the population at all ages. Detecting the development of epilepsy, i.e., epileptogenesis (EPG), before any seizures occur could allow for early interventions and potentially more effective treatments. Here, we investigate if modern machine learning (ML) techniques can detect EPG from intra-cranial electroencephalography (EEG) recordings prior to the occurrence of any seizures. For this we use a rodent model of epilepsy where EPG is triggered by electrical stimulation of the brain. We propose a ML framework for EPG identification, which combines a deep convolutional neural network (CNN) with a prediction aggregation method to obtain the final classification decision. Specifically, the neural network is trained to distinguish five second segments of EEG recordings taken from either the pre-stimulation period or the post-stimulation period. Due to the gradual development of epilepsy, there is enormous overlap of the EEG patterns before and after the stimulation. Hence, a prediction aggregation process is introduced, which pools predictions over a longer period. By aggregating predictions over one hour, our approach achieves an area under the curve (AUC) of 0.99 on the EPG detection task. This demonstrates the feasibility of EPG prediction from EEG recordings.


title: Training Generative Adversarial Networks with Limited Data

abstract: Training generative adversarial networks (GAN) using too little data typically leads to discriminator overfitting, causing training to diverge. We propose an adaptive discriminator augmentation mechanism that significantly stabilizes training in limited data regimes. The approach does not require changes to loss functions or network architectures, and is applicable both when training from scratch and when fine-tuning an existing GAN on another dataset. We demonstrate, on several datasets, that good results are now possible using only a few thousand training images, often matching StyleGAN2 results with an order of magnitude fewer images. We expect this to open up new application domains for GANs. We also find that the widely used CIFAR-10 is, in fact, a limited data benchmark, and improve the record FID from 5.59 to 2.42.


title: Frontiers in Mortar Methods for Isogeometric Analysis

abstract: Complex geometries as common in industrial applications consist of multiple patches, if spline based parametrizations are used. The requirements for the generation of analysis-suitable models are increasing dramatically since isogeometric analysis is directly based on the spline parametrization and nowadays used for the calculation of higher-order partial differential equations. The computational, or more general, the engineering analysis necessitates suitable coupling techniques between the different patches. Mortar methods have been successfully applied for coupling of patches and for contact mechanics in recent years to resolve the arising issues within the interface. We present here current achievements in the design of mortar technologies in isogeometric analysis within the Priority Program SPP 1748, Reliable Simulation Techniques in Solid Mechanics. Development of Non-standard Discretisation Methods, Mechanical and Mathematical Analysis.


title: Constraints from gravitational wave detections of binary black hole mergers on the $^{12}\rm{C}\left(\alpha,\gamma\right)^{16}\!\rm{O}$ rate

abstract: Gravitational wave detections are starting to allow us to probe the physical processes in the evolution of very massive stars through the imprints they leave on their final remnants. Stellar evolution theory predicts the existence of a gap in the black hole mass distribution at high mass due to the effects of pair-instability. Previously, we showed that the location of the gap is robust against model uncertainties, but it does depend sensitively on the uncertain $^{12}\rm{C}\left(\alpha,\gamma\right)^{16}\!\rm{O}$ rate. This rate is of great astrophysical significance and governs the production of oxygen at the expense of carbon. We use the open source MESA stellar evolution code to evolve massive helium stars to probe the location of the mass gap. We find that the maximum black hole mass below the gap varies between $40\rm{M}_\odot$ to $90\rm{M}_\odot$, depending on the strength of the uncertain $^{12}\rm{C}\left(\alpha,\gamma\right)^{16}\!\rm{O}$ reaction rate. With the first ten gravitational-wave detections of black holes, we constrain the astrophysical S-factor for $^{12}\rm{C}\left(\alpha,\gamma\right)^{16}\!\rm{O}$, at $300\rm{keV}$, to $S_{300}>175\rm{\,keV\, barns}$ at 68% confidence. With $\mathcal{O}(50)$ detected binary black hole mergers, we expect to constrain the S-factor to within $\pm10$-$30\rm{\,keV\, barns}$. We also highlight a role for independent constraints from electromagnetic transient surveys. The unambiguous detection of pulsational pair instability supernovae would imply that $S_{300}>79\rm{\,keV\, barns}$. Degeneracies with other model uncertainties need to be investigated further, but probing nuclear stellar astrophysics poses a promising science case for the future gravitational wave detectors.


title: Wide binaries are rare in open clusters

abstract: The population statistics of binary stars are an important output of star formation models. However populations of wide binaries evolve over time due to interactions within a system's birth environment and the unfolding of wide, hierarchical triple systems. Hence the wide binary populations observed in star forming regions or OB associations may not accurately reflect the wide binary populations that will eventually reach the field. We use Gaia DR2 data to select members of three open clusters, Alpha~Per, the Pleiades and Praesepe and to flag cluster members that are likely unresolved binaries due to overluminosity or elevated astrometric noise. We then identify the resolved wide binary population in each cluster, separating it from coincident pairings of unrelated cluster members. We find that these clusters have an average wide binary fraction in the 300-3000\,AU projected separation range of 2.1$\pm^{0.4}_{0.2}$\% increasing to 3.0$\pm^{0.8}_{0.7}$\% for primaries with masses in the 0.5-1.5\,$M_{\odot}$ range. This is significantly below the observed field wide binary fraction, but shows some wide binaries survive in these dynamically highly processed environments. We compare our results with another open cluster (the Hyades) and two populations of young stars that likely originated in looser associations (Young Moving Groups and the Pisces-Eridanus stream). We find that the Hyades also has a deficit of wide binaries while the products of looser associations have wide binary fractions at or above field level.


title: Assessing the Impact of Space School UK

abstract: Space School UK (SSUK) is a series of summer residential programmes for secondary school aged students, held at the University of Leicester over 3 weeks each year. Each programme involves space-related activities run by a team of mentors - currently including university students, graduates, teachers and young professionals associated with the space sector - all of whom attended SSUK as students themselves. It includes the 6-day Space School UK and the 8-day Senior Space School UK (collectively SSUK) which are for 13-15 and 16-18 year olds respectively. This paper seeks to evaluate and present the benefits of SSUK to individuals who participate in the programme, organisations involved in the running of SSUK, and to highlight and promote these benefits to the wider UK and global space community. We also address which facets of SSUK make for such an engaging and encouraging experience for the students, that are missing from students' traditional education. We seek to show how SSUK acts as an excellent example of how to bridge the gap between secondary and tertiary space education. Through an analysis of our alumni survey results, we show that attending SSUK has a significant impact on career choices and prepares students for Higher Education, regardless of background. Some groups, such as women, and those from non-selective schools, reported a higher impact in some of these areas than others. Metrics such as skills learned, goals achieved, alongside knowledge of careers and Higher Education possibilities are discussed for various demographics.


title: Simulating the "hidden giant" in cold and self-interacting dark matter models

abstract: We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial orbit has a small pericenter. After evolving in the Milky Way's tidal field, the simulated galaxies experience significant mass loss and their stellar distributions expand accordingly. These tidal effects are more prominent if the initial halo concentration is lower and if the self-scattering cross section is larger. Our results show that Antlia 2-like galaxies could be realized in CDM if the halo concentration is low and the stellar distribution is diffuse at the infall time, while these conditions could be relaxed in SIDM. We also find all the simulated galaxies predict approximately the same stellar velocity dispersion after imposing selection criteria for stellar particles. This has important implications for testing dark matter models using tidally disturbed systems.


title: Power System Disturbance Classification with Online Event-Driven Neuromorphic Computing

abstract: Accurate online classification of disturbance events in a transmission network is an important part of wide-area monitoring. Although many conventional machine learning techniques are very successful in classifying events, they rely on extracting information from PMU data at control centers and processing them through CPU/GPUs, which are highly inefficient in terms of energy consumption. To solve this challenge without compromising accuracy, this paper presents a novel methodology based on event-driven neuromorphic computing architecture for classification of power system disturbances. A Spiking Neural Network (SNN)-based computing framework is proposed, which exploits sparsity in disturbances and promotes local event driven operation for unsupervised learning and inference from incoming data. Spatio-temporal information of PMU signals is first extracted and encoded into spike trains and classification is achieved with SNN-based supervised and unsupervised learning framework. Moreover, a QR decomposition-based selection technique is proposed to identify signals participating in the low rank subspace of multiple disturbance events. Performance of the proposed method is validated on data collected from a 16-machine, 5-area New England-New York system.


title: Evolution and observational signatures of the cosmic ray electron spectrum in SN 1006

abstract: Supernova remnants (SNRs) are believed to be the source of Galactic cosmic rays (CRs). SNR shocks accelerate CR protons and electrons which reveal key insights into the non-thermal physics by means of their synchrotron and $\gamma$-ray emission. The remnant SN 1006 is an ideal particle acceleration laboratory because it is observed across all electromagnetic wavelengths from radio to $\gamma$-rays. We perform three-dimensional (3D) magnetohydrodynamics (MHD) simulations where we include CR protons and follow the CR electron spectrum. By matching the observed morphology and non-thermal spectrum of SN 1006 in radio, X-rays and $\gamma$-rays, we gain new insight into CR electron acceleration and magnetic field amplification. 1. We show that a mixed leptonic-hadronic model is responsible for the $\gamma$-ray radiation: while leptonic inverse-Compton emission and hadronic pion-decay emission contribute equally at GeV energies observed by Fermi, TeV energies observed by imaging air Cherenkov telescopes are hadronically dominated. 2. We show that quasi-parallel acceleration (i.e., when the shock propagates at a narrow angle to the upstream magnetic field) is preferred for CR electrons and that the electron acceleration efficiency of radio-emitting GeV electrons at quasi-perpendicular shocks is suppressed at least by a factor ten. This precludes extrapolation of current one-dimensional plasma particle-in-cell simulations of shock acceleration to realistic SNR conditions. 3. To match the radial emission profiles and the $\gamma$-ray spectrum, we require a volume-filling, turbulently amplified magnetic field and that the Bell-amplified magnetic field is damped in the immediate post-shock region. Our work connects micro-scale plasma physics simulations to the scale of SNRs.


title: Fainter harder brighter softer: a correlation between $\alpha_{\rm ox}$ , X-ray spectral state and Eddington ratio in tidal disruption events

abstract: We explore the accretion states of tidal disruption events (TDEs) using a sample of 7 X- ray bright sources. To this end, we estimate the relative contribution of the disk and corona to the observed X-ray emission through spectral modeling, and assess the X-ray brightness (through ${\alpha}$_ox, L$_{2 keV}$ and f$_{Edd,X}$) as a function of the Eddington ratio. We report strong positive correlations between ${\alpha}$_ox and f$_{Edd,bol}$; f$_{Edd,X}$ and f$_{Edd,UV}$ ; and an anti-correlation for L$_{2 keV}$ and f$_{Edd,UV}$. TDEs at high f$_{Edd,bol}$ have thermal dominated X-ray spectra and high (soft) ${\alpha}$_ox, whereas those at low f$_{Edd,bol}$ show a significant power-law contribution and low (hard) ${\alpha}$_ox. Similar to X-ray binaries and active galactic nuclei, the transition between X-ray spectral states occurs around f$_{Edd,bol}$ ${\approx}$ 0.03, although the uncertainty is large due to the small sample size. Our results suggest that X-ray surveys are more likely to discover TDEs at low f$_{Edd,bol}$, whereas optical surveys are more sensitive to TDEs at high Eddington ratios. The X-ray and optical selected TDEs have different UV and X-ray properties, which should be taken into account when deriving rates, luminosity and black hole mass functions. TDEs around the most massive supermassive black holes are observed in the hard state; this could indicate that TDE evolution is faster around more massive BHs.


title: Invertible Networks or Partons to Detector and Back Again

abstract: For simulations where the forward and the inverse directions have a physics meaning, invertible neural networks are especially useful. A conditional INN can invert a detector simulation in terms of high-level observables, specifically for ZW production at the LHC. It allows for a per-event statistical interpretation. Next, we allow for a variable number of QCD jets. We unfold detector effects and QCD radiation to a pre-defined hard process, again with a per-event probabilistic interpretation over parton-level phase space.


title: Resolving the Hubble Tension with New Early Dark Energy

abstract: New Early Dark Energy (NEDE) is a component of vacuum energy at the electron volt scale, which decays in a first-order phase transition shortly before recombination [arXiv:1910.10739]. The NEDE component has the potential to resolve the tension between recent local measurements of the expansion rate of the Universe using supernovae (SN) data and the expansion rate inferred from the early Universe through measurements of the cosmic microwave background (CMB) when assuming $\Lambda$CDM. We discuss in depth the two-scalar field model of the NEDE phase transition including the process of bubble percolation, collision, and coalescence. We also estimate the gravitational wave signal produced during the collision phase and argue that it can be searched for using pulsar timing arrays. In a second step, we construct an effective cosmological model, which describes the phase transition as an instantaneous process, and derive the covariant equations that match perturbations across the transition surface. Fitting the cosmological model to CMB, baryonic acoustic oscillations and SN data, we report $H_0 = 69.6^{+1.0}_{-1.3} \, \textrm{km}\, \textrm{s}^{-1}\, \textrm{Mpc}^{-1}$ $(68 \%$ C.L.) without the local measurement of the Hubble parameter, bringing the tension down to $2.5\, \sigma$. Including the local input, we find $H_0 = 71.4 \pm 1.0 \, \textrm{km}\, \textrm{s}^{-1}\, \textrm{Mpc}^{-1}$ $(68 \%$ C.L.) and strong evidence for a non-vanishing NEDE component with a $\simeq 4\, \sigma$ significance.


title: On the asymptotics of wide networks with polynomial activations

abstract: We consider an existing conjecture addressing the asymptotic behavior of neural networks in the large width limit. The results that follow from this conjecture include tight bounds on the behavior of wide networks during stochastic gradient descent, and a derivation of their finite-width dynamics. We prove the conjecture for deep networks with polynomial activation functions, greatly extending the validity of these results. Finally, we point out a difference in the asymptotic behavior of networks with analytic (and non-linear) activation functions and those with piecewise-linear activations such as ReLU.


title: The imprint of arms and bars on rotation curves: in-plane and off-plane

abstract: Within Rotation Curves (RC) is encoded the kinematical state of the stellar disc as well as information about the dynamical mechanisms driving the secular evolution of galaxies. To explain the characteristic features of RCs that arise by the influence of spiral patterns and bar, we study the kinematics of the stellar disc in a set of spiral galaxy models specifically tailored for this purpose. We find that, for our models, the induced non-circular motions are more prominent for spirals with larger pitch angle, the ones typical in late type galaxies. Moreover, inside corotation, stars rotate slower along the spiral arms than along the inter-arm, which translates into a local minima or maxima in the RC, respectively. We also see, from off-plane RC, that the rotation is faster for stars that at observed closer to the plane, and diminishes as one looks farther off plane; this trend is more noticeable in our Sa galaxy model than our Sc galaxy model. Additionally, in a previous work we found that the diagonal ridges in the $V_{\phi}-R$ plane, revealed through the Gaia DR2, have a resonant origin due to the spiral arms and bar and that these ridges project themselves as wiggles in the RC; here, we further notice that the development of these ridges, and the development of high orbital eccentricities in the stellar disc are the same. Hence, we conclude that, the following explanations of bumps and wiggles in RCs are equivalent: they are manifestations of diagonal ridges in the $V_{\phi}-R$ plane, or of the rearrangement of the orbital eccentricities in the stellar disc.


title: Radio Emission from Ultra-stripped Supernovae as Diagnostics for Properties of the Remnant Double Neutron Star Binaries

abstract: An ultra-stripped supernova (SN) is an explosion of a helium or C+O star whose outer envelope has been stripped away by a companion neutron star. A double neutron star (DNS) binary is believed to be left after the explosion, which will emit the gravitational wave later at the coalescence. Recent detections of a few candidates for the ultra-stripped SN have constrained the properties of the explosion and the progenitor, but little information is given as to whether the remnant DNS binary will merge within the cosmic age. A large fraction of the material stripped away from the helium star through the binary interaction is expected to escape from the system and form circumstellar material (CSM). The CSM should be traced by radio emission induced by the collision with the SN ejecta. Based on the stellar evolution models previously developed, we calculate the expected radio luminosities from ultra-stripped SNe. We find that high radio luminosity at its maximum can be an indicator of small separation of a DNS binary leading to its merger within the cosmic age. Our results can be used to optimize the strategy for the radio follow-up observations such as observational epochs and frequencies.


title: Discovery and Follow-up of ASASSN-19dj: An X-ray and UV Luminous TDE in an Extreme Post-Starburst Galaxy

abstract: We present observations of ASASSN-19dj, a nearby tidal disruption event (TDE) discovered in the post-starburst galaxy KUG 0810+227 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of d $\simeq98$ Mpc. We observed ASASSN-19dj from $-$21 to 392 d relative to peak ultraviolet (UV)/optical emission using high-cadence, multiwavelength spectroscopy and photometry. From the ASAS-SN $g$-band data, we determine that the TDE began to brighten on 2019 February 6.8 and for the first 16 d the rise was consistent with a flux $\propto t^2$ power-law. ASASSN-19dj peaked in the UV/optical on 2019 March 6.5 (MJD = 58548.5) at a bolometric luminosity of $L = (6.2 \pm 0.2) \times 10^{44} \text{ erg s}^{-1}$. Initially remaining roughly constant in X-rays and slowly fading in the UV/optical, the X-ray flux increased by over an order of magnitude $\sim$225 d after peak, resulting from the expansion of the X-ray emitting region. The late-time X-ray emission is well fitted by a blackbody with an effective radius of $\sim1 \times 10^{12} \text{ cm}$ and a temperature of $\sim6 \times 10^{5} \text{ K}$. The X-ray hardness ratio becomes softer after brightening and then returns to a harder state as the X-rays fade. Analysis of Catalina Real-Time Transient Survey images reveals a nuclear outburst roughly 14.5 yr earlier with a smooth decline and a luminosity of $L_V\geq1.4 \times 10^{43}$ erg s$^{-1}$, although the nature of the flare is unknown. ASASSN-19dj occurred in the most extreme post-starburst galaxy yet to host a TDE, with Lick H$\delta_{A}$ = $7.67\pm0.17$ \AA.


title: SPISEA: A Python-Based Simple Stellar Population Synthesis Code for Star Clusters

abstract: We present SPISEA (Stellar Population Interface for Stellar Evolution and Atmospheres), an open-source Python package that simulates simple stellar populations. The strength of SPISEA is its modular interface which offers the user control of 13 input properties including (but not limited to) the Initial Mass Function, stellar multiplicity, extinction law, and the metallicity-dependent stellar evolution and atmosphere model grids used. The user also has control over the Initial-Final Mass Relation in order to produce compact stellar remnants (black holes, neutron stars, and white dwarfs). We demonstrate several outputs produced by the code, including color-magnitude diagrams, HR-diagrams, luminosity functions, and mass functions. SPISEA is object-oriented and extensible, and we welcome contributions from the community. The code and documentation are available on GitHub and ReadtheDocs, respectively.


title: Topological Correlators and Surface Defects from Equivariant Cohomology

abstract: We find a one-dimensional protected subsector of $\mathcal{N}=4$ matter theories on a general class of three-dimensional manifolds. By means of equivariant localization, we identify a dual quantum mechanics computing BPS correlators of the original model in three dimensions. Specifically, applying the Atiyah-Bott-Berline-Vergne formula to the original action demonstrates that this localizes on a one-dimensional action with support on the fixed-point submanifold of suitable isometries. We first show that our approach reproduces previous results obtained on $S^3$. Then, we apply it to the novel case of $S^2 \times S^1$ and show that the theory localizes on two noninteracting quantum mechanics with disjoint support. We prove that the BPS operators of such models are naturally associated with a noncommutative star product, while their correlation functions are essentially topological. Finally, we couple the three-dimensional theory to general $\mathcal{N}=(2,2)$ surface defects and extend the localization computation to capture the full partition function and BPS correlators of the mixed-dimensional system.


title: Characterizing Quantum Networks: Insights from Coherence Theory

abstract: Networks based on entangled quantum systems enable interesting applications in quantum information processing and the understanding of the resulting quantum correlations is essential for advancing the technology. We show that the theory of quantum coherence provides powerful tools for analyzing this problem. For that, we demonstrate that a recently proposed approach to network correlations based on covariance matrices can be improved and analytically evaluated for the most important cases.


title: UV Shadows in EFTs: Accidental Symmetries, Robustness and No-Scale Supergravity

abstract: We argue that accidental approximate scaling symmetries are robust predictions of weakly coupled string vacua, and show that their interplay with supersymmetry and other (generalised) internal symmetries underlies the ubiquitous appearance of no-scale supergravities in low-energy 4D EFTs. We identify 4 nested types of no-scale supergravities, and show how leading quantum corrections can break scale invariance while preserving some no-scale properties (including non-supersymmetric flat directions). We use these ideas to classify corrections to the low-energy 4D supergravity action in perturbative 10D string vacua, including both bulk and brane contributions. Our prediction for the K\"ahler potential at any fixed order in $\alpha'$ and string loops agrees with all extant calculations. p-form fields play two important roles: they spawn many (generalised) shift symmetries; and space-filling 4-forms teach 4D physics about higher-dimensional phenomena like flux quantisation. We argue that these robust symmetry arguments suffice to understand obstructions to finding classical de Sitter vacua, and suggest how to get around them in UV complete models.


title: And In The Darkness Unbind Them: High-Resolution Simulations of Dark Matter Subhalo Disruption in a Milky Way-like Tidal Field

abstract: We compare the results of high-resolution simulations of individual dark matter subhalos evolving in external tidal fields with and without baryonic bulge and disk components, where the average dark matter particle mass is three orders of magnitude smaller than cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation is used to setup our initial conditions and provides a basis for our simulations of subhalos in a dark matter-only tidal field, while an observationally motivated model for the Milky Way is used for the tidal field that is comprised of a dark matter halo, a stellar disk, and a stellar bulge. Our simulations indicate that including stellar components in the tidal field results in the number of subhalos in Milky Way-like galaxies being only $65\%$ of what is predicted by $\Lambda$ Cold Dark Matter ($\Lambda$CDM). For subhalos with small pericentres $(r_p \lesssim 25$ kpc), the subhalo abundance is reduced further to $40\%$, with the surviving subhalos being less dense and having a tangentially-anisotropic orbital distribution. Conversely, subhalos with larger pericentres are minimally affected by the inclusion of a stellar component in the tidal field, with the total number of outer subhalos $\approx 75\%$ of the $\Lambda$CDM prediction. The densities of outer subhalos are comparable to predictions from $\Lambda$CDM, with the subhalos having an isotropic distribution of orbits. These ratios are higher than those found in previous studies that include the effects baryonic matter, which are affected by spurious disruption caused by low resolution.


title: Duality and Supersymmetry Constraints on the Weak Gravity Conjecture

abstract: Positivity bounds coming from consistency of UV scattering amplitudes are in general insufficient to prove the weak gravity conjecture for theories beyond Einstein-Maxwell. Additional ingredients about the UV may be necessary to exclude those regions of parameter space which are na\"ively in conflict with the predictions of the weak gravity conjecture. In this paper we explore the consequences of imposing additional symmetries inherited from the UV theory on higher-derivative operators for Einstein-Maxwell-dilaton-axion theory. Using black hole thermodynamics, for a preserved SL($2,\mathbb{R}$) symmetry we find that the weak gravity conjecture then does follow from positivity bounds. For a preserved O($d,d;\mathbb{R}$) symmetry we find a simple condition on the two Wilson coefficients which ensures the positivity of corrections to the charge-to-mass ratio and that follows from the null energy condition alone. We find that imposing supersymmetry on top of either of these symmetries gives corrections which vanish identically, as expected for BPS states.


title: Diffuse Ionized Gas in Simulations of Multiphase, Star-Forming Galactic Disks

abstract: It has been hypothesized that photons from young, massive star clusters are responsible for maintaining the ionization of diffuse warm ionized gas seen in both the Milky Way and other disk galaxies. For a theoretical investigation of the warm ionized medium (WIM), it is crucial to solve radiation transfer equations where the ISM and clusters are modeled self-consistently. To this end, we employ a Solar neighborhood model of TIGRESS, a magnetohydrodynamic simulation of the multiphase, star-forming ISM, and post-process the simulation with an adaptive ray tracing method to transfer UV radiation from star clusters. We find that the WIM volume filling factor is highly variable, and sensitive to the rate of ionizing photon production and ISM structure. The mean WIM volume filling factor rises to ~0.15 at |z|~1 kpc. Approximately half of ionizing photons are absorbed by gas and half by dust; the cumulative ionizing photon escape fraction is 1.1%. Our time-averaged synthetic H$\alpha$ line profile matches WHAM observations on the redshifted (outflowing) side, but has insufficient intensity on the blueshifted side. Our simulation matches the Dickey-Lockman neutral density profile well, but only a small fraction of snapshots have high-altitude WIM density consistent with Reynolds Layer estimates. We compute a clumping correction factor C = <n_e>/sqrt<n_e^2>~0.2 that is remarkably constant with distance from the midplane and time; this can be used to improve estimates of ionized gas mass and mean electron density from observed H$\alpha$ surface brightness profiles in edge-on galaxies.


title: Non-minimal Tinges of Unimodular Gravity

abstract: Unimodular Gravity is normally assumed to be equivalent to General Relativity for all matters but the character of the Cosmological Constant. Here we discuss this equivalence in the presence of a non-minimally coupled scalar field. We show that when we consider gravitation to be dynamical in a QFT sense, quantum corrections can distinguish both theories if the non-minimal coupling is non-vanishing. In order to show this, we construct a path integral formulation of Unimodular Gravity, fixing the complicated gauge invariance of the theory and computing all one-loop divergences. We find a combination of the couplings in the Lagrangian to which we can assign a physical meaning. It tells whether quantum gravitational phenomena can be ignored or not at a given energy scale. Its renormalization group flow differs depending on if it is computed in General Relativity or Unimodular Gravity.


title: Mechanical oscillator thermometry in the nonlinear optomechanical regime

abstract: Optomechanical systems are promising platforms for controlled light-matter interactions. They are capable of providing several fundamental and practical novel features when the mechanical oscillator is cooled down to nearly reach its ground state. In this framework, measuring the effective temperature of the oscillator is perhaps the most relevant step in the characterization of those systems. In conventional schemes, the cavity is driven strongly, and the overall system is well-described by a linear (Gaussian preserving) Hamiltonian. Here, we depart from this regime by considering an undriven optomechanical system via non-Gaussian radiation-pressure interaction. To measure the temperature of the mechanical oscillator, initially in a thermal state, we use light as a probe to coherently interact with it and create an entangled state. We show that the optical probe gets a nonlinear phase, resulting from the non-Gaussian interaction, and undergoes an incoherent phase diffusion process. To efficiently infer the temperature from the entangled light-matter state, we propose using a nonlinear Kerr medium before a homodyne detector. Remarkably, placing the Kerr medium enhances the precision to nearly saturate the ultimate quantum bound given by the quantum Fisher information. Furthermore, it also simplifies the thermometry procedure as it makes the choice of the homodyne local phase independent of the temperature, which avoids the need for adaptive sensing protocols.


title: Observational constraints on the slope of the radial acceleration relation at low accelerations

abstract: The radial acceleration relation (RAR) locally relates the `observed' acceleration inferred from the dynamics of a system to the acceleration implied by its baryonic matter distribution. The relation as traced by galaxy rotation curves is one-to-one with remarkably little scatter, implying that the dynamics of a system can be predicted simply by measuring its density profile as traced by e.g. stellar light or gas emission lines. Extending the relation to accelerations below those usually probed by practically observable kinematic tracers is challenging, especially once accounting for faintly emitting baryons, such as the putative warm-hot intergalactic medium, becomes important. We show that in the low-acceleration regime, the (inverted) RAR predicts an unphysical, declining enclosed baryonic mass profile for systems with `observed' acceleration profiles steeper than $g_{\rm obs}\propto r^{-1}$ (corresponding to density profiles steeper than isothermal - $\rho(r)\propto r^{-2}$). If the RAR is tantamount to a natural law, such acceleration profiles cannot exist. We apply this argument to test the compatibility of an extrapolation of the rotation curve-derived RAR to low accelerations with data from galaxy-galaxy weak lensing, dwarf spheroidal galaxy stellar kinematic, and outer Milky~Way dynamical measurements, fully independent of the uncertainties inherent in direct measurements of the baryonic matter distribution. In all cases we find that the data weakly favour a break to a steeper low-acceleration slope. Improvements in measurements and modelling of the outer Milky~Way, and weak lensing, seem like the most promising path toward stronger constraints on the low-acceleration behaviour of the RAR.


title: Probing the Planck scale: The modification of the time evolution operator due to the quantum structure of spacetime

abstract: The propagator which evolves the wave-function in NRQM, can be expressed as a matrix element of a time evolution operator: i.e $ G_{\rm NR}(x)= \langle{\mathbf{x}_2}|{U_{\rm NR}(t)}|{\mathbf{x}_1}\rangle$ in terms of the orthonormal eigenkets $|{\mathbf{x}}\rangle$ of the position operator. In QFT, it is not possible to define a conceptually useful single-particle position operator or its eigenkets. It is also not possible to interpret the relativistic (Feynman) propagator $G_R(x)$ as evolving any kind of single-particle wave-functions. In spite of all these, it is indeed possible to express the propagator of a free spinless particle, in QFT, as a matrix element $\langle{\mathbf{x}_2}|{U_{\rm R}(t)}|{\mathbf{x}_1}\rangle$ for a suitably defined time evolution operator and (non-orthonormal) kets $|{\mathbf{x}}\rangle$ labeled by spatial coordinates. At mesoscopic scales, which are close but not too close to Planck scale, one can incorporate quantum gravitational corrections to the propagator by introducing a zero-point-length. It turns out that even this QG corrected propagator can be expressed as a matrix element $\langle{\mathbf{x}_2}|{U_{\rm QG}(t)}|{\mathbf{x}_1}\rangle$. I describe these results and explore several consequences. It turns out that the evolution operator $U_{\rm QG}(t)$ becomes non-unitary for sub-Planckian time intervals while remaining unitary for time interval is larger than Planck time. The results can be generalised to any ultrastatic curved spacetime.


title: Demographics of neutron stars in young massive and open clusters

abstract: Star clusters appear to be the ideal environment for the assembly of neutron star-neutron star (NS-NS) and black hole-neutron star (BH-NS) binaries. These binaries are among the most interesting astrophysical objects, being potential sources of gravitational waves (GWs) and gamma-ray bursts. We use for the first time high-precision N-body simulations of young massive and open clusters to study the origin and dynamical evolution of NSs, within clusters with different initial masses, metallicities, primordial binary fractions, and prescriptions for the compact object natal kicks at birth. We find that the radial profile of NSs is shaped by the BH content of the cluster, which partially quenches the NS segregation due to the BH-burning process. This leaves most of the NSs out of the densest cluster regions, where NS-NS and BH-NS binaries could potentially form. Due to a large velocity kick that they receive at birth, most of the NSs escape the host clusters, with the bulk of their retained population made up of NSs of $\sim 1.3$ M$_\odot$ coming from the electron-capture supernova process. The details of the primordial binary fraction and pairing can smear out this trend. Finally, we find that a subset of our models produce NS-NS mergers, leading to a rate of $\sim 0.01$--$0.1$ Gpc$^{-3}$ yr$^{-1}$ in the local Universe, and compute an upper limit of $\sim 3\times 10^{-2}$--$3\times 10^{-3}$ Gpc$^{-3}$ yr$^{-1}$ for the BH-NS merger rate. Our estimates are several orders of magnitude smaller than the current empirical merger rate from LIGO/Virgo, in agreement with the recent rate estimates for old globular clusters.


title: The Massive and Distant Clusters of WISE Survey X: Initial Results from a Sunyaev-Zeldovich Effect Study of Massive Galaxy Clusters at z>1 using MUSTANG2 on the GBT

abstract: The properties of galaxy clusters as a function of redshift can be utilized as an important cosmological tool. We present initial results from a program of follow-up observations of the Sunyaev-Zeldovich effect (SZE) in high redshift galaxy clusters detected at infrared wavelengths in the Massive and Distant Clusters of WISE Survey (MaDCoWS). Using typical on-source integration times of 3-4 hours per cluster, MUSTANG2 on the Green Bank Telescope was able to measure strong detections of SZE decrements and statistically significant masses on 14 out of 16 targets. On the remaining two, weaker (3.7 sigma) detections of the SZE signal and strong upper limits on the masses were obtained. In this paper we present masses and pressure profiles of each target and outline the data analysis used to recover these quantities. Of the clusters with strong detections, three show significantly flatter pressure profiles while, from the MUSTANG2 data, five others show signs of disruption at their cores. However, outside of the cores of the clusters, we were unable to detect significant amounts of asymmetry. Finally, there are indications that the relationship between optical richness used by MaDCoWS and SZE-inferred mass may be significantly flatter than indicated in previous studies.


title: End-to-end Sinkhorn Autoencoder with Noise Generator

abstract: In this work, we propose a novel end-to-end sinkhorn autoencoder with noise generator for efficient data collection simulation. Simulating processes that aim at collecting experimental data is crucial for multiple real-life applications, including nuclear medicine, astronomy and high energy physics. Contemporary methods, such as Monte Carlo algorithms, provide high-fidelity results at a price of high computational cost. Multiple attempts are taken to reduce this burden, e.g. using generative approaches based on Generative Adversarial Networks or Variational Autoencoders. Although such methods are much faster, they are often unstable in training and do not allow sampling from an entire data distribution. To address these shortcomings, we introduce a novel method dubbed end-to-end Sinkhorn Autoencoder, that leverages sinkhorn algorithm to explicitly align distribution of encoded real data examples and generated noise. More precisely, we extend autoencoder architecture by adding a deterministic neural network trained to map noise from a known distribution onto autoencoder latent space representing data distribution. We optimise the entire model jointly. Our method outperforms competing approaches on a challenging dataset of simulation data from Zero Degree Calorimeters of ALICE experiment in LHC. as well as standard benchmarks, such as MNIST and CelebA.


title: Optimizing generalization on the train set: a novel gradient-based framework to train parameters and hyperparameters simultaneously

abstract: Generalization is a central problem in Machine Learning. Most prediction methods require careful calibration of hyperparameters carried out on a hold-out \textit{validation} dataset to achieve generalization. The main goal of this paper is to present a novel approach based on a new measure of risk that allows us to develop novel fully automatic procedures for generalization. We illustrate the pertinence of this new framework in the regression problem. The main advantages of this new approach are: (i) it can simultaneously train the model and perform regularization in a single run of a gradient-based optimizer on all available data without any previous hyperparameter tuning; (ii) this framework can tackle several additional objectives simultaneously (correlation, sparsity,...) $via$ the introduction of regularization parameters. Noticeably, our approach transforms hyperparameter tuning as well as feature selection (a combinatorial discrete optimization problem) into a continuous optimization problem that is solvable via classical gradient-based methods ; (iii) the computational complexity of our methods is $O(npK)$ where $n,p,K$ denote respectively the number of observations, features and iterations of the gradient descent algorithm. We observe in our experiments a significantly smaller runtime for our methods as compared to benchmark methods for equivalent prediction score. Our procedures are implemented in PyTorch (code is available for replication).


title: Extreme data compression while searching for new physics

abstract: Bringing a high-dimensional dataset into science-ready shape is a formidable challenge that often necessitates data compression. Compression has accordingly become a key consideration for contemporary cosmology, affecting public data releases, and reanalyses searching for new physics. However, data compression optimized for a particular model can suppress signs of new physics, or even remove them altogether. We therefore provide a solution for exploring new physics \emph{during} data compression. In particular, we store additional agnostic compressed data points, selected to enable precise constraints of non-standard physics at a later date. Our procedure is based on the maximal compression of the MOPED algorithm, which optimally filters the data with respect to a baseline model. We select additional filters, based on a generalised principal component analysis, which are carefully constructed to scout for new physics at high precision and speed. We refer to the augmented set of filters as MOPED-PC. They enable an analytic computation of Bayesian evidences that may indicate the presence of new physics, and fast analytic estimates of best-fitting parameters when adopting a specific non-standard theory, without further expensive MCMC analysis. As there may be large numbers of non-standard theories, the speed of the method becomes essential. Should no new physics be found, then our approach preserves the precision of the standard parameters. As a result, we achieve very rapid and maximally precise constraints of standard and non-standard physics, with a technique that scales well to large dimensional datasets.


title: Learning to Learn Kernels with Variational Random Features

abstract: In this work, we introduce kernels with random Fourier features in the meta-learning framework to leverage their strong few-shot learning ability. We propose meta variational random features (MetaVRF) to learn adaptive kernels for the base-learner, which is developed in a latent variable model by treating the random feature basis as the latent variable. We formulate the optimization of MetaVRF as a variational inference problem by deriving an evidence lower bound under the meta-learning framework. To incorporate shared knowledge from related tasks, we propose a context inference of the posterior, which is established by an LSTM architecture. The LSTM-based inference network can effectively integrate the context information of previous tasks with task-specific information, generating informative and adaptive features. The learned MetaVRF can produce kernels of high representational power with a relatively low spectral sampling rate and also enables fast adaptation to new tasks. Experimental results on a variety of few-shot regression and classification tasks demonstrate that MetaVRF delivers much better, or at least competitive, performance compared to existing meta-learning alternatives.


title: Discovery of rotation axis alignments in Milky Way globular clusters

abstract: There is an increasing number of recent observational results which show that some globular clusters exhibit internal rotation while they travel along their orbital trajectories around the Milky Way center. Based on these findings, we looked for any relationship between the inclination angles of the globular clusters' orbits with respect to the Milky Way plane and those of their rotation. We discovered that the relative inclination, in the sense rotation axis inclination - orbit axis inclination, is a function of the globular cluster's orbit inclination. Rotation and orbit axes are aligned for an inclination of ~ 56deg, while the rotation axis inclination is far from the orbit's one between ~ 20deg and -20deg when the latter increases from 0deg up to 90deg. We further investigated the origin of such a linear relationship and found no correlation with the semimajor axes and eccentricities of the globular clusters' orbits, nor with the internal rotation strength, the globular clusters' sizes, actual and tidally disrupted masses, half-mass relaxation times, among others. The uncovered relationship will impact on the development of numerical simulations of the internal rotation of globular clusters, on our understanding about the interaction of the globular clusters with the Milky Way gravitational field, and on the observational campaigns for increasing the number of studied globular clusters with detected internal rotation.


title: Low-frequency divergence and quantum geometry of the bulk photovoltaic effect in topological semimetals

abstract: We study the low-frequency properties of the bulk photovoltaic effect in topological semimetals. The bulk photovoltaic effect is a nonlinear optical effect that generates DC photocurrents under uniform irradiation, allowed by noncentrosymmetry. It is a promising mechanism for a terahertz photodetection based on topological semimetals. Here, we systematically investigate the low-frequency behavior of the second-order optical conductivity in point-node semimetals. Through symmetry and power-counting analysis, we show that Dirac and Weyl points with tilted cones show the leading low-frequency divergence. In particular, we find new divergent behaviors of the conductivity of Dirac and Weyl points under circularly polarized light, where the conductivity scales as $\omega^{-2}$ and $\omega^{-1}$ near the gap-closing point in two and three dimensions, respectively. We provide a further perspective on the low-frequency bulk photovoltaic effect by revealing the complete quantum geometric meaning of the second-order optical conductivity tensor. The bulk photovoltaic effect has two origins, which are the transition of electron position and the transition of electron velocity during the optical excitation, and the resulting photocurrents are respectively called the shift current and the injection current. Based on an analysis of two-band models, we show that the injection current is controlled by the quantum metric and Berry curvature, whereas the shift current is governed by the Christoffel symbols near the gap-closing points in semimetals. Finally, for further demonstrations of our theory beyond simple two-band models, we perform first-principles calculations on magnetic Dirac semimetal MnGeO$_3$and Weyl semimetal PrGeAl. Our work brings out new insights into the structure of nonlinear optical responses as well as for the design of semimetal-based terahertz photodetectors.


title: Wave-particle interactions in a long traveling wave tube with upgraded helix

abstract: We investigate the interaction of electromagnetic waves and electron beams in a 4 meters long traveling wave tube (TWT). The device is specially designed to simulate beam-plasma experiments without appreciable noise. This TWT presents an upgraded slow wave structure (SWS) that results in more precise measurements and makes new experiments possible. We introduce a theoretical model describing wave propagation through the SWS and validated by the experimental dispersion relation, impedance, phase and group velocities. We analyze nonlinear effects arising from the beam-wave interaction, such as the modulation of the electron beam and the wave growth and saturation process. When the beam current is low, the wave growth coefficient and saturation amplitude follow the linear theory predictions. However, for high values of current, nonlinear space charge effects become important and these parameters deviate from the linear predictions, tending to a constant value. After saturation, we also observe trapping of the beam electrons, which alters the wave amplitude along the TWT.


title: The cost of approximate controllability of heat equation with dynamical boundary conditions

abstract: We consider the heat equation with dynamic bounary conditions involving gradient terms in a bounded domain. In this paper we study the cost of approximate controllability for this equation. Combining new developed Carleman estimates and some optimization techniques, we obtain explicit bounds of the minimal norm control. We consider the linear and the semilinear cases.


title: A low Hubble Constant from galaxy distribution observations

abstract: An accurate determination of the Hubble constant remains a puzzle in observational cosmology. The possibility of a new physics has emerged with a significant tension between the current expansion rate of our Universe measured from the cosmic microwave background by the Planck satellite and from local methods. In this paper, new tight estimates on this parameter are obtained by considering two data sets from galaxy distribution observations: galaxy cluster gas mass fractions and baryon acoustic oscillation measurements. Priors from the Big Bang nucleosynthesis (BBN) were also considered. By considering the flat $\Lambda$CDM and XCDM models, and the non-flat $\Lambda$CDM model, our main results are: $H_0=65.9^{+1.5}_{-1.5}$ km s$^{-1}$ Mpc$^{-1}$, $H_0=65.9^{+4.4}_{-4.0}$ km s$^{-1}$ Mpc$^{-1}$ and $H_0=64.3^{+ 4.5}_{- 4.4}$ km s$^{-1}$ Mpc$^{-1}$ in $2\sigma$ c.l., respectively. These estimates are in full agreement with the Planck satellite results. Our analyses in these cosmological scenarios also support a negative value for the deceleration parameter at least in 3$\sigma$ c.l..


title: Estimation of temperature of cosmological apparent horizons: a new approach

abstract: We consider radiation from cosmological apparent horizon in Friedmann-Lemaitre-Robertson-Walker(FLRW) model in a double-null coordinate setting. As the spacetime is dynamic, there is no timelike Killing vector, instead we have Kodama vector which acts as dynamical time. We construct the positive frequency modes of the Kodama vector across the horizon. The conditional probability that a signal reaches the central observer when it is crossing from the outside gives the temperature associated with the horizon.


title: A system to test 2D optoelectronic devices in high vacuum

abstract: The exploration of electronic and optoelectronic properties of two-dimensional (2D) materials has become one of the most attractive line of research since the isolation of graphene. Such 'all-surface materials' present a strong sensitivity to environmental conditions and thus characterization of the devices based on these materials usually requires measurement systems operating in high-vacuum. However, conventional optoelectronic probe-station testing systems are are not compatible with high vacuum operation and vacuum-compatible versions are rather expensive. Here, we present a high-vacuum system specifically designed to test electronic and optoelectronic devices based on 2D materials. This system can be implemented with low budget and it is mostly based on the assembly of commercially available standard vacuum and optic components. Despite the simplicity of this system we demonstrate full capabilities to characterize optoelectronic devices in a broad range of wavelengths with fast pumping/venting speed and possibility of modulating the device temperature (room temperature to ~150deg).


title: Data Driven Prediction Architecture for Autonomous Driving and its Application on Apollo Platform

abstract: Autonomous Driving vehicles (ADV) are on road with large scales. For safe and efficient operations, ADVs must be able to predict the future states and iterative with road entities in complex, real-world driving scenarios. How to migrate a well-trained prediction model from one geo-fenced area to another is essential in scaling the ADV operation and is difficult most of the time since the terrains, traffic rules, entities distributions, driving/walking patterns would be largely different in different geo-fenced operation areas. In this paper, we introduce a highly automated learning-based prediction model pipeline, which has been deployed on Baidu Apollo self-driving platform, to support different prediction learning sub-modules' data annotation, feature extraction, model training/tuning and deployment. This pipeline is completely automatic without any human intervention and shows an up to 400\% efficiency increase in parameter tuning, when deployed at scale in different scenarios across nations.


title: Bounds on the Ricci curvature and solutions to the Einstein equations for weighted graphs

abstract: This is a preliminary study of the equation of motion of Euclidean classical gravity on a graph, based on the Lin-Lu-Yau Ricci curvature on graphs. We observe that the constant edge weights configuration gives the unique solution on an infinite tree w.r.t. the asymptotically constant boundary condition. We study the minimum and maximum of the action w.r.t. certain boundary conditions, on several types of graphs of interest. We also exhibit a new class of solutions to the equations of motion on the infinite regular tree.


title: Quantum-classical correspondence for gl(1|1) supersymmetric Gaudin magnet with boundary

abstract: We extend duality between the quantum integrable Gaudin models with boundary and the classical Calogero-Moser systems associated with root systems of classical Lie algebras $B_N$, $C_N$, $D_N$ to the case of supersymmetric ${\rm gl}(m|n)$ Gaudin models with $m+n=2$. Namely, we show that the spectra of quantum Hamiltonians for all such magnets being identified with the classical particles velocities provide the zero level of the classical action variables.


title: Light-wave Control of Non-equilibrium Correlated States using Quantum Femtosecond Magnetism and Time-Periodic Modulation of Coherent Transport

abstract: Lightwave quantum electronics utilizes the oscillating carrier wave of intense laser fields to control quantum materials properties. Using quantum kinetic equations of motion, we describe lightwave-driven nonlinear quantum transport of electronic spin and charge with simultaneous quantum fluctuations of non-collinear local spins. During cycles of field oscillations, spin-charge inter-atomic quantum excitations trigger non-adiabatic time evolution of an antiferromagnetic insulator state into a metallic non-equilibrium state with transient magnetization. Lightwave modulation of electronic hopping changes the energy landscape and establishes a non-thermal pathway to laser-induced transitions in correlated systems with strong local magnetic exchange interactions.


title: Quantum Algorithm for Smoothed Particle Hydrodynamics

abstract: We present a quantum computing algorithm for the smoothed particle hydrodynamics (SPH) method. We use a normalization procedure to encode the SPH operators and domain discretization in a quantum register. We then perform the SPH summation via an inner product of quantum registers. Using a one-dimensional function, we test the approach in a classical sense for the kernel sum and first and second derivatives of a one-dimensional function, using both the Gaussian and Wendland kernel functions, and compare various register sizes against analytical results. Error convergence is exponentially fast in the number of qubits. We extend the method to solve the one-dimensional advection and diffusion partial differential equations, which are commonly encountered in fluids simulations. This work provides a foundation for a more general SPH algorithm, eventually leading to highly efficient simulations of complex engineering problems on gate-based quantum computers.


title: Generalized Cline's formula for G-Drazin inverse in a ring

abstract: In this paper, we give a generalized Cline's formula for the generalized Drazin inverse. Let $R$ be a ring, and let $a,b,c,d\in R$ satisfying $$\begin{array}{c} (ac)^2 = (db)(ac), (db)^2 = (ac)(db);\\ b(ac)a = b(db)a, c(ac)d = c(db)d.\end{array}$$ Then $ac\in R^d$ if and only if $bd\in R^d$. In this case, $(bd)^d = b((ac)^d)^2 d$. We also present generalized Cline's formulas for Drazin and group inverses. Some weaker conditions in a Banach algebra are also investigated. These extend the main results of Cline's formula on g-Drazin inverse of Liao, Chen and Cui (Bull. Malays. Math. Soc., 37 (2014), 37-42), Lian and Zeng (Turk. J. Math., 40 (2016), 161-165) and Miller and Zguitti (Rend. Circ. Mat. Palermo, II. Ser., 67 (2018), 105-114). As an application, new common spectral property of bounded linear operators over Banach spaces are obtained.


title: Backdoor Smoothing: Demystifying Backdoor Attacks on Deep Neural Networks

abstract: Backdoor attacks mislead machine-learning models to output an attacker-specified class when presented a specific trigger at test time. These attacks require poisoning the training data to compromise the learning algorithm, e.g., by injecting poisoning samples containing the trigger into the training set, along with the desired class label. Despite the increasing number of studies on backdoor attacks and defenses, the underlying factors affecting the success of backdoor attacks, along with their impact on the learning algorithm, are not yet well understood. In this work, we aim to shed light on this issue by unveiling that backdoor attacks induce a smoother decision function around the triggered samples -- a phenomenon which we refer to as \textit{backdoor smoothing}. To quantify backdoor smoothing, we define a measure that evaluates the uncertainty associated to the predictions of a classifier around the input samples. Our experiments show that smoothness increases when the trigger is added to the input samples, and that this phenomenon is more pronounced for more successful attacks. We also provide preliminary evidence that backdoor triggers are not the only smoothing-inducing patterns, but that also other artificial patterns can be detected by our approach, paving the way towards understanding the limitations of current defenses and designing novel ones.


title: Search for Axionlike-Particle-Induced Prompt Gamma-Ray Emission from Extragalactic Core-Collapse Supernovae with the Fermi Large Area Telescope

abstract: During a core-collapse supernova (SN), axionlike particles (ALPs) could be produced through the Primakoff process and subsequently convert into $\gamma$ rays in the magnetic field of the Milky Way. We do not find evidence for such a $\gamma$-ray burst in observations of extragalactic SNe with the Fermi Large Area Telescope (LAT). The SN explosion times are estimated from optical light curves and we find a probability of about $\sim$90% that the LAT observed at least one SN at the time of the core collapse. Under the assumption that at least one SN was contained within the LAT field of view, we exclude photon-ALP couplings $\gtrsim 2.6\times10^{-11}$GeV$^{-1}$ for ALP masses $m_a \lesssim 3\times 10^{-10}$ eV, within a factor of $\sim 5$ of previous limits from SN1987A.


title: Double single-channel Kondo coupling in graphene with Fe molecules

abstract: We study the interaction between graphene and a single-molecule-magnet, [Fe4(L)2(dpm)6]. Focusing on the closest Iron ion in a hollow position with respect to the graphene sheet, we derive a channel selective tunneling Hamiltonian, that couples different d orbitals of the Iron atom to precise independent combinations of sublattice and valley degrees of freedom of the electrons in graphene. When looking at the spin-spin interaction between the molecule and the graphene electrons, close to the Dirac point the channel selectivity results in a channel decoupling of the Kondo interaction, with two almost independent Kondo systems weakly interacting among themselves. The formation of magnetic moments and the development of a full Kondo effect depends on the charge state of the graphene layer.


title: Tree/Endofunction Bijections and Concentration Inequalities

abstract: We demonstrate a method for proving precise concentration inequalities in uniformly random trees on $n$ vertices, where $n\geq1$ is a fixed positive integer. The method uses a bijection between mappings $f\colon\{1,\ldots,n\}\to\{1,\ldots,n\}$ and doubly rooted trees on $n$ vertices. The main application is a concentration inequality for the number of vertices connected to an independent set in a uniformly random tree, which is then used to prove partial unimodality of its independent set sequence. So, we give probabilistic arguments for inequalities that often use combinatorial arguments.


title: Designing Experiments: Student Learning Experience and Behaviour in Undergraduate Physics Laboratories

abstract: We investigated physics students' learning experience and behaviour in a second-year laboratory by analyzing transcribed audio recordings of laboratory sessions. One student group was given both a problem and procedure and asked to analyze and explain their results. Another was provided with only the problem and asked to design and execute the experiment, interpret the data, and draw conclusions. These two approaches involved different levels of student inquiry and they have been described as guided and open inquiry respectively. The latter gave students more opportunities to practice "designing experiments," one of the six major learning outcomes in the recommendations for the undergraduate physics laboratory curriculum by the American Association of Physics Teachers (AAPT). Qualitative analysis was performed of the audio transcripts to identify emergent themes and it was augmented by quantitative analysis for a richer understanding of students' experiences. An important finding is that significant improvements can be made to undergraduate laboratories impacting both student learning experience and behaviour by increasing the level of inquiry in laboratory experiments. This is most easily achieved by requiring students to design their own experimental procedures.


title: On next to soft corrections to Drell-Yan and Higgs Boson productions

abstract: We present a framework that resums threshold enhanced large logarithms to all orders in perturbation theory for the production of a pair of leptons in Drell-Yan process and of Higgs boson in gluon fusion as well as in bottom quark annihilation. We restrict ourselves to contributions from diagonal partonic channels. These logarithms include the distributions $((1-z)^{-1} \ln^i(1-z))_+$ resulting from soft plus virtual (SV) and the logarithms $\ln^i(1-z)$ from next-to-SV (NSV) contributions. We use collinear factorisation and renormalisation group invariance to achieve this. The former allows one to define a Soft-Collinear (SC) function which encapsulates soft and collinear dynamics of the perturbative results to all orders in strong coupling constant. The logarithmic structure of these results are governed by universal infrared anomalous dimensions and process dependent functions of Sudakov differential equation that the SC satisfies. The solution to the differential equation is obtained by proposing an all-order ansatz in dimensional regularisation, owing to several state-of-the-art perturbative results available to third order. The $z$ space solutions thus obtained provide an integral representation to sum up large logarithms originating from both soft and collinear configurations, conveniently in Mellin $N$ space. We show that in $N$ space, tower of logarithms $a_s^n/N^\alpha \ln^{2n-\alpha} (N), a_s^n/N^\alpha \ln^{2n-1-\alpha}(N) \cdots $ etc for $\alpha =0,1$ are summed to all orders in $a_s$.


title: Image-Based Model Predictive Control via Dynamic Mode Decomposition

abstract: We present a data-driven model predictive control (MPC) framework for systems with high state-space dimensionalities. This work is motivated by the need to exploit sensor data that appears in the form of images (e.g., 2D or 3D spatial fields reported by thermal cameras). We propose to use dynamic mode decomposition (DMD) to directly build a low-dimensional model from image data and we use such model to obtain a tractable MPC controller. We demonstrate the scalability of this approach (which we call DMD-MPC) by using a 2D thermal diffusion system. Here, we assume that the evolution of the thermal field is captured by 50x50 pixel images, which results in a 2500-dimensional state-space. We show that the dynamics of this high-dimensional space can be accurately predicted by using a 40-dimensional DMD model and we show that the field can be manipulated satisfactorily by using an MPC controller that embeds the low-dimensional DMD model. We also show that the DMD-MPC controller significantly outperforms a standard MPC controller that uses data from a finite set of spatial locations (proxy locations) to manipulate the high-dimensional thermal field. This comparison illustrates the value of information embedded in image data.


title: Deep Reinforcement Learning for Electric Transmission Voltage Control

abstract: Today, human operators primarily perform voltage control of the electric transmission system. As the complexity of the grid increases, so does its operation, suggesting additional automation could be beneficial. A subset of machine learning known as deep reinforcement learning (DRL) has recently shown promise in performing tasks typically performed by humans. This paper applies DRL to the transmission voltage control problem, presents open-source DRL environments for voltage control, proposes a novel modification to the "deep Q network" (DQN) algorithm, and performs experiments at scale with systems up to 500 buses. The promise of applying DRL to voltage control is demonstrated, though more research is needed to enable DRL-based techniques to consistently outperform conventional methods.


title: Weak variable step-size schemes for stochastic differential equations based on controlling conditional moments

abstract: We address the weak numerical solution of stochastic differential equations driven by independent Brownian motions (SDEs for short). This paper develops a new methodology to design adaptive strategies for determining automatically the step-sizes of the numerical schemes that compute the mean values of smooth functions of the solutions of SDEs. First, we introduce a general method for constructing variable step-size weak schemes for SDEs, which is based on controlling the match between the first conditional moments of the increments of the numerical integrator and the ones corresponding to an additional weak approximation. To this end, we use certain local discrepancy functions that do not involve sampling random variables. Precise directions for designing suitable discrepancy functions and for selecting starting step-sizes are given. Second, we introduce a variable step-size Euler scheme, together with a variable step-size second order weak scheme via extrapolation. Finally, numerical simulations are presented to show the potential of the introduced variable step-size strategy and the adaptive scheme to overcome known instability problems of the conventional fixed step-size schemes in the computation of diffusion functional expectations.


title: Is deep learning necessary for simple classification tasks?

abstract: Automated machine learning (AutoML) and deep learning (DL) are two cutting-edge paradigms used to solve a myriad of inductive learning tasks. In spite of their successes, little guidance exists for when to choose one approach over the other in the context of specific real-world problems. Furthermore, relatively few tools exist that allow the integration of both AutoML and DL in the same analysis to yield results combining both of their strengths. Here, we seek to address both of these issues, by (1.) providing a head-to-head comparison of AutoML and DL in the context of binary classification on 6 well-characterized public datasets, and (2.) evaluating a new tool for genetic programming-based AutoML that incorporates deep estimators. Our observations suggest that AutoML outperforms simple DL classifiers when trained on similar datasets for binary classification but integrating DL into AutoML improves classification performance even further. However, the substantial time needed to train AutoML+DL pipelines will likely outweigh performance advantages in many applications.


title: Bandits with Partially Observable Confounded Data

abstract: We study linear contextual bandits with access to a large, confounded, offline dataset that was sampled from some fixed policy. We show that this problem is closely related to a variant of the bandit problem with side information. We construct a linear bandit algorithm that takes advantage of the projected information, and prove regret bounds. Our results demonstrate the ability to take advantage of confounded offline data. Particularly, we prove regret bounds that improve current bounds by a factor related to the visible dimensionality of the contexts in the data. Our results indicate that confounded offline data can significantly improve online learning algorithms. Finally, we demonstrate various characteristics of our approach through synthetic simulations.


title: Preconditioned accelerated gradient descent methods for locally Lipschitz smooth objectives with applications to the solution of nonlinear PDEs

abstract: We develop a theoretical foundation for the application of Nesterov's accelerated gradient descent method (AGD) to the approximation of solutions of a wide class of partial differential equations (PDEs). This is achieved by proving the existence of an invariant set and exponential convergence rates when its preconditioned version (PAGD) is applied to minimize locally Lipschitz smooth, strongly convex objective functionals. We introduce a second-order ordinary differential equation (ODE) with a preconditioner built-in and show that PAGD is an explicit time-discretization of this ODE, which requires a natural time step restriction for energy stability. At the continuous time level, we show an exponential convergence of the ODE solution to its steady state using a simple energy argument. At the discrete level, assuming the aforementioned step size restriction, the existence of an invariant set is proved and a matching exponential rate of convergence of the PAGD scheme is derived by mimicking the energy argument and the convergence at the continuous level. Applications of the PAGD method to numerical PDEs are demonstrated with certain nonlinear elliptic PDEs using pseudo-spectral methods for spatial discretization, and several numerical experiments are conducted. The results confirm the global geometric and mesh size-independent convergence of the PAGD method, with an accelerated rate that is improved over the preconditioned gradient descent (PGD) method.


title: IDEAL: Inexact DEcentralized Accelerated Augmented Lagrangian Method

abstract: We introduce a framework for designing primal methods under the decentralized optimization setting where local functions are smooth and strongly convex. Our approach consists of approximately solving a sequence of sub-problems induced by the accelerated augmented Lagrangian method, thereby providing a systematic way for deriving several well-known decentralized algorithms including EXTRA arXiv:1404.6264 and SSDA arXiv:1702.08704. When coupled with accelerated gradient descent, our framework yields a novel primal algorithm whose convergence rate is optimal and matched by recently derived lower bounds. We provide experimental results that demonstrate the effectiveness of the proposed algorithm on highly ill-conditioned problems.


title: Analytical investigation of wave absorption by a rotating black hole analogue

abstract: Perturbations in a draining vortex can be described analytically in terms of confluent Heun functions. In the context of analogue models of gravity in ideal fluids, we investigate analytically the absorption length of waves in a draining bathtub, a rotating black hole analogue, using confluent Heun functions. We compare our analytical results with the corresponding numerical ones, obtaining excellent agreement.


title: Hairy black holes by gravitational decoupling

abstract: Black holes with hair represented by generic fields surrounding the central source of the vacuum Schwarzschild metric are examined under the minimal set of requirements consisting of i) the existence of a well defined event horizon and ii) the strong or dominant energy condition for the hair outside the horizon. We develop our analysis by means of the gravitational decoupling approach. We find that trivial deformations of the seed Schwarzschild vacuum preserve the energy conditions and provide a new mechanism to evade the no-hair theorem based on a primary hair associated with the charge generating these transformations. Under the above conditions i) and ii), this charge consistently increases the entropy from the minimum value given by the Schwarzschild geometry. As a direct application, we find a non-trivial extension of the Reissner-Nordstrom black hole showing a surprisingly simple horizon. Finally, the non-linear electrodynamics generating this new solution is fully specified.


title: Generalized Jacobson's lemma in a Banach algebra

abstract: Let A be a Banach algebra, and let a; b; c 2 A satisfying a(ba)^2 = abaca = acaba = (ac)^2a: We prove that 1 - ba\in A^d if and only if 1 - ac \in A^d. In this case, (1-ac)^d =1-a(1-ba)^{\pi}(1-\alpha(1+ba))^{-1}bac (1+ac)+a((1-ba)^d)bac. This extends the main result on g-Drazin inverse of Corach (Comm. Algebra, 41(2013), 520{531).


title: Every 7-Dimensional Abelian Variety over the p-adic Numbers has a Reducible $\ell$-adic Galois Representation

abstract: Let $K$ be a complete, discretely valued field with finite residue field and $G_K$ its absolute Galois group. The subject of this note is the study of the set of positive integers $d$ for which there exists an absolutely irreducible $\ell$-adic representation of $G_K$ of dimension $d$ with rational traces on inertia. Our main result is that non-Sophie Germain primes are not in this set when the residue characteristic of $K$ is $> 3$. The result stated in the title is a special case.


title: Eikonal amplitudes and non-global logarithms from the BMS equation

abstract: The Banfi-Marchesini-Smye (BMS) equation accounts for non-global logarithms to all orders in perturbation theory in the large-Nc approximation. We show that the squared amplitudes for the emission of soft energy-ordered gluons are correctly embedded in this equation, and explicitly verify that they coincide with those derived in our previous work in the large-Nc limit up to sixth order in the strong coupling. We perform analytical calculations for the non-global logarithms up to fourth order for the specific hemisphere mass distribution in e+ e- collisions, thus confirming our previous semi-numerical results. We show that the solution to the BMS equation may be cast into a product of an infinite number of exponentials each of which resums a class of Feynman diagrams that manifest a symmetry pattern, and explicitly carry out the computation of the first of these exponentials.


title: A Bayesian response-adaptive dose finding and comparative effectiveness trial

abstract: Aims: Combinations of treatments can offer additional benefit over the treatments individually. However, trials of these combinations are lower priority than the development of novel therapies, which can restrict funding, timelines and patient availability. This paper develops a novel trial design to facilitate the evaluation of novel combination therapies that combines elements of phase II and phase III trials. Methods: This trial uses response adaptive randomisation to increase the information collected about successful novel drug combinations and Bayesian dose-response modelling to undertake a comparative-effectiveness analysis for the most successful dose combination against a relevant comparator. We used simulation methods to evaluate the probability of selecting the correct optimal dose combination, the operating characteristics and predictive power of this design for a trial in pain management and sedation in paediatric emergency departments. Results: With 410 participants, 5 interim updates of the randomisation ratio and a probability of effectiveness of 0.93, 0.88 and 0.83 for the three dose combinations, we have an 83% chance of randomising the largest number of patients to the drug with the highest probability of effectiveness. Based on this adaptive randomisation procedure, the comparative effectiveness analysis has a type I error of less than 5% and a 93% chance of correctly concluding non-inferiority when the probability of effectiveness for the optimal combination therapy is 0.9. In this case, the trial has a 77% chance of meeting its dual aims of dose finding and comparative effectiveness. Finally, the Bayesian predictive power of the trial is over 90%. Conclusion: The proposed trial has high potential to meet the dual study objectives within a feasible level of recruitment, minimising the administrative burden and recruitment time for a trial.


title: Gaze estimation problem tackled through synthetic images

abstract: In this paper, we evaluate a synthetic framework to be used in the field of gaze estimation employing deep learning techniques. The lack of sufficient annotated data could be overcome by the utilization of a synthetic evaluation framework as far as it resembles the behavior of a real scenario. In this work, we use U2Eyes synthetic environment employing I2Head datataset as real benchmark for comparison based on alternative training and testing strategies. The results obtained show comparable average behavior between both frameworks although significantly more robust and stable performance is retrieved by the synthetic images. Additionally, the potential of synthetically pretrained models in order to be applied in user's specific calibration strategies is shown with outstanding performances.


title: The imprint of dark subhaloes on the circumgalactic medium

abstract: The standard model of cosmology, the LCDM model, robustly predicts the existence of a multitude of dark matter 'subhaloes' around galaxies like the Milky Way. A wide variety of observations have been proposed to look for the gravitational effects such subhaloes would induce in observable matter. Most of these approaches pertain to the stellar or cool gaseous phases of matter. Here we propose a new approach, which is to search for the perturbations that such dark subhaloes would source in the warm/hot circumgalactic medium (CGM) around normal galaxies. With a combination of analytic theory, carefully-controlled high-resolution idealised simulations, and full cosmological hydrodynamical simulations (the ARTEMIS simulations), we calculate the expected signal and how it depends on important physical parameters (subhalo mass, CGM temperature, and relative velocity). We find that dark subhaloes enhance both the local CGM temperature and density and, therefore, also the pressure. For the pressure and density, the fluctuations can vary in magnitude from tens of percent (for subhaloes with M_sub=10^10 Msun) to a few percent (for subhaloes with M_sub=10^8 Msun), although this depends strongly on the CGM temperature. The subhaloes also induce fluctuations in the velocity field ranging in magnitude from a few km/s up to 25 km/s. We propose that X-ray, Sunyaev-Zel'dovich effect, radio dispersion measure, and quasar absorption line observations can be used to measure these fluctuations and place constraints on the abundance and distribution of dark subhaloes, thereby placing constraints on the nature of dark matter.


title: Non-Convex SGD Learns Halfspaces with Adversarial Label Noise

abstract: We study the problem of agnostically learning homogeneous halfspaces in the distribution-specific PAC model. For a broad family of structured distributions, including log-concave distributions, we show that non-convex SGD efficiently converges to a solution with misclassification error $O(\opt)+\eps$, where $\opt$ is the misclassification error of the best-fitting halfspace. In sharp contrast, we show that optimizing any convex surrogate inherently leads to misclassification error of $\omega(\opt)$, even under Gaussian marginals.


title: Faster DBSCAN via subsampled similarity queries

abstract: DBSCAN is a popular density-based clustering algorithm. It computes the $\epsilon$-neighborhood graph of a dataset and uses the connected components of the high-degree nodes to decide the clusters. However, the full neighborhood graph may be too costly to compute with a worst-case complexity of $O(n^2)$. In this paper, we propose a simple variant called SNG-DBSCAN, which clusters based on a subsampled $\epsilon$-neighborhood graph, only requires access to similarity queries for pairs of points and in particular avoids any complex data structures which need the embeddings of the data points themselves. The runtime of the procedure is $O(sn^2)$, where $s$ is the sampling rate. We show under some natural theoretical assumptions that $s \approx \log n/n$ is sufficient for statistical cluster recovery guarantees leading to an $O(n\log n)$ complexity. We provide an extensive experimental analysis showing that on large datasets, one can subsample as little as $0.1\%$ of the neighborhood graph, leading to as much as over 200x speedup and 250x reduction in RAM consumption compared to scikit-learn's implementation of DBSCAN, while still maintaining competitive clustering performance.


title: Percolation in random graphs with higher-order clustering

abstract: Percolation theory can be used to describe the structural properties of complex networks using the generating function formulation. This mapping assumes that the network is locally tree-like and does not contain short-range loops between neighbours. In this paper we use the generating function formulation to examine clustered networks that contain simple cycles and cliques of any order. We use the natural generalisation to the Molloy-Reed criterion for these networks to describe their critical properties and derive an analytical description of the size of the giant component, providing solutions for Poisson and power-law networks. We find that networks comprising larger simple cycles behave increasingly more tree-like. Conversley, clustering comprised of larger cliques increasingly deviate from the tree-like solution, although the behaviour is strongly dependent on the degree-assortativity.


title: deepSIP: Linking Type Ia Supernova Spectra to Photometric Quantities with Deep Learning

abstract: We present {\tt deepSIP} (deep learning of Supernova Ia Parameters), a software package for measuring the phase and -- for the first time using deep learning -- the light-curve shape of a Type Ia supernova (SN~Ia) from an optical spectrum. At its core, {\tt deepSIP} consists of three convolutional neural networks trained on a substantial fraction of all publicly-available low-redshift SN~Ia optical spectra, onto which we have carefully coupled photometrically-derived quantities. We describe the accumulation of our spectroscopic and photometric datasets, the cuts taken to ensure quality, and our standardised technique for fitting light curves. These considerations yield a compilation of 2754 spectra with photometrically characterised phases and light-curve shapes. Though such a sample is significant in the SN community, it is small by deep-learning standards where networks routinely have millions or even billions of free parameters. We therefore introduce a data-augmentation strategy that meaningfully increases the size of the subset we allocate for training while prioritising model robustness and telescope agnosticism. We demonstrate the effectiveness of our models by deploying them on a sample unseen during training and hyperparameter selection, finding that Model~I identifies spectra that have a phase between $-10$ and 18\,d and light-curve shape, parameterised by $\Delta m_{15}$, between 0.85 and 1.55\,mag with an accuracy of 94.6\%. For those spectra that do fall within the aforementioned region in phase--$\Delta m_{15}$ space, Model~II predicts phases with a root-mean-square error (RMSE) of 1.00\,d and Model~III predicts $\Delta m_{15}$ values with an RMSE of 0.068\,mag.


title: Deep Convolutional Likelihood Particle Filter for Visual Tracking

abstract: We propose a novel particle filter for convolutional-correlation visual trackers. Our method uses correlation response maps to estimate likelihood distributions and employs these likelihoods as proposal densities to sample particles. Likelihood distributions are more reliable than proposal densities based on target transition distributions because correlation response maps provide additional information regarding the target's location. Additionally, our particle filter searches for multiple modes in the likelihood distribution, which improves performance in target occlusion scenarios while decreasing computational costs by more efficiently sampling particles. In other challenging scenarios such as those involving motion blur, where only one mode is present but a larger search area may be necessary, our particle filter allows for the variance of the likelihood distribution to increase. We tested our algorithm on the Visual Tracker Benchmark v1.1 (OTB100) and our experimental results demonstrate that our framework outperforms state-of-the-art methods.


title: First-Order Methods for Large-Scale Market Equilibrium Computation

abstract: Market equilibrium is a solution concept with many applications such as digital ad markets, fair division, and resource sharing. For many classes of utility functions, equilibria can be captured by convex programs. We develop simple first-order methods suitable for solving these convex programs for large-scale markets. We focus on three practically-relevant utility classes: linear, quasilinear, and Leontief utilities. Using structural properties of market equilibria under each utility class, we show that the corresponding convex programs can be reformulated as optimization of a structured smooth convex function over a polyhedral set, for which projected gradient achieves linear convergence. To do so, we utilize recent linear convergence results under weakened strong-convexity conditions, and further refine the relevant constants in existing convergence results. Then, we show that proximal gradient (a generalization of projected gradient) with a practical linesearch scheme achieves linear convergence under the Proximal-PL condition, a recently developed error bound condition for convex composite problems. For quasilinear utilities, we show that Mirror Descent applied to a new convex program achieves sublinear last-iterate convergence and yields a form of Proportional Response dynamics, an elegant, interpretable algorithm for computing market equilibria originally developed for linear utilities. Numerical experiments show that Proportional Response dynamics is highly efficient for computing approximate market equilibria, while projected gradient with linesearch can be much faster when higher-accuracy solutions are needed.


title: Curvature of planar aesthetic curves

abstract: Farin proposed a method for designing Bezier curves with monotonic curvature and torsion. Such curves are relevant in design due to their aesthetic shape. The method relies on applying a matrix M to the first edge of the control polygon of the curve in order to obtain by iteration the remaining edges. With this method, sufficient conditions on the matrix $M$ are provided, which lead to the definition of Class A curves, generalising a previous result by Mineur et al for plane curves with M being the composition of a dilatation and a rotation. However, Cao and Wang have shown counterexamples for such conditions. In this paper, we revisit Farin's idea of using the subdivision algorithm to relate the curvature at every point of the curve to the curvature at the initial point in order to produce a closed formula for the curvature of planar curves in terms of the eigenvalues of the matrix M and the seed vector for the curve, the first edge of the control polygon. Moreover, we give new conditions in order to produce planar curves with monotonic curvature. The main difference is that we do not require our conditions on the eigenvalues to be preserved under subdivision of the curve. This facilitates giving a unified derivation of the existing results and obtain more general results in the planar case.


title: New models with independent dynamical torsion and nonmetricity fields

abstract: We propose a gravitational model which allows the independent dynamical behaviour of the torsion and nonmetricity fields to be displayed in the framework of Metric-Affine gauge theory of gravity. For this task, we derive a new exact black hole solution referred to this model which extends the role of torsion of the main well-known exact solutions based on Weyl-Cartan geometry and constitutes the first known isolated gravitational system characterized by a metric tensor with independent spin and dilation charges.


title: Solar System Tests in Modified Teleparallel Gravity

abstract: In this paper, we study different Solar System tests in a modified Teleparallel gravity theory based on an arbitrary function $f(T,B)$ which depends on the scalar torsion $T$ and the boundary term $B$. To do this, we first find new perturbed spherically symmetric solutions around Schwarzschild for different power-law forms of the arbitrary Lagrangian. Then, for each model we calculated the photon sphere, perihelion shift, deflection of light, Cassini experiment, Shapiro delay and the gravitational redshift. Finally, we confront these computations with different known experiments from these Solar System tests to put different bounds on the mentioned models. We then conclude that $f(T,B)$ is compatible with these Solar System experiments with a wide range of parameters which are relevant for cosmology.


title: Dimensional Reduction of Helium-4 Inside Argon Plated MCM-41 Nanopores

abstract: The angstrom-scale coherence length describing the superfluid wavefunction of helium-4 at low temperatures has prevented its preparation in a truly one-dimensional geometry. Mesoporous ordered silica-based structures, such as the molecular sieve MCM-41, offer a promising avenue towards physical confinement, but the minimal pore diameters that can be chemically synthesized have proven to be too large to reach the quasi-one-dimensional limit. We present an active nano-engineering approach to this problem by pre-plating MCM-41 with a single, well controlled layer of Ar gas before filling the pores with helium. The structure inside the pore is investigated via experimental adsorption isotherms and neutron scattering measurements that are in agreement with large scale quantum Monte Carlo simulations. The results demonstrate angstrom and kelvin scale tunability of the effective confinement potential experienced by 4He atoms inside the MCM-41, with the Ar layer reducing the diameter of the confining media into a regime where a number of solid layers surround a one-dimensional quantum liquid.


title: An Unsupervised Information-Theoretic Perceptual Quality Metric

abstract: Tractable models of human perception have proved to be challenging to build. Hand-designed models such as MS-SSIM remain popular predictors of human image quality judgements due to their simplicity and speed. Recent modern deep learning approaches can perform better, but they rely on supervised data which can be costly to gather: large sets of class labels such as ImageNet, image quality ratings, or both. We combine recent advances in information-theoretic objective functions with a computational architecture informed by the physiology of the human visual system and unsupervised training on pairs of video frames, yielding our Perceptual Information Metric (PIM). We show that PIM is competitive with supervised metrics on the recent and challenging BAPPS image quality assessment dataset and outperforms them in predicting the ranking of image compression methods in CLIC 2020. We also perform qualitative experiments using the ImageNet-C dataset, and establish that PIM is robust with respect to architectural details.