Fast Simulation of Cosmological Neutral Hydrogen based on the Halo Model

Abstract

Cosmological neutral hydrogen (HI) surveys provide a promising tomographic probe of the post-reionization era and of the standard model of cosmology. Simulations of this signal are crucial for maximizing the utility of these surveys. We present a fast method for simulating the cosmological distribution of HI based on a halo model approach. Employing the approximate $\text{PINOCCHIO}$ code, we generate the past light cone of dark matter halos. Subsequently, the halos are populated with HI according to a HI-halo mass relation. The nature of 21 cm intensity mapping demands large-volume simulations with a high halo mass resolution. To fulfill both requirements, we simulate a past light cone for declinations between -15° and -35° in the frequency range from 700 to 800 MHz, matching HIRAX, the Hydrogen Intensity and Real-time Analysis eXperiment. We run $\text{PINOCCHIO}$ for a 1 h${}^{-3}$Gpc${}^3$ box with 6700${}^3$ simulation particles. With this configuration, halos with masses as low as M${}_{\text{min}}$ = 4.3 $\times$ 10${}^9$M${}_{\odot }$ are simulated, resulting in the recovery of more than 97% of the expected HI density. From the dark matter and HI past light cone, maps with a width of 5 MHz are created. To validate the simulations, we have implemented and present here an analytical dark matter and HI halo model in $\text{PyCosmo}$, a Python package tailored for theoretical cosmological predictions. We perform extensive comparisons between analytical predictions and the simulations for the mass function, mass density, power spectrum, and angular power spectrum for dark matter and HI. We find close agreement in the mass function and mass densities, with discrepancies within a few percent. For the three-dimensional power spectra and angular power spectra, we observe an agreement better than 10%.

32 pages, 13 figures. Updated to match published version in JCAP. The PyCosmo package is available on https://cosmology.ethz.ch/research/software-lab/PyCosmo.html and can be installed from the Python Package Index at https://pypi.org. The simulation dataset is available at https://cosmology.ethz.ch/research/software-lab/cosmological-neutral-hydrogen-simulation.html

References in corpus (20)

• Array Programming with NumPy
• DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations
• The Five-Hundred-Meter Aperture Spherical Radio Telescope (FAST) Project
• Baryon Acoustic Oscillation Intensity Mapping as a Test of Dark Energy
• Measurement of 21 cm brightness fluctuations at z ~ 0.8 in cross-correlation
• Late-time cosmology with 21cm intensity mapping experiments
• DESI 2024 IV: Baryon Acoustic Oscillations from the Lyman Alpha Forest
• Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder
• HIRAX: A Probe of Dark Energy and Radio Transients
• HI intensity mapping with MeerKAT: power spectrum detection in cross-correlation with WiggleZ galaxies
• Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment
• Forecasts on the Dark Energy and Primordial Non-Gaussianity Observations with the Tianlai Cylinder Array
• SKAO HI Intensity Mapping: Blind Foreground Subtraction Challenge
• MeerKLASS: MeerKAT Large Area Synoptic Survey
• CoLoRe: fast cosmological realisations over large volumes with multiple tracers
• The Hubble Sphere Hydrogen Survey
• Simulating cosmologies beyond $\Lambda$CDM with PINOCCHIO
• Symbolic Implementation of Extensions of the $\text{PyCosmo}$ Boltzmann Solver
• The HI intensity mapping power spectrum: insights from recent measurements
• A cross-bispectrum estimator for CMB-HI intensity mapping correlations

Cited by in corpus (3)

• Large, fast and accurate HI intensity maps with latent overlap diffusion
• SHAM-OT: Rapid Subhalo Abundance Matching with Optimal Transport
• Analytical modeling of the one-dimensional power spectrum of 21-cm forest based on a halo model method