skip to main content

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] more » 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). « less
Authors:
; ; ; ; ; ; ; ;
Award ID(s):
1909841
Publication Date:
NSF-PAR ID:
10378456
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
508
Issue:
3
Page Range or eLocation-ID:
4484 to 4511
ISSN:
0035-8711
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 andmore »[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.

    « less
  2. 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 loopsmore »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/>.« 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 extensionmore »of the low-α sequence. Our observational results are important constraints to Galactic simulations and models of chemical evolution.« less
  4. ABSTRACT Using a sample of red giant stars from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16, we infer the conditional distribution $p([\alpha /{\rm Fe}]\, |\, [{\rm Fe}/{\rm H}])$ in the Milky Way disk for the α-elements Mg, O, Si, S, and Ca. In each bin of [Fe/H] and Galactocentric radius R, we model p([α/Fe]) as a sum of two Gaussians, representing ‘low-α’ and ‘high-α’ populations with scale heights $z_1=0.45\, {\rm kpc}$ and $z_2=0.95\, {\rm kpc}$, respectively. By accounting for age-dependent and z-dependent selection effects in APOGEE, we infer the [α/Fe] distributions that would be found for a fair sample of long-lived stars covering all z. Near the Solar circle, this distribution is bimodal at sub-solar [Fe/H], with the low-α and high-α peaks clearly separated by a minimum at intermediate [α/Fe]. In agreement with previous results, we find that the high-α population is more prominent at smaller R, lower [Fe/H], and larger |z|, and that the sequence separation is smaller for Si and Ca than for Mg, O, and S. We find significant intrinsic scatter in [α/Fe] at fixed [Fe/H] for both the low-α and high-α populations, typically ∼0.04-dex. The means, dispersions, and relative amplitudes of thismore »two-Gaussian description, and the dependence of these parameters on R, [Fe/H], and α-element, provide a quantitative target for chemical evolution models and a test for hydrodynamic simulations of disk galaxy formation. We argue that explaining the observed bimodality will probably require one or more sharp transitions in the disk’s gas accretion, star formation, or outflow history in addition to radial mixing of stellar populations.« less
  5. Abstract Type Ia supernovae are critical for feedback and elemental enrichment in galaxies. Recent surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific supernova Ia rate at z ∼ 0 may be ≲ 20 − 50 × higher in lower-mass galaxies than at Milky Way-mass. Independently, observations show that the close-binary fraction of solar-type Milky Way stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of metallicity-dependent rate models on galaxies of $M_* \sim 10^7\, \rm {M}_{\odot }-10^{11}\, \rm {M}_{\odot }$. First, we benchmark our simulated star-formation histories (SFHs) against observations, and show that assumed stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent rate models, potentially causing ASAS-SN and DES observations to agree more than might appear. Models in which the supernova Ia rate increases with decreasing metallicity ($\propto Z^{-0.5 \; \rm {to} \; -1}$) provide significantly better agreement with observations. Encouragingly, these rate increases (≳ 10 × in low-mass galaxies) do not significantly impact galaxy masses and morphologies, which remain largely unaffected except for our most extreme models.more »We explore implications for both [Fe/H] and [$\alpha /\rm {Fe}$] enrichment; metallicity-dependent rate models can improve agreement with the observed stellar mass-metallicity relations in low-mass galaxies. Our results demonstrate that a range of metallicity-dependent rate models are viable for galaxy formation and motivate future work.« less