Simulating Cosmic Reionization at Large Scales I: the Geometry of Reionization

TLDR: The first large-scale radiative transfer simulation of cosmic reionization was presented in this paper, in a simulation volume of (100 h 1 Mpc) 3, which is more than a 2 orders of magnitude improvement over previous simulations.

Abstract: We present the first large-scale radiative transfer simulat ions of cosmic reionization, in a simulation volume of (100 h 1 Mpc) 3 . This is more than a 2 orders of magnitude improvement over previous simulations. We achieve this by combining the results from extremely large, cosmological, N-body simulations with a new, fast and effici ent code for 3D radiative transfer, C 2 -Ray, which we have recently developed. These simulations allow us to do the first numerical studies of the large-scale structure of reionization w hich at the same time, and crucially, properly take account of the dwarf galaxy ionizing sources which are primarily responsible for reionization. In our realization, reionization starts around z � 21, and final overlap occurs by z � 11. The resulting electron-scattering optical depth is in goo d agreement with the firstyear WMAP polarization data. We show that reionization clearly proceeded in an inside-out fashion, with the high-density regions being ionized earli er, on average, than the voids. Ionization histories of smaller-size (5 to 10 comoving Mpc) subregions exabit a large scatter about the mean and do not describe the global reionization history well. This is true even when these subregions are at the mean density of the universe, which shows that small-box simulations of reionization have little predictive power for the evolut ion of the mean ionized fraction. The minimum reliable volume size for such predictions is � 30 Mpc. We derive the power-spectra of the neutral, ionized and total gas density fields and show t hat there is a significant boost of the density fluctuations in both the neutral and the ionized c omponents relative to the total at arcminute and larger scales. We find two populations of H II re gions according to their size, numerous, mid-sized (� 10 Mpc) regions and a few, rare, very large regions tens of Mpc in size. Thus, local overlap on fairly large scales of tens of Mp c is reached by z � 13, when our volume is only about 50% ionized, and well before the global overlap. We derive the statistical distributions of the ionized fraction and ionized gas densi ty at various scales and for the first time show that both distributions are clearly non-Gaussian. All these quantities are critical for predicting and interpreting the observational signals from reionization from a variety of observations like 21-cm emission, Ly-α emitter statistics, Gunn-Peterson optical depth and small-scale CMB secondary anisotropies due to patchy reionization.