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%.

arXiv comment32 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 (65)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• A Universal Density Profile from Hierarchical Clustering
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sample
• Large scale bias and the peak background split
• Profiles of dark haloes: evolution, scatter, and environment
• Statistical Properties of X-ray Clusters: Analytic and Numerical Comparisons
• Halo Models of Large Scale Structure
• Ellipsoidal collapse and an improved model for the number and spatial distribution of dark matter haloes
• An analytic model for the spatial clustering of dark matter haloes
• The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological Implications from two Decades of Spectroscopic Surveys at the Apache Point observatory
• Analytic model for galaxy and dark matter clustering
• DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations
• Halo occupation numbers and galaxy bias
• The Large Scale Bias of Dark Matter Halos: Numerical Calibration and Model Tests
• Power Spectra for Cold Dark Matter and its Variants
• The Five-Hundred-Meter Aperture Spherical Radio Telescope (FAST) Project
• Probing Dark Energy with Baryonic Acoustic Oscillations from Future Large Galaxy Redshift Surveys
• Hydrogen Epoch of Reionization Array (HERA)
• An accurate halo model for fitting non-linear cosmological power spectra and baryonic feedback models
• 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
• Cosmology with Phase 1 of the Square Kilometre Array; Red Book 2018: Technical specifications and performance forecasts
• The halo mass function through the cosmic ages
• Accurate halo-model matter power spectra with dark energy, massive neutrinos and modified gravitational forces
• DESI 2024 IV: Baryon Acoustic Oscillations from the Lyman Alpha Forest
• Cosmic Complementarity: H_0 and Omega_m from Combining CMB Experiments and Redshift Surveys
• Data Analysis for Precision 21 cm Cosmology
• PINOCCHIO: pinpointing orbit-crossing collapsed hierarchical objects in a linear density field
• Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder
• All-Sky Interferometry with Spherical Harmonic Transit Telescopes
• Coaxing Cosmic 21cm Fluctuations from the Polarized Sky using m-mode Analysis
• HIRAX: A Probe of Dark Energy and Radio Transients
• Improving lognormal models for cosmological fields
• Neutral hydrogen surveys for high redshift galaxy clusters and proto-clusters
• An accurate tool for the fast generation of dark matter halo catalogs
• An observational and theoretical view of the radial distribution of HI gas in galaxies
• Strong Gravitational Lensing and Velocity Function as Tools to Probe Cosmological Parameters: Current Constraints and Future Predictions
• Cosmological Perturbation Theory Using the FFTLog: Formalism and Connection to QFT Loop Integrals
• HI intensity mapping with MeerKAT: power spectrum detection in cross-correlation with WiggleZ galaxies
• HI constraints from the cross-correlation of eBOSS galaxies and Green Bank Telescope intensity maps
• Predicting the number, spatial distribution and merging history of dark matter haloes
• A Virtual Sky with Extragalactic HI and CO Lines for the SKA and ALMA
• LRP 2020 Whitepaper: The Canadian Hydrogen Observatory and Radio-transient Detector (CHORD)
• 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
• Fluctuations in 21cm Emission After Reionization
• A halo model for cosmological neutral hydrogen: abundances and clustering
• Improving fast generation of halo catalogs with higher-order Lagrangian perturbation theory
• SKAO HI Intensity Mapping: Blind Foreground Subtraction Challenge
• Tomographic galaxy clustering with the Subaru Hyper Suprime-Cam first year public data release
• MeerKLASS: MeerKAT Large Area Synoptic Survey
• The atomic Hydrogen content of the post-reionization Universe
• PINOCCHIO and the hierarchical build-up of dark matter haloes
• The bias of DLAs at z ~ 2.3: contraining stellar feedback in shallow potential wells
• The BINGO telescope: a new instrument exploring the new 21-cm cosmology window
• Angpow: a software for the fast computation of accurate tomographic power spectra
• 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