---
title: 'Thamnos: Stellar-Halo Substructure'
url: https://www.emergentmind.com/topics/thamnos
type: topic
---

# Thamnos: Stellar-Halo Substructure

Thamnos is a stellar-halo substructure of the Milky Way, identified in the Gaia era as a low-energy retrograde component of the local halo and generally interpreted as debris from an ancient accretion event [1909.08924]. Across the subsequent literature, Thamnos has been recovered in multiple dynamical and chemo-dynamical searches, but its astrophysical status remains contested: some studies treat it as a distinct low-mass progenitor, while others argue that the usual “Thamnos” selection region in integral-of-motion space is a superposition of genuine accreted debris, Gaia-Sausage-Enceladus (GSE), heated in situ halo or disc stars, and possibly $\omega$ Cen-related material [2509.13408].

## 1. Discovery and basic dynamical identity

Thamnos was first isolated in a local halo sample built from Gaia DR2, supplemented with RAVE, APOGEE, and LAMOST, by applying HDBSCAN in the four-dimensional space of specific orbital energy $E$, angular momentum $L_z$, orbital eccentricity, and metallicity [1909.08924]. In that work, two overlapping clumps were linked and collectively named Thamnos. The defining dynamical quantities were
$$
E = \tfrac12\,(v_r^2+v_\theta^2+v_\phi^2)+\Phi(r),
\qquad
L_z = R\,v_\phi.
$$

In the original characterization, Thamnos occupied a more tightly bound region than Sequoia, with Thamnos 1 at $E/10^5\,(\mathrm{km^2\,s^{-2}})\in[-1.65,-1.45]$ and circularity $< -0.75$, and Thamnos 2 at $E/10^5\,(\mathrm{km^2\,s^{-2}})\in[-1.80,-1.60]$ and $-0.75<\mathrm{circ}< -0.40$ [1909.08924]. The substructure was described as retrograde, low-inclination, and mildly eccentric, with $\langle v_\phi\rangle \simeq -150\,\mathrm{km\,s^{-1}}$ for Thamnos 2 and $\simeq -200\,\mathrm{km\,s^{-1}}$ for Thamnos 1, and with eccentricities $\sim0.3$–$0.6$ [1909.08924].

Later analyses recovered the same component in larger local-halo samples. In a very metal-poor Gaia DR3 + LAMOST DR9 sample, Thamnos was identified in action-energy space with mean orbital parameters $\langle e\rangle = 0.51\pm0.15$, $\langle r_{\rm apo}\rangle = 9.87\pm1.8\,\mathrm{kpc}$, $\langle r_{\rm peri}\rangle = 3.31\pm0.9\,\mathrm{kpc}$, and $\langle z_{\max}\rangle = 3.70\pm1.3\,\mathrm{kpc}$ [2312.07825]. A separate LAMOST + Gaia study found a larger Thamnos sample with mean eccentricity $0.47$ and mean $|Z|_{\max}=2.73\,\mathrm{kpc}$, again reinforcing the picture of a moderately eccentric, low-energy retrograde population [2407.14508].

## 2. Selection functions and identification methodologies

Although the physical interpretation of Thamnos differs across papers, its operational definition is usually tied to a bounded region of orbital phase space. A common APOGEE-based selection adopts
$$
-1.8 < E < -1.6 \quad (\times 10^5\,\mathrm{km^2\,s^{-2}}),\qquad L_z<0,\qquad e<0.7,
$$
which yielded $N_{\rm Thamnos}=121$ stars in one chemical-characterization study and was reused in later contamination analyses [2204.04233; 2509.13408].

Different groups have recovered Thamnos with distinct clustering formalisms. These include Shared Nearest Neighbor clustering in $(J_R,J_\phi,J_z,E)$ space [2312.07825], the StarGO self-organizing-map algorithm in the seven-dimensional space $\{[\mathrm{Fe/H}],[\mathrm{Mg/Fe}],[\mathrm{Al/Fe}],L_z,J_R,J_z,E\}$ [2504.01536], t-SNE chemo-kinematic tagging based on $\{E,J_r,J_z,L_z,r_{\rm apo},r_{\rm peri},e,[\mathrm{Fe/H}]\}$ [2510.26952], wavelet transforms in $(L_n,\sqrt{J_r})$ space [2602.19647], CLiMB novelty detection for RR Lyrae in $(E,L_z,L_\perp)$ [2602.19858], and GS3 Hunter, which combines Siamese neural networks with K-means [2512.18693].

Representative definitions illustrate both the convergence and the non-uniformity of the literature:

| Study | Phase-space definition | Reported outcome |
|---|---|---|
| Koppelman et al. | Thamnos 1: $E/10^5\in[-1.65,-1.45]$, $\mathrm{circ}<-0.75$; Thamnos 2: $E/10^5\in[-1.80,-1.60]$, $-0.75<\mathrm{circ}<-0.40$ | Two low-energy retrograde clumps [1909.08924] |
| APOGEE chemical characterization | $-1.8<E<-1.6$, $L_z<0$, $e<0.7$ | $121$ stars [2204.04233] |
| GALAH DR4 wavelet analysis | Thamnos 1 centered at $\sqrt{J_r}=13.5\pm4.0$, $L_n=-600\pm350$; Thamnos 2 at $\sqrt{J_r}=8.6\pm3.5$, $L_n=-1184\pm175$ | Two retrograde groups [2602.19647] |
| RR Lyrae CLiMB analysis | Semi-supervised ellipsoidal assignment in $(E,L_z,L_\perp)$ around labeled Thamnos centroids | Sample increased from $26$ to $48$ RR Lyrae [2602.19858] |

The sign convention for retrograde motion is not uniform across all summaries. Some papers use the standard $L_z<0$ or $v_\phi<0$ for retrograde motion, whereas others report the same substructure with opposite-sign conventions for $J_\phi$ or $L_z$ [1909.08924; 2409.04529]. This suggests that cross-study comparisons of Thamnos must track the adopted convention explicitly rather than rely on sign alone.

## 3. Chemical-abundance patterns and metallicity structure

The chemical characterization of Thamnos is one of the main reasons it has remained astrophysically significant. In APOGEE-based analyses, Thamnos typically appears as a metal-poor, $\alpha$-enhanced population distinct from GES and from disc-like structures. One study reported $\langle[\mathrm{Fe/H}]\rangle\approx -1.4$ with a range from $-1.8$ to $-1.2$, together with $\langle[\mathrm{Mg/Fe}]\rangle \approx +0.28\pm0.05$, $\langle[\mathrm{Si/Fe}]\rangle \approx +0.24\pm0.06$, $\langle[\mathrm{Ca/Fe}]\rangle \approx +0.20\pm0.06$, $\langle[\mathrm{Al/Fe}]\rangle\approx -0.15\pm0.08$, and $\langle[\mathrm{Ce/Fe}]\rangle\approx -0.05\pm0.10$ [2204.04233]. In the same work, comparison with GES at $[\mathrm{Fe/H}]_{\rm comp}=-1.3$ gave $\chi^2_{\rm Thamnos-GES}=74.9$ with $p_{\chi^2}<10^{-3}$, supporting chemical distinctness [2204.04233].

A separate high-resolution study of low-$E$ retrograde groups overlapping Thamnos found a more metal-poor dominant population. Its metallicity distribution was fit by two Gaussians with peaks at $[\mathrm{Fe/H}]_1=-2.17\pm0.03$ and $[\mathrm{Fe/H}]_2=-1.49\pm0.05$, where the lower-metallicity component was identified as the main Thamnos population [2601.09796]. The same analysis reported a flat $\alpha$-plateau with $\langle[\mathrm{Mg/Fe}]\rangle \approx +0.40$, $\langle[\mathrm{Ca/Fe}]\rangle \approx +0.35$, $\langle[\mathrm{Si/Fe}]\rangle \approx +0.43$, and $\langle[\mathrm{Ti/Fe}]\rangle \approx +0.38$, together with no evidence for an $\alpha$ knee [2601.09796].

Recent GALAH DR4 work further resolved the retrograde low-energy region into Thamnos 1 and Thamnos 2. In that analysis, kernel-density estimates gave metallicity peaks of $-1.17\pm0.18$ for Thamnos 1 and $-1.28\pm0.79$ for Thamnos 2, while median abundance differences indicated higher $[\alpha/\mathrm{Fe}]$ for Thamnos 2, with approximately $+0.30$ dex versus $+0.22$ dex for Thamnos 1 [2602.19647]. Thamnos 2 also showed stronger $r$-process signatures, with $[\mathrm{Sm/Fe}]\simeq +0.18$ versus $\simeq +0.10$ for Thamnos 1 [2602.19647].

Taken together, these results suggest two overlapping chemical descriptions of Thamnos. One is a kinematically selected envelope centered near $[\mathrm{Fe/H}] \sim -1.4$ to $-1.5$ with moderate $\alpha$ enhancement; the other is a very metal-poor component near $[\mathrm{Fe/H}] \sim -2.1$ that several later studies interpret as the genuine accreted core. This interpretation is strengthened by the repeated claim that no clear $\alpha$ knee is seen in the chemically cleaner, metal-poor component [2204.04233; 2601.09796].

## 4. Stellar ages, star-formation history, and progenitor chronology

Chronological studies consistently place Thamnos among the oldest identified halo building blocks, though the exact age scale depends on tracer selection and fitting methodology. In CMD-fitting work based on a local 5D Gaia DR3 sample, the stricter “Thamnos B” selection yielded a metallicity distribution extending from approximately $-2.5$ to $-1.3$ dex and a half-mass lookback time
$$
t_{50}^{\rm Thamnos}=12.3^{+0.3}_{-0.3}\,\mathrm{Gyr},
$$
with a steep early star-formation history that rapidly declined for younger ages [2408.13763].

A complementary H3 survey analysis of MSTO and subgiant stars reached an even earlier chronology. For 35 stars classified as Thamnos members, nearly all inferred ages exceeded $12\,\mathrm{Gyr}$, with a sharp peak near $13\,\mathrm{Gyr}$ [2409.04529]. Fitting a truncated-Gaussian star-formation history gave $\mu_{\rm SFH}=13.79\,^{+0.15}_{-0.23}\,\mathrm{Gyr}$, $\sigma_{\rm SFH}=0.31\,^{+0.24}_{-0.15}\,\mathrm{Gyr}$, and $f_{\rm out}=0.17\pm0.08$, with $84\%$ of the stellar mass assembled by $\tau_{84}\approx13.4\pm0.3\,\mathrm{Gyr}$ [2409.04529]. Under the assumption that SFH truncation marks accretion, that study inferred $t_{\rm merge}\approx13.4\pm0.3\,\mathrm{Gyr\ ago}$ and described Thamnos as the earliest identified accretion event in the metal-poor halo [2409.04529].

A broader age compilation based on LAMOST and Gaia found somewhat younger but still ancient values: mean age $11.60\,\mathrm{Gyr}$, mode $12.89\,\mathrm{Gyr}$, and peak metallicity near $-1.36$ dex for a 3006-star Thamnos sample [2407.14508]. This does not contradict the older determinations so much as indicate method dependence in age recovery and in sample contamination.

The progenitor mass inferred for Thamnos is generally small. The discovery paper estimated $M_\star \lesssim 5\times10^6\,M_\odot$ from mock dwarf contours in $E$–$L_z$ space [1909.08924], and later high-resolution work continued to describe the genuine accreted component as a tiny, metal-poor dwarf with $M_\star\lesssim5\times10^6\,M_\odot$ [2510.06332]. A plausible implication is that Thamnos represents an early, low-mass accretion event whose star formation was truncated quickly after infall.

## 5. Internal complexity, contamination, and relation to other systems

The main controversy surrounding Thamnos is whether it is a chemically coherent relic of a single minor merger or a mixed region in integral-of-motion space. Several recent studies favor the latter interpretation. A chemically driven Gaussian-mixture-model comparison of kinematically selected halo substructures found that Thamnos likely contains GSE and heated disc stars in significant amounts [2509.13408]. In that analysis, the fraction chemically compatible with GSE was reported as $48\pm5\%$ for one definition, corresponding to an absolute $\simeq60\%$, or $59\pm4\%$ under an alternative definition; the metal-poor-disc-compatible fraction was $19\pm4\%$; and among the remaining outliers, $\sim80\%$ aligned chemically with $\omega$ Cen [2509.13408].

Dedicated UVES spectroscopy sharpened the same point. In a sample of 140 Thamnos candidates, the metallicity distribution was best fit by two Gaussians: a “true Thamnos” component with mean $-2.15$ and $\sigma=0.45$, comprising $22\%$ of the area, and a contaminant component with mean $-1.35$ and $\sigma=0.25$, comprising $78\%$ of the area [2510.06332]. That work concluded that the dynamically selected Thamnos region is dominated by in situ halo stars, with a smaller metal-poor accreted component and residual GSE contamination below $20\%$ [2510.06332].

RR Lyrae analyses point in the same direction. In a CLiMB-based Gaia DR3 study, the Thamnos RR Lyrae metallicity distribution was explicitly bimodal: a metal-poor peak with $23$ stars at $\langle[\mathrm{Fe/H}]\rangle=-1.94$ dex and $\sigma=0.20$ dex, and a metal-rich peak with $10$ stars at $\langle[\mathrm{Fe/H}]\rangle=-1.35$ dex and $\sigma=0.12$ dex [2602.19858]. The metal-poor peak was interpreted as the genuine Thamnos population, whereas the metal-rich peak was attributed to contamination from GSE and Splash/Aurora-like stars [2602.19858].

Relations to $\omega$ Cen have also become part of the modern discussion. A chemo-kinematic tagging study recovered a small Thamnos cluster containing NGC 288, NGC 5139 ($\omega$ Cen), the NGC 288 stream, and the Fimbulthul stream, and argued that the chemo-dynamic properties of $\omega$ Cen are consistent with a common accretion with Thamnos [2510.26952]. This does not establish identity between the two systems, but it does reinforce the idea that the standard Thamnos region may host multiple lineages of accreted debris.

An earlier chemodynamical study had already hinted at a link between Thamnos and GSE. In a broad retrograde sample, stars selected to probe Thamnos and GSE were different in overall metallicity and $[\mathrm{Mg/Fe}]$, but shared the same radial and vertical metallicity gradients, the same positive $\partial v_\phi/\partial[\mathrm{M/H}]\simeq+25\pm5\,\mathrm{km\,s^{-1}\,dex^{-1}}$, and the same eccentricity-metallicity trend [2009.07192]. This suggests that even if Thamnos is a distinct structure, its observed local manifestation overlaps dynamically with debris from other early halo-building events.

## 6. Specialized tracers and uncommon abundance signatures

Beyond ordinary giant-star samples, Thamnos has been traced with chemically distinctive populations. In the R-Process Alliance extended sample of $r$-process-enhanced halo stars, Thamnos was re-identified as CDTG-22, containing $8$ member stars with $\langle[\mathrm{Fe/H}]\rangle=-2.06\pm0.58$, $\langle[\mathrm{Mg/Fe}]\rangle=+0.54\pm0.16$, and $\langle[\mathrm{Eu/Fe}]\rangle=+0.54\pm0.21$ [2208.09712]. That work interpreted Thamnos as a chemically old, $r$-process self-enriched system, consistent with early enrichment before disruption [2208.09712].

An even more unusual result concerns Thamnos-2 and beryllium. A recent abundance study of nine stars associated with Thamnos found four new Be-rich members and showed that the previously known Be-rich stars HD 106038 and HD 132475 are also consistent with Thamnos membership [2607.08179]. In that sample, all currently known Be-rich stars were associated with Thamnos-2, and the Be enhancement correlated tightly with silicon and neutron-capture elements, including
$$
A(\mathrm{Be})=(1.05\pm0.15)\,[\mathrm{Si/Fe}] + 0.40,\qquad r=+0.93,
$$
and
$$
A(\mathrm{Be})=(0.98\pm0.20)\,[\mathrm{Ba/Fe}] + 0.35,\qquad r=+0.88
$$
[2607.08179]. The authors interpreted this pattern as the imprint of a rare spallation event, plausibly associated with a hypernova in the Thamnos-2 progenitor [2607.08179].

Thamnos has also been extended into the substellar regime. A survey of cold subdwarfs reported that the extreme L subdwarf 2MASS J053253.46+824646.5 and the mild T subdwarf CWISE J113010.07+313944.7 may be part of the Thamnos population on the basis of full 3D LSR velocities: for J0532+8246, $U=-55\pm2$, $V=-328\pm2$, $W=+66\pm2\,\mathrm{km\,s^{-1}}$; for J1130+3139, $U=-82\pm8$, $V=-265\pm24$, $W=+32\pm8\,\mathrm{km\,s^{-1}}$ [2411.01378]. This broadens the empirical scope of Thamnos beyond traditional stellar tracers.

The cumulative literature therefore presents Thamnos as both a named substructure and an evolving interpretive problem. Its repeated recovery in independent surveys supports the reality of a low-energy retrograde halo component, while its chemically resolved bimodality and cross-contamination with GSE, in situ halo populations, and possibly $\omega$ Cen-related debris indicate that “Thamnos” often denotes a structured region of phase space rather than a single mono-chemical population [1909.08924; 2510.06332]. The strongest common ground across recent work is that a genuinely ancient, metal-poor accreted component is present within that region, and that this component records one of the earliest low-mass mergers contributing to the Milky Way’s inner halo [2408.13763].

Source: https://www.emergentmind.com/topics/thamnos