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

Abstract

The first measurements of the 21 cm brightness temperature power spectrum from the epoch of reionization will very likely be achieved in the near future by radio interferometric array experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the Square Kilometre Array (SKA). Standard MCMC analyses use an explicit likelihood approximation to infer the reionization parameters from the 21 cm power spectrum. In this paper, we present a new Bayesian inference of the reionization parameters where the likelihood is implicitly defined through forward simulations using density estimation likelihood-free inference (DELFI). Realistic effects including thermal noise and foreground avoidance are also applied to the mock observations from the HERA and SKA. We demonstrate that this method recovers accurate posterior distributions for the reionization parameters, and outperforms the standard MCMC analysis in terms of the location and size of credible parameter regions. With the minutes-level processing time once the network is trained, this technique is a promising approach for the scientific interpretation of future 21 cm power spectrum observation data. Our code 21cmDELFI-PS is publicly available at this link.

arXiv comment15 pages, 8 figures, 4 tables. Accepted for publication in ApJ. Comments welcome

References in corpus (49)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• Astropy: A Community Python Package for Astronomy
• Planck 2015 results. XIII. Cosmological parameters
• Binary companions of evolved stars in APOGEE DR14: Search method and catalog of ~5,000 companions
• LOFAR: The LOw-Frequency ARray
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• The Murchison Widefield Array: the Square Kilometre Array Precursor at low radio frequencies
• The frontier of simulation-based inference
• The GMRT 150 MHz All-sky Radio Survey: First Alternative Data Release TGSS ADR1
• 21cmFAST: A Fast, Semi-Numerical Simulation of the High-Redshift 21-cm Signal
• Hydrogen Epoch of Reionization Array (HERA)
• Efficient Simulations of Early Structure Formation and Reionization
• The Precision Array for Probing the Epoch of Reionization: 8 Station Results
• What Next-Generation 21 cm Power Spectrum Measurements Can Teach Us About the Epoch of Reionization
• Reionization and the Cosmic Dawn with the Square Kilometre Array
• A simulation calibrated limit on the HI power spectrum from the GMRT Epoch of Reionization experiment
• Improved upper limits on the 21-cm signal power spectrum of neutral hydrogen at $\mathbf{z}\approx 9.1$ from LOFAR
• New Limits on 21cm EoR From PAPER-32 Consistent with an X-Ray Heated IGM at z=7.7
• 21CMMC: an MCMC analysis tool enabling astrophysical parameter studies of the cosmic 21 cm signal
• Introducing the THESAN project: radiation-magnetohydrodynamic simulations of the epoch of reionization
• Deep multi-redshift limits on Epoch of Reionisation 21cm Power Spectra from Four Seasons of Murchison Widefield Array Observations
• First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum
• The Baryon Acoustic Oscillation Broadband and Broad-beam Array: Design Overview and Sensitivity Forecasts
• Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology
• Redshift Space Distortion of the 21cm Background from the Epoch of Reionization I: Methodology Re-examined
• Nuisance hardened data compression for fast likelihood-free inference
• Lyα emission line reconstruction for high-$z$ QSOs
• Emulating Simulations of Cosmic Dawn for 21cm Power Spectrum Constraints on Cosmology, Reionization, and X-ray Heating
• Approximate Bayesian Computation for Astronomical Model Analysis: A Case Study in Galaxy Demographics and Morphological Transformation at High Redshift
• Deep learning from 21-cm tomography of the Cosmic Dawn and Reionization
• 21CMMC with a 3D light-cone: the impact of the co-evolution approximation on the astrophysics of reionisation and cosmic dawn
• Analysing the 21cm signal from the Epoch of Reionization with artificial neural networks
• Emulation of reionization simulations for Bayesian inference of astrophysics parameters using neural networks
• PAPER-64 Constraints On Reionization II: The Temperature Of The z=8.4 Intergalactic Medium
• Constraining the astrophysics and cosmology from 21cm tomography using deep learning with the SKA
• Simulation-Based Inference of Reionization Parameters From 3D Tomographic 21 cm Lightcone Images
• Statistics of the epoch of reionization (EoR) 21-cm signal -- II. The evolution of the power spectrum error-covariance
• Evaluating machine learning techniques for predicting power spectra from reionization simulations
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• A machine learning approach to galaxy properties: joint redshift-stellar mass probability distributions with Random Forest
• Machine learning astrophysics from 21 cm lightcones: impact of network architectures and signal contamination
• Adding Context to JWST Surveys with Current and Future 21cm Radio Observations
• Epoch of reionization parameter estimation with the 21-cm bispectrum
• The impact of non-Gaussianity on the error covariance for observations of the Epoch of Reionization 21-cm power spectrum
• Improved supervised learning methods for EoR parameters reconstruction
• The impact of non-Gaussianity on the Epoch of Reionization parameter forecast using 21-cm power spectrum measurements
• Extracting the 21-cm Power Spectrum and the reionization parameters from mock datasets using Artificial Neural Networks
• Understanding the Impact of Semi-numeric Reionization Models when using CNNs

Cited by in corpus (24)

• Exploring the likelihood of the 21-cm power spectrum with simulation-based inference
• 21cmEMU: an emulator of 21cmFAST summary observables
• Constraining the X-ray heating and reionization using 21-cm power spectra with Marginal Neural Ratio Estimation
• Efficient parameter inference for gravitational wave signals in the presence of transient noises using temporal and time-spectral fusion normalizing flow
• Calibrating cosmological simulations with implicit likelihood inference using galaxy growth observables
• How informative are summaries of the cosmic 21-cm signal?
• CoLFI: Cosmological Likelihood-free Inference with Neural Density Estimators
• Scalable inference with Autoregressive Neural Ratio Estimation
• Extracting the Global 21-cm signal from Cosmic Dawn and Epoch of Reionization in the presence of Foreground and Ionosphere
• Efficient exploration of reionization parameters for the upcoming 21 cm observations using a photon conserving semi-numerical model SCRIPT
• Efficient simulation of discrete galaxy populations and associated radiation fields during the first billion years
• Deep learning-driven likelihood-free parameter inference for 21-cm forest observations
• Analyzing the 21cm forest with Wavelet Scattering Transform: Insight into non-Gaussian features of the 21cm forest
• Comparing sampling techniques to chart parameter space of 21 cm Global signal with Artificial Neural Networks
• Comparison of Bayesian inference methods using the Loreli II database of hydro-radiative simulations of the 21-cm signal
• Reionization Parameter Inference from 3D Minkowski Functionals of the 21 cm Signals
• Cosmological Constraints with Void Lensing I: the Simulation-Based Inference Framework
• Reionization relics in the cross-correlation between the Ly$\alpha$ forest and 21 cm intensity mapping in the post-reionization era
• Reionisation time fields reconstruction from 21 cm signal maps
• Map Reconstruction of radio observations with Conditional Invertible Neural Networks
• From ANN to BNN: Inferring Reionization Parameters using Uncertainty-aware Emulators of 21-cm Summaries
• 3D ScatterNet: Inference from 21 cm Light-cones
• A bubble size distribution model for the Epoch of Reionization
• Nonlinear reconstruction of 21cm global signal from 21cm power spectrum with artificial neural networks