Nitrogen evolution in the halo, thick disc, thin disc and bulge of the Galaxy

We study the evolution of nitrogen (N) in the Galactic halo, thick disc, thin disc, and bulge by comparing detailed chemical evolution models with recent observations. The models used in this work have already been constrained to explain the abundance patterns of α-elements and the metallicity distribution functions of halo, disc, and bulge stars; here, we adopt them to investigate the origin and evolution of N in the different Galactic components. First, we consider different sets of yields and study the importance of the various channels proposed for N production. Secondly, we apply the reference models to study the evolution of both the Galactic discs and bulge. We conclude that: i) primary N produced by rotating massive stars is required to reproduce the plateau in log(N/O) and [N/Fe] ratios at low metallicity, as well as the secondary and primary production from low- and intermediate-mass stars to reproduce the data of the thin disc; ii) the parallel model can provide a good explanation of the evolution of N abundance in the thick and thin discs, and we confirm that the thick disc has evolved much faster than the thin disc, in agreement with the results from the abundance patterns of other chemical elements; and iii) finally, we present new model predictions for N evolution in the Galactic bulge, and we show that the observations in bulge stars can be explained if massive stars rotate fast during the earliest phases of Galactic evolution, in agreement with findings from the abundance pattern of carbon

Abundance gradients along the Galactic disc from chemical evolution models

In this paper, we study the formation and chemical evolution of the Milky Way disc with particular focus on the abundance patterns ([\u3b1/Fe] versus [Fe/H]) at different Galactocentric distances, the present-time abundance gradients along the disc, and the time evolution of abundance gradients. We consider the chemical evolution models for the Galactic disc developed by Grisoni et al. for the solar neighbourhood, both the two-infall and the one-infall ones, and we extend our analysis to the other Galactocentric distances. In particular, we examine the processes that mainly influence the formation of the abundance gradients: the inside-out scenario, a variable star formation efficiency, and radial gas flows. We compare our model results with recent abundance patterns obtained along the Galactic disc from the APOGEE survey and with abundance gradients observed from Cepheids, open clusters, H II regions, and PNe. We conclude that the inside-out scenario is a key ingredient but cannot be the only one to explain abundance patterns at different Galactocentric distances and abundance gradients. Further ingredients, such as radial gas flows and variable star formation efficiency, are needed to reproduce the observed features in the thin disc. The evolution of abundance gradients with time is also shown, although firm conclusions cannot still be drawn

Chemical evolution of the Milky Way: constraints on the formation of the thick and thin discs

We study the evolution of Milky Way thick and thin discs in the light of the most recent observational data. In particular, we analyse abundance gradients of O, N, Fe, and Mg along the thin disc as well as the [Mg/Fe] versus [Fe/H] relations and the metallicity distribution functions at different Galactocentric distances. We run several models starting from the two-infall paradigm, assuming that the thick and thin discs formed by means of two different infall episodes, and we explore several physical parameters, such as radial gas flows, variable efficiency of star formation, different times for the maximum infall on to the disc, different distributions of the total surface mass density of the thick disc, and enriched gas infall. Our best model suggests that radial gas flows and variable efficiency of star formation should be acting together with the inside-out mechanism for the thin disc formation. The time-scale for maximum infall on to the thin disc, which determines the gap between the formation of the two discs, should be tmax ≃ 3.25 Gyr. The thick disc should have an exponential, small-scale length density profile and gas infall on the inner thin disc should be enriched. We also compute the evolution of Gaia–Enceladus system and study the effects of possible interactions with the thick and thin discs. We conclude that the gas lost by Enceladus or even part of it could have been responsible for the formation of the thick disc but not the thin disc

APOGEE DR16: A multi-zone chemical evolution model for the Galactic disc based on MCMC methods

Context. The analysis of the latest release of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE DR16) data suggests the existence of a clear distinction between two sequences of disc stars at different Galactocentric distances in the [α/Fe] versus [Fe/H] abundance ratio space: the so-called high-α sequence, classically associated with an old population of stars in the thick disc with high average [α/Fe], and the low-α sequence, which mostly comprises relatively young stars in the thin disc with low average [α/Fe].Aims. We aim to constrain a multi-zone two-infall chemical evolution model designed for regions at different Galactocentric distances using measured chemical abundances from the APOGEE DR16 sample.Methods. We performed a Bayesian analysis based on a Markov chain Monte Carlo method to fit our multi-zone two-infall chemical evolution model to the APOGEE DR16 data.Results. An inside-out formation of the Galaxy disc naturally emerges from the best fit of our two-infall chemical-evolution model to APOGEE-DR16: Inner Galactic regions are assembled on shorter timescales compared to the external ones. In the outer disc (with radii R > 6 kpc), the chemical dilution due to a late accretion event of gas with a primordial chemical composition is the main driver of the [Mg/Fe] versus [Fe/H] abundance pattern in the low-α sequence. In the inner disc, in the framework of the two-infall model, we confirm the presence of an enriched gas infall in the low-α phase as suggested by chemo-dynamical models. Our Bayesian analysis of the recent APOGEE DR16 data suggests a significant delay time, ranging from ∼3.0 to 4.7 Gyr, between the first and second gas infall events for all the analysed Galactocentric regions. The best fit model reproduces several observational constraints such as: (i) the present-day stellar and gas surface density profiles; (ii) the present-day abundance gradients; (iii) the star formation rate profile; and (iv) the solar abundance values.Conclusions. Our results propose a clear interpretation of the [Mg/Fe] versus [Fe/H] relations along the Galactic discs. The signatures of a delayed gas-rich merger which gives rise to a hiatus in the star formation history of the Galaxy are impressed in the [Mg/Fe] versus [Fe/H] relation, determining how the low-α stars are distributed in the abundance space at different Galactocentric distances, which is in agreement with the finding of recent chemo-dynamical simulations

Galactic Archaeology with asteroseismic ages: Evidence for delayed gas infall in the formation of the Milky Way disc

Context. Precise stellar ages from asteroseismology have become available and can help to set stronger constraints on the evolution of the Galactic disc components. Recently, asteroseismology has confirmed a clear age difference in the solar annulus between two distinct sequences in the [α/Fe] versus [Fe/H] abundance ratios relation: the high-α and low-α stellar populations. Aims: We aim to reproduce these new data with chemical evolution models including different assumptions for the history and number of accretion events. Methods: We tested two different approaches: a revised version of the "two-infall" model where the high-α phase forms by a fast gas accretion episode and the low-α sequence follows later from a slower gas infall rate, and the parallel formation scenario where the two disc sequences form coevally and independently. Results: The revised two-infall model including uncertainties in age and metallicity is capable of reproducing: i) the [α/Fe] versus [Fe/H] abundance relation at different Galactic epochs, ii) the age-metallicity relation and the time evolution [α/Fe]; iii) the age distribution of the high-α and low-α stellar populations, iv) the metallicity distribution function. The parallel approach is not capable of properly reproducing the stellar age distribution, in particular at old ages. Conclusions: The best chemical evolution model is the revised two-infall one, where a consistent delay of ∼4.3 Gyr in the beginning of the second gas accretion episode is a crucial assumption to reproduce stellar abundances and ages

The AMBRE project: chemical evolution models for the Milky Way thick and thin discs

We study the chemical evolution of the thick and thin discs of the Galaxy by comparing detailed chemical evolution models with recent data from the Archéologie avec Matisse Basée sur les aRchives de l'ESO project. The data suggest that the stars in the thick and thin discs form two distinct sequences with the thick disc stars showing higher [α/Fe] ratios. We adopt two different approaches to model the evolution of thick and thin discs. In particular, we adopt (i) a two-infall approach where the thick disc forms fast and before the thin disc and by means of a fast gas accretion episode, whereas the thin disc forms by means of a second accretion episode on a longer time-scale; (ii) a parallel approach, where the two discs form in parallel but at different rates. By comparing our model results with the observed [Mg/Fe] versus [Fe/H] and the metallicity distribution functions in the two Galactic components, we conclude that the parallel approach can account for a group of α-enhanced metal-rich stars present in the data, whereas the two-infall approach cannot explain these stars unless they are the result of stellar migration. In both approaches, the thick disc has formed on a time-scale of accretion of 0.1 Gyr, whereas the thin disc formed on a time-scale of 7 Gyr in the solar region. In the two-infall approach, a gap in star formation between the thick and thin disc formation of several hundreds of Myr should be present, at variance with the parallel approach where no gap is present

Impact of very massive stars on the chemical evolution of extremely metal-poor galaxies

Context. In recent observations of extremely metal-poor, low-mass, starburst galaxies, almost solar Fe/O ratios are reported, despite N/O ratios consistent with the low metallicity. Aims: We aim to investigate if the peculiar Fe/O ratios can be a distinctive signature of an early enrichment produced by very massive objects dying as pair-instability supernova (PISN). Methods: We ran chemical evolution models with yields that account for the contribution by PISN. We used both the non-rotating stellar yields from a recent study and new yields from rotating very massive stars calculated specifically for this work. We also searched for the best initial mass function (IMF) that is able to reproduce the observations. Results: We can reproduce the observations by adopting a bi-modal IMF and by including an initial burst of rotating very massive stars. Only with a burst of very massive stars can we reproduce the almost solar Fe/O ratios at the estimated young ages. We also confirm that rotation is absolutely needed to concomitantly reproduce the observed N/O ratios. Conclusions: These results stress the importance of very massive stars in galactic chemical evolution studies and strongly support a top-heavy initial mass function in the very early evolutionary stages of metal-poor starburst galaxies

The effects of the initial mass function on Galactic chemical enrichment

Context. We have been seeing mounting evidence that the stellar initial mass function (IMF) might extend far beyond the canonical Mi ∼ 100 M⊙ limit, but the impact of such a hypothesis on the chemical enrichment of galaxies is yet to be clarified. Aims. We aim to address this question by analysing the observed abundances of thin- and thick-disc stars in the Milky Way with chemical evolution models that account for the contribution of very massive stars dying as pair instability supernovae. Methods. We built new sets of chemical yields from massive and very massive stars up to Mi ∼ 350 M⊙ by combining the wind ejecta extracted from our hydrostatic stellar evolution models with explosion ejecta from the literature. Using a simple chemical evolution code, we analysed the effects of adopting different yield tables by comparing predictions against observations of stars in the solar vicinity. Results. After several tests, we set our focus on the [O/Fe] ratio that best separates the chemical patterns of the two Milky Way components. We find that with a standard IMF, truncated at Mi ∼ 100 M⊙, we can reproduce various observational constraints for thin-disc stars; however, the same IMF fails to account for the [O/Fe] ratios of thick-disc stars. The best results are obtained by extending the IMF up to Mi = 350 M⊙, while including the chemical ejecta of very massive stars in the form of winds and pair instability supernova (PISN) explosions. Conclusions. Our study indicates that PISN may have played a significant role in shaping the chemical evolution of the thick disc of the Milky Way. Including their chemical yields makes it easier to reproduce not only the level of the α-enhancement, but also the observed slope of thick-disc stars in the [O/Fe] vs. [Fe/H] diagram. The bottom line is that the contribution of very massive stars to the chemical enrichment of galaxies is potentially quite important and should not be neglected in models of chemical evolution

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