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

We present a new method for modeling inhomogeneous cosmic reionization on large scales. Utilizing high-resolution radiation-hydrodynamic simulations with 2048^3 dark matter particles, 2048^3 gas cells, and 17 billion adaptive rays in a L = 100 Mpc/h box, we show that the density and reionization-redshift fields are highly correlated on large scales (>~ 1 Mpc/h). This correlation can be statistically represented by a scale-dependent linear bias. We construct a parametric funct...

The cosmic microwave background provides an image of the Universe 0.4 million years after the big bang, when atomic hydrogen formed out of free electrons and protons. One of the primary goals of observational cosmology is to obtain follow-up images of the Universe during the epoch of reionization, hundreds of millions of years later, when cosmic hydrogen was ionized once again by the UV photons emitted from the first galaxies. To achieve this goal, new observatories are being...

Redshifted 21cm measurements of the structure of ionised regions that grow during reionization promise to provide a new probe of early galaxy and structure formation. One of the challenges of modelling reionization is to account both for the sub-halo scale physics of galaxy formation and the regions of ionization on scales that are many orders of magnitude larger. To bridge this gap we first calculate the statistical relationship between ionizing luminosity and Mpc-scale over...

The transformation of cold neutral intergalactic hydrogen into a highly ionized warm plasma marks the end of the cosmic dark ages and the beginning of the age of galaxies. The details of this process reflect the nature of the early sources of radiation and heat, the statistical characteristics of the large-scale structure of the Universe, the thermodynamics and chemistry of cosmic baryons, and the histories of star formation and black hole accretion. A number of massive data ...

We introduce the THESAN project, a suite of large volume (L = 95.5 cMpc) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the IGM during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter (DM) and baryonic mass resolutions of $3.1 \times 10^6$ $M_\odot$ and $5.8 \times 10^5 $ $M_\odot$, respectively. The gravitational forces are softened on scale...

The universe goes through several phase transitions during its formative stages. Cosmic reionization is the last of them, where ultraviolet and X-ray radiation escape from the first generations of galaxies heating and ionizing their surroundings and subsequently the entire intergalactic medium. There is strong observational evidence that cosmic reionization ended approximately one billion years after the Big Bang, but there are still uncertainties that will be clarified with ...

In this first of several application papers, we investigate the mechanics of reionization from stellar sources in high-z galaxies, the utility of various clumping factors on estimating the recombination time in the IGM, and the photon budget required to achieve reionization. We test the accuracy of the static and time-dependent models of Madau et al. as predictors of reionization completion/maintenance. We simulate a WMAP7 LCDM cosmological model in a 20 Mpc comoving cube wit...

We study the evolution of ionization fronts around the first proto-galaxies by using high resolution numerical cosmological (Lambda+CDM model) simulations and Monte Carlo radiative transfer methods. We present the numerical scheme in detail and show the results of test runs from which we conclude that the scheme is both fast and accurate. As an example of interesting cosmological application, we study the reionization produced by a stellar source of total mass M=2 10^8 M_\odo...

Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models that are pre-computed via the initial density field, or post-processed where feedb...

Large-scale polarization of the cosmic microwave background measured by the WMAP satellite requires a mean optical depth to Thomson scattering, tau_e ~ 0.17. The reionization of the universe must therefore have begun at relatively high redshift. We have studied the reionization process using supercomputer simulations of a large and representative region of a universe which has cosmological parameters consistent with the WMAP results. Our simulations follow both the radiative ...