Showing 1–50 of 55 results for author: Ocvirk, P

First upper limits on the 21 cm signal power spectrum from cosmic dawn from one night of observations with NenuFAR

Authors: S. Munshi, F. G. Mertens, L. V. E. Koopmans, A. R. Offringa, B. Semelin, D. Aubert, R. Barkana, A. Bracco, S. A. Brackenhoff, B. Cecconi, E. Ceccotti, S. Corbel, A. Fialkov, B. K. Gehlot, R. Ghara, J. N. Girard, J. M. Grießmeier, C. Höfer, I. Hothi, R. Mériot, M. Mevius, P. Ocvirk, A. K. Shaw, G. Theureau, S. Yatawatta , et al. (2 additional authors not shown)

Abstract: The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 2… ▽ More The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift of $z$ = 20.3. NenuFAR is a new low-frequency radio interferometer, operating in the 10-85 MHz frequency range, currently under construction at the Nançay Radio Observatory in France. It is a phased array instrument with a very dense uv coverage at short baselines, making it one of the most sensitive instruments for 21 cm cosmology analyses at these frequencies. Our analysis adopts the foreground subtraction approach, in which sky sources are modeled and subtracted through calibration and residual foregrounds are subsequently removed using Gaussian process regression. The final power spectra are constructed from the gridded residual data cubes in the uv plane. Signal injection tests are performed at each step of the analysis pipeline, the relevant pipeline settings are optimized to ensure minimal signal loss, and any signal suppression is accounted for through a bias correction on our final upper limits. We obtain a best 2$σ$ upper limit of $2.4\times 10^7$ $\text{mK}^{2}$ at $z$ = 20.3 and $k$ = 0.041 $h\,\text{cMpc}^{-1}$. We see a strong excess power in the data, making our upper limits two orders of magnitude higher than the thermal noise limit. We investigate the origin and nature of this excess power and discuss further improvements to the analysis pipeline that can potentially mitigate it and consequently allow us to reach thermal noise sensitivity when multiple nights of observations are processed in the future. △ Less

Submitted 30 April, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 27 pages, 21 figures, and 6 tables; accepted for publication in Astronomy and Astrophysics (A&A); language edits implemented; typos corrected

Journal ref: A&A 681, A62 (2024)

The Sphinx Public Data Release: Forward Modelling High-Redshift JWST Observations with Cosmological Radiation Hydrodynamics Simulations

Authors: Harley Katz, Joki Rosdahl, Taysun Kimm, Jeremy Blaizot, Nicholas Choustikov, Marion Farcy, Thibault Garel, Martin G. Haehnelt, Leo Michel-Dansac, Pierre Ocvirk

Abstract: The recent launch of JWST has ushered in a new era of high-redshift astronomy by providing detailed insights into the gas and stellar populations of galaxies in the epoch of reionization. Interpreting these observations and translating them into constraints on the physics of early galaxy formation is a complex challenge that requires sophisticated models of star formation and the interstellar medi… ▽ More The recent launch of JWST has ushered in a new era of high-redshift astronomy by providing detailed insights into the gas and stellar populations of galaxies in the epoch of reionization. Interpreting these observations and translating them into constraints on the physics of early galaxy formation is a complex challenge that requires sophisticated models of star formation and the interstellar medium (ISM) in high-redshift galaxies. To this end, we present Version 1 of the Sphinx$^{20}$ public data release. Sphinx$^{20}$ is a full box cosmological radiation hydrodynamics simulation that simultaneously models the large-scale process of cosmic reionization and the detailed physics of a multiphase ISM, providing a statistical sample of galaxies akin to those currently being observed by JWST. The data set contains $\sim14,000$ mock images and spectra of the stellar continuum, nebular continuum, and 52 nebular emission lines, including Ly$α$, for each galaxy in Sphinx$^{20}$ with a star formation rate $\geq0.3\ {\rm M_{\odot}\ yr^{-1}}$. All galaxy emission has been processed with dust radiative transfer and/or resonant line radiative transfer, and data is provided for ten viewing angles for each galaxy. Additionally, we provide a comprehensive set of intrinsic galaxy properties, including halo masses, stellar masses, star formation histories, and ISM characteristics (e.g., metallicity, ISM gas densities, LyC escape fractions). This paper outlines the data generation methods, presents a comparative analysis with JWST ERS and Cycle 1 observations, and addresses data set limitations. The Sphinx$^{20}$ data release can be downloaded at the following URL: https://github.com/HarleyKatz/SPHINX-20-data △ Less

Submitted 3 December, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: 33 pages, 34 figures, Accepted to The Open Journal of Astrophysics, data can be downloaded at the following URL: https://github.com/HarleyKatz/SPHINX-20-data

New insight on the nature of cosmic reionizers from the CEERS survey

Authors: S. Mascia, L. Pentericci, A. Calabrò, P. Santini, L. Napolitano, P. Arrabal Haro, M. Castellano, M. Dickinson, P. Ocvirk, J. S. W. Lewis, R. Amorín, M. Bagley, R. N. J. Cleri, L. Costantin, A. Dekel, S. L. Finkelstein, A. Fontana, M. Giavalisco, N. A. Grogin, N. P. Hathi, M. Hirschmann, B. W. Holwerda, I. Jung, J. S. Kartaltepe, A. M. Koekemoer , et al. (7 additional authors not shown)

Abstract: The Epoch of Reionization (EoR) began when galaxies grew in abundance and luminosity, so their escaping Lyman continuum (LyC) radiation started ionizing the surrounding neutral intergalactic medium (IGM). Despite significant recent progress, the nature and role of cosmic reionizers are still unclear: in order to define them, it would be necessary to directly measure their LyC escape fraction (… ▽ More The Epoch of Reionization (EoR) began when galaxies grew in abundance and luminosity, so their escaping Lyman continuum (LyC) radiation started ionizing the surrounding neutral intergalactic medium (IGM). Despite significant recent progress, the nature and role of cosmic reionizers are still unclear: in order to define them, it would be necessary to directly measure their LyC escape fraction ($f_{esc}$). However, this is impossible during the EoR due to the opacity of the IGM. Consequently, many efforts at low and intermediate redshift have been made to determine measurable indirect indicators in high-redshift galaxies so that their $f_{esc}$ can be predicted. This work presents the analysis of the indirect indicators of 62 spectroscopically confirmed star-forming galaxies at $6 \leq z \leq 9$ from the Cosmic Evolution Early Release Science (CEERS) survey, combined with 12 sources with public data from other JWST-ERS campaigns. From the NIRCam and NIRSpec observations, we measured their physical and spectroscopic properties. We discovered that on average $6<z<9$ star-forming galaxies are compact in the rest-frame UV ($r_e \sim $ 0.4 kpc), are blue sources (UV-$β$ slope $\sim $ -2.17), and have a predicted $f_{esc}$ of about 0.13. A comparison of our results to models and predictions as well as an estimation of the ionizing budget suggests that low-mass galaxies with UV magnitudes fainter than $M_{1500} = -18$ that we currently do not characterize with JWST observations probably played a key role in the process of reionization. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 14 pages, 11 figures, submitted to A&A

Reionisation time fields reconstruction from 21 cm signal maps

Authors: Julien Hiegel, Emilie Thélie, Dominique Aubert, Jonathan Chardin, Nicolas Gillet, Pierre Galois, Nicolas Mai, Pierre Ocvirk, Rodrigo Ibata

Abstract: During the Epoch of reionisation, the intergalactic medium is reionised by the UV radiation from the first generation of stars and galaxies. One tracer of the process is the 21 cm line of hydrogen that will be observed by the Square Kilometre Array (SKA) at low frequencies, thus imaging the distribution of ionised and neutral regions and their evolution. To prepare for these upcoming observations,… ▽ More During the Epoch of reionisation, the intergalactic medium is reionised by the UV radiation from the first generation of stars and galaxies. One tracer of the process is the 21 cm line of hydrogen that will be observed by the Square Kilometre Array (SKA) at low frequencies, thus imaging the distribution of ionised and neutral regions and their evolution. To prepare for these upcoming observations, we investigate a deep learning method to predict from 21 cm maps the reionisation time field (treion(r)), i.e. the time at which each location has been reionised. treion(r) encodes the propagation of ionisation fronts in a single field, gives access to times of local reionisation or to the extent of the radiative reach of early sources. Moreover it gives access to the time evolution of ionisation on the plane of sky, when such evolution is usually probed along the line-of-sight direction. We trained a convolutional neural network (CNN) using simulated 21 cm maps and reionisation times fields produced by the simulation code 21cmFAST . We also investigate the performance of the CNN when adding instrumental effects. Globally, we find that without instrumental effects the 21 cm maps can be used to reconstruct the associated reionisation times field in a satisfying manner: the quality of the reconstruction is dependent on the redshift at which the 21 cm observation is being made and in general it is found that small scale (<10cMpc/h) features are smoothed in the reconstructed field, while larger scale features are well recovered. When instrumental effects are included, the scale dependance of reconstruction is even further pronounced, with significant smoothing on small and intermediate scales. △ Less

Submitted 2 October, 2023; v1 submitted 2 July, 2023; originally announced July 2023.

Journal ref: A&A 679, A125 (2023)

Cosmic Variance and the Inhomogeneous UV Luminosity Function of Galaxies During Reionization

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Joseph S. W. Lewis, Dominique Aubert, Jenny G. Sorce, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Romain Teyssier, Gustavo Yepes

Abstract: When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- a… ▽ More When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from $z\simeq6$-15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same $z$. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by $z=6$ (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, `rolling-over' at $M_\text{UV}\gtrsim-17$. CoDa II's UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at $z=6$-8, including the faint end, except for the faintest data point at $z=6$, based on one galaxy at $M_\text{UV}=-12.5$. According to CoDa II, the probability of observing the latter is $\sim5\%$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. We find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2-4, primarily by underestimating the variance in halo abundance when the sample volume is small. △ Less

Submitted 8 August, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: 16 pages, 16 figures, accepted by MNRAS 07/20/23, comments welcome

Accurate modelling of the Lyman-$α$ coupling for the 21-cm signal, observability with NenuFAR and SKA

Authors: Benoit Semelin, Romain Mériot, Florent Mertens, Léon V. E. Koopmans, Dominique Aubert, Rennan Barkana, Anastasia Fialkov, Satyapan Munshi, Pierre Ocvirk

Abstract: The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluct… ▽ More The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluctuations of the spin temperature of hydrogen. To interpret future measurements, controlling the level of accuracy of the Lyman-$α$ flux modelling is mandatory. In this work, we evaluate the impact of various approximations that exist in the main fast modelling approach compared to the results of a costly full radiative transfer simulation. The fast SPINTER code, presented in this work, computes the Lyman-$α$ flux including the effect of wing scatterings for an inhomogeneous emissivity field, but assuming an otherwise homogeneous expanding universe. The LICORICE code computes the full radiative transfer in the Lyman-$α$ line without any substantial approximation. We find that the difference between homogeneous and inhomogeneous gas density and temperature is very small for the computed flux. On the contrary, neglecting the effect of gas velocities produces a significant change in the computed flux. We identify the causes (mainly Doppler shifts due to velocity gradients) and quantify the magnitude of the effect in both an idealised setup and a realistic cosmological situation. We find that the amplitude of the effect, up to a factor of $\sim 2$ on the $21$ cm signal power spectrum on some scales (depending on both other model parameters and the redshift), can be easily discriminated with an SKA-like survey and already be approached, particularly for exotic signals, by the ongoing NenuFAR Cosmic Dawn Key Science Program. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 14 pages, accepted for publication in A&A

Journal ref: A&A 672, A162 (2023)

Topology of Reionisation times: concepts, measurements and comparisons to gaussian random field predictions

Authors: Emilie Thélie, Dominique Aubert, Nicolas Gillet, Julien Hiegel, Pierre Ocvirk

Abstract: In the next decade, radio telescopes like the Square Kilometer Array (SKA) will explore the Universe at high redshift, and particularly during the Epoch of Reionisation (EoR). The first structures emerged during this epoch, and their radiations have reionised the previously cold and neutral gas of the Universe creating ionised bubbles that percolate at the end of the EoR (at a redshift of approxim… ▽ More In the next decade, radio telescopes like the Square Kilometer Array (SKA) will explore the Universe at high redshift, and particularly during the Epoch of Reionisation (EoR). The first structures emerged during this epoch, and their radiations have reionised the previously cold and neutral gas of the Universe creating ionised bubbles that percolate at the end of the EoR (at a redshift of approximately 6). SKA will produce 2D images of the distribution of the neutral gas at many redshifts, pushing us to develop tools and simulations to understand its properties. This paper aims at measuring topological statistics of the EoR in the "reionisation times" fields from both cosmological and semi-analytical simulations. This field informs us about the time of reionisation of the gas at each position, is used to probe the inhomogeneities of reionisation histories and can possibly be extracted from 21 cm maps. We also compare these measurements with analytical predictions from the gaussian random field (GRF) theory. The GRF theory allows us to compute many statistics of a field: PDFs of the field or its gradient, isocontour length, critical point distributions, and skeleton length. We compare these theoretical predictions to measurements made on reionisation time fields extracted from an EMMA and a 21cmFAST simulations at 1 a cMpc/h resolution. We also compared our results to GRFs generated from the fitted power spectra of the simulation maps. Both EMMA and 21cmFAST reionisation time fields (treion(r)) are close to be gaussian fields, in contrast with the 21 cm, density or ionisation fraction that are all proven to be non-gaussian. Only accelerating ionisation fronts at the end of the EoR seem to be a cause of small non-gaussianities in treion(r). Overall our results indicate that an analytical description of the reionisation percolation can be reasonably made within the framework of GRF theory. △ Less

Submitted 22 February, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: 22 pages, 20 figures, 3 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 672, A184 (2023)

Reionization time of the Local Group and Local-Group-like halo pairs

Authors: Jenny G. Sorce, Pierre Ocvirk, Dominique Aubert, Stefan Gottloeber, Paul R. Shapiro, Taha Dawoodbhoy, Gustavo Yepes, Kyungjin Ahn, Ilian T. Iliev, Joseph S. W. Lewis

Abstract: Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper… ▽ More Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired halos within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully-coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation and isolation criteria but in other environments, alongside isolated halos within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out given our reionization model. Within this model, halos in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ~625 Myr against 660+/-4 Myr (z~8.25 against 7.87+/-0.02) on average. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS, 11 pages, 5 figures

LyC escape from SPHINX galaxies in the Epoch of Reionization

Authors: Joakim Rosdahl, Jérémy Blaizot, Harley Katz, Taysun Kimm, Thibault Garel, Martin Haehnelt, Laura C. Keating, Sergion Martin-Alvarez, Léo Michel-Dansac, Pierre Ocvirk

Abstract: We measure escape fractions, $f_{\rm esc}$, of ionizing radiation from galaxies in the SPHINX suite of cosmological radiation-hydrodynamical simulations of reionization, resolving halos with $M_{\rm vir} \gtrapprox 7.5 \times 10^7 \ M_{\odot}$ with a minimum cell width of $\approx 10$ pc. Our new and largest $20$ co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with ha… ▽ More We measure escape fractions, $f_{\rm esc}$, of ionizing radiation from galaxies in the SPHINX suite of cosmological radiation-hydrodynamical simulations of reionization, resolving halos with $M_{\rm vir} \gtrapprox 7.5 \times 10^7 \ M_{\odot}$ with a minimum cell width of $\approx 10$ pc. Our new and largest $20$ co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times $10^{11} \ M_{\odot}$. The simulated galaxies agree well with observational constraints of the UV luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over timescales of a few Myrs, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the inter-galactic medium. Statistically, $f_{\rm esc}$ peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies ($M_{*} \approx 10^7 \ M_{\odot}$, $M_{1500} \approx -17$, $Z\lesssim 5 \times 10^{-3} \ Z_{\odot}$), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to SPHINX, galaxies too dim to be yet observed, with $M_{1500} \gtrapprox -17$, provide about $55$ percent of the photons contributing to reionization. The global averaged $f_{\rm esc}$ naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 20 figures, 21 pages plus appendices. Accepted for publication in MNRAS

DUSTiER (DUST in the Epoch of Reionization): dusty galaxies in cosmological radiation-hydrodynamical simulations of the Epoch of Reionization with RAMSES-CUDATON

Authors: Joseph S. W. Lewis, Pierre Ocvirk, Yohan Dubois, Dominique Aubert, Jonathan Chardin, Nicolas Gillet, Émilie Thélie

Abstract: In recent years, interstellar dust has become a crucial topic in the study of the high and very high redshift Universe. Evidence points to the existence of high dust masses in massive star forming galaxies already during the Epoch of Reionization, potentially affecting the escape of ionising photons into the intergalactic medium. Moreover, correctly estimating dust extinction at UV wavelengths is… ▽ More In recent years, interstellar dust has become a crucial topic in the study of the high and very high redshift Universe. Evidence points to the existence of high dust masses in massive star forming galaxies already during the Epoch of Reionization, potentially affecting the escape of ionising photons into the intergalactic medium. Moreover, correctly estimating dust extinction at UV wavelengths is essential for precise ultra-violet luminosity function (UVLF) prediction and interpretation. In this paper, we investigate the impact of dust on the observed properties of high redshift galaxies, and cosmic reionization. To this end, we couple a physical model for dust production to the fully coupled radiation-hydrodynamics cosmological simulation code RAMSES-CUDATON, and perform a $16^3$, $2048^3$, simulation, that we call DUSTiER for DUST in the Epoch of Reionization. It yields galaxies with dust masses and UV slopes compatible with constraints at z $\geq 5$. We find that extinction has a dramatic impact on the bright end of the UVLF, even as early as $\rm z=8$, and our dusty UVLFs are in better agreement with observations than dust-less UVLFs. The fraction of obscured star formation rises up to 45% at $\rm z=5$, in agreement with some of the latest results from ALMA. Finally, we find that dust reduces the escape of ionising photons from galaxies more massive than $10^{10} M_\odot$ (brighter than $\approx -18$ MAB1600) by >10%, and possibly up to 80-90% for our most massive galaxies. Nevertheless, we find that the ionising escape fraction is first and foremost set by neutral Hydrogen in galaxies, as the latter produces transmissions up to 100 times smaller than through dust alone. △ Less

Submitted 9 January, 2023; v1 submitted 8 April, 2022; originally announced April 2022.

Comments: First submitted to MNRAS on 03/11/21; accepted on 30/12/22 after two revisions

Scattering of Lyα Photons through the Reionizing Intergalactic Medium: I. Spectral Energy Distribution

Authors: Hyunbae Park, Hyo Jeong Kim, Kyungjin Ahn, Hyunmi Song, Intae Jung, Pierre Ocvirk, Paul R. Shapiro, Taha Dawoodbhoy, Jenny G. Sorce, Ilian T. Iliev

Abstract: During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previ… ▽ More During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previous literature. Then, we run our code on the simulated IGM of the CoDaII simulation. We find that the scattered light can leave an observable imprint on the emergent spectrum if collected over several square arcminutes. Scattering can redden the emission by increasing the path lengths of photons, but it can also make the photons bluer by upscattering them according to the peculiar motion of the scatterer. The photons emitted on the far blue side of the resonance appear more extended in both frequency and space compared to those emitted near the resonance. This provides a discriminating feature for the blueward emission, which cannot be constrained from the unscattered light coming directly from the source. The ionization state of the IGM also affects the scattered light spectrum. When the source is in a small HII region, the emission goes through more scatterings in the surrounding HI region regardless of the initial frequency and ends up more redshifted and spatially extended. This can result in a weakening of the scattered light toward high $z$ during reionization. Our results provide a framework for interpreting the scattered light to be measured by high-$z$ integral-field-unit surveys. △ Less

Submitted 16 May, 2022; v1 submitted 13 February, 2022; originally announced February 2022.

Comments: Accepted for publication in the ApJ

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Authors: Joseph S. W. Lewis, Pierre Ocvirk, Jenny G. Sorce, Yohan Dubois, Dominique Aubert, Luke Conaboy, Paul R. Shapiro, Taha Dawoodbhoy, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Yann Rasera, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Émilie Thélie

Abstract: Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) a… ▽ More Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$). △ Less

Submitted 22 August, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: Submitted to MNRAS Accepted 2022 August 08. Received 2022 August 05; in original form 2022 February 11

A first look at the topology of reionization redshifts in models of the Epoch of the Reionization

Authors: Emilie Thélie, Dominique Aubert, Nicolas Gillet, Pierre Ocvirk

Abstract: During the EoR, the first stars and galaxies appear while creating ionized bubbles that will eventually percolate near z=6. These ionized bubbles and percolation process are nowadays under a lot of scrutiny since observations of the HI gas will be carried on in the next decade with e.g. the SKA radiotelescope. Studies of the EoR are performed on semi-analytical and fully numerical cosmological sim… ▽ More During the EoR, the first stars and galaxies appear while creating ionized bubbles that will eventually percolate near z=6. These ionized bubbles and percolation process are nowadays under a lot of scrutiny since observations of the HI gas will be carried on in the next decade with e.g. the SKA radiotelescope. Studies of the EoR are performed on semi-analytical and fully numerical cosmological simulations to investigate e.g. the topology of the process. We analyse the topology of EoR models through regions that are under the radiative influence of ionization sources. They are associated with peak patches of reionization redshift (zreion) field, for which we measure their general properties (e.g. number, size, shape, orientation). We aim at gaining insights on the geometry of the reionization process and how it relates to the matter distribution for example. We also assess how such measurements can be used to quantify the influence of physical parameters on the reionization models or the differences between fully numerical simulations and semi-analytical models. We use the DisPerSE framework (which applies the Morse theory and the persistent homology) on different EoR scenarios through gas density and zreion maps. We find that we can distinguish between EoR models with different sources using simple analyses on the number, shape and size distributions of the peak patches. For every model, we statistically show that they are rather prolate and aligned with the gas filaments. We briefly highlight that the percolation process can be followed studying zreion fields with different persistence thresholds. We show that fully numerical EMMA simulations can be made consistent with 21cmFAST models in this topological framework as long as the source distribution is diffuse enough. △ Less

Submitted 23 November, 2021; originally announced November 2021.

Comments: 13 pages, 20 figures, 1 table. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 658, A139 (2022)

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

Authors: Harley Katz, Joakim Rosdahl, Taysun Kimm, Thibault Garel, Jérémy Blaizot, Martin G. Haehnelt, Léo Michel-Dansac, Sergio Martin-Alvarez, Julien Devriendt, Adrianne Slyz, Romain Teyssier, Pierre Ocvirk, Nicolas Laporte, Richard Ellis

Abstract: ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenom… ▽ More ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63μm}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged) △ Less

Submitted 5 January, 2022; v1 submitted 2 August, 2021; originally announced August 2021.

Comments: 20 pages, 18 figures, Accepted by MNRAS

Best Practices for Data Publication in the Astronomical Literature

Authors: Tracy X. Chen, Marion Schmitz, Joseph M. Mazzarella, Xiuqin Wu, Julian C. van Eyken, Alberto Accomazzi, Rachel L. Akeson, Mark Allen, Rachael Beaton, G. Bruce Berriman, Andrew W. Boyle, Marianne Brouty, Ben Chan, Jessie L. Christiansen, David R. Ciardi, David Cook, Raffaele D'Abrusco, Rick Ebert, Cren Frayer, Benjamin J. Fulton, Christopher Gelino, George Helou, Calen B. Henderson, Justin Howell, Joyce Kim , et al. (20 additional authors not shown)

Abstract: We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to… ▽ More We present an overview of best practices for publishing data in astronomy and astrophysics journals. These recommendations are intended as a reference for authors to help prepare and publish data in a way that will better represent and support science results, enable better data sharing, improve reproducibility, and enhance the reusability of data. Observance of these guidelines will also help to streamline the extraction, preservation, integration and cross-linking of valuable data from astrophysics literature into major astronomical databases, and consequently facilitate new modes of science discovery that will better exploit the vast quantities of panchromatic and multi-dimensional data associated with the literature. We encourage authors, journal editors, referees, and publishers to implement the best practices reviewed here, as well as related recommendations from international astronomical organizations such as the International Astronomical Union (IAU) for publication of nomenclature, data, and metadata. A convenient Checklist of Recommendations for Publishing Data in the Literature is included for authors to consult before the submission of the final version of their journal articles and associated data files. We recommend that publishers of journals in astronomy and astrophysics incorporate a link to this document in their Instructions to Authors. △ Less

Submitted 16 April, 2022; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 19 pages, 1 figure, 3 tables, accepted for publication in ApJS

Crucial Factors for Lyα Transmission in the Reionizing Intergalactic Medium: Infall Motion, HII Bubble Size, and Self-shielded Systems

Authors: Hyunbae Park, Intae Jung, Hyunmi Song, Pierre Ocvirk, Paul R. Shapiro, Taha Dawoodbhoy, Ilian T. Iliev, Kyungjin Ahn, Michele Bianco, Hyo Jeong Kim

Abstract: Using the CoDa II simulation, we study the Ly$α$ transmissivity of the intergalactic medium (IGM) during reionization. At $z>6$, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Ly$α$, suppressing the tran… ▽ More Using the CoDa II simulation, we study the Ly$α$ transmissivity of the intergalactic medium (IGM) during reionization. At $z>6$, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Ly$α$, suppressing the transmission up to roughly the circular velocity of the galaxy. In some sight lines, an optically thin blob generated by a supernova in a neighboring galaxy results in a peak feature, which can be mistaken for a blue peak. Redward of the resonance absorption, the damping-wing opacity correlates with the global IGM neutral fraction and the UV magnitude of the source galaxy. Brighter galaxies tend to suffer lower opacity because they tend to reside in larger HII regions, and the surrounding IGM transmits redder photons, which are less susceptible to attenuation, owing to stronger infall velocity. The HII regions are highly nonspherical, causing both sight-line-to-sight-line and galaxy-to-galaxy variation in opacity. Also, self-shielded systems within HII regions strongly attenuate the emission for certain sight lines. All these factors add to the transmissivity variation, requiring a large sample size to constrain the average transmission. The variation is largest for fainter galaxies at higher redshift. The 68\% range of the transmissivity is similar to or greater than the median for galaxies with $M_{\rm UV}\ge-21$ at $z\ge7$, implying that more than a hundred galaxies would be needed to measure the transmission to 10\% accuracy. △ Less

Submitted 6 December, 2021; v1 submitted 22 May, 2021; originally announced May 2021.

Comments: 18 pages, 16 figures, 4 tables, Published in ApJ

Lyman-alpha opacities at z=4-6 require low mass, radiatively-suppressed galaxies to drive cosmic reionization

Authors: Pierre Ocvirk, Joseph S. W. Lewis, Nicolas Gillet, Jonathan Chardin, Dominique Aubert, Nicolas Deparis, Emilie Thelie

Abstract: The high redshift Lyman-alpha forest, in particular the Gunn-Peterson trough, is the most unambiguous signature of the neutral to ionized transition of the intergalactic medium (IGM) taking place during the Epoch of Reionization (EoR). Recent studies, e.g. Kulkarni et al. (2019a) and Keating et al. (2019), showed that reproducing the observed Lyman-alpha opacities after overlap required a non-mono… ▽ More The high redshift Lyman-alpha forest, in particular the Gunn-Peterson trough, is the most unambiguous signature of the neutral to ionized transition of the intergalactic medium (IGM) taking place during the Epoch of Reionization (EoR). Recent studies, e.g. Kulkarni et al. (2019a) and Keating et al. (2019), showed that reproducing the observed Lyman-alpha opacities after overlap required a non-monotonous evolution of cosmic emissivity: rising, peaking at z=6, and then decreasing onwards to z=4. Such an evolution is puzzling considering galaxy build-up and the cosmic star formation rate are still continously on the rise at these epochs. Here, we use new RAMSES-CUDATON simulations to show that such a peaked evolution may occur naturally in a fully coupled radiation-hydrodynamical framework. In our fiducial run, cosmic emissivity at z>6 is dominated by a low mass (M$_{\rm DM}<2.10^9$ M$_{\odot}$), high escape fraction halo population, driving reionization, up to overlap. Approaching z=6, this population is radiatively suppressed due to the rising ionizing UV background, and its emissivity drops. In the meantime, the high mass halo population builds up and its emissivity rises, but not fast enough to compensate the dimming of the low mass haloes, because of low escape fractions. The combined ionizing emissivity of these two populations therefore naturally results in a rise and fall of the cosmic emissivity, from z=12 to z=4, with a peak at z=6. An alternative run, which features higher escape fractions for the high mass haloes and later suppression at low mass, leads to overshooting the ionizing rate, over-ionizing the IGM and therefore too low Lyman-alpha opacities. △ Less

Submitted 1 September, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

21cm signal predictions at Cosmic Dawn and Reionization with coupled radiative-hydrodynamics

Authors: Nicolas J. F. Gillet, Dominique Aubert, Florent G. Mertens, Pierre Ocvirk

Abstract: The process of heating and reionization of the Universe at high redshift links small scale structures/galaxy formation and large scale inter-galactic medium properties. Even if the first one is difficult to observe, an observation window is opening on the second one, with the promising development of current and future radio telescopes. They will permit to observe the 21cm brightness temperature g… ▽ More The process of heating and reionization of the Universe at high redshift links small scale structures/galaxy formation and large scale inter-galactic medium properties. Even if the first one is difficult to observe, an observation window is opening on the second one, with the promising development of current and future radio telescopes. They will permit to observe the 21cm brightness temperature global signal and fluctuations. The need of large scale simulations is therefore strong to understand the properties of the IGM that will be observed. But at the same time the urge to resolve the structures responsible of those process is important. We introduce in this study, a coupled hydro-radiative transfer simulations of the Cosmic Dawn and Reionization with a simple sub-grid star formation process developed and calibrated on the state of the art simulation CoDaII. This scheme permits to follow consistently dark matter, hydrodynamics and radiative transfer evolution's on large scales, while the sub-grid models bridges to the galaxy formation scale. We process the simulation to produce 21cm signal as close as possible to the observations. △ Less

Submitted 4 March, 2021; originally announced March 2021.

Cosmological Magnetogenesis: The Biermann Battery during the Epoch of Reionization

Authors: Omar Attia, Romain Teyssier, Harley Katz, Taysun Kimm, Sergio Martin-Alvarez, Pierre Ocvirk, Joakim Rosdahl

Abstract: We investigate the effect of the Biermann battery during the Epoch of Reionization (EoR) using cosmological Adaptive Mesh Refinement simulations within the framework of the SPHINX project. We develop a novel numerical technique to solve for the Biermann battery term in the Constrained Transport method, preserving both the zero divergence of the magnetic field and the absence of Biermann battery fo… ▽ More We investigate the effect of the Biermann battery during the Epoch of Reionization (EoR) using cosmological Adaptive Mesh Refinement simulations within the framework of the SPHINX project. We develop a novel numerical technique to solve for the Biermann battery term in the Constrained Transport method, preserving both the zero divergence of the magnetic field and the absence of Biermann battery for isothermal flows. The structure-preserving nature of our numerical method turns out to be very important to minimise numerical errors during validation tests of the propagation of a Strömgren sphere and of a Sedov blast wave. We then use this new method to model the evolution of a 2.5 and 5 co-moving Mpc cosmological box with a state-of-the-art galaxy formation model within the RAMSES code. Contrary to previous findings, we show that three different Biermann battery channels emerge: the first one is associated with linear perturbations before the EoR, the second one is the classical Biermann battery associated with reionization fronts during the EoR, and the third one is associated with strong, supernova-driven outflows. While the two former channels generate spontaneously volume-filling magnetic fields with a strength on the order or below $10^{-20}$ G, the latter, owing to the higher plasma temperature and a marginally-resolved turbulent dynamo, reaches a field strength as high as $10^{-18}$ G in the intergalactic medium around massive haloes. △ Less

Submitted 19 May, 2021; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: 14 pages, 18 figures, accepted in MNRAS

Lyman-alpha transmission properties of the intergalactic medium in the CoDaII simulation

Authors: Max Gronke, Pierre Ocvirk, Charlotte Mason, Jorryt Matthee, Sarah E. I. Bosman, Jenny G. Sorce, Joseph Lewis, Kyungjin Ahn, Dominique Aubert, Taha Dawoodbhoy, Ilian T. Iliev, Paul R. Shapiro, Gustavo Yepes

Abstract: The decline in abundance of Lyman-$α$ (Ly$α$) 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$α$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$.… ▽ More The decline in abundance of Lyman-$α$ (Ly$α$) 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$α$ 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$α$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$α$ 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 1\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$α$ peak. We discuss our results in the context of the interpretation of observations. △ Less

Submitted 21 October, 2021; v1 submitted 29 April, 2020; originally announced April 2020.

Comments: 13 pages, 12 figures. Published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 3, December 2021, Pages 3697-3709

Galactic ionising photon budget during the Epoch of Reionisation in the Cosmic Dawn II simulation

Authors: Joseph S. W. Lewis, Pierre Ocvirk, Dominique Aubert, Jenny G. Sorce, Paul R. Shapiro, Nicolas Deparis, Taha Dawoodbhoy, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Kyungjin Ahn, Ilian T. Iliev, Jonathan Chardin

Abstract: 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 transmissio… ▽ More 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. △ Less

Submitted 18 June, 2020; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: Submitted to MNRAS on 09/01/2020. Accepted on 12/06/2020

A deep learning model to emulate simulations of cosmic reionization

Authors: Jonathan Chardin, Grégoire Uhlrich, Dominique Aubert, Nicolas Deparis, Nicolas Gillet, Pierre Ocvirk, Joseph Lewis

Abstract: 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}}$ a… ▽ More 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. △ Less

Submitted 18 October, 2019; v1 submitted 16 May, 2019; originally announced May 2019.

Comments: 11 pages, 8 figures, accepted in MNRAS (official release 21st of october)

Cosmic Dawn II (CoDa II): a new radiation-hydrodynamics simulation of the self-consistent coupling of galaxy formation and reionization

Authors: Pierre Ocvirk, Dominique Aubert, Jenny G. Sorce, Paul R. Shapiro, Nicolas Deparis, Taha Dawoodbhoy, Joseph Lewis, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Kyungjin Ahn, Ilian T. Iliev, Yehuda Hoffman

Abstract: 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 RA… ▽ More 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 $Λ$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. △ Less

Submitted 14 October, 2020; v1 submitted 27 November, 2018; originally announced November 2018.

Comments: 20 pages, 14 figures, accepted for publication in MNRAS

Suppression of Star Formation in Low-Mass Galaxies Caused by the Reionization of their Local Neighborhood

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Dominique Aubert, Nicolas Gillet, Jun-Hwan Choi, Ilian T. Iliev, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Anson D'Aloisio, Hyunbae Park, Yehuda Hoffman

Abstract: 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 s… ▽ More 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. △ Less

Submitted 25 July, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 15 pages, 15 figures, accepted by MNRAS 07/17/18, comments welcome

Journal ref: MNRAS, 480, 1740-1753 (2018)

Impact of the reduced speed of light approximation on the post-overlap neutral hydrogen fraction in numerical simulations of the epoch of reionization

Authors: P. Ocvirk, D. Aubert, J. Chardin, N. Deparis, J. Lewis

Abstract: 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 som… ▽ More 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. △ Less

Submitted 3 May, 2019; v1 submitted 6 March, 2018; originally announced March 2018.

Comments: Includes languages corrections for UK english conventions. Accepted for publication in A&A, 9 pages, 9 figures

Journal ref: A&A 626, A77 (2019)

Impact of the reduced speed of light approximation on ionization front velocities in cosmological simulations of the epoch of reionization

Authors: Nicolas Deparis, Dominique Aubert, Pierre Ocvirk, Jonathan Chardin, Joseph Lewis

Abstract: 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 g… ▽ More 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. △ Less

Submitted 8 January, 2019; v1 submitted 5 March, 2018; originally announced March 2018.

Comments: 11 pages, 13 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 622, A142 (2019)

The Inhomogeneous Reionization Times of Present-day Galaxies

Authors: Dominique Aubert, Nicolas Deparis, Pierre Ocvirk, Paul R. Shapiro, Ilian T. Iliev, Gustavo Yepes, Stefan Gottloeber, Yehuda Hoffman, Romain Teyssier

Abstract: 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 t… ▽ More 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). △ Less

Submitted 20 February, 2018; v1 submitted 5 February, 2018; originally announced February 2018.

Comments: 8 pages, 5 figures, updated version accepted for publication in ApJL

The SPHINX Cosmological Simulations of the First Billion Years: the Impact of Binary Stars on Reionization

Authors: Joakim Rosdahl, Harley Katz, Jeremy Blaizot, Taysun Kimm, Leo Michel-Dansac, Thibault Garel, Martin Haehnelt, Pierre Ocvirk, Romain Teyssier

Abstract: 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… ▽ More 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. △ Less

Submitted 25 June, 2018; v1 submitted 22 January, 2018; originally announced January 2018.

Comments: 18 pages, 15 figures. Accepted for publication in MNRAS

French SKA White Book - The French Community towards the Square Kilometre Array

Authors: F. Acero, J. -T. Acquaviva, R. Adam, N. Aghanim, M. Allen, M. Alves, R. Ammanouil, R. Ansari, A. Araudo, E. Armengaud, B. Ascaso, E. Athanassoula, D. Aubert, S. Babak, A. Bacmann, A. Banday, K. Barriere, F. Bellossi, J. -P. Bernard, M. G. Bernardini, M. Béthermin, E. Blanc, L. Blanchet, J. Bobin, S. Boissier , et al. (153 additional authors not shown)

Abstract: 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 p… ▽ More 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. △ Less

Submitted 28 March, 2018; v1 submitted 19 December, 2017; originally announced December 2017.

Comments: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferrière, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Deschênes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

Gaia Data Release 1. Testing the parallaxes with local Cepheids and RR Lyrae stars

Authors: Gaia Collaboration, G. Clementini, L. Eyer, V. Ripepi, M. Marconi, T. Muraveva, A. Garofalo, L. M. Sarro, M. Palmer, X. Luri, R. Molinaro, L. Rimoldini, L. Szabados, I. Musella, R. I. Anderson, T. Prusti, J. H. J. de Bruijne, A. G. A. Brown, A. Vallenari, C. Babusiaux, C. A. L. Bailer-Jones, U. Bastian, M. Biermann, D. W. Evans, F. Jansen , et al. (566 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 1 May, 2017; originally announced May 2017.

Comments: 29 pages, 25 figures. Accepted for publication by A&A

Journal ref: A&A 605, A79 (2017)

Gaia Data Release 1. Open cluster astrometry: performance, limitations, and future prospects

Authors: Gaia Collaboration, F. van Leeuwen, A. Vallenari, C. Jordi, L. Lindegren, U. Bastian, T. Prusti, J. H. J. de Bruijne, A. G. A. Brown, C. Babusiaux, C. A. L. Bailer-Jones, M. Biermann, D. W. Evans, L. Eyer, F. Jansen, S. A. Klioner, U. Lammers, X. Luri, F. Mignard, C. Panem, D. Pourbaix, S. Randich, P. Sartoretti, H. I. Siddiqui, C. Soubiran , et al. (567 additional authors not shown)

Abstract: 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 ast… ▽ More 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. △ Less

Submitted 3 March, 2017; originally announced March 2017.

Comments: Accepted for publication by A&A. 21 pages main text plus 46 pages appendices. 34 figures main text, 38 figures appendices. 8 table in main text, 19 tables in appendices

Journal ref: A&A 601, A19 (2017)

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Authors: Pierre Ocvirk, Nicolas Gillet, Paul R. Shapiro, Dominique Aubert, Ilian T. Iliev, Romain Teyssier, Gustavo Yepes, Jun-Hwan Choi, David Sullivan, Alexander Knebe, Stefan Gottloeber, Anson D'Aloisio, Hyunbae Park, Yehuda Hoffman, Timothy Stranex

Abstract: 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 reso… ▽ More 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) △ Less

Submitted 6 September, 2016; v1 submitted 30 October, 2015; originally announced November 2015.

Comments: 26 pages, 16 figures, accepted in MNRAS

EMMA: an AMR cosmological simulation code with radiative transfer

Authors: Dominique Aubert, Nicolas Deparis, Pierre Ocvirk

Abstract: 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.… ▽ More 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. △ Less

Submitted 17 August, 2015; originally announced August 2015.

Comments: Accepted for publication in MNRAS, 28 pages, 27 figures

Vast planes of satellites in a high resolution simulation of the Local Group: comparison to Andromeda

Authors: N. Gillet, P. Ocvirk, D. Aubert, A. Knebe, N. Libeskind, G. Yepes, S. Gottlöber, Y. Hoffman

Abstract: 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, rou… ▽ More 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. △ Less

Submitted 5 January, 2015; v1 submitted 9 December, 2014; originally announced December 2014.

The reionization of galactic satellite populations

Authors: P. Ocvirk, N. Gillet, D. Aubert, A. Knebe, N. Libeskind, J. Chardin, S. Gottlöber, G. Yepes, Y. Hoffman

Abstract: 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… ▽ More 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. △ Less

Submitted 29 September, 2014; originally announced September 2014.

Comments: 19 pages, 6 figures, accepted in ApJ

Local reionizations histories with merger tree of HII regions

Authors: Jonathan Chardin, Dominique Aubert, Pierre Ocvirk

Abstract: 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.… ▽ More 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. △ Less

Submitted 23 June, 2014; originally announced June 2014.

Comments: 16 pages, 9 figures, Accepted for publication in A&A

MSC Class: 83F05

Journal ref: A&A 568, A52 (2014)

High resolution simulations of the reionization of an isolated Milky Way - M31 galaxy pair

Authors: P. Ocvirk, D. Aubert, J. Chardin, A. Knebe, N. Libeskind, S. Gottlöber, G. Yepes, Y. Hoffman

Abstract: 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 o… ▽ More 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. △ Less

Submitted 30 September, 2013; originally announced October 2013.

Comments: Accepted for publication in ApJ. 14 pages, 4 figures, 1 table

A study of simulated histories of reionization with merger trees of HII regions

Authors: Jonathan Chardin, Dominique Aubert, Pierre Ocvirk

Abstract: 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… ▽ More 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. △ Less

Submitted 2 October, 2012; originally announced October 2012.

Comments: Accepted for publication in A&A

MSC Class: 85A25; 85A40 ACM Class: J.2.3; J.2.9

Constraining local UV field geometry at reionization using Milky Way satellites

Authors: P. Ocvirk, D. Aubert

Abstract: 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 propagat… ▽ More 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. △ Less

Submitted 7 November, 2011; originally announced November 2011.

Comments: Proceedings of the conference "Assembling the puzzle of the Milky Way" held at Le Grand Bornand, 2011

Fitting galaxy spectra with STECKMAP: a user guide

Authors: P. Ocvirk

Abstract: 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… ▽ More 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 https://github.com/pocvirk/STECKMAP. 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. △ Less

Submitted 6 March, 2023; v1 submitted 23 August, 2011; originally announced August 2011.

Comments: 27 pages, updated with github repository url

A signature of the internal reionisation of the Milky Way?

Authors: Pierre Ocvirk, Dominique Aubert

Abstract: 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… ▽ More 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. △ Less

Submitted 4 August, 2011; originally announced August 2011.

Comments: Accepted for publication in Letters of the Monthly Notices of the Royal Astronomical Society. 6 pages, 3 figures

Star formation history of barred disc galaxies

Authors: Patricia Sanchez-Blazquez, Pierre Ocvirk, Brad K. Gibson, Isabel Perez, Reynier F. Peletier

Abstract: 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 d… ▽ More 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) △ Less

Submitted 19 March, 2011; originally announced March 2011.

Comments: 24 pages, 23 figures, accepted for publication in MNRAS

Pinning down the ram-pressure-induced halt of star formation in the Virgo cluster spiral galaxy NGC 4388. A joint inversion of spectroscopic and photometric data

Authors: C. Pappalardo, A. Lancon, B. Vollmer, P. Ocvirk, S. Boissier, A. Boselli

Abstract: 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 m… ▽ More 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. △ Less

Submitted 18 February, 2010; originally announced February 2010.

Comments: 12 pages, 10 figures, accepted for publication in A&A

Fake star formation bursts: blue horizontal branch stars masquerade as young stars in optical integrated light spectroscopy

Authors: P. Ocvirk

Abstract: 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 Schi… ▽ More 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. △ Less

Submitted 16 November, 2009; originally announced November 2009.

Comments: 14 pages, 6 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.709:88-96,2010

The impact of binary star yields on the spectra of galaxies

Authors: A. E. Sansom, R. G. Izzard, P. Ocvirk

Abstract: 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 bin… ▽ More 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. △ Less

Submitted 2 July, 2009; originally announced July 2009.

Comments: 15 pages, 5 figures, 4 tables, accepted for publication in MNRAS

Cosmic star-formation history from a non-parametric inversion of infrared galaxy counts

Authors: Damien Le Borgne, David Elbaz, Pierre Ocvirk, Christophe Pichon

Abstract: [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… ▽ More [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. △ Less

Submitted 2 July, 2009; v1 submitted 23 January, 2009; originally announced January 2009.

Comments: 20 pages, 15 figures, 1 table. Accepted in A&A. The model is available at http://www.iap.fr/users/leborgne/irmodel09/index.html

Bimodal gas accretion in the Horizon-MareNostrum galaxy formation simulation

Authors: P. Ocvirk, C. Pichon, R. Teyssier

Abstract: (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 hi… ▽ More (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. △ Less

Submitted 7 August, 2008; v1 submitted 31 March, 2008; originally announced March 2008.

Comments: 15 pages, 6 figures, accepted for publication in MNRAS, 29th of july 2008

Spectroscopic ages and metallicities of stellar populations: validation of full spectrum fitting

Authors: M. Koleva, Ph. Prugniel, P. Ocvirk, D. Le Borgne, C. Soubiran

Abstract: 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 wit… ▽ More 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. △ Less

Submitted 6 January, 2008; originally announced January 2008.

Comments: accepted in MNRAS; Full version available at http://www-obs.univ-lyon1.fr/labo/perso/prugniel/mina/koleva.pdf

Analysis of stellar populations with large empirical libraries at high spectral resolution

Authors: Philippe Prugniel, Mina Koleva, Pierre Ocvirk, Damien Le Borgne, Caroline Soubiran

Abstract: 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 libra… ▽ More 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 [$α$/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 ($α$-enhanced) globular clusters. △ Less

Submitted 6 March, 2007; originally announced March 2007.

Comments: 4 pages; proceedings of IAU Symposium No. 241, "Stellar Populations as Building Blocks of Galaxies", editors A. Vazdekis and R. Peletier

Comparison of different spectral population models

Authors: Mina Koleva, Philippe Prugniel, Pierre Ocvirk, Damien Le Borgne, Igor Chilingarian, Caroline Soubiran

Abstract: 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 i… ▽ More 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). △ Less

Submitted 6 March, 2007; originally announced March 2007.

Comments: 2 pages; proceedings of IAU Symposium No. 241, "Stellar Populations as Building Blocks of Galaxies", editors A. Vazdekis and R. Peletier