Abstract
We estimate the 3D density profile of the Galactic dark matter (DM) halo within r less than or similar to 30kpc from the Galactic centre by using the astrometric data for halo RR Lyrae stars from Gaia DR2. We model both the stellar halo distribution function and the Galactic potential, fully taking into account the survey selection function, the observational errors, and the missing line-of-sight velocity data for RR Lyrae stars. With a Bayesian method, we infer the model parameters, including the density flattening of the DM halo q, which is assumed to be constant as a function of radius. We find that 99 per cent of the posterior distribution of q is located at q > 0.963, which strongly disfavours a flattened DM halo. We cannot draw any conclusions as to whether the Galactic DM halo at is prolate, because we restrict ourselves to axisymmetric oblate halo models with q <= 1. Our DM density profile might be biased especially in the inner few kpc, due to the uncertainty in the baryonic distribution. Our result is in tension with predictions from cosmological hydrodynamical simulations that advocate more oblate (q similar to 0.8 +/- 0.15) DM haloes within of the virial radius for Milky-Way-sized galaxies. An alternative possibility, based on our validation tests with a cosmological simulation, is that the true value q of the Galactic halo could be consistent with cosmological simulations but that disequilibrium in the Milky Way potential is inflating our measurement of q by 0.1-0.2. As a by-product, our model constrains the DM density in the Solar neighbourhood to be rho(DM,circle dot) = (9.01(-2.0)(+0.18)) x 10(-3) M(circle dot)pc(-3) = 0.342(-0.007)(+0.007) M GeVcm(-3), consistent with other recent measurements.