The local PNG bias of neutral Hydrogen, HI

Abstract

We use separate universe simulations with the IllustrisTNG galaxy formation model to predict the local PNG bias parameters bΦ and bΦδ of atomic neutral hydrogen, HI. These parameters and their relation to the linear density bias parameter b 1 play a key role in observational constraints of the local PNG parameter f NL using the HI power spectrum and bispectrum. Our results show that the popular calculation based on the universality of the halo mass function overpredicts the bΦ (b 1) and bΦδ (b 1) relations measured in the simulations. In particular, our results show that at z ≲ 1 the HI power spectrum is more sensitive to f NL compared to previously thought (bΦ is more negative), but is less sensitive at other epochs (bΦ is less positive). We discuss how this can be explained by the competition of physical effects such as that large-scale gravitational potentials with local PNG (i) accelerate the conversion of hydrogen to heavy elements by star formation, (ii) enhance the effects of baryonic feedback that eject the gas to regions more exposed to ionizing radiation, and (iii) promote the formation of denser structures that shield the HI more efficiently. Our numerical results can be used to revise existing forecast studies on f NL using 21 cm line-intensity mapping data. Despite this first step towards predictions for the local PNG bias parameters of HI, we emphasize that more work is needed to assess their sensitivity on the assumed galaxy formation physics and HI modeling strategy.