Simulating cosmic reionization at large scales - I. The geometry of reionization

We present the first large-scale radiative transfer simulations of cosmic reionization, in a simulation volume of (100 h(-1) Mpc)(3). This is more than a two 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 efficient 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 which 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 similar to 21, and final overlap occurs by z similar to 11. The resulting electron-scattering optical depth is in good agreement with the first-year Wilkinson Microwave Anisotropy Probe (WMAP) polarization data. We show that reionization clearly proceeded in an inside-out fashion, with the high-density regions being ionized earlier, on average, than the voids. Ionization histories of smaller-size (5-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 evolution of the mean ionized fraction. The minimum reliable volume size for such predictions is similar to 30 Mpc. We derive the power spectra of the neutral, ionized and total gas density fields and show that there is a significant boost of the density fluctuations in both the neutral and the ionized components relative to the total at arcmin and larger scales. We find two populations of H II regions according to their size, numerous, mid-sized (similar to 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 Mpc is reached by z similar to 13, when our volume is only about 50 per cent ionized, and well before the global overlap. We derive the statistical distributions of the ionized fraction and ionized gas density 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 alpha emitter statistics, Gunn-Peterson optical depth and small-scale cosmic microwave background secondary anisotropies due to patchy reionization.