Weighing Neutrinos with Cosmic Neutral Hydrogen

We investigate the signatures left by massive neutrinos on the spatial distribution of neutral hydrogen (H I) in the post-reionization era by running hydrodynamic simulations that include massive neutrinos as additional collisionless particles. We find that halos in massive/massless neutrino cosmologies host a similar amount of neutral hydrogen, although for a fixed halo mass, on average, the H I mass increases with the sum of the neutrino masses. Our results show that H I is more strongly clustered in cosmologies with massive neutrinos, while its abundance, Omega(H I) (z), is lower. These effects arise mainly from the impact of massive neutrinos on cosmology: they suppress both the amplitude of the matter power spectrum on small scales and the abundance of dark matter halos. Modeling the H I distribution with hydrodynamic simulations at z > 3 and a simple analytic model at z < 3, we use the Fisher matrix formalism to conservatively forecast the constraints that Phase 1 of the Square Kilometre Array will place on the sum of neutrino masses, M-nu = Sigma m(nu). We find that with 10,000 hr of interferometric observations at 3 less than or similar to z less than or similar to 6 from a deep and narrow survey with SKA1-LOW, the sum of the neutrino masses can be measured with an error sigma(M-nu) less than or similar to 0.3 eV (95% CL). Similar constraints can be obtained with a wide and deep SKA1-MID survey at z less than or similar to 3, using the single-dish mode. By combining data from MID, LOW, and Planck, plus priors on cosmological parameters from a Stage IV spectroscopic galaxy survey, the sum of the neutrino masses can be determined with an error sigma(M-nu) similar or equal to 0.06 eV (95% CL).