REVIEW
Weak Evolution of Cosmic Atomic Hydrogen over the Past 4.5 Billion Years
Astrophysics of Galaxies · Submitted 6 Jul 2026 · Published version · 93 references
Pith's verdict: The central qualitative result — that Ω_HI evolves much more weakly than the CSFRD — is robust and well-supported. The quantitative headline (factor 1.35, or 1.12 after systematics) is sensitive to an extrapolation below 10^9 M⊙ that the paper's own forward model cannot validate.
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
Cosmic atomic hydrogen barely declined while star formation plummeted
The reading
What carries the argument
The paper introduces a luminosity-bias correction method that models the HI gas fraction as a function of both r-band luminosity and stellar mass, f_HI(L_r|M_*), using an 8-parameter double-power-law fit. This correction accounts for the fact that flux-limited optical samples are biased toward more luminous (and gas-richer) galaxies at fixed mass. By measuring how f_HI depends on luminosity within each stellar-mass bin, the method reconstructs the intrinsic gas fraction that would be observed in a complete sample, without requiring mass-completeness cuts that would sacrifice the low-mass galaxy population.
Extended reading notes
Core claim
The cosmic atomic hydrogen density has remained nearly constant over the past 4.5 billion years, declining by at most a factor of 1.35 (or 1.12 after systematic corrections), while the cosmic star formation rate density fell by a factor of 2.46 over the same interval. This mismatch, present at fixed stellar mass across the galaxy population, rules out rapid HI depletion as the cause of declining star formation and points to the HI-to-H2 phase conversion or gas accretion regulation as the operative bottleneck.
Load-bearing premise
The luminosity bias correction assumes that the intrinsic distribution of r-band luminosity at fixed stellar mass does not evolve over 0 < z < 0.41, using the local (z < 0.05) sample as the reference for all redshifts. If this distribution actually evolves with redshift, the corrected HI gas fractions at z > 0.25 could be systematically offset.
Editorial extensions
If this is right
- Galaxy formation models must reproduce a baryon cycle where the atomic hydrogen reservoir stays nearly invariant while molecular gas and star formation decline substantially — a tighter constraint than matching Omega_HI alone.
- The bottleneck driving late-time cosmic star formation decline lies downstream of the atomic gas reservoir, in the atomic-to-molecular phase conversion or in the regulation of gas accretion onto galactic disks, not in the depletion of HI itself.
- Future 21-cm intensity mapping surveys can use the measured Omega_HI(z) relation, Omega_HI proportional to (1+z)^0.9, as a benchmark over 0 < z < 0.41.
- The f_HI(M_*) relation's near-invariance with redshift provides a stable calibration target for semi-analytic and hydrodynamical models aiming to predict gas partitioning across cosmic time.
Reading between the lines
- If the HI reservoir is stable while H2 declines, the star formation efficiency per unit molecular gas may itself be evolving, or the molecular gas depletion timescale may be lengthening — a question the paper does not directly address but that follows from the decoupling it establishes.
- The near-invariance of f_HI(M_*) suggests that gas accretion and gas loss are finely balanced across all stellar masses, which may point to a self-regulating equilibrium in the atomic gas phase that is insensitive to the processes driving the molecular phase decline.
- Extending the same f_HI(L_r|M_*) correction framework to higher redshifts (z > 0.4) with deeper surveys would test whether the weak evolution breaks down at earlier epochs when cosmic gas accretion rates were higher.
Original abstract
The cosmic star formation rate density (CSFRD) has declined sharply toward the present day, but the roles of the atomic and molecular gas reservoirs remain uncertain. We measure the cosmic HI density, $\Omega_{\mathrm{HI}}$, over $0<z<0.41$ by combining HI spectra from the Five-hundred-meter Aperture Spherical Telescope with optical spectroscopy from the Dark Energy Spectroscopic Instrument for $\sim2.5$ million galaxies across $\sim12,000\,{\rm deg}^2$. We measure a raw decrease in $\Omega_{\mathrm{HI}}$ by a factor of $1.35\pm0.10$ over the past 4.5 Gyr. Even after applying the conservative systematic corrections from our forward model, the inferred decline is only $1.12\pm0.10$ -- still far weaker than the CSFRD decline (a factor of 2.46). The molecular gas density, in contrast, is known to evolve more closely with star formation. At fixed stellar mass, the average HI gas fraction evolves by less than 0.2 dex, showing that the weak evolution is present across the galaxy population. These quantitative differences rule out rapid depletion of galaxy HI as the primary driver of the late-time CSFRD decline, and provide a stringent benchmark for models of gas accretion, phase conversion and star-formation regulation.
Figures
Figures reproduced from arXiv:2607.05326 by the authors.
References
How the literature uses it
1 arXiv paper cite this in 1 sentence
Papers citing it
Cites it as background 1 paper
-
Neutral atomic hydrogen () plays a fundamental role in galaxy evolution by providing the primary reservoir from which molecular gas and stars ultimately form (e.g., [ref] ).
FASHI DR2: A Catalog of 132 Low-Redshift HI 21 cm Absorption Systems · Introduction · 0.97
What it cites · 4 sentences without a confident label
Uses data or values from 3 papers
-
The 21 cm spectral data used in this work are drawn from the FAST All Sky survey (FASHI [ref] ), which employs FAST’s 19-beam receiver to scan the sky over $-14^∘< DEC < 64^∘$.
The FAST All Sky HI Survey DR2: the FASHI Catalog and the HI Mass Function · FASHI Observations and data · 0.83
-
The input data for CIGALE include optical and infrared photometry in the $g$, $r$, $z$, $W1$, and $W2$ bands, obtained from Data Release 9 of the DESI Legacy Imaging Surveys (LS-DR9 [ref] ), along with spectroscopic redshifts from the loa release of DESI DR2 [ref] .
DESI DR2 Results II: Measurements of Baryon Acoustic Oscillations and Cosmological Constraints · Calculations for $M_∗$, $L_r$, $M_HI$, and $f_HI$ · 0.99
-
The input data for CIGALE include optical and infrared photometry in the $g$, $r$, $z$, $W1$, and $W2$ bands, obtained from Data Release 9 of the DESI Legacy Imaging Surveys (LS-DR9 [ref] ), along with spectroscopic redshifts from the loa release of DESI DR2 [ref] .
DESI DR2 Results I: Baryon Acoustic Oscillations from the Lyman Alpha Forest · Calculations for $M_∗$, $L_r$, $M_HI$, and $f_HI$ · 0.99
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.