Fully coupled simulation of cosmic reionization. II. Recombinations, clumping factors, and the photon budget for reionization
- CASS, University of California, San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0424 (United States)
- Southern Methodist University, 6425 Boaz Lane, Dallas, TX 75205 (United States)
- Center for Relativistic Astrophysics, Georgia Institute of Technology, 837 State St, Atlanta, GA 30332 (United States)
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 800{sup 3} 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 n-dot{sub ion}×(1 Gyr)/n{sub H} is ∼9 at z = 5, in rough agreement with measurements of the ionizing emissivity by Becker and 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 Q{sub HII} is more predictive, but overestimates the redshift of reionization completion z{sub reion} 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-bar{sub esc}(CGM)≈0.7.
- OSTI ID:
- 22365685
- Journal Information:
- Astrophysical Journal, Vol. 789, Issue 2; Other Information: Country of input: International Atomic Energy Agency (IAEA); ISSN 0004-637X
- Country of Publication:
- United States
- Language:
- English
Similar Records
HIGH-RESOLUTION SIMULATIONS OF THE REIONIZATION OF AN ISOLATED MILKY WAY-M31 GALAXY PAIR
THE OPACITY OF THE INTERGALACTIC MEDIUM DURING REIONIZATION: RESOLVING SMALL-SCALE STRUCTURE