Search

The formation mechanism of globular clusters (GCs) has long been debated by astronomers. It was recently proposed that Supersonically Induced Gas Objects (SIGOs), which formed in the early Universe due to the supersonic relative motion of baryons and dark matter at recombination, could be the progenitors of early globular clusters. In order to become GCs, SIGOs must form stars relatively efficiently despite forming outside of dark matter halos. We investigate the potential for star formation in SIGOs using cosmological hydrodynamic simulations, including the aforementioned relative motions of baryons and dark matter, molecular hydrogen cooling in primordial gas clouds, and including explicit star formation. We find that SIGOs do form stars and that the nascent star clusters formed through this process are accreted by dark matter halos on short timescales (a few hundreds of Myr). Thus, SIGOs may be found as intact substructures within these halos, analogous to many present-day GCs. From this result, we conclude that SIGOs are capable of forming star clusters with similar properties to globular clusters in the early Universe and we discuss their detectablity by upcoming JWST surveys., Comment: 12 pages, 5 figures

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback and magnetic fields, to investigate the effectiveness of different fitting methods for volume and column density probability distribution functions (PDFs). The first method fits a piecewise lognormal and power-law (PL) function to recover PDF parameters such as the PL slope and transition density. The second method fits a polynomial spline function and examines the first and second derivatives of the spline to determine the PL slope and the functional transition density. We demonstrate that fitting a spline allows us to directly determine if the data has multiple PL slopes. The first PL (set by the transition between lognormal and PL function) can also be visualized in the derivatives directly. In general, the two methods produce fits that agree reasonably well for volume density but vary for column density, likely due to the increased statistical noise in column density maps as compared to volume density. We test a well-known conversion for estimating volume density PL slopes from column density slopes and find that the spline method produces a better match (\c{hi}2 of 2.38 vs \c{hi}2 of 5.92), albeit with a significant scatter. Ultimately, we recommend the use of both fitting methods on column density data to mitigate the effects of noise., Comment: 14 pages, 8 figures, Submitted to ApJ. Comments welcome

The emission of $\gamma$-rays after a neutron capture in a cryogenic detector can generate mono-energetic nuclear recoils in the sub-keV regime, of direct interest for the calibration of Dark Matter and Coherent Elastic Neutrino Nucleus Scattering experiments. Here we show that accurate predictions of the nuclear recoil spectra induced by neutron captures require taking into account the interplay between the development in time of the de-excitation $\gamma$-cascade of the target nucleus and that of the associated atomic collisions in matter. We present detailed simulations coupling the FIFRELIN code for the description of the $\gamma$-cascades and the IRADINA code for the modelling of the fast atomic movements in matter. Nuclear recoil spectra are predicted, and made available to the community, for concrete cases of Al$_2$O$_3$, Si, Ge and CaWO$_4$ crystals exposed to a low intensity beam of thermal neutrons. We find that timing effects cause new calibration peaks to emerge in the recoil spectra and also impact the shape of the continuous recoil distribution. We discuss how they could give access to a rich physics program, spanning the accurate study of the response of cryogenic detectors in the sub-keV range, tests of solid state physics simulations and tests of nuclear models.

An unprecedented array of new observational capabilities are starting to yield key constraints on models of the epoch of first light in the Universe. In this Letter we discuss the implications of the UV radiation background at cosmic dawn inferred by recent JWST observations for radio experiments aimed at detecting the redshifted 21-cm hyperfine transition of diffuse neutral hydrogen. Under the basic assumption that the 21-cm signal is activated by the Ly$\alpha$ photon field produced by metal-poor stellar systems, we show that a detection at the low frequencies of the EDGES and SARAS3 experiments may be expected from a simple extrapolation of the declining UV luminosity density inferred at $z\lesssim 14$ from JWST early galaxy data. Accounting for an early radiation excess above the CMB suggests a shallower or flat evolution to simultaneously reproduce low and high-$z$ current UV luminosity density constraints, which cannot be entirely ruled out, given the large uncertainties from cosmic variance and the faint-end slope of the galaxy luminosity function at cosmic dawn. Our findings raise the intriguing possibility that a high star formation efficiency at early times may trigger the onset of intense Ly$\alpha$ emission at redshift $z\lesssim 20$ and produce a cosmic 21-cm absorption signal 200 Myr after the Big Bang., Comment: 8 pages, 3 figures, accepted to ApJL

New-generation direct searches for low mass dark matter feature detection thresholds at energies well below 100 eV, much lower than the energies of commonly used X-ray calibration sources. This requires new calibration sources with sub-keV energies. When searching for nuclear recoil signals, the calibration source should ideally cause mono-energetic nuclear recoils in the relevant energy range. Recently, a new calibration method based on the radiative neutron capture on $^{182}$W with subsequent de-excitation via single $\gamma$-emission leading to a nuclear recoil peak at 112 eV was proposed. The CRESST-III dark matter search operated several CaWO$_{4}$-based detector modules with detection thresholds below 100 eV in the past years. We report the observation of a peak around the expected energy of 112 eV in the data of three different detector modules recorded while irradiated with neutrons from different AmBe calibration sources. We compare the properties of the observed peaks with Geant-4 simulations and assess the prospects of using this for the energy calibration of CRESST-III detectors., Comment: 8 pages, 4 figures; submitted to Phys. Rev. D

The Electron, Proton and Alpha Monitor, EPAM, located at the L1 Position approximately 1-million miles from the earth in the direction of the sun, was designed to detect fluctuations in solar output through counting the numbers of various particles hitting the detector. The EPAM detector is part of an early warning system that can alert the earth to coronal mass ejection events that can damage our electronic grids and satellite equipment. EPAM gives a real-time estimate of changes in the local solar magnetic field directed towards the earth, recorded in the fluctuations of solar particles being ejected. This paper presents an analysis of fluctuations in data taken by the Geological Survey of Israel, GSI, compared to the changes in detected numbers of protons as seen by EPAM. Surprisingly, the GSI and EPAM detectors show an unexpected correlation between the variation in count rate detected by the GSI detectors and an increased numbers of protons seen at EPAM; well above statistical significance of 5-sigma, indicating a non-random connection between the data sets. The statistically significant overlap between data taken by these two detectors, subject to very different conditions, may hint at a Primakoff mechanism whereby exotic particles, e.g. galactic Dark Matter, couple through magnetic fields to both photons and even nuclei. This work builds on an earlier paper on the observations of Radon decay and their implications for particle physics., Comment: N/A

Machine learning techniques are effective for building predictive models because they identify patterns in large datasets. Development of a model for complex real-life problems often stop at the point of publication, proof of concept or when made accessible through some mode of deployment. However, a model in the medical domain risks becoming obsolete as patient demographics, systems and clinical practices change. The maintenance and monitoring of predictive model performance post-publication is crucial to enable their safe and effective long-term use. We will assess the infrastructure required to monitor the outputs of a machine learning algorithm, and present two scenarios with examples of monitoring and updates of models, firstly on a breast cancer prognosis model trained on public longitudinal data, and secondly on a neurodegenerative stratification algorithm that is currently being developed and tested in clinic., Comment: 12 pages, 3 figures

Formal verification of cyber-physical and robotic systems requires that we can accurately model physical quantities that exist in the real-world. The use of explicit units in such quantities can allow a higher degree of rigour, since we can ensure compatibility of quantities in calculations. At the same time, improper use of units can be a barrier to safety and therefore it is highly desirable to have automated sanity checking in physical calculations. In this paper, we contribute a mechanisation of the International System of Quantities (ISQ) and the associated SI unit system in Isabelle/HOL. We show how Isabelle can be used to provide a type system for physical quantities, and automated proof support. Quantities are parameterised by dimension types, which correspond to base vectors, and thus only quantities of the same dimension can be equated. Since the underlying "algebra of quantities" induces congruences on quantity and SI types, specific tactic support is developed to capture these. Our construction is validated by a test-set of known equivalences between both quantities and SI units. Moreover, the presented theory can be used for type-safe conversions between the SI system and others, like the British Imperial System (BIS)., Comment: 10 pages, submitted to ICECCS 2023

We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC-pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion ($\sigma_v$) versus gas surface density ($\Sigma_{\mathrm{mol}}$) to the one observed in local spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For GMC-pixels in simulated mergers, we see a significant increase of factor of 5 - 10 in both $\Sigma_{\mathrm{mol}}$ and $\sigma_v$, which puts these pixels above the trend of PHANGS galaxies in the $\sigma_v$ vs $\Sigma_{\mathrm{mol}}$ diagram. This deviation may indicate that GMCs in the simulated mergers are much less gravitationally bound compared with simulated control galaxies with virial parameter ($\alpha_{\mathrm{vir}}$) reaching 10 - 100. Furthermore, we find that the increase in $\alpha_{\mathrm{vir}}$ happens at the same time as the increase in global star formation rate (SFR), which suggests stellar feedback is responsible for dispersing the gas. We also find that the gas depletion time is significantly lower for high $\alpha_{\mathrm{vir}}$ GMCs during a starburst event. This is in contrast to the simple physical picture that low $\alpha_{\mathrm{vir}}$ GMCs are easier to collapse and form stars on shorter depletion times. This might suggest that some other physical mechanisms besides self-gravity are helping the GMCs in starbursting mergers collapse and form stars., Comment: 22 pages, 11 figures. Accepted to ApJ. Link to animation updated

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high and low density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations., Comment: 22 pages, 13 figures. Submitted to ApJ. Comments welcome

Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the ISM/CGM transition region which also provides the necessary level of refinement at the base of the wind. In the absence of our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the SFR of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model., Comment: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

We present limits on the spin-independent interaction cross section of dark matter particles with silicon nuclei, derived from data taken with a cryogenic calorimeter with 0.35 g target mass operated in the CRESST-III experiment. A baseline nuclear recoil energy resolution of $(1.36\pm 0.05)$ eV$_{\text{nr}}$, currently the lowest reported for macroscopic particle detectors, and a corresponding energy threshold of $(10.0\pm 0.2)$ eV$_{\text{nr}}$ have been achieved, improving the sensitivity to light dark matter particles with masses below 160 MeV/c$^2$ by a factor of up to 20 compared to previous results. We characterize the observed low energy excess, and we exclude noise triggers and radioactive contaminations on the crystal surfaces as dominant contributions., Comment: 8 pages, 5 figures; precised the position of the calibration source in Fig. 1; extended the discussion about the observed energy spectrum; added the DM limit curve to ancillary files. Published in Phys. Rev. D

This paper presents a new analysis of the observations of radon decay in an enclosed environment by the Geological Survey of Israel (GSI) between 2007 and 2012 [4]; for a more complete list of experiments performed by GSI on radon see also references [1-3]. The data exhibit a large peak around local noon followed by an abrupt drop, and by a second peak around 6PM local time. Additionally, there is also a very low amplitude peak occurring before daybreak. The salient features of the GSI radon decay data can be modeled as arising from a change in the radon decay rate, rather than due to the changes in the local concentration of radon (N0). Such a model may provide a clue to long theorized axionic, dark matter, interactions. Finally, new experimentation is suggested that can distinguish between changes in N0 versus changes in decay rate. Should a follow-up experiment show an effect similar to that seen in GSI, this could have significant implications for elementary particle physics.

We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful far-ultraviolet (FUV) diagnostic. Molecular hydrogen (H2) is the most abundant molecule in the universe and a key ingredient for star and planet formation, but is typically not observed directly because its symmetric atomic structure and lack of a dipole moment mean there are no spectral lines at visible wavelengths and few in the infrared. Hyperion uses molecular hydrogen's wealth of FUV emission lines to achieve three science objectives: (1) determining how star formation is related to molecular hydrogen formation and destruction at the boundaries of molecular clouds; (2) determining how quickly and by what process massive star feedback disperses molecular clouds; and (3) determining the mechanism driving the evolution of planet-forming disks around young solar-analog stars. Hyperion conducts this science using a straightforward, highly-efficient, single-channel instrument design. Hyperion's instrument consists of a 48 cm primary mirror, with an f/5 focal ratio. The spectrometer has two modes, both covering 138.5-161.5 nm bandpasses. A low resolution mode has a spectral resolution of R>10,000 with a slit length of 65 arcmin, while the high resolution mode has a spectral resolution of R>50,000 over a slit length of 5 armin. Hyperion occupies a 2 week long, high-earth, Lunar resonance TESS-like orbit, and conducts 2 weeks of planned observations per orbit, with time for downlinks and calibrations. Hyperion was reviewed as Category I, which is the highest rating possible, but was not selected., Comment: Accepted to JATIS, 9 Figures

We demonstrate that two supersoluble complements of an abelian base in a finite split extension are conjugate if and only if, for each prime $p$, a Sylow $p$-subgroup of one complement is conjugate to a Sylow $p$-subgroup of the other. As a corollary, we find that any two supersoluble complements of an abelian subgroup $N$ in a finite split extension $G$ are conjugate if and only if, for each prime $p$, there exists a Sylow $p$-subgroup $S$ of $G$ such that any two complements of $S\cap N$ in $S$ are conjugate in $G$. In particular, restricting to supersoluble groups allows us to ease D. G. Higman's stipulation that the complements of $S\cap N$ in $S$ be conjugate within $S$. We then consider group actions and obtain a fixed point result for non-coprime actions analogous to Glauberman's lemma., Comment: to appear in the Proceedings of the Edinburgh Mathematical Society (2022)

Known as the "Missing Baryon Problem", about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations used a novel stacking analysis and detected an OVII absorption line toward the sightline of a luminous quasar, hinting that the missing baryons may reside in the WHIM. To explore how the properties of the OVII absorption line depend on feedback physics, we compare the observational results with predictions obtained from the Cosmology and Astrophysics with MachinE Learning (CAMEL) Simulation suite. CAMELS consists of cosmological simulations with state-of-the-art supernova (SN) and active galactic nuclei (AGN) feedback models from the IllustrisTNG and SIMBA simulations, with varying strengths. We find that the simulated OVII column densities are higher in the outskirts of galaxies than in the large-scale WHIM, but they are consistently lower than those obtained in the Chandra observations, for all feedback runs. We establish that the OVII distribution is primarily sensitive to changes in the SN feedback prescription, whereas changes in the AGN feedback prescription have minimal impact. We also find significant differences in the OVII column densities between the IllustrisTNG and SIMBA runs. We conclude that the tension between the observed and simulated OVII column densities cannot be explained by the wide range of feedback models implemented in CAMELS., Comment: 11 pages, 9 figures. Matching to the published version on MNRAS

Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a $50$ Myr timescale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in halos, without a need for jet precession or wide opening angles. Comparing the jet feedback with a `kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star forming phase in place, while the wind acts more locally on the star forming phase, which manifests itself in different time-variability properties., Comment: 22 pages, 20 figures, accepted for publication in MNRAS, comments welcome

The mass distribution of dense cores is a potential key to understand the process of star formation. Applying dendrogram analysis to the CARMA-NRO Orion C$^{18}$O ($J$=1--0) data, we identify 2342 dense cores, about 22 \% of which have virial ratios smaller than 2, and can be classified as gravitationally bound cores. The derived core mass function (CMF) for bound starless cores which are not associate with protostars has a slope similar to Salpeter's initial mass function (IMF) for the mass range above 1 $M_\odot$, with a peak at $\sim$ 0.1 $M_\odot$. We divide the cloud into four parts based on the declination, OMC-1/2/3, OMC-4/5, L1641N/V380 Ori, and L1641C, and derive the CMFs in these regions. We find that starless cores with masses greater than 10 $M_\odot$ exist only in OMC-1/2/3, whereas the CMFs in OMC-4/5, L1641N, and L1641C are truncated at around 5--10 $M_\odot$. From the number ratio of bound starless cores and Class II objects in each subregion, the lifetime of bound starless cores is estimated to be 5--30 free-fall times, consistent with previous studies for other regions. In addition, we discuss core growth by mass accretion from the surrounding cloud material to explain the coincidence of peak masses between IMFs and CMFs. The mass accretion rate required for doubling the core mass within a core lifetime is larger than that of Bondi-Hoyle accretion by a factor of order 2. This implies that more dynamical accretion processes are required to grow cores., Comment: 58 pages, 33 figures, 21 tables, accepted by ApJS

A supersonic relative velocity between dark matter (DM) and baryons (the stream velocity) at the time of recombination induces the formation of low mass objects with anomalous properties in the early Universe. We widen the scope of the `Supersonic Project' paper series to include objects we term Dark Matter + Gas Halos Offset by Streaming (DM GHOSts)--diffuse, DM-enriched structures formed because of a physical offset between the centers of mass of DM and baryonic overdensities. We present an updated numerical investigation of DM GHOSts and Supersonically Induced Gas Objects (SIGOs), including the effects of molecular cooling, in high resolution hydrodynamic simulations using the AREPO code. Supplemented by an analytical understanding of their ellipsoidal gravitational potentials, we study the population-level properties of these objects, characterizing their morphology, spin, radial mass, and velocity distributions in comparison to classical structures in non-streaming regions. The stream velocity causes deviations from sphericity in both the gas and DM components and lends greater rotational support to the gas. Low mass ($<\sim 10^{5.5}$ M$_\odot$) objects in regions of streaming demonstrate core-like rotation and mass profiles. Anomalies in the rotation and morphology of DM GHOSts could represent an early Universe analogue to observed ultra-faint dwarf galaxies with variations in DM content and unusual rotation curves., Comment: 27 pages, 20 figures. Accepted ApJ

The CRESST experiment employs cryogenic calorimeters for the sensitive measurement of nuclear recoils induced by dark matter particles. The recorded signals need to undergo a careful cleaning process to avoid wrongly reconstructed recoil energies caused by pile-up and read-out artefacts. We frame this process as a time series classification task and propose to automate it with neural networks. With a data set of over one million labeled records from 68 detectors, recorded between 2013 and 2019 by CRESST, we test the capability of four commonly used neural network architectures to learn the data cleaning task. Our best performing model achieves a balanced accuracy of 0.932 on our test set. We show on an exemplary detector that about half of the wrongly predicted events are in fact wrongly labeled events, and a large share of the remaining ones have a context-dependent ground truth. We furthermore evaluate the recall and selectivity of our classifiers with simulated data. The results confirm that the trained classifiers are well suited for the data cleaning task., Comment: 12 pages, 8 figures, 6 tables

Active galactic nuclei (AGN) feedback models are generally calibrated to reproduce galaxy observables such as the stellar mass function and the bimodality in galaxy colors. We use variations of the AGN feedback implementations in the IllustrisTNG (TNG) and Simba cosmological hydrodynamic simulations to show that the low redshift Lyman-$\alpha$ forest can provide constraints on the impact of AGN feedback. We show that TNG over-predicts the number density of absorbers at column densities $N_{\rm HI} < 10^{14}$ cm$^{-2}$ compared to data from the Cosmic Origins Spectrograph (in agreement with previous work), and we demonstrate explicitly that its kinetic feedback mode, which is primarily responsible for galaxy quenching, has a negligible impact on the column density distribution (CDD) of absorbers. In contrast, we show that the fiducial Simba model which includes AGN jet feedback is the preferred fit to the observed CDD of the $z = 0.1$ Lyman-$\alpha$ forest across five orders of magnitude in column density. We show that the Simba results with jets produce a quantitatively better fit to the observational data than the Simba results without jets, even when the UVB is left as a free parameter. AGN jets in Simba are high speed, collimated, weakly-interacting with the interstellar medium (via brief hydrodynamic decoupling) and heated to the halo virial temperature. Collectively these properties result in stronger long-range impacts on the IGM when compared to TNG's kinetic feedback mode, which drives isotropic winds with lower velocities at the galactic radius. Our results suggest that the low redshift Lyman-$\alpha$ forest provides plausible evidence for long-range AGN jet feedback., Comment: 15 pages, 5 figures, 1 table, accepted to ApJL