Gaia-Enceladus-Sausage (GES) Members

• GES members are the accreted halo stars, defined by high orbital eccentricities, specific [Fe/H] ranges, and low [Al/Fe] signatures via dynamical and chemical selection.
• They exhibit highly radial orbits with pericenters around 3–5 kpc and apocenters of 10–30 kpc, consistent with a major dwarf-galaxy merger 8–12 Gyr ago.
• Age and abundance trends in both field stars and associated globular clusters trace bursty star formation episodes, refining our reconstruction of the Milky Way’s merger history.

Gaia-Enceladus-Sausage (GES) members constitute the dominant, highly radial, metal-intermediate accreted debris population in the Galactic halo, originating from one or more major dwarf-galaxy mergers 8–12 Gyr ago. Robust identification and characterization of GES members have been central to reconstructing the Milky Way’s merger history, with recent studies revealing subtle chemical, dynamical, and chronological features, as well as a population of associated globular clusters.

1. Dynamical and Chemodynamical Selection of GES Members

GES member identification relies on integrals-of-motion (IoM), action/energy space, and chemical tagging. Canonical dynamical identifiers include:

• Integrals of motion: Orbital energy $$E = \frac{1}{2}(v_R^2 + v_\phi^2 + v_z^2) + \Phi(R)$$, vertical angular momentum $L_z = R\,v_\phi$, and $L_\perp = \sqrt{L_x^2 + L_y^2}$.
• Action space cuts: $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$ and $|L_z|\leq 500$ kpc km s$^{-1}$ select the high-eccentricity, low-rotation population (Feuillet et al., 2020, González-Koda et al., 27 Feb 2025).
• Orbital eccentricity: $e > 0.7$, typically $e>0.8$ for the "sausage" component (Perottoni et al., 2022, Liu et al., 2024).
• Energy and angular momentum box: $-2.0 \times 10^5 < E < -1.2 \times 10^5$ km$^2$/s$L_z = R\,v_\phi$0, $L_z = R\,v_\phi$1–1500 kpc km s$L_z = R\,v_\phi$2 for both field stars and globular clusters (Aguado-Agelet et al., 27 Feb 2025, Perottoni et al., 2022).

Chemodynamical selection is further refined by:

• Metallicity: $L_z = R\,v_\phi$3; median values span $L_z = R\,v_\phi$4 to $L_z = R\,v_\phi$5 dex, with observed dispersions $L_z = R\,v_\phi$6–$L_z = R\,v_\phi$7 dex (Myeong et al., 2022, Feuillet et al., 2021).
• Alpha and odd-Z elements: GSE has uniformly low $L_z = R\,v_\phi$8 (typically $L_z = R\,v_\phi$9 dex), a distinct $L_\perp = \sqrt{L_x^2 + L_y^2}$0–$L_\perp = \sqrt{L_x^2 + L_y^2}$1 track with a "knee" at $L_\perp = \sqrt{L_x^2 + L_y^2}$2, and $L_\perp = \sqrt{L_x^2 + L_y^2}$3 (Feuillet et al., 2021, Myeong et al., 2022).
• Probabilistic mixture modeling: Gaussian Mixture Models and unsupervised algorithms applied to abundance and orbital parameters assign GSE memberships with high purity (Myeong et al., 2022, Perottoni et al., 2022).

Recent studies employing HDBSCAN clustering in $L_\perp = \sqrt{L_x^2 + L_y^2}$4 space have yielded member lists of over 1,300 stars from Gaia DR3, with eccentricity and angular momentum cuts consistent with previous methodology (Liu et al., 2024).

2. Kinematic, Spatial, and Orbital Properties

GES members are defined by highly radial orbits and low net angular momentum:

• Eccentricity: Median $L_\perp = \sqrt{L_x^2 + L_y^2}$5, majority of orbits satisfy $L_\perp = \sqrt{L_x^2 + L_y^2}$6; for RR Lyrae, GSE-fraction stars have $L_\perp = \sqrt{L_x^2 + L_y^2}$7 (Kunder et al., 15 Jul 2025, Liu et al., 2024).
• Pericenter and apocenter: $L_\perp = \sqrt{L_x^2 + L_y^2}$8–5 kpc; $L_\perp = \sqrt{L_x^2 + L_y^2}$9–30 kpc (Perottoni et al., 2022, González-Koda et al., 27 Feb 2025).
• Spatial extent: Members are distributed across the halo with a density profile $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$0 and flattening axis ratio $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$1 increasing from $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$2 at $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$3 kpc to $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$4 at $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$5 kpc (Wu et al., 2022).
• Phase-space signature: Vertical bar in $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$6–$$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$7, "sausage"-shaped locus in $$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$8–$$\sqrt{J_R} \in [30,50]\ \mathrm{(kpc\,km\,s^{-1})^{1/2}}$$9, and spread in Galactocentric radii to at least 20 kpc (Ding et al., 9 Sep 2025).
• Bulk rotation: $|L_z|\leq 500$0 near zero or slight retrograde; $|L_z|\leq 500$1 kpc km/s.

3. Chemical Abundance Patterns and Metallicity Distribution

GES members show distinctive chemical evolution:

• Metallicity distribution: Median $|L_z|\leq 500$2 in APOGEE+Gaia, $|L_z|\leq 500$3 from APOGEE-based samples, $|L_z|\leq 500$4 from SkyMapper–Gaia; MDF dispersions $|L_z|\leq 500$5–$|L_z|\leq 500$6 dex (Myeong et al., 2022, Feuillet et al., 2021).
• Alpha–elements: $|L_z|\leq 500$7, declining from a plateau at low $|L_z|\leq 500$8 to lower values at higher metallicity ("knee" at $|L_z|\leq 500$9) (Myeong et al., 2022).
• Odd-Z and Fe-peak elements: $^{-1}$0 is systematically low (–0.24 to –0.25 dex), $^{-1}$1 enhanced ($^{-1}$2), and $^{-1}$3 subsolar (–0.3 dex) (Myeong et al., 2022, Feuillet et al., 2021).
• r-process enrichment: GSE is unusually rich in r-process elements; $^{-1}$4 and $^{-1}$5 (Myeong et al., 2022, Aguado-Agelet et al., 27 Feb 2025).

4. Age Distributions, Star Formation History, and the Age–Metallicity Relation

GES member ages and SFH have been characterized using isochrone fitting and CMD-based methods:

• Dominant age: Mean/median stellar age $^{-1}$6–$^{-1}$7 Gyr, with the main population spanning $^{-1}$8–$^{-1}$9 Gyr and full width at half maximum $e > 0.7$0 Gyr (Feuillet et al., 2021, González-Koda et al., 27 Feb 2025).
• SFH and AMR: Multiple studies identify two major episodes of star formation:
  • "Isolated" evolution phase at $e > 0.7$1 Gyr, with $e > 0.7$2.
  • "Merger-induced" phase at $e > 0.7$3–$e > 0.7$4 Gyr, with $e > 0.7$5 increasing to $e > 0.7$6 (González-Koda et al., 27 Feb 2025).
• Recent population: A fourth, younger subpopulation at $e > 0.7$7 Gyr and $e > 0.7$8 is observed, but its association with GSE is unresolved (González-Koda et al., 27 Feb 2025).
• AMR functional form: For globular clusters,

$e > 0.7$9

with 1$e>0.8$0 scatter of 0.15 Gyr (Aguado-Agelet et al., 27 Feb 2025).

• Cluster and field AMR: Field star and globular cluster AMRs coincide, confirming that cluster formation tracks global GSE chemical evolution (González-Koda et al., 27 Feb 2025, Aguado-Agelet et al., 27 Feb 2025).

5. GSE-Associated Globular Clusters

Selection of globular clusters (GCs) as GSE members is based on integrals-of-motion, specifically:

• Selection box: $e>0.8$1 km$e>0.8$2/s$e>0.8$3, $e>0.8$4 kpc km s$e>0.8$5 (Aguado-Agelet et al., 27 Feb 2025).
• Outliers: Four clusters (NGC 288, NGC 6205, NGC 5286, NGC 7099) are dynamical outliers or show chemical/age deviation; e.g., NGC 288 and NGC 6205 are likely in-situ "Splash" clusters, NGC 7099 is younger and retrograde (Sequoia member), NGC 5286 is older with high Eu/Si (Aguado-Agelet et al., 27 Feb 2025).
• Two formation epochs: Nine bona-fide GSE clusters exhibit bimodal cluster ages:
  • First cohort formed at $e>0.8$6 Gyr (NGC 2298, 5897, 6341, 6779, 7089),
  • Second at $e>0.8$7 Gyr (NGC 362, 1261, 1851, 2808).
  • Duration of both episodes is short ($e>0.8$80.3 Gyr each), separation $e>0.8$92 Gyr, with the timing attributed to pericentric passage and final coalescence (Aguado-Agelet et al., 27 Feb 2025).
• Age–Si–Eu relation: $-2.0 \times 10^5 < E < -1.2 \times 10^5$0 rises with decreasing age, while $-2.0 \times 10^5 < E < -1.2 \times 10^5$1 declines and $-2.0 \times 10^5 < E < -1.2 \times 10^5$2 decreases, tracing the onset of r-process production and SNIa enrichment; these trends provide additional membership discrimination (Aguado-Agelet et al., 27 Feb 2025).

6. Multiple Progenitor Events and Substructure

Several analyses support a composite origin for the classical GSE structure:

• Multiple accretion events: Action-inclination decomposition reveals at least three distinct major merger components—low-inclination prograde and retrograde "Sausage" events, and a high-inclination, nearly radial "Enceladus" component—plus additional minor outer-halo retrograde dwarfs (Kim et al., 2021).
• Distinct metallicity peaks: Kinematic/chemical separation finds $-2.0 \times 10^5 < E < -1.2 \times 10^5$3 peaks at $-2.0 \times 10^5 < E < -1.2 \times 10^5$4 (LOI, radial), $-2.0 \times 10^5 < E < -1.2 \times 10^5$5 (LOI, tangential), $-2.0 \times 10^5 < E < -1.2 \times 10^5$6 (HOI, radial), $-2.0 \times 10^5 < E < -1.2 \times 10^5$7 (HOI, tangential) that map onto different mass progenitors (Kim et al., 2021).
• Simulation results: TNG50 and VINTERGATAN-GM cosmological analogues show that both single- and double-merger GSE analogues yield similar present-day chemical–dynamical signatures; explicit mass ratios and star-formation burst features can only be teased apart with additional chronochemical constraints (Folsom et al., 2024, Rey et al., 2022).

7. Summary Table: Core Properties of GES Members

Property (Field Stars)	Value or Range	Source(s)
Median $-2.0 \times 10^5 < E < -1.2 \times 10^5$8	–1.20 to –1.15 (APOGEE), –1.17 (SM)	(Myeong et al., 2022, Feuillet et al., 2021, Feuillet et al., 2020)
Dispersion $-2.0 \times 10^5 < E < -1.2 \times 10^5$9	0.23–0.34 dex	(Myeong et al., 2022, Feuillet et al., 2021)
Median $^2$0	0.18 ($^2$1)	(Myeong et al., 2022)
Median $^2$2	–0.24 to –0.25 ($^2$3)	(Myeong et al., 2022, Feuillet et al., 2021)
Median Age	10–12 Gyr	(Feuillet et al., 2021, González-Koda et al., 27 Feb 2025)
Eccentricity $^2$4	$^2$50.8 (majority), selection $^2$6	(Perottoni et al., 2022, Liu et al., 2024)
Pericenter $^2$7	$^2$83–5 kpc	(Perottoni et al., 2022, Perottoni et al., 2022)
Apocenter $^2$9	$L_z = R\,v_\phi$0010–30 kpc	(Perottoni et al., 2022, Perottoni et al., 2022)
$L_z = R\,v_\phi$01 Range	$L_z = R\,v_\phi$02 or $L_z = R\,v_\phi$03 kpc km/s	(Feuillet et al., 2020, Perottoni et al., 2022)

*SM: SkyMapper *Full membership selection and additional covariance information are provided in the referenced works.

8. Implications and Prospects

GES membership definitions have matured from broad orbital cuts to multi-dimensional probabilistic models incorporating kinematics, actions, and precise chemical tagging. The population encompasses a tightly correlated system of field stars and globular clusters, tracing a bursty, merger-driven star-formation history that dominates the local stellar halo. Detailed study of age–abundance trends (e.g., [Eu/Fe], [Eu/Si]), chemodynamical substructure, and the full distribution of metallicity and age enables separation of genuine GSE members from in-situ contaminants and allows for nuanced distinctions between single versus multi-progenitor origin scenarios. Simulation analogues confirm that similar present-day properties can originate from a range of early assembly pathways, underscoring the need for high-precision ages and multi-element abundance diagnostics to fully reconstruct the Milky Way’s merger sequence (Rey et al., 2022, Folsom et al., 2024, Aguado-Agelet et al., 27 Feb 2025, González-Koda et al., 27 Feb 2025).