Exploring the Model Dependence of MCMC-based 21 cm Power Spectrum Parameter Constraints

Abstract Detection and analysis of the cosmic 21 cm signal of neutral hydrogen has long been considered the most promising route towards exploration of the Epoch of Reionization (EoR). 21CMMC, a Markov Chain Monte Carlo sampler of the seminumerical simulation code 21cmFAST, has already been used in conjunction with published upper limits on the 21 cm signal from the Murchison Widefield Array, the Low Frequency Array, and the Hydrogen Epoch of Reionization Array to constrain the astrophysics of the EoR. Here, we investigate the extent to which analysis of the EoR performed using 21CMMC is dependent on the underlying seminumerical model. We used 21cmFAST to simulate two datasets of 21 cm light-cones that differ only in the algorithm used to identify ionized regions (the so-called “bubble-finding” algorithm). We then tested 21CMMC’s ability to return key astrophysical parameters when using the different bubble-finding algorithms. We find that the performance of 21CMMC depends sensitively on the agreement between the astrophysical model of our mock data and the model used for sampling. This result has important implications for the analysis of the 21 cm signal performed using 21CMMC and further motivates investigation into model-independent analysis techniques for 21 cm EoR data.