Abstract:The detection of the 21-cm signal at $z\gtrsim6$ will reveal insights into the properties of the first galaxies responsible for driving reionisation. To extract this information, we perform parameter inference which requires embedding 3D simulations of the 21-cm signal within a Bayesian inference pipeline. Presently, when performing inference we must choose which sources of uncertainty to sample and which to hold fixed. Since the astrophysics of galaxies are much more uncertain than those of the underlying halo-mass function (HMF), we usually parameterise and model the former while fixing the latter. However, in doing so we may bias our inference of the properties of these first galaxies. In this work, we explore the consequences of assuming an incorrect choice of HMF and quantify the relative biases in our inferred astrophysical model parameters when considering the wrong HMF. We then relax this assumption by constructing a generalised five parameter model for the HMF and simultaneously recover these parameters along with our underlying astrophysical model. For this analysis, we use 21cmFAST and perform Simulation-Based Inference by applying marginal neural ratio estimation to learn the likelihood-to-evidence ratio using Swyft. Using a mock 1000 hour observation of the 21-cm power spectrum from the forthcoming Square Kilometre Array, conservatively assuming foreground wedge avoidance, we find assuming the incorrect HMF can bias the recovered astrophysical parameters by up to $\sim3-4\sigma$ even when including independent information from observed luminosity functions. When considering our generalised HMF model, we recover constraints on our astrophysical parameters with a factor of $\sim2-4$ larger marginalised uncertainties. Importantly, these constraints are unbiased, agnostic to the underlying HMF and therefore more conservative.

What problem does this paper attempt to address?

This paper explores the impact of the Halo Mass Function (HMF) on inferring astrophysical parameters during the Epoch of Reionisation (EoR). The authors point out that when analyzing the 21cm signal to understand the properties of early galaxies, the HMF is often fixed, while more uncertainties are attributed to the astrophysical properties of galaxies. However, this approach may lead to biased inferences of these early galaxy properties. The paper uses simulations of the 21cm signal and employs a Bayesian inference pipeline to investigate the biases introduced by incorrect HMF assumptions and quantifies the relative biases in the astrophysical model parameters when considering incorrect HMFs. They utilize the 21CMFAST simulation and the Swyft tool for the foundational inference, learning the likelihood-to-evidence ratio to avoid direct sampling in the high-dimensional parameter space. The study finds that incorrect HMF assumptions can introduce biases of up to 3-4 standard deviations in the recovered astrophysical parameters, even when incorporating independent observed luminosity function information. Subsequently, they employ a five-parameter generic HMF model to simultaneously recover the astrophysical parameters and HMF parameters, finding that the resulting constraints are unbiased and insensitive to the specific form of the HMF, although the uncertainties increase by 2-4 times, making the results more conservative. The paper also discusses the computational costs of using traditional methods such as 21CMMC for Bayesian inference and highlights the advantages of using machine learning methods like foundational inference, such as handling non-Gaussian statistics and improving efficiency. Finally, the paper investigates the impact of HMF on parameter inference using simulated 1000-hour observations of the 21cm power spectrum with the Square Kilometre Array (SKA) and demonstrates the differences in the global history and mean brightness temperature of the 21cm signal under different HMF models.

Exploring the role of the halo mass function for inferring astrophysical parameters during reionisation

Exploring the role of the halo-mass function for inferring astrophysical parameters during reionization

Inferring astrophysical parameters using the 2D cylindrical power spectrum from reionisation

Inferring astrophysical parameters using the 2D cylindrical power spectrum from reionization

Simultaneously constraining the astrophysics of reionisation and the epoch of heating with 21CMMC

Reionization Parameter Inference from 3D Minkowski Functionals of the 21 cm Signals

Fisher matrix forecasts on the astrophysics of galaxies during the epoch of reionisation from the 21-cm power spectra

Implicit Likelihood Inference of Reionization Parameters from the 21 Cm Power Spectrum

21cmfish: Fisher-matrix framework for fast parameter forecasts from the cosmic 21-cm signal

Predicting the 21 cm field with a Hybrid Effective Field Theory approach

Exploring the likelihood of the 21-cm power spectrum with simulation-based inference

How Accurately Can 21 Cm Tomography Constrain Cosmology?

Cosmic variance suppression in radiation-hydrodynamic modeling of the reionization-era 21-cm signal

Dark matter halo properties from spatially integrated HI flux profiles

Cosmic variance suppression in radiation-hydrodynamic modelling of the reionization-era 21-cm signal

Constraining cosmology and ionization history with combined 21 cm power spectrum and global signal measurements

Constraining the X-ray heating and reionization using 21-cm power spectra with Marginal Neural Ratio Estimation

21-cm fluctuations from primordial magnetic fields

Testing common approximations to predict the 21cm signal at the Epoch of Reionization and Cosmic Dawn

Impact of astrophysical scatter on the Epoch of Reionization [H I]$_{\rm 21cm}$ bispectrum