2025ApJ...994..128O

Recent measurements of the Milky Way rotation curve found a sharp decline at around 15-20 kpc from the center of the Galaxy, suggesting that the Galactic dark matter halo is much less massive than predicted by other dynamical tracers. To address this tension, we study the validity of the assumptions made in calculating the Milky Way's rotation curve. To do so, we apply the Jeans equation, the current standard approach of measuring rotation curves, to three cosmological zoom-in simulations of Milky Way-like galaxies from the FIRE-2 Latte suite. Using synthetic Gaia surveys, we replicate the sample selection process and calculation employed in measuring the Milky Way rotation curve. We examine four failure modes of this calculation and find that the measured curves deviate from the true curve by 5%-20% rather than below 5%, as estimated by previous works. Interestingly, there is a large galaxy-to-galaxy variance, and different systematics dominate different galaxies. We rederive the Milky Way's dark matter density profile with the rotation curve while incorporating systematics from the simulations. The posterior distribution of the density profiles is consistent with a fiducial Navarro-Frenk-White (NFW) profile when assuming a generalized NFW profile for dark matter. We find that the virial mass is $7.3{2}_{-1.53}^{+1.98}\times 1{0}^{11}$M_⊙, consistent with other probes of the Milky Way's mass. However, we recommend that the field move away from relying solely on the rotation curve when studying the dark matter profile and adopt methods that incorporate additional probes and/or do not heavily depend on assumptions described in this study.