---
title: Two-Infall Scenario in the Milky Way Disc
url: https://www.emergentmind.com/topics/two-infall-scenario
type: topic
---

# Two-Infall Scenario in the Milky Way Disc

The Two-Infall Scenario is a Galactic chemical-evolution framework in which the Milky Way disc is assembled through two main gas-accretion episodes rather than through one smooth, continuous supply. In its modern delayed formulation, the first infall is rapid and produces the old high-\(\alpha\) population, identified with the thick disc, whereas the second infall is slower and delayed by roughly \(4\)–\(4.6\) Gyr, producing the younger low-\(\alpha\) population, identified with the thin disc. The delayed second episode dilutes a previously enriched interstellar medium (ISM), generates a characteristic loop in abundance space, and thereby explains the observed high-\(\alpha\)/low-\(\alpha\) bimodality in the solar neighbourhood [1809.00914] [2204.07597].

## 1. Conceptual structure and historical reformulation

In the modern literature, the scenario is not merely the statement that there were “two phases” of disc growth. Its defining claim is that the two chemically distinct disc sequences are the imprint of two distinct gas-accretion episodes separated by a substantial lull in accretion and star formation. The first episode forms the high-\(\alpha\), thick-disc-like population in a rapid, intense star-forming phase dominated initially by core-collapse supernova enrichment. The second episode begins much later, brings in low-metallicity gas, lowers \([\mathrm{Fe/H}]\) at nearly fixed \([\alpha/\mathrm{Fe}]\), and then forms the low-\(\alpha\), thin-disc-like sequence as Type Ia supernova iron accumulates [1809.00914].

A major revision relative to the classical Chiappini-style implementation is the timing of the second infall. The classical version typically placed the second onset after \(\sim 1\) Gyr, whereas the revised model constrained by asteroseismic ages requires a much longer delay, \(t_{\max}=4.3\) Gyr, and identifies that delay as the crucial assumption for reproducing both abundances and ages [1809.00914]. In this reformulation, the scenario is better described as a delayed-accretion model than as a threshold-driven star-formation-gap model: the key physical ingredient is renewed gas supply after a prolonged reduction in accretion, not necessarily a strictly imposed cessation of star formation.

This formulation was tested directly against the “parallel” alternative, in which the high-\(\alpha\) and low-\(\alpha\) sequences form coevally in independent gas reservoirs. The parallel model can reproduce the bimodality in chemical space qualitatively, but it fails once stellar ages are imposed, especially because it predicts a low-\(\alpha\) population that is too old and lacks the delayed-infall dilution signature [1809.00914].

## 2. Mathematical formulation in Galactic chemical evolution

The standard implementation is a one-zone chemical-evolution model for the solar annulus or solar vicinity, usually centered at \(R=8\) kpc. A representative governing equation for the surface density of species \(i\) is
$$
\dot \Sigma_i(R,t)= -X_i(R,t)\,\mathrm{SFR}(R,t) + \mathcal{R}_i(R,t) + \mathcal{B}_i(R,t),
$$
where \(X_i\) is the ISM mass fraction, \(\mathcal{R}_i\) the stellar return term, and \(\mathcal{B}_i\) the gas-accretion term [1809.00914].

In the delayed two-infall formulation used for the solar vicinity, the infall law is written as
$$
\mathcal{I}_i(t,z=0) \equiv (\mathcal{X}_i)_{\rm inf} \left[ \mathcal{A}_1 \, e^{-t/t_1} + \theta(t-t_{\rm max}) \mathcal{A}_2 \, e^{-(t-t_{\rm max})/t_2} \right],
$$
with \(t_1\) and \(t_2\) the infall timescales of the first and second episodes, \(t_{\rm max}\) the onset time of the second accretion episode, and \(\theta\) the Heaviside step function [2204.07597]. The normalizations are fixed so as to reproduce the present-day total surface mass density of \(47.1 \pm 3.4~M_\odot\,\mathrm{pc}^{-2}\) in the solar neighbourhood [2204.07597].

Star formation follows a Kennicutt–Schmidt law,
$$
\psi \propto \nu \sigma_g^k, \qquad k=1.5,
$$
with \(\nu\) the star-formation efficiency and \(\sigma_g\) the gas surface density. In the APOGEE-constrained ES21 realization, the adopted efficiencies are \(\nu_1=2~\mathrm{Gyr}^{-1}\) and \(\nu_2=1~\mathrm{Gyr}^{-1}\) for the high-\(\alpha\) and low-\(\alpha\) phases, respectively [2204.07597]. In the APOKASC age-constrained revised two-infall model, the preferred parameters are \(\tau_{D1}=0.1~\mathrm{Gyr}\), \(\tau_{D2}=8~\mathrm{Gyr}\), \(t_{\max}=4.3~\mathrm{Gyr}\), \(\nu=1.3~\mathrm{Gyr}^{-1}\), \(\Sigma_{\mathrm{tot1}}=8~M_\odot\,\mathrm{pc}^{-2}\), and \(\Sigma_{\mathrm{tot2}}=64~M_\odot\,\mathrm{pc}^{-2}\) [1809.00914].

These models adopt the Scalo (1986) IMF, constant in time and space, and standard delayed enrichment from Type Ia supernovae in the Matteucci formalism [1809.00914] [2204.07597]. In the local disc realizations discussed here, galactic winds are not included, and radial migration is not included in the model proper [1809.00914] [2204.07597].

The abundance-space morphology follows directly from this setup. After the second infall begins, low-metallicity gas dilutes the ISM and produces a nearly horizontal excursion in \([\alpha/\mathrm{Fe}]\)–\([\mathrm{Fe}/\mathrm{H}]\), because dilution lowers metallicity more directly than the abundance ratio. Renewed star formation then raises \([\alpha/\mathrm{Fe}]\) through core-collapse enrichment, and later Type Ia supernovae lower \([\alpha/\mathrm{Fe}]\) again, creating the loop that populates the low-\(\alpha\) branch [1809.00914]. For comparison with APOKASC, one implementation uses
$$
[\mathrm{M/H}] = [\mathrm{Fe/H}] + \log\left(0.638\times 10^{[\alpha/\mathrm{Fe}]} + 0.362\right),
$$
following Salaris et al. (1993) [1809.00914].

## 3. Empirical support from ages, abundance planes, and metallicity distributions

The decisive observational shift in favour of the revised delayed scenario came from asteroseismic ages. In the APOKASC sample analyzed in this context, the high-\(\alpha\) population peaks at about \(11\) Gyr, while the low-\(\alpha\) population peaks at about \(2\) Gyr [1809.00914]. After kinematic pruning and selection-function corrections, the final APOKASC sample used for the age test contains \(1197\) red giants [1809.00914].

The revised two-infall model reproduces several observables simultaneously once age and metallicity uncertainties are convolved into the synthetic predictions. These include the \([\alpha/\mathrm{Fe}]\)–\([\mathrm{Fe}/\mathrm{H}]\) relation at different Galactic epochs, the age–metallicity relation, the time evolution of \([\alpha/\mathrm{Fe}]\), the age distributions of the high-\(\alpha\) and low-\(\alpha\) populations, and the metallicity distribution function [1809.00914]. Quantitatively, the observed median age of the high-\(\alpha\) population is \(10.40^{+1.86}_{-2.71}\) Gyr, compared with \(10.53^{+2.23}_{-2.14}\) Gyr in the synthetic revised two-infall model; for the low-\(\alpha\) population, the observed median is \(3.76^{+3.32}_{-2.17}\) Gyr, compared with \(4.31^{+3.16}_{-2.79}\) Gyr in the model [1809.00914].

The same age information disfavors the parallel scenario. Its synthetic low-\(\alpha\) age distribution has median \(7.06^{+4.13}_{-3.76}\) Gyr, substantially older than the observed \(3.76^{+3.32}_{-2.17}\) Gyr, and it lacks the horizontal dilution signature that the delayed second infall produces in the abundance plane [1809.00914].

A recurrent misconception is that the delayed two-infall picture requires an explicitly imposed star-formation threshold or a strict hiatus. The revised local model tested against APOKASC does not include a star-formation threshold; rather, it uses a delayed second gas-accretion event to mimic a period of substantially lowered gas accretion, which is sufficient to generate the observed age-separated chemical structure [1809.00914].

## 4. Vertical chemo-dynamical signatures

A further development tests whether the same delayed-infall mechanism that explains the abundance bimodality also predicts the vertical distribution of \([\mathrm{Mg/Fe}]\). In this extension, the chemistry is first computed in the solar vicinity and then linked to vertical structure by assigning each simple stellar population (SSP) an age, chemical composition, mass, vertical action, and maximum height \(|z_{\max}|\) above the plane [2204.07597].

The SSPs are defined as “an assembly of coeval and chemically homogeneous stars,” formed with time-step \(\Delta \tau = 0.005~\mathrm{Gyr}\), yielding \(2740\) SSPs in total [2204.07597]. The age–vertical action relation is taken from Ting & Rix (2019),
$$
\widehat{J_z}(R_\odot,\tau) = 0.91 + 1.81 \cdot \tau^{1.09} \quad [\mathrm{kpc\ km\ s}^{-1}],
$$
and each SSP orbit is integrated in the static Galactic potential `MWPotential2014` from `galpy`, consisting of bulge, Miyamoto–Nagai disc, and NFW halo [2204.07597]. Scatter in the empirical action relation is modeled with Gaussian perturbations, using \(\sigma_{J_z}=0.5\,J_z\) and \(\sigma_{J_z}=1.0\,J_z\) as test cases [2204.07597].

This construction yields a nontrivial result: the delayed-infall model reproduces not only the high-\(\alpha\)/low-\(\alpha\) bimodality in the \([\mathrm{Mg/Fe}]\)–\([\mathrm{Fe/H}]\) plane, but also a corresponding vertical dichotomy in \(|z_{\max}|\)–\([\mathrm{Mg/Fe}]\) [2204.07597]. For stars younger than \(8\) Gyr, which are mostly low-\(\alpha\) thin-disc stars and for which the Ting relation is directly applicable, the predicted \(|z_{\max}|\)–\([\mathrm{Mg/Fe}]\) relation agrees with the APOGEE DR16 plus astroNN data, especially for the \(|z|<2\) kpc and \(|z|<1\) kpc samples. When high-\(\alpha\) stars are included, the model reproduces the existence of the vertical dichotomy for the \(|z|<2\) kpc sample [2204.07597].

The main tension concerns the high-\(\alpha\) regime. The model predicts too flat a growth of \(|z_{\max}|\) with \([\mathrm{Mg/Fe}]\) for high-\(\alpha\) stars compared with the observed median trend. Proposed explanations in the study include halo contamination in the \(|z|<2\) kpc sample, external perturbations such as minor mergers, dynamical heating associated with events such as Gaia-Enceladus, the inadequacy of extrapolating the Ting relation to old high-\(\alpha\) stars, the limitations of a static-potential approximation, and sizable observational uncertainties in vertical action estimates [2204.07597]. If the data are dissected chemically into high-\(\alpha\) and low-\(\alpha\) populations, the observed \(|z_{\max}|\) distributions agree much better with the model, especially for \(\sigma_{J_z}=1\cdot J_z\) [2204.07597].

## 5. Gaia DR3 and the three-infall refinement

Gaia DR3 introduced a sharper test of the scenario by exposing a population of young, often massive, low-\(\alpha\) disc stars with unexpectedly low metallicity and low \([\mathrm{X/Fe}]\) values for several \(\alpha\)-elements. In the solar vicinity, the comparison sample is restricted to giant stars with \(\log g<3.5\) and guiding radii \(R_g\in[8.1,8.4]\) kpc, with ages from the Gaia DR3 catalogue of Kordopatis et al. (2022) and relative age errors \(<0.5\) [2206.12436]. The abundance planes examined are \([\mathrm{X/Fe}]\) versus \([\mathrm{M/H}]\) for \(X=\mathrm{Mg}, \mathrm{Si}, \mathrm{Ca}, \mathrm{Ti}, \alpha\) [2206.12436].

The modern delayed two-infall model remains broadly successful for the classical high-\(\alpha\)/low-\(\alpha\) dichotomy, but it fails specifically for the youngest age bins. Gaia DR3 shows many young stars with under-solar metallicity and low \([\mathrm{Mg/Fe}]\), \([\mathrm{Ca/Fe}]\), \([\mathrm{Ti/Fe}]\), and similar low-\(\alpha\) behaviour; Recio-Blanco et al. (2022b), as summarized in this context, identified many massive stars with under-solar \([\mathrm{M/H}]\), especially Ca-poor objects with \([\mathrm{Ca/Fe}]\) as low as \(\sim -0.3\) dex over \([\mathrm{M/H}] \in [-0.5,0]\) [2206.12436]. The older two-infall model, when projected into young age bins, forms too many young stars at too high metallicity.

The proposed response is not a rejection of the Two-Infall Scenario but an explicit extension of it. The three-infall model keeps the first two infalls essentially identical to the delayed two-infall baseline and splits the low-\(\alpha\) disc formation into two distinct gas-accretion episodes [2206.12436]. The infall law becomes a sum of three exponentials with Heaviside onsets, corresponding to high-\(\alpha\), low-\(\alpha\) part I, and low-\(\alpha\) part II [2206.12436]. The first infall remains very rapid, with \(\tau_1=0.103~\mathrm{Gyr}\); the second has \(\tau_2=4.110~\mathrm{Gyr}\) and onset \(t_{\rm max1}=4.085~\mathrm{Gyr}\); the third begins at \(t_{\rm max2}=11.000~\mathrm{Gyr}\), i.e. about \(2.7\) Gyr ago for \(t_G=13.7\) Gyr, and has \(\tau_3=0.150~\mathrm{Gyr}\) [2206.12436]. The star-formation efficiencies are \(\nu_1=2.0~\mathrm{Gyr}^{-1}\), \(\nu_2=1.3~\mathrm{Gyr}^{-1}\), and \(\nu_3=0.5~\mathrm{Gyr}^{-1}\), with \(\sigma_2/\sigma_3=2.33\), so the low-\(\alpha\) mass split is about \(70\%\) in low-\(\alpha\) part I and \(30\%\) in low-\(\alpha\) part II [2206.12436].

Physically, the third infall is interpreted as a recent dilution event: fresh primordial or mildly enriched gas lowers the ISM metallicity and forms chemically impoverished young stars. The model also uses mild pre-enrichment for the two low-\(\alpha\) infalls, setting their gas composition to the high-\(\alpha\)-phase model abundance at \([\mathrm{Fe/H}] = -0.75\) dex divided by a factor of five [2206.12436]. In abundance space, the first two infalls preserve the standard explanation for the main bifurcation, while the third introduces a second, smaller loop inside the low-\(\alpha\) regime. The central conclusion is therefore precise: Gaia DR3 pushes the field beyond the two-infall model by extending its logic, not by replacing its foundation [2206.12436].

## 6. Current challenges from large stellar age catalogues

Large age catalogues across the disc have tightened the allowed parameter space of the scenario considerably. In a multi-zone VICE implementation spanning \(0 \leq R_{\rm gal} < 20\) kpc in annuli of width \(100\) pc, with stellar migration, APOGEE DR17 abundances, and age catalogues from Leung et al. (2023) and Roberts et al. (2025, submitted), the Two-Infall Scenario still reproduces the existence of the \([\alpha/{\rm Fe}]\) bimodality, but it struggles with several features of the age–abundance structure [2508.00988].

The first challenge is the age–metallicity relation of the thin disc. Because the second infall introduces a large quantity of low-metallicity gas, the model generically predicts a substantial dilution event followed by several Gyr of re-enrichment. The data, however, show that the mode of the stellar metallicity distribution is remarkably constant over much of the past \(6\)–\(8\) Gyr, and more broadly that the thin-disc epoch is much flatter in age–metallicity than a strong dilution-reset event would predict [2508.00988]. In the solar annulus, the fiducial models underpredict the metallicity of \(5\)–\(8\) Gyr old stars by up to \(\sim 0.3\)–\(0.5\) dex [2508.00988].

The second challenge is the age distribution of local metal-rich stars. In the two-infall picture, the most metal-rich stars at fixed radius should arise either from the metal-rich tip of the first sequence or from the very late endpoint of the second, so the model predicts a bimodal age distribution. By contrast, in APOGEE with neural-network ages, local metal-rich stars with \([{\rm Fe/H}] \ge +0.1\) have a broad, mostly unimodal age distribution peaking around \(\sim 5\) Gyr, with very few older than \(\sim 10\) Gyr [2508.00988].

A third generic issue is an extra intermediate-\(\alpha\) overdensity near the turnover after the second infall, around \([{\rm Fe/H}] \approx -0.3\) and \([{\rm O/Fe}] \approx 0.1\)–0.2 depending on parameters. The models predict this feature, whereas APOGEE shows a trough there [2508.00988].

Several modifications mitigate but do not remove these tensions. Longer \(\tau_2\) smooths the low-\(\alpha\) loop, larger thick/thin mass ratios reduce dilution, stronger migration broadens local distributions, and pre-enrichment of the second infall to \({\rm [X/H]}_{\rm CGM}=-0.5\) or \(-0.7\) weakens dilution to about the \(\sim0.3\)-dex level [2508.00988]. Yet the age structure of local metal-rich stars remains fundamentally incorrect in the tested models. The resulting assessment is not that the scenario is ruled out, but that it remains viable only in modified and tightly restricted form [2508.00988].

## 7. Extensions, reinterpretations, and terminological scope

The two-infall logic has also been extended beyond the local Galactic disc. In a bulge-specific one-zone OMEGA++ implementation constrained by more than \(30{,}000\) Galactic chemical-evolution models and machine learning, the preferred solution for the Milky Way bulge consists of an early rapid collapse with \(t_1\sim 0.1\) Gyr, \(\tau_1\sim 0.09\) Gyr, and \({\rm SFE}_1\sim 3~\mathrm{Gyr}^{-1}\), followed by a delayed second episode at \(t_2\sim 5.1\) Gyr with \(\tau_2\sim 1.7\) Gyr, \(\sigma_2\sim 0.69\), and reduced efficiency \(\Delta{\rm SFE}\sim 0.72\) [2512.08090]. In that reinterpretation, the first phase builds the majority of the bulge and the second is chemically required to reproduce the metal-rich peak of the bulge metallicity distribution function and the decline in \([\alpha/{\rm Fe}]\) at high \([\mathrm{Fe/H}]\). The study explicitly emphasizes that the posterior supports more than one gas-supply episode, but not necessarily exactly two sharply separated physical events; in a one-zone model, multiple minor late inflows would be absorbed into the effective second episode [2512.08090].

The phrase “two infall” also appears in other astrophysical settings, but these uses are conceptually distinct. In circumstellar-disk hydrodynamics, late inclined infall onto an already formed star–disk system can create a second-generation outer disk that is misaligned with respect to the primordial inner disk, or even counter-rotating in retrograde cases; this is a physically concrete two-stage accretion picture, but it is not the Galactic chemical-evolution Two-Infall Scenario [2110.04309]. Conversely, methanol studies of high-mass star-forming regions use a modified two-layer radiative-transfer model to infer inward motions from redshifted absorption; those results support large-scale inward collapse and, in W31C, possibly hierarchical or multi-layer infall, but they do not identify two temporally separate infall episodes and should not be conflated with the Galactic chemical-evolution usage of the term [2201.01792].

Within Galactic archaeology, the Two-Infall Scenario therefore remains best understood as an effective description of Milky Way enrichment history in which the main high-\(\alpha\)/low-\(\alpha\) dichotomy arises from two major gas-accretion episodes, with the modern debate centered less on whether two episodes are useful than on how delayed, how chemically disruptive, and how literally discrete those episodes must be.

Source: https://www.emergentmind.com/topics/two-infall-scenario