skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 16, 2026

Title: Semianalytic Modeling of Dark Matter Subhalo Encounters with Thin Stellar Streams: Statistical Predictions for GD-1-like Streams in Cold Dark Matter
Abstract Stellar streams from disrupted globular clusters are dynamically cold structures that are sensitive to perturbations from dark matter subhalos, allowing them in principle to trace the dark matter substructure in the Milky Way. We model, within the context of Λ cold dark matter, the likelihood of dark matter subhalos to produce a significant feature in a GD-1-like stream and analyze the properties of such subhalos. We generate many realizations of the subhalo population within a Milky Way mass host halo using the semianalytic codeSatGen, accounting for effects such as tidal stripping and dynamical friction. The subhalo distributions are combined with a GD-1-like stream model, and the impact of subhalos that pass close to the stream are modeled withGala. We find that subhalos with masses in the range 2 × 106M–108Mat the time of the stream–subhalo encounter, corresponding to masses of about 2 × 107M–109Mat the time of infall, are the likeliest to produce gaps in a GD-1-like stream. We find that gaps occur on average ∼3 times per realization of the host system. These gaps have typical widths of ∼(5–27)° and fractional underdensities of ∼(10–30)%, with larger gaps being caused by heavier subhalos. The stream–subhalo encounters responsible for these have impact parameters (0.1–1.5) kpc and relative velocities ∼(200–410) km s−1. We also investigate the effects of increasing the host-halo mass on the gap properties and formation rate.  more » « less
Award ID(s):
2210533
PAR ID:
10646498
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
American Astronomical Society
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
991
Issue:
1
ISSN:
0004-637X
Page Range / eLocation ID:
66
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT A variety of observational campaigns seek to test dark matter models by measuring dark matter subhaloes at low masses. Despite their predicted lack of stars, these subhaloes may be detectable through gravitational lensing or via their gravitational perturbations on stellar streams. To set measurable expectations for subhalo populations within Lambda cold dark matter, we examine 11 Milky Way (MW)-mass haloes from the FIRE-2 baryonic simulations, quantifying the counts and orbital fluxes for subhaloes with properties relevant to stellar stream interactions: masses down to $$10^{6}\, \text{M}_\odot$$, distances ≲50 kpc of the galactic centre, across z = 0 − 1 (tlookback = 0–8 Gyr). We provide fits to our results and their dependence on subhalo mass, distance, and lookback time, for use in (semi)analytical models. A typical MW-mass halo contains ≈16 subhaloes $$\gt 10^{7}\, \text{M}_\odot$$ (≈1 subhalo $$\gt 10^{8}\, \text{M}_\odot$$) within 50 kpc at z ≈ 0. We compare our results with dark matter-only versions of the same simulations: because they lack a central galaxy potential, they overpredict subhalo counts by 2–10×, more so at smaller distances. Subhalo counts around a given MW-mass galaxy declined over time, being ≈10× higher at z = 1 than at z ≈ 0. Subhaloes have nearly isotropic orbital velocity distributions at z ≈ 0. Across our simulations, we also identified 4 analogues of Large Magellanic Cloud satellite passages; these analogues enhance subhalo counts by 1.4–2.1 times, significantly increasing the expected subhalo population around the MW today. Our results imply an interaction rate of ∼5 per Gyr for a stream like GD-1, sufficient to make subhalo–stream interactions a promising method of measuring dark subhaloes. 
    more » « less
  2. Abstract TheFermiLarge Area Telescope (Fermi-LAT) has been widely used to search for Weakly Interacting Massive Particle (WIMP) dark matter signals due to its unparalleled sensitivity in the GeV energy band. The leading constraints for WIMP byFermi-LAT are obtained from the analyses of dwarf spheroidal galaxies within the Local Group, which are compelling targets for dark matter searches due to their relatively low astrophysical backgrounds and high dark matter content. In the meantime, the search for heavy dark matter with masses above TeV remains a compelling and relatively unexplored frontier. In this study, we utilize 14-yearFermi-LAT data to search for dark matter annihilation and decay signals in 8 classical dwarf spheroidal galaxies within the Local Group. We consider secondary emission caused by electromagnetic cascades of prompt gamma rays and electrons/positrons from dark matter, which enables us to extend the search withFermi-LAT to heavier dark matter cases. We also update the dark matter subhalo model with informative priors respecting the fact that they reside in subhalos of our Milky Way halo aiming to enhance the robustness of our results. We place constraints on dark matter annihilation cross section and decay lifetime for dark matter masses ranging from 103GeV to 1011GeV, where our limits are more stringent than those obtained by many other high-energy gamma-ray instruments. 
    more » « less
  3. Abstract The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)–mass cosmological zoom-in simulations and an empirical galaxy–halo connection model, and assess sensitivity to galaxy formation and DM signals when marginalizing over galaxy–halo connection uncertainties. We find that a survey of all satellites around one MW-mass host can constrain a galaxy formation cutoff at peak virial masses of M 50 = 10 8 M at the 1σlevel; however, a tail toward low M 50 prevents a 2σmeasurement. In this scenario, combining hosts with differing bright satellite abundances significantly reduces uncertainties on M 50 at the 1σlevel, but the 2σtail toward low M 50 persists. We project that observations of one (two) complete satellite populations can constrain warm DM models withmWDM≈ 10 keV (20 keV). Subhalo mass function (SHMF) suppression can be constrained to ≈70%, 60%, and 50% that in cold dark matter (CDM) at peak virial masses of 108, 109, and 1010M, respectively; SHMF enhancement constraints are weaker (≈20, 4, and 2 times that in CDM, respectively) due to galaxy–halo connection degeneracies. These results motivate searches for faint dwarf galaxies beyond the MW and indicate that ongoing missions like Euclid and upcoming facilities including the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope will probe new galaxy formation and DM physics. 
    more » « less
  4. Abstract We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (fdiss= 0.75). In this paper, we examine the properties of dwarf galaxies withM*∼ 105–109Min both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm2g−1. On the contrary, models with (σ/m) = 10 cm2g−1produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm2g−1are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hirotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm2g−1,fdiss= 0.75 are in tension in massive dwarfs (Mhalo∼ 1011M) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures. 
    more » « less
  5. A<sc>bstract</sc> Cosmologically stable, light particles that came into thermal contact with the Standard Model in the early universe may persist today as a form of hot dark matter. For relics with masses in the eV range, their role in structure formation depends critically on their mass. We trace the evolution of such hot relics and derive their density profiles around cold dark matter halos, introducing a framework for theirindirect detection. Applying this framework to axions — a natural candidate for a particle that can reach thermal equilibrium with the Standard Model in the early universe and capable of decaying into two photons — we establish stringent limits on the axion-photon couplinggusing current observations of dwarf galaxies, the Milky Way halo, and galaxy clusters. Our results set new bounds on hot axions in the$$ \mathcal{O}\left(1-10\right) $$ O 1 10 eV range. 
    more » « less