Geometry of the Milky Way's dark matter from dynamical models of the tilted stellar halo

The shape and orientation of the Milky Way's dark matter halo remain poorly constrained. Observations of the accreted stellar halo show that it is triaxial and tilted with respect to the disc. If this configuration is long-lived, it can be used to place constraints on the shape and orientation of the dark matter halo that can support it close to steady state. We fit equilibrium orbit-superposition (Schwarzschild) models to the stellar halo in a realistic Milky Way potential with a tilted dark matter halo, assuming that the long axes of each halo and the disc normal are coplanar. These models are matched to parametric density fits and velocity anisotropy measurements of Gaia Sausage-Enceladus (GSE) stars at radii r∈[6,60] kpc. The observations are consistent with a (near-)prolate dark matter halo whose density has a short-to-long axis ratio of qdm=0.87-0.09+0.05. The long axis is inclined at an angle of βdm=43₋₈⁺²²\,^∘ to the disc plane, which exceeds the stellar halo tilt by ≈18^∘. Spherical haloes cannot support the observed structure of the GSE in equilibrium. The best-fitting dynamical GSE model has a radius-dependent shape and orientation; between radii of 6 and 60 kpc the tilt increases from β_*(r)≈10^∘ to ≈35^∘. Our model provides a good fit to the observed triaxial structure and dynamics of the GSE. It is therefore an excellent source of realistic initial conditions for simulations of the halo, such as for investigating perturbations from satellites or the Galactic bar.