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: The Rapid Formation of the Metal-poor Milky Way
Abstract Our understanding of the assembly timeline of the Milky Way has been transforming along with the dramatic increase in astrometric and spectroscopic data available over the past several years. Many substructures in chemo-dynamical space have been discovered and identified as the remnants of various galactic mergers. To investigate the timeline of these mergers, we select main-sequence turnoff and subgiant stars (MSTOs) from the H3 survey, finding members in seven metal-poor components of the halo: Gaia-Sausage/Enceladus (GSE), the Helmi Streams, Thamnos, Sequoia, Wukong/LMS-1, Arjuna, and I’itoi. We also select out a metal-poor in situ population to facilitate comparison to the evolution of the Milky Way itself at these early epochs. We fit individual isochrone ages to the MSTOs in each of these substructures and use the resulting age distributions to infer simple star formation histories (SFHs). For GSE, we resolve an extended SFH that truncates ≈10 Gyr ago, as well as a clear age–metallicity relation. From this age distribution and measured SFH we infer that GSE merged with the Milky Way at a time 9.5–10.2 Gyr ago, in agreement with previous estimates. We infer that the other mergers occurred at various times ranging from 9 to 13 Gyr ago, and that the metal-poor in situ Galaxy built up within only a few billion years. These results reinforce the emerging picture that both the disk and halo of the Milky Way experienced a rapid assembly.  more » « less
Award ID(s):
2406729
PAR ID:
10612804
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
The Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
978
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
152
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The star formation histories (SFHs) of galactic stellar haloes offer crucial insights into the merger history of the galaxy and the effects of those mergers on their hosts. Such measurements have revealed that while the Milky Way’s most important merger was 8–10 Gyr ago, M31’s largest merger was more recent, within the last few Gyr. Unfortunately, the required halo SFH measurements are extremely observationally expensive outside of the Local Group. Here, we use asymptotic giant branch (AGB) stars brighter than the tip of the red giant branch (RGB) to constrain stellar halo SFHs. Both stellar population models and archival data sets show that the AGB/RGB ratio constrains the time before which 90 per cent of the stars formed, t90. We find AGB stars in the haloes of three highly inclined roughly Milky Way-mass galaxies with resolved star measurements from the Hubble Space Telescope; this population is most prominent in the stellar haloes of NGC 253 and NGC 891, suggesting that their stellar haloes contain stars born at relatively late times, with inferred t90 ∼ 6 ± 1.5 Gyr. This ratio also varies from region to region, tending towards higher values along the major axis and in tidal streams or shells. By combining our measurements with previous constraints, we find a tentative anticorrelation between halo age and stellar halo mass, a trend that exists in models of galaxy formation but has never been elucidated before, i.e. the largest stellar haloes of Milky Way-mass galaxies were assembled more recently. 
    more » « less
  2. The vast majority of Milky Way stellar halo stars were likely accreted from a small number (<~3) of relatively large dwarf galaxy accretion events. However, the timing of these events is poorly constrained and predominantly relies on indirect dynamical mixing arguments or imprecise age measurements of stars associated with debris structures. Here, we aim to infer robust stellar ages for stars associated with galactic substructures to more directly constrain the merger history of the Galaxy. By combining kinematic, asteroseismic, and spectroscopic data where available, we infer stellar ages for a sample of 10 red giant stars that were kinematically selected to be within the stellar halo, a subset of which are associated with the Gaia–Enceladus–Sausage halo substructure, and compare their ages to 3 red giant stars in the Galactic disk. Despite systematic differences in both absolute and relative ages determined here, age rankings of stars in this sample are robust. Passing the same observable inputs to multiple stellar age determination packages, we measure a weighted average age for the Gaia–Enceladus–Sausage stars in our sample of 8+/-3 (stat.)+/-1 (sys.) Gyr. We also determine hierarchical ages using isochrones for the populations of Gaia–Enceladus–Sausage, in situ halo and disk stars, finding a Gaia–Enceladus–Sausage population age of 8.0+2.3-3.2 Gyr. Although we cannot distinguish hierarchical population ages of halo or disk structures with our limited data and sample of stars, this framework should allow a distinct characterization of Galactic substructures using larger stellar samples and additional data available in the near future 
    more » « less
  3. Abstract The majority of the Milky Way’s stellar halo consists of debris from our galaxy’s last major merger, the Gaia-Sausage-Enceladus (GSE). In the past few years, stars from the GSE have been kinematically and chemically studied in the inner 30 kpc of our galaxy. However, simulations predict that accreted debris could lie at greater distances, forming substructures in the outer halo. Here we derive metallicities and distances using Gaia DR3 XP spectra for an all-sky sample of luminous red giant stars, and map the outer halo with kinematics and metallicities out to 100 kpc. We obtain follow-up spectra of stars in two strong overdensities—including the previously identified outer Virgo Overdensity—and find them to be relatively metal rich and on predominantly retrograde orbits, matching predictions from simulations of the GSE merger. We argue that these are apocentric shells of GSE debris, forming 60–90 kpc counterparts to the 15–20 kpc shells that are known to dominate the inner stellar halo. Extending our search across the sky with literature radial velocities, we find evidence for a coherent stream of retrograde stars encircling the Milky Way from 50 to 100 kpc, in the same plane as the Sagittarius Stream but moving in the opposite direction. These are the first discoveries of distant and structured imprints from the GSE merger, cementing the picture of an inclined and retrograde collision that built up our galaxy’s stellar halo. 
    more » « less
  4. Abstract Several lines of evidence suggest that the Milky Way underwent a major merger at z ∼ 2 with the Gaia-Sausage-Enceladus (GSE) galaxy. Here we use H3 Survey data to argue that GSE entered the Galaxy on a retrograde orbit based on a population of highly retrograde stars with chemistry similar to the largely radial GSE debris. We present the first tailored N -body simulations of the merger. From a grid of ≈500 simulations we find that a GSE with M ⋆ = 5 × 10 8 M ⊙ , M DM = 2 × 10 11 M ⊙ best matches the H3 data. This simulation shows that the retrograde stars are stripped from GSE’s outer disk early in the merger. Despite being selected purely on angular momenta and radial distributions, this simulation reproduces and explains the following phenomena: (i) the triaxial shape of the inner halo, whose major axis is at ≈35° to the plane and connects GSE’s apocenters; (ii) the Hercules-Aquila Cloud and the Virgo Overdensity, which arise due to apocenter pileup; and (iii) the 2 Gyr lag between the quenching of GSE and the truncation of the age distribution of the in situ halo, which tracks the lag between the first and final GSE pericenters. We make the following predictions: (i) the inner halo has a “double-break” density profile with breaks at both ≈15–18 kpc and 30 kpc, coincident with the GSE apocenters; and (ii) the outer halo has retrograde streams awaiting discovery at >30 kpc that contain ≈10% of GSE’s stars. The retrograde (radial) GSE debris originates from its outer (inner) disk—exploiting this trend, we reconstruct the stellar metallicity gradient of GSE (−0.04 ± 0.01 dex r 50 − 1 ). These simulations imply that GSE delivered ≈20% of the Milky Way’s present-day dark matter and ≈50% of its stellar halo. 
    more » « less
  5. ABSTRACT We model the stellar abundances and ages of two disrupted dwarf galaxies in the Milky Way stellar halo: Gaia-Sausage Enceladus (GSE) and Wukong/LMS-1. Using a statistically robust likelihood function, we fit one-zone models of galactic chemical evolution with exponential infall histories to both systems, deriving e-folding time-scales of τin = 1.01 ± 0.13 Gyr for GSE and $$\tau _\text{in} = 3.08^{+3.19}_{-1.16}$$ Gyr for Wukong/LMS-1. GSE formed stars for $$\tau _\text{tot} = 5.40^{+0.32}_{-0.31}$$ Gyr, sustaining star formation for ∼1.5–2 Gyr after its first infall into the Milky Way ∼10 Gyr ago. Our fit suggests that star formation lasted for $$\tau _\text{tot} = 3.36^{+0.55}_{-0.47}$$ Gyr in Wukong/LMS-1, though our sample does not contain any age measurements. The differences in evolutionary parameters between the two are qualitatively consistent with trends with stellar mass M⋆ predicted by simulations and semi-analytic models of galaxy formation. Our inferred values of the outflow mass-loading factor reasonably match $$\eta \propto M_\star ^{-1/3}$$ as predicted by galactic wind models. Our fitting method is based only on Poisson sampling from an evolutionary track and requires no binning of the data. We demonstrate its accuracy by testing against mock data, showing that it accurately recovers the input model across a broad range of sample sizes (20 ≤ N ≤ 2000) and measurement uncertainties (0.01 ≤ σ[α/Fe], σ[Fe/H] ≤ 0.5; $$0.02 \le \sigma _{\log _{10}(\text{age})} \le 1$$). Due to the generic nature of our derivation, this likelihood function should be applicable to one-zone models of any parametrization and easily extensible to other astrophysical models which predict tracks in some observed space. 
    more » « less