21CMMC: an MCMC analysis tool enabling astrophysical parameter studies of the cosmic 21 cm signal

Abstract

We introduce 21CMMC: a parallelized, Monte Carlo Markov Chain analysis tool, incorporating the epoch of reionization (EoR) seminumerical simulation 21CMFAST. 21CMMC estimates astrophysical parameter constraints from 21 cm EoR experiments, accommodating a variety of EoR models, as well as priors on model parameters and the reionization history. To illustrate its utility, we consider two different EoR scenarios, one with a single population of galaxies (with a mass-independent ionizing efficiency) and a second, more general model with two different, feedback-regulated populations (each with mass-dependent ionizing efficiencies). As an example, combining three observations (z=8, 9 and 10) of the 21 cm power spectrum with a conservative noise estimate and uniform model priors, we find that interferometers with specifications like the Low Frequency Array/Hydrogen Epoch of Reionization Array (HERA)/Square Kilometre Array 1 (SKA1) can constrain common reionization parameters: the ionizing efficiency (or similarly the escape fraction), the mean free path of ionizing photons and the log of the minimum virial temperature of star-forming haloes to within 45.3/22.0/16.7, 33.5/18.4/17.8 and 6.3/3.3/2.4 per cent, ~$1\sigma$ fractional uncertainty, respectively. Instead, if we optimistically assume that we can perfectly characterize the EoR modelling uncertainties, we can improve on these constraints by up to a factor of ~few. Similarly, the fractional uncertainty on the average neutral fraction can be constrained to within $\lesssim 10$ per cent for HERA and SKA1. By studying the resulting impact on astrophysical constraints, 21CMMC can be used to optimize (i) interferometer designs; (ii) foreground cleaning algorithms; (iii) observing strategies; (iv) alternative statistics characterizing the 21 cm signal; and (v) synergies with other observational programs.

Accepted to MNRAS, matches online version. 20 pages, 4 tables, 11 colour figures

References in corpus (20)

• A luminous quasar at a redshift of z = 7.085
• Cosmology at Low Frequencies: The 21 cm Transition and the High-Redshift Universe
• Efficient Simulations of Early Structure Formation and Reionization
• The morphology of HII regions during reionization
• The observed ionization rate of the intergalactic medium and the ionizing emissivity at z >5: Evidence for a photon starved and extended epoch of reionization
• On the (Non)Evolution of HI Disks over Cosmic Time
• The Epoch of Reionization Window: I. Mathematical Formalism
• New observations of z~7 galaxies: evidence for a patchy reionization
• Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests
• The Epoch of Reionization Window: II. Statistical Methods for Foreground Wedge Reduction
• Can the intergalactic medium cause a rapid drop in Lyman alpha emission at z>6?
• Detecting the Rise and Fall of 21 cm Fluctuations with the Murchison Widefield Array
• The Keck+Magellan Survey for Lyman Limit Absorption I: The Frequency Distribution of Super Lyman Limit Systems
• The Physics and Early History of the Intergalactic Medium
• Was Star-Formation Suppressed in High-Redshift Minihalos?
• UV Radiative Feedback on High-Redshift Proto-Galaxies
• Higher Order Contributions to the 21 cm Power Spectrum
• Early galaxy formation in warm dark matter cosmologies
• Lyman continuum galaxies and the escape fraction of Lyman break galaxies
• Adding Context to JWST Surveys with Current and Future 21cm Radio Observations

Cited by in corpus (24)

• Inferences on the Timeline of Reionization at z~8 From the KMOS Lens-Amplified Spectroscopic Survey
• Line-Intensity Mapping: 2017 Status Report
• Improving the Epoch of Reionization Power Spectrum Results from Murchison Widefield Array Season 1 Observations
• Constraining the intergalactic medium at $z\approx$ 9.1 using LOFAR Epoch of Reionization observations
• Emulation of reionization simulations for Bayesian inference of astrophysics parameters using neural networks
• Mitigating Internal Instrument Coupling for 21 cm Cosmology I: Temporal and Spectral Modeling in Simulations
• Mitigating Internal Instrument Coupling II: A Method Demonstration with the Hydrogen Epoch of Reionization Array
• Bubble size statistics during reionization from 21-cm tomography
• The impact of interference excision on 21-cm Epoch of Reionization power spectrum analyses
• Reionization inference from the CMB optical depth and E-mode polarization power spectra
• 21SSD: a public database of simulated 21-cm signals from the epoch of reionization
• Constraining the reionization history using deep learning from 21cm tomography with the Square Kilometre Array
• Epoch of Reionisation 21cm Forecasting From MCMC-Constrained Semi-Numerical Models
• Contamination of the Epoch of Reionization power spectrum in the presence of foregrounds
• The impact of non-Gaussianity on the error covariance for observations of the Epoch of Reionization 21-cm power spectrum
• The Escape Fraction of Ionizing Photons During the Epoch of Reionization: observability with the Square Kilometre Array
• The impact of non-Gaussianity on the Epoch of Reionization parameter forecast using 21-cm power spectrum measurements
• A unified framework for 21cm tomography sample generation and parameter inference with Progressively Growing GANs
• Fundamental uncertainty levels of 21cm power spectra from a delay analysis
• Bayesian Model Selection with Future 21cm Observations of The Epoch of Reionisation
• Constraining the Reionization History using Bayesian Normalizing Flows
• Upper Limits on the 21 cm Power Spectrum at z = 5.9 from Quasar Absorption Line Spectroscopy
• Parameters Estimation from the 21 cm signal using Variational Inference
• It's Always Darkest Before the Cosmic Dawn: Early Results from Novel Tools and Telescopes for 21 cm Cosmology