The effects of galactic fountains on the chemical evolution of galaxies

2008/2009In this thesis we study the effect of galactic fountains, namely gas and flows from the disk of galaxies produced by multiple supernova explosions, on the chemical evolution of galaxies. Sequential supernova explosions create a superbubble, whereas the swept up interstellar medium is concentrated in a supershell which can break out a stratified medium, producing bipolar outflows. The gas of the supershells can fragment into clouds which eventually fall toward the disk producing so-called galactic fountains. Many works in literature have dealt with superbubble expansion in stratified media. However, very few papers in the past have taken into account the chemical evolution of the superbubble and how the supershell get polluted from the metals produced by supernova explosions. With this thesis for the first time the effect of galactic fountains we consider in a detailed chemical evolution model for the Milky Way. In the first part of our work we study the expansion law and chemical enrichment of a supershell powered by the energetic feedback of a typical Galactic OB association at various galactocentric radii. We follow the orbits of the fragments created when the supershell breaks out and we compare their kinetic and chemical properties with the available observations of high - and intermediate - velocity clouds. We use the Kompaneets (1960) approximation for the evolution of the superbubble driven by sequential supernova explosions and we compute the abundances of oxygen and iron residing in the thin cold supershell. Due to Rayleigh-Taylor instabilities we assume that supershells are fragmented and we follow the orbit of the clouds either ballistically or by means of a hybrid model considering viscous interaction between the clouds and the extra-planar gas. We find that if the initial metallicity is solar, the pollution from the dying stars of the OB association has a negligible effect on the chemical composition of the clouds. The maximum height reached by the clouds above the plane seldom exceeds 2 kpc and when averaging over different throwing angles, the landing coordinate differs from the throwing coordinate by only 1 kpc. Therefore, it is unlikely that galactic fountains can affect abundance gradients on large scales. The range of heights and [O/Fe] ratios spanned by our clouds suggest that the high velocity clouds cannot have a Galactic origin, whereas intermediate velocity clouds have kinematic properties similar to our predicted clouds but have observed overabundances of the [O/Fe] ratios that can be reproduced only with initial metallicities which are too low compared to those of the Galaxy disk. Even if it is unlikely that galactic fountains can affect abundance gradients on large scales, they can still affect the chemical enrichment of the interstellar medium (ISM) because of the time-delay due to the non-negligible time taken by fountains to orbit around and fall back into the Galaxy. This implies a delay in the mixing of metals in ISM which conflicts with the instantaneous mixing approximation usually assumed in all models in literature. We test whether relaxing this approximation in a detailed chemical evolution model can improve or worsen the agreement with observations. To do that, we investigate two possible causes for relaxing of the instantaneous mixing: i) the ``galactic fountain time delay effect'' and ii) the ``metal cooling time delay effect''. We find that the effect of galactic fountains is negligible if an average time delay of 0.1 Gyr, as suggested by our model, is assumed. Longer time delays produce differences in the results but they are not realistic. The metal cooling time delays produce strong effects on the evolution of the chemical abundances only if we adopt stellar yields depending on metallicity. If, instead, the yields computed for the solar chemical composition are adopted, negligible effects are produced, as in the case of the galactic fountain delay. The relaxation of the IMA by means of the galactic fountain model, where the delay is considered only for massive stars and only in the disk, does not affect the chemical evolution results. The combination of metal dependent yields and time delay in the chemical enrichment from all stars starting from the halo phase, instead, produces results at variance with observations.XXII Ciclo197

The Fall of a Giant. Chemical evolution of Enceladus, alias the Gaia Sausage

We present the first chemical evolution model for Enceladus, alias the Gaia Sausage, to investigate the star formation history of one of the most massive satellites accreted by the Milky Way during a major merger event. Our best chemical evolution model for Enceladus nicely fits the observed stellar [$\alpha$/Fe]-[Fe/H] chemical abundance trends, and reproduces the observed stellar metallicity distribution function, by assuming low star formation efficiency, fast infall time scale, and mild outflow intensity. We predict a median age for Enceladus stars $12.33^{+0.92}_{-1.36}$ Gyr, and - at the time of the merger with our Galaxy ($\approx10$ Gyr ago from Helmi et al.) - we predict for Enceladus a total stellar mass $M_{\star} \approx 5 \times 10^{9}\,\text{M}_{\odot}$. By looking at the predictions of our best model, we discuss that merger events between the Galaxy and systems like Enceladus may have inhibited the gas accretion onto the Galaxy disc at high redshifts, heating up the gas in the halo. This scenario could explain the extended period of quenching in the star formation activity of our Galaxy about 10 Gyr ago, which is predicted by Milky Way chemical evolution models, in order to reproduce the observed bimodality in [$\alpha$/Fe]-[Fe/H] between thick- and thin-disc stars.Comment: Accepted for publication in MNRAS Letter

ISM simulations: an overview of models

Until recently the dynamical evolution of the interstellar medium (ISM) was simu- lated using collisional ionization equilibrium (CIE) conditions. However, the ISM is a dynamical system, in which the plasma is naturally driven out of equilibrium due to atomic and dynamic processes operating on different timescales. A step forward in the field comprises a multi-fluid approach taking into account the joint thermal and dynamical evolutions of the ISM gas

Chemical clocks and their time zones: understanding the [s/Mg]--age relation with birth radii

The relative enrichment of s-process to $\alpha$-elements ([s/$\alpha$]) has been linked with age, providing a potentially useful avenue in exploring the Milky Way's chemical evolution. However, the age--[s/$\alpha$] relationship is non-universal, with dependencies on metallicity and current location in the Galaxy. In this work, we examine these chemical clock tracers across birth radii ($\rm \text{R}_\text{birth}$), recovering the inherent trends between the variables. We derive $\rm \text{R}_\text{birth}$ and explore the [s/$\alpha$]--age--$\rm \text{R}_\text{birth}$ relationship for 36,652 APOGEE DR17 red giant and 24,467 GALAH DR3 main sequence turnoff and subgiant branch disk stars using [Ce/Mg], [Ba/Mg], and [Y/Mg]. We discover that the age--[s/Mg] relation is strongly dependent on birth location in the Milky Way, with stars born in the inner disk having the weakest correlation. This is congruent with the Galaxy's initially weak, negative [s/Mg] radial gradient, which becomes positive and steep with time. We show that the non-universal relations of chemical clocks is caused by their fundamental trends with $\rm \text{R}_\text{birth}$ over time, and suggest that the tight age--[s/Mg] relation obtained with solar-like stars is due to similar $\rm \text{R}_\text{birth}$ for a given age. Our results are put into context with a Galactic chemical evolution model, where we demonstrate the need for data-driven nucleosynthetic yields.Comment: submitted to MNRA

Origin of neutron capture elements with the Gaia-ESO survey: the evolution of s- and r-process elements across the Milky Way

We study the abundance patterns and the radial gradients of s-process elements (Y, Zr, Ba, La and Ce), r-process elements (Eu) and mixed-process elements (Mo, Nd and Pr) in the Galactic thin disc by means of a detailed two-infall chemical evolution model for the Milky Way with state-of-the-art nucleosynthesis prescriptions. We consider r-process nucleosynthesis from merging neutron stars (MNS), magneto-rotational supernovae (MR-SNe) and s-process synthesis from low- and intermediate- mass stars (LIMS) and rotating massive stars. The predictions of our model are compared with data from the sixth data release of the Gaia-ESO survey, from which we consider 62 open clusters with age > 0.1 Gyr and 1300 Milky Way disc field stars. We conclude that: i) the [Eu/Fe] vs. [Fe/H] is reproduced by both a prompt and a delayed source, but the quick source completely dominates the Eu production; ii) rotation in massive stars contribute substantially to the s-process elements of the first peak, but MNS and MR-SNe are necessary in order to reproduce the observations; iii) due to the adopted yields, our model overpredicts Pr and underpredicts Nd, while the [Mo/Fe] vs. [Fe/H] is nicely reproduced. For the radial gradients, we conclude that: i) our predicted slope of the [Fe/H] gradient is in agreement with the one observed in open clusters by Gaia-ESO and other high-resolution spectroscopic surveys. ii) The predicted slope of the [Eu/H] radial gradient is steeper than the observed one, independently on how quick the production of Eu is. We discuss the possible causes of this discrepancy in terms of both different Galaxy formation scenarios and stellar radial migration effects. iii) For all the elements belonging to the second s-process peak (Ba, La, Ce) as well as for Pr, we predict a plateau at low Galactocentric distances, which is probably due to the enhanced enrichment from LIMS in the inner regions.Comment: 17 pages, 11 figures. Submitted to MNRA

Electron Distribution in the Galactic Disk - Results From a Non-Equilibrium Ionization Model of the ISM

Using three-dimensional non-equilibrium ionization (NEI) hydrodynamical simulation of the interstellar medium (ISM), we study the electron density, $n_{e}$, in the Galactic disk and compare it with the values derived from dispersion measures towards pulsars with known distances located up to 200 pc on either side of the Galactic midplane. The simulation results, consistent with observations, can be summarized as follows: (i) the DMs in the simulated disk lie between the maximum and minimum observed values, (ii) the log derived from lines of sight crossing the simulated disk follows a Gaussian distribution centered at \mu=-1.4 with a dispersion \sigma=0.21, thus, the Galactic midplane =0.04\pm 0.01$cm$^{-3}$, (iii) the highest electron concentration by mass (up to 80%) is in the thermally unstable regime (200<T<10^{3.9} K), (iv) the volume occupation fraction of the warm ionized medium is 4.9-6%, and (v) the electrons have a clumpy distribution along the lines of sight.Comment: Letter accepted for publication in Monthly Notices of the Royal Astronomical Societ

The cerium content of the Milky Way as revealed by Gaia DR3 GSP-Spec abundances

[Abstract]: The recent Gaia third data release contains a homogeneous analysis of millions of high-quality Radial Velocity Spectrometer (RVS) stellar spectra by the GSP-Spec module. This led to the estimation of millions of individual chemical abundances and allows us to chemically map the Milky Way. The published GSP-Spec abundances include three heavy elements produced by neutron-captures in stellar interiors: Ce, Zr, and Nd. Aims. We study the Galactic content in cerium based on these Gaia/RVS data and discuss the chemical evolution of this element. Methods. We used a sample of about 30 000 local thermal equilibrium Ce abundances, selected after applying different combinations of GSP-Spec flags. Based on the Gaia DR3 astrometric data and radial velocities, we explore the cerium content in the Milky Way and, in particular, in its halo and disc components. Results. The high quality of the Ce GSP-Spec abundances is quantified through literature comparisons. We found a rather flat [Ce/Fe] versus [M/H] trend. We also found a flat radial gradient in the disc derived from field stars and, independently, from about 50 open clusters. This agrees with previous studies. The [Ce/Fe] vertical gradient was also estimated. We also report an increasing [Ce/Ca] versus [Ca/H] in the disc, illustrating the late contribution of asymptotic giant branch stars with respect to supernovae of type II. Our cerium abundances in the disc, including the young massive population, are well reproduced by a new three-infall chemical evolution model. In the halo population, the M 4 globular cluster is found to be enriched in cerium. Moreover, 11 stars with cerium abundances belonging to the Thamnos, Helmi Stream, and Gaia-Sausage-Enceladus accreted systems were identified from chemo-dynamical diagnostics. We found that the Helmi Stream might be slightly underabundant in cerium compared to the two other systems. Conclusions. This work illustrates the high quality of the GSP-Spec chemical abundances, which significantly contribute to unveiling the heavy-element evolution history of the Milky Way.We thank the referee for their valuable comments. ES received funding from the European Union’s Horizon 2020 research and innovation program under SPACE-H2020 grant agreement number 101004214 (EXPLORE project). ARB also acknowledges support from this Horizon program. PAP and EP thanks the Centre National d’Etudes Spatiales (CNES) for funding support. VG acknowledges support from the European Research Council Consolidator Grant funding scheme (project ASTEROCHRONOMETRY, G.A. n. 772293, http://www.asterochronometry.eu ). Special thanks to Niels Nieuwmunster and Botebar for grateful comments on figures. This work has made use of data from the European Space Agency (ESA) mission Gaia ( https://www.cosmos.esa.int/gaia ), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium ). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement