2024ApJ...976...54P

Thin stellar streams originating from globular clusters (GCs) are among the most sensitive tracers of low-mass dark matter subhalos. Joint analysis of the entire population of stellar streams will place the most robust constraints on the dark matter subhalo mass function, and therefore the nature of dark matter. Here we use a hierarchical model of GC formation to forecast the total number, masses, and radial distribution of dissolved GC in Milky Way-like galaxies. Furthermore, we generate mock stellar streams from these progenitors' orbital histories taking into account the clusters' formation and accretion times, mass, and metallicity. Out of ∼10,000 clusters more massive than 10⁴M_⊙, ∼9000 dissolved in the central bulge and are fully phase mixed at the present, while the remaining ∼1000 survive as coherent stellar streams. This suggests that the current census of ∼80 GC streams in the Milky Way is severely incomplete. Beyond 15 kpc from the Galactic center we are missing ∼100 streams, of which the vast majority are from accreted GCs. Deep Rubin photometry (g ≲ 27.5) would be able to detect these streams, even the most distant ones beyond >75 kpc. We also find that M31 will have an abundance of streams at galactocentric radii of 30-100 kpc. We conclude that future surveys will find a multitude of stellar streams from GCs, which can be used for dark matter subhalo searches.