To unravel the assembly history of our Galaxy we need to learn how many stars were born when, from what material and on what orbits. This requires precise age determinations for a large sample of stars that extend to the oldest possible ages (around 14 Gyr)9,12. Subgiant stars, which are stars sustained by hydrogen shell fusion, can be unique tracers for such purposes, as they exist in the brief stellar evolutionary phase that permits the most precise and direct age determination, because their luminosity is a direct measure of their age. Moreover, the chemical element compositions determined from the spectra of their photosphere surfaces accurately reflect their birth material composition billions of years ago. This makes subgiants the best practical tracers of Galactic archaeology, even compared to main-sequence turn-off stars, whose surface abundances may be altered by atomic diffusion effects13. However, because of the short lifetime of their evolutionary phase, subgiant stars are relatively rare, and large surveys are essential to build a large sample of these objects with good spectra, which have not been available in the past.

With the recent data release (eDR3) of the Gaia mission14,15 and the recent data release (DR7) of the LAMOST spectroscopic survey16,17, we identify a set of approximately 250,000 subgiant stars based on their position in the effective temperatures (Teff)–absolute magnitude (MK) diagram (Fig. 1a). The ages (τ) of these subgiant stars are estimated by fitting to the Yonsei–Yale (YY) stellar isochrones18 with a Bayesian approach, which draws on the astrometric distances (parallaxes), apparent magnitudes (fluxes), spectroscopic chemical abundances ([Fe/H], [α/Fe] where α refers to α elements Mg, Si, Ca, Ti), Teff and MK. As summarized in Fig. 1b, the sample stars have a median relative age uncertainty of only 7.5% across the age range from 1.5 Gyr to the age of the Universe (13.8 Gyr; ref. 19). The lower age limit of our sample is inherent to our approach: younger and hence more luminous subgiants can be confused with a different stellar evolutionary phase, the horizontal branch phase for far older stars, which would cause serious sample contamination. This sample constitutes a 100-fold leap in sample size for stars with comparably precise and consistent age estimates20,21. In addition, it is a large sample that covers a large spatial volume across the Milky Way (Fig. 1c) and most of the pertinent range in age and in metallicity (1.5 Gyr < τ < 13.8 Gyr, and −2.5 < [Fe/H] < 0.4). The sample also has a straightforward spatial selection function that allows us to estimate the space density of the tracers. These ingredients enable an alternative view of the Milky Way’s assembly history, especially the early formation history.

Fig. 1: The subgiant star sample with precise ages.
figure 1

a, Illustration of the subgiant selection in the TeffMK diagram, shown for the solar metallicity bin of −0.1 < [Fe/H] < 0.1. In total, the subgiant sample contains 247,104 stars. The solid curves are isochrones from the YY stellar evolution models18 for solar metallicity ([Fe/H] = 0, [α/Fe] = 0) for ages of 1, 2, 3, 4, 6, 8, 10, 12, 14, 16, 18 and 20 Gyr, illustrating how stellar ages can be determined from the position in the TeffMK diagram if [Fe/H] is known. The two straight lines bracket the region within which we define our subgiant star sample. b, Distribution in the relative age precision as a function of age: the mode of this precision distribution is at 6% and the median at 7.5%. For the subsequent analysis we will only use stars with a relative age precision of less than 15% (horizontal dashed line). Histograms in the top and right are normalized to the peak value Nmax. c, Spatial distribution of our subgiant sample stars in the RZ plane of Galactic cylindrical coordinates. The full extent of the Galactocentric radius in the sample is 6 kpc R 14 kpc and that of the distance from the Galactic mid-plane is −5 kpc Z 6 kpc. The bulk of the sample (90%) covers 7.2  kpc R 10.4 kpc and −1.2  kpc Z 2 kpc, as illustrated by the dashed lines.

Our Galaxy’s stellar age–metallicity distribution

The photospheric metallicity of any subgiant star of age τ reflects the element composition of the gas from which it formed at the epoch τ Gyr ago. The overall distribution of these stellar metallicities at different epochs, p(τ, [Fe/H]), thus encodes the chemical enrichment history of our Milky Way galaxy. Figure 2a presents this distribution for our data. It shows that the age–metallicity distribution exhibits a number of prominent and distinct sequences, including at least two age-separated sequences with [Fe/H] > −1, and a sequence of exclusively old stars at low metallicity, [Fe/H] < −1. The density of p(τ, [Fe/H]) may change with stellar orbit or Galactocentric radius, in the range our sample covers (6–14 kpc; Fig. 1). Yet, the ‘morphology’ of the distribution varies only slightly, enabling us to focus on the radially averaged distribution p(τ, [Fe/H]) here.

Fig. 2: Stellar age–metallicity relation revealed by our subgiant star sample.
figure 2

a, Stellar distribution in the age–[Fe/H] plane for the whole subgiant star sample, colour-coded by the stellar number density, N. b, Stellar density distribution in the plane of the azimuthal action Jϕ (equivalent to angular momentum LZ) versus radial action JR. The vertical line delineates Jϕ = 1,500 kpc km s–1, which separates the sample into high angular momentum (yellow background) and low angular momentum regimes. c, Stellar density distribution in the [Fe/H]–[α/Fe] plane. The red solid line separates the sample into high-α and low-α (yellow background) regimes. d, Probability distribution of stellar age p(τ | [Fe/H]), normalized to the peak value for each [Fe/H], for stars with high angular momentum and low [α/Fe] (yellow background regimes in b and c). e, Similar to d but for stars with low angular momentum or high [α/Fe]. The two regimes exhibit a sharp distinction at τ 8 Gyr. Prominent structures are shown for both regimes, such as the V-shaped structure in the late phase (d), and the metal-poor ([Fe/H]  −1) ‘halo’ and metal-rich ([Fe/H]  −1) ‘disk’ sequences in the early phase (e). In the early phase, the two sequences merge at [Fe/H]  −1, but the metal-rich sequence is older than the metal-poor sequence by around 2 Gyr at this metallicity, leading to a Z-shaped structure in p(τ | [Fe/H]).

It turns out that the complexity of p(τ, [Fe/H]) (Fig. 2a) can be unravelled by dividing the sample into two subsamples using stellar quantities that are neither τ nor [Fe/H]: the angular momentum Jϕ (also denoted as LZ) and the ‘α-enhancement’, [α/Fe]. Extensive observations indicate that the majority of stars in the Milky Way formed from gradually enriched gas on high-angular momentum orbits, or the extended (‘thin’) disk4,22, at high Jϕ and low [α/Fe]. It is also well established that the distribution of Galactic stars in the [α/Fe]–[Fe/H] plane is bimodal, with a high-α sequence reflecting rapid enrichment and a low-α sequence reflecting gradual enrichment, which indicates a natural way to divide any sample in the [α/Fe]–[Fe/H] plane8. This inspired our approach to divide our sample into two, separating the dominant sample portion of gradually enriched disk stars with high angular momentum from the rest. Specifically, we used the cut

$${begin{array}{c}begin{array}{cc}{J}_{varphi } > 1500,{rm{kpc}}.{rm{km}}/{rm{s}} & {rm{and}}end{array}\ {begin{array}{cc}[alpha /Fe] > 0.16, & {rm{if}},[{rm{Fe}}/{rm{H}}], > -0.5,\ ,[alpha /Fe] < -0.16[{rm{Fe}}/{rm{H}}],+0.08, & {rm{if}},[{rm{Fe}}/{rm{H}}], > -0.5,end{array}end{array}$$


which is illustrated as a yellow shaded area in Fig. 2b, c. The resulting subsamples in the τ–[Fe/H] plane are shown in Fig. 2d, e, where it is crucial to recall that the sample split involved neither of the quantities on the two axes, τ and [Fe/H]. As we want to focus first on the Milky Way’s elemental enrichment history, rather than its star-formation history, we normalize the distribution p(τ, [Fe/H]) at each [Fe/H] to yield p(τ | [Fe/H]), the age distribution at a given [Fe/H].

Figure 2d, e shows that this cut in angular momentum and [α/Fe] separates the Milky Way’s enrichment history neatly into two distinct age regimes, with a rather sharp transition at τ 8 Gyr. We will therefore refer to these two portions, not clearly apparent in earlier data, as (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{late}}}) and (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{early}}}). The distribution of (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{late}}}) clearly exhibits a V-shape23. This shape is presumably a consequence of the secular evolution of the dynamically quiescent disk; the metal-rich ([Fe/H]  −0.1) branch arises from stars that have migrated from the inner disk to near the Solar radius. The slope of that branch in (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{late}}}) then results from the (negative) radial metallicity gradient in the disk1 and the fact that the stars that have migrated more needed more time to do so, and are hence older. Analogously, we presume the lower branch of (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{late}}}) at [Fe/H]  −0.1 to arise from stars that were born further out and have migrated inwards6. A quantitative comparison with secular evolution models of the Galactic disk4,22 is part of separate ongoing work.

The older stars, reflected in (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{early}}}), show two prominent sequences with distinct [Fe/H](τ) relations. The stars with −2.5 < [Fe/H] < −1.0 reflect the well-established stellar halo population of our Milky Way, whereas the more metal-rich sequence ([Fe/H]  −1) reflects the Milky Way’s inner, high-α (thick) disk24; this designation as an old disk component is also justified by the stars’ angular momentum, as we will show below.

The morphology of the old disk sequence in (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{early}}}) is the most striking feature in Fig. 2e; it reveals an exceptionally clear, continuous and tight age–metallicity relation from [Fe/H]  −1 at 13 Gyr ago all the way to [Fe/H]  0.5 at 7 Gyr ago. A simple model for p(τ | [Fe/H]) of this sequence (Supplementary Information) finds an intrinsic age dispersion of less than 0.82  Gyr at a given [Fe/H] across this 6 Gyr interval (Extended Data Fig. 1). Given the sequence’s slope, this implies that the [Fe/H] dispersion at a given age is smaller than 0.22 dex across the 1.5 dex range in [Fe/H].

Both the halo and old disk sequences extend to [Fe/H]  −1. However, at that [Fe/H] value, the old disk sequence is approximately 2 Gyr older than the halo sequence, leading to a Z-shaped structure in (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{early}}}). This feature is a second aspect of the distribution that has not, to our knowledge, been seen before21.

Formation and enrichment of the Milky Way’s old disk

Tentative hints for some of these features in p(τ | [Fe/H]) have been seen in earlier work24,25 (see the discussion in the Supplementary Information) but these studies lacked the sample size or precision for definitive inferences about the Galactic formation history. Figure 2 shows clearly that the old, high-α ‘thick’ disk of our Milky Way started to form approximately 13 Gyr ago, which is only 0.8 Gyr after the Big Bang19, and extended over 5–6 Gyr, and the interstellar stellar medium (ISM) forming the stars was continually enriched by more than 1 dex, from [Fe/H]  −1 to 0.5. The tightness of this [Fe/H]–age sequence implies that the ISM must have remained spatially mixed thoroughly during this entire period. Had there been any radial (or azimuthal) [Fe/H] variations (or gradients) in excess of 0.2 dex in the star-forming ISM at any time, this would have increased the resulting [Fe/H]–age scatter beyond what is seen. Such gradients, along with orbital migration, are the main reason that the later Galactic disk shows a considerably higher [Fe/H] dispersion at a given age4,26. The results also show that the formation of the Milky Way’s old, α-enhanced disk overlapped in time with the formation of the halo stars: the earliest disk stars are 1–2 Gyr older than the major halo populations at [Fe/H]  −1 (see the Z-shaped structure).

In Fig. 3 we examine the (p{(tau |[{rm{Fe}}/{rm{H}}])}_{{rm{early}}}) distribution more closely by separating stars with at least modest angular momentum, Jϕ > 500 kpc km s–1, from those stars on nearly radial or even retrograde orbits, Jϕ < 500 kpc km s–1. This further sample differentiation by angular momentum leads again to two nearly disjoint p(τ | [Fe/H]) distributions. The first (Fig. 3, upper panel), with mostly [Fe/H] > −1, is dominated by the tight p(τ | [Fe/H]) sequence that we we have already attributed to the old disk. The second, predominately [Fe/H] < −1.2, reflects the halo.

Fig. 3: Probability of stellar distribution in the Jϕ versus [Fe/H] plane, p(τ, [Fe/H]), for stars formed in the early phase.
figure 3

The stars formed in the early phase are divided into Jϕ > 500 kpc km s–1 (upper) and Jϕ < 500 kpc km s–1 (lower). The stellar distribution probability is normalized to the peak value so that the colour from blue to red represents a value from 0 to unity. Note that this is different from p(τ | [Fe/H]) in Fig. 2, which is normalized for each [Fe/H]. The histograms show the distribution integrated over [Fe/H] (top panel) or age (right panels). In the top panel, the age distribution p(τ) is a measure of the relative star-formation history. The dashed curve in red is the result after correcting for the volume selection effect. The vertical dashed line delineates a constant age of 11.2 Gyr, when the star-formation rate reaches its maximum.

Note that Fig. 3, lower panel shows a distinct set of stars with Jϕ < 500 kpc km s–1, for which the p(τ | [Fe/H]) locus indicates that they are the oldest and most metal-poor part of the old disk sequence (see also Extended Data Fig. 2). These stars indicate that some of the oldest members of the old disk sequence were present during an early merger event, by which they were ‘splashed’ to low-angular-momentum orbits27,28. This ancient merger event is presumably the merger with the Gaia-Enceladus satellite galaxy11 (also known as Gaia Sausage10; hereafter Gaia-Sausage-Enceladus), which has contributed most of the Milky Way’s halo stars7,29. The fact that the splashed old disk stars with very little angular momentum are exclusively seen at τ 11 Gyr constitutes strong evidence that the major merger process between the old disk and the Gaia-Sausage-Enceladus satellite galaxy was largely completed 11 Gyr ago. This epoch is 1 Gyr earlier than previous estimates that were based on the lower age limit of the halo stars, 10 Gyr (refs. 11,21,30).

Figure 3 shows the volume-corrected two-dimensional distribution p(τ, [Fe/H]) (see the Supplementary Information for the correction of the volume selection effect), rather than the p(τ | [Fe/H]) of Fig. 2. Figure 3 reveals a remarkable feature, namely that the star-formation rate of the old disk reached a prominent maximum at  around 11.2 Gyr ago, apparently just when the merger with the Gaia-Sausage-Enceladus satellite galaxy was completed, and then continuously declined with time. The most obvious interpretation of this coincidence is that the perturbation from the Gaia-Sausage-Enceladus satellite galaxy greatly enhanced the star formation of the old disk. Note that this star-formation peak among the old disk stars ~11 Gyr ago is very consistent with earlier indications of such a peak based on abundances only31.

To put our results into the bigger picture of galaxy formation and evolution, the multiple assembly phases are seen to be universal among present-day star-forming galaxies. Using the IllustriesTNG simulation, Wang et al.32 showed that galaxy mergers and interactions have played a crucial role in inducing gas inflow, resulting in multiple star formation episodes, intermitted by quiescent phases. Observationally, the best testbed for this theoretical picture would be here at home within our Galaxy. Our study has demonstrated the power of such tests for galactic assembly and enrichment history in the full cosmic timeline, from the very early epoch (τ 13 Gyr or redshift z > 10) to the current time.

Source link

Invest In Films & Earn Yearly

Invest in films & earn yearly

initial deposit returned + 25% min on top

+ 50% return on royalties for life

This will close in 20 seconds