Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate ("SFR") between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8.

Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RAMSES-CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star-formation algorithm, making reionization end earlier, at $z \gtrsim 6$, thereby better matching the observations of intergalactic Lyman-alpha opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at $z = 6$. As in CoDa I, reionization feedback suppresses star formation in haloes below $\sim 2 \times 10^9$ M$_{\odot}$, though suppression here is less severe, a possible consequence of modifying the star-formation algorithm. Suppression is environment-dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of $\Lambda$CDM constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, nor by nearby groups, nor by each other.

We present a deep learning model trained to emulate the radiative transfer during the epoch of cosmological reionization. CRADLE (Cosmological Reionization And Deep LEarning) is an autoencoder convolutional neural network that uses two-dimensional maps of the star number density and the gas density field at z=6 as inputs and that predicts 3D maps of the times of reionization $\mathrm{t_{reion}}$ as outputs. These predicted single fields are sufficient to describe the global reionization history of the intergalactic medium in a given simulation. We trained the model on a given simulation and tested the predictions on another simulation with the same paramaters but with different initial conditions. The model is successful at predicting $\mathrm{t_{reion}}$ maps that are in good agreement with the test simulation. We used the power spectrum of the $\mathrm{t_{reion}}$ field as an indicator to validate our model. We show that the network predicts large scales almost perfectly but is somewhat less accurate at smaller scales. While the current model is already well-suited to get average estimates about the reionization history, we expect it can be further improved with larger samples for the training, better data pre-processing and finer tuning of hyper-parameters. Emulators of this kind could be systematically used to rapidly obtain the evolving HII regions associated with hydro-only simulations and could be seen as precursors of fully emulated physics solvers for future generations of simulations.

The reduced speed of light approximation is used in a variety of simulations of the epoch of reionization and galaxy formation. Its popularity stems from its ability to drastically reduce the computing cost of a simulation, by allowing the use of larger, and therefore fewer timesteps to reach a solution. It is physically motivated by the fact that ionization fronts rarely propagate faster than some fraction of the speed of light. However, no global proof of the physical validity of this approach is available, and possible artefacts resulting from this approximation therefore need to be identifited and characterized to allow its proper use. In this paper we investigate the impact of the reduced speed of light approximation on the predicted properties of the intergalactic medium. To this end we use fully coupled radiation-hydrodynamics RAMSES-CUDATON simulations of the epoch of reionization. We find that reducing the speed of light by a factor 5 (20, 100) leads to overestimating the post-reionization volume-weighted $x_{HI}$ by a similar factor ~5 (20, 100) with respect to full speed of light simulations. We show that the error is driven by the hydrogen - photon chemistry. In photo-ionization equilibrium, reducing the speed of light has the same effect as artificially reducing the photon density or the reaction cross-section and leads to an underestimated ionizing flux. We confirm this interpretation by running additional simulations using a reduced speed of light in the photon propagation module, but keeping this time the full speed of light in the chemistry module. With this setup, the post-reionization neutral hydrogen fractions converge to the full speed of light value, which validates our explanation. Increasing spatial resolution beyond a cell size of 1 kpc physical, so as to better resolve Lyman-limit systems, does not significantly affect our conclusions.

Coupled radiative-hydrodynamics simulations of the epoch of reionization aim to reproduce the propagation of ionization fronts during the transition before the overlap of HII regions. Many of these simulations use moment-based methods to track radiative transfer processes using explicit solvers and are therefore subject to strict stability conditions regarding the speed of light, which implies a great computational cost. It can be reduced by assuming a reduced speed of light, and this approximation is now widely used to produce large-scale simulations of reionization. We introduce a new method for estimating and comparing the ionization front speeds based on maps of the reionization redshifts. We applied it to a set of cosmological simulations of the reionization using a set of reduced speeds of light, and measured the evolution of the ionization front speeds during the reionization process. We find that ionization fronts progress via a two-stage process, the first stage at low velocity as the fronts emerge from high density regions and a second later stage just before the overlap, during which front speeds increase close to the speed of light. Using a set of small 8Mpc/h^3 simulations, we find that a minimal velocity of 0.3c is able to model these two stages in this specific context without significant impact. Values as low as 0.05c can model the first low velocity stage, but limit the acceleration at later times. Lower values modify the distribution of front speeds at all times. Using larger 64Mpc/h^3 volumes that better account for distant sources, we find that reduced speed of light has a greater impact on reionization times and front speeds in underdense regions that are reionized at late times. The same quantities measured in dense regions with slow fronts are less sensitive to c values.

Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen to form the well-known structures of the local universe, including the Local Group and Virgo, in a (91 Mpc)$^3$ volume large enough to model both global and local reionization. Reionization simulation CoDa I-AMR, by CPU-GPU code EMMA, used (2048)$^3$ particles and (2048)$^3$ initial cells, adaptively-refined, while N-body simulation CoDa I-DM2048, by Gadget2, used (2048)$^3$ particles, to find reionization times for all galaxies at z = 0 with masses $M(z=0)\ge 10^8 M_\odot$. Galaxies with $M(z=0) \gtrsim 10^{11} M_\odot$ reionized earlier than the universe as a whole, by up to $\sim$ 500 Myrs, with significant scatter. For Milky-Way-like galaxies, $z_R$ ranged from 8 to 15. Galaxies with $M(z=0) \lesssim 10^{11} M_\odot$ typically reionized as late or later than globally-averaged $50\%$ reionization at $\langle z_R\rangle =7.8$, in neighborhoods where reionization was completed by external radiation. The spread of reionization times within galaxies was sometimes as large as the galaxy-to-galaxy scatter. The Milky Way and M31 reionized earlier than global reionization but later than typical for their mass, neither dominated by external radiation. Their most massive progenitors at $z>6$ had $z_R$ = 9.8 (MW) and 11 (M31), while their total masses had $z_R$ = 8.2 (both).

EMMA is a cosmological simulation code aimed at investigating the reionization epoch. It handles simultaneously collisionless and gas dynamics, as well as radiative transfer physics using a moment-based description with the M1 approximation. Field quantities are stored and computed on an adaptive 3D mesh and the spatial resolution can be dynamically modified based on physically-motivated criteria. Physical processes can be coupled at all spatial and temporal scales. We also introduce a new and optional approximation to handle radiation : the light is transported at the resolution of the non-refined grid and only once the dynamics have been fully updated, whereas thermo-chemical processes are still tracked on the refined elements. Such an approximation reduces the overheads induced by the treatment of radiation physics. A suite of standard tests are presented and passed by EMMA, providing a validation for its future use in studies of the reionization epoch. The code is parallel and is able to use graphics processing units (GPUs) to accelerate hydrodynamics and radiative transfer calculations. Depending on the optimizations and the compilers used to generate the CPU reference, global GPU acceleration factors between x3.9 and x16.9 can be obtained. Vectorization and transfer operations currently prevent better GPU performances and we expect that future optimizations and hardware evolution will lead to greater accelerations.