We present an analysis of the evolution of the Lyman-series forest into the epoch of reionization using cosmological radiative transfer simulations in a scenario where reionization ends late. We explore models with different midpoints of reionization and gas temperatures. We find that once the simulations have been calibrated to match the mean flux of the observed Lyman-$\alpha$ forest at $4 < z < 6$, they also naturally reproduce the distribution of effective optical depths of the Lyman-$\beta$ forest in this redshift range. We note that the tail of the largest optical depths that is most challenging to match corresponds to the long absorption trough of ULAS J0148+0600, which we have previously shown to be rare in our simulations. We consider the evolution of the Lyman-series forest out to higher redshifts, and show that future observations of the Lyman-$\beta$ forest at $z>6$ will discriminate between different reionization histories. The evolution of the Lyman-$\alpha$ and Lyman-$\gamma$ forests are less promising as a tool for pushing studies of reionization to higher redshifts due to the stronger saturation and foreground contamination, respectively.

We show that, in the context of patchy reionisation, an accurate description of the angular power spectrum of the kinetic Sunyaev-Zel'dovich (kSZ) effect is not possible with simple scaling relations between the amplitude of the spectrum and global parameters, such as the reionisation midpoint and its duration. We introduce a new parameterisation of this spectrum, based on a novel description of the power spectrum of the free electrons density contrast Pee (k,z) in terms of the reionisation global history and morphology. We directly relate features of the spectrum to the typical ionised bubble size at different stages in the process and, subsequently, to the angular scale at which the patchy kSZ power spectrum reaches its maximum. We successfully calibrated our results on a custom set of advanced radiative hydrodynamical simulations and later found our parameterisation to be a valid description of a wide range of other simulations and, therefore, reionisation physics. In the end, and as long as the global reionisation history is known, two parameters are sufficient to derive the angular power spectrum. Such an innovative framework applied to cosmic microwave background data and combined with 21cm intensity mapping will allow a first consistent detection of the amplitude and shape of the patchy kSZ signal, giving in turn access to the physics of early light sources.

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.

The decline in abundance of Lyman-$\alpha$ (Ly$\alpha$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$\alpha$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$\alpha$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$\alpha$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 5\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$\alpha$ peak. We discuss our results in the context of the interpretation of observations.

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.

A long (110 cMpc/$h$) and deep absorption trough in the Ly$\alpha$ forest has been observed extending down to redshift 5.5 in the spectrum of ULAS J0148+0600. Although no Ly$\alpha$ transmission is detected, Ly$\beta$ spikes are present which has led to claims that the gas along this trough must be ionized. Using high resolution cosmological radiative transfer simulations in large volumes, we show that in a scenario where reionization ends late ($z \sim 5.2$), our simulations can reproduce troughs as long as observed. In this model, we find that the troughs are caused by islands of neutral hydrogen. Small ionized holes within the neutral islands allow for the transmission of Ly$\beta$. We have also modelled the Ly$\alpha$ emitter population around the simulated troughs, and show that there is a deficit of Ly$\alpha$ emitters close to the trough as is observed.

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.

High-redshift QSO spectra show large spatial fluctuations in the Ly-alpha opacity of the intergalactic medium on surprisingly large scales at z>~5.5. We present a radiative transfer simulation of cosmic reionization driven by galaxies that reproduces this large scatter and the rapid evolution of the Ly-alpha opacity distribution at 5<z<6. The simulation also reproduces the low Thomson scattering optical depth reported by the latest CMB measurement and is consistent with the observed short near-zones and strong red damping wings in the highest-redshift QSOs. It also matches the rapid disappearance of observed Ly-alpha emission by galaxies at z>~6. Reionization is complete at z=5.3 in our model, and 50% of the volume of the Universe is ionized at z=7. Agreement with the Ly-alpha forest data in such a late reionization model requires a rapid evolution of the ionizing emissivity of galaxies that peaks at z~6.8. The late end of reionization results in a large scatter in the photoionisation rate and the neutral hydrogen fraction at redshifts as low as z<~5.5 with large residual neutral 'islands' that can produce very long Gunn-Peterson troughs resembling those seen in the data.

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.

Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical simulation of reionization and galaxy formation in the Local Universe, in a volume large enough to model reionization globally but with enough resolving power to follow all atomic-cooling galactic halos in that volume. For every halo identified at a given time, we find the redshift at which the surrounding IGM reionized, along with its instantaneous star formation rate (`SFR') and baryonic gas-to-dark matter ratio ($M_\text{gas}/M_\text{DM}$). The average SFR per halo with $M < 10^9 \text{ M}_\odot$ was steady in regions not yet reionized, but declined sharply following local reionization. For $M > 10^{10} \text{ M}_\odot$, this SFR continued through local reionization, increasing with time, instead. For $10^9 < M < 10^{10} \text{ M}_\odot$, the SFR generally increased modestly through reionization, followed by a modest decline. In general, halo SFRs were higher for regions that reionized earlier. A similar pattern was found for $M_\text{gas}/M_\text{DM}$, which declined sharply following local reionization for $M < 10^9 \text{ M}_\odot$. Local reionization time correlates with local matter overdensity, which determines the local rates of structure formation and ionizing photon consumption. The earliest patches to develop structure and reionize ultimately produced more stars than they needed to finish and maintain their own reionization, exporting their `surplus' starlight to help reionize regions that developed structure later.

We calibrate here cosmological radiative transfer simulation with ATON/RAMSES with a range of measurements of the Lyman alpha opacity from QSO absorption spectra. We find the Lyman alpha opacity to be very sensitive to the exact timing of hydrogen reionisation. Models reproducing the measured evolution of the mean photoionisation rate and average mean free path reach overlap at z ~ 7 and predict an accelerated evolution of the Lyman alpha opacity at z > 6 consistent with the rapidly evolving luminosity function of Lyman alpha emitters in this redshift range. Similar to "optically thin" simulations our full radiative transfer simulations fail, however, to reproduce the high-opacity tail of the Lyman alpha opacity PDF at z > 5. We argue that this is due to spatial UV fluctuations in the post-overlap phase of reionisation on substantially larger scales than predicted by our source model, where the ionising emissivity is dominated by large numbers of sub-L* galaxies. We further argue that this suggests a significant contribution to the ionising UV background by much rarer bright sources at high redshift.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

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).

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

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.

We search for vast planes of satellites (VPoS) in a high resolution simulation of the Local Group performed by the CLUES project, which improves significantly the resolution of former similar studies. We use a simple method for detecting planar configurations of satellites, and validate it on the known plane of M31. We implement a range of prescriptions for modelling the satellite populations, roughly reproducing the variety of recipes used in the literature, and investigate the occurence and properties of planar structures in these populations. The structure of the simulated satellite systems is strongly non-random and contains planes of satellites, predominantly co-rotating, with, in some cases, sizes comparable to the plane observed in M31 by Ibata et al.. However the latter is slightly richer in satellites, slightly thinner and has stronger co-rotation, which makes it stand out as overall more exceptional than the simulated planes, when compared to a random population. Although the simulated planes we find are generally dominated by one real structure, forming its backbone, they are also partly fortuitous and are thus not kinematically coherent structures as a whole. Provided that the simulated and observed planes of satellites are indeed of the same nature, our results suggest that the VPoS of M31 is not a coherent disc and that one third to one half of its satellites must have large proper motions perpendicular to the plane.

We use high resolution simulations of the formation of the local group post-processed by a radiative transfer code for UV photons, to investigate the reionization of the satellite populations of an isolated Milky Way-M31 galaxy pair in a variety of scenarios. We use an improved version of ATON which includes a simple recipe for radiative feedback. In our baseline models, reionization is initiated by low mass, radiatively regulated haloes at high redshift, until more massive haloes appear, which then dominate and complete the reionization process. We investigate the relation between reionization history and present-day positions of the satellite population. We find that the average reionization redshift (zr) of satellites is higher near galaxy centers (MW and M31). This is due to the inside-out reionization patterns imprinted by massive haloes within the progenitor during the EoR, which end up forming the center of the galaxy. Thanks to incomplete dynamical mixing during galaxy assembly, these early patterns survive down to present day, resulting in a a clear radial gradient in the average satellites reionization redshift, up to the virial radius of MW and M31 and beyond. In the lowest emissivity scenario, the outer satellites are reionized about 180 Myr later than the inner satellites. This delay decreases with increasing source model emissivity, or in the case of external reionization by Virgo or M31, because reionization happens faster overall, and becomes spatially quasi-uniform at the highest emissivity.

We constrain the initial stage of the reionization process around progenitors of galaxies, such as the extent of the initial HII region before its fusion with the UV background and the duration of its propagation. We use a set of reionisation simulations with different resolutions and ionizing source recipes. A catalog of the HII regions properties is built thanks to a merger tree of HII regions. We draw local reionization histories as a function of time and investigate variations according to the halo mass progenitors of the regions. We then extrapolate the halo mass inside the region from high z to z=0 to make predictions about the reionization histories of z=0 galaxies. We found that the later an HII region appears, the smaller will be its related lifetime and volume before it sees the global UV background. Quantitatively the duration and the extent of the initial growth of an HII region is strongly dependent on the mass of the inner halo and can be as long as 50 % of the reionization epoch. We found that the most massive is a halo today, the earlier it appears and the larger are the extension and the duration of propagation of its HII region. Quantitative predictions differ depending on the box size or the source model: small simulated volumes are affected by proximity effects between HII regions and halo-based source models predict smaller regions and slower I-front expansion than in models using star particles as ionizing sources. Our results suggests that Milky Way-type halos have a maximal extent of 1.1 Mpc/h for the initial HII region that established itself in 150-200$\pm 20$ Myrs. This is consistent with prediction made using constrained Local Group simulation. Considering halos with masses comparable to those of the Local Group (MW+M31), our result suggests that statistically it has not been influenced by an external front coming from a Virgo-like cluster.

We present the results of a set of numerical simulations aimed at studying reionization at galactic scale. We use a high resolution simulation of the formation of the Milky Way-M31 system to simulate the reionization of the local group. The reionization calculation was performed with the post-processing radiative transfer code ATON and the underlying cosmological simulation was performed as part of the CLUES project. We vary the source models to bracket the range of source properties used in the literature. We investigate the structure and propagation of the galatic ionization fronts by a visual examination of our reionization maps. Within the progenitors we find that reionization is patchy, and proceeds locally inside out. The process becomes patchier with decreasing source photon output. It is generally dominated by one major HII region and 1-4 additional isolated smaller bubbles, which eventually overlap. Higher emissivity results in faster and earlier local reionization. In all models, the reionization of the Milky Way and M31 are similar in duration, i.e. between 203 Myr and 22 Myr depending on the source model, placing their zreion between 8.4 and 13.7. In all models except the most extreme, the MW and M31 progenitors reionize internally, ignoring each other, despite being relatively close to each other even during the epoch of reionization. Only in the case of strong supernova feedback suppressing star formation in haloes less massive than 10^9 M_sun, and using our highest emissivity, we find that the MW is reionized by M31.

We present a new implementation of radiation hydrodynamics (RHD) in the adaptive mesh refinement (AMR) code RAMSES. The multi-group radiative transfer (RT) is performed on the AMR grid with a first-order Godunov method using the M1 closure for the Eddington tensor, and is coupled to the hydrodynamics via non-equilibrium thermochemistry of hydrogen and helium. This moment-based approach has the large advantage that the computational cost is independent of the number of radiative sources - it can even deal with continuous regions of emission such as bound-free emission from gas. As it is built directly into RAMSES, the RT takes natural advantage of the refinement and parallelization strategies already in place. Since we use an explicit advection solver for the radiative transport, the time step is restricted by the speed of light - a severe limitation that can be alleviated using the so--called "reduced speed of light" approximation. We propose a rigorous framework to assess the validity of this approximation in various conditions encountered in cosmology and galaxy formation. We finally perform with our newly developed code a complete suite of RHD tests, comparing our results to other RHD codes. The tests demonstrate that our code performs very well and is ideally suited for exploring the effect of radiation on current scenarios of structure and galaxy formation.

We describe a new methodology to analyze the reionization process in numerical simulations: The evo- lution of the reionization is investigated by focusing on the merger histories of individual HII regions. From the merger tree of ionized patches, one can track the individual evolution of the regions such as e.g. their size, or investigate the properties of the percolation process by looking at the formation rate, the frequency of mergers and the number of individual HII regions involved in the mergers. By applying this technique to cosmological simulations with radiative transfer, we show how this methodology is a good candidate to quantify the impact of the star formation adopted on the history of the reionization. As an application we show how different models of sources result in different evolutions and geometry of the reionization even though they produce e.g. similar ionized fraction or optical depth.

We describe a new methodology to analyze the reionization process in numerical simulations: the chronology and the geometry of reionization is investigated by means of merger histories of individual HII regions. From the merger tree of ionized patches, one can track the individual evolution of the regions properties such as e.g. their size, or the intensity of the percolation process by looking at the formation rate, the frequency of mergers and the number of individual HII regions involved in the mergers. We apply the merger tree technique to simulations of reionization with three different kinds of ionizing source models and two resolutions. Two of them use star particles as ionizing sources. In this case we confront two emissivity evolutions for the sources in order to reach the reionization at z ~ 6. As an alternative we built a semi-analytical model where the dark matter halos extracted from the density fields are assumed as ionizing sources. We then show how this methodology is a good candidate to quantify the impact of the adopted star formation on the history of the observed reionization. The semi-analytical model shows a homogeneous reionization history with 'local' hierarchical growth steps for individual HII regions. On the other hand auto-consistent models for star formation tend to present fewer regions with a dominant region in size which governs the fusion process early in the reionization at the expense of the 'local' reionizations. The differences are attenuated when the resolution of the simulation is increased.

We present a new semi-analytical model of the population of satellite galaxies of the Milky Way, aimed at estimating the effect of the geometry of reionization at galaxy scale on the properties of the satellites. In this model reionization can be either: (A) externally-driven and uniform, or (B) internally-driven, by the most massive progenitor of the Milky Way. In the latter scenario the propagation of the ionisation front and photon dilution introduce a delay in the photo-evaporation of the outer satellites' gas with respect to the inner satellites. As a consequence, outer satellites experience a longer period of star formation than those in the inner halo. We use simple models to account for star formation, the propagation of the ionisation front, photo-evaporation and observational biases. Both scenarios yield a model satellite population that matches the observed luminosity function and mass-to-light ratios. However, the predicted population for scenario (B) is significantly more extended spatially than for scenario (A), by about 0.3 dex in distance, resulting in a much better match to the observations. The survival of the signature left by the local UV field during reionization on the radial distribution of satellites makes it a promising tool for studying the reionization epoch at galaxy scale in the Milky Way and nearby galaxies resolved in stars with forthcoming large surveys.

We present a new semi-analytical model of the population of satellite galaxies of the Milky Way, aimed at estimating the effect of the spatial structure of reionisation at galaxy scale on the properties of the satellites. In this model reionisation can be either: (A) externally-driven and uniform, or (B) internally-driven, by the most massive progenitor of the Milky-Way. In the latter scenario the propagation of the ionisation front and photon dilution introduce a delay in the photo-evaporation of the outer satellites' gas with respect to the inner satellites. As a consequence, outer satellites experience a longer period of star formation than those in the inner halo. We use simple models to account for star formation, the propagation of the ionisation front, photo-evaporation and observational biases. Both scenarios yield a model satellite population at z=0 that matches the observed luminosity function and mass-to-light ratios. However, the predicted population for scenario (B) is significantly more extended spatially than for scenario (A), by about 0.3 dex in distance, resulting in a much better match to the observations. The survival of this structural signature imprinted by the local UV field during reionisation on the radial distribution of satellites makes it a promising tool for studying the reionisation epoch at galaxy scale in the Milky Way and nearby galaxies resolved in stars with forthcoming large surveys. However, more work is needed to determine how the effect reported here can be disentangled from that of cosmic variance between different realisations of Milky Way haloes.

We present a set of cosmological simulations with radiative transfer in order to model the reionization history of the Universe. Galaxy formation and the associated star formation are followed self-consistently with gas and dark matter dynamics using the RAMSES code, while radiative transfer is performed as a post-processing step using a moment-based method with M1 closure relation in the ATON code. The latter has been ported to a multiple Graphics Processing Units (GPU) architecture using CUDA + MPI, resulting in an overall acceleration (x80) that allows us to tackle radiative transfer problems at resolution of 1024^3 + 2 levels of refinement for the hydro adaptive grid and 1024^3 for the RT cartesian grid. We observe a good convergence between our different resolution runs as long as the effects of finite resolution on the star formation history are properly taken into account. We also show that the neutral fraction depends on the total mass density, in a way close to the predictions of photoionization equilibrium, as long as the effect of self-shielding is included in the background radiation model. However we still fail at reproducing the z=6 constraints on the H neutral fraction and the intensity of the UV background. In order to account for unresolved density fluctuations, we added a simple clumping factor model. Using our most spatially resolved simulation (12.5 Mpc/h-1024^3) to calibrate our subgrid model, we have resimulated our largest box (100 Mpc/h 1024^3), successfully reproducing the observed level of H neutral fraction at z=6. We don't reproduce the photoionization rate inferred from the same observations. We argue that this discrepancy could be explained by the fact that the average radiation intensity and the average neutral fraction depends on different regions of the gas density distribution, so that one quantity cannot be simply deduced from the other.

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

We look for observational signatures that could discriminate between Newtonian and modified Newtonian (MOND) dynamics in the Milky Way, in view of the advent of large astrometric and spectroscopic surveys. Indeed, a typical signature of MOND is an apparent disk of "phantom" dark matter, which is uniquely correlated with the visible disk-density distribution. Due to this phantom dark disk, Newtonian models with a spherical halo have different signatures from MOND models close to the Galactic plane. The models can thus be differentiated by measuring dynamically (within Newtonian dynamics) the disk surface density at the solar radius, the radial mass gradient within the disk, or the velocity ellipsoid tilt angle above the Galactic plane. Using the most realistic possible baryonic mass model for the Milky Way, we predict that, if MOND applies, the local surface density measured by a Newtonist will be approximately 78 Msun/pc2 within 1.1 kpc of the Galactic plane, the dynamically measured disk scale-length will be enhanced by a factor of 1.25 with respect to the visible disk scale-length, and the local vertical tilt of the velocity ellipsoid at 1 kpc above the plane will be approximately 6 degrees. None of these tests can be conclusive for the present-day accuracy of Milky Way data, but they will be of prime interest with the advent of large surveys such as GAIA.

We have performed a 70 billion dark-matter particles N-body simulation in a 2 $h^{-1}$ Gpc periodic box, using the concordance, cosmological model as favored by the latest WMAP3 results. We have computed a full-sky convergence map with a resolution of $\Delta \theta \simeq 0.74$ arcmin$^{2}$, spanning 4 orders of magnitude in angular dynamical range. Using various high-order statistics on a realistic cut sky, we have characterized the transition from the linear to the nonlinear regime at $\ell \simeq 1000$ and shown that realistic galactic masking affects high-order moments only below $\ell < 200$. Each domain (Gaussian and non-Gaussian) spans 2 decades in angular scale. This map is therefore an ideal tool for testing map-making algorithms on the sphere. As a first step in addressing the full map reconstruction problem, we have benchmarked in this paper two denoising methods: 1) Wiener filtering applied to the Spherical Harmonics decomposition of the map and 2) a new method, called MRLens, based on the modification of the Maximum Entropy Method on a Wavelet decomposition. While the latter is optimal on large spatial scales, where the signal is Gaussian, MRLens outperforms the Wiener method on small spatial scales, where the signal is highly non-Gaussian. The simulated full-sky convergence map is freely available to the community to help the development of new map-making algorithms dedicated to the next generation of weak-lensing surveys.

We present a statistical study of the effects induced by substructures on the deflection potential of dark matter halos in the strong lensing regime. This investigation is based on the pertubative solution around the Einstein radius (Alard 2007) in which all the information on the deflection potential is specified by only a pair of one-dimensional functions on this ring. Using direct comparison with ray-tracing solutions, we found that the iso-contours of lensed images predicted by the pertubative solution is reproduced with a mean error on their radial extension of less than 1% - in units of the Einstein radius, for reasonable substructure masses. It demonstrates the efficiency of the approximation to track possible signatures of substructures. We have evaluated these two fields and studied their properties for different lens configurations modelled either through massive dark matter halos from a cosmological N-body simulation, or via toy models of Monte Carlo distribution of substructures embedded in a triaxial Hernquist potential. As expected, the angular power spectra of these two fields tend to have larger values for larger harmonic numbers when substructures are accounted for and they can be approximated by power-laws, whose values are fitted as a function of the profile and the distribution of the substructures.

This technical paper describes a software package that was designed to produce initial conditions for large cosmological simulations in the context of the Horizon collaboration. These tools generalize E. Bertschinger's Grafic1 software to distributed parallel architectures and offer a flexible alternative to the Grafic2 software for ``zoom'' initial conditions, at the price of large cumulated cpu and memory usage. The codes have been validated up to resolutions of 4096^3 and were used to generate the initial conditions of large hydrodynamical and dark matter simulations. They also provide means to generate constrained realisations for the purpose of generating initial conditions compatible with, e.g. the local group, or the SDSS catalog.

A radiative transfer scheme is presented, based on a moment description of the equation of radiative transfer and the so-called ``M1 closure model'' for the Eddington tensor. This model features a strictly hyperbolic transport step for radiation: it has been implemented using standard Godunov--like techniques in a new code called ATON. Coupled to simple models of ionization chemistry and photo-heating, ATON is able to reproduce the results of other schemes on a various set of standard tests such as the expansion of a HII region, the shielding of the radiation by dense clumps and cosmological ionization by multiple sources. Being simple yet robust, such a scheme is intended to be naturally and easily included in grid--based cosmological fluid solvers.

In a serie of three papers, the dynamical interplay between environments and dark matter haloes is investigated, while focussing on the dynamical flows through their virial sphere. Our method relies on both cosmological simulations, to constrain the environments, and an extension to the classical matrix method to derive the response of the halo (see Pichon & Aubert (2006), paper I). The current paper focuses on the statistical characterisation of the environments surrounding haloes, using a set of large scale simulations. Our description relies on a `fluid' halocentric representation where the interactions between the halo and its environment are investigated in terms of a time dependent external tidal field and a source term characterizing the infall. The method is applied to 15000 haloes, with masses between 5 x 10^12 Ms and 10^14 Ms evolving between z = 1 and z = 0. The net accretion at the virial radius is found to decrease with time, resulting from both an absolute decrease of infall and from a growing contribution of outflows. Infall is found to be mainly radial and occurring at velocities ~ 0.75 V200. Outflows are also detected through the virial sphere and occur at lower velocities ~ 0.6 V200 on more circular orbits. The external tidal field is found to be strongly quadrupolar and mostly stationnary, possibly reflecting the distribution of matter in the halo's near environment. The coherence time of the small scale fluctuations of the potential hints a possible anisotropic distribution of accreted satellites. The flux density of mass on the virial sphere appears to be more clustered than the potential while the shape of its angular power spectrum seems stationnary.

We present a numerical technique to compute the gravitational lensing induced by simulated haloes. It relies on a 2D-Tree domain decomposition in the lens plane combined with a description of N-Body particles as extended clouds with a non-singular density. This technique is made fully adaptive by the use of a density-dependent smoothing which allows one to probe the lensing properties of haloes from the densest regions in the center or in substructures to the low-density regions in the outskirts. 'Smooth Particle Lensing' has some promising features. First, the deflection potential, the deflection angles, the convergence and the shear are direct and separate end-products of the SPL calculation and can be computed at an arbitrary distribution of points on the lens plane. Second, this flexibility avoids the use of interpolation or a finite differentiation procedure on a grid, does not require padding the region with zeros and focuses the computing power on relevant regions. The SPL algorithm is tested by populating isothermal spheres and ellipsoids with particles and then comparing the lensing calculations to the classical FFT-based technique and analytic solutions. We assess issues related to the resolution of the lensing code and the limitations set by the simulations themselves. We conclude by discussing how SPL can be used to predict the impact of substructures on strong lensing and how it can be generalized to weak lensing and cosmic shear simulations.

We investigate statistically the dynamical consequences of cosmological fluxes of matter and related moments on progenitors of today's dark matter haloes. Their dynamics is described via canonical perturbation theory which accounts for two types of perturbations: the tidal field corresponding to fly-bys and accretion of dark matter through the halo's outer boundary. he dynamical equations are solved linearly, order by order, projecting on a biorthogonal basis to consistently satisfy the field equation. Since our solution of the Boltzmann Poisson equations is explicit, it allows statistical predictions for the ensemble distribution of the inner dynamical features of haloes. The secular evolution of open galactic haloes is investigated: we derive the kinetic equation which governs the quasi-linear evolution of dark matter profile induced by infall and its corresponding gravitational correlations. This yields a Fokker Planck-like equation for the angle-averaged underlying distribution function. We show how these extensions to the classical theory could be used to (i) observationally constrain the statistical nature of the infall (ii) predict the observed distribution and correlations of substructures in upcoming surveys, (iii) predict the past evolution of the observed distribution of clumps, and finally (iv) weight the relative importance of the intrinsic (via the unperturbed distribution function) and external (tidal and/ or infall) influence of the environment in determining the fate of galaxies.

We measure the anisotropy of dark matter flows on small scales in the near environment of haloes using a large set of simulations. We rely on two different approaches to quantify the anisotropy of the cosmic infall: we measure the flows at the haloes' virial radius while describing the infalling matter via fluxes through a spherical shell; we measure the spatial and kinematical distributions of satellites and substructures around haloes detected by the subclump finder ADAPTAHOP first described in the appendix. The two methods are found to be in agreement both qualitatively and quantitatively via one and two points statistics.The infall takes place preferentially in the plane perpendicular to the direction defined by the halo's spin. We computed the excess of equatorial accretion both through rings and via a harmonic expansion of the infall. The level of anisotropy of infalling matter is found to be ~15 %. The substructures have their spin orthogonal to their velocity vector in the halo's rest frame at a level of about 5%, suggestive of an image of a flow along filamentary structures which provides an explanation for the measured anisotropy. We conclude that a halo does not see its environment as an isotropic perturbation, investigate how the anisotropy is propagated inwards using perturbation theory, and discuss briefly implications for weak lensing, warps and the thickness of galactic disks.

We measure the spatial clustering of metals in the intergalactic medium from z = 1.7 to 3.0, as traced by 643 CIV and 104 SiIV N >= 10^12 cm^-2 absorption systems in 19 high signal-to-noise (40-80) and high resolution (R = 45000) quasar spectra. The number densities and two-point correlation functions of both these species are largely constant with redshift, suggesting the bulk of metal ejection occurred at z >= 3. However, at z <= 1.9 some additional signature appears in the CIV correlation function at 500 km/s, associated with four strong and peculiar systems. At all redshifts, the CIV and SiIV correlation functions exhibit a steep rise at large separations and a flatter profile at small separations, with an elbow occurring at ~150 km/s. We show that these properties are consistent with metals confined within bubbles with a typical radius R_s about sources of mass >= M_s, and use numerical simulations to derive best-fit values of R_s ~ 2 comoving Mpc and M_s ~ 5x10^11 solar masses at z = 3. This does not exclude that metals could have been produced at higher redshifts in smaller, but equally rare, objects. At the level of detection of this survey, IGM enrichment is likely to be incomplete and inhomogeneous, with a filling factor ~ 10%.