2025A&A...700A.131R

The accretion and processing of neutral and ionized gas play substantial roles in the evolution of the Milky Way. From the position of the Sun, circumgalactic gas flows in the Milky Way halo are known to span a large range of radial velocities, but the complex kinematics of the circumgalactic medium (CGM) cannot be fully reconstructed from observations because of the blending with foreground interstellar gas in the Milky Way disk. For this paper we used three zoom-in magnetohydrodynamic simulations of the Milky Way and the Local Group from the HESTIA project to systematically investigate the radial velocity distribution of neutral hydrogen (H I) clouds in the CGM in the (simulated) Local Standard of Rest (LSR) velocity frame. Our three simulations, which exhibit substantial differences in their global CGM properties, reveal that 48-65 percent of the extraplanar H I at z > 2 kpc above the plane is confined to a velocity range |v_LSR| ≤ 100 km s^–1, implying that the gas is (at least partly) corotating with the underlying disk. In the two most realistic Milky Way realizations, the CGM velocity distribution is skewed toward negative velocities (in particular for H I clouds at vertical distances z > 10 kpc), indicating a net accretion of neutral gas. These results are in line with the statistics from UV absorption-line measurements of the Milky Way CGM, and we also find broad agreement with the Illustris TNG50 simulation. Our study supports a scenario in which a substantial fraction of the Milky Way's CGM resides close to the disk at |v_LSR| ≤ 100 km s^–1, where it is hiding from observations as its spectral signatures are covered by foreground interstellar gas features. We furthermore find that 97 percent of the clumps live in the Milky Way halo and are not associated with satellite galaxies. The clumps are magnetized with a magnetic pressure often dominating over the thermal pressure.