Implications of a contracted dark matter halo for the Milky Way's inferred virial mass

TL;DR

This study assesses how well dark matter halo properties of the Milky Way can be inferred from limited radial data, revealing biases in common models and proposing a contracted halo model for unbiased estimates.

Contribution

The paper demonstrates that standard parametric models underestimate halo mass when fitting only inner regions and shows that a contracted halo model provides more accurate, unbiased mass estimates across all radii.

Findings

NFW profiles underestimate virial mass by ~2 when fitted to inner data.

Contracted halo model yields unbiased mass estimates across all radii.

Biases decrease with increasing radial coverage, especially beyond 50 kpc.

Abstract

We investigate how reliably the global properties of Milky Way-mass dark matter haloes can be recovered from dynamical data over a limited radial range, particularly $\lesssim 30~\mathrm{kpc}$ where observations are most sensitive but baryonic processes modify the halo structure. Using the ARTEMIS simulations, which produce varying degrees of baryon-induced contraction, we fit dark matter profiles over restricted radial ranges using commonly adopted parametric models. Assuming negligible observational uncertainties allows the systematic errors from these choices to be isolated. When fits are confined to inner radii, an NFW profile underestimates the virial mass by a factor of $\approx 2$ on average ($\approx 4$ for some systems), and the concentration by a factor of $\approx 2$. Einasto and generalised-NFW models provide excellent local fits but retain similar global biases. In contrast, the contracted halo prescription from Cautun et al. (2020) yields stable extrapolations and recovers unbiased halo mass estimates over all radii. The inferred mass improves systematically with increasing radial coverage, and tracers beyond $\gtrsim 50~\mathrm{kpc}$ largely eliminate the mean bias for all models. The local dark matter density at the Solar radius is recovered to within $\lesssim 5\%$ for all profiles other than NFW. These biases are sufficient to reconcile recent low Milky Way mass estimates derived from inner rotation-curve analyses with the canonical $\approx 10^{12}~\mathrm{M}_\odot$. We additionally find a halo-to-halo scatter of $\gtrsim 0.1$ dex ($\approx 25\%$) persists even under idealised conditions, setting a likely lower limit for the precision of halo mass estimates.