Comparison of Bayesian inference methods using the Loreli II database of hydro-radiative simulations of the 21-cm signal
arXiv:2411.03093 · doi:10.1051/0004-6361/202452901
Abstract
While the observation of the 21 cm signal from the Cosmic Dawn and Epoch of Reionization is an instrumental challenge, the interpretation of a prospective detection is still open to questions regarding the modelling of the signal and the Bayesian inference techniques that bridge the gap between theory and observations. To address some of these questions, we present Loreli II, a database of nearly 10 000 simulations of the 21 cm signal run with the Licorice 3D radiative transfer code. With Loreli II, we explore a 5-dimensional astrophysical parameter space where star formation, X-ray emissions, and UV emissions are varied. We then use this database to train neural networks and perform Bayesian inference on 21 cm power spectra affected by thermal noise at the level of 100 hours of observation with the Square Kilometer Array. We study and compare three inference techniques : an emulator of the power spectrum, a Neural Density Estimator that fits the implicit likelihood of the model, and a Bayesian Neural Network that directly fits the posterior distribution. We measure the performances of each method by comparing them on a statistically representative set of inferences, notably using the principles of Simulation-Based Calibration. We report errors on the 1-D marginalized posteriors (biases and over/under confidence) below of the standard deviation for the emulator and below for the other methods. We conclude that at our noise level and our sampling density of the parameter space, an explicit Gaussian likelihood is sufficient. This may not be the case at lower noise level or if a denser sampling is used to reach higher accuracy. We then apply the emulator method to recent HERA upper limits and report weak constraints on the X-ray emissivity parameter of our model.
21 pages, 16 figures, submitted to A&A. Processed simulated data (21-cm power spectra and global signals) available at https://21ssd.obspm.fr/
References in corpus (46)
- Planck 2015 results. XIII. Cosmological parameters
- The cosmological simulation code GADGET-2
- Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
- An absorption profile centred at 78 megahertz in the sky-averaged spectrum
- 21-cm cosmology
- Hands-on Bayesian Neural Networks -- a Tutorial for Deep Learning Users
- The Evolution of the Galaxy Rest-Frame Ultraviolet Luminosity Function Over the First Two Billion Years
- Cosmological Parameter Estimation Using 21 cm Radiation from the Epoch of Reionization
- New Determinations of the UV Luminosity Functions from z~9 to z~2 show a remarkable consistency with halo growth and a constant star formation efficiency
- The Dearth of z~10 Galaxies in all HST Legacy Fields -- The Rapid Evolution of the Galaxy Population in the First 500 Myr
- Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at from LOFAR
- Full-Data Results of Hubble Frontier Fields: UV Luminosity Functions at and a Consistent Picture of Cosmic Reionization
- ALMA Spectroscopic Survey in the Hubble Ultra Deep Field: The Infrared Excess of UV-selected z=2-10 galaxies as a function of UV-continuum Slope and Stellar Mass
- 21CMMC: an MCMC analysis tool enabling astrophysical parameter studies of the cosmic 21 cm signal
- The z = 9-10 galaxy population in the Hubble Frontier Fields and CLASH surveys: The z=9 LF and further evidence for a smooth decline in UV luminosity at z >= 8
- The Extreme Faint End of the UV Luminosity Function at Through Gravitational Telescopes: a comprehensive assessment of strong lensing uncertainties
- C^2-Ray: A new method for photon-conserving transport of ionizing radiation
- Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations
- Inhomogeneous recombinations during cosmic reionization
- HERA Phase I Limits on the Cosmic 21-cm Signal: Constraints on Astrophysics and Cosmology During the Epoch of Reionization
- Simulating IGM Reionization
- UV Luminosity Density Results at z>8 from the First JWST/NIRCam Fields: Limitations of Early Data Sets and the Need for Spectroscopy
- Nuisance hardened data compression for fast likelihood-free inference
- Light cone effect on the reionization 21-cm power spectrum
- Generalized massive optimal data compression
- 21CMMC with a 3D light-cone: the impact of the co-evolution approximation on the astrophysics of reionisation and cosmic dawn
- The Scattering of Lyman-series Photons in the Intergalactic Medium
- The subtlety of Ly-a photons: changing the expected range of the 21-cm signal
- Simulation-Based Inference of Reionization Parameters From 3D Tomographic 21 cm Lightcone Images
- A new MWA limit on the 21 cm Power Spectrum at Redshifts 13 17
- The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization
- First upper limits on the 21 cm signal power spectrum from cosmic dawn from one night of observations with NenuFAR
- Evaluating machine learning techniques for predicting power spectra from reionization simulations
- Prediction of the 21-cm signal from reionization: comparison between 3D and 1D radiative transfer schemes
- 21SSD: a public database of simulated 21-cm signals from the epoch of reionization
- Implicit Likelihood Inference of Reionization Parameters from the 21 cm Power Spectrum
- Exploring the likelihood of the 21-cm power spectrum with simulation-based inference
- 21cmEMU: an emulator of 21cmFAST summary observables
- 21cmVAE: A Very Accurate Emulator of the 21-cm Global Signal
- Constraining the X-ray heating and reionization using 21-cm power spectra with Marginal Neural Ratio Estimation
- BEoRN: A fast and flexible framework to simulate the epoch of reionisation and cosmic dawn
- How informative are summaries of the cosmic 21-cm signal?
- Accurate modelling of the Lyman- coupling for the 21-cm signal, observability with NenuFAR and SKA
- The collective X-ray luminosity of HMXB as a SFR indicator
- A prescription for the conditional mass function of dark matter haloes
- Understanding the Impact of Semi-numeric Reionization Models when using CNNs