Search

We describe observations of low energy excess (LEE) events (background events observed in all light dark matter direct detection calorimeters) and noise in a two-channel silicon athermal phonon detector with 375 meV baseline energy resolution. We measure two distinct LEE populations: ``shared'' multichannel events with a pulse shape consistent with athermal phonon events, and sub-eV events which couple nearly exclusively to a single channel with a significantly faster pulse shape. These ``singles'' are consistent with events occurring within the aluminum athermal phonon collection fins. Similarly, our measured detector noise is higher than the theoretical expectation. Measured noise can be split into an uncorrelated component, consistent with shot noise from small energy depositions within the athermal phonon sensor itself, and a correlated component, consistent with shot noise from energy depositions within the silicon crystal's phonon system., Comment: 6 pages, 5 figures

Mature radio galaxies such as M87 belong to a specific subclass of active galaxies (AGN) whose evolution in time endows them with five distinguishing characteristics, including (1) low excitation emission, (2) low star formation rates, (3) high bulge stellar-velocity dispersion, (4) bright stellar nuclei, and (5) weak or nonexistent merger signatures. We show how to understand these seemingly disparate characteristics as originating from the time evolution of powerful radio quasars and describe a new model prediction that tilted accretion disks in AGN are expected to occur in bright quasars but not in other subclasses of AGN. The picture we present should be understood as the most compelling evidence for counter-rotation as a key element in feedback from accreting black holes., Comment: 11 pages, 6 Figures, Published in the Universe journal

The CODEX-$\beta$ apparatus is a demonstrator for the proposed future CODEX-b experiment, a long-lived-particle detector foreseen for operation at IP8 during HL-LHC data-taking. The demonstrator project, intended to collect data in 2025, is described, with a particular focus on the design, construction, and installation of the new apparatus.

The first direct detection of gravitational waves from binary neutron stars on the 17th of August, 2017, (GW170817) heralded the arrival of a new messenger for probing neutron star astrophysics and provided the first constraints on neutron star equation of state from gravitational wave observations. Significant computational effort was expended to obtain these first results and therefore, as observations of binary neutron star coalescence become more routine in the coming observing runs, there is a need to improve the analysis speed and flexibility. Here, we present a rapid approach for inferring the neutron star equation of state based on Normalising Flows. As a demonstration, using the same input data, our approach, ASTREOS, produces results consistent with those presented by the LIGO-Virgo collaboration but requires < 1 sec to generate neutron star equation of state confidence intervals. Furthermore, ASTREOS allows for non-parametric equation of state inference. This rapid analysis will not only facilitate neutron star equation of state studies but can potentially enhance future alerts for electromagnetic follow-up observations of binary neutron star mergers., Comment: 4 + 2 pages, 2 figures

We apply and extend standard tools for void statistics to cosmological simulations that solve Einstein's equations with numerical relativity (NR). We obtain a simulated void catalogue without Newtonian approximations, using a new watershed void finder which operates on fluid-based NR simulations produced with the Einstein Toolkit. We compare and contrast measures of void size and void fraction, and compare radial stacked density profiles to empirically-derived Hamaus-Sutter-Wandelt (HSW) density profiles and profiles based on distance to void boundaries. We recover statistics roughly consistent with Newtonian N-body simulations where such a comparison is meaningful. We study variation of dynamical spatial curvature and local expansion explicitly demonstrating the spatial fluctuations of these quantities in void regions. We find that voids in our simulations expand ~10-30% faster than the global average and the kinetic curvature density parameter in the centre of voids reaches ~60-80%. Within the limits of resolution of the simulations, the results are consistent with the parameters of the Timescape model of cosmological backreaction., Comment: 15 pages, 14 figures

The counter-rotation between black holes and accretion disk configuration was introduced over a decade ago to elucidate the nature of the radio loud/radio-quiet dichotomy and the jet-disk connection, but has since been applied to a plethora of observations across space and time. We briefly review the paradigm in which counter-rotation is key for the triggering of radio galaxies and its observational support, then apply it to a series of observations concerning FR0 radio galaxies. FR0 radio galaxies appear to be radio galaxies in transition, with low-spinning black holes and thus weaker but tilted jets with respect to an earlier radio quasar phase. As a result, FR0 radio galaxies are prescribed to be in an earlier phase of star formation suppression in radio galaxies, compared to a later phase that is unlikely to be less than tens of millions of years in the future if they have enough accretion fuel to evolve into more powerful FRI radio galaxies. FR0 radio galaxies will have a greater or lesser star formation suppression feedback effect depending on how long they live. Tilted jets also enhance stellar velocities in the bulge. Because FR0 jet lengths are of the same order of magnitude as the radius of the stellar bulge, FR0 jets are prescribed to have begun, more or less recently depending on their age, to affect stellar velocity dispersions as well. As a result, they will be associated with dispersion values that tend to be larger than for characteristically non-jetted active galaxies, but smaller than giant radio galaxies such as M87 that have experienced a long-term tilted and more powerful FRI jet. With these ideas it is possible to make a coarse-grained prediction for the slope of the M-{\sigma} plane for FR0 radio galaxies with values between 4 and 8., Comment: 13 pages, 4 Figures, Accepted for publication in the Journal of High Energy Astrophysics

The discovery of the optical counterpart, along with the gravitational waves from GW170817, of the first binary neutron star merger, opened up a new era for multi-messenger astrophysics. Combining the GW data with the optical counterpart, also known as AT2017gfo, classified as a kilonova, has revealed the nature of compact binary merging systems by extracting enriched information about the total binary mass, the mass ratio, the system geometry, and the equation of state. Even though the detection of kilonova brought about a revolution in the domain of multi-messenger astronomy, since there has been only one kilonova from a gravitational wave detected binary neutron star merger event so far, this limits the exact understanding of the origin and propagation of the kilonova. Here, we use a conditional variational autoencoder trained on light curve data from two kilonova models having different temporal lengths, and consequently, generate kilonova light curves rapidly based on physical parameters of our choice with good accuracy. Once trained, the time scale for light curve generation is of the order of a few milliseconds, thus speeding up generating light curves by $1000$ times compared to the simulation. The mean squared error between the generated and original light curves is typically $0.015$ with a maximum of $0.08$ for each set of considered physical parameter; while having a maximum of $\approx0.6$ error across the whole parameter space. Hence, implementing this technique provides fast and reliably accurate results., Comment: 19 pages, 7 figures (3 additional figures in appendix), accepted to ApJ

The monotonic dependence of the outputs of a neural network on some of its inputs is a crucial inductive bias in many scenarios where domain knowledge dictates such behavior. This is especially important for interpretability and fairness considerations. In a broader context, scenarios in which monotonicity is important can be found in finance, medicine, physics, and other disciplines. It is thus desirable to build neural network architectures that implement this inductive bias provably. In this work, we propose a weight-constrained architecture with a single residual connection to achieve exact monotonic dependence in any subset of the inputs. The weight constraint scheme directly controls the Lipschitz constant of the neural network and thus provides the additional benefit of robustness. Compared to currently existing techniques used for monotonicity, our method is simpler in implementation and in theory foundations, has negligible computational overhead, is guaranteed to produce monotonic dependence, and is highly expressive. We show how the algorithm is used to train powerful, robust, and interpretable discriminators that achieve competitive performance compared to current state-of-the-art methods across various benchmarks, from social applications to the classification of the decays of subatomic particles produced at the CERN Large Hadron Collider., Comment: 9 pages, 4 figures, ICLR 2023 final submission

We present a novel fully Bayesian analysis to constrain short gamma-ray burst jet structures associated with cocoon, wide-angle and simple top-hat jet models, as well as the binary neutron star merger rate. These constraints are made given the distance and inclination information from GW170817, observed flux of GRB170817A, observed rate of short gamma-ray bursts detected by Swift, and the neutron star merger rate inferred from LIGO's first and second observing runs. A separate analysis is conducted where a fitted short gamma-ray burst luminosity function is included to provide further constraints. The jet structure models are further constrained using the observation of GW190425 and we find that the assumption that it produced a GRB170817-like short gamma-ray burst that went undetected due to the jet geometry is consistent with previous observations. We find and quantify evidence for low luminosity and wide-angled jet structuring in the short gamma-ray burst population, independently from afterglow observations, with log Bayes factors of $0.45{-}0.55$ for such models when compared to a classical top-hat jet. Slight evidence is found for a Gaussian jet structure model over all others when the fitted luminosity function is provided, producing log Bayes factors of $0.25{-}0.9\pm0.05$ when compared to the other models. However without considering GW190425 or the fitted luminosity function, the evidence favours a cocoon-like model with log Bayes factors of $0.14\pm0.05$ over the Gaussian jet structure. We provide new constraints to the binary neutron star merger rates of $1{-}1300$Gpc$^{-3}$yr$^{-1}$ or $2{-}680$Gpc$^{-3}$yr$^{-1}$ when a fitted luminosity function is assumed., Comment: 18 pages, 8 figures

We present an improved version of the nested sampling algorithm nessai in which the core algorithm is modified to use importance weights. In the modified algorithm, samples are drawn from a mixture of normalising flows and the requirement for samples to be independently and identically distributed (i.i.d.) according to the prior is relaxed. Furthermore, it allows for samples to be added in any order, independently of a likelihood constraint, and for the evidence to be updated with batches of samples. We call the modified algorithm i-nessai. We first validate i-nessai using analytic likelihoods with known Bayesian evidences and show that the evidence estimates are unbiased in up to 32 dimensions. We compare i-nessai to standard nessai for the analytic likelihoods and the Rosenbrock likelihood, the results show that i-nessai is consistent with nessai whilst producing more precise evidence estimates. We then test i-nessai on 64 simulated gravitational-wave signals from binary black hole coalescence and show that it produces unbiased estimates of the parameters. We compare our results to those obtained using standard nessai and dynesty and find that i-nessai requires 2.68 and 13.3 times fewer likelihood evaluations to converge, respectively. We also test i-nessai of an 80-second simulated binary neutron star signal using a Reduced-Order-Quadrature (ROQ) basis and find that, on average, it converges in 24 minutes, whilst only requiring $1.01\times10^6$ likelihood evaluations compared to $1.42\times10^6$ for nessai and $4.30\times10^7$ for dynesty. These results demonstrate the i-nessai is consistent with nessai and dynesty whilst also being more efficient., Comment: 41 pages, 12 figures, 2 tables, published in Machine Learning: Science and Technology, code available at https://doi.org/10.5281/zenodo.8124198, updated with accepted version

Strain-engineering is a powerful means to tune the polar, structural, and electronic instabilities of incipient ferroelectrics. KTaO3 is near a polar instability and shows anisotropic superconductivity in electron-doped samples. Here, we demonstrate growth of high quality KTaO3 thin films by molecular-beam epitaxy. Tantalum was provided by both a suboxide source emanating a TaO2 flux from Ta2O5 contained in a conventional effusion cell as well as an electron-beam-heated tantalum source. Excess potassium and a combination of ozone and oxygen (10 \% O3 + 90 \% O2) were simultaneously supplied with the TaO$_2$ (or tantalum) molecular beams to grow the KTaO$_3$ films. Laue fringes suggest that the films are smooth with an abrupt film/substrate interface. Cross-sectional scanning transmission electron microscopy does not show any extended defects and confirms that the films have an atomically abrupt interface with the substrate. Atomic force microscopy reveals atomic steps at the surface of the grown films. Reciprocal space mapping demonstrates that the films, when sufficiently thin, are coherently strained to the SrTiO$_3$ (001) and GdScO$_3$ (110) substrates.

We report the use of suboxide molecular-beam epitaxy (S-MBE) to grow $\beta$-Ga$_2$O$_3$ at a growth rate of ~1 ${\mu}$m/h with control of the silicon doping concentration from 5x10$^{16}$ to 10$^{19}$ cm$^{-3}$. In S-MBE, pre-oxidized gallium in the form of a molecular beam that is 99.98\% Ga$_2$O, i.e., gallium suboxide, is supplied. Directly supplying Ga2O to the growth surface bypasses the rate-limiting first step of the two-step reaction mechanism involved in the growth of $\beta$-Ga$_2$O$_3$ by conventional MBE. As a result, a growth rate of ~1 ${\mu}$m/h is readily achieved at a relatively low growth temperature (T$_{sub}$ = 525 $^\circ$C), resulting in films with high structural perfection and smooth surfaces (rms roughness of < 2 nm on ~1 ${\mu}$m thick films). Silicon-containing oxide sources (SiO and SiO$_2$) producing an SiO suboxide molecular beam are used to dope the $\beta$-Ga$_2$O$_3$ layers. Temperature-dependent Hall effect measurements on a 1 ${\mu}$m thick film with a mobile carrier concentration of 2.7x10$^{17}$ cm$^{-3}$ reveal a room-temperature mobility of 124 cm$^2$ V$^{-1}$ s$^{-1}$ that increases to 627 cm$^2$ V$^{-1}$ s$^{-1}$ at 76 K; the silicon dopants are found to exhibit an activation energy of 27 meV. We also demonstrate working MESFETs made from these silicon-doped $\beta$-Ga$_2$O$_3$ films grown by S-MBE at growth rates of ~1 ${\mu}$m/h., Comment: 19 pages, 7 figures, 2 tables, 2 pages supplementary materials

The performance of superconducting qubits is degraded by a poorly characterized set of energy sources breaking the Cooper pairs responsible for superconductivity, creating a condition often called ``quasiparticle poisoning". Both superconducting qubits and low threshold dark matter calorimeters have observed excess bursts of quasiparticles or phonons that decrease in rate with time. Here, we show that a silicon crystal glued to its holder exhibits a rate of low-energy phonon events that is more than two orders of magnitude larger than in a functionally identical crystal suspended from its holder in a low-stress state. The excess phonon event rate in the glued crystal decreases with time since cooldown, consistent with a source of phonon bursts which contributes to quasiparticle poisoning in quantum circuits and the low-energy events observed in cryogenic calorimeters. We argue that relaxation of thermally induced stress between the glue and crystal is the source of these events., Comment: 19 pages (main + supplementary), 6 figures. W. A. Page and R. K. Romani contributed equally to this work. Correspondence should be addressed to R. K. Romani

In this Snowmass contribution we present a comprehensive status update on the progress and plans for the proposed CODEX-b detector, intended to search for long-lived particles beyond the Standard Model. We review the physics case for the proposal and present recent progress on optimization strategies for the detector and shielding design, as well as the development of new fast and full simulation frameworks. A summary of the technical design for a smaller demonstrator detector (CODEX-$\beta$) for the upcoming Run~3 of the LHC is also discussed, alongside the road towards realization of the full experiment at the High-Luminosity LHC., Comment: 13 pages, 13 figures, Snowmass contribution

Owing to the computational complexity of electronic structure algorithms running on classical digital computers, the range of molecular systems amenable to simulation remains tightly circumscribed even after many decades of work. Quantum computers hold the promise of transcending such limitations although in the current era the size and noise of these devices militates against significant progress. Here we describe a new and chemically intuitive approach that permits a subdomain of the electronic structure of a molecule to be calculated accurately on a quantum device, while the rest of the molecule is described at a lower level of accuracy using density functional theory running on a classical computer. We demonstrate that our method produces improved results for molecules that cannot be simulated fully on quantum computers but which can be resolved classically at a lower level of approximation. Our algorithm is tunable, so that the size of the quantum simulation can be adjusted to run on available quantum resources. Therefore, as quantum devices become larger, our method will enable increasingly large subdomains to be studied accurately., Comment: 8 pages, 3 figures, Supplementary Information pdf available in anc folder. Updated pre-print includes changes to acknowledgements, and minor changes to the density matrix definition and orbital localization discussion

We settle in the affirmative the Graham-Sloane conjecture., Comment: arXiv admin note: text overlap with arXiv:2111.12812

In the field of multi-messenger astronomy, Bayesian inference is commonly adopted to compare the compatibility of models given the observed data. However, to describe a physical system like neutron star mergers and their associated gamma-ray burst (GRB) events, usually more than ten physical parameters are incorporated in the model. With such a complex model, likelihood evaluation for each Monte Carlo sampling point becomes a massive task and requires a significant amount of computational power. In this work, we perform quick parameter estimation on simulated GRB X-ray light curves using an interpolated physical GRB model. This is achieved by generating a grid of GRB afterglow light curves across the parameter space and replacing the likelihood with a simple interpolation function in the high-dimensional grid that stores all light curves. This framework, compared to the original method, leads to a $\sim$90$\times$ speedup per likelihood estimation. It will allow us to explore different jet models and enable fast model comparison in the future., Comment: 9 pages, 4 figures, accepted to the special issue of Universe, "Waiting for GODOT -- Present and Future of Multi-Messenger Astronomy"

In this paper, we describe the potential of the LHCb experiment to detect Stealth physics. This refers to dynamics beyond the Standard Model that would elude searches that focus on energetic objects or precision measurements of known processes. Stealth signatures include long-lived particles and light resonances that are produced very rarely or together with overwhelming backgrounds. We will discuss why LHCb is equipped to discover this kind of physics at the Large Hadron Collider and provide examples of well-motivated theoretical models that can be probed with great detail at the experiment., Comment: 82 pages, 31 figures. This is the version of the article before peer review or editing, as submitted by an author to Reports on Progress in Physics. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at doi:10.1088/1361-6633/ac4649

We present a novel method for sampling iso-likelihood contours in nested sampling using a type of machine learning algorithm known as normalising flows and incorporate it into our sampler nessai. Nessai is designed for problems where computing the likelihood is computationally expensive and therefore the cost of training a normalising flow is offset by the overall reduction in the number of likelihood evaluations. We validate our sampler on 128 simulated gravitational wave signals from compact binary coalescence and show that it produces unbiased estimates of the system parameters. Subsequently, we compare our results to those obtained with dynesty and find good agreement between the computed log-evidences whilst requiring 2.07 times fewer likelihood evaluations. We also highlight how the likelihood evaluation can be parallelised in nessai without any modifications to the algorithm. Finally, we outline diagnostics included in nessai and how these can be used to tune the sampler's settings., Comment: 20 pages, 11 figures, accepted to PRD, typos corrected, added new figure, other minor corrections, code available at https://doi.org/10.5281/zenodo.4550693