21-cm Cosmology Research Articles

Soft photon heating: a semi-analytic framework and applications to 21-cm cosmology

ABSTRACT The presence of an abundant population of low-frequency photons at high redshifts (such as a radio background) can source leading order effects on the evolution of the matter and spin temperatures through rapid free–free absorptions. This effect, known as soft photon heating, can have a dramatic impact on the differential brightness temperature, $\Delta T_{\rm b}$, a central observable in 21-cm cosmology. Here, we introduce a semi-analytic framework to describe the dynamics of soft photon heating, providing a simplified set of evolution equations and a useful numerical scheme which can be used to study this generic effect. We also perform quasi-instantaneous and continuous soft photon injections to elucidate the different regimes in which soft photon heating is expected to impart a significant contribution to the global 21-cm signal and its fluctuations. We find that soft photon backgrounds produced after recombination with spectral index $\gamma \gt 3.0$ undergo significant free–free absorption, and therefore this heating effect cannot be neglected. The effect becomes stronger with steeper spectral index, and in some cases the injection of a synchrotron-like spectrum ($\gamma = 3.6$) can suppress the amplitude of $\Delta T_{\rm b}$ relative to the standard model prediction (where an additional radio background is absent), making the global 21-cm signal even more difficult to detect in these scenarios.

Understanding spectral artefacts in SKA-Low 21-cm cosmology experiments: the impact of cable reflections

Understanding spectral artefacts in SKA-Low 21-cm cosmology experiments: the impact of cable reflections

Minimal mass of thermal dark matter and the viability of millicharged particles affecting 21-cm cosmology

Thermal freeze-out offers an attractive explanation of the dark matter density free from fine-tuning of initial conditions. For dark matter with a mass below tens of MeV, photons, electrons, and neutrinos are the only available direct Standard Model annihilation products. Using a full three-sector abundance calculation, we determine the minimal mass of dark matter, allowing for an arbitrary branching into electrons/photons and neutrinos that is compatible with current cosmological observations. The analysis takes into account the heat transfer between the various sectors from annihilation elastic scattering, representing the first fully self-consistent analysis that tracks the respective sectors’ temperatures. We thereby provide accurate thermal annihilation cross sections, particularly for velocity-dependent cases, and deduce the sensitivity of current and upcoming CMB experiments to MeV thermal dark matter. In the latter context, we also establish the fine-tuned parameter region where a tiny admixture of neutrinos in the final states rules in MeV-scale p-wave annihilating DM into electrons. Finally, we show that a sub-% millicharged dark matter with an interaction strength that interferes with 21-cm cosmology is still allowed when freeze-out is supplemented with annihilation into neutrinos. For all cases considered, we provide concrete particle physics models and supplement our findings with a discussion of other relevant experimental results. Published by the American Physical Society 2024

Probing new physics at cosmic dawn with 21-cm cosmology

21-cm cosmology provides an exciting opportunity to probe new physics dynamics in the early universe. In particular, a tiny sub-component of dark matter that interacts strongly with the visible sector may cool the gas in the intergalactic medium and significantly alter the expected absorption signal at Cosmic Dawn. However, the information about new physics in this observable is obscured by astrophysical systematic uncertainties. In the absence of a microscopic framework describing the astrophysical sources, these uncertainties can be encoded in a bottom up effective theory for the 21-cm observables in terms of unconstrained astrophysical fluxes. In this paper, we take a first step towards a careful assessment of the degeneracies between new physics effects and the uncertainties in these fluxes. We show that the latter can be constrained by combining measurements of the UV luminosity function, the Planck measurement of the CMB optical depth to reionization, and an upper bound on the unresolved X-ray flux. Leveraging those constraints, we demonstrate how new physics signatures can be disentangled from astrophysical effects. Focusing on the case of millicharged dark matter, we find sharp predictions, with small uncertainties within the viable parameter space.1

On the constraints on superconducting cosmic strings from 21-cm cosmology

ABSTRACT Constraints on the potential properties of superconducting cosmic strings provide an indirect probe of physics beyond the standard model at energies inaccessible to terrestrial particle colliders. In this study, we perform the first joint Bayesian analysis to extract constraints on superconducting cosmic strings from current 21-cm signal measurements while accounting rigorously for the uncertainties in foregrounds and high redshift astrophysics. We include the latest publicly available 21-cm power spectrum upper limits from HERA, 21-cm global signal data from SARAS 3, and the synergistic probe of the unresolved X-ray background in our final analysis. This paper thus constitutes the first attempt to use 21-cm power spectrum data to probe cosmic strings. In contrast to previous works, we find no strong constraints can be placed on superconducting cosmic strings from current 21-cm measurements. This is because of uncertainties in the X-ray emission efficiency of the first galaxies, with X-ray emissivities greater than 3 × 1040 erg s−1 M$_{\odot }^{-1}$ yr able to mask the presence of cosmic strings in the 21-cm signal. We conclude by discussing the prospects for future constraints from definitive 21-cm signal measurements and argue that the recently proposed soft photon heating should be cause for optimism due to its potential to break degeneracies that would have otherwise made the signatures of cosmic strings difficult to distinguish from those of astrophysical origin.

New tool for 21-cm cosmology. II. Investigating the effect of early linear fluctuations

New tool for 21-cm cosmology. II. Investigating the effect of early linear fluctuations

A Bayesian method to mitigate the effects of unmodelled time-varying systematics for 21-cm cosmology experiments

ABSTRACT Radio observations of the neutral hydrogen signal from the Cosmic Dawn and Epoch of Reionization have helped to provide constraints on the properties of the first stars and galaxies. Since this global 21-cm cosmological signal from the Cosmic Dawn is effectively constant on observing time-scales and since effects resulting from systematics will vary with time, the effects of these systematics can be mitigated without the need for a model of the systematic. We present a method to account for unmodelled time-varying systematics in 21-cm radio cosmology experiments using a squared exponential Gaussian process kernel to account for correlations between time bins in a fully Bayesian way. We find by varying the model parameters of a simulated systematic that the Gaussian process method improves our ability to recover the signal parameters by widening the posterior in the presence of a systematic and reducing the bias in the mean fit parameters. When varying the amplitude of a model sinusoidal systematic between 0.25 and 2.00 times the 21-cm signal amplitude and the period between 0.5 and 4.0 times the signal width, we find on average a 5 per cent improvement in the root mean squared error of the fitted signal. We can use the fitted Gaussian process hyperparameters to identify the presence of a systematic in the data, demonstrating the method’s utility as a diagnostic tool. Furthermore, we can use Gaussian process regression to calculate a mean fit to the residuals over time, providing a basis for producing a model of the time-varying systematic.

Marginal post-processing of Bayesian inference products with normalizing flows and kernel density estimators

ABSTRACT Bayesian analysis has become an indispensable tool across many different cosmological fields, including the study of gravitational waves, the cosmic microwave background, and the 21-cm signal from the Cosmic Dawn, among other phenomena. The method provides a way to fit complex models to data describing key cosmological and astrophysical signals and a whole host of contaminating signals and instrumental effects modelled with ‘nuisance parameters’. In this paper, we summarize a method that uses masked autoregressive flows and kernel density estimators to learn marginal posterior densities corresponding to core science parameters. We find that the marginal or ‘nuisance-free’ posteriors and the associated likelihoods have an abundance of applications, including the calculation of previously intractable marginal Kullback–Leibler divergences and marginal Bayesian model dimensionalities, likelihood emulation, and prior emulation. We demonstrate each application using toy examples, examples from the field of 21-cm cosmology, and samples from the Dark Energy Survey. We discuss how marginal summary statistics like the Kullback–Leibler divergences and Bayesian model dimensionalities can be used to examine the constraining power of different experiments and how we can perform efficient joint analysis by taking advantage of marginal prior and likelihood emulators. We package our multipurpose code up in the pip-installable code margarine for use in the wider scientific community.

Constraints on dark matter–neutrino interaction from 21-cm cosmology and forecasts on SKA1-Low

ABSTRACT In this article, we have done a thorough investigation of the possible effects of interaction between dark matter (DM) and neutrinos on reionization history. We have constrained the interaction strength using 21-cm cosmology and found out possible deviations from standard, non-interacting Lambda cold dark matter scenario. Comparing the results with the existing constraints from present cosmological observations reveals that 21-cm observations are more competent to constrain the interaction strength by a few orders of magnitude. We have also searched for prospects of detecting any such interaction in the upcoming 21-cm mission SKA1-Low by doing a forecast analysis and error estimation.

Interpreting the Hi 21-cm cosmology maps through Largest Cluster Statistics. Part I. Impact of the synthetic SKA1-Low observations

We analyse the evolution of the largest ionized region using the topological and morphological evolution of the redshifted 21-cm signal coming from the neutral hydrogen distribution during the different stages of reionization. For this analysis, we use the “Largest Cluster Statistics” — LCS. We mainly study the impact of the array synthesized beam on the LCS analysis of the 21-cm signal considering the upcoming low-frequency Square Kilometer Array (SKA1-Low) observations using a realistic simulation for such observation based on the 21cmE2E-pipeline using OSKAR. We find that bias in LCS estimation is introduced in synthetic observations due to the array beam. This in turn shifts the apparent percolation transition point towards the later stages of reionization. The biased estimates of LCS, occurring due to the effect of the lower resolution (lack of longer baselines) and the telescope synthesized beam will lead to a biased interpretation of the reionization history. This is important to note while interpreting any future 21-cm signal images from upcoming or future telescopes like the SKA, HERA, etc. We conclude that one may need denser uv-coverage at longer baselines for a better deconvolution of the array synthesized beam from the 21-cm images and a relatively unbiased estimate of LCS from such images.

Use of time dependent data in Bayesian global 21-cm foreground and signal modelling

ABSTRACT Global 21-cm cosmology aims to investigate the cosmic dawn and epoch of reionization by measuring the sky averaged H i absorption signal, which requires, accurate modelling of, or correction for, the bright radio foregrounds and distortions arising from chromaticity of the antenna beam. We investigate the effect of improving foreground modelling by fitting data sets from many observation times simultaneously in a single Bayesian analysis, fitting for the same parameter set by performing these fits on simulated data. We find that for a hexagonal dipole antenna, this simultaneous fitting produces a significant improvement in the accuracy of the recovered 21-cm signal, relative to fitting a time average of the data. Furthermore, the recovered models of the foreground are also seen to become more accurate by up to a factor of ∼2–3 relative to time averaged fitting. For a less chromatic log spiral antenna, no significant improvement in signal recovery was found by this process. However, the modelling of the foregrounds was still significantly improved. We also investigate extending this technique to fit multiple data sets from different antennas simultaneously for the same parameters. This is also found to improve both 21-cm signal and foreground modelling, to a higher degree than fitting data set from multiple times from the same antenna.

Effective bias expansion for 21-cm cosmology in redshift space

A near-future detection of the 21-cm signal from the epoch of reionization will provide unique opportunities to probe the underlying cosmology, provided that such cosmological information can be extracted with precision. To this end, we further develop effective field theory (EFT) inspired techniques for the 21-cm brightness temperature field during the epoch of reionization, incorporating renormalized bias and a treatment of redshift space distortions. Notably, we confirm that in redshift space, measures of the 21-cm brightness, e.g., the power spectrum, should have irreducible contributions that lack a bias coefficient and therefore contain direct, astrophysics-free information about the cosmological density field; in this work, we study this effect beyond linear order. To validate our theoretical treatment, we fit the predicted EFT Fourier-space shapes to the thesan suite of hydrodynamical simulations of reionization at the field level, where the considerable number of modes prevents overfitting. We find agreement at the level of a few percent between the 21-cm power spectrum from the EFT fits and simulations over the wave number range $k\ensuremath{\lesssim}0.8\text{ }\text{ }\mathrm{h}/\mathrm{Mpc}$ and neutral fraction ${x}_{\mathrm{HI}}\ensuremath{\gtrsim}0.4$, which is imminently measurable by the Hydrogen Epoch of Reionization Array and future experiments. The ability of the EFT to describe the 21-cm signal extends to simulations that have different astrophysical prescriptions for reionization as well as simulations with interacting dark matter.

Removing systematics-induced 21-cm foreground residuals by cross-correlating filtered data

Observations of the redshifted 21-cm signal emitted by neutral hydrogen represent a promising probe of large-scale structure in the universe. However, cosmological 21-cm signal is challenging to observe due to astrophysical foregrounds which are several orders of magnitude brighter. Traditional linear foreground removal methods can optimally remove foregrounds for a known telescope response but are sensitive to telescope systematic errors such as antenna gain and delay errors, leaving foreground contamination in the recovered signal. Non-linear methods such as principal component analysis, on the other hand, have been used successfully for foreground removal, but they lead to signal loss that is difficult to characterize and requires careful analysis. In this paper, we present a systematics-robust foreground removal technique which combines both linear and non-linear methods. We first obtain signal and foreground estimates using a linear filter. Under the assumption that the signal estimate is contaminated by foreground residuals induced by parameterizable systematic effects, we infer the systematics-induced contamination by cross-correlating the initial signal and foreground estimates. Correcting for the inferred error, we are able to subtract foreground contamination from the linearly filtered signal up to the first order in the amplitude of the telescope systematics. In simulations of an interferometric 21-cm survey, our algorithm removes foreground leakage induced by complex gain errors by one to two orders of magnitude in the power spectrum. Our technique thus eases the requirements on telescope characterization for modern and next-generation 21-cm cosmology experiments.

Bayesian data analysis for sky-averaged 21-cm experiments in the presence of ionospheric effects

ABSTRACT The ionosphere introduces chromatic distortions on low frequency radio waves, and thus poses a hurdle for 21-cm cosmology. In this paper, we introduce time-varying chromatic ionospheric effects on simulated antenna temperature data of a global 21-cm data analysis pipeline, and try to detect the injected global signal. We demonstrate that given turbulent ionospheric conditions, more than 5 per cent error in our knowledge of the ionospheric parameters could lead to comparatively low evidence and high root-mean-square error (RMSE), suggesting a false or null detection. When using a constant antenna beam for cases that include data at different times, the significance of the detection lowers as the number of time samples increases. It is also shown that for observations that include data at different times, readjusting beam configurations according to the time-varying ionospheric conditions should greatly improve the significance of a detection, yielding higher evidences and lower RMSE, and that it is a necessary procedure for a successful detection when the ionospheric conditions are not ideal.

Bayesian evidence-driven diagnosis of instrumental systematics for sky-averaged 21-cm cosmology experiments

AbstractWe demonstrate the effectiveness of a Bayesian evidence -based analysis for diagnosing and disentangling the sky-averaged 21-cm signal from instrumental systematic effects. As a case study, we consider a simulated REACH pipeline with an injected systematic. We demonstrate that very poor performance or erroneous signal recovery is achieved if the systematic remains unmodelled. These effects include sky-averaged 21-cm posterior estimates resembling a very deep or wide signal. However, when including parameterised models of the systematic, the signal recovery is dramatically improved in performance. Most importantly, a Bayesian evidence-based model comparison is capable of determining whether or not such a systematic model is needed as the true underlying generative model of an experimental dataset is in principle unknown. We, therefore, advocate a pipeline capable of testing a variety of potential systematic errors with the Bayesian evidence acting as the mechanism for detecting their presence.

Bounds on dark matter annihilation cross-sections from inert doublet model in the context of 21-cm cosmology of dark ages

We study the fluctuations in the brightness temperature of 21-cm signal [Formula: see text] at the dark ages ([Formula: see text]) with a dark matter (DM) candidate in Inert Doublet Model (IDM). We then explore the effects of different fractions of IDM DM on [Formula: see text] signal. The IDM DM masses are chosen in few tens of GeV region as well as in the high mass region beyond 500 GeV. It has been observed that the [Formula: see text] signal is more sensitive in the DM mass range of 70–80 GeV. A lower bound on annihilation cross-section for this DM is also obtained by analyzing the [Formula: see text] signal. This is found to lie within the range [Formula: see text] cm3/s for the IDM DM mass range 10 GeV[Formula: see text] GeV.

Peering into the dark (ages) with low-frequency space interferometers

The Dark Ages and Cosmic Dawn are largely unexplored windows on the infant Universe (z ~ 200–10). Observations of the redshifted 21-cm line of neutral hydrogen can provide valuable new insight into fundamental physics and astrophysics during these eras that no other probe can provide, and drives the design of many future ground-based instruments such as the Square Kilometre Array (SKA) and the Hydrogen Epoch of Reionization Array (HERA). We review progress in the field of high-redshift 21-cm Cosmology, in particular focussing on what questions can be addressed by probing the Dark Ages at z > 30. We conclude that only a space- or lunar-based radio telescope, shielded from the Earth’s radio-frequency interference (RFI) signals and its ionosphere, enable the 21-cm signal from the Dark Ages to be detected. We suggest a generic mission design concept, CoDEX, that will enable this in the coming decades.

Maxsmooth: rapid maximally smooth function fitting with applications in Global 21-cm cosmology

ABSTRACTMaximally Smooth Functions (MSFs) are a form of constrained functions in which there are no inflection points or zero crossings in high-order derivatives. Consequently, they have applications to signal recovery in experiments where signals of interest are expected to be non-smooth features masked by larger smooth signals or foregrounds. They can also act as a powerful tool for diagnosing the presence of systematics. The constrained nature of MSFs makes fitting these functions a non-trivial task. We introduce maxsmooth, an open-source package that uses quadratic programming to rapidly fit MSFs. We demonstrate the efficiency and reliability of maxsmooth by comparison to commonly used fitting routines and show that we can reduce the fitting time by approximately two orders of magnitude. We introduce and implement with maxsmooth Partially Smooth Functions, which are useful for describing elements of non-smooth structure in foregrounds. This work has been motivated by the problem of foreground modelling in 21-cm cosmology. We discuss applications of maxsmooth to 21-cm cosmology and highlight this with examples using data from the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) and the Large-aperture Experiment to Detect the Dark Ages (LEDA) experiments. We demonstrate the presence of a sinusoidal systematic in the EDGES data with a log-evidence difference of 86.19 ± 0.12 when compared to a pure foreground fit. MSFs are applied to data from LEDA for the first time in this paper and we identify the presence of sinusoidal systematics. maxsmooth is pip installable and available for download at https://github.com/htjb/maxsmooth.

Redundant-baseline calibration of the hydrogen epoch of reionization array

ABSTRACT In 21-cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21-cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions – both in data and in simulations – and present strategies for mitigating that structure.

Emulating the global 21-cm signal from Cosmic Dawn and Reionization

ABSTRACT The 21-cm signal of neutral hydrogen is a sensitive probe of the Epoch of Reionization (EoR), Cosmic Dawn, and the Dark Ages. Currently, operating radio telescopes have ushered in a data-driven era of 21-cm cosmology, providing the first constraints on the astrophysical properties of sources that drive this signal. However, extracting astrophysical information from the data is highly non-trivial and requires the rapid generation of theoretical templates over a wide range of astrophysical parameters. To this end emulators are often employed, with previous efforts focused on predicting the power spectrum. In this work, we introduce 21cmgem– the first emulator of the global 21-cm signal from Cosmic Dawn and the EoR. The smoothness of the output signal is guaranteed by design. We train neural networks to predict the cosmological signal using a database of ∼30 000 simulated signals which were created by varying seven astrophysical parameters: the star formation efficiency and the minimal mass of star-forming haloes; the efficiency of the first X-ray sources and their spectrum parametrized by spectral index and the low-energy cut-off; the mean-free path of ionizing photons, and the cosmic microwave background optical depth. We test the performance with a set of ∼2000 simulated signals, showing that the relative error in the prediction has an rms of 0.0159. The algorithm is efficient, with a running time per parameter set of 0.16 s. Finally, we use the database of models to check the robustness of relations between the features of the global signal and the astrophysical parameters that we previously reported.