A major goal in the field of galaxy formation is to understand the formation of the Milky Way’s disk. The first step toward doing this is to empirically describe its present state. We use the new high-dimensional data set of 19 abundances from 27,135 red clump Apache Point Observatory Galactic Evolution Experiment stars to examine the distribution of clusters defined using abundances. We explore different dimension reduction techniques and implement a nonparametric agglomerate hierarchical clustering method. We see that groups defined using abundances are spatially separated, as a function of age. Furthermore, the abundance groups represent different distributions in the [Fe/H]–age plane. Ordering our clusters by age reveals patterns suggestive of the sequence of chemical enrichment in the disk over time. Our results indicate that a promising avenue to trace the details of the disk’s assembly is via a full interpretation of the empirical connections we report.

To understand the formation and evolution of the Milky Way disk, we must connect its current properties to its past. We explore hydrodynamical cosmological simulations to investigate how the chemical abundances of stars might be linked to their origins. Using hierarchical clustering of abundance measurements in two Milky Way–like simulations with distributed and steady star formation histories, we find that groups of chemically similar stars comprise different groups in birth place (R birth) and time (age). Simulating observational abundance errors (0.05 dex), we find that to trace distinct groups of (R birth, age) requires a large vector of abundances. Using 15 element abundances (Fe, O, Mg, S, Si, C, P, Mn, Ne, Al, N, V, Ba, Cr, Co), up to ≈10 groups can be defined with ≈25% overlap in (R birth, age). We build a simple model to show that in the context of these simulations, it is possible to infer a star’s age and R birth from abundances with precisions of ±0.06 Gyr and ±1.17 kpc, respectively. We find that abundance clustering is ineffective for a third simulation, where low-α stars form distributed in the disk and early high-α stars form more rapidly in clumps that sink toward the Galactic center as their constituent stars evolve to enrich the interstellar medium. However, this formation path leads to large age dispersions across the [α/Fe]–[Fe/H] plane, which is inconsistent with the Milky Way’s observed properties. We conclude that abundance clustering is a promising approach toward charting the history of our Galaxy.

Chemical abundances of Milky Way disk stars are empirical tracers of its enrichment history. However, they capture joint-information that is valuable to disentangle. In this work, we quantify how individual abundances evolve across the present-day Galactic radius, at fixed supernovae contribution ([Fe/H], [Mg/Fe]). We use 18,135 Apache Point Observatory Galactic Evolution Experiment Data Release 17 red clump stars and 7943 GALactic Archaeology with HERMES Data Release 3 main-sequence stars to compare the abundance distributions conditioned on ([Fe/H], [Mg/Fe]) across 3–13 kpc and 6.5–9.5 kpc, respectively. We examine 15 elements: C, N, Al, K (light), O, Si, S, Ca, (α), Mn, Ni, Cr, Cu, (iron-peak) Ce, Ba (s-process) and Eu (r-process). We find that the conditional neutron-capture and light elements most significantly trace variations in the disk’s enrichment history, with absolute conditional radial gradients ≤0.03 dex kpc−1. The other elements studied have absolute conditional gradients ≲0.01 dex kpc−1. We uncover structured conditional abundance variations with [Fe/H] for the low-α, but not the high-α , sequence. The average scatter between the mean conditional abundances at different radii is σ intrinsic ≈ 0.02 dex (Ce, Eu, Ba σ intrinsic > 0.05 dex). These results serve as a measure of the magnitude via which different elements trace Galactic radial enrichment history once fiducial supernovae correlations are accounted for. Furthermore, we uncover subtle systematic variations in moments of the conditional abundance distributions and bimodal differences in [Al/Fe]. These suggest a nonuniform enrichment of each chemical cell, and will presumably constrain chemical evolution models of the Galaxy.

Chemical abundances are an essential tool in untangling the Milky Way’s enrichment history. However, the evolution of the interstellar medium abundance gradient with cosmic time is lost as a result of radial mixing processes. For the first time, we quantify the evolution of many observational abundances across the Galactic disk as a function of lookback time and birth radius, $\rm \text{R}_\text{birth}$. Using an empirical approach, we derive $\rm \text{R}_\text{birth}$ estimates for 145,447 APOGEE DR17 red giant disk stars, based solely on their ages and $\rm [Fe/H]$. We explore the detailed evolution of 6 abundances (Mg, Ca (α), Mn (iron-peak), Al, C (light), Ce (s-process)) across the Milky Way disk using 87,426 APOGEE DR17 red giant stars. We discover that the interstellar medium had three fluctuations in the metallicity gradient ∼9, ∼6, and ∼4 Gyr ago. The first coincides with the end of high-α sequence formation around the time of the Gaia-Sausage-Enceladus disruption, while the others are likely related to passages of the Sagittarius dwarf galaxy. A clear distinction is found between present-day observed radial gradients with age and the evolution with lookback time for both [X/Fe] and [X/H], resulting from the significant flattening and inversion in old populations due to radial migration. We find the $\rm [Fe/H]$–$\rm [\alpha /Fe]$ bimodality is also seen as a separation in the $\rm \text{R}_\text{birth}$–$\rm [X/Fe]$ plane for the light and α-elements. Our results recover the chemical enrichment of the Galactic disk over the past 12 Gyr, providing tight constraints on Galactic disk chemical evolution models.

Recent observations of the Milky Way (MW) found an unexpected steepening of the star-forming gas metallicity gradient around the time of the Gaia-Sausage-Enceladus (GSE) merger event. Here we investigate the influence of early ($t_{\rm {merger}}\lesssim 5$ Gyr) massive ($M_{\rm {gas}}^{\rm {merger}}/M_{\rm {gas}}^{\rm {main}}(t_{\rm {merger}})\gtrsim 10~{{\%}}$) merger events such as the Gaia-Sausage Enceladus merger in the MW on the evolution of the cold gas metallicity gradient. We use the NIHAO-UHD suite of cosmological hydrodynamical simulations of MW-mass galaxies to study the frequency of massive early mergers and their detailed impact on the morphology and chemistry of the gaseous disks. We find a strong steepening of the metallicity gradient at early times for all four galaxies in our sample which is caused by a sudden increase in the cold gas disk size (up to a factor of 2) in combination with the supply of un-enriched gas (∼0.75 dex lower compared to the main galaxy) by the merging dwarf galaxies. The mergers mostly affect the galaxy outskirts and lead to an increase in cold gas surface density of up to 200% outside of ∼8 kpc. The addition of un-enriched gas breaks the self-similar enrichment of the inter-stellar-medium and causes a dilution of the cold gas in the outskirts of the galaxies. The accreted stars and the ones formed later out of the accreted gas inhabit distinct tracks offset to lower [α/Fe] and [Fe/H] values compared to the main galaxy’s stars. We find that such mergers can contribute significantly to the formation of a second, low-α sequence as is observed in the MW.