Fully coupled simulation of cosmic reionization. II. recombinations, clumping factors, and the photon budget for reionization

Abstract

We use a fully coupled cosmological simulation including dark matter dynamics, multispecies hydrodynamics, nonequilibrium chemical ionization, flux-limited diffusion radiation transport, and a parameterized model of star formation and feedback (thermal and radiative) to investigate the epoch of hydrogen reionization in detail. In this paper, the first of several application papers, we investigate the mechanics of reionization from stellar sources forming in high-z galaxies, the utility of various formulations for the gas clumping factor on accurately estimating the effective recombination time in the intergalactic medium (IGM), and the photon budget required to achieve reionization. We also test the accuracy of the static and time-dependent models of Madau et al. as predictors of reionization completion/maintenance. We simulate a WMAP7 ΛCDM cosmological model in a 20 comoving Mpc cube, resolved with 8003 uniform fluid cells and dark matter particles. By tuning our star formation recipe to approximately match the observed high-redshift star formation rate density and galaxy luminosity function, we have created a fully coupled radiation hydrodynamical realization of hydrogen reionization, which begins to ionize at z ≈ 10 and is completed at z ≈ 5.8 without further tuning. We find that roughly two ionizing photons per H atom are required to convert the neutral IGM to a highly ionized state. After reionization concludes, we find that the quantity ṅion×(1 Gyr)/nH is ∼9 at z = 5, in rough agreement with measurements of the ionizing emissivity by Becker & Bolton. The complicated events during reionization that lead to this number can be generally described as inside-out, but in reality, the narrative depends on the level of ionization of the gas one attributes as being ionized. We find that the formula for the ionizing photon production rate needed to maintain the IGM in an ionized state derived by Madau et al. should not be used to predict the epoch of reionization completion because it ignores history-dependent terms in the global ionization balance which are not ignorable. We find that the time-dependent model for the ionized volume fraction QH II is more predictive, but overestimates the redshift of reionization completion zreion by Δz ≈ 1. We propose a revised formulation of the time-dependent model that agrees with our simulation to a few percent accuracy. Finally, we use our simulation data to measure the absorption of ionizing radiation due to circumgalactic gas resolved on our mesh and find f̄esc(CGM) ≈ 0.7. © 2014. The American Astronomical Society. All rights reserved.