Gaia-Sausage/Enceladus Accretion Event

• Gaia-Sausage/Enceladus is a significant ancient merger where a massive dwarf galaxy merged with the Milky Way 9–11 Gyr ago, leaving distinct kinematic and chemical signatures.
• Its debris, marked by high radial anisotropy and a narrow metallicity distribution, underpins studies of galactic assembly and the evolution of the Milky Way’s halo and disk.
• Integrated chemo-dynamical models and simulations reveal detailed progenitor properties and star formation histories that explain chemical bimodality and the disk transition process.

The Gaia-Sausage/Enceladus (GSE) accretion event refers to a substantial ancient merger between the Milky Way and a massive dwarf galaxy, occurring approximately 8–11 Gyr ago. This event is responsible for the radially biased stellar halo component discovered through Gaia data and constitutes the last major galactic collision in the Milky Way’s history. The GSE remnant is identified through its distinct kinematics—highly eccentric, low-rotation orbits—and unique chemical abundance patterns, which together offer critical insight into both the merger itself and the broader context of early galactic assembly.

1. Progenitor Properties and Assembly Timeline

Multiple dynamical and chemical approaches converge on the following progenitor characteristics:

• Stellar Mass: Estimates from chemo-dynamical modeling and scaling relations yield $M_\star \sim 5\times10^8$–$2\times10^9\,M_\odot$ (Ernandes et al., 10 May 2025, Feuillet et al., 2021), comparable to the present-day LMC.
• Halo (Dark Matter) Mass: Abundance matching with globular cluster population implies $M_{\rm halo} \sim 10^{10}$–$10^{11}\,M_\odot$ (1904.03185, Lane et al., 2023).
• Orbit: Gaia proper motions reveal highly radial plunges ($e\gtrsim0.7$), no net rotation ($L_z \sim 0$), and velocity anisotropy $\beta \sim 0.8$–$0.9$ at large radii (Lane et al., 4 Sep 2025, Folsom et al., 2024).
• Accretion Epoch: The event occurred at redshift $z\simeq2$ (lookback $t_{\rm merge} \sim 9$–$2\times10^9\,M_\odot$0 Gyr) (Ernandes et al., 10 May 2025, González-Koda et al., 27 Feb 2025, Merrow et al., 2023), and contributed up to 50% of the inner stellar halo’s mass.

Recent star formation history analyses (CMD fitting, globular cluster isochrones) reveal two distinct episodes: isolated formation ($2\times10^9\,M_\odot$113.5–12 Gyr, [M/H] $2\times10^9\,M_\odot$2 –1.6) followed by merger-induced quenching at $2\times10^9\,M_\odot$311–10 Gyr ([M/H] $2\times10^9\,M_\odot$4 –0.8), with a minor younger population ([M/H] $2\times10^9\,M_\odot$5 –0.4, age $2\times10^9\,M_\odot$68.5 Gyr) of uncertain association (González-Koda et al., 27 Feb 2025, Aguado-Agelet et al., 27 Feb 2025).

2. Kinematics, Density Profile, and Structure

The GSE debris is best described as a prolate–triaxial, radially anisotropic halo component:

• Velocity Anisotropy:

$2\times10^9\,M_\odot$7 beyond $2\times10^9\,M_\odot$8 kpc, declining to $2\times10^9\,M_\odot$9 at $M_{\rm halo} \sim 10^{10}$0–$M_{\rm halo} \sim 10^{10}$1 kpc; two-component Osipkov-Merritt DFs yield best matches (scale radii $M_{\rm halo} \sim 10^{10}$2 kpc, $M_{\rm halo} \sim 10^{10}$3 kpc, mixture $M_{\rm halo} \sim 10^{10}$4) (Lane et al., 4 Sep 2025).

• Density Profile:

Triaxial broken power law, with inner slope $M_{\rm halo} \sim 10^{10}$5, outer slope $M_{\rm halo} \sim 10^{10}$6, and break radius $M_{\rm halo} \sim 10^{10}$7 kpc (Lane et al., 2023).

• Mass:

Recent high-purity kinematic samples yield $M_{\rm halo} \sim 10^{10}$8, about 15–25% of halo stellar mass in the range $M_{\rm halo} \sim 10^{10}$9 kpc (Lane et al., 2023, Lane et al., 4 Sep 2025).

• Shape:

Axis ratios $10^{11}\,M_\odot$0, major axis inclined $10^{11}\,M_\odot$116° below the plane, oriented toward Galactic rotation (Lane et al., 2023).

3. Chemical Abundance Signatures

Distinct abundance tracks as a result of the GSE's star formation history:

• [Fe/H] Distribution:

Narrow MDF centered at [Fe/H] $10^{11}\,M_\odot$2 –1.15, sigma $10^{11}\,M_\odot$30.3 dex (Feuillet et al., 2021). The metal-weak tail is steeper than that of the full halo, indicating a lack of extremely metal-poor stars ([M/H] $10^{11}\,M_\odot$4 –2.0) compared to less massive accreted systems (Bonifacio et al., 2021).

• Heavy Elements:

[Eu/Ba] rises monotonically with increasing [Fe/H], revealing a $10^{11}\,M_\odot$52 Gyr period of star formation and a sharp quenching at [Fe/H] $10^{11}\,M_\odot$6 –0.5 that prevented the late s-process resurgence (Ernandes et al., 10 May 2025).

• Alpha Elements:

[Mg/Fe] presents a high-$10^{11}\,M_\odot$7 plateau for [Fe/H] $10^{11}\,M_\odot$8 –1.4, followed by a "knee" (Type Ia onset) leading to a low-$10^{11}\,M_\odot$9 track—fundamental to the disk bimodality (Feuillet et al., 2021, Johnson et al., 9 Oct 2025).

Figurative abundance—age relations from globular cluster studies indicate a [Si/Fe] decline and [Eu/Fe] rise correlating ($e\gtrsim0.7$0[Si/Fe]/$e\gtrsim0.7$1 $e\gtrsim0.7$2 –0.11 dex/Gyr; $e\gtrsim0.7$3[Eu/Fe]/$e\gtrsim0.7$4 $e\gtrsim0.7$5 +0.15 dex/Gyr), with [Eu/Si] increasing at +0.26 dex/Gyr (Aguado-Agelet et al., 27 Feb 2025).

4. Metallicity Gradients and Internal Structure

Simulations and data reveal that the GSE progenitor had well-ordered, negative metallicity gradients prior to accretion; these were substantially blurred by merger:

• Pre-infall Gradients:

Radial: $e\gtrsim0.7$6[Fe/H] $e\gtrsim0.7$7 –0.09 to –0.03 dex/kpc; energy-space: $e\gtrsim0.7$8[Fe/H] $e\gtrsim0.7$9 –1.99 to –0.41 dex/($L_z \sim 0$0km$L_z \sim 0$1s$L_z \sim 0$2) (Carrillo et al., 29 Sep 2025).

• Present-day Debris:

The observed halo gradient, mapping mean orbital radius and energy, is –0.014 dex/kpc and –0.28 dex/($L_z \sim 0$3 km$L_z \sim 0$4 s$L_z \sim 0$5) respectively (Khoperskov et al., 2023); the original progenitor likely had a true gradient $L_z \sim 0$6–0.1 dex/kpc.

• Chemical Tagging:

Stars' current energy or mean orbital radius act as proxies for their birth radius in the progenitor, reflecting the central metal-rich enrichment and the extended metal-poor outskirts (Khoperskov et al., 2023).

5. Star Formation History and Quenching

Detailed reconstruction using CMD fitting, elemental ratios, and globular cluster ages yields:

• Two Principal Epochs:

Initial isolated star formation ($L_z \sim 0$713.5–12 Gyr), followed by a merger-induced burst and rapid quenching at $L_z \sim 0$811–10 Gyr (González-Koda et al., 27 Feb 2025, Aguado-Agelet et al., 27 Feb 2025).

• Duration:

Total star formation over $L_z \sim 0$93–4 Gyr, with each globular cluster formation burst lasting $\beta \sim 0.8$0 Gyr and separated by $\beta \sim 0.8$12 Gyr (Aguado-Agelet et al., 27 Feb 2025).

• SFR Quenching Mechanisms:

Merger-driven dynamical heating, gas consumption, and transition from cold to hot gas accretion suppressed SFR, marked by a dip in the gas disc scale length and an abrupt truncation of chemical enrichment at [Fe/H] $\beta \sim 0.8$2 –0.5 (Funakoshi et al., 30 Jul 2025).

6. Role in Disk Structure and Chemical Bimodality

The GSE event exerted profound influence on the Milky Way’s disk and halo structure:

• Disk Transition:

The merger coincided with the thick-to-thin disk switch at $\beta \sim 0.8$310 Gyr, ending the "Great Galactic Starburst" and triggering gas-disc shrinkage and subsequent inside-out thin disk growth (Funakoshi et al., 30 Jul 2025).

• Chemical Bimodality:

Retrograde, radial merger models robustly generate the MW $\beta \sim 0.8$4–Fe bimodality by temporally suppressing SFR, leaving a gap between high-$\beta \sim 0.8$5 and low-$\beta \sim 0.8$6 tracks that depend sensitively on orbit and feedback prescriptions (Johnson et al., 9 Oct 2025, Amarante et al., 1 Dec 2025). Prograde mergers fail to imprint such bimodality.

• Bar Formation:

Cosmological Auriga simulations show that GSE-like mergers (mass ratio $\beta \sim 0.8$7–0.1, high eccentricity) directly seed bar formation in $\beta \sim 0.8$8 Gyr through tidal torque and induced central starbursts (Merrow et al., 2023).

7. Multiplicity of the GSE Event

Recent cosmological simulations suggest that not all GSE-like halos or debris structures are mono-genic:

• Multiple Progenitors:

IllustrisTNG50 analogues reveal %%%%69$10^{11}\,M_\odot$70%%%%32 GSE-like halos are constructed from two mergers, distinguished by star formation histories and chemical patterns (Folsom et al., 2024).

• Observational Discrimination:

Single-merger analogues infall later ($0.9$16 Gyr ago) with higher [Fe/H] and lower [Mg/Fe], whereas two-merger analogues infall earlier ($0.9$211 Gyr ago) and are more $0.9$3-enhanced. However, $0.9$4 radial anisotropy as observed in the MW is rare—favoring a single, massive progenitor origin (Folsom et al., 2024, Kim et al., 2021).

• Outer Halo Debris:

Some studies argue for additional low-mass, high-inclination accretions (e.g., Sequoia-like events), contributing to outer-halo retrograde and tangential debris (Kim et al., 2021, 1904.03185).

---

Summary Table: GSE Progenitor and Remnant Properties

Property	Value / Range	Reference
Stellar Mass	$0.9$5–$0.9$6	(Ernandes et al., 10 May 2025Feuillet et al., 2021Vincenzo et al., 2019)
Halo Mass	$0.9$7–$0.9$8	(1904.03185Lane et al., 2023)
Orbit	$0.9$9, $z\simeq2$0, $z\simeq2$1	(Lane et al., 4 Sep 2025Folsom et al., 2024)
MDF Peak	[Fe/H] $z\simeq2$2 –1.15 to –1.3, FWHM $z\simeq2$30.6 dex	(Feuillet et al., 2021Bonifacio et al., 2021)
Pre-infall Gradient	$z\simeq2$4 –0.1 dex/kpc (radial)	(Carrillo et al., 29 Sep 2025Khoperskov et al., 2023)
Merger Epoch	$z\simeq2$5 ($z\simeq2$69–11 Gyr ago)	(Ernandes et al., 10 May 2025González-Koda et al., 27 Feb 2025Merrow et al., 2023)
Disk Transition	Thick-to-thin at $z\simeq2$710 Gyr	(Funakoshi et al., 30 Jul 2025Merrow et al., 2023)
Key Event Types	Single/dual major merger, retrograde orbit	(Folsom et al., 2024Johnson et al., 9 Oct 2025)

---

• (Ernandes et al., 10 May 2025, Feuillet et al., 2021, Lane et al., 2023, Lane et al., 4 Sep 2025, Carrillo et al., 29 Sep 2025, Funakoshi et al., 30 Jul 2025, González-Koda et al., 27 Feb 2025, Aguado-Agelet et al., 27 Feb 2025, Johnson et al., 9 Oct 2025, Amarante et al., 1 Dec 2025, Khoperskov et al., 2023, Bonifacio et al., 2021, 1904.03185, Folsom et al., 2024, Merrow et al., 2023, Kim et al., 2021, Vincenzo et al., 2019).

The GSE merger provides a unique laboratory for dissecting the interplay between orbital dynamics, chemical enrichment, star formation shutdown, disk structure evolution, and the origin of key stellar populations in the Milky Way. Its chemical and dynamical signatures remain benchmarks for accretion-driven galaxy assembly scenarios.