Recognition: unknown
Cosmic evolution of the [CII]-to-molecular gas relation
Pith reviewed 2026-05-10 02:16 UTC · model grok-4.3
The pith
No single conversion factor from [CII] luminosity to molecular gas mass works across all redshifts.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The global L_[CII]-M_mol relation evolves from a steep, [CII]-deficient regime at very low metallicity to an almost linear behaviour, similar to calibrations at z approximately 2, once the ISM reaches Z greater than or equal to 0.05-0.1 solar at z less than or equal to 5. Over this evolution, alpha_[CII] spans nearly three orders of magnitude, from greater than 10,000 down to approximately 10 solar masses per solar luminosity, even though the [CII] emission remains spatially correlated with the molecular gas. A unique, redshift-independent alpha_[CII] therefore cannot recover molecular gas masses across the regimes explored. [CII] remains a viable tracer of molecular gas at very high redshif
What carries the argument
The redshift-dependent [CII]-to-molecular gas conversion factor alpha_[CII], measured globally and at 100 parsec scales from synthetic data cubes, which varies with the evolving metallicity, ISM phase mix, and merger activity in the simulation.
If this is right
- At very high redshifts, [CII] traces molecular gas only when conversion factors explicitly account for metallicity, ISM phase mix, and merger events.
- The L_[CII]-M_mol and L_[CII]-SFR relations approach linear behaviour once metallicity exceeds roughly 0.05 solar, resembling z approximately 2 calibrations.
- Spatially resolved (100 pc) [CII] emission stays correlated with molecular gas even while the global conversion factor varies by orders of magnitude.
- Merger-driven changes in ISM phase structure introduce additional offsets in alpha_[CII] beyond the metallicity trend.
Where Pith is reading between the lines
- High-redshift surveys that apply local alpha_[CII] values may systematically underestimate molecular gas masses in metal-poor systems.
- The metallicity dependence found here could be tested directly by pairing [CII] maps with metallicity diagnostics in future ALMA or JWST observations.
- Similar conversion-factor evolution may apply to other tracers such as CO lines, affecting gas-mass estimates in the same high-z regime.
- Galaxy-formation models that assume a fixed [CII] conversion would mispredict star-formation efficiencies at z greater than 5.
Load-bearing premise
The Vintergatan zoom-in simulation plus Skirt post-processing faithfully reproduces the ISM phase structure, metallicity distribution, and [CII] excitation conditions in real high-redshift galaxies.
What would settle it
Independent measurements of molecular gas mass (for example via CO or dust) and [CII] luminosity in a statistical sample of z greater than 5 galaxies that yield a constant alpha_[CII] near the local value would falsify the predicted strong redshift evolution.
Figures
read the original abstract
The [CII] 158 $\mu$m line is widely used to trace star formation and the gas contents of high-redshift galaxies. However, it remains unclear under which physical conditions it reliably traces the molecular reservoir, and whether a unique conversion factor $\alpha_{\rm [CII]}$ can be applied across cosmic time. We investigate the evolution of the relation between the [CII] luminosity and molecular gas mass from $z\simeq10$ to $z\simeq0.2$ using the Vintergatan simulation, a high-resolution cosmological zoom-in of a Milky Way-like galaxy. We post-process the snapshots with the Skirt radiative transfer code to generate synthetic [CII] data cubes. We measure global and spatially resolved (100 pc) relations between [CII] luminosity ($L_{\rm [CII]}$), star formation rate (SFR), and molecular gas mass ($M_{\rm mol}$). We follow the redshift evolution of the [CII]-to-molecular gas conversion factor $\alpha_{\rm [CII]}$, and link these trends to the evolution of the interstellar medium (ISM) phases. The global $L_{\rm [CII]}$-$M_{\rm mol}$ and $L_{\rm [CII]}$-SFR relations evolve from a steep, [CII]-deficient regime at very low metallicity to an almost linear behaviour, similar to calibrations at $z\approx2$, once the ISM reaches $Z \gtrsim 0.05$-$0.1\,Z_\odot$ at $z\lesssim5$. Over this evolution, $\alpha_{\rm [CII]}$ spans nearly three orders of magnitude, from $\gtrsim 10^4$ down to $\approx10 \,\rm{M_\odot\,L_\odot^{-1}}$, even though the [CII] emission remains spatially correlated with the molecular gas. A unique, redshift-independent $\alpha_{\rm [CII]}$ therefore cannot recover molecular gas masses across the regimes we explore. [CII] remains a viable tracer of molecular gas at very high redshifts, but only when used with conversion factors that explicitly account for metallicity, ISM phase mix, and merger events.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses the Vintergatan cosmological zoom-in simulation of a Milky Way-like halo, post-processed with the SKIRT radiative transfer code, to track the redshift evolution (z≈10 to z≈0.2) of the [CII] 158 μm luminosity–molecular gas mass relation. Synthetic [CII] data cubes are generated to measure global and 100-pc resolved L_[CII]–M_mol and L_[CII]–SFR relations, from which the conversion factor α_[CII] is derived and linked to ISM phase evolution, metallicity, and merger activity. The key result is that α_[CII] varies by nearly three orders of magnitude, precluding a unique redshift-independent value, while [CII] remains a viable tracer only when metallicity-, phase-, and merger-dependent corrections are applied.
Significance. If the reported trends hold, the work supplies a physically motivated framework for interpreting high-redshift [CII] observations, demonstrating that forward-modeling from hydrodynamics plus line transfer can reveal when and why [CII] fails as a direct molecular-gas tracer. The absence of circularity—trends emerge from the simulation rather than being tuned to observations—is a clear strength, as is the explicit connection between α_[CII] evolution and ISM phase partitioning. The result would affect ALMA-based gas-mass estimates at z≳5 and motivate metallicity-dependent conversion factors in future surveys.
major comments (2)
- [abstract and results] The central claim that no unique, redshift-independent α_[CII] recovers M_mol across regimes rests on trends extracted from a single Vintergatan zoom-in of one Milky Way-like halo. Because the reported three-order-of-magnitude variation in α_[CII] is driven by this halo’s specific metallicity build-up, merger sequence, and ISM phase partitioning, it is unclear whether the same evolution applies to galaxies with different stellar masses, assembly histories, or feedback implementations at z≳5 (abstract and results sections).
- [methods and discussion] The assumption that the Vintergatan + SKIRT pipeline faithfully reproduces high-redshift ISM conditions is load-bearing for the conclusion that [CII] remains viable only with explicit corrections. No quantitative assessment is provided of how changes in sub-grid physics, spatial resolution, or [CII] excitation assumptions would alter the derived α_[CII](z) trends (methods and discussion sections).
minor comments (2)
- [results] Clarify the exact spatial scale and aperture used for the 'global' versus '100 pc resolved' relations; the transition between them is not numerically specified.
- [results] Add a brief table or figure panel showing the redshift evolution of mean metallicity and molecular gas fraction alongside α_[CII] to make the causal links more transparent.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments, which have prompted us to clarify several aspects of our analysis. We respond to each major comment below and indicate the revisions we will make to the manuscript.
read point-by-point responses
-
Referee: [abstract and results] The central claim that no unique, redshift-independent α_[CII] recovers M_mol across regimes rests on trends extracted from a single Vintergatan zoom-in of one Milky Way-like halo. Because the reported three-order-of-magnitude variation in α_[CII] is driven by this halo’s specific metallicity build-up, merger sequence, and ISM phase partitioning, it is unclear whether the same evolution applies to galaxies with different stellar masses, assembly histories, or feedback implementations at z≳5 (abstract and results sections).
Authors: We agree that the study is based on a single high-resolution zoom-in simulation of a Milky Way-like halo. This setup was selected to enable a self-consistent, high-fidelity tracking of one galaxy’s ISM evolution from z≈10 to z≈0.2. The reported variation in α_[CII] is physically tied to the secular increase in gas-phase metallicity and the associated shift in ISM phase partitioning, both of which are generic processes expected to operate across galaxies of different masses. The precise redshift at which the transition to a near-linear L_[CII]–M_mol relation occurs will, however, depend on a galaxy’s specific assembly history and feedback implementation. In the revised manuscript we will add an explicit paragraph in the discussion section that states this limitation, cautions against over-generalizing the quantitative α_[CII](z) curve, and notes that the qualitative conclusion (no unique redshift-independent conversion factor) follows from the metallicity dependence itself, which is a robust feature of the underlying physics. revision: partial
-
Referee: [methods and discussion] The assumption that the Vintergatan + SKIRT pipeline faithfully reproduces high-redshift ISM conditions is load-bearing for the conclusion that [CII] remains viable only with explicit corrections. No quantitative assessment is provided of how changes in sub-grid physics, spatial resolution, or [CII] excitation assumptions would alter the derived α_[CII](z) trends (methods and discussion sections).
Authors: We recognize that a quantitative sensitivity study varying sub-grid physics, resolution, or the [CII] excitation treatment is not presented. Performing such a study would require a new suite of simulations and is beyond the computational scope of the present work. The Vintergatan model has been validated against a range of observational diagnostics at both low and high redshift, and SKIRT is a standard, well-tested radiative-transfer code. In the revised manuscript we will expand the methods and discussion sections with a qualitative assessment of how plausible changes in these ingredients could affect the results, drawing on existing literature that explores [CII] emission under different ISM conditions and resolutions. We will emphasize that the dominant driver of the α_[CII] evolution is the large-scale metallicity build-up rather than small-scale details of the sub-grid modeling. revision: partial
Circularity Check
No circularity: forward-modeling yields emergent α_[CII] trends
full rationale
The paper runs the Vintergatan zoom-in simulation, post-processes snapshots with Skirt to produce synthetic [CII] cubes, then directly computes global and resolved L_[CII]–M_mol relations and the resulting α_[CII](z) from those outputs. No observational data are used to fit α_[CII], no self-referential definitions equate inputs to outputs, and no uniqueness theorems or prior self-citations are invoked as load-bearing premises for the central claim. The reported three-order-of-magnitude variation and the conclusion that a redshift-independent factor fails both follow from the hydrodynamics, metallicity build-up, and radiative transfer inside the model rather than being presupposed by construction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The Vintergatan simulation and its sub-grid physics accurately capture the metallicity evolution and ISM phase structure of high-redshift galaxies.
- domain assumption The Skirt radiative transfer code correctly computes the [CII] 158 micron emission given the simulated gas and radiation field.
Reference graph
Works this paper leans on
-
[1]
2025, A&A, 702, A206
Accard, C., Béthermin, M., Boquien, M., et al. 2025, A&A, 702, A206
2025
-
[2]
2017, MNRAS, 470, 4750
Accurso, G., Saintonge, A., Catinella, B., et al. 2017, MNRAS, 470, 4750
2017
-
[3]
2021, MNRAS, 503, 5826
Agertz, O., Renaud, F., Feltzing, S., et al. 2021, MNRAS, 503, 5826
2021
-
[4]
Algera, H. S. B., Rowland, L., Stefanon, M., et al. 2026, MNRAS, 545, staf1897 Astropy Collaboration, Price-Whelan, A. M., Sip˝ocz, B. M., et al. 2018, AJ, 156, 123
2026
-
[5]
2011, ApJS, 196, 22 Béthermin, M., Fudamoto, Y ., Ginolfi, M., et al
Baes, M., Verstappen, J., De Looze, I., et al. 2011, ApJS, 196, 22 Béthermin, M., Fudamoto, Y ., Ginolfi, M., et al. 2020, A&A, 643, A2
2011
-
[6]
D., Wolfire, M., & Leroy, A
Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207 Bouché, N., Dekel, A., Genzel, R., et al. 2010, ApJ, 718, 1001
2013
-
[7]
J., Smit, R., Schouws, S., et al
Bouwens, R. J., Smit, R., Schouws, S., et al. 2022, ApJ, 931, 160
2022
-
[8]
& Charlot, S
Bruzual, G. & Charlot, S. 2003, MNRAS, 344, 1000
2003
-
[9]
& Baes, M
Camps, P. & Baes, M. 2015, Astronomy and Computing, 9, 20
2015
-
[10]
& Baes, M
Camps, P. & Baes, M. 2020, Astronomy and Computing, 31, 100381
2020
-
[11]
Carilli, C. L. & Walter, F. 2013, ARA&A, 51, 105
2013
-
[12]
2003, PASP, 115, 763
Chabrier, G. 2003, PASP, 115, 763
2003
-
[13]
J., et al
Choustikov, N., Katz, H., Cameron, A. J., et al. 2026, The Open Journal of As- trophysics, 9, 58199 da Cunha, E., Groves, B., Walter, F., et al. 2013, ApJ, 766, 13 De Looze, I., Cormier, D., Lebouteiller, V ., et al. 2014, A&A, 568, A62
2026
-
[14]
Draine, B. T. & Salpeter, E. E. 1979, ApJ, 231, 77
1979
-
[15]
L., Schaerer, D., Lemaux, B
Faisst, A. L., Schaerer, D., Lemaux, B. C., et al. 2020, ApJS, 247, 61
2020
-
[16]
2019, MNRAS, 489, 1
Ferrara, A., Vallini, L., Pallottini, A., et al. 2019, MNRAS, 489, 1
2019
-
[17]
W., Lang, D., & Goodman, J
Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306
2013
-
[18]
J., et al
Freundlich, J., Combes, F., Tacconi, L. J., et al. 2013, A&A, 553, A130
2013
-
[19]
Fudamoto, Y ., Smit, R., Bowler, R. A. A., et al. 2022, ApJ, 934, 144
2022
-
[20]
J., Gracia-Carpio, J., et al
Genzel, R., Tacconi, L. J., Gracia-Carpio, J., et al. 2010, MNRAS, 407, 2091
2010
-
[21]
Glover, S. C. O., Federrath, C., Mac Low, M.-M., & Klessen, R. S. 2010, MN- RAS, 404, 2
2010
-
[22]
K., et al
Harikane, Y ., Ouchi, M., Inoue, A. K., et al. 2020, ApJ, 896, 93
2020
-
[23]
2025, A&A, 699, A80
Herrera-Camus, R., González-López, J., Förster Schreiber, N., et al. 2025, A&A, 699, A80
2025
-
[24]
F., Kereš, D., Oñorbe, J., et al
Hopkins, P. F., Kereš, D., Oñorbe, J., et al. 2014, MNRAS, 445, 581
2014
-
[25]
F., Quataert, E., & Murray, N
Hopkins, P. F., Quataert, E., & Murray, N. 2012, MNRAS, 421, 3488
2012
-
[26]
Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90
2007
-
[27]
Evaluating star formation rates at z = 5
Ismail, D., Kraljic, K., Béthermin, M., et al. 2026, arXiv e-prints, arXiv:2601.05916
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[28]
P., Köhler, M., Ysard, N., Bocchio, M., & Verstraete, L
Jones, A. P., Köhler, M., Ysard, N., Bocchio, M., & Verstraete, L. 2017, A&A, 602, A46
2017
-
[29]
U., Baes, M., van der Wel, A., et al
Kapoor, A. U., Baes, M., van der Wel, A., et al. 2023, MNRAS, 526, 3871
2023
-
[30]
U., Baes, M., van der Wel, A., et al
Kapoor, A. U., Baes, M., van der Wel, A., et al. 2024, A&A, 692, A79
2024
-
[31]
Katz, H., Rey, M. P., Cadiou, C., et al. 2025, arXiv e-prints, arXiv:2510.05201
-
[32]
R., Burkhart, B., Forbes, J
Krumholz, M. R., Burkhart, B., Forbes, J. C., & Crocker, R. M. 2018, MNRAS, 477, 2716
2018
-
[33]
R., Dekel, A., & McKee, C
Krumholz, M. R., Dekel, A., & McKee, C. F. 2012, ApJ, 745, 69
2012
-
[34]
R., McKee, C
Krumholz, M. R., McKee, C. F., & Tumlinson, J. 2009, ApJ, 693, 216
2009
-
[35]
2018, A&A, 609, A130 Le Fèvre, O., Béthermin, M., Faisst, A., et al
Lagache, G., Cousin, M., & Chatzikos, M. 2018, A&A, 609, A130 Le Fèvre, O., Béthermin, M., Faisst, A., et al. 2020, A&A, 643, A1
2018
-
[36]
L., Faucher-Giguère, C.-A., & Lidz, A
Liu, L.-J., Sun, G., Faisst, A. L., Faucher-Giguère, C.-A., & Lidz, A. 2026, arXiv e-prints, arXiv:2601.00959
-
[37]
L., Satyapal, S., Fischer, J., et al
Luhman, M. L., Satyapal, S., Fischer, J., et al. 2003, ApJ, 594, 758
2003
-
[38]
2020, MNRAS, 496, 5160
Lupi, A., Pallottini, A., Ferrara, A., et al. 2020, MNRAS, 496, 5160
2020
-
[39]
F., Faucher-Giguère, C.-A., et al
Ma, X., Hopkins, P. F., Faucher-Giguère, C.-A., et al. 2016, MNRAS, 456, 2140
2016
-
[40]
& Dickinson, M
Madau, P. & Dickinson, M. 2014, ARA&A, 52, 415
2014
-
[41]
C., Cormier, D., Hony, S., et al
Madden, S. C., Cormier, D., Hony, S., et al. 2020, A&A, 643, A141
2020
-
[42]
J., Hollenbach, D., et al
Malhotra, S., Kaufman, M. J., Hollenbach, D., et al. 2001, ApJ, 561, 766
2001
-
[43]
& Krumholz, M
Narayanan, D. & Krumholz, M. R. 2017, MNRAS, 467, 50
2017
-
[44]
P., Burkhart, B., Mac Low, M.-M., et al
Olsen, K. P., Burkhart, B., Mac Low, M.-M., et al. 2021, The Astrophysical Jour- nal, 922, 88
2021
-
[45]
Pabst, C. H. M., Goicoechea, J. R., Hacar, A., et al. 2022, A&A, 658, A98
2022
-
[46]
S., Werk, J
Peeples, M. S., Werk, J. K., Tumlinson, J., et al. 2014, ApJ, 786, 54
2014
-
[47]
L., Langer, W
Pineda, J. L., Langer, W. D., & Goldsmith, P. F. 2014, A&A, 570, A121
2014
-
[48]
2019, A&A, 621, A104
Renaud, F., Bournaud, F., Daddi, E., & Weiß, A. 2019, A&A, 621, A104
2019
-
[49]
2025, A&A, 694, A56
Renaud, F., Ratcliffe, B., Minchev, I., et al. 2025, A&A, 694, A56
2025
-
[50]
2024, A&A, 689, A273
Rizzo, F., Bacchini, C., Kohandel, M., et al. 2024, A&A, 689, A273
2024
-
[51]
2024, A&A, 687, A35
Roman-Oliveira, F., Rizzo, F., & Fraternali, F. 2024, A&A, 687, A35
2024
-
[52]
E., Hodge, J., Bouwens, R., et al
Rowland, L. E., Hodge, J., Bouwens, R., et al. 2024, MNRAS, 535, 2068
2024
-
[53]
2020, A&A, 643, A3 Segovia Otero, Á., Renaud, F., & Agertz, O
Schaerer, D., Ginolfi, M., Béthermin, M., et al. 2020, A&A, 643, A3 Segovia Otero, Á., Renaud, F., & Agertz, O. 2022, MNRAS, 516, 2272
2020
-
[54]
Somerville, R. S. & Davé, R. 2015, ARA&A, 53, 51
2015
-
[55]
J., Genzel, R., Saintonge, A., et al
Tacconi, L. J., Genzel, R., Saintonge, A., et al. 2018, ApJ, 853, 179
2018
-
[56]
J., Neri, R., Genzel, R., et al
Tacconi, L. J., Neri, R., Genzel, R., et al. 2013, ApJ, 768, 74
2013
-
[57]
2002, A&A, 385, 337
Teyssier, R. 2002, A&A, 385, 337
2002
-
[58]
2015, ApJ, 813, 36
Vallini, L., Gallerani, S., Ferrara, A., Pallottini, A., & Yue, B. 2015, ApJ, 813, 36
2015
-
[59]
2025, A&A, 700, A117
Vallini, L., Pallottini, A., Kohandel, M., et al. 2025, A&A, 700, A117
2025
-
[60]
E., et al
Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261
2020
-
[61]
R., Olsen, K
Vizgan, D., Greve, T. R., Olsen, K. P., et al. 2022, ApJ, 929, 92
2022
-
[62]
2020, ApJ, 902, 111
Walter, F., Carilli, C., Neeleman, M., et al. 2020, ApJ, 902, 111
2020
-
[63]
G., Vallini, L., & Chevance, M
Wolfire, M. G., Vallini, L., & Chevance, M. 2022, ARA&A, 60, 247
2022
-
[64]
2018, MNRAS, 481, 1976 Article number, page 12 C
Zanella, A., Daddi, E., Magdis, G., et al. 2018, MNRAS, 481, 1976 Article number, page 12 C. Accard et al.: Cosmic evolution of the [Cii]-to-molecular gas relation Appendix A: Orientation, extraction radius, timescale, and transparency effects In this appendix we quantify how viewing angle, extraction aper- ture, and line opacity affect the inferred [Cii]...
2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.