skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: 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).  more » « less
Award ID(s):
1909841
PAR ID:
10378456
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
508
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
4484 to 4511
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We investigate the impact of bursts in star formation on the predictions of one-zone chemical evolution models, adopting oxygen (O), iron (Fe), and strontium (Sr), as representative α, iron-peak, and s-process elements, respectively. To this end, we develop and make use of the Versatile Integrator for Chemical Evolution (VICE), a python package designed to handle flexible user-specified evolutionary parameters. Starbursts driven by a temporary boost of gas accretion rate create loops in [O/Fe]–[Fe/H] evolutionary tracks and a peak in the stellar [O/Fe] distribution at intermediate values. Bursts driven by a temporary boost of star formation efficiency have similar effects, and they also produce a population of α-deficient stars during the depressed star formation phase following the burst. This α-deficient population is more prominent if the outflow rate is tied to a time-averaged star formation rate (SFR) instead of the instantaneous SFR. Theoretical models of Sr production predict a strong metallicity dependence of supernova and asymptotic giant branch star yields, though comparison to data suggests an additional, nearly metallicity-independent source. Evolution of [Sr/Fe] and [Sr/O] during a starburst is complex because of this metallicity dependence and the multiple time-scales at play. Moderate amplitude (10–20 per cent) sinusoidal oscillations in SFR produce loops in [O/Fe]–[Fe/H] tracks and multiple peaks in [O/Fe] distributions, a potential source of intrinsic scatter in observed sequences. We investigate the impact of a factor ∼2 enhancement of Galactic star formation ∼2 Gyr ago, as suggested by some recent observations. VICE is publicly available at <http://pypi.org/project/vice/>. 
    more » « less
  2. ABSTRACT We derive empirical constraints on the nucleosynthetic yields of nitrogen by incorporating N enrichment into our previously developed and empirically tuned multizone galactic chemical evolution model. We adopt a metallicity-independent (‘primary’) N yield from massive stars and a metallicity-dependent (‘secondary’) N yield from AGB stars. In our model, galactic radial zones do not evolve along the observed [N/O]–[O/H] relation, but first increase in [O/H] at roughly constant [N/O], then move upward in [N/O] via secondary N production. By t ≈ 5 Gyr, the model approaches an equilibrium [N/O]–[O/H] relation, which traces the radial oxygen gradient. Reproducing the [N/O]–[O/H] trend observed in extragalactic systems constrains the ratio of IMF-averaged N yields to the IMF-averaged O yield of core-collapse supernovae. We find good agreement if we adopt $$y_\text{N}^\text{CC}/y_\text{O}^\text{CC}=0.024$$ and $$y_\text{N}^\text{AGB}/y_\text{O}^\text{CC} = 0.062(Z/Z_\odot)$$. For the theoretical AGB yields we consider, simple stellar populations release half their N after only ∼250 Myr. Our model reproduces the [N/O]–[O/H] relation found for Milky Way stars in the APOGEE survey, and it reproduces (though imperfectly) the trends of stellar [N/O] with age and [O/Fe]. The metallicity-dependent yield plays the dominant role in shaping the gas-phase [N/O]–[O/H] relation, but the AGB time-delay is required to match the stellar age and [O/Fe] trends. If we add ∼40 per cent oscillations to the star formation rate, the model reproduces the scatter in the gas phase [N/O]–[O/H] relation observed in external galaxies by MaNGA. We discuss implications of our results for theoretical models of N production by massive stars and AGB stars. 
    more » « less
  3. ABSTRACT Stellar ages are a crucial component to studying the evolution of the Milky Way. Using Gaia DR2 distance estimates, it is now possible to estimate stellar ages for a larger volume of evolved stars through isochrone matching. This work presents [M/H]–age and [α/M]–age relations derived for different spatial locations in the Milky Way disc. These relations are derived by hierarchically modelling the star formation history of stars within a given chemical abundance bin. For the first time, we directly observe that significant variation is apparent in the [M/H]–age relation as a function of both Galactocentric radius and distance from the disc mid-plane. The [M/H]–age relations support claims that radial migration has a significant effect in the plane of the disc. Using the [M/H] bin with the youngest mean age at each radial zone in the plane of the disc, the present-day metallicity gradient is measured to be −0.059 ± 0.010 dex kpc−1, in agreement with Cepheids and young field stars. We find a vertically flared distribution of young stars in the outer disc, confirming predictions of models and previous observations. The mean age of the [M/H]–[α/M] distribution of the solar neighbourhood suggests that the high-[M/H] stars are not an evolutionary extension of the low-α sequence. Our observational results are important constraints to Galactic simulations and models of chemical evolution. 
    more » « less
  4. Abstract We present new maps of the Milky Way disk showing the distribution of metallicity ([Fe/H]),α-element abundances ([Mg/Fe]), and stellar age, using a sample of 66,496 red giant stars from the final data release (DR17) of the Apache Point Observatory Galactic Evolution Experiment survey. We measure radial and vertical gradients, quantify the distribution functions for age and metallicity, and explore chemical clock relations across the Milky Way for the low-αdisk, high-αdisk, and total population independently. The low-αdisk exhibits a negative radial metallicity gradient of −0.06 ± 0.001 dex kpc−1, which flattens with distance from the midplane. The high-αdisk shows a flat radial gradient in metallicity and age across nearly all locations of the disk. The age and metallicity distribution functions shift from negatively skewed in the inner Galaxy to positively skewed at large radius. Significant bimodality in the [Mg/Fe]–[Fe/H] plane and in the [Mg/Fe]–age relation persist across the entire disk. The age estimates have typical uncertainties of ∼0.15 in log(age) and may be subject to additional systematic errors, which impose limitations on conclusions drawn from this sample. Nevertheless, these results act as critical constraints on galactic evolution models, constraining which physical processes played a dominant role in the formation of the Milky Way disk. We discuss how radial migration predicts many of the observed trends near the solar neighborhood and in the outer disk, but an additional more dramatic evolution history, such as the multi-infall model or a merger event, is needed to explain the chemical and age bimodality elsewhere in the Galaxy. 
    more » « less
  5. ABSTRACT Stellar radial migration plays an important role in reshaping a galaxy’s structure and the radial distribution of stellar population properties. In this work, we revisit reported observational evidence for radial migration and quantify its strength using the age–[Fe/H] distribution of stars across the Milky Way with APOGEE data. We find a broken age–[Fe/H] relation in the Galactic disc at r > 6 kpc, with a more pronounced break at larger radii. To quantify the strength of radial migration, we assume stars born at each radius have a unique age and metallicity, and then decompose the metallicity distribution function (MDF) of mono-age young populations into different Gaussian components that originated from various birth radii at rbirth < 13 kpc. We find that, at ages of 2 and 3 Gyr, roughly half the stars were formed within 1 kpc of their present radius, and very few stars (<5 per cent) were formed more than 4 kpc away from their present radius. These results suggest limited short-distance radial migration and inefficient long-distance migration in the Milky Way during the last 3 Gyr. In the very outer disc beyond 15 kpc, the observed age–[Fe/H] distribution is consistent with the prediction of pure radial migration from smaller radii, suggesting a migration origin of the very outer disc. We also estimate intrinsic metallicity gradients at ages of 2 and 3 Gyr of −0.061 and −0.063 dex kpc−1, respectively. 
    more » « less