HI within and around observed and simulated galaxy discs -- Comparing MeerKAT observations with mock data from TNG50 and FIRE-2

Atomic hydrogen (HI) is an ideal tracer of gas flows in and around galaxies, and it is uniquely observable in the nearby Universe. Here we make use of wide-field (~1 square degree), spatially resolved (down to 22"), high-sensitivity (~$10^{18}$ cm$^{-2}$) HI observations of 5 nearby galaxies with stellar mass of $5\times10^{10}$ M$_\odot$, taken with the MeerKAT radio telescope. Four of these were observed as part of the MHONGOOSE survey. We characterise their main HI properties and compare these with synthetic HI data from a sample of 25 similarly massive star-forming galaxies from the TNG50 (20) and FIRE-2 (5) suites of cosmological hydrodynamical simulations. Globally, the simulated systems have HI and molecular hydrogen (H$_2$) masses in good agreement with the observations, but only when the H$_2$ recipe of Blitz & Rosolowsky (2006) is employed. The other recipes that we tested overestimate the H$_2$-to-HI mass fraction by up to an order of magnitude. On a local scale, we find two main discrepancies between observed and simulated data. First, the simulated galaxies show a more irregular HI morphology than the observed ones due to the presence of HI with column density $<10^{20}$ cm$^{-2}$ up to ~100 kpc from the galaxy centre, in spite of the fact that they inhabit more isolated environments than the observed targets. Second, the simulated galaxies and in particular those from the FIRE-2 suite, feature more complex and overall broader HI line profiles than the observed ones. We interpret this as being due to the combined effect of stellar feedback and gas accretion, which lead to a large-scale gas circulation that is more vigorous than in the observed galaxies. Our results indicate that, with respect to the simulations, gentler processes of gas inflows and outflows are at work in the nearby Universe, leading to more regular and less turbulent HI discs.