Dark Matter Constraints from a Unified Analysis of Strong Gravitational Lenses and Milky Way Satellite Galaxies

TL;DR

This paper combines strong gravitational lensing and Milky Way satellite data to improve constraints on warm dark matter properties, accounting for model degeneracies and covariances, and demonstrates a significant tightening of existing limits.

Contribution

It presents a novel multidimensional analysis that jointly infers dark matter substructure properties from heterogeneous data sources, advancing the methodology for small-scale structure studies.

Findings

Sets new upper limit on the half-mode mass, $M_{hm}<10^{7.0} M_{\odot}$

Disfavors warm dark matter particle mass below 9.7 keV at 95% confidence

Improves existing limits on $m_{WDM}$ by approximately 30%

Abstract

Joint analyses of small-scale cosmological structure probes are relatively unexplored and promise to advance measurements of microphysical dark matter properties using heterogeneous data. Here, we present a multidimensional analysis of dark matter substructure using strong gravitational lenses and the Milky Way (MW) satellite galaxy population, accounting for degeneracies in model predictions and using covariances in the constraining power of these individual probes for the first time. We simultaneously infer the projected subhalo number density and the half-mode mass describing the suppression of the subhalo mass function in thermal relic warm dark matter (WDM), $M_{\mathrm{hm}}$, using the semianalytic model $\mathrm{\texttt{Galacticus}}$ to connect the subhalo population inferred from MW satellite observations to the strong lensing host halo mass and redshift regime. Combining MW satellite and strong lensing posteriors in this parameter space yields $M_{\mathrm{hm}}<10^{7.0}\ M_{\mathrm{\odot}}$ (WDM particle mass $m_{\mathrm{WDM}}>9.7\ \mathrm{keV}$) at $95\%$ confidence and disfavors $M_{\mathrm{hm}}=10^{7.4}\ M_{\mathrm{\odot}}$ ($m_{\mathrm{WDM}}=7.4\ \mathrm{keV}$) with a 20:1 marginal likelihood ratio, improving limits on $m_{\mathrm{WDM}}$ set by the two methods independently by $\sim 30\%$. These results are marginalized over the line-of-sight contribution to the strong lensing signal, the mass of the MW host halo, and the efficiency of subhalo disruption due to baryons and are robust to differences in the disruption efficiency between the MW and strong lensing regimes at the $\sim 10\%$ level. This work paves the way for unified analyses of next-generation small-scale structure measurements covering a wide range of scales and redshifts.