Disentangling the Distant Stellar Halo Using K Giants in DESI Year 3 Data

We present a sample of 88,959 K giants from Dark Energy Spectroscopic Instrument Milky Way Survey Year 3 data, which we use to characterize the chemo-dynamical properties of the stellar halo at Galactocentric distances of 12 to ∼100 kpc. Using HDBSCAN, we identify five prominent stellar halo substructures: Aleph, the Sagittarius stream, Gaia-Sausage-Enceladus (GSE), Cetus-Palca, and the Orphan–Chenab stream. We present the properties of each of these structures as they appear in our catalog, and examine how uncertainties on distance affect the characterization of substructure with this approach. We also examine regions associated with previously reported overdensities (such as the Virgo Overdensities and the Sagittarius spur) that we do not recover with HDBSCAN. The size and distance range of our catalog allows us to explore in detail the residual stellar halo, comprising stars that we do not associate with any substructure. We find that samples of ∼2000 outer halo stars with both highly prograde and highly retrograde angular momenta have similar metallicity distribution functions (MDFs), which do not resemble the MDFs of either GSE or Sagittarius. Both the prograde and retrograde residual halo MDFs are bimodal, with a metal-poor peak at [Fe/H] ∼ −2 and a metal-rich peak at [Fe/H] ∼ −1.3 (prograde) or −1.5 (retrograde). The MDF for lower-angular-momentum residual halo K giants does not show clear evidence for a metal-poor peak, and broadly resembles the MDF of GSE, even at much lower binding energies than GSE itself. We discuss possible interpretations of these findings for GSE accretion scenarios.