Cosmic Reionization Redux

Nickolay Y. Gnedin,Xiaohui Fan

DOI: https://doi.org/10.1086/505790

2006-03-29

Abstract:We show that numerical simulations of reionization that resolve the Lyman Limit systems (and, thus, correctly count absorptions of ionizing photons) have converged to about 10% level for 5<z<6.2 and are in reasonable agreement (within 10%) with the SDSS data in this redshift interval. The SDSS data thus constraint the redshift of overlap of cosmic HII regions to z_{OVL} = 6.1+-0.15. At higher redshifts, the simulations are far from convergence on the mean Gunn-Peterson optical depth, but achieve good convergence for the mean neutral hydrogen fraction. The simulations that fit the SDSS data, however, do not have nearly enough resolution to resolve the earliest episodes of star formation, and are very far from converging on the precise value of the optical depth to Thompson scattering - any value between 6 and 10% is possible, depending on the convergence rate of the simulations and the fractional contribution of PopIII stars. This is generally consistent with the third-year WMAP results, but much higher resolution simulation are required to come up with the sufficiently precise value for the Thompson optical depth that can be statistically compared with the WMAP data.

Astrophysics

What problem does this paper attempt to address?

This paper attempts to solve some key problems in the cosmic reionization period, especially the degree of agreement between the simulation of the reionization process and the observational data. Specifically, the paper mainly focuses on the following aspects:

1. **Convergence of numerical simulations**: The author shows that numerical simulations capable of resolving the Lyman Limit system have converged to approximately the 10% level in the redshift range \(5 < z < 6.2\), and these simulation results are basically consistent with SDSS (Sloan Digital Sky Survey) data (with an error within 10%). This indicates that in this redshift range, the simulation can describe the reionization process relatively accurately.

2. **Determination of the overlap redshift**: By comparing with SDSS data, the paper determines that the redshift of the overlap of HII regions in the universe is \(z_{\text{OVL}} = 6.1 \pm 0.15\). This redshift marks the time point when the reionization process is nearing its end.

3. **Simulation challenges in the high - redshift range**: In the higher redshift range (\(z > 6.2\)), the simulation has not yet converged in terms of the average Gunn - Peterson optical depth. Nevertheless, the simulation has better convergence in terms of the average neutral hydrogen fraction, which indicates that the simulation is still reliable in terms of certain physical quantities.

4. **Estimation of Thompson scattering optical depth**: The paper points out that the current simulation resolution is not sufficient to resolve the earliest star - formation events, and therefore it is unable to accurately calculate the Thompson scattering optical depth \(\tau_T\). According to the simulation results, the value of \(\tau_T\) may be between 6% and 10%, depending on the convergence rate of the simulation and the contribution ratio of Pop III stars.

### Summary
The main objective of the paper is to evaluate whether the existing numerical simulations can accurately describe the cosmic reionization process and to verify the validity of these simulations by comparing them with observational data (such as SDSS and WMAP). The study finds that in the redshift range \(5 < z < 6.2\), the simulation is basically consistent with the observational data, but there are still challenges in the higher redshift range and in the accurate calculation of the Thompson scattering optical depth. Higher - resolution simulations are required in the future to further improve the understanding of the reionization process.