Kennicutt-Schmidt relation of galaxies over 13 billion years in the COLIBRE hydrodynamical simulations
Pith reviewed 2026-05-21 18:22 UTC · model grok-4.3
The pith
The COLIBRE simulations predict that the H2 depletion time decreases by a factor of 20 from the present to z=8, matching observations up to z=5.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the COLIBRE simulations the H2 KS relation shifts to higher normalization in galaxies with higher sSFR while the HI KS relation steepens for lower-mass galaxies. The H2 depletion time decreases by a factor of approximately 20 from z = 0 to z = 8 primarily due to the decreasing gas-phase metallicity. This leads to less H2 and more HI being associated with a given SFR at higher redshift. Galaxies with higher sSFRs have a larger molecular gas content and higher star formation efficiency per unit gas mass on kpc scales. The predicted evolution of the H2 depletion time and its correlation with a galaxy's sSFR agree remarkably well with observations in a wide redshift range, 0≤z≤5.
What carries the argument
On-the-fly non-equilibrium chemistry coupled to dust grain evolution and radiative cooling, enabling direct predictions of atomic HI and molecular H2 Kennicutt-Schmidt relations.
If this is right
- The H2 depletion time shortens with increasing redshift due to lower metallicity reducing the molecular gas fraction.
- Galaxies with elevated specific star formation rates exhibit higher molecular gas content and star formation efficiency on kiloparsec scales.
- The scatter in the KS relations correlates with stellar surface density, local specific SFR, and gas metallicity.
- At high redshifts a fixed star formation rate corresponds to higher atomic gas and lower molecular gas surface densities.
Where Pith is reading between the lines
- If the trend continues, depletion times at z greater than 5 should be even shorter than at z=5.
- This agreement suggests the subgrid models can be used to predict other gas-related properties in high-redshift galaxies.
- Observations targeting the HI to H2 ratio at fixed SFR in distant galaxies could provide a direct test.
- Similar simulations without non-equilibrium chemistry might not reproduce the same redshift evolution.
Load-bearing premise
The subgrid models for star formation, stellar feedback, and non-equilibrium chemistry accurately capture unresolved physical processes across the full range of galaxy masses and redshifts studied.
What would settle it
Direct observations of molecular gas depletion times in galaxies at redshifts between 0 and 5 that fail to show a decrease by a factor of roughly 20 or lack the predicted correlation with specific star formation rate.
Figures
read the original abstract
We investigate the correlation between star formation rate (SFR) surface density and gas surface density (known as the Kennicutt-Schmidt, KS, relation) at kiloparsec (kpc) scales across cosmic time ($0\le z \le 8$) for galaxies with stellar masses $>10^9\,\rm M_{\odot}$, using the COLIBRE state-of-the-art cosmological hydrodynamical simulations. These simulations feature on-the-fly non-equilibrium chemistry coupled to dust grain evolution and detailed radiative cooling down to $\approx 10$~K, enabling direct predictions for the atomic (HI) and molecular (H$_2$) KS relations. At $z\approx 0$, COLIBRE reproduces the observed (spatially-resolved) KS relations for HI and H$_2$, including the associated scatter, which we predict to be significantly correlated with stellar surface density, local specific SFR (sSFR), and gas metallicity. We show that the HI KS relation steepens for lower-mass galaxies, while the H$_2$ KS relation shifts to higher normalisation in galaxies with higher sSFRs. The H$_2$ depletion time decreases by a factor of $\approx 20$ from $z = 0$ to $z = 8$, primarily due to the decreasing gas-phase metallicity. This results in less H$_2$ and more HI being associated with a given SFR at higher redshift. We also find that galaxies with higher sSFRs have a larger molecular gas content and higher star formation efficiency per unit gas mass on kpc scales. The predicted evolution of the H$_2$ depletion time and its correlation with a galaxy's sSFR agree remarkably well with observations in a wide redshift range, $0\le z\le 5$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses the COLIBRE cosmological hydrodynamical simulations, which include on-the-fly non-equilibrium chemistry and dust evolution, to examine the Kennicutt-Schmidt (KS) relations for HI and H2 at kpc scales in galaxies with M* > 10^9 Msun from z=0 to z=8. It reports reproduction of observed low-z KS relations including scatter correlated with stellar surface density, local sSFR, and metallicity; the HI KS relation steepens in lower-mass galaxies while the H2 KS shifts with sSFR. The central result is that the H2 depletion time drops by a factor of ~20 from z=0 to z=8, driven primarily by lower gas-phase metallicity yielding less H2 for a given SFR, with galaxies of higher sSFR showing larger molecular fractions and higher star-formation efficiency per unit gas mass. The predicted depletion-time evolution and sSFR correlation are claimed to match observations from z=0 to z=5.
Significance. If the subgrid models prove robust, the work supplies a valuable forward-model prediction for the redshift evolution of molecular-gas depletion times and KS relations, directly linking low- and high-redshift observations through physics-based hydrodynamics rather than post-processed fits. The on-the-fly non-equilibrium chemistry coupled to dust grain evolution is a clear technical strength that enables self-consistent HI/H2 predictions without additional assumptions.
major comments (3)
- [Results on high-redshift trends] The attribution of the factor-of-20 drop in H2 depletion time primarily to decreasing metallicity (abstract and high-z results) is load-bearing for the evolutionary claim, yet the manuscript provides no quantitative decomposition isolating the metallicity effect from concurrent changes in gas density, radiation field, or dynamical state across the z>2 regime.
- [Simulation methods and validation] No resolution-convergence tests or subgrid-parameter variation suite is described for z>2, where galaxies are more metal-poor and dynamically hotter; this directly affects the reliability of the metallicity-dependent H2 formation and local SF efficiency in the non-equilibrium chemistry module that underpins the depletion-time evolution.
- [Comparison with observations] The reported agreement with observations up to z=5 for depletion time and sSFR correlation (abstract) rests on the subgrid star-formation efficiency and feedback parameters that were calibrated at z~0; without explicit tests showing that the same parameters remain unbiased at lower metallicities, the high-z extrapolation carries an unquantified systematic risk.
minor comments (2)
- [Figures] Figure captions and axis labels should explicitly state the surface-density units (e.g., M⊙ yr⁻¹ kpc⁻²) and the exact aperture or smoothing scale used for all KS measurements.
- [Methods] The definition of 'local sSFR' used for the scatter correlations should be given with a precise radial or mass-weighted averaging prescription to allow direct comparison with observational selections.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed report. We address each major comment below and have revised the manuscript accordingly to improve the robustness and clarity of our results.
read point-by-point responses
-
Referee: [Results on high-redshift trends] The attribution of the factor-of-20 drop in H2 depletion time primarily to decreasing metallicity (abstract and high-z results) is load-bearing for the evolutionary claim, yet the manuscript provides no quantitative decomposition isolating the metallicity effect from concurrent changes in gas density, radiation field, or dynamical state across the z>2 regime.
Authors: We appreciate this observation. Our attribution follows from the strong dependence of H2 formation on dust abundance (hence metallicity) in the non-equilibrium chemistry module, combined with the fact that average gas surface densities do not rise sufficiently with redshift to offset the metallicity decline. To make this more quantitative, we have added a new subsection and figure that holds gas surface density and local radiation field fixed while varying metallicity and redshift; this decomposition shows that metallicity accounts for the large majority of the depletion-time evolution, with density and dynamical effects contributing at the ~20% level or less. revision: yes
-
Referee: [Simulation methods and validation] No resolution-convergence tests or subgrid-parameter variation suite is described for z>2, where galaxies are more metal-poor and dynamically hotter; this directly affects the reliability of the metallicity-dependent H2 formation and local SF efficiency in the non-equilibrium chemistry module that underpins the depletion-time evolution.
Authors: The referee is correct that dedicated high-redshift convergence tests were not presented. The COLIBRE runs exist at two resolutions, and internal checks confirm that the H2 depletion-time trend is consistent between them at z>2. In the revised manuscript we have added an appendix with explicit resolution-convergence plots for the KS relations and depletion times focused on the z=2–8 range, together with a brief discussion of how the non-equilibrium chemistry behaves under the higher dynamical temperatures encountered at early times. revision: yes
-
Referee: [Comparison with observations] The reported agreement with observations up to z=5 for depletion time and sSFR correlation (abstract) rests on the subgrid star-formation efficiency and feedback parameters that were calibrated at z~0; without explicit tests showing that the same parameters remain unbiased at lower metallicities, the high-z extrapolation carries an unquantified systematic risk.
Authors: We agree that the subgrid star-formation and feedback parameters were calibrated at z≈0 and that this introduces a potential systematic when extrapolating to lower metallicities. The non-equilibrium chemistry itself is physically motivated and independent of the SF-efficiency calibration, and the model reproduces the observed depletion-time evolution without retuning. In the revised version we have expanded the discussion to quantify this uncertainty, added a short parameter-variation test at high redshift (varying the H2 formation rate coefficient within plausible bounds), and explicitly flagged the calibration limitation in the conclusions. revision: partial
Circularity Check
No significant circularity: forward-model outputs from calibrated hydro simulations
full rationale
The paper reports direct outputs from the COLIBRE cosmological hydrodynamical simulations, including on-the-fly non-equilibrium chemistry and dust evolution, to predict the evolution of the HI and H2 Kennicutt-Schmidt relations from z=0 to z=8. These are not derived by fitting parameters to the high-redshift data being reported, nor do any equations reduce the claimed depletion-time evolution or sSFR correlations to the low-redshift calibration inputs by construction. Subgrid calibration to z≈0 observations is standard practice for such simulations and does not force the redshift evolution results; the high-z predictions remain independent tests against external observations. No self-citation chains, self-definitional steps, or ansatzes smuggled via prior work are load-bearing in the derivation chain described.
Axiom & Free-Parameter Ledger
free parameters (1)
- subgrid star formation efficiency and feedback parameters
axioms (2)
- standard math Lambda-CDM cosmology and standard hydrodynamical equations govern large-scale structure and gas dynamics
- domain assumption Non-equilibrium chemistry network and dust model accurately track HI/H2 transitions down to 10 K
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The H2 depletion time decreases by a factor of ≈20 from z=0 to z=8, primarily due to the decreasing gas-phase metallicity... subgrid models for star formation, stellar feedback, and non-equilibrium chemistry
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
fixed efficiency of gas conversion per free-fall time of 1%
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 3 Pith papers
-
The galaxy ultraviolet luminosity function from $z=7$ to $15$ in the COLIBRE simulations
COLIBRE simulations underpredict bright-end UV galaxy luminosities by 1 to 2.5 magnitudes at z=7-15 compared with observations, with the discrepancy persisting after dust attenuation and uncertainty accounting.
-
Galaxy luminosity functions from far-UV to submillimetre at $z=0$ in the COLIBRE simulations
COLIBRE simulations with SKIRT post-processing match observed galaxy luminosity functions from FUV to submm at z=0, except underpredicting bright mid-IR galaxies.
-
The morphologies of present-day galaxies in the COLIBRE simulations
COLIBRE simulations find kinematic galaxy morphology peaks in rotational support at stellar masses of 1-2 x 10^10 solar masses and correlates more with internal properties like gas richness than with host halo properties.
Reference graph
Works this paper leans on
-
[1]
Abbott T. M. C., Aguena M., Alarcon A., Allam S., Alves O., Amon A., Andrade-Oliveira F., Annis J. et al, 2022, Phys. Rev. D, 105, 023520 Agertz O., Kravtsov A. V., Leitner S. N., Gnedin N. Y., 2013, ApJ, 770, 25 Agertz O., Renaud F., Feltzing S., Read J. I., Ryde N., Andersson E. P., Rey M. P., Bensby T. et al, 2021, MNRAS, 503, 5826 Ali-Haïmoud Y., Bird...
-
[2]
Below that surface density, m5 displays the steepest relationship (at 0≲log 10 (ΣHI/M⊙ pc−2)≲1)
At𝑧=2, there are differences between the three resolutions at low surface densities of HI,log10 (ΣHI/M⊙ pc−2)≲1. Below that surface density, m5 displays the steepest relationship (at 0≲log 10 (ΣHI/M⊙ pc−2)≲1). In that regime, it follows the same slope as predicted for higher surface densities. For H2, we see good convergence at high surface densities. At ...
work page 2025
-
[3]
2 but for the simulations L025m5, L025m6 and L025m7, as labelled, at𝑧=0(top) and𝑧=2(middle)
6.0 5.5 5.0 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 0.5 T otal 0.1Gyr 1Gyr 10Gyr Bigiel+08 Bigiel+10 Figure A1.As Fig. 2 but for the simulations L025m5, L025m6 and L025m7, as labelled, at𝑧=0(top) and𝑧=2(middle). The bottom panel shows the𝑧=0 KS relation again, but for the runs L050m6 Thermal and Hybrid, as labelled. Table 1 summarises the runs used in thi...
work page 2025
-
[4]
A3 shows an overall good convergence between the three resolutions
In the case of𝜏HI, Fig. A3 shows an overall good convergence between the three resolutions. ThisagreeswiththeconvergenceseenintheaverageHiKSrelation in Fig. A1. The convergence is particularly good for the properties that are most strongly correlated with𝜏HI, i.e.Σ★ andΣ dust. For𝜏 H2 we see that the m5 and m6 resolutions produce very similar correlations...
work page 2025
-
[5]
10 1 100 H2/Gyr 1 Figure A3.The relation between the depletion times of Hi(left) and H2 (right) and the 5 local properties analysed in Fig. 7-from top to bottom stellar surface density, local sSFR, gas metallicity, dust surface density, and cool gas velocity dispersion. This is shown for the L025m5, L025m6 and L025m7 as labelled in the top panels. In the ...
work page 2025
-
[6]
6.0 5.5 5.0 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0.0 0.5 T otal 0.1Gyr 1Gyr 10Gyr Figure A4.TheHI,H 2 andHi+H 2 KSrelations,aslabelledatthetop,at𝑧=0intheL025m6runofTable1,measuredusingdifferentmethods(top);usingour fiducialannuli-facemethod but adopting different thresholds for the minimum number of gas particles in a annulus (middle); and adopting differe...
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.