Galactic Chemical Evolution Models Research Articles

Impact of the new 12C+12C reaction rate on presupernova nucleosynthesis* *Supported by the National Natural Science Foundation of China (11988101, 11890694), and the National Key R&D Program of China (2019YFA0405502). K. Nomoto is supported by the World Premier International Research Center Initiative (WPI), MEXT, Japan, and the Japan Society for the Promotion of Science (JSPS) KAKENHI grant (JP17K05382, JP20K04024, JP21H04499)

The 12C+12C reaction rate plays an essential role in stellar evolution and nucleosynthesis. Nevertheless, the uncertainties of this reaction rate are still large. We calculate a series of stellar evolution models with the near solar abundance from the zero-age main-sequence through presupernova stages for initial masses of 20 M to 40 M . The 12C+12C reaction rates from two different studies are used in our investigation. One is the rate obtained using the Trojan Horse Method (THM) by Tumino et al. [Nature 557(7707), 687 (2018)], and the other was obtained by Mukhamedzhanov et al. [Physical Review C 99(6), 064618 (2019)] (Muk19). Then, comparisons of the nucleosynthesis and presupernova isotopic abundances are conducted. In particular, we find that in the C burning shell, models with the THM produce a smaller amount of 23Na and some neutron-rich isotopes than Muk19. The difference in the abundance ratios of Na/Mg, S/Mg, Ar/Mg, and K/Mg between the two models are apparent. We compare Na/Mg obtained from our theoretical presupernovae models with Na/Mg in stellar atmospheres observed with high-resolution spectra as well as from the latest galactic chemical evolution model. Although Na/Mg obtained using the THM is within 2σ of the observed stellar ratio, the theoretical uncertainty on Na/Mg introduced by the uncertainty of the 12C+12C reaction rate is almost equivalent to the standard deviation of astronomical observations. Therefore, a more accurate 12C+12C reaction rate is crucial.

CHEMOUT: CHEMical complexity in star-forming regions of the OUTer Galaxy

Context. Nitrogen isotopic ratios are a key tool for tracing Galactic stellar nucleosynthesis. Aims. We present the first study of the 14N/15N abundance ratio in the outer regions of the Milky Way (namely, for galactocentric distances, RGC, from 12 up to 19 kpc), with the aim to study the stellar nucleosynthesis effects in the global Galactic trend. Methods. We analysed IRAM 30 m observations towards a sample of 35 sources in the context of the CHEMical complexity in star-forming regions of the OUTer Galaxy (CHEMOUT) project. We derived the 14N/15N ratios from HCN and HNC for 14 and 3 sources, respectively, using the J = 1–0 rotational transition of HN13C, H15NC, H13CN, and HC15N. Results. The results found in the outer Galaxy have been combined with previous measurements obtained in the inner Galaxy. We find an overall linear decreasing H13CN/HC15N ratio with increasing RGC. This translates to a parabolic 14N/15N ratio with a peak at 11 kpc. Updated Galactic chemical evolution models have been taken into account and compared with the observations. The parabolic trend of the 14N/15N ratio with RGC can be naturally explained (i) by a model that assumes novae as the main 15N producers on long timescales (≥1 Gyr) and (ii) by updated stellar yields for low- and intermediate-mass stars.

The C/M Ratio of AGB Stars in the Local Group Galaxies

The number ratio of carbon-rich to oxygen-rich asymptotic giant branch (AGB) stars (the so-called C/M ratio) is closely related to the evolution environment of the host galaxy. This work studies the C/M ratio in 14 galaxies within the Local Group with the most complete and clean sample of member stars identified in our previous works. The borderlines between carbon-rich AGB and oxygen-rich AGB stars as well as red supergiants are defined by Gaussian mixture model fitting to the number density in the (J−K)/K diagram for the member stars of the LMC and M33, and then applied to the other galaxies by shifting the difference in the position of tip red giant branch (TRGB). The C/M ratios are obtained after precise and consistent categorization. Although for galaxies with larger distance modulo there is greater uncertainty, the C/M ratio is clearly found to decrease with the color index (J−K)0 of TRGB as the indicator of metallicity, which agrees with previous studies and can be explained by the fact that carbon stars are more easily formed in a metal-poor environment. Furthermore, the C/M ratio within M33 is found to increase with galactocentric distance, which coincides with this scenario and the galactic chemical evolution model. On the other hand, the C/M ratio within M31 is found to decrease with galactocentric radius, which deserves further study.

Non-LTE abundance corrections for late-type stars from 2000 Å to 3 µm

Context. It is well known that cool star atmospheres depart from local thermodynamic equilibrium (LTE). Making an accurate abundance determination requires taking those effects into account, but the necessary non-LTE (hereafter NLTE) calculations are often lacking. Aims. Our goal is to provide detailed estimates of NLTE effects for FGK type stars for all spectral lines from the ultraviolet (UV) to the near infrared (NIR) that are potentially useful as abundance diagnostics. The first paper in this series focusses on the light elements Na, Mg, and Al. Methods. The code PySME was used to compute curves of growth for 2158 MARCS model atmospheres in the parameter range 3800 < Teff < 8000 K, 0.0 < log(g) < 5.5, and −5 < [Fe/H] < +0.5. Two microturbulence values, 1 and 2 kms−1, and nine abundance points spanning −1 < [X/Fe] < 1 for element X, are used to construct individual line curves of growth by calculating the equivalent widths of 35 Na lines, 134 Mg lines, and 34 Al lines. The lines were selected in the wavelength range between 2000 Å and 3 µm. Results. We demonstrate the power of the new grids with LTE and NLTE abundance analysis by means of equivalent width measurements of five benchmark stars; the Sun, Arcturus, HD 84937, HD 140283 and HD 122563. For Na, the NLTE abundances are lower than in LTE and show markedly reduced line-to-line scatter in the metal-poor stars. For Mg, we confirm previous reports of a significant ~0.25 dex LTE ionisation imbalance in metal-poor stars that is only slightly improved in NLTE (~0.18 dex). LTE abundances based on Mg II lines agree better with models of Galactic chemical evolution. For Al, NLTE calculations strongly reduce an ~0.6 dex ionisation imbalance seen in LTE for the metal-poor stars. The abundance corrections presented in this work are in good agreement with previous studies for the subset of lines that overlap, with the exception of strongly saturated lines. Conclusions. A consensus between different abundance diagnostics is the most powerful tool available to stellar spectroscopists to assess the accuracy of the models. Here we report that NLTE abundance analysis in general leads to improved agreement, in particular for metal-poor stars. The residual scatter is believed to be caused mainly by unresolved blends and/or poor atomic data, with the notable exception of Mg, which calls for further investigation.

Chemical Tracing and the Origin of Carbon in the Galactic Disk

A basic problem in studies of the evolution of chemical elements in galaxies is the uncertainties in the yields of elements produced by different types of stars. The possibilities of tracing the sites producing chemical elements and corresponding yields in stellar populations by studying ratios of abundances in stars of different ages and metallicities, with an approach with minimal assumptions concerning the yields, is explored by means of simple models of Galactic chemical evolution. Elemental abundances of carbon and oxygen, obtained by recent observations of samples of solar-type stars with estimated ages in the thin disk of the Galaxy, are analysed. Constraints on the yields from winds of intermediate-mass stars and of hot massive stars, including core-collapse supernovae, are derived. It is found that a dominating contribution of carbon from massive stars is most probable, although stars in the mass interval of two to three solar masses may have provided some amounts of carbon in the Sun. The results are consistent with those obtained by using theoretical yields and more elaborate models of Galactic evolution. The uncertainties as regards the mixing of stellar populations due to migration of stars in the Galactic disk may be important for the conclusions. Variations in the star formation rates, lack of chemical homogeneity in the Galactic gas, the inflow of gas from the intergalactic space and possible variations in the Initial mass function may also limit conclusions about the sites and their yields. Very accurate abundance ratios and the determination of stellar ages provide further important constraints on the yields.

The origin of elements: the need for UV spectra

Thanks to the long-term collaborations between nuclear and astrophysics, we have good understanding on stellar nucleosynthesis, except for the elements around Ti and some neutron-capture elements. From the comparison between observations and Galactic chemical evolution models, it is necessary to have the rapid neutron-capture process associated with core-collapse supernovae, although the explosion mechanism is unknown. The impact of rotating massive stars is also shown in this paper. Many of the key elements can be exclusively obtained in the UV, and therefore without UV spectra it would not be possible to fully understand the origin of elements in the universe.

Viability of Novae as Sources of Galactic Lithium

Of all the light elements, the evolution of lithium (Li) in the Milky Way is perhaps the most difficult to explain. Li is difficult to synthesize and is easily destroyed, making most stellar sites unsuitable for producing Li in sufficient quantities to account for the protosolar abundance. For decades, novae have been proposed as a potential explanation for this “Galactic Li problem,” and the recent detection of 7Be in the ejecta of multiple nova eruptions has breathed new life into this theory. In this work, we assess the viability of novae as dominant producers of Li in the Milky Way. We present the most comprehensive treatment of novae in a galactic chemical evolution code to date, testing theoretically and observationally derived nova Li yields by integrating metallicity-dependent nova ejecta profiles computed using the binary population synthesis code binary_c with the galactic chemical evolution code OMEGA+. We find that our galactic chemical evolution models which use observationally derived Li yields account for the protosolar Li abundance very well, while models relying on theoretical nova yields cannot reproduce the protosolar observation. A brief exploration of physical uncertainties including single-stellar yields, the metallicity resolution of our nova treatment, common-envelope physics, and nova accretion efficiencies indicates that this result is robust to physical assumptions. Scatter within the observationally derived Li yields in novae is identified as the primary source of uncertainty, motivating further observations of 7Be in nova ejecta.

Neutron-capture elements record the ordered chemical evolution of the disc over time

ABSTRACT An ensemble of chemical abundances probing different nucleosynthetic channels can be leveraged to build a comprehensive understanding of the chemical and structural evolution of the Galaxy. Using GALAH DR3 data, we seek to trace the enrichment by the supernovae Ia, supernovae II, asymptotic giant branch stars, and neutron-star mergers and/or collapsars nucleosynthetic sources by studying the [Fe/H], [α/Fe], [Ba/Fe], and [Eu/Fe] chemical compositions of ∼50 000 red giant stars, respectively. Employing small [Fe/H]–[α/Fe] cells, which serve as an effective reference-frame of supernovae contributions, we characterize the abundance-age profiles for [Ba/Fe] and [Eu/Fe]. Our results disclose that these age–abundance relations vary across the [Fe/H]–[α/Fe] plane. Within cells, we find negative age–[Ba/Fe] relations and flat age–[Eu/Fe] relations. Across cells, we see the slope of the age–[Ba/Fe] relations evolve smoothly and the [Eu/Fe] relations vary in amplitude. We subsequently model our empirical findings in a theoretical setting using the flexible Chempy Galactic chemical evolution (GCE) code, using the mean [Fe/H], [Mg/Fe], [Ba/Fe], and age values for stellar populations binned in [Fe/H], [Mg/Fe], and age space. We find that within a one-zone framework, an ensemble of GCE model parameters vary to explain the data. Using present day orbits from Gaia EDR3 measurements we infer that the GCE model parameters, which set the observed chemical abundance distributions, vary systematically across mean orbital radii. Under our modelling assumptions, the observed chemical abundances are consistent with a small gradient in the high-mass end of the initial mass function (IMF) across the disc, where the IMF is more top heavy towards the inner disc and more bottom heavy in the outer disc.

Mantle Degassing Lifetimes through Galactic Time and the Maximum Age Stagnant-lid Rocky Exoplanets Can Support Temperate Climates

The ideal exoplanets to search for life are those within a star’s habitable zone. However, even within the habitable zone, planets can still develop uninhabitable climate states. Sustaining a temperate climate over geologic (∼gigayear) timescales requires a planet to contain sufficient internal energy to power a planetary-scale carbon cycle. A major component of a rocky planet’s energy budget is the heat produced by the decay of radioactive elements, especially 40K, 232Th, 235U, and 238U. As the planet ages and these elements decay, this radiogenic energy source dwindles. Here we estimate the probability distribution of the amount of these heat-producing elements that enter into rocky exoplanets through Galactic history by combining the system-to-system variation seen in stellar abundance data with the results from Galactic chemical evolution models. From this, we perform Monte Carlo thermal evolution models that maximize the mantle cooling rate, thus allowing us to create a pessimistic estimate of lifetime a rocky, stagnant-lid exoplanet can support a global carbon cycle through Galactic history. We apply this framework to a sample of 17 likely rocky exoplanets with measured ages, seven of which we predict are likely to be actively degassing today, despite our pessimistic assumptions. For the remaining planets, including those orbiting TRAPPIST-1, we cannot confidently assume that they currently contain sufficient internal heat to support mantle degassing at a rate sufficient to sustain a global carbon cycle or temperate climate without additional tidal heating or undergoing plate tectonics.

Common Envelope Jet Supernova r-process Yields Can Reproduce [Eu/Fe] Abundance Evolution in the Galaxy

Abstract We use a numerical Galactic chemical evolution model and find that the common envelope jet supernova (CEJSN) r-process scenario can account for both the very early average ratio of europium to iron and its evolution at later times in the the Milky Way. In the CEJSN scenario a neutron star (NS) spirals in inside a red supergiant (RSG) star all the way to the core and destroys it. According to this scenario r-process isotopes are nucleosynthesized inside neutron-rich jets that the accretion disk around the NS launches inside the core. The merger of an NS with an RSG core already takes place in the very young Galaxy. We conclude that CEJSNe can be a major contributor to r-process nucleosynthesis.

The K2 Galactic Archaeology Program Data Release 3: Age-abundance Patterns in C1–C8 and C10–C18

Abstract We present the third and final data release of the K2 Galactic Archaeology Program (K2 GAP) for Campaigns C1–C8 and C10–C18. We provide asteroseismic radius and mass coefficients, κ R and κ M , for ∼19,000 red giant stars, which translate directly to radius and mass given a temperature. As such, K2 GAP DR3 represents the largest asteroseismic sample in the literature to date. K2 GAP DR3 stellar parameters are calibrated to be on an absolute parallactic scale based on Gaia DR2, with red giant branch and red clump evolutionary state classifications provided via a machine-learning approach. Combining these stellar parameters with GALAH DR3 spectroscopy, we determine asteroseismic ages with precisions of ∼20%–30% and compare age-abundance relations to Galactic chemical evolution models among both low- and high-α populations for α, light, iron-peak, and neutron-capture elements. We confirm recent indications in the literature of both increased Ba production at late Galactic times as well as significant contributions to r-process enrichment from prompt sources associated with, e.g., core-collapse supernovae. With an eye toward other Galactic archeology applications, we characterize K2 GAP DR3 uncertainties and completeness using injection tests, suggesting that K2 GAP DR3 is largely unbiased in mass/age, with uncertainties of 2.9% (stat.) ± 0.1% (syst.) and 6.7% (stat.) ± 0.3% (syst.) in κ R and κ M for red giant branch stars and 4.7% (stat.) ± 0.3% (syst.) and 11% (stat.) ± 0.9% (syst.) for red clump stars. We also identify percent-level asteroseismic systematics, which are likely related to the time baseline of the underlying data, and which therefore should be considered in TESS asteroseismic analysis.

Galactic Chemical Evolution of Radioactive Isotopes with an s-process Contribution

Abstract Analysis of inclusions in primitive meteorites reveals that several short-lived radionuclides (SLRs) with half-lives of 0.1–100 Myr existed in the early solar system (ESS). We investigate the ESS origin of 107Pd, 135Cs, and 182Hf, which are produced by slow neutron captures (the s-process) in asymptotic giant branch (AGB) stars. We modeled the Galactic abundances of these SLRs using the OMEGA+ galactic chemical evolution (GCE) code and two sets of mass- and metallicity-dependent AGB nucleosynthesis yields (Monash and FRUITY). Depending on the ratio of the mean-life τ of the SLR to the average length of time between the formations of AGB progenitors γ, we calculate timescales relevant for the birth of the Sun. If τ/γ ≳ 2, we predict self-consistent isolation times between 9 and 26 Myr by decaying the GCE predicted 107Pd/108Pd, 135Cs/133Cs, and 182Hf/180Hf ratios to their respective ESS ratios. The predicted 107Pd/182Hf ratio indicates that our GCE models are missing 9%–73% of 107Pd and 108Pd in the ESS. This missing component may have come from AGB stars of higher metallicity than those that contributed to the ESS in our GCE code. If τ/γ ≲ 0.3, we calculate instead the time (T LE) from the last nucleosynthesis event that added the SLRs into the presolar matter to the formation of the oldest solids in the ESS. For the 2 M ⊙, Z = 0.01 Monash model we find a self-consistent solution of T LE = 25.5 Myr.

Scaling relations and baryonic cycling in local star-forming galaxies

Gas accretion and stellar feedback processes link metal content, star formation, and gas and stellar mass (and the potential depth) in star-forming galaxies. Constraining this hypersurface has been challenging because of the need for measurements of H I and H2 gas masses spanning a broad parameter space. A recent step forward has been achieved through the Metallicity And Gas for Mass Assembly sample of local star-forming galaxies, which consists of homogeneously determined parameters and a significant quantity of dwarf galaxies, with stellar masses as low as ∼105 − 106 M⊙. Here, in the third paper of a series, we adopt a standard galactic chemical evolution model, with which we can quantify stellar-driven outflows. In particular, we constrain the difference between the mass-loading in accretion and outflows and the wind metal-loading factor. The resulting model reproduces very well the local mass–metallicity relation, and the observed trends of metallicity with gas fraction. Although the difference in mass loading between accreted and expelled gas is extremely difficult to constrain, we find indications that, on average, the amount of gas acquired through accretion is roughly the same as the gas lost through bulk stellar outflows, a condition roughly corresponding to a “gas equilibrium” scenario. In agreement with previous work, the wind metal-loading factor shows a steep increase toward lower mass and circular velocity, indicating that low-mass galaxies are more efficient at expelling metals, thus shaping the mass–metallicity relation. Effective yields are found to increase with mass up to an inflection mass threshold, with a mild decline at larger masses and circular velocities. A comparison of our results for metal loading in outflows with the expectations for their mass loading favors momentum-driven winds at low masses, rather than energy-driven ones.

The bi-modal 7Li distribution of the Milky Way’s thin-disk dwarf stars

Context. The lithium abundance, A(Li), in stellar atmospheres suffers from various enhancement and depletion processes during the star’s lifetime. While several studies have demonstrated that these processes are linked to the physics of stellar formation and evolution, the role that Galactic-scale events play in the galactic A(Li) evolution is not yet well understood. Aims. We aim to demonstrate that the observed A(Li) bi-modal distribution, in particular in the FGK-dwarf population, is not a statistical artefact and that the two populations connect through a region with a low number of stars. We also want to investigate the role that Galactic-scale events play in shaping the A(Li) distribution of stars in the thin disk. Methods. We use statistical techniques along with a Galactic chemical evolution model for A(Li) that includes most of the well-known 7Li production and depletion channels. Results. We confirm that the FGK main-sequence stars belonging to the Milky Way’s thin disk present a bi-modal A(Li) distribution. We demonstrate that this bi-modality can be generated by a particular Milky Way star formation history profile combined with the stellar evolution’s 7Li depletion mechanisms. We show that A(Li) evolution can be used as an additional proxy for the star formation history of our Galaxy.

Metallicity-dependent nucleosynthetic yields of Type Ia supernovae originating from double detonations of sub-MCh white dwarfs

Double detonations in sub-Chandrasekhar mass carbon-oxygen white dwarfs (WD) with helium shells ares potential explosion mechanisms for Type Ia supernovae. The mechanism consists of a shell detonation and subsequent core detonation. The focus of our study is the effect of the progenitor metallicity on the nucleosynthetic yields. For this, we computed and analyzed a set of 11 different models with varying core and shell masses at four different metallicities each. This results in a total of 44 models at metallicities between 0.01 Z⊙ and 3 Z⊙. Our models show a strong impact of the metallicity in the high-density regime. The presence of 22Ne causes a neutron-excess that shifts the production from 56Ni to stable isotopes such as 54Fe and 58Ni in the α-rich freeze-out regime. The isotopes of the metallicity implementation further serve as seed nuclei for additional reactions in the shell detonation. The production of 55Mn increases with metallicity, confirming the results of previous work. A comparison of elemental ratios relative to iron shows a good match to solar values for some models. Super-solar values are reached for Mn at 3 Z⊙ and solar values in some models at Z⊙. This indicates that the required contribution of Type Ia supernovae originating from Chandrasekhar-mass WDs can be lower than estimated in previous work to reach solar values of [Mn/Fe] at [Fe/H] = 0. Our galactic chemical evolution models suggest that Type Ia supernovae from sub-Chandrasekhar mass white dwarfs, along with core-collapse supernovae, could account for more than 80% of the solar Mn abundance. Using metallicity-dependent Type Ia supernova yields helps to reproduce the upward trend of [Mn/Fe] as a function of metallicity for the solar neighborhood. These chemical evolution predictions, however, depend on the massive star yields adopted in the calculations.

Metallicity-suppressed collapsars cannot be the dominant r-process source in the milky way

ABSTRACT We develop a high-performance analytical model of Galactic Chemical Evolution, which accounts for delay time distributions and lock-up of stellar yields in a thermal-phased ISM. The model is capable of searching, for the first time, through the high-dimensional parameter space associated with the r-process enrichment of the Milky Way by its possible sources: Neutron Star Mergers and Collapsar events. Their differing formation mechanisms give these two processes different time dependencies, a property which has frequently been used to argue in favour of collapsars as the dominant r-process source. However, we show that even with large degrees of freedom in the allowed thermal, structural, and chemical properties of the galaxy, large regions of parameter space are in strong tension with the data. In particular, whilst we are able to find models in which neutron star mergers produce the majority of r-process material, the data rule out all models with dominant collapsar yields. With no other identified source, we conclude that Neutron Star Mergers must be the dominant contributors to the modern Milky Way r-process budget.

Chemical abundances of 1111 FGK stars from the HARPS GTO planet search program

Context. To understand the formation and composition of planetary systems, it is essential to have insights into the chemical composition of their host stars. In particular, C/O elemental ratios are useful for constraining the density and bulk composition of terrestrial planets. Aims. We study the carbon abundances with a twofold objective. On the one hand, we want to evaluate the behaviour of carbon in the context of Galactic chemical evolution. On the other hand, we focus on the possible dependence of carbon abundances on the presence of planets and on the impact of various factors (such as different oxygen lines) on the determination of C/O elemental ratios. Methods. We derived chemical abundances of carbon from two atomic lines for 757 FGK stars in the HARPS-GTO sample, observed with high-resolution (R ~ 115 000) and high-quality spectra. The abundances were derived using a standard Local Thermodynamic Equilibrium analysis with automatically measured Equivalent Widths injected into the code MOOG and a grid of Kurucz ATLAS9 atmospheres. Oxygen abundances, derived using different lines, were taken from previous papers in this series and updated with the new stellar parameters. Results. We find that thick- and thin-disk stars are chemically disjunct for [C/Fe] across the full metallicity range that they have in common. Moreover, the population of high-α metal-rich stars also presents higher and clearly separated [C/Fe] ratios than thin-disk stars up to [Fe/H] ~ 0.2 dex. The [C/O] ratios present a general flat trend as a function of [O/H] but becomes negative at [O/H] ≳ 0dex. This trend is more clear when considering stars of similar metallicity. We find tentative evidence that stars with low-mass planets at lower metallicities have higher [C/Fe] ratios than stars without planets at the same metallicity, in the same way as has previously been found for α elements. Finally, the elemental C/O ratios for the vast majority of our stars are below 0.8 when using the oxygen line at 6158 Å, however, the forbidden oxygen line at 6300 Å provides systematically higher C/O values (going above 1.2 in a few cases) which also show a dependence on Teff. Moreover, by using different atmosphere models the C/O ratios can have a non-negligible difference for cool stars. Therefore, C/O ratios should be scaled to a common solar reference in order to correctly evaluate its behaviour. We find no significant differences in the distribution of C/O ratios for the different populations of planet hosts, except when comparing the stars without detected planets with the stars hosting Jupiter-type planets. However, we note that this difference might be caused by the different metallicity distributions of both populations. Conclusions. The derivation of homogeneous abundances from high-resolution spectra in samples that are modest in size is of great utility in constraining models of Galactic chemical evolution. The combination of these high-quality data with the long-term study of planetary presence in our sample is crucial for achieving an accurate understanding of the impact of stellar chemical composition on planetary formation mechanisms.

Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars

Context. Next to H and He, carbon is, together with oxygen, the most abundant element in the Universe and widely used when modelling the formation and evolution of galaxies and their stellar populations. For the Milky Way bulge, there are currently essentially no measurements of carbon in un-evolved stars, hampering our abilities to properly compare Galactic chemical evolution models to observational data for this still enigmatic stellar population. Aims. We aim to determine carbon abundances for our sample of 91 microlensed dwarf and subgiant stars in the Galactic bulge. Together with new determinations for oxygen this forms the first statistically significant sample of bulge stars that have C and O abundances measured, and for which the C abundances have not been altered by the nuclear burning processes internal to the stars. Methods. Our analysis is based on high-resolution spectra for a sample of 91 dwarf and subgiant stars that were obtained during microlensing events when the brightnesses of the stars were highly magnified. Carbon abundances were determined through spectral line synthesis of six C I lines around 9100 Å, and oxygen abundances using the three O I lines at about 7770 Å. One-dimensional (1D) MARCS model stellar atmospheres calculated under the assumption of local thermodynamic equilibrium (LTE) were used, and non-LTE corrections were applied when calculating the synthetic spectra for both C and O. Results. Carbon abundances was possible to determine for 70 of the 91 stars in the sample and oxygen abundances for 88 of the 91 stars in the sample. The [C/Fe] ratio evolves essentially in lockstep with [Fe/H], centred around solar values at all [Fe/H]. The [O/Fe]–[Fe/H] trend has an appearance very similar to that observed for other α-elements in the bulge, with the exception of a continued decrease in [O/Fe] at super-solar [Fe/H], where other α-elements tend to level out. When dividing the bulge sample into two sub-groups, one younger than 8 Gyr and one older than 8 Gyr, the stars in the two groups follow exactly the elemental abundance trends defined by the solar neighbourhood thin and thick disks, respectively. Comparisons with recent models of Galactic chemical evolution in the [C/O]–[O/H] plane show that the models that best match the data are the ones that have been calculated with the Galactic thin and thick disks in mind. Conclusions. We conclude that carbon, oxygen, and the combination of the two support the idea that the majority of the stars in the Galactic bulge have a secular origin; that is, they are formed from disk material. We cannot exclude that a fraction of stars in the bulge could be classified as a classical bulge population, but it would have to be small. More dedicated and advanced models of the inner region of the Milky Way are needed to make more detailed comparisons to the observations.

Estimation of ages and masses via carbon and nitrogen abundances for 556 007 giants from LAMOST

Estimating ages for a large number of giants is of great importance for studying Galactic evolution. In this work, we determine stellar ages and masses for 556 007 giants from LAMOST Data Release 5 (DR5) with empirical relations estimated from chemical [C/N] abundance ratios. Our sample reveals the two well-known sequences in the age–[α/M] relation. The high-α sequence is composed of stars older than 8 Gyr and low-α sequence is composed of stars with age ranging from 0 Gyr to 13.8 Gyr. Our sample also shows a flat age–[M/H] relation up until 12 Gyr. We compare these distributions with Galactic Chemical Evolution models for reference. When looking at the spatial distribution of stars in 2 Gyr age bins, we find that young stars are concentrated towards the Galactic plane and older stars extend to higher height above and below the disk. We find a smooth transition of median Galactic height for different age bins, which suggests a strong age-dependence on Galactic scale height.

Stellar migration and chemical enrichment in the milky way disc: a hybrid model

ABSTRACT We develop a hybrid model of galactic chemical evolution that combines a multiring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc’s abundance structure: (i) the dependence of the [O/Fe]–[Fe/H] distribution on radius Rgal and mid-plane distance |z|; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with Rgal and |z|; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with Rgal, |z|, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] > 0.1; (v) a population of young and intermediate-age α-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age–[O/H] and age–[Fe/H] relations, with large scatter and a median age of ∼4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate ∼2 Gyr ago improve agreement with the observed age–[Fe/H] and age–[O/H] relations, but worsen agreement with the observed age–[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice).