Several authors have recently applied Jeans modelling to -based datasets to infer the circular velocity curve for the Milky Way. These works have consistently found evidence for a continuous decline in the rotation curve beyond $ kpc $, which may indicate the existence of a light dark matter (DM) halo. Using a large sample of Gaia DR3 data, we aim to derive the rotation curve of the Milky Way using the Jeans equations, and to quantify the role of systematic effects, both in the data and those inherent to the Jeans methodology under the assumptions of axisymmetry and time independence. We used data from the Gaia DR3 radial velocity spectrometer sample, supplemented with distances inferred through Bayesian frameworks, to determine the radial variation of the second moments of the velocity distribution for stars close to the Galactic plane. We used these profiles to determine the rotation curve using the Jeans equations under the assumption of axisymmetry and explored how they vary with azimuth and position above and below the plane of the Galactic disc. We applied the same methodology to an N-body simulation of a Milky Way-like galaxy impacted by a satellite akin the Sagittarius dwarf, and to the Auriga suite of cosmological simulations. The circular velocity curve we infer for the Milky Way is consistent with previous findings out to $ 15$ kpc, where our statistics are robust. Due to the larger number of stars in our sample, we are able to reveal evidence of disequilibrium and deviations from axisymmetry closer in. For example, we find that the second moment of $V_R$ flattens out at $R kpc, and that the second moment of $V_ is different above and below the plane for $R kpc. Our exploration of the simulations indicates that these features are typical of galaxies that have been perturbed by external satellites. From the simulations, we also estimate that the difference between the true circular velocity curve and that inferred from Jeans equations can be as high as $15<!PCT!>$, but that it is likely of the order of $10<!PCT!>$ for the Milky Way. This is higher than the systematic uncertainties associated with the observations or those linked to most modelling assumptions when using the Jeans equations. However, if the density of the tracer population were truncated at large radii instead of being exponential as often assumed, this could lead to the erroneous conclusion of a steeply declining rotation curve. We find that steady-state axisymmetric Jeans modelling becomes less robust at large radii, indicating that particular caution must be exercised when interpreting the rotation curve inferred in those regions. A more careful and sophisticated approach may be necessary for precision measurements of the DM content of our Galaxy.
Read full abstract