The density of the Milky Way’s corona at z ≈ 1.6 through ram pressure stripping of the Draco dSph galaxy

Satellite galaxies within the Milky Way’s (MW’s) virial radius Rvir are typically devoid of cold gas due to ram pressure stripping by the MW’s corona. The density of this corona is poorly constrained today and essentially unconstrained in the past, but can be estimated using ram pressure stripping. In this paper, we probe the MW’s corona at z ≈ 1.6 using the Draco dwarf spheroidal galaxy. We assume that (i) Draco’s orbit is determined by its interaction with the MW, whose dark matter halo we evolve in time following cosmologically motivated prescriptions, (ii) Draco’s star formation was quenched by ram pressure stripping and (iii) the MW’s corona is approximately smooth, spherical, and in hydrostatic equilibrium. We used Gaia proper motions to set the initial conditions and Draco’s star formation history to estimate its past gas content. We found indications that Draco was stripped of its gas during the first pericentric passage. Using 3D hydrodynamical simulations at a resolution that enables us to resolve individual supernovae and assuming no tidal stripping, which we estimate to be a minor effect, we find a density of the MW corona ≥8 × 10−4 cm−3 at a radius ≈0.72Rvir . This provides evidence that the MW’s corona was already in place at z ≈ 1.6 and with a higher density than today. If isothermal, this corona would have contained all the baryons expected by the cosmological baryon fraction. Extrapolating to today shows good agreement with literature constraints if feedback has removed <~30 per cent of baryons accreted on to the halo.