Does it matter? A more careful treatment of density fluctuations in 21-cm simulations

The cosmological 21-cm signal is sourced from hyperfine transitions in neutral hydrogen atoms. Yet, although hydrogen atoms belong to the sub-class of baryonic matter, it is customary to simplify their treatment and analyze the 21-cm signal as if all the matter in the Universe was in the form of collisionless cold dark matter (CDM). This is usually done by evolving the density field via a scale-independent growth factor (SIGF). In this work, we separate the baryons from CDM and evolve the two species with a proper scale-dependent growth factor (SDGF). By incorporating the SDGF in the 21cmFirstCLASS code, we demonstrate the effect that baryons and CDM have on the 21-cm signal at the linear dark ages epoch and the subsequent non-linear epochs of cosmic dawn and reionization. Our analysis shows that the baryonic nature of hydrogen cannot be ignored during the dark ages, and that non-linear effects in density-field evolution must be accounted for after stars have formed. Furthermore, we discuss how the 21-cm signal is modified at lower redshifts, where ground-based 21-cm interferometers are mostly sensitive, due to the choice of working with either the "linear" or "non-linear" matter over-density (that is, the over-density as computed from linear perturbation theory, and non-linear perturbation theory, respectively) in the extended Press-Schechter formalism. Our code is publicly available at https://github.com/jordanflitter/21cmFirstCLASS.