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.

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.

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 present the SPHINX suite of cosmological adaptive mesh refinement simulations, the first radiation-hydrodynamical simulations to simultaneously capture large-scale reionization and the escape of ionizing radiation from thousands of resolved galaxies. Our $5$ and $10$ co-moving Mpc volumes resolve haloes down to the atomic cooling limit and model the inter-stellar medium with better than $\approx10$ pc resolution. The project has numerous goals in improving our understanding of reionization and making predictions for future observations. In this first paper we study how the inclusion of binary stars in computing stellar luminosities impacts reionization, compared to a model that includes only single stars. Owing to the suppression of galaxy growth via strong feedback, our galaxies are in good agreement with observational estimates of the galaxy luminosity function. We find that binaries have a significant impact on the timing of reionization: with binaries, our boxes are $99.9$ percent ionized by volume at $z\approx 7$, while without them our volumes fail to reionize by $z=6$. These results are robust to changes in volume size, resolution, and feedback efficiency. The escape of ionizing radiation from individual galaxies varies strongly and frequently. On average, binaries lead to escape fractions of $\approx 7-10$ percent, about $3.5$ times higher than with single stars only. The higher escape fraction is a result of a shallower decline in ionizing luminosity with age, and is the primary reason for earlier reionization, although the higher integrated luminosity with binaries also plays a sub-dominant role.

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

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

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

STECKMAP stands for STEllar Content and Kinematics via Maximum A Posteriori likelihood. It is a tool for interpreting galaxy spectra in terms of their stellar populations, through the derivation of their star formation history, age-metallicity relation, kinematics and extinction. To do so, the observed spectrum is projected onto a temporal sequence of models of single stellar populations, so as to determine a linear combination of these models, that fit the observed spectrum best. The weights of the various components of this linear combination indicate the stellar content of the population. This procedure is regularized using various penalizing functions. The principles of the method are detailed in Ocvirk et al. 2006a,b. The STECKMAP software package is public and freely available at http://astro.u-strasbg.fr/~ocvirk/. A number of authors have already adopted it and use it in their daily research. This user guide aims at accompanying the user through the setup and first runs of STECKMAP. The last chapter will help the user to understand and improve his results and experience with the code.

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 the first results of a pilot study aimed at understanding the influence of bars on the evolution of galaxy discs through the study of their stellar content. We examine here the kinematics, star formation history, mass-weighted, luminosity-weighted, and single stellar population (SSP) equivalent ages and metallicities for four galaxies ranging from lenticulars to late-type spirals. The data employed extends to 2-3 disc scalelengths, with S/N(A)>50. Several techniques are explored to derive star formation histories and SSP-equivalent parameters, each of which are shown to be compatible. We demostrate that the age-metallicity degeneracy is highly reduced by using spectral fitting techniques --instead of indices-- to derive these parameters. We found that the majority of the stellar mass in our sample is composed of old (~10 Gyr) stars. This is true in the bulge and the disc region, even beyond two disc scalelengths. In the bulge region, we find that the young, dynamically cold, structures produced by the presence of the bar (e.g., nuclear discs or rings) are responsible for shaping the bulges' age and metallicity gradients. In the disc region, a larger fraction of young stars is present in the external parts of the disc compared with the inner disc. The disc growth is, therefore, compatible with a moderate inside-out formation scenario, where the luminosity weighted age changes from ~10 Gyrs in the centre, to ~4 Gyrs at two disc scalelengths, depending upon the galaxy. For two galaxies, we compare the metallicity and age gradients of the disc major axis with that of the bar, finding very important differences. In particular, the stellar population of the bar is more similar to the bulge than to the disc, indicating that, at least in those two galaxies, bars formed long ago and have survived to the present day. (abridged)

In a galaxy cluster, the evolution of spiral galaxies depends on their cluster environment. Ram pressure due to the rapid motion of a spiral galaxy within the hot intracluster medium removes the galaxy's interstellar medium from the outer disk. Once the gas has left the disk, star formation stops. The passive evolution of the stellar populations should be detectable in optical spectroscopy and multi-wavelength photometry. The goal of our study is to recover the stripping age of the Virgo spiral galaxy NGC 4388, i.e. the time elapsed since the halt of star formation in the outer galactic disk using a combined analysis of optical spectra and photometry. We performed VLT FORS2 long-slit spectroscopy of the inner star-forming and outer gas-free disk of NGC 4388. We developed a non-parametric inversion tool that allows us to reconstruct the star formation history of a galaxy from spectroscopy and photometry. The tool was tested on a series of mock data using Monte Carlo simulations. The results from the non-parametric inversion were refined by applying a parametric inversion method. The star formation history of the unperturbed galactic disk is flat. The non-parametric method yields a rapid decline of star formation < 200 Myr ago in the outer disk. The parametric method is not able to distinguish between an instantaneous and a long-lasting star formation truncation. The time since the star formation has dropped by a factor of two from its pre-stripping value is 190 +- 30 Myr. We are able to give a precise stripping age that is consistent with revised dynamical models.

Model color magnitude diagrams of low-metallicity globular clusters usually show a deficit of hot evolved stars with respect to observations. We investigate quantitatively the impact of such modelling inaccuracies on the significance of star formation history reconstructions obtained from optical integrated spectra. To do so, we analyse the sample of spectra of galactic globular clusters of Schiavon et al. with STECKMAP (Ocvirk et al.) and the stellar population models Vazdekis et al. and Bruzual & Charlot, and focus on the reconstructed stellar age distributions. Firstly, we show that background/foreground contamination correlates with E(B-V), which allows us to define a clean subsample of uncontaminated GCs, on the basis of a E(B-V) filtering. We then identify a "confusion zone" where fake young bursts of star formation pop up in the star formation history although the observed population is genuinely old. These artifacts appear for 70-100% of cases depending on the population model used, and contribute up to 12% of the light in the optical. Their correlation with the horizontal branch ratio indicates that the confusion is driven by HB morphology: red horizontal branch clusters are well fitted by old stellar population models while those with a blue HB require an additional hot component. The confusion zone extends over [Fe/H]=[-2,-1.2], although we lack the data to probe extreme high and low metallicity regimes. As a consequence, any young starburst superimposed on an old stellar population in this metallicity range could be regarded as a modeling artifact, if it weighs less than 12% of the optical light, and if no emission lines typical of an HII region are present.

[Abridged] This paper aims at providing new conservative constraints to the cosmic star-formation history from the empirical modeling of mid- and far-infrared data. We perform a non-parametric inversion of galaxy counts at 15, 24, 70, 160, and 850 microns simultaneously. It is a "blind" search (no redshift information is required) of all possible evolutions of the infrared luminosity function of galaxies, from which the evolution of the star-formation rate density and its uncertainties are derived. The cosmic infrared background (CIRB) measurements are used a posteriori to tighten the range of solutions. The inversion relies only on two hypotheses: (1) the luminosity function remains smooth both in redshift and luminosity, (2) a set of infrared spectral energy distributions (SEDs) of galaxies must be assumed. The range of star-formation histories that we derive is well constrained and consistent with redshift-based measurements from deep surveys. The redshift decompositions of the counts are also recovered successfully. Therefore, multi-wavelength counts and CIRB (both projected observations) alone seem to contain enough information to recover the cosmic star-formation history with quantifiable errors. A peak of the SFRD at z~2 is preferred, although higher redshifts are not excluded. We also find a good consistency between the observed evolution of the stellar mass density and the prediction from our model. Finally, the inability of the inversion to model perfectly and simultaneously all the multi-wavelength infrared counts (especially at 160 microns where an excess is seen around 20 mJ) implies either (i) the existence of a sub-population of colder galaxies, (ii) a larger dispersion of dust temperatures among local galaxies than expected, (iii) or a redshift evolution of the infrared SEDs of galaxies.

One of the complexities in modelling integrated spectra of stellar populations is the effect of interacting binary stars besides type Ia supernovae (SNeIa). These include common envelope systems, CVs, novae, and are usually ignored in models predicting the chemistry and spectral absorption line strengths in galaxies. In this paper predictions of chemical yields from populations of single and binary stars are incorporated into a galactic chemical evolution model to explore the significance of the effects of these other binary yields. Effects on spectral line strengths from different progenitor channels of SNeIa are also explored. Small systematic effects are found when the yields from binaries, other than SNeIa, are included, for a given star formation history. These effects are, at present, within the observational uncertainties on the line strengths. More serious differences can arise in considering different types of SNIa models, their rates and contributions.

(abridged) The physics of diffuse gas accretion and the properties of the cold and hot modes of accretion onto proto-galaxies between z=2 and z=5.4 is investigated using the large cosmological simulation performed with the RAMSES code on the MareNostrum supercomputing facility. Galactic winds, chemical enrichment, UV background heating and radiative cooling are taken into account in this very high resolution simulation. Using accretion--weighted temperature histograms, we have perfomed two different measurements of the thermal state of the gas accreted towards the central galaxy. The first measurement, performed using accretion--weighted histograms on a spherical surface of radius 0.2 Rvir centred on the densest gas structure near the halo centre of mass, is a good indicator of the presence of an accretion shock in the vicinity of the galactic disc. We define the hot shock mass, Mshock, as the typical halo mass separating cold dominated from hot dominated accretion in the vicinity of the galaxy. The second measurement is performed by radially averaging histograms between 0.2 Rvir and Rvir, in order to detect radially extended structures such as gas filaments: this is a good proxy for detecting cold streams feeding the central galaxy. We define Mstream as the transition mass separating cold dominated from hot dominated accretion in the outer halo, marking the disappearance of these cold streams. We find a hot shock transition mass of Mshock=10^{11.6} Msun (dark matter), with no significant evolution with redshift. Conversely, we find that Mstream increases sharply with z. Our measurements are in agreement with the analytical predictions of Birnboim & Dekel (2003) and Dekel & Birnboim (2006), if we correct their model by assuming low metallicity (< 10^{-3} Zsun) for the filaments, correspondingly to our measurements.

Fitting whole spectra at intermediate spectral resolution (R = 1000 -- 3000), to derive physical properties of stellar populations, appears as an optimized alternative to methods based on spectrophotometric indices: it uses all the redundant information contained in the signal. This paper addresses the validation of the method and it investigates the quality of the population models together with the reliability of the fitting procedures. We are using two algorithms: STECKMAP, a non-parametric regularized program and NBURSTS a parametric non-linear minimization. We compare three spectral synthesis models for single stellar populations: Pegase-HR, Galaxev (BC03) and Vazdekis/Miles, and we analyse spectra of Galactic clusters whose populations are known from studies of color-magnitude diagrams (CMD) and spectroscopy of individual stars. We find that: (1) The quality of the models critically depends on the stellar library they use. Pegase-HR and Vazdekis/Miles are consistent, while the comparison between Pegase-HR and BC03 shows some systematics reflecting the limitations of the stellar library (STELIB) used to generate the latter models; (2) The two fitting programs are consistent; (3) For globular clusters and M67 spectra, the method restitutes metallicities in agreement with spectroscopy of stars within 0.14 dex; (4) The spectroscopic ages are very sensitive to the presence of a blue horizontal branch (BHB) or of blue stragglers. A BHB morphology results in a young SSP-equivalent age. Fitting a free amount of blue stars in addition to the SSP model to mimic the BHB improves and stabilizes the fit and restores ages in agreement with CMDs studies. This method is potentially able to disentangle age or BHB effects in extragalactic clusters.

The stellar population models dramatically progressed with the arrival of large and complete libraries, ELODIE, CFLIB (=Indo-US) and MILES at a relatively high resolution. We show that the quality of the fits is not anymore limited by the size of the stellar libraries in a large range of ages (0.1 to 10 Gyrs) and metallicities (-2 to +0.4 dex). The main limitations of the empirical stellar libraries are (i) the coverage of the parameters space (lack of hot stars of low metallicity), (ii) the precision and homogeneity of the atmospheric parameters and (iii) the non-resolution of individual element abundances (in particular [$\alpha$/Fe]). Detailed abundance measurements in the large libraries, and usage of theoretical libraries are probably the next steps, and we show that a combination between an empirical (ELODIE) and a theoretical library (Coelho et al. 2005) immediately improves the modeling of ($\alpha$-enhanced) globular clusters.

We have compared simple stellar populations (SSPs) generated with different population synthesis tools: BC03, Vazdekis and Pegase.HR and different stellar libraries: ELODIE3.1, SteLib and MILES. We find that BC03/SteLib SSPs are biased toward solar metallicity, however Pegase.HR/ELODIE3.1 and Vazdekis/MILES are extremely consistent. The extensive coverage of the space of atmospheric parameters in the large stellar libraries allows precise synthesis for a large range of ages (0.1 .. 10 Gyr) and metallicities (-2 .. +0.4 dex) limited by the quality of the determination of stellar parameters (like temperature scale of the giants).

We are using full spectrum fitting to determine the ages and metallicities of Galactic clusters (M67 and globular clusters). We find that the method is very accurate to measure the metallicity. Blue horizontal branches are well identified as a 'young' sub-component, and the age of the 'old' component is in fair agreement with CMD determinations.

We introduce STECKMAP (STEllar Content and Kinematics via Maximum A Posteriori), a method to recover the kinematical properties of a galaxy simultaneously with its stellar content from integrated light spectra. It is an extension of STECMAP (astro-ph/0505209) to the general case where the velocity distribution of the underlying stars is also unknown. %and can be used as is for the analysis of large sets of data. The reconstructions of the stellar age distribution, the age-metallicity relation, and the Line-Of-Sight Velocity Distribution (LOSVD) are all non-parametric, i.e. no specific shape is assumed. The only a propri we use are positivity and the requirement that the solution is smooth enough. The smoothness parameter can be set by GCV according to the level of noise in the data in order to avoid overinterpretation. We use single stellar populations (SSP) from PEGASE-HR (R=10000, lambda lambda = 4000-6800 Angstrom, Le Borgne et al. 2004) to test the method through realistic simulations. Non-Gaussianities in LOSVDs are reliably recovered with SNR as low as 20 per 0.2 Angstrom pixel. It turns out that the recovery of the stellar content is not degraded by the simultaneous recovery of the kinematic distribution, so that the resolution in age and error estimates given in Ocvirk et al. 2005 remain appropriate when used with STECKMAP. We also explore the case of age-dependent kinematics (i.e. when each stellar component has its own LOSVD). We separate the bulge and disk components of an idealized simplified spiral galaxy in integrated light from high quality pseudo data (SNR=100 per pixel, R=10000), and constrain the kinematics (mean projected velocity, projected velocity dispersion) and age of both components.

This paper describes STECMAP (STEllar Content via Maximum A Posteriori), a flexible, non-parametric inversion method for the interpretation of the integrated light spectra of galaxies, based on synthetic spectra of single stellar populations (SSPs). We focus on the recovery of a galaxy's star formation history and stellar age-metallicity relation. We use the high resolution SSPs produced by PEGASE-HR to quantify the informational content of the wavelength range 4000 - 6800 Angstroms. A detailed investigation of the properties of the corresponding simplified linear problem is performed using singular value decomposition. It turns out to be a powerful tool for explaining and predicting the behaviour of the inversion. We provide means of quantifying the fundamental limitations of the problem considering the intrinsic properties of the SSPs in the spectral range of interest, as well as the noise in these models and in the data. We performed a systematic simulation campaign and found that, when the time elapsed between two bursts of star formation is larger than 0.8 dex, the properties of each episode can be constrained with a precision of 0.04 dex in age and 0.02 dex in metallicity from high quality data (R=10 000, signal-to-noise ratio SNR=100 per pixel), not taking model errors into account. The described methods and error estimates will be useful in the design and in the analysis of extragalactic spectroscopic surveys.

One way to extract more information from the integrated light of galaxies is to improve the spectral resolution at which observations and analysis are carried out. The population synthesis code currently providing the highest spectral resolution is Pegase-HR, which was made available by D. Le Borgne et al. in 2004. Based on an empirical stellar library, it provides synthetic spectra between 4000 and 6800 A at lambda/d(lambda)=10000 for any star formation history, with or without chemical evolution. Such a resolution is particularly useful for the study of low mass galaxies, massive star clusters, or other galaxy regions with low internal velocity dispersions. After a summary of the main features of Pegase-HR and comparisons with other population synthesis codes, this paper focuses on the inversion of optical galaxy spectra. We explore the limits of what information can or can not be recovered, based on theoretical principles and extensive simulations. First applications to extragalactic objects are shown.

We present inversion techniques aiming at recovering the stellar content (age, metallicity, extinction) and kinematics of a stellar population from its absorption line spectrum. These techniques use new synthetic high resolution spectra produced by PEGASE_HR and minimization algorithms. Application to SDSS data and to a simulated disc-bulge system are presented as examples.