Anomalously high nitrogen-to-oxygen abundance ratios [N/O] are observed in globular clusters (GCs), among the field stars of the Milky Way (MW), and even in the gas in a z ≈ 11 galaxy. Using data from the APOGEE Data Release 17 and the Gaia Data Release 3, we present several independent lines of evidence that most of the MW’s high-[N/O] stars were born in situ in massive bound clusters during the early, pre-disc evolution of the Galaxy. Specifically, we show that distributions of metallicity [Fe/H], energy, the angular momentum Lz, and distance of the low-metallicity high-[N/O] stars match the corresponding distributions of stars of the Aurora population and of the in situ GCs. We also show that the fraction of in situ field high-[N/O] stars, fN/O, increases rapidly with decreasing metallicity. During epochs when metallicity evolves from $\rm [Fe/H]=-1.5$ to $\rm [Fe/H]=-0.9$, the Galaxy spins up and transitions from a turbulent Aurora state to a coherently rotating disc. This transformation is accompanied by many qualitative changes. In particular, we show that high N/O abundances similar to those observed in GN-z11 were common before the spin-up ($\rm [Fe/H]\lesssim -1.5$) when up to $\approx 50~{{\ \rm per\ cent}}-70~{{\ \rm per\ cent}}$ of the in situ stars formed in massive bound clusters. The dramatic drop of fN/O at $\rm [Fe/H]\gtrsim -0.9$ indicates that after the disc emerges the fraction of stars forming in massive bound clusters decreases by two orders of magnitude.
We use accurate estimates of aluminium abundance from the APOGEE Data Release 17 and Gaia Early Data Release 3 astrometry to select a highly pure sample of stars with metallicity −1.5 ≲ [Fe/H] ≲ 0.5 born in-situ in the Milky Way proper. The low-metallicity ([Fe/H] ≲ −1.3) in-situ component we dub Aurora is kinematically hot with an approximately isotropic velocity ellipsoid and a modest net rotation. Aurora stars exhibit large scatter in metallicity and in many element abundance ratios. The median tangential velocity of the in-situ stars increases sharply with metallicity between [Fe/H] = −1.3 and −0.9, the transition that we call the spin-up. The observed and theoretically expected age–metallicity correlations imply that this increase reflects a rapid formation of the MW disc over ≈1–2 Gyr. The transformation of the stellar kinematics as a function of [Fe/H] is accompanied by a qualitative change in chemical abundances: the scatter drops sharply once the Galaxy builds up a disc during later epochs corresponding to [Fe/H] > −0.9. Results of galaxy formation models presented in this and other recent studies strongly indicate that the trends observed in the MW reflect generic processes during the early evolution of progenitors of MW-sized galaxies: a period of chaotic pre-disc evolution, when gas is accreted along cold narrow filaments and when stars are born in irregular configurations, and subsequent rapid disc formation. The latter signals formation of a stable hot gaseous halo around the MW progenitor, which changes the mode of gas accretion and allows development of coherently rotating disc.
more » « less- Award ID(s):
- 1911111
- NSF-PAR ID:
- 10412781
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 514
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 689-714
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT We characterize the 3D spatial variations of [Fe/H], [Mg/H], and [Mg/Fe] in stars at the time of their formation, across 11 simulated Milky Way (MW)- and M31-mass galaxies in the FIRE-2 simulations, to inform initial conditions for chemical tagging. The overall scatter in [Fe/H] within a galaxy decreased with time until $\approx 7 \, \rm {Gyr}$ ago, after which it increased to today: this arises from a competition between a reduction of azimuthal scatter and a steepening of the radial gradient in abundance over time. The radial gradient is generally negative, and it steepened over time from an initially flat gradient $\gtrsim 12 \, \rm {Gyr}$ ago. The strength of the present-day abundance gradient does not correlate with when the disc ‘settled’; instead, it best correlates with the radial velocity dispersion within the galaxy. The strength of azimuthal variation is nearly independent of radius, and the 360 deg scatter decreased over time, from $\lesssim 0.17 \, \rm {dex}$ at $t_{\rm lb} = 11.6 \, \rm {Gyr}$ to $\sim 0.04 \, \rm {dex}$ at present-day. Consequently, stars at $t_{\rm lb} \gtrsim 8 \, \rm {Gyr}$ formed in a disc with primarily azimuthal scatter in abundances. All stars formed in a vertically homogeneous disc, Δ[Fe/H]$\le 0.02 \, \rm {dex}$ within $1 \, \rm {kpc}$ of the galactic mid-plane, with the exception of the young stars in the inner $\approx 4 \, \rm {kpc}$ at z ∼ 0. These results generally agree with our previous analysis of gas-phase elemental abundances, which reinforces the importance of cosmological disc evolution and azimuthal scatter in the context of stellar chemical tagging. We provide analytic fits to our results for use in chemical-tagging analyses.
-
We investigate the inner regions of the Milky Way using data from APOGEE and Gaia EDR3. Our inner Galactic sample has more than 26 500 stars within | X Gal |< 5 kpc, | Y Gal |< 3.5 kpc, | Z Gal |< 1 kpc, and we also carry out the analysis for a foreground-cleaned subsample of 8000 stars that is more representative of the bulge–bar populations. These samples allow us to build chemo-dynamical maps of the stellar populations with vastly improved detail. The inner Galaxy shows an apparent chemical bimodality in key abundance ratios [ α /Fe], [C/N], and [Mn/O], which probe different enrichment timescales, suggesting a star formation gap (quenching) between the high- and low- α populations. Using a joint analysis of the distributions of kinematics, metallicities, mean orbital radius, and chemical abundances, we can characterize the different populations coexisting in the innermost regions of the Galaxy for the first time. The chemo-kinematic data dissected on an eccentricity–| Z | max plane reveal the chemical and kinematic signatures of the bar, the thin inner disc, and an inner thick disc, and a broad metallicity population with large velocity dispersion indicative of a pressure-supported component. The interplay between these different populations is mapped onto the different metallicity distributions seen in the eccentricity–| Z | max diagram consistently with the mean orbital radius and V ϕ distributions. A clear metallicity gradient as a function of | Z | max is also found, which is consistent with the spatial overlapping of different populations. Additionally, we find and chemically and kinematically characterize a group of counter-rotating stars that could be the result of a gas-rich merger event or just the result of clumpy star formation during the earliest phases of the early disc that migrated into the bulge. Finally, based on 6D information, we assign stars a probability value of being on a bar orbit and find that most of the stars with large bar orbit probabilities come from the innermost 3 kpc, with a broad dispersion of metallicity. Even stars with a high probability of belonging to the bar show chemical bimodality in the [ α /Fe] versus [Fe/H] diagram. This suggests bar trapping to be an efficient mechanism, explaining why stars on bar orbits do not show a significant, distinct chemical abundance ratio signature.more » « less
-
null (Ed.)ABSTRACT In hierarchical structure formation, metal-poor stars in and around the Milky Way (MW) originate primarily from mergers of lower mass galaxies. A common expectation is therefore that metal-poor stars should have isotropic, dispersion-dominated orbits that do not correlate strongly with the MW disc. However, recent observations of stars in the MW show that metal-poor ($\rm {[Fe/H]}\lesssim -2$) stars are preferentially on prograde orbits with respect to the disc. Using the Feedback In Realistic Environments 2 (FIRE-2) suite of cosmological zoom-in simulations of MW/M31-mass galaxies, we investigate the prevalence and origin of prograde metal-poor stars. Almost all (11 of 12) of our simulations have metal-poor stars on preferentially prograde orbits today and throughout most of their history: we thus predict that this is a generic feature of MW/M31-mass galaxies. The typical prograde-to-retrograde ratio is ∼2:1, which depends weakly on stellar metallicity at $\rm {[Fe/H]}\lesssim -1$. These trends predicted by our simulations agree well with MW observations. Prograde metal-poor stars originate largely from a single Large/Small Magellanic Cloud (LMC/SMC)-mass gas-rich merger $7\!-\!12.5\, \rm {Gyr}$ ago, which deposited existing metal-poor stars and significant gas on an orbital vector that sparked the formation of and/or shaped the orientation of a long-lived stellar disc, giving rise to a prograde bias for all low-metallicity stars. We find subdominant contributions from in situ stars formed in the host galaxy before this merger, and in some cases, additional massive mergers. We find few clear correlations between any properties of our MW/M31-mass galaxies at z = 0 and the degree of this prograde bias as a result of diverse merger scenarios.more » « less
-
more » « less
ABSTRACT The r-process-enhanced (RPE) stars provide fossil records of the assembly history of the Milky Way (MW) and the nucleosynthesis of the heaviest elements. Observations by the R-Process Alliance (RPA) and others have confirmed that many RPE stars are associated with chemo-dynamically tagged groups, which likely came from accreted dwarf galaxies of the MW. However, we do not know how RPE stars are formed. Here, we present the result of a cosmological zoom-in simulation of an MW-like galaxy with r-process enrichment, performed with the highest resolution in both time and mass. Thanks to this advancement, unlike previous simulations, we find that most highly RPE (r-II; [Eu/Fe] > +0.7) stars are formed in low-mass dwarf galaxies that have been enriched in r-process elements for [Fe/H] $\lt -2.5$, while those with higher metallicity are formed in situ, in locally enhanced gas clumps that were not necessarily members of dwarf galaxies. This result suggests that low-mass accreted dwarf galaxies are the main formation site of r-II stars with [Fe/H] $\, \lt -2.5$. We also find that most low-metallicity r-II stars exhibit halo-like kinematics. Some r-II stars formed in the same halo show low dispersions in [Fe/H] and somewhat larger dispersions of [Eu/Fe], similar to the observations. The fraction of simulated r-II stars is commensurate with observations from the RPA, and the distribution of the predicted [Eu/Fe] for halo r-II stars matches that observed. These results demonstrate that RPE stars can be valuable probes of the accretion of dwarf galaxies in the early stages of their formation.
