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
This content will become publicly available on July 1, 2026
Reconstructing the Milky Way chemical map with the galactic chemical evolution tool OMEGA+ from SDSS-MWM
Context.Although current observations indicate that there are two distinct sequences of disk stars in the [α/M] versus [M/H] parameter space, further complexity is evident in the chemical makeup of the Milky Way and consequently suggests a complicated evolutionary history. Aims.We developed two-infall galactic chemical evolution (GCE) models consistent with the Galactic chemical map. Methods.We obtained new GCE models simulating the chemical evolution of the Milky Way, as constrained by a golden sample of 394 000 stellar abundances of the Milky Way Mapper survey from data release 19 of SDSS-V. The separation between the chemical thin and thick disks was defined using [Mg/M]. We used the chemical evolution environmentOMEGA+combined with Levenberg-Marquardt (LM) and bootstrapping algorithms for the optimization and error estimation. We simulated the entire Galactic disk and considered six galactocentric regions, allowing for a more detailed analysis of the formation of the inner, middle, and outer Galaxy. We investigated the evolution ofα, odd-Z, and iron-peak elements, covering 15 species altogether. Results.The chemical thin and thick disks are separated by Mg observations, which the otherα-elements show similar trends with, while odd-Z species demonstrate different patterns as functions of metallicity. In the inner Galactic disk regions, the locus of the low-Mg sequence is gradually shifted toward higher metallicity, while the high-Mg phase is less populated. The best-fit GCE models show a well-defined peak in the rate of the infalling matter as a function of the Galactic age, confirming a merger event about 10 Gyr ago. We show that the timescale of gas accretion, the exact time of the second infall and the ratio between the surface mass densities associated with the second infall event and the formation event vary with the distance from the Galactic center. According to the models, the disk was assembled within a timescale of (0.32±0.02) Gyr during a primary formation phase, followed by an increasing accretion rate over a (0.55±0.06) Gyr-timescale and a relaxation phase that lasted (2.86±0.70) Gyr, with a second peak seen for the infall rate at (4.13±0.19) Gyr. Conclusions.Our best Galaxy evolution models are consistent with an inside-out formation scenario of the Milky Way disk and in agreement with the findings of recent chemodynamical simulations.
more »
« less
- Award ID(s):
- 1927130
- PAR ID:
- 10649112
- Publisher / Repository:
- EDP Sciences
- Date Published:
- Journal Name:
- Astronomy & Astrophysics
- Volume:
- 699
- ISSN:
- 0004-6361
- Page Range / eLocation ID:
- A293
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Spectroscopic studies of elliptical galaxies show that their stellar population ages, mean metallicity, andαenhancement traced by [Mg/Fe] all increase with galaxy stellar mass or velocity dispersion. We use one-zone galactic chemical evolution (GCE) models with a flexible star formation history (SFH) to model the age, [Mg/H], and [Mg/Fe] inferred from simple stellar population (SSP) fits to observed ellipticals atz∼ 0 andz∼ 0.7. We show that an SSP fit to the spectrum computed from a full GCE model gives ages and abundances close to the light-weighted, logarithmically averaged values of the composite stellar population, 〈age〉, 〈[Mg/H]〉, and 〈[Mg/Fe]〉. With supernova Mg and Fe yields fixed to values motivated by Milky Way stellar populations, we find that predicted 〈[Mg/H]〉–〈age〉 and 〈[Mg/Fe]〉–〈age〉 relations are surprisingly insensitive to SFH parameters: Older galaxies have higher 〈[Mg/Fe]〉, but the detailed form of the SFH has limited impact. The star formation efficiency (SFE) and outflow efficiency affect the early and late evolution of 〈[Mg/H]〉, respectively; explaining observed trends requires higher SFE and lower outflows in more massive galaxies. With core-collapse supernova yields calibrated to the plateau [Mg/Fe]cc≈ 0.45 observed in many Milky Way studies, our models underpredict the observed 〈[Mg/Fe]〉 ratios of ellipticals by 0.05–0.1 dex. Increasing the core-collapse yield ratio to [Mg/Fe]cc= 0.55 improves the agreement, though the models remain below the data. We discuss potential resolutions of this discrepancy, including the possibility that many ellipticals terminate their star formation with a self-enriching, terminating burst that reduces the light-weighted age and boosts 〈[Mg/Fe]〉.more » « less
-
Abstract Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan–Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the parent of the OC stream, we use astrometry from the early third data release of ESA’s Gaia mission (Gaia EDR3) and radial velocity information from the Sloan Digital Sky Survey (SDSS)-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to find up to 13 stars that are likely members of the OC stream. We use the APOGEE survey to study the chemical nature (for up to 10 stars) of the OC stream in theα(O, Mg, Ca, Si, Ti, and S), odd-Z(Al, K, and V), Fe-peak (Fe, Ni, Mn, Co, and Cr), and neutron-capture (Ce) elemental groups. We find that the stars that make up the OC stream are not consistent with a monometallic population and have a median metallicity of −1.92 dex with a dispersion of 0.28 dex. Our results also indicate that the α elements are depleted compared to the known Milky Way populations and that its [Mg/Al] abundance ratio is not consistent with second-generation stars from globular clusters. The detailed chemical pattern of these stars, namely the [α/Fe]–[Fe/H] plane and the metallicity distribution, indicates that the OC stream progenitor is very likely to be a dwarf spheroidal galaxy with a mass of ∼106M⊙.more » « less
-
Abstract The Milky Way (MW) stellar disk has both a thin and a thick component. The thin disk is composed mostly of younger stars (≲8 Gyr) with a lower abundance ofα-elements, while the thick disk contains predominantly older stars (≳8–12 Gyr) with a higherαabundance, giving rise to anα-bimodality most prominent at intermediate metallicities. A proposed explanation for the bimodality is an episode of clumpy star formation, where high-αstars form in massive clumps that appear in the first few billion years of the MW’s evolution, while low-αstars form throughout the disk and over a longer time span. To better understand the evolution of clumps, we track them and their constituent stars in two clumpy MW simulations that reproduce theα-abundance bimodality, one with 10% and the other with 20% supernova feedback efficiency. We investigate the paths that these clumps take in the chemical space ([O/Fe]–[Fe/H]) as well as their mass, star formation rate (SFR), formation location, lifetime, and merger history. The clumps in the simulation with lower feedback last longer on average, with several lasting hundreds of millions of years. Some of the clumps do not reach high-α, but the ones that do on average have a higher SFR, longer lifetime, greater mass, and form closer to the Galactic center than the ones that do not. Most clumps that reach high-αmerge with others and eventually spiral into the Galactic center, but shed stars along the way to form most of the thick-disk component.more » « less
-
Abstract Many nucleosynthetic channels create the elements, but two-parameter models characterized byαand Fe nonetheless predict stellar abundances in the Galactic disk to accuracies of 0.02–0.05 dex for most measured elements, near the level of current abundance uncertainties. It is difficult to make individual measurements more precise than this to investigate lower-amplitude nucleosynthetic effects, but population studies of mean abundance patterns can reveal more subtle abundance differences. Here, we look at the detailed abundances for 67,315 stars from the Apache Point Observatory Galactic Evolution Experiment (or APOGEE) Data Release 17, but in abundance residuals away from a best-fit two-parameter, data-driven nucleosynthetic model. We find that these residuals show complex structures with respect to age, guiding radius, and vertical action that are not random and are also not strongly correlated with sources of systematic error such as ,Teff, and radial velocity. The residual patterns, especially in Na, C+N, Mn, and Ce, trace kinematic structures in the Milky Way, such as the inner disk, thick disk, and flared outer disk. A principal component analysis suggests that most of the observed structure is low-dimensional and can be explained by a few eigenvectors. We find that some, but not all, of the effects in the low-αdisk can be explained by dilution with fresh gas, so that the abundance ratios resemble those of stars with higher metallicity. The patterns and maps we provide can be combined with accurate forward models of nucleosynthesis, star formation, and gas infall to provide a more detailed picture of star and element formation in different Milky Way components.more » « less
An official website of the United States government
