Search

Recent JWST observations at $z > 6$ may imply galactic ionizing photon production in excess of prior expectations. Under observationally motivated assumptions about escape fractions, these suggest a $z \sim 8-9$ end to reionization, in strong tension with the $z < 6$ end required by the Ly$\alpha$ forest. In this work, we use radiative transfer simulations to understand what different observations tell us about when reionization ended and when it started. We consider a model that ends too early (at $z \approx 8$) alongside two more realistic scenarios that end late at $z \approx 5$: one that starts late ($z \sim 9$) and another that starts early ($z \sim 13$). We find that the latter requires up to an order-of-magnitude evolution in galaxy ionizing properties at $6 < z < 12$, perhaps in tension with recent measurements of $\xi_{\rm ion}$ by JWST, which indicate little evolution. We also study how these models compare to recent measurements of the Ly$\alpha$ forest opacity, mean free path, IGM thermal history, visibility of $z > 8$ Ly$\alpha$ emitters, and the patchy kSZ signal from the CMB. We find that neither of the late-ending scenarios is conclusively disfavored by any single data set. However, a majority of these observables, spanning several distinct types of observations, prefer a late start. Not all probes agree with this conclusion, hinting at a possible lack of concordance between observables. Observations by multiple experiments (including JWST, Roman, and CMB-S4) in the coming years will either establish a concordance picture of reionization's early stages or reveal systematics in data and/or theoretical modeling., Comment: 21+7 pages, 13+2 figures, submitted to ApJ. Comments welcome

The connection between galaxies and their host dark matter (DM) halos is critical to our understanding of cosmology, galaxy formation, and DM physics. To maximize the return of upcoming cosmological surveys, we need an accurate way to model this complex relationship. Many techniques have been developed to model this connection, from Halo Occupation Distribution (HOD) to empirical and semi-analytic models to hydrodynamic. Hydrodynamic simulations can incorporate more detailed astrophysical processes but are computationally expensive; HODs, on the other hand, are computationally cheap but have limited accuracy. In this work, we present NeHOD, a generative framework based on variational diffusion model and Transformer, for painting galaxies/subhalos on top of DM with an accuracy of hydrodynamic simulations but at a computational cost similar to HOD. By modeling galaxies/subhalos as point clouds, instead of binning or voxelization, we can resolve small spatial scales down to the resolution of the simulations. For each halo, NeHOD predicts the positions, velocities, masses, and concentrations of its central and satellite galaxies. We train NeHOD on the TNG-Warm DM suite of the DREAMS project, which consists of 1024 high-resolution zoom-in hydrodynamic simulations of Milky Way-mass halos with varying warm DM mass and astrophysical parameters. We show that our model captures the complex relationships between subhalo properties as a function of the simulation parameters, including the mass functions, stellar-halo mass relations, concentration-mass relations, and spatial clustering. Our method can be used for a large variety of downstream applications, from galaxy clustering to strong lensing studies., Comment: Submitted to ApJ; 30 + 6 pages; 11 + 4 figures; Comments welcomed

Analyses of the cosmic 21-cm signal are hampered by astrophysical foregrounds that are far stronger than the signal itself. These foregrounds, typically confined to a wedge-shaped region in Fourier space, often necessitate the removal of a vast majority of modes, thereby degrading the quality of the data anisotropically. To address this challenge, we introduce a novel deep generative model based on stochastic interpolants to reconstruct the 21-cm data lost to wedge filtering. Our method leverages the non-Gaussian nature of the 21-cm signal to effectively map wedge-filtered 3D lightcones to samples from the conditional distribution of wedge-recovered lightcones. We demonstrate how our method is able to restore spatial information effectively, considering both varying cosmological initial conditions and astrophysics. Furthermore, we discuss a number of future avenues where this approach could be applied in analyses of the 21-cm signal, potentially offering new opportunities to improve our understanding of the Universe during the epochs of cosmic dawn and reionization., Comment: 12 pages, 5 figures

Observations of the 21-cm signal are opening a window to the cosmic-dawn epoch, when the first stars formed. These observations are usually interpreted with semi-numerical or hydrodynamical simulations, which are often computationally intensive and inflexible to changes in cosmological or astrophysical effects. Here, we present an effective, fully analytic model for the impact of the first stars on the 21-cm signal, using the modular code Zeus21. Zeus21 employs an analytic prescription of the star formation rate density (SFRD) to recover the fully nonlinear and nonlocal correlations of radiative fields that determine the 21-cm signal. We introduce the earliest Population III (Pop III) stars residing in low-mass molecular-cooling galaxies in Zeus21, with distinct spectra from later Pop II stars. We also self-consistently model feedback in the form of $H_2$-dissociating Lyman-Werner (LW) radiation, as well as dark matter-baryon relative velocities, both of which suppress star formation in the lowest-mass halos. LW feedback produces a scale-dependence on the SFRD fluctuations, due to the long mean free path of LW photons. Relative velocities give rise to "wiggles" in the spatial distribution of the 21-cm signal; we present an improved calculation of the shape of these velocity-induced acoustic oscillations, showing they remain a standard ruler at cosmic dawn. Our improved version of Zeus21 predicts the 21-cm global signal and power spectra in agreement with simulations at the $\sim 10\%$ level, yet is at least three orders of magnitude faster. This public code represents a step towards efficient and flexible parameter inference at cosmic dawn, allowing us to predict the first billion years of the universe in mere seconds., Comment: 37 pages, 17 figures, submitted to PRD

We study the cosmological impact of warm, dark-sector relic particles produced as Hawking radiation in a primordial-black-hole-dominated universe before big bang nucleosynthesis. If those dark-sector particles are stable, they would survive to the present day as "Hawking relics" and modify the subsequent growth of cosmological structure. We show that such relics are produced with much larger momenta, but in smaller quantities than the familiar thermal relics considered in standard cosmology. Consequently, Hawking relics with keV-MeV masses affect the growth of large-scale structure in a similar way to eV-scale thermal relics like massive neutrinos. We model their production and evolution, and show that their momentum distributions are broader than comparable relics with thermal distributions. Warm Hawking relics affect the growth of cosmological perturbations and we constrain their abundance to be less than $2\%$ of the dark matter over a broad range of their viable parameter space. Finally, we examine how future measurements of the matter power spectrum can distinguish Hawking relics from thermal particles., Comment: 34 pages, 13 figures, prepared for submission to JCAP

We introduce the DREAMS project, an innovative approach to understanding the astrophysical implications of alternative dark matter models and their effects on galaxy formation and evolution. The DREAMS project will ultimately comprise thousands of cosmological hydrodynamic simulations that simultaneously vary over dark matter physics, astrophysics, and cosmology in modeling a range of systems -- from galaxy clusters to ultra-faint satellites. Such extensive simulation suites can provide adequate training sets for machine-learning-based analyses. This paper introduces two new cosmological hydrodynamical suites of Warm Dark Matter, each comprised of 1024 simulations generated using the Arepo code. One suite consists of uniform-box simulations covering a $(25~h^{-1}~{\rm M}_\odot)^3$ volume, while the other consists of Milky Way zoom-ins with sufficient resolution to capture the properties of classical satellites. For each simulation, the Warm Dark Matter particle mass is varied along with the initial density field and several parameters controlling the strength of baryonic feedback within the IllustrisTNG model. We provide two examples, separately utilizing emulators and Convolutional Neural Networks, to demonstrate how such simulation suites can be used to disentangle the effects of dark matter and baryonic physics on galactic properties. The DREAMS project can be extended further to include different dark matter models, galaxy formation physics, and astrophysical targets. In this way, it will provide an unparalleled opportunity to characterize uncertainties on predictions for small-scale observables, leading to robust predictions for testing the particle physics nature of dark matter on these scales., Comment: 28 pages, 8 figures, DREAMS website: https://www.dreams-project.org

The formation redshift and abundance of the first stars and galaxies is highly sensitive to the build up of low mass dark matter halos as well as astrophysical feedback effects which modulate star formation in these low mass halos. The 21-cm signal at cosmic dawn will depend strongly on the formation of these first luminous sources and thus can be used to constrain unknown astrophysical and dark matter properties in the early universe. In this paper, we explore how well we could measure properties of dark matter using the 21-cm power spectrum at $z>10$, given unconstrained astrophysical parameters. We create a generalizable form of the dark matter halo mass function for models with damped and/or oscillatory linear power spectra, finding a single "smooth-k" window function which describes a broad range of models including CDM. We use this to make forecasts for structure formation using the Effective Theory of Structure Formation (ETHOS) framework to explore a broad parameter space of dark matter models. We make predictions for the 21-cm power spectrum observed by HERA varying both cosmological ETHOS parameters as well as astrophysical parameters. Using a Markov Chain Monte Carlo forecast we find that the ETHOS dark matter parameters are degenerate with astrophysical parameters linked to star formation in low mass dark matter halos but not with X-ray heating produced by the first generation of stars. After marginalizing over uncertainties in astrophysical parameters we demonstrate that with just 540 days of HERA observations it should be possible to distinguish between CDM and a broad range of dark matter models with suppression at wavenumbers $k\lesssim 200\,h$Mpc$^{-1}$ assuming a moderate noise level. These results demonstrate the potential of 21-cm observations to constrain the matter power spectrum on scales smaller than current probes., Comment: 13 pages, 9 figures, submitted to PRD

New James Webb Space Telescope (JWST) observations are revealing the first galaxies to be prolific producers of ionizing photons, which we argue gives rise to a tension between different probes of reionization. Over the last two decades a consensus has emerged where star-forming galaxies are able to generate enough photons to drive reionization, given reasonable values for their number densities, ionizing efficiencies $\xi_{\rm ion}$ (per unit UV luminosity), and escape fractions $f_{\rm esc}$. However, some new JWST observations infer high values of $\xi_{\rm ion}$ during reionization and an enhanced abundance of earlier ($z\gtrsim 9$) galaxies, dramatically increasing the number of ionizing photons produced at high $z$. Simultaneously, recent low-$z$ studies predict significant escape fractions for faint reionization-era galaxies. Put together, we show that the galaxies we have directly observed ($M_{\rm UV} < -15$) not only can drive reionization, but would end it too early. That is, our current galaxy observations, taken at face value, imply an excess of ionizing photons and thus a process of reionization in tension with the cosmic microwave background (CMB) and Lyman-$\alpha$ forest. Considering galaxies down to $M_{\rm UV}\approx -11$, below current observational limits, only worsens this tension. We discuss possible avenues to resolve this photon budget crisis, including systematics in either theory or observations., Comment: 6+3 pages, 3+2 figures, 1 table, updated to match accepted version

The 21-cm signal provides a novel avenue to measure the thermal state of the universe during cosmic dawn and reionization (redshifts $z\sim 5-30$), and thus to probe energy injection from decaying or annihilating dark matter (DM). These DM processes are inherently inhomogeneous: both decay and annihilation are density dependent, and furthermore the fraction of injected energy that is deposited at each point depends on the gas ionization and density, leading to further anisotropies in absorption and propagation. In this work, we develop a new framework for modeling the impact of spatially inhomogeneous energy injection and deposition during cosmic dawn, accounting for ionization and baryon density dependence, as well as the attenuation of propagating photons. We showcase how this first completely inhomogeneous treatment affects the predicted 21-cm power spectrum in the presence of exotic sources of energy injection, and forecast the constraints that upcoming HERA measurements of the 21-cm power spectrum will set on DM decays to photons and to electron/positron pairs. These projected constraints considerably surpass those derived from CMB and Lyman-$\alpha$ measurements, and for decays to electron/positron pairs they exceed all existing constraints in the sub-GeV mass range, reaching lifetimes of $\sim 10^{28}\,\mathrm{s}$. Our analysis demonstrates the unprecedented sensitivity of 21-cm cosmology to exotic sources of energy injection during the cosmic dark ages. Our code, $\mathtt{DM21cm}$, includes all these effects and is publicly available in an accompanying release., Comment: 33 pages, 22 figures, public code at https://github.com/yitiansun/DM21cm

The first stars are expected to form through molecular-hydrogen (H$_2$) cooling, a channel that is especially sensitive to the thermal and ionization state of gas, and can thus act as a probe of exotic energy injection from decaying or annihilating dark matter (DM). Here, we use a toy halo model to study the impact of DM-sourced energy injection on the H$_2$ content of the first galaxies, and thus estimate the threshold mass required for a halo to form stars at high redshifts. We find that currently allowed DM models can significantly change this threshold, producing both positive and negative feedback. In some scenarios, the extra heating of the gas raises the halo mass required for collapse, whereas in others, energy injection lowers the threshold by increasing the free-electron fraction and catalyzing H$_2$ formation. The direction of the effect can be redshift-dependent. We also bracket the uncertainties from self-shielding of halos from Lyman-Werner radiation. Hence, exotic energy injection can both delay and accelerate the onset of star formation; we show how this can impact the timing of 21cm signals at cosmic dawn. We encourage detailed simulation follow-ups in the most promising regions of parameter space identified in this work., Comment: 32 pages, 13 figures, 3 tables

The high-redshift galaxy UV luminosity function (UVLF) has become essential for understanding the formation and evolution of the first galaxies. Yet, UVLFs only measure galaxy abundances, giving rise to a degeneracy between the mean galaxy luminosity and its stochasticity. Here, we show that upcoming clustering measurements with the James Webb Space Telescope (JWST), as well as with Roman, will be able to break this degeneracy, even at redshifts $z \gtrsim 10$. First, we demonstrate that current Subaru Hyper Suprime-Cam (HSC) measurements of the galaxy bias at $z\sim 4-6$ point to a relatively tight halo-galaxy connection, with low stochasticity. Then, we show that the larger UVLFs observed by JWST at $z\gtrsim 10$ can be explained with either a boosted average UV emission or an enhanced stochasticity. These two models, however, predict different galaxy biases, which are potentially distinguishable in JWST and Roman surveys. Galaxy-clustering measurements, therefore, will provide crucial insights into the connection between the first galaxies and their dark-matter halos, and identify the root cause of the enhanced abundance of $z \gtrsim 10$ galaxies revealed with JWST during its first year of operations., Comment: 6+3 pages, 3+3 figures. Published in MNRAS Lett

The early-science observations made by the James Webb Space Telescope (JWST) have revealed an excess of ultra-massive galaxy candidates that appear to challenge the standard cosmological model ($\Lambda$CDM). Here, we argue that any modifications to $\Lambda$CDM that can produce such ultra-massive galaxies in the early Universe would also affect the UV galaxy luminosity function (UV LF) inferred from the Hubble Space Telescope (HST). The UV LF covers the same redshifts ($z\approx 7-10$) and host-halo masses $(M_\mathrm{h}\approx 10^{10}-10^{12}\, M_\odot$) as the JWST candidates, but tracks star-formation rate rather than stellar mass. We consider beyond-$\Lambda$CDM power-spectrum enhancements and show that any departure large enough to reproduce the abundance of ultra-massive JWST candidates is in conflict with the HST data. Our analysis, therefore, severely disfavors a cosmological explanation for the JWST abundance problem. Looking ahead, we determine the maximum allowable stellar-mass function and provide projections for the high-$z$ UV LF given our constraints on cosmology from current HST data., Comment: 10 pages, 5 figures - v2: Matches version accepted for publication in PRL

We test the impact of gravitational lensing on the lifetime estimates of seven high-redshift quasars at redshift $z\gtrsim6$. The targeted quasars are identified by their small observed proximity zone sizes, which indicate extremely short quasar lifetimes $(t_Q\lesssim10^5 \text{ yrs})$. However, these estimates of quasar lifetimes rely on the assumption that the observed luminosities of the quasars are intrinsic and not magnified by gravitational lensing, which would bias the lifetime estimates towards younger ages. In order to test possible effects of gravitational lensing, we obtain high-resolution images of the seven quasars with the {\em Hubble Space Telescope (HST)} and look for signs of strong lensing. We do not find any evidence of strong lensing, i.e., all quasars are well-described by point sources, and no foreground lensing galaxy is detected. We estimate that the strong lensing probabilities for these quasars are extremely small $(\sim1.4\times10^{-5})$, and show that weak lensing changes the estimated quasar lifetimes by only $\lesssim0.2$ dex. We thus confirm that the short lifetimes of these quasars are intrinsic. The existence of young quasars indicates a high obscured fraction, radiatively inefficient accretion, and/or flickering light curves for high-redshift quasars. We further discuss the impact of lensing magnification on measurements of black hole masses and Eddington ratios of quasars., Comment: 14 pages, 4 figures, 3 tables. Accepted by the Astrophysical Journal

The 21-cm signal holds the key to understanding the first structure formation during cosmic dawn. Theoretical progress over the last decade has focused on simulations of this signal, given the nonlinear and nonlocal relation between initial conditions and observables (21-cm or reionization maps). Here, instead, we propose an {\it effective} and fully analytic model for the 21-cm signal during cosmic dawn. We take advantage of the exponential-like behavior of the local star-formation rate density (SFRD) against densities at early times to analytically find its correlation functions including nonlinearities. The SFRD acts as the building block to obtain the statistics of radiative fields (X-ray and Lyman-$\alpha$ fluxes), and therefore the 21-cm signal. We implement this model as the public Python package Zeus21. This code can fully predict the 21-cm global signal and power spectrum in $\sim 1$ s, with negligible memory requirements. When comparing against state-of-the-art semi-numerical simulations from 21CMFAST we find agreement to $\sim10\%$ precision in both the 21-cm global signal and power spectra, after accounting for a (previously missed) underestimation of adiabatic fluctuations in 21CMFAST. Zeus21 is modular, allowing the user to vary the astrophysical model for the first galaxies, and interfaces with the cosmological code CLASS, which enables searches for beyond standard-model cosmology in 21-cm data. This represents a step towards bringing 21-cm to the era of precision cosmology., Comment: 20+4 pages, 21+4 figures. Updated to match accepted version. Code available at https://github.com/JulianBMunoz/Zeus21

The 21-cm signal from neutral hydrogen in the early universe will provide unprecedented information about the first stars and galaxies. Extracting this information, however, requires accounting for many unknown astrophysical processes. Semi-numerical simulations are key for exploring the vast parameter space of said processes. These simulations use approximate techniques such as excursion-set and perturbation theory to model the 3D evolution of the intergalactic medium, at a fraction of the computational cost of hydrodynamic and/or radiative transfer simulations. However, exploring the enormous parameter space of the first galaxies can still be computationally expensive. Here we introduce 21cmfish, a Fisher-matrix wrapper for the semi-numerical simulation 21cmFAST. 21cmfish facilitates efficient parameter forecasts, scaling to significantly higher dimensionalities than MCMC approaches, assuming a multi-variate Gaussian posterior. Our method produces comparable parameter uncertainty forecasts to previous MCMC analyses but requires ~10$^4$x fewer simulations. This enables a rapid way to prototype analyses adding new physics and/or additional parameters. We carry out a forecast for HERA using the largest astrophysical parameter space to-date, with 10 free parameters, spanning both population II and III star formation. We find X-ray parameters for the first galaxies could be measured to sub-percent precision, and, though they are highly degenerate, the stellar-to-halo mass relation and ionizing photon escape fraction for population II and III galaxies can be constrained to ~10% precision (logarithmic quantities). Using a principal component analysis we find HERA is most sensitive to the product of the ionizing escape fraction and the stellar-to-halo mass fraction for population II galaxies., Comment: 18 pages, accepted for publication in MNRAS. Fisher-matrix wrapper for 21cmFAST available here: https://github.com/charlottenosam/21cmfish

We present a Lagrangian model of galaxy clustering bias in which we train a neural net using the local properties of the smoothed initial density field to predict the late-time mass-weighted halo field. By fitting the mass-weighted halo field in the AbacusSummit simulations at z=0.5, we find that including three coarsely spaced smoothing scales gives the best recovery of the halo power spectrum. Adding more smoothing scales may lead to 2-5% underestimation of the large-scale power and can cause the neural net to overfit. We find that the fitted halo-to-mass ratio can be well described by two directions in the original high-dimension feature space. Projecting the original features into these two principal components and re-training the neural net either reproduces the original training result, or outperforms it with a better match of the halo power spectrum. The elements of the principal components are unlikely to be assigned physical meanings, partly owing to the features being highly correlated between different smoothing scales. Our work illustrates a potential need to include multiple smoothing scales when studying galaxy bias, and this can be done easily with machine-learning methods that can take in high dimensional input feature space., Comment: submitted to JCAP

We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $\Delta^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $\Delta^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars., Comment: 57 pages, 37 figures. Updated to match the accepted ApJ version. Corresponding author: Joshua S. Dillon

Cosmological and astrophysical observations currently provide the only robust, positive evidence for dark matter. Cosmic probes of dark matter, which seek to determine the fundamental properties of dark matter through observations of the cosmos, have emerged as a promising means to reveal the nature of dark matter. This report summarizes the current status and future potential of cosmic probes to inform our understanding of the fundamental nature of dark matter in the coming decade., Comment: Report of the CF3 Topical Group for Snowmass 2021; 35 pages, 10 figures, many references. V3 updates Fig 3-2 and the author list

We use Swampland principles to theoretically disfavor regions of the parameter space of dark matter and other darkly charged particles that may exist. The Festina Lente bound, the analogue of the Weak-Gravity conjecture in de Sitter, places constraints on the mass and charge of dark particles, which here we show cover regions in parameter space that are currently allowed by observations. As a consequence, a broad set of new ultra-light particles are in the Swampland, independently of their cosmic abundance, showing the complementarity of Quantum Gravity limits with laboratory and astrophysical studies. In parallel, a Swampland bound on the UV cutoff associated to the axion giving a St\"{u}ckelberg photon its longitudinal mode translates to a new constraint on the kinetic mixings and masses of dark photons. This covers part of the parameter space targeted by upcoming dark-photon direct-detection experiments. Moreover, it puts astrophysically interesting models in the Swampland, including freeze-in dark matter through an ultra-light dark photon, as well as radio models invoked to explain the 21-cm EDGES anomaly., Comment: 31 pages, 5 figures

Primordial Black Holes (PBHs) are a viable candidate to comprise some or all of the dark matter and provide a unique window into the high-energy physics of the early universe. This white paper discusses the scientific motivation, current status, and future reach of observational searches for PBHs. Future observational facilities supported by DOE, NSF, and NASA will provide unprecedented sensitivity to PBHs. However, devoted analysis pipelines and theoretical modeling are required to fully leverage these novel data. The search for PBHs constitutes a low-cost, high-reward science case with significant impact on the high energy physics community., Comment: 22 pages, 4 figures, submission to the Snowmass 2021 process. Update to add references

Many well-motivated extensions of the Standard Model predict the existence of new light species that may have been produced in the early universe. Prominent examples include axions, sterile neutrinos, gravitinos, dark photons, and more. The gravitational influence of light relics leaves imprints in the cosmic microwave background fluctuations, the large-scale structure of the universe and the primordial element abundances. In this paper, we detail the physics of cosmological light relics, and describe how measurements of their relic density and mass serve as probes of physics beyond the Standard Model. A measurement of the light relic density at the precision of upcoming cosmological surveys will point the way toward new physics or severely constrain the range of viable extensions to the Standard Model., Comment: Contribution to Snowmass 2021. 59 pages, 12 figures

The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determine the underlying properties of dark matter across the vast landscape of dark matter theories. This white paper summarizes the ongoing rapid development of theoretical and experimental methods, as well as new opportunities, to use dark matter halo measurements as a pillar of dark matter physics., Comment: White paper submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). 88 pages, 9 figures. Comments welcome

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged), Comment: Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)

While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade., Comment: Contribution to Snowmass 2021, 24 pages, 4 figures, comments welcome

Upper limits from the current generation of interferometers targeting the 21-cm signal from high redshifts have recently begun to rule out physically realistic, though still extreme, models of the Epoch of Reionization (EoR). While inferring the detailed properties of the first galaxies is one of the most important motivations for measuring the high-$z$ 21-cm signal, they can also provide useful constraints on the properties of the intergalactic medium (IGM). Motivated by this, we build a simple, phenomenological model for the 21-cm power spectrum that works directly in terms of IGM properties, which bypasses the computationally expensive 3-D semi-numerical modeling generally employed in inference pipelines and avoids explicit assumptions about galaxy properties. The key simplifying assumptions are that (i) the ionization field is binary, and composed of spherical bubbles with an abundance described well by a parametric bubble size distribution, and (ii) that the spin temperature of the ``bulk'' IGM outside bubbles is uniform. Despite the simplicity of the model, the mean ionized fraction and spin temperature of the IGM recovered from mock 21-cm power spectra generated with \textsc{21cmfast} are generally in good agreement with the true input values. This suggests that it is possible to obtain comparable constraints on the IGM using models with very different assumptions, parameters, and priors. Our approach will thus be complementary to semi-numerical models as upper limits continue to improve in the coming years., Comment: 23 pages, 14 figures, accepted for publication in MNRAS, code publicly available at https://github.com/mirochaj/micro21cm

The formation of the first galaxies during cosmic dawn and reionization (at redshifts $z=5-30$), triggered the last major phase transition of our universe, as hydrogen evolved from cold and neutral to hot and ionized. The 21-cm line of neutral hydrogen will soon allow us to map these cosmic milestones and study the galaxies that drove them. To aid in interpreting these observations, we upgrade the publicly available code {\tt 21cmFAST}. We introduce a new, flexible parametrization of the additive feedback from: an inhomogeneous, $H_2$-dissociating (Lyman-Werner; LW) background; and dark matter -- baryon relative velocities; which recovers results from recent, small-scale hydrodynamical simulations with both effects. We perform a large, "best-guess" simulation as the 2021 installment of the Evolution of 21-cm Structure (EOS) project. This improves the previous release with a galaxy model that reproduces the observed UV luminosity functions (UVLFs), and by including a population of molecular-cooling galaxies. The resulting 21-cm global signal and power spectrum are significantly weaker, primarily due to a more rapid evolution of the star-formation rate density required to match the UVLFs. Nevertheless, we forecast high signal-to-noise detections for both HERA and the SKA. We demonstrate how the stellar-to-halo mass relation of the unseen, first galaxies can be inferred from the 21-cm evolution. Finally, we show that the spatial modulation of X-ray heating due to relative velocities provides a unique acoustic signature that is detectable at $z \approx 10-15$ in our fiducial model. Ours are the first public simulations with joint inhomogeneous LW and relative-velocity feedback across the entire cosmic dawn and reionization., Comment: 23+3 pages, 23+3 figures. Matches accepted version. Simulation output available at https://www.dropbox.com/sh/dqh9r6wb0s68jfo/AACc9ZCqsN0SQ_JJN7GRVuqDa?dl=0

Observations of high-redshift galaxies have provided us with a rich tool to study the physics at play during the epoch of reionisation. The luminosity function (LF) of these objects is an indirect tracer of the complex processes that govern galaxy formation, including those of the first dark-matter structures. In this work, we present an extensive analysis of the UV galaxy LF at high redshifts to extract cosmological and astrophysical parameters. We provide a number of phenomenological approaches in modelling the UV LF and take into account various sources of uncertainties and systematics in our analysis, including cosmic variance, dust extinction, scattering in the halo-galaxy connection, and the Alcock-Paczy\'{n}ski effect. Using UV LF measurements from the Hubble Space Telescope together with external data on the matter density, we derive the large-scale matter clustering amplitude to be $\sigma_8=0.76^{+0.12}_{-0.14}$, after marginalising over the unknown astrophysical parameters. We find that with current data this result is only weakly sensitive to our choice of astrophysical modelling, as well as the calibration of the underlying halo mass function. As a cross check, we run our analysis pipeline with mock data from the IllustrisTNG hydrodynamical simulations and find consistent results with their input cosmology. In addition, we perform a simple forecast for future space telescopes, where an improvement of roughly 30% upon our current result is expected. Finally, we obtain constraints on astrophysical parameters and the halo-galaxy connection for the models considered here. All methods discussed in this work are implemented in the form of a versatile likelihood code, GALLUMI, which we make public., Comment: 25 pages, 12 figures - v3: Updated funding information, matches version accepted for publication in PRD

The epochs of cosmic dawn and reionisation present promising avenues for understanding the role of dark matter (DM) in our cosmos. The first galaxies that populated the Universe during these eras resided in DM halos that were much less massive than their counterparts today. Consequently, observations of such galaxies can provide us with a handle on the clustering of DM in an otherwise currently inaccessible regime. In this work, we use high-redshift UV galaxy luminosity-function (UV LF) data from the Hubble Space Telescope to study the clustering properties of DM at small scales. In particular, we present new measurements of the matter power spectrum at wavenumbers $0.5\,\mathrm{Mpc}^{-1} < k < 10\,\mathrm{Mpc}^{-1}$ to roughly 30% precision, obtained after marginalising over the unknown astrophysics. These new data points cover the uncharted redshift range $4\leq z\leq 10$ and encompass scales beyond those probed by Cosmic-Microwave-Background and large-scale-structure observations. This work establishes the UV LF as a powerful tool to probe the nature of DM in a different regime than other cosmological and astrophysical data sets., Comment: 8 pages, 4 figures - v2: Expanded on discussion of results, matches version accepted for publication in ApJL

We present a non-parametric Lagrangian biasing model and fit the ratio of the halo and mass densities at the field level using the mass-weighted halo field in the AbacusSummit simulations at $z=0.5$. Unlike the perturbative halo bias model widely used in interpreting the observed large-scale structure traced by galaxies, we find a non-negative halo-to-mass ratio that increases monotonically with the linear overdensity $\delta_1$ in the initial Lagrangian space. The bias expansion, however, does not guarantee non-negativity of the halo counts, and may give rising halo number counts at negative overdensities. The shape of the halo-to-mass ratio is unlikely to be described by a polynomial of $\delta_1$ and other quantities. Especially for massive halos with $6\times10^{12}\ h^{-1}\ M_\odot$, the halo-to-mass ratio starts soaring up at $\delta_1>0$, substantially different from the predictions of the bias expansion. We show that for the halo masses we consider ($M>3\times10^{11}\ h^{-1}\ M_\odot$) a non-parametric halo-to-mass ratio as a function of $\delta_1$ and its local derivative $\nabla^2\delta_1$ can recover the halo power spectra to sub-percent level for wavenumbers $k=0.01-0.1\ h\ {\rm Mpc}^{-1}$ given a proper smoothing scale to filter the initial density field, even though we do not fit the power spectrum directly. However, there is mild dependence of the recovery of the halo power spectrum on the smoothing scale and other input parameters. At $k<0.01\ h\ {\rm Mpc}^{-1}$ and for massive halos with $M>6\times10^{12}\ h^{-1}\ M_\odot$, our non-parametric model leads to a few percent overestimation of the halo power spectrum, indicating the need for larger or multiple smoothing scales. The halo-to-mass ratios obtained qualitatively agree with intuitions from extended Press-Schechter theory. We compare our framework to the bias expansion and discuss possible extensions., Comment: JCAP accepted

Despite the success of the standard $\Lambda$CDM model of cosmology, recent data improvements have made tensions emerge between low- and high-redshift observables, most importantly in determinations of the Hubble constant $H_0$ and the (rescaled) clustering amplitude $S_8$. The high-redshift data, from the cosmic microwave background (CMB), crucially relies on recombination physics for its interpretation. Here we study how small-scale baryon inhomogeneities (i.e., clumping) can affect recombination and consider whether they can relieve both the $H_0$ and $S_8$ tensions. Such small-scale clumping, which may be caused by primordial magnetic fields or baryon isocurvature below kpc scales, enhances the recombination rate even when averaged over larger scales, shifting recombination to earlier times. We introduce a flexible clumping model, parametrized via three spatial zones with free densities and volume fractions, and use it to study the impact of clumping on CMB observables. We find that increasing $H_0$ decreases both $\Omega_m$ and $S_8$, which alleviates the $S_8$ tension. On the other hand, the shift in $\Omega_m$ is disfavored by the low-$z$ baryon-acoustic-oscillations measurements. We find that the clumping parameters that can change the CMB sound horizon enough to explain the $H_0$ tension also alter the damping tail, so they are disfavored by current Planck 2018 data. We test how the CMB damping-tail information rules out changes to recombination by first removing $\ell>1000$ multipoles in Planck data, where we find that clumping could resolve the $H_0$ tension. Furthermore, we make predictions for future CMB experiments, as their improved damping-tail precision can better constrain departures from standard recombination. Both the Simons Observatory and CMB-S4 will provide decisive evidence for or against clumping as a resolution to the $H_0$ tension., Comment: 21 page, 20 figures. Code available: CLASS fork at https://github.com/misharash/class_public and mock CMB likelihood for Cobaya at https://github.com/misharash/cobaya_mock_cmb. Updated to match the published version

Many scenarios of physics beyond the standard model predict new light, weakly coupled degrees of freedom, populated in the early universe and remaining as cosmic relics today. Due to their high abundances, these relics can significantly affect the evolution of the universe. For instance, massless relics produce a shift $\Delta N_{\rm eff}$ to the cosmic expectation of the effective number of active neutrinos. Massive relics, on the other hand, additionally become part of the cosmological dark matter in the later universe, though their light nature allows them to freely stream out of potential wells. This produces novel signatures in the large-scale structure (LSS) of the universe, suppressing matter fluctuations at small scales. We present the first general search for such light (but massive) relics (LiMRs) with cosmic microwave background (CMB) and LSS data, scanning the 2D parameter space of their masses $m_X$ and temperatures $T_X^{(0)}$ today. In the conservative minimum-temperature ($T_X^{(0)}=0.91$ K) scenario, we rule out Weyl (and higher-spin) fermions -- such as the gravitino -- with $m_X\geq 2.26$ eV at 95% C.L., and set analogous limits of $m_X\leq 11.2, 1.06, 1.56$ eV for scalar, vector, and Dirac-fermion relics. This is the first search for LiMRs with joint CMB, weak-lensing, and full-shape galaxy data; we demonstrate that weak-lensing data is critical for breaking parameter degeneracies, while full-shape information presents a significant boost in constraining power relative to analyses with only baryon acoustic oscillation parameters. Under the combined strength of these datasets, our constraints are the tightest and most comprehensive to date., Comment: 13 pages, 7 figures, updated to match version appearing in PRD

Upcoming measurements of the 21-cm line of neutral hydrogen will open a new observational window into the early stages of structure growth, providing a unique opportunity for probing large-scale cosmological signatures using the small-scale signals from the first stars. In this paper we evaluate the detection significance of compensated isocurvature perturbations (CIPs) from observations of the 21-cm hydrogen-line during the cosmic-dawn era. CIPs are modulations of the relative baryon and dark-matter density that leave the total matter density unchanged. We find that, under different assumptions for feedback and foregrounds, the ongoing HERA and upcoming SKA1-low experiments will provide constraints on uncorrelated CIPs at the level of $\sigma(A_{\rm CIP})= 10^{-3}-10^{-4}$, comparable to the sensitivity of upcoming CMB experiments, and potentially exceeding the constraints from cosmic-variance limited BAO surveys., Comment: 13 pages, 5 figures, comments welcome

All the news that's (un)fit to publish., Comment: v2: Fixed typo in crossword solutions

Primordial black holes (PBHs), formed out of large overdensities in the early Universe, are a viable dark matter (DM) candidate over a broad range of masses. Ultra-light, asteroid-mass PBHs with masses around $10^{17}$ g are particularly interesting as current observations allow them to constitute the entire DM density. PBHs in this mass range emit $\sim$ MeV photons via Hawking radiation which can directly be detected by the gamma ray telescopes, such as the upcoming AMEGO. In this work we forecast how well an instrument with the sensitivity of AMEGO will be able to detect, or rule out, PBHs as a DM candidate, by searching for their evaporating signature when marginalizing over the Galactic and extra-Galactic gamma-ray backgrounds. We find that an instrument with the sensitivity of AMEGO could exclude non-rotating PBHs as the only DM component for masses up to $7 \times 10^{17}$ g at 95% confidence level (C.L.) for a monochromatic mass distribution, improving upon current bounds by nearly an order of magnitude. The forecasted constraints are more stringent for PBHs that have rotation, or which follow extended mass distributions., Comment: v2: 10 pages, 5 figures. Results Unchanged. Matches version published in Physical Review D

Upcoming data of the 21-cm hydrogen line during cosmic dawn ($z\sim 10-30$) will revolutionize our understanding of the astrophysics of the first galaxies. Here we present a case study on how to exploit those same measurements to learn about the nature of dark matter (DM) at small scales. Focusing on the Effective Theory of Structure Formation (ETHOS) paradigm, we run a suite of simulations covering a broad range of DM microphysics, connecting the output of $N$-body simulations to dedicated 21-cm simulations to predict the evolution of the 21-cm signal across the entire cosmic dawn. We find that observatories targeting both the global signal and the 21-cm power spectrum are sensitive to all ETHOS models we study, and can distinguish them from CDM if the suppression wavenumber is smaller than $k\approx 300\, h/$Mpc, even when accounting for feedback with a phenomenological model. This is an order of magnitude smaller comoving scales than currently constrained by other data sets, including the Lyman-$\alpha$ forest. Moreover, if a prospective 21-cm detection confirmed a deficiency of power at small scales, we show that ETHOS models with strong dark acoustic oscillations can be discriminated from the pure suppression of warm dark matter, showing the power of 21-cm data to understand the behavior of DM at the smallest physical scales., Comment: 16 pages, 16 figures. Comments are welcome

This brief code paper presents a new Python-wrapped version of the popular 21cm cosmology simulator, 21cmFAST. The new version, v3+, maintains the same core functionality of previous versions of 21cmFAST, but features a simple and intuitive interface, and a great deal more flexibility. This evolution represents the work of a formalized collaboration, and the new version, available publicly on GitHub, provides a single point-of-reference for all future upgrades and community-added features. In this paper, we describe simple usage of 21cmFAST, some of its new features, and provide a simple performance benchmark., Comment: 6 pages, 3 figures, published in JOSS

UV luminosity functions provide a wealth of information on the physics of galaxy formation in the early Universe. Given that this probe indirectly tracks the evolution of the mass function of dark matter halos, it has the potential to constrain alternative theories of structure formation. One of such scenarios is the existence of primordial non-Gaussianity at scales beyond those probed by observations of the Cosmic Microwave Background. Through its impact on the halo mass function, such small-scale non-Gaussianity would alter the abundance of galaxies at high redshifts. In this work we present an application of UV luminosity functions as measured by the Hubble Space Telescope to constrain the non-Gaussianity parameter $f_\mathrm{NL}$ for wavenumbers above a cut-off scale $k_{\rm cut}$. After marginalizing over the unknown astrophysical parameters and accounting for potential systematic errors, we arrive at a $2\sigma$ bound of $f_{\rm NL}=71^{+426}_{-237}$ for a cut-off scale $k_{\rm cut}=0.1\,\mathrm{Mpc}^{-1}$ in the bispectrum of the primordial gravitational potential. Moreover, we perform forecasts for the James Webb Space Telescope and the Nancy Grace Roman Space Telescope, finding an expected improvement of a factor $3-4$ upon the current bound., Comment: 18 pages, 13 figures, 2 tables - v2: journal version

The advanced LIGO O3a run catalog has been recently published, and it includes several events with unexpected mass properties, including mergers with individual masses in the lower and upper mass gaps, as well as mergers with unusually small mass ratios between the binary components. Here we entertain the possibility that these outliers are the outcome of hierarchical mergers of black holes or neutron stars in the dense environments of globular clusters. We use the coagulation equation to study the evolution of the black hole mass function within a typical cluster. Our prescription allows us to monitor how various global quantities change with time, such as the total mass and number of compact objects in the cluster, its overall merger rate, and the probability to form intermediate-mass black holes via a runaway process. By accounting for the LIGO observational bias, we predict the merger event distributions with respect to various variables such as the individual masses M1 and M2, their ratio q, and redshift z, and we compare our predictions with the published O3a data. We study how these distributions depend on the merger-rate and ejections parameters and produce forecasts for the (tight) constraints that can be placed on our model parameters using the future dataset of the O5 run. Finally, we also consider the presence of a static channel with no coagulation producing merger events alongside the dynamic channel, finding that the two can be distinguished based solely on the merger mass distribution with future O5 data., Comment: 20 pages, 16 figures, 1 table; Included recoil and delay-time effects and added a comparison with LIGO O3a results; Version accepted by MNRAS

A promising avenue to measure the total, and potentially individual, mass of neutrinos consists of leveraging cosmological datasets, such as the cosmic microwave background and surveys of the large-scale structure of the universe. In order to obtain unbiased estimates of the neutrino mass, however, many effects ought to be included. Here we forecast, via a Markov Chain Monte Carlo likelihood analysis, whether measurements by two galaxy surveys: DESI and {\it Euclid}, when added to the CMB-S4 experiment, are sensitive to two effects that can alter neutrino-mass measurements. The first is the slight difference in the suppression of matter fluctuations that each neutrino-mass hierarchy generates, at fixed total mass. The second is the growth-induced scale-dependent bias (GISDB) of haloes produced by massive neutrinos. We find that near-future surveys can distinguish hierarchies with the same total mass only at the $1\sigma$ level; thus, while these are poised to deliver a measurement of the sum of neutrino masses, they cannot significantly discern the mass of each individual neutrino in the foreseeable future. We further find that neglecting the GISDB induces up to a $1\sigma$ overestimation of the total neutrino mass, and we show how to absorb this effect via a redshift-dependent parametrization of the scale-independent bias., Comment: 12 pages, 5 figures. Matches PRD version

Cosmological data provide a powerful tool in the search for physics beyond the Standard Model (SM). An interesting target are light relics, new degrees of freedom which decoupled from the SM while relativistic. Nearly massless relics contribute to the radiation energy budget, and are commonly parametrized as variations in the effective number $N_{\rm eff}$ of neutrino species. Additionally, relics with masses greater than $10^{-4}$ eV become non-relativistic before today, and thus behave as matter instead of radiation. This leaves an imprint in the clustering of the large-scale structure of the universe, as light relics have important streaming motions, mirroring the case of massive neutrinos. Here we forecast how well current and upcoming cosmological surveys can probe light massive relics (LiMRs). We consider minimal extensions to the SM by both fermionic and bosonic relic degrees of freedom. By combining current and upcoming cosmic-microwave-background and large-scale-structure surveys, we forecast the significance at which each LiMR, with different masses and temperatures, can be detected. We find that a very large coverage of parameter space will be attainable by upcoming experiments, opening the possibility of exploring uncharted territory for new physics beyond the SM., Comment: 19 pages, 13 figures. v2: matches the published version

Cosmological measurements of the 21-cm line of neutral hydrogen are poised to dramatically enhance our understanding of the early universe. In particular, both the epochs of reionization and cosmic dawn remain largely uncharted, and the 21-cm signal is one of the few probes to reach them. The simplest 21-cm measurement is the global signal (GS), which corresponds to the averaged absorption or emission of 21-cm photons across the entire sky. While bright radio foregrounds swamp the cosmic signal over the entire frequency range observable, presenting a formidable hurdle, they can in principle be subtracted, given enough sensitivity. Here, however, we point out an additional---and irreducible---source of uncertainty for the 21-cm GS: cosmic variance. The cosmic-variance noise arises from the finite volume of the universe accessible to 21-cm experiments. Due to the cosmological redshifting of 21-cm photons, each observed frequency probes our universe during a particular cosmic age, corresponding to a narrow redshift slice. The presence of large 21-cm fluctuations makes the GS within each slice different than the GS averaged over the entire universe. We estimate the size of this cosmic-variance noise, and find that for a standard scenario it has a size of $\sim 0.1$ mK, which is $\sim 10\%$ of the size of the expected instrumental noise of a year-long experiment. Interestingly, cosmic variance can overtake instrumental noise for scenarios with extreme 21-cm fluctuations, such as those suggested to explain the sharpness of the claimed EDGES detection. Moreover, as large-scale 21-cm fluctuations are coherent over long distances, cosmic variance correlates the measurements of the GS at nearby redshifts, leading to off-diagonal uncertainties that have so far been neglected., Comment: 5+9 pages, 3+5 figures

The International Gamma-Ray Astrophysics Laboratory (INTEGRAL) satellite has yielded unprecedented measurements of the soft gamma-ray spectrum of our Galaxy. Here we use those measurements to set constraints on dark matter (DM) that decays or annihilates into photons with energies $E\approx 0.02-2$ MeV. First, we revisit the constraints on particle DM that decays or annihilates to photon pairs. In particular, for decaying DM, we find that previous limits were overstated by roughly an order of magnitude. Our new, conservative analysis finds that the DM lifetime must satisfy $\tau\gtrsim 5\times 10^{26}\,{\rm s}\times (m_{\chi}/\rm MeV)^{-1}$ for DM masses $m_{\chi}=0.054-3.6$ MeV. For MeV-scale DM that annihilates into photons INTEGRAL sets the strongest constraints to date. Second, we target ultralight primordial black holes (PBHs) through their Hawking radiation. This makes them appear as decaying DM with a photon spectrum peaking at $E\approx 5.77/(8\pi G M_{\rm PBH})$, for a PBH of mass $M_{\rm PBH}$. We use the INTEGRAL data to demonstrate that, at 95\% C.L., PBHs with masses less than $1.2\times 10^{17}$ g cannot comprise all of the DM, setting the tightest bound to date on ultralight PBHs., Comment: 9 pages, 5 figures. Updated to match accepted version

The very first galaxies that started the cosmic dawn likely resided in so-called "minihaloes", with masses of $\sim10^5$-$10^8\mathrm{M}_\odot$, accreting their gas from the intergalactic medium through H$_2$ cooling. Such molecularly cooled galaxies (MCGs) mostly formed in pristine environments, hosted massive, metal-free stars, and were eventually sterilized by the build-up of a disassociating (Lyman-Werner; LW) background. Therefore, their properties might be very different from the galaxies we see in the later Universe. Although MCGs are probably too faint to be observed directly, we could nevertheless infer their properties from the imprint they leave in the cosmic 21-cm signal. Here we quantify this imprint by extending the public simulation code 21cmFAST to allow for a distinct population of MCGs. We allow MCGs to have different properties from other galaxies, including unique scaling relations for their stellar-to-halo mass ratios, ionizing escape fractions, and spectral energy distributions. We track inhomogeneous recombinations, disassociative LW feedback, and photoheating from reionization. After demonstrating how MCGs can shape the 21-cm signal, we explore to what extent current observations can already place constraints on their properties. The cosmic microwave background optical depth from Planck sets an upper limit on the product of the ionizing escape fraction and the stellar mass in MCGs. When including also the timing of the putative EDGES absorption signal, we find an additional strong degeneracy between the stellar mass and the X-ray luminosity of MCGs. If proven to be of cosmic origin, the timing of the EDGES signal would have been set by MCGs., Comment: 19 pages, 8 figures, 1 table; Updated to match the published version. Only minor changes

The distribution of matter fluctuations in our universe is key for understanding the nature of dark matter and the physics of the early cosmos. Different observables have been able to map this distribution at large scales, corresponding to wavenumbers $k\lesssim 10$ Mpc$^{-1}$, but smaller scales remain much less constrained. The 21-cm line is a promising tracer of early stellar formation, which took place in small haloes (with masses $M\sim 10^6-10^8M_\odot$), formed out of matter overdensities with wavenumbers as large as $k\approx100$ Mpc$^{-1}$. Here we forecast how well both the 21-cm global signal, and its fluctuations, could probe the matter power spectrum during cosmic dawn ($z=12-25$). We find that the long-wavelength modes (with $k\lesssim40$ Mpc$^{-1}$) are highly degenerate with astrophysical parameters, whereas the modes with $k= (40-80)$ Mpc$^{-1}$ are more readily observable. This is further illustrated in terms of the principal components of the matter power spectrum, which peak at $k\sim 50$ Mpc$^{-1}$ both for a typical experiment measuring the 21-cm global signal and its fluctuations. We find that, imposing broad priors on astrophysical parameters, a global-signal experiment can measure the amplitude of the matter power spectrum integrated over $k= (40-80)$ Mpc$^{-1}$ with a precision of tens of percent. A fluctuation experiment, on the other hand, can constrain the power spectrum to a similar accuracy over both the $k=(40-60)$ Mpc$^{-1}$ and $(60-80)$ Mpc$^{-1}$ ranges even without astrophysical priors. The constraints outlined in this work would be able to test the behavior of dark matter at the smallest scales yet measured, for instance probing warm-dark matter masses up to $m_{\rm WDM}=8$ keV for the global signal and $14$ keV for the 21-cm fluctuations. This could shed light on the nature of dark matter beyond the reach of other cosmic probes., Comment: 23 pages, 16 figures. Matches accepted version

Fast radio bursts (FRBs) are highly energetic radio pulses from cosmological origins. Despite an abundance of detections, their nature remains elusive. At least a subset of FRBs is expected to repeat, as the daily FRB rate surpasses that of any known cataclysmic event, which has been confirmed by observations. One of the proposed mechanisms to generate repeating FRBs is supergiant pulses from young and highly spinning NSs, in which case FRBs could inherit the periodicity of their parent NS. Here we examine the consequences of such a population of periodic fast radio bursts (PFRBs). We calculate the rate and lifetime of PFRB progenitors, and find that each newly born highly spinning NS has to emit a number $N_{\rm PFRB}\sim 10^2$ of bursts during its active lifetime of $\tau\sim 100$ years, after which it becomes too dim and crosses a PFRB "death line" analogous to the pulsar one. We propose several tests of this hypothesis. First, the period of PFRBs would increase over time, and their luminosity would decrease, due to the NS spin-down. Second, PFRBs may show modest amounts of rotation measure, given the lack of expelled matter from the pulsar, as opposed to the magnetar-sourced FRBs proposed to explain the first repeater FRB 121102. As an example, we study whether the second confirmed repeater (FRB 180814) is a PFRB, given the preference for an inter-pulse separation of 13 ms within its sub-bursts. We show that, if confirmed, this period would place FRB 180814 in a different category as FRB 121102. We develop tests that would identify---and characterize---the prospective population of PFRBs., Comment: 8 pages, 1 figure. Comments are welcome

We devise a test of nonlinear departures from general relativity (GR) using time delays in strong gravitational lenses. We use a phenomenological model of gravitational screening as a step discontinuity in the measure of curvature per unit mass, at a radius $\Lambda$. The resulting slip between two scalar gravitational potentials leads to a shift in the apparent positions and time delays of lensed sources, relative to the GR predictions, of size $\gamma_{\rm PN}-1$. As a proof of principle, we use measurements of two lenses, RXJ1121-1231 and B1608+656, to constrain deviations from GR to be below $|\gamma_{\rm PN}-1| \leq 0.2 \times (\Lambda/100\, \rm kpc)$. These constraints are complementary to other current probes, and are the tightest in the range $\Lambda=10-200$ kpc, showing that future measurements of strong-lensing time delays have great promise to seek departures from general relativity on kpc-Mpc scales., Comment: 9 pages, 3 figures. Expanded discussion. Published in Phys Rev D

The redshifted 21-cm line of hydrogen holds great potential for the study of cosmology, as it can probe otherwise unobservable cosmic epochs. In particular, measurements of the 21-cm power spectrum during cosmic dawn---the era when the first stars were formed---will provide us with a wealth of information about the astrophysics of stellar formation, as well as the origin of fluctuations in our Universe. In addition to their usually considered density fluctuations, dark matter and baryons possess large relative velocities with respect to each other, due to the baryon acoustic oscillations (BAOs) suffered by the latter, which suppress the formation of stars during cosmic dawn, leaving an imprint on 21-cm observables during this era. Here we present 21cmvFAST, a version of the publicly available code 21cmFAST modified to account for this effect. Previous work has shown that the inclusion of relative velocities produces an acoustic modulation on the large-scale 21-cm power spectrum during cosmic dawn. By comparing analytic calculations with simulations from 21cmvFAST, here we demonstrate that this modulation takes the form of robust velocity-induced acoustic oscillations (VAOs), during both the Lyman-$\alpha$ coupling era and the subsequent epoch of heating. The unique shape of these VAOs, which is determined by the acoustic physics that generated the relative velocities, can be analytically computed and is easily distinguishable from the usual density-sourced fluctuations. We find that, under a broad range of astrophysical assumptions, the VAOs are detectable at high significance by the upcoming HERA interferometer, which could therefore confirm the presence of acoustic oscillations at cosmic dawn., Comment: 25 pages, 22 figures. Code available at https://github.com/JulianBMunoz/21cmvFAST. Updated to match accepted version

The matter in our Universe comes in two flavors: dark and baryonic. Of these, only the latter couples to photons, giving rise to the well-known baryon acoustic oscillations and, in the process, generating supersonic relative velocities between dark matter and baryons. These velocities---imprinted with the acoustic scale in their genesis---impede the formation of the first stars during cosmic dawn ($z\sim 20$), modulating the expected 21-cm signal from this era. In a companion paper we showed, combining numerical simulations and analytic models, that this modulation takes the form of robust velocity-induced acoustic oscillations (VAOs), with a well-understood shape that is frozen at recombination, and unaffected by the unknown astrophysics of star formation. Here we propose using these VAOs as a standard ruler at cosmic dawn. We find that three years of 21-cm power-spectrum data from the upcoming HERA interferometer should be able to measure the Hubble expansion rate $H(z)$ at $z=15-20$ to percent-level precision, ranging from $0.3\%$ to $11\%$ depending on the strength of astrophysical feedback processes and foregrounds. This would provide a new handle on the expansion rate of our Universe during an otherwise unprobed epoch, opening a window to the mysterious cosmic-dawn era., Comment: 7 pages, 3 figures. Updated to match accepted version

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal)., Comment: Astro2020 Decadal Survey science white paper; list of endorsers available at https://github.com/veragluscevic/Astro2020-DM-Cosmology-Endorsers