HIDES -- I. The population and diversity of HI-rich 'dark' galaxies in the Hestia and Auriga simulations

TL;DR

This study investigates the population and diversity of HI-rich 'dark' galaxies in the Hestia and Auriga simulations, revealing their characteristics, environmental effects, and implications for galaxy formation theories.

Contribution

It provides the first detailed analysis of HI-rich 'dark' galaxies in these simulations, highlighting their demographics, environmental influences, and a new empirical density fit.

Findings

89 HI-rich 'dark' galaxies identified, mostly with low star formation

Environmental effects like ram pressure stripping influence HI content

Large scatter in HI properties explained by environment and history

Abstract

We present our investigation of HI-rich 'Dark' galaxiEs in Simulations (HIDES), specifically using the Hestia and Auriga simulations in this work. We select galaxies that are faint ($M_g > -10$) and contain sufficient HI ($M_\mathrm{HI} > 10^5\,M_\odot$), and identify 89 such objects, only one of which is completely starless. Their demographics generally converge across simulations of different resolution, with $M_{200} \sim 10^{9.5}\,M_\odot$, $M_\mathrm{gas} \sim 10^{7.4}\,M_\odot$, $M_\mathrm{HI} \sim 10^{6.5}\,M_\odot$, $M_\mathrm{*} \sim 10^{5.6}\,M_\odot$, low gas metallicity, little or no current star formation, and a mean stellar age of $\sim$ 11 Gyr, and with some of them can survive in dense environments as close as $\sim$ 300 kpc from a Milky-Way mass neighbor. We find a large scatter in their HI density profiles and $M_\mathrm{HI} - M_\mathrm{*}$ relation, which cannot be fully explained by current halo mass or concentration, but can be attributed to ram pressure stripping in dense environments, past mergers, and stellar feedback. In particular, close encounters with massive halos and dense environments can reshape the HI content, which may explain the asymmetric HI map of an intriguing observed analogue, Cloud-9. An empirical fit, $n = 0.25 \left(d_\mathrm{MW}/{1\,\mathrm{Mpc}}\right)^{-1.4}\, \mathrm{Mpc}^{-3}$, based on their number density extended to 3.7 Mpc in constrained local volume simulations, is also provided to aid observational forecasts. We conclude that both mass assembly history and environmental history play a crucial role in the formation and subsequent diversity of these galaxies.