Detailed Chemical Evolution Models Research Articles

(Re)mind the gap: A hiatus in star formation history unveiled by APOGEE DR17

Context. Analysis of several spectroscopic surveys indicates the presence of a bimodality between the disc stars in the abundance ratio space of [α/Fe] versus [Fe/H]. The two stellar groups are commonly referred to as the high-α and low-α sequences. Some models capable of reproducing such a bimodality invoke the presence of a hiatus in the star formation history in our Galaxy, whereas other models explain the two sequences by means of stellar migration. Aims. Our aim is to show that the existence of the gap in the star formation rate between high-α and low-α is evident in the stars of APOGEE DR17, if one plots [Fe/α] versus [α/H], confirming previous suggestions. We then try to interpret the data by means of detailed chemical models. Methods. We compare the APOGEE DR17 red giant stars with the predictions of a detailed chemical evolution model based on the two-infall paradigm, taking into account also the possible accretion of dwarf satellites. Results. The APOGEE DR17 abundance ratios [Fe/α] versus [α/H] exhibit a sharp increase in [Fe/α] at a nearly constant [α/H] (where α elements considered are Mg, Si, O) during the transition between the two disc phases. This observation strongly supports the hypothesis that a hiatus in star formation occurred during this evolutionary phase. Notably, the most pronounced growth in the [Fe/α] versus [α/H] relation is observed for oxygen, as this element is exclusively synthesised in core-collapse supernovae. The revised version of the two-infall chemical evolution model proposed in this study reproduces the APOGEE DR17 abundance ratios better than before. Particularly noteworthy is the model’s ability to predict the hiatus in the star formation between the two infalls of gas, which form the thick and thin disc, respectively, and thus generate abundance ratios compatible with APOGEE DR17 data. Conclusions. We show that the signature of a hiatus in the star formation is imprinted in the APOGEE DR17 abundance ratios. A chemical model predicting a pause in the star formation of a duration of roughly 3.5 Gyr, and in which the high-α disc starts forming from pre-enriched gas by a previous encounter with a dwarf galaxy, could well explain the observations

Mapping radial abundance gradients with Gaia -ESO open clusters. Evidence of recent gas accretion in the Milky Way disk

Recent evidence from spectroscopic surveys points towards the presence of a metal-poor, young stellar population in the low-alpha , chemically thin disk. In this context, the investigation of the spatial distribution and time evolution of precise, unbiased abundances is fundamental to disentangle the scenarios of formation and evolution of the Galaxy. We study the evolution of abundance gradients in the Milky Way by taking advantage of a large sample of open star clusters, which are among the best tracers for this purpose. In particular, we used data from the last release of the Gaia -ESO survey. We performed a careful selection of open cluster member stars, excluding those members that may be affected by biases in spectral analysis. We compared the cleaned open cluster sample with detailed chemical evolution models for the Milky Way, using well-tested stellar yields and prescription for radial migration. We tested different scenarios of Galaxy evolution to explain the data, namely, the two-infall and the three-infall frameworks, which suggest the chemical thin disk is formed by one or two subsequent gas accretion episodes, respectively. With the performed selection in cluster member stars, we still find a metallicity decrease between intermediate-age ($1<$ Age/Gyr $<3$) and young (Age $<1$ Gyr) open clusters. This decrease cannot be explained in the context of the two-infall scenario, even by accounting for the effect of migration and yield prescriptions. The three-infall framework, with its late gas accretion in the last 3 Gyr, is able to explain the low metallic content in young clusters. However, we have invoked a milder metal dilution for this gas infall episode relative to previous findings. To explain the observed low metallic content in young clusters, we propose that a late gas accretion episode triggering a metal dilution would have taken place, extending the framework of the three-infall model for the first time to the entire Galactic disk.

Galactic archaeology with [Mg/Mn] versus [Al/Fe] abundance ratios

Context. The diagram depicting the abundance ratios [Mg/Mn] versus [Al/Fe] has gained significant attention in recent literature as a valuable tool for exploring fundamental aspects of the evolution of the Milky Way (MW) and the nearby dwarf galaxies. In particular, this combination of elements is supposed to be highly sensitive to the history of star formation, unveiled by the imprints left on those abundances. Unfortunately, a complete discussion on the uncertainties associated with these nuclei is still missing, making it difficult to know how reliable the associated results are. Aims. In this work, we aim to analyse, by means of detailed chemical evolution models, the uncertainties related to the nucleosynthesis of Mg, Al, Mn, and Fe to show how different yield prescriptions can substantially affect the trends in the [Mg/Mn] versus [Al/Fe] plane. In fact, if different nucleosynthesis assumptions produce conflicting results, then the [Mg/Mn] versus [Al/Fe] diagram does not represent a strong diagnostic for the star formation history (SFH) of a galaxy. Methods. We discuss the results on the [Mg/Mn] versus [Al/Fe] diagram, as predicted by several MW and Large Magellanic Cloud (LMC) chemical evolution models adopting different nucleosynthesis prescriptions. Results. The results show that the literature nucleosynthesis prescriptions require some corrective factors to reproduce the APOGEE DR17 abundances of Mg, Al, and Mn in the MW and that the same factors can also improve the results for the LMC. In particular, we show that by modifying the massive star yields of Mg and Al, the behaviour of the [Mg/Mn] versus [Al/Fe] plot changes substantially. Conclusions. In conclusion, by changing the chemical yields within their error bars, one obtains trends that differ significantly, making it very difficult to draw any reliable conclusion on the SFH of galaxies. The proposed diagram is therefore very uncertain from a theoretical point of view, and it could represent a good diagnostic for star formation, only if the uncertainties on the nucleosynthesis of the above-mentioned elements (Mg, Mn, Al, and Fe) could be reduced by future stellar calculations.

Metal and dust evolution in ALMA REBELS galaxies: insights for future JWST observations

ABSTRACT Atacama Large Millimeter/submillimeter Array (ALMA) observations revealed the presence of significant amounts of dust in the first Gyr of Cosmic time. However, the metal and dust build-up picture remains very uncertain due to the lack of constraints on metallicity. JWST has started to reveal the metal content of high-redshift targets, which may lead to firmer constraints on high-redshift dusty galaxies evolution. In this work, we use detailed chemical and dust evolution models to explore the evolution of galaxies within the ALMA Reionization Era Bright Emission Line Survey (REBELS) survey, testing different metallicity scenarios that could be inferred from JWST observations. In the models, we track the build-up of stellar mass using non-parametric star formation histories for REBELS galaxies. Different scenarios for metal and dust evolution are simulated by allowing different prescriptions for gas flows and dust processes. The model outputs are compared with measured dust scaling relations, by employing metallicity-dependent calibrations for the gas mass based on the [C ii] 158 μm line. Independently of the galaxies metal content, we found no need for extreme dust prescriptions to explain the dust masses revealed by ALMA. However, different levels of metal enrichment will lead to different dominant dust production mechanisms, with stardust production dominant over other interstellar medium dust processes only in the metal-poor case. This points out how metallicity measurements from JWST will significantly improve our understanding of the dust build-up in high-redshift galaxies. We also show that models struggle to reproduce observables such as dust-to-gas and dust-to-stellar ratios simultaneously, possibly indicating an overestimation of the gas mass through current calibrations, especially at high metallicities.

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

ABSTRACT We investigate the origin of neutron-capture elements by analysing their abundance patterns and radial gradients in the Galactic thin disc. We adopt a detailed two-infall chemical evolution model for the Milky Way, including state-of-the-art nucleosynthesis prescriptions for neutron-capture elements. We consider r-process nucleosynthesis from merging neutron stars (MNS) and 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] versus [Fe/H] diagram is reproduced by both prompt and delayed sources, with the prompt source dominating Eu production; (ii) rotation in massive stars significantly contributes to the first peak s-process elements, but MNS and MR-SNe are necessary to match the observations; and (iii) our model slightly underpredicts Mo and Nd, while accurately reproducing the [Pr/Fe] versus [Fe/H] trend. Regarding the radial gradients, we find that: (i) our predicted [Fe/H] gradient slope agrees with observations from Gaia-ESO and other high-resolution spectroscopic surveys; (ii) the predicted [Eu/H] radial gradient slope is steeper than the observed one, regardless of how quick the production of Eu is, prompting discussion on different Galaxy-formation scenarios and stellar radial migration effects; and (iii) elements in the second s-process peak as well as Nd and Pr exhibit a plateau at low-Galactocentric distances, likely due to enhanced enrichment from LIMS in the inner regions.

GalCEM. I. An Open-source Detailed Isotopic Chemical Evolution Code

This is the first of a series of papers that will introduce a user-friendly, detailed, and modular Galactic Chemical Evolution Model, GalCEM, that tracks isotope masses as a function of time in a given galaxy. The list of tracked isotopes automatically adapts to the complete set provided by the input yields. The present iteration of GalCEM tracks 86 elements broken down into 451 isotopes. The prescription includes massive stars, low-to-intermediate-mass stars, and Type Ia supernovae as enrichment channels. We have developed a preprocessing tool that extracts multidimensional interpolation curves from the input yield tables. These interpolation curves improve the computation speeds of the full convolution integrals, which are computed for each isotope and for each enrichment channel. We map the integrand quantities onto consistent array grids in order to perform the numerical integration at each time step. The differential equation is solved with a fourth-order Runge–Kutta method. We constrain our analysis to the evolution of all light and intermediate elements from carbon to zinc, and lithium. Our results are consistent up to the extremely metal-poor regime with Galactic abundances. We provide tools to track the mass rate change of individual isotopes on a typical spiral galaxy with a final baryonic mass of 5 × 1010 M ⊙. Future iterations of the work will extend to the full periodic table by including the enrichment from neutron-capture channels as well as spatially dependent treatments of galaxy properties. GalCEM is publicly available at https://github.com/egjergo/GalCEM.

Mg/Fe] ratios in the solar neighbourhood: Stellar yields and chemical evolution scenarios

Context. The [Mg/Fe] abundance ratios are a fundamental fossil signature used to trace the chemical evolution of the disc and to divide it into low-α and high-α populations. Despite the huge observational and theoretical efforts, discrepancies between models and data are still present and several explanations have been put forward to explain the [α/Fe] bimodality. Aims. In this work we take advantage of a new AMBRE:HARPS dataset, which provides new more precise [Mg/Fe] estimations and reliable stellar ages for a subsample of stars, to study the [α/Fe] bimodality and the evolution of the solar neighbourhood. Methods. The data are compared with detailed chemical evolution models for the Milky Way, exploring the most used prescriptions for stellar yields and different formation scenarios for the Galactic disc (i.e. the delayed two-infall and the parallel models), including prescriptions for stellar radial migration. Results. We see that most of the stellar yield prescriptions struggle to reproduce the observed trend of the data and that semi-empirical yields describe best the [Mg/Fe] evolution in the thick and thin discs. In particular, most of the yields still predict a steeper decrease of the [Mg/Fe] ratio at high metallicity than shown by the data. The bulk of the data are well reproduced by the parallel and two-infall scenarios, but both scenarios have problems in explaining the most metal-rich and metal-poor tails of the low-α data. These tails can be explained in light of radial migration from the inner and outer disc regions, respectively. Conclusions. Despite the evidence of stellar migration, it is difficult to estimate the actual contribution of stars from other parts of the disc to the solar vicinity in the data we adopt. However, the comparison between data and models suggests that peculiar histories of star formation, such as that of the two-infall model, are still needed to reproduce the observed distribution of stars.

Detailed Chemical Abundances of Star Clusters in the Large Magellanic Cloud

We derive the first detailed chemical abundances of three star clusters in the Large Magellanic Cloud (LMC), NGC 1831 (436 ± 22 Myr), NGC 1856 (350 ± 18 Myr), and [SL63]268 (1230 ± 62 Myr) using integrated-light spectroscopic observations obtained with the Magellan Echelle spectrograph on Magellan Baade telescope. We derive [Fe/H], [Mg/Fe], [Ti/Fe], [Ca/Fe], [Ni/Fe], [Mn/Fe], [Cr/Fe], and [Na/Fe] for the three clusters. Overall, our results match the LMC abundances obtained in the literature as well as those predicted by detailed chemical evolution models. For clusters NGC 1831 and NGC 1856, the [Mg/Fe] ratios appear to be slightly depleted compared with [Ca/Fe] and [Ti/Fe]. This could be hinting at the well-known Mg–Al abundance anti-correlation observed in several Milky Way globular clusters. We note, however, that higher signal-to-noise observations are needed to confirm such a scenario, particularly for NGC 1831. We also find a slightly enhanced integrated-light [Na/Fe] ratio for cluster [SL63]268 compared with those from the LMC field stars, possibly supporting a scenario of intracluster abundance variations. We stress that detailed abundance analysis of individual stars in these LMC clusters is required to confirm the presence or absence of multiple stellar populations.

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

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

Evolution of neutron capture elements in dwarf galaxies

ABSTRACT We study the evolution of europium (Eu) and barium (Ba) abundances in Local Group dwarf spheroidal and ultrafaint dwarf galaxies by means of detailed chemical evolution models and compare our results with new sets of homogeneous abundances. The adopted models include gas infall and outflow and have been previously tested. We investigate several production scenarios for r-process elements: merging neutron stars and magnetorotational-driven supernovae. Production of Ba through the main s-process acting in low- and intermediate-mass stars is considered as well. We also test different sets of nucleosynthesis yields. For merging neutron stars we adopt either a constant and short delay time for merging or a delay time distribution function. Our simulations show that (i) if r-process elements are produced only by a quick source, it is possible to reproduce the [Eu/Fe] versus [Fe/H], but those models fail in reproducing the [Ba/Fe] versus [Fe/H]. (ii) If r-process elements are produced only with longer delays the opposite happens. (iii) If both a quick source and a delayed one are adopted, such as magnetorotational-driven supernovae and merging neutron stars with a delay time distribution, the [Eu/Fe] abundance pattern is successfully reproduced, but models still fail in reproducing the [Ba/Fe]. (iv) On the other hand, the characteristic abundances of Reticulum II can be reproduced only if both the Eu and the r-process fraction of Ba are produced on short and constant time delays during a single merging event. We discuss also other possible interpretations, including an inhomogeneous mixing of gas that might characterize this galaxy.

The evolution of Lithium: implications of a universal Spite plateau

ABSTRACT The cosmological 7Li problem consists in explaining why the primordial Li abundance, as predicted by the standard Big Bang nucleosynthesis theory with constraints from WMAP and Planck, is a factor of 3 larger than the Li abundance measured in the stars of the Spite plateau defined by old, warm dwarf stars of the Milky Way halo. Several explanations have been proposed to explain this difference, including various Li depletion processes as well as non-standard Big Bang nucleosynthesis, but the main question remains unanswered. In this paper, we present detailed chemical evolution models for dwarf spheroidal and ultra faint galaxies, compute the galactic evolution of 7Li abundance in these objects, and compare it with observations of similar objects. In our models, Li is mainly produced by novae and cosmic rays, and to a minor extent, by low and intermediate mass stars. We adopt the yield combination that best fits the Li abundances in the Milky Way stars. It is evident that the observations of dwarf objects define a Spite plateau, identical to that observed in the Milky Way, thus suggesting that the Spite plateau could be a universal feature and its meaning should be discussed. The predictions of our models for dwarf galaxies are obtained by assuming as Li primordial abundance either the one detected in the atmospheres of the oldest halo stars (Spite plateau; A(Li) ∼ 2.2 dex), or the one from cosmological observations (WMAP; A(Li) ∼ 2.66 dex). Finally, we discuss the implications of the universality of the Spite plateau results.

Detailed chemical evolution models of the Tucana dSph galaxy

ABSTRACT We use semi-analytical methods to obtain detailed chemical evolution models (CEMs) for the dwarf spheroidal (dSph) galaxy Tucana. Published star formation rates and the age–metallicity relationship are used to constrain the observables. The results show that Tucana: (i) behaved like a closed box for 75 per cent of its life, (ii) had either a primordial-gas accretion or a metal-rich wind during 15 per cent of its life (between 0.5 and 2.0 Gyr), and (iii) lost 95 per cent of its gas through a well-mixed wind at t ∼ 4.5 Gyr. Specifically, we find two CEMs: the metal-dilution model and the metal-loss model, which differ mainly during the range 0.5–2.0 Gyr. In order to discriminate between these CEMs, we compare the predicted [Xi/Fe]–[Fe/H] trends, which differ less than the average error of the observed trends for other dSphs of the Local Group. Furthermore, the models predict very different metallicity distribution functions. Therefore, an observational metallicity distribution function for Tucana is essential in order to discriminate between the metal-dilution and the metal-loss scenarios. In addition, because the difference of [Xi/Fe] between the two models is less than the average of the errors observed for others dSph and ultra-faint dwarf galaxies of the Local Group, greater precision is required in future observations, so that the errors are less than the difference between the models, thus enabling discrimination between them.

Fluorine in the solar neighbourhood: modelling the Galactic thick and thin discs

ABSTRACT We investigate the evolution of the abundance of fluorine in the Milky Way thick and thin discs by means of detailed chemical evolution models compared with recent observational data. The chemical evolution models adopted here have already been shown to fit the observed abundance patterns of CNO and α-elements as well as the metallicity distribution functions for the Galactic thick and thin disc stars. We apply them here to the study of the origin and evolution of fluorine, which is still a matter of debate. First, we study the importance of the various sites proposed for the production of fluorine. Then, we apply the reference models to follow the evolution of the two different Galactic components. We conclude that rotating massive stars are important producers of F and they can set a plateau in F abundance below [Fe/H] = −0.5 dex, though its existence for [Fe/H]<−1 has yet to be confirmed by extensive observations of halo stars. In order to reproduce the F abundance increase in the discs at late times, instead, a contribution from lower mass stars – single asymptotic giant branch stars and/or novae – is required. The dichotomy between the thick and thin discs is more evident in the [F/O] versus [O/H] plot than in the [F/Fe] versus [Fe/H] one, and we confirm that the thick disc has evolved much faster than the thin disc, in agreement with findings from the abundance patterns of other chemical elements.

Heavy element evolution in the inner regions of the Milky Way

ABSTRACT We present results for the evolution of the abundances of heavy elements (O, Mg, Al, Si, K, Ca, Cr, Mn, Ni, and Fe) in the inner Galactic regions (RGC ≤ 4 kpc). We adopt a detailed chemical evolution model already tested for the Galactic bulge and compare the results with Apache Point Observatory Galactic Evolution Experiment data. We start with a set of yields from the literature that are considered the best to reproduce the abundance patterns in the solar vicinity. We find that, in general, the predicted trends nicely reproduce the data but in some cases either the trend or the absolute values of the predicted abundances need to be corrected, even by large factors, in order to reach the best agreement. We suggest how the current stellar yields should be modified to reproduce the data and we discuss whether such corrections are reasonable in the light of the current knowledge of stellar nucleosynthesis. However, we also critically discuss the observations. Our results suggest that Si, Ca, Cr, and Ni are the elements for which the required corrections are the smallest, while for Mg and Al moderate modifications are necessary. On the other hand, O and K need the largest corrections to reproduce the observed patterns, a conclusion already reached for solar vicinity abundance patterns, with the exception of oxygen. For Mn, we apply corrections already suggested in previous works.

On the origin of dust in galaxy clusters at low-to-intermediate redshift

ABSTRACT Stacked analyses of galaxy clusters at low-to-intermediate redshift show signatures attributable to dust, but the origin of this dust is uncertain. We test the hypothesis that the bulk of cluster dust derives from galaxy ejecta. To do so, we employ dust abundances obtained from detailed chemical evolution models of galaxies. We integrate the dust abundances over cluster luminosity functions (one-slope and two-slope Schechter functions). We consider both a hierarchical scenario of galaxy formation and an independent evolution of the three main galactic morphologies: elliptical/S0, spiral and irregular. We separate the dust residing within galaxies from the dust ejected in the intracluster medium. To the latter, we apply thermal sputtering. The model results are compared to low-to-intermediate redshift observations of dust masses. We find that in any of the considered scenarios, elliptical/S0 galaxies contribute negligibly to the present-time intracluster dust, despite producing the majority of gas-phase metals in galaxy clusters. Spiral galaxies, instead, provide both the bulk of the spatially unresolved dust and of the dust ejected into the intracluster medium. The total dust-to-gas mass ratio in galaxy clusters amounts to 10−4, while the intracluster medium dust-to-gas mass ratio amounts to 10−6 at most. These dust abundances are consistent with the estimates of cluster observations at 0.2 < z < 1. We propose that galactic sources, spiral galaxies in particular, are the major contributors to the cluster dust budget.

Galactic Archaeology at High Redshift: Inferring the Nature of GRB Host Galaxies from Abundances

Abstract We identify the nature of high-redshift long gamma-ray burst (LGRB) host galaxies by comparing the observed abundance ratios in the interstellar medium with detailed chemical evolution models accounting for the presence of dust. We compare abundance data from LGRB afterglow spectra to abundance patterns as predicted by our models for different galaxy types. We analyze [X/Fe] abundance ratios (where X is C, N, O, Mg, Si, S, Ni, Zn) as functions of [Fe/H]. Different galaxies (irregulars, spirals, spheroids) are, in fact, characterized by different star formation histories, which produce different [X/Fe] versus [Fe/H] relations (“time-delay model”). This allows us to identify the star formation history of the host galaxies and to infer their age (i.e., the time elapsed from the beginning of star formation) at the time of the GRB events. Unlike previous works, we use newer models in which we adopt updated stellar yields and prescriptions for dust production, accretion, and destruction. We consider a sample of seven LGRB host galaxies. Our results suggest that two of them (GRB 050820, GRB 120815A) are star-forming spheroids, two (GRB 081008, GRB 161023A) are spirals, and three (GRB 090926A, GRB 050730, GRB 120327A) are irregulars. The inferred ages of the considered host galaxies span from 10 Myr to slightly more than 1 Gyr.

Modelling the chemical evolution of Zr, La, Ce, and Eu in the Galactic discs and bulge

ABSTRACT We study the chemical evolution of Zr, La, Ce, and Eu in the Milky Way discs and bulge by means of chemical evolution models compared with spectroscopic data. We consider detailed chemical evolution models for the Galactic thick disc, thin disc, and bulge, which have been already tested to reproduce the observed [α/Fe] versus [Fe/H] diagrams and metallicity distribution functions for the three different components, and we apply them to follow the evolution of neutron capture elements. In the [Eu/Fe] versus [Fe/H] diagram, we observe and predict three distinct sequences corresponding to the thick disc, thin disc, and bulge, similar to what happens for the α-elements. We can nicely reproduce the three sequences by assuming different time-scales of formation and star formation efficiencies for the three different components, with the thin disc forming on a longer time-scale of formation with respect to the thick disc and bulge. On the other hand, in the [X/Fe] versus [Fe/H] diagrams for Zr, La, and Ce, the three populations are mixed and also from the model point of view there is an overlapping between the predictions for the different Galactic components, but the observed behaviour can be also reproduced by assuming different star formation histories in the three components. In conclusions, it is straightforward to see how different star formation histories can lead to different abundance patterns and also looking at the abundance patterns of neutron capture elements can help in constraining the history of formation and evolution of the major Galactic components.

Evolution of lithium in the Milky Way halo, discs, and bulge

Abstract In this work, we study the Galactic evolution of lithium by means of chemical evolution models in the light of the most recent spectroscopic data from Galactic stellar surveys. We consider detailed chemical evolution models for the Milky Way halo, discs, and bulge, and we compare our model predictions with the most recent spectroscopic data for these different Galactic components. In particular, we focus on the decrease of lithium at high metallicity observed by the AMBRE Project, the Gaia-ESO Survey, and other spectroscopic surveys, which still remains unexplained by theoretical models. We analyse the various lithium producers and confirm that novae are the main source of lithium in the Galaxy, in agreement with other previous studies. Moreover, we show that, by assuming that the fraction of binary systems giving rise to novae is lower at higher metallicity, we can suggest a novel explanation to the lithium decline at super-solar metallicities: the aforementioned assumption is based on independent constraints on the nova system birth rate, which have been recently proposed in the literature. As regards the thick disc, it is less lithium enhanced due to the shorter time-scale of formation and higher star formation efficiency with respect to the thin disc; therefore, we have a faster evolution and the ‘reverse knee’ in the A(Li) versus [Fe/H] relation is shifted towards higher metallicities. Finally, we present our predictions about lithium evolution in the Galactic bulge, which, however, still need further data to be confirmed or disproved.

The origin of stellar populations in the Galactic bulge from chemical abundances

ABSTRACT In this work, we study the formation and chemical evolution of the Galactic bulge with particular focus on the abundance pattern ([Mg/Fe] versus [Fe/H]), metallicity, and age distribution functions. We consider detailed chemical evolution models for the Galactic bulge and inner disc, with the aim of shedding light on the connection between these components and the origin of bulge stars. In particular, we first present a model assuming a fast and intense star formation, with the majority of bulge stars forming on a time-scale less than 1 Gyr. Then we analyse the possibility of two distinct stellar populations in the bulge, as suggested by Gaia-ESO and APOGEE data. These two populations, one metal poor and the other metal rich, can have had two different origins: (i) the metal rich formed after a stop of ∼250 Myr in the star formation rate of the bulge or (ii) the metal-rich population is made of stars formed in the inner disc and brought into the bulge by the early secular evolution of the bar. We also examine the case of multiple starbursts in the bulge with consequent formation of multiple populations, as suggested by studies of microlensed stars. After comparing model results and observations, we suggest that the most likely scenario is that there are two main stellar populations, both made mainly by old stars (>10 Gyr), with the metal-rich and younger one formed from inner thin disc stars, in agreement with kinematical arguments. However, on the basis of dynamical simulations, we cannot completely exclude that the second population formed after a stop in the star formation during the bulge evolution, so that all the stars formed in situ.

A new delay time distribution for merging neutron stars tested against Galactic and cosmic data

The merging of two neutron stars (MNS) is thought to be the source of short gamma-ray bursts (SGRB) and gravitational wave transients, as well as the main production site of r-process elements like Eu. We have derived a new delay time distribution (DTD) for MNS systems from theoretical considerations and we have tested it against (i) the SGRB redshift distribution and (ii) the Galactic evolution of Eu and Fe, in particular the [Eu/Fe] vs [Fe/H] relation. For comparison, we also tested other DTDs, as proposed in the literature. To address the first item, we have convolved the DTD with the cosmic star formation rate, while for the second we have employed a detailed chemical evolution model of the Milky Way. We have also varied the DTD of Type Ia SNe (the main Fe producers), the contribution to Eu production from core-collapse SNe, as well as explored the effect of a dependence on the metallicity of the occurrence probability of MNS. Our main results can be summarized as follows: (i) the SGRB redshift distribution can be fitted using DTDs for MNS that produce average timescales of 300-500 Myr; (ii) if the MNS are the sole producers of the Galactic Eu and the occurrence probability of MNS is constant the Eu production timescale must be on the order of <30 Myr; (iii) allowing for the Eu production in core-collapse SNe, or adopting a metallicity-dependent occurrence probability, allow us to reproduce both observational constraints, but many uncertainties are still present in both assumptions.