• ### French SKA White Book - The French Community towards the Square Kilometre Array(1712.06950)

March 28, 2018 astro-ph.IM
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.
• ### The impact of the reduced speed of light approximation on the post-overlap neutral hydrogen fraction in numerical simulations of the epoch of reionization(1803.02434)

March 6, 2018 astro-ph.CO, astro-ph.GA
Context: 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. Aims: In this paper we investigate the impact of the reduced speed of light approximation on the predicted properties of the intergalactic medium. Methods: To this end we use fully coupled radiation-hydrodynamics RAMSES-CUDATON simulations of the epoch of reionization. Results: We find that reducing the speed of light by a factor 5 (20, 100) leads to overestimating the post-reionization average volume-weighted neutral hydrogen fraction 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 by considering the analytical solution for a strongly ionized hydrogen gas in photo-ionization equilibrium. In this regime, 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 true speed of light in the chemistry module. With this setup, the post-reionization neutral hydrogen fractions converge to the true value, which validates our explanation.
• 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.
• ### Vast planes of satellites in a high resolution simulation of the Local Group: comparison to Andromeda(1412.3110)

Jan. 5, 2015 astro-ph.GA
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.
• ### The reionization of galactic satellite populations(1409.8285)

Sept. 29, 2014 astro-ph.CO, astro-ph.GA
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.
• ### High resolution simulations of the reionization of an isolated Milky Way - M31 galaxy pair(1310.0009)

Sept. 30, 2013 astro-ph.CO, astro-ph.GA
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.
• ### Constraining local UV field geometry at reionization using Milky Way satellites(1111.1663)

Nov. 7, 2011 astro-ph.CO, astro-ph.GA
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.
• ### Fitting galaxy spectra with STECKMAP: a user guide(1108.4631)

Aug. 23, 2011 astro-ph.GA, astro-ph.IM
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.
• ### Fake star formation bursts: blue horizontal branch stars masquerade as young stars in optical integrated light spectroscopy(0911.3156)

Nov. 16, 2009 astro-ph.CO
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.
• ### Bimodal gas accretion in the Horizon-MareNostrum galaxy formation simulation(0803.4506)

Aug. 7, 2008 astro-ph
(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.
• ### Spectroscopic ages and metallicities of stellar populations: validation of full spectrum fitting(0801.0871)

Jan. 6, 2008 astro-ph
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.