Galactic Substructure and Dark Matter Annihilation in the Milky Way Halo

We study the effects of substructure on the rate of dark-matter annihilation in the Galactic halo. We use an analytic model for substructure that can extend numerical simulation results to scales too small to be resolved by the simulations. We first calibrate the analytic model to numerical simulations, and then determine the annihilation boost factor, for standard WIMP models as well as those with Sommerfeld (or other) enhancements, as a function of Galactocentric radius in the Milky Way. We provide an estimate of the dependence of the gamma-ray intensity of WIMP annihilation as a function of angular distance from the Galactic center. This methodology, coupled with future numerical simulation results can be a powerful tool that can be used to constrain WIMP properties using Fermi all-sky data.

💡 Research Summary

The paper investigates how dark‑matter substructure within the Milky Way halo influences the annihilation rate of weakly interacting massive particles (WIMPs). Because state‑of‑the‑art N‑body simulations cannot resolve subhalos below roughly 10⁶ M⊙, the authors develop an analytic framework that extrapolates the subhalo mass function, spatial distribution, and internal density profiles down to the theoretically motivated cutoff of ~10⁻⁶ M⊙. The model is calibrated against high‑resolution simulations such as Via Lactea II and Aquarius, ensuring that the mass‑function slope (α≈1.9), concentration‑mass relation, and tidal stripping near the Galactic center are faithfully reproduced.

With the calibrated model in hand, the authors compute the “boost factor” – the ratio of the annihilation luminosity including subhalos to that of a smooth halo – as a function of Galactocentric radius. Two annihilation scenarios are considered: (i) the standard velocity‑independent s‑wave cross‑section σ₀, and (ii) a Sommerfeld‑enhanced cross‑section that scales as 1/v or 1/v² at low relative velocities. To evaluate the velocity distribution inside each subhalo, they adopt a virial relation linking subhalo mass and radius to a characteristic velocity dispersion, and they incorporate the effect of tidal stripping, which reduces subhalo densities in the inner ≲1 kpc.

The results show that for a standard WIMP the boost factor is modest: it is close to unity within the inner kiloparsec (where subhalos are heavily disrupted) and rises to 5–20 at distances of 10–30 kpc, where subhalos survive relatively intact. When Sommerfeld enhancement is included, the boost can become dramatically larger in the outer halo because the typical subhalo velocity dispersion drops to ≲10⁻⁴ c, driving the cross‑section upward. In this case the boost factor reaches 30–200 depending on the exact velocity scaling and the assumed saturation of the enhancement.

The authors then translate the radial boost profile into an observable gamma‑ray intensity map. By integrating the annihilation emissivity along the line of sight for different angular separations from the Galactic Center, they produce predicted gamma‑ray fluxes as a function of viewing angle. The angular profile exhibits a steep rise toward the center, reflecting the combined effect of the smooth halo density cusp and the increasing contribution of subhalos at larger radii. These predictions are directly comparable to the all‑sky data from the Fermi Large Area Telescope. The paper demonstrates that, with a reliable substructure model, one can use the angular distribution of the diffuse gamma‑ray background to place constraints on WIMP mass, annihilation cross‑section, and any velocity‑dependent enhancement such as Sommerfeld effects.

Finally, the authors argue that their analytic methodology is flexible: as future simulations push to ever‑smaller mass scales and better resolve tidal disruption, the model parameters can be updated, sharpening the connection between substructure physics and indirect dark‑matter searches. This synergy between theory, simulation, and observation offers a powerful avenue for probing the particle nature of dark matter using existing and forthcoming gamma‑ray data.