Pith. sign in

REVIEW 2 major objections 4 minor 2 cited by

Predicting the resolved CO emission of $z=1-3$ star-forming galaxies

T0 review · 2 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read This paper predicts a steep radial rise in the CO-to-H2 conversion factor in z=1–3 galaxies, from ~1–5 in the centers to >100 at ~15 kpc, so CO emission cannot be read as a uniform tracer of molecular gas across these disks.

desk verdict Solid statistical advance in predicting resolved CO at z=1-3; the central gradient is likely robust, but the magnitude hinges on an unscaled CR prescription the authors themselves flag. read the letter →

arxiv 2506.13899 v2 pith:TMQRGUCX submitted 2025-06-16 astro-ph.GA

classification astro-ph.GA
keywords COemissionmoleculargasH2massCO-to-H2conversionfactorhigh-redshiftgalaxiescosmicnoonsyntheticobservationsgalaxysizes
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper asks how reliably resolved CO emission traces molecular hydrogen in the peak star-forming epoch, z=1–3, by generating synthetic CO(1–0) through CO(5–4) maps for hundreds of main-sequence galaxies drawn from a cosmological simulation. It argues that the CO-to-H2 conversion factor $\alpha_{\rm CO}$ is not constant across these disks: it stays near 1–5 $\mathrm{M_\odot}\,(\mathrm{K\,km\,s^{-1}\,pc^2})^{-1}$ in the central 1–3 kpc, then climbs steeply to values above 100 at roughly 15 kpc. The consequence is that using a single $\alpha_{\rm CO}$, or reading a CO half-light radius as the H2 size, systematically underweights the extended molecular gas: more than a quarter of H2 lies beyond twice the CO(1–0) half-light radius, and CO half-light radii are on average about 27–29% smaller than H2 half-mass radii. If correct, this matters because most resolved CO studies at these redshifts are effectively seeing only the dense centers of galaxies, not their full molecular reservoirs.

What carries the argument

The carrying mechanism is the SLICK modeling pipeline: each gas particle in the cosmological simulation is treated as a virialized, pressure-confined cloud with a power-law density profile, split into concentric zones; a chemical network solves for H2 and CO abundances under an equilibrium between heating, from photoelectric, dust, cosmic-ray, and CMB processes, and cooling, and a spectral-line code computes CO line luminosities using escape-probability radiative transfer, with the interstellar radiation field and cosmic-ray field set locally from the 64 nearest neighboring clouds. This sub-resolution machinery is what converts the simulation's unresolved gas particles into resolved, roughly 500 pc CO maps, and the paper's radial $\alpha_{\rm CO}$ gradients, compact CO sizes, and excitation trends all follow from how these cloud properties vary with galactocentric radius.

What would settle it

Measure resolved CO(1–0) and CO(3–2) emission together with an independent, dust-based molecular-gas surface density map for a sample of z~2 main-sequence galaxies. If $\alpha_{\rm CO}$ stays near 3–5 $\mathrm{M_\odot}\,(\mathrm{K\,km\,s^{-1}\,pc^2})^{-1}$ beyond 10 kpc, or if the CO(1–0) half-light radius matches the H2 half-mass radius instead of sitting about 28% smaller, the predicted steep gradient and size bias are ruled out.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that resolved CO emission in z=1–3 star-forming galaxies is not a uniform tracer of H2. Using sub-resolution cloud modeling plus an equilibrium line-radiation calculation on galaxies from a cosmological simulation, the authors predict that $\alpha_{\rm CO(1-0)}$ rises by one to two orders of magnitude from galaxy centers to about 15 kpc, that higher-$J$ CO transitions trace even more centrally concentrated gas, and that CO excitation, as measured by ratios like CO(3–2)/CO(1–0), declines from suprathermal values near 0.9–2 in the centers to roughly 0.1 in the outskirts, with line ratios increasing toward higher redshift and lower stellar mass. They attribute the $\alpha_{\rm CO}$ gradient primarily to a declining CO abundance per H2 molecule and to falling density, radiation field, and turbulence at large radius, rather than to metallicity alone. The predicted CO half-light radii are compact, around 1–5 kpc, comparable to current observations, but about 27–29% smaller than the H2 half-mass radii, with about 30% of H2 located outside twice the CO(1–0) half-light radius.

Load-bearing premise

Everything rests on treating each unresolved gas particle in the simulation as a single virialized, pressure-confined molecular cloud with equilibrium chemistry and a cosmic-ray field that scales linearly with local star-formation surface density; if z=1–3 gas is not organized this way, the predicted radial $\alpha_{\rm CO}$ gradient, compact CO sizes, and excitation trends are model artifacts.

Editorial extensions

If this is right

  • Observed CO half-light radii at z=1–3 should be interpreted as sizes of the dense, star-forming molecular cores, not the full H2 disk; on average they are about 27–29% smaller than H2 half-mass radii.
  • Resolved studies that assume a constant $\alpha_{\rm CO}$ will systematically underestimate molecular gas surface densities in galaxy outskirts, by factors of tens to hundreds beyond roughly 10 kpc.
  • Tracing H2 beyond 3–5 kpc in cosmic-noon galaxies is predicted to be very difficult with current facilities because the CO surface brightness drops faster than the H2 density.
  • Higher-$J$ transitions such as CO(5–4) are progressively more compact and excitation-biased toward dense central gas, so different CO lines cannot be used interchangeably for individual galaxies.
  • CO excitation gradients, with line ratios peaking in centers and increasing with redshift, provide a way to infer gas physical state: suprathermal ratios require warm, roughly 30–100 K, and dense, above 100 cm$^{-3}$, gas.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A direct observational test could come from combining resolved CO(1–0) with dust-based or other independent H2 maps in the same z~2 galaxies; if those maps confirm that CO half-light radii understate H2 sizes by about 30%, then integrated CO-based gas masses at high redshift may also be biased, not just resolved profiles.
  • Because the predicted $\alpha_{\rm CO}$ gradient is driven more by the local interstellar radiation field and density than by metallicity, galaxy-wide $\alpha_{\rm CO}$–metallicity scaling relations may hide large radial scatter, and resolved local-galaxy calibrations are the nearest analogues that could test the sub-grid assumptions.
  • A testable extension is to rerun the same pipeline with alternative cosmic-ray prescriptions, for instance with a flat rather than star-formation-proportional cosmic-ray field; the paper flags this as a key uncertainty, and since cosmic-ray heating boosts CO(1–0) brightness and high-$J$ excitation, the predicted $\alpha_{\rm CO}$ gradient and suprathermal ratios would weaken or strengthen accordin
  • The result implies that future high-resolution CO surveys of cosmic noon should either include low-$J$ CO(1–0) to capture the extended diffuse gas or explicitly model a radially varying $\alpha_{\rm CO}$; otherwise stacked CO sizes will be biased toward compact centers even when individual galaxies show a wide scatter.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. This paper uses the Simba cosmological simulation combined with the SLICK post-processing pipeline (Garcia et al. 2024) to generate synthetic, spatially resolved CO(1–0) to CO(5–4) maps for star-forming galaxies at z = 0.8–3.5, and compares the CO emission with the underlying H2 mass distribution. The central claims are that the CO-to-H2 conversion factor αCO rises from about 1–5 M⊙ (K km s−1 pc2)−1 in the central 1–3 kpc of massive galaxies to values exceeding 100 at ~15 kpc, that CO half-light radii are on average 27–29% smaller than H2 half-mass radii with roughly 30% of the H2 mass lying beyond two CO half-light radii, and that CO excitation ratios increase toward galaxy centers and with redshift, reaching suprathermal values in warm, dense gas. The paper argues that these results imply that a single αCO cannot be used to trace H2 across z = 1–3 disks and that combined CO plus H2 corrections will be needed for resolved observations.

Significance. If these predictions hold, they are directly relevant to the interpretation of current and upcoming resolved CO observations at z = 1–3: the paper provides a physically motivated, spatially resolved αCO prescription, predicts quantitative size-ratio corrections, and makes falsifiable predictions for CO line-ratio gradients. The work is carefully presented, with percentile spreads, a localized ISRF treatment that improves agreement with z = 0 GMC data, clear comparisons to observed sizes and line ratios, and an explicit list of caveats. The target z = 1–3 outputs are forward predictions rather than fits to the same data, so the exercise is not circular. The main weakness is that the headline quantitative claims depend strongly on the assumed cosmic-ray ionization rate scaling, which is both ambiguously written in Eq. (9) and acknowledged in Sec. 4.2 to be potentially overestimated in galaxy centers; no robustness test is provided. This needs to be addressed before the magnitude of the αCO gradient and the size ratios can be fully endorsed.

major comments (2)
  1. [Sec. 4.2; Eq. (9); Fig. 4; Table 1] The central quantitative claims—the αCO gradient in Fig. 4 and the CO half-light to H2 half-mass radius ratios in Table 1—are directly affected by the assumed cosmic-ray ionization rate. Sec. 4.1 states that the higher gas temperature caused by increased ISRF and CR strength raises the CO(1–0) brightness and thus lowers αCO in the centers. The authors acknowledge in Sec. 4.2 that the linear Σ_SFR–ξ_CR scaling 'may not hold in all environments', and, citing Krumholz et al. (2023), that starburst systems may have only moderately enhanced CR rates, which would imply that central CR heating is overestimated in the present model. Yet no alternative CR prescription or sensitivity test is presented. I request a robustness test with a plausible alternative CR treatment (e.g., a capped CR enhancement or the Krumholz et al. prescription) recomputing at least the αCO(r) profiles, the r31(r) profiles, and the Table 1 size ratios. Without such a test, the magnitude of the predicted αCO range and the 27–29% size ratios are not demonstrated to be robust to an uncertainty the authors themselves flag.
  2. [Sec. 2.2.3; Eq. (9); Sec. 4.2] The description of the CR normalization is internally inconsistent and invites a factor-of-10 misreading. The text says 'assuming ζ−16 = 0.1 is the CR ionization rate and ξ_CR,MW = 10−16 s−1 is the CR field in the solar neighborhood'. If ζ−16 is a dimensionless prefactor, calling it the ionization rate is wrong; if ζ−16 is already a rate, multiplying it by ξ_CR,MW double-counts the normalization. The product of the two values is 10−17 s−1, consistent with the '10^-17 s^-1' stated in Sec. 4.2, but a reader implementing Eq. (9) literally from Sec. 2.2.3 could easily adopt 10−16 s−1. Please rewrite the sentence and Eq. (9) so that ζ−16 is unambiguously a dimensionless scaling relative to ξ_CR,MW, or state the solar-neighborhood CR ionization rate directly and use only one symbol for it.
minor comments (4)
  1. [Sec. 2.3; Figs. 3–7; Table 1] The number of galaxies entering each redshift–stellar mass bin after the merger and convergence cuts is not reported anywhere. Because some bins are strongly affected by the merger cut (the text notes only one high-mass z = 2.5–3.5 galaxy and omits that bin), the 16th–84th percentile spreads would be much more interpretable if N per bin were given in each panel or in Table 1.
  2. [Sec. 2.2.2; Sec. 3.1] The definition of the molecular gas mass used for αCO includes the helium contribution (0.4×M_C) and is used interchangeably with M_H2. Since αCO is then effectively a molecular-gas-to-CO ratio rather than a pure H2-to-CO ratio, the comparison to literature αCO values (e.g., the observed values discussed in Sec. 4.1) should state explicitly whether the observational values include helium (i.e., the usual 1.36 correction).
  3. [Fig. 8] The observational data points ('diamonds') in Fig. 8 are not individually identified or assigned uncertainties. Please add a legend or a caption note specifying which source (Kaasinen, Ikeda, Rizzo, Tadaki, and where applicable individual galaxies) corresponds to each symbol, so that the comparison is reproducible.
  4. [Sec. 2.3] The prose in Sec. 2.3 says the three bins are 'around z=1,2,3', while the figures use z=0.8–1.2, z=1.8–2.2, and z=2.5–3.5. This is clear in the figures, but the text should use the wider labels consistently to avoid confusion between the snapshot centers and the bin edges.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the predicted radial alpha_CO, excitation, and size trends are forward-model outputs from a pipeline validated on external z=0 data, not fits to the quantities being predicted.

full rationale

The predicted radial alpha_CO profiles, CO excitation gradients, and half-light-to-half-mass size ratios are not fitted to the quantities they are used to claim. SLICK's sub-resolution cloud properties (Eqs. 2-4), the chemical network, the ISRF and CR prescriptions (Eqs. 8-9), and Despotic's line solver are fixed inputs adopted from prior work; alpha_CO = M_H2 / L'_CO is computed as an output of that forward model, not imposed as a target. The only calibration against data is the z=0 M_H2-L'_CO comparison from Garcia et al. (2024), which is an external benchmark outside this paper's target (resolved z=1-3 CO emission); using it as a validation anchor does not make the z=1-3 predictions fitted. The CR-Sigma_SFR scaling in Eq. (9) is a clearly stated modeling assumption rather than a hidden restatement of the results, and Sec. 4.2 explicitly concedes that overestimating the CR ionization rate could overpredict the line ratios and that the assumed scaling 'may not hold in all environments'; this is an acknowledged uncertainty to be explored with alternative CR prescriptions, not a circular reduction. No equation or fitting step makes the headline alpha_CO range, the r1/2,CO/r1/2,H2 ~ 0.7 values, or the suprathermal r31 values equal to any input by construction.

Assumptions & free parameters 7 free parameters · 6 assumptions · 0 invented entities

The model depends on several calibrated normalizations (CR rate, ISRF, GMC pressure-radius relation, dust scaling) and domain assumptions about unresolved ISM physics. None of these is a new physical entity; they are extrapolations of local calibrations to z=1-3. The CR normalization is the most fragile because the paper uses two different values and states the line ratios are sensitive to it.

free parameters (7)
  • CR ionization rate normalization (ζ_-16, ξ_CR,MW) = ζ_-16=0.1, ξ_CR,MW=10^-16 s^-1 (Sec 2.2.3); 10^-17 s^-1 (Sec 4.2)
    Sets CR heating in DESPOTIC; drives high-J CO excitation and suprathermal line ratios. The paper states two different values.
  • Solar-neighborhood SFR surface density Σ_SFR,MW = 790 M_sun Myr^-1 kpc^-2
    Normalizes the localized UV ISRF strength chi in Eq. 8; controls cloud heating and CO luminosity.
  • GMC pressure-radius normalization (P_ext/k_B, M_C) = 10^4 cm^-3 K and 290 M_sun
    Eq. 2 sets cloud radius from external pressure and cloud mass; determines density and line excitation.
  • Dust absorption coefficient κ_abs and β_UV,MW = κ_abs=1.078e5 cm^2 g^-1; β_UV,MW≈0.5
    Eqs. 6-7 set UV attenuation and ISRF; calibrated to MW extinction curve and gas surface density.
  • Dust-to-gas ratio relation coefficients = log DGR = 2.445 log(Z/Zsun) - 2.029; DMR_MW=0.44
    Scales dust abundance and dust-gas coupling (Eqs. 10-11), affecting CO formation and gas heating.
  • Virial parameter α_vir = 1
    Assumes clouds are virialized to set velocity dispersion and clumping; affects CO line brightness.
  • Nearest-neighbor count for ISRF = 64
    Number of neighboring particles used to compute local SFR density and ISRF; chosen after comparing number vs distance thresholds.
assumptions (6)
  • domain assumption Each SIMBA gas particle represents a virialized molecular cloud with a power-law density profile (Eq. 3) and pressure-radius relation (Eq. 2).
    Sub-resolution ISM model; if z=1-3 cloud populations differ from local GMCs, all CO maps and alpha_CO gradients change.
  • domain assumption DESPOTIC equilibrium chemistry, thermal balance, and escape-probability radiative transfer correctly predict CO line luminosities for such clouds.
    Required to compute CO(1-0) through CO(5-4); extrapolated from z=0 validation to z=1-3.
  • ad hoc to paper The CR ionization rate scales linearly with local SFR surface density (Eq. 9).
    Boosts high-J excitation in centers; authors note starbursts may have only moderately enhanced CR rates.
  • ad hoc to paper The ISRF strength is set by the 64 nearest-neighbor particles with dust attenuation (Eqs. 6-8).
    Choice of neighbor count affects chi and hence CO luminosity; motivated by agreement with z=0 GMC data.
  • domain assumption SIMBA's sub-resolution H2 prescription (Krumholz & Gnedin 2011) and dust evolution correctly trace molecular gas at z=1-3.
    Input H2 masses and metallicities set the denominator of alpha_CO and the chemistry.
  • domain assumption X-ray/AGN heating is negligible for J≤5 CO lines in the selected main-sequence galaxies.
    Authors exclude AGN; X-rays can dominate high-J CO excitation in AGN hosts.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Predicting the resolved CO emission of $z=1-3$ star-forming galaxies." pith.science (2026). https://pith.science/paper/TMQRGUCX

@misc{pith2026250613899,
  author       = {Pith},
  title        = {Pith review of: Predicting the resolved CO emission of $z=1-3$ star-forming galaxies},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TMQRGUCX}},
  note         = {Machine review of arXiv:2506.13899}
}
abstract

(Abridged) Resolved observations of the CO emission from $z=1-3$ star-forming galaxies are becoming increasingly common, with new high-resolution surveys on the horizon. We aim to inform the interpretation of this resolved CO emission by creating synthetic observations and testing to what extent routinely observed CO transitions can be used to trace H$_2$ across galaxy disks. To this end, we extract $z=1-3$ massive star-forming galaxies (on and above the main sequence) from the SIMBA cosmological simulation and predict their spatially resolved CO(1$-$0)-to-CO(5$-$4) emission using the $\texttt{SLICK}$ pipeline, which combines sub-resolution modeling of the cloud population with the DESPOTIC spectral line calculation code. We find that the CO(1$-$0)-to-H$_2$ ratio ($\alpha_{\rm CO}$) varies significantly within these galaxy disks$-$from values of $\sim1-5$ $\mathrm{M_\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ in the central 1-3 kpc of the most massive galaxies to $>100$ $\mathrm{M_\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$ at $\sim$ 15 kpc. Thus, the use of a single $\alpha_{\rm CO}$ to derive the H$_2$ surface density leads to severe underestimates of the H$_2$ contribution in its outskirts. As expected, higher-$J$ CO transitions trace molecular gas in the centers at higher densities, whereas CO(1$-$0) better traces the more diffuse, extended molecular gas. We see significant variations in the CO excitation, with CO(3$-$2)/CO(1$-$0) line luminosity ratios of the most massive galaxies at $z\sim2$ declining from $\sim$ 0.9 in the galaxy centers to $\sim$ 0.1 in the outskirts. On average, line ratios increase substantially toward higher redshifts and lower galaxy stellar masses. We predict that tracing molecular gas with CO beyond 3-5 kpc of cosmic noon galaxies will be challenging with current facilities due to the drastic increase in $\alpha_{\rm CO}$.

Figures

Figures reproduced from arXiv: 2506.13899 by the authors.

Figure 1
Figure 1. Different treatments of the interstellar radiation field, as shown through the correlation between the clouds’ H2 mass vs. CO(1–0) luminosity: the left-hand panel shows a non-localized implementation (i.e., same ISRF for all clouds in a galaxy) whereas the right-hand panel shows the localized treatment (considering the 64 nearest-neighboring clouds). Slick’s clouds are indicated by the hexagonal bins, where the colo… view at source ↗
Figure 2
Figure 2. Example of the resolved CO emission and H [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Average radial profiles of CO surface brightness and H [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (10 more)
Figure 4
Figure 4. Figure 4: Average radial profile of the CO-to-H2 conversion factor (αCO), for different transitions (shown in the legend on the right). The x axis is plotted in log scale. The solid lines represent the median αCO, averaged across the entire galaxy sample shown in each panel. The…
Figure 5
Figure 5. Figure 5: Average radial profile of the CO line luminosity ratios, in the same bins of redshift (columns) and stellar mass (rows) as the [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: Histogram of the size ratio of CO(1–0) half-light radii [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]
Figure 7
Figure 7. Figure 7: Histograms of the size ratios of (a) CO(3–2) half-light [PITH_FULL_IMAGE:figures/full_fig_p011_7.png]
Figure 8
Figure 8. Figure 8: CO half-light radii as a function of galaxy stellar mass, in three redshift bins (one for each panel, as labelled at top), for [PITH_FULL_IMAGE:figures/full_fig_p012_8.png]
Figure 9
Figure 9. Figure 9: Average radial profile of abundances of HI, H [PITH_FULL_IMAGE:figures/full_fig_p012_9.png]
Figure 10
Figure 10. Figure 10: αCO vs. gas-phase metallicity per pixel, including the pixels of all synthetic CO(1–0) maps for the entire z=1–3 sample. The color-coding represents the number of pixels per hexbin as shown via a color bar on the right of the figure. The columns again show different r…
Figure 11
Figure 11. Figure 11: αCO vs. UV field strength, χ, per pixel. The column and row separation, and color-coding are the same as for [PITH_FULL_IMAGE:figures/full_fig_p013_11.png]
Figure 12
Figure 12. Figure 12: αCO vs. hydrogen number density per pixel. The column and row separation, and color-coding are the same as for [PITH_FULL_IMAGE:figures/full_fig_p014_12.png]
Figure 13
Figure 13. Figure 13: αCO vs. non-thermal velocity dispersion per pixel. The column and row separation, and color-coding are the same as for [PITH_FULL_IMAGE:figures/full_fig_p014_13.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Vz-GAL Dusty Star-Forming Galaxies: Revisiting the CO-H2 Conversion Factor Tension

    astro-ph.GA 2026-07 conditional novelty 6.0 of 10

    The CO-to-H2 conversion factor in high-z dusty star-forming galaxies is not required to be 0.8; intermediate to near-Galactic values (median ~4-5) remain dynamically viable.

  2. The Arizona Molecular ISM Survey with the SMT: The Diverse Carbon Monoxide Line Ratios and Spectral Line Energy Distributions of Star Forming Galaxies

    astro-ph.GA 2025-07 accept novelty 6.0 of 10

    CO line ratios in 47 nearby star-forming galaxies rise with star formation rate, leading to empirical prescriptions for converting higher-J CO luminosities to molecular gas masses.

Reference graph

Works this paper leans on

90 extracted references · 48 canonical work pages · cited by 2 Pith papers

  1. [1]

    B., Narayanan, D., Whitaker, K

    Akins, H. B., Narayanan, D., Whitaker, K. E., et al. 2022, ApJ, 929, 94 ALMA Partnership, Vlahakis, C., Hunter, T. R., et al. 2015, ApJ, 808, L4

  2. [2]

    A., Wagg, J., et al

    Aravena, M., Hodge, J. A., Wagg, J., et al. 2014, MNRAS, 442, 558

  3. [3]

    G., van Dishoeck, E

    Bisbas, T. G., van Dishoeck, E. F., Hu, C.-Y ., & Schruba, A. 2023, MNRAS, 519, 729

  4. [4]

    D., Leroy, A

    Bolatto, A. D., Leroy, A. K., Rosolowsky, E., Walter, F., & Blitz, L. 2008, ApJ, 686, 948

  5. [5]

    D., Warren, S

    Bolatto, A. D., Warren, S. R., Leroy, A. K., et al. 2015, ApJ, 809, 175

  6. [6]

    D., Wolfire, M., & Leroy, A

    Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207

  7. [7]

    & Bica, E

    Bonatto, C. & Bica, E. 2011, MNRAS, 415, 2827

  8. [8]

    A., van der Werf, P., Weiss, A., et al

    Boogaard, L. A., van der Werf, P., Weiss, A., et al. 2020, ApJ, 902, 109

Show all 90 references
  1. [9]

    2015, A&A, 575, A56 Calistro Rivera, G., Hodge, J

    Bournaud, F., Daddi, E., Weiß, A., et al. 2015, A&A, 575, A56 Calistro Rivera, G., Hodge, J. A., Smail, I., et al. 2018, ApJ, 863, 56

  2. [10]

    2002, Chinese J

    Chang, R.-X., Shu, C.-G., & Hou, J.-L. 2002, Chinese J. Astron. Astrophys., 2, 226

  3. [11]

    A., Smail, I., et al

    Chen, C.-C., Hodge, J. A., Smail, I., et al. 2017, ApJ, 846, 108

  4. [12]

    M., Chastenet, J., et al

    Chiang, I.-D., Sandstrom, K. M., Chastenet, J., et al. 2024, ApJ, 964, 18

  5. [13]

    K., Hayward, C

    Cochrane, R. K., Hayward, C. C., Anglés-Alcázar, D., et al. 2019, MNRAS, 488, 1779 Davé, R., Anglés-Alcázar, D., Narayanan, D., et al. 2019, MNRAS, 486, 2827 Davé, R., Crain, R. A., Stevens, A. R. H., et al. 2020, MNRAS, 497, 146

  6. [14]

    2001, A&A, 373, 853

    Dumke, M., Nieten, C., Thuma, G., Wielebinski, R., & Walsh, W. 2001, A&A, 373, 853

  7. [15]

    2024, MNRAS, 527, 8727 Frias Castillo, M., Hodge, J., Rybak, M., et al

    Esposito, F., Vallini, L., Pozzi, F., et al. 2024, MNRAS, 527, 8727 Frias Castillo, M., Hodge, J., Rybak, M., et al. 2023, ApJ, 945, 128

  8. [16]

    2024, ApJ, 974, 197

    Garcia, K., Narayanan, D., Popping, G., et al. 2024, ApJ, 974, 197

  9. [17]

    C., & Wolfire, M

    Gong, M., Ostriker, E. C., & Wolfire, M. G. 2017, ApJ, 843, 38 González-López, J., Barrientos, L. F., Gladders, M. D., et al. 2017, ApJ, 846, L22

  10. [18]

    Greve, T. R. & Sommer-Larsen, J. 2008, A&A, 480, 335

  11. [19]

    C., Weiss, A., Yun, M

    Harrington, K. C., Weiss, A., Yun, M. S., et al. 2021, ApJ, 908, 95

  12. [20]

    Hopkins, P. F. 2015, MNRAS, 450, 53

  13. [21]

    K., Belfiore, F., Lelli, F., et al

    Hunt, L. K., Belfiore, F., Lelli, F., et al. 2023, A&A, 675, A64

  14. [22]

    2022, ApJ, 933, 11

    Ikeda, R., Tadaki, K.-i., Iono, D., et al. 2022, ApJ, 933, 11

  15. [23]

    A., Gerin, M., et al

    Indriolo, N., Neufeld, D. A., Gerin, M., et al. 2012, ApJ, 758, 83

  16. [24]

    2020, MNRAS, 498, 5960

    Inoue, S., Yoshida, N., & Yajima, H. 2020, MNRAS, 498, 5960

  17. [25]

    2020, ApJ, 899, 37

    Kaasinen, M., Walter, F., Novak, M., et al. 2020, ApJ, 899, 37

  18. [26]

    2017, MNRAS, 471, 2917

    Kamenetzky, J., Rangwala, N., & Glenn, J. 2017, MNRAS, 471, 2917

  19. [27]

    C., Richings, A

    Keating, L. C., Richings, A. J., Murray, N., et al. 2020, MNRAS, 499, 837

  20. [28]

    Krumholz, M. R. 2014, MNRAS, 437, 1662

  21. [29]

    R., Crocker, R

    Krumholz, M. R., Crocker, R. M., & Offner, S. S. R. 2023, MNRAS, 520, 5126

  22. [30]

    Krumholz, M. R. & Gnedin, N. Y . 2011, ApJ, 729, 36

  23. [31]

    Lagos, C. d. P., Bayet, E., Baugh, C. M., et al. 2012, MNRAS, 426, 2142

  24. [32]

    G., et al

    Lelli, F., Zhang, Z.-Y ., Bisbas, T. G., et al. 2023, A&A, 672, A106

  25. [33]

    K., Rosolowsky, E., Usero, A., et al

    Leroy, A. K., Rosolowsky, E., Usero, A., et al. 2022, ApJ, 927, 149

  26. [34]

    2019, MNRAS, 490, 1425

    Li, Q., Narayanan, D., & Davé, R. 2019, MNRAS, 490, 1425

  27. [35]

    Lower, S., Narayanan, D., Hu, C.-Y ., & Privon, G. C. 2024, ApJ, 965, 123

  28. [36]

    2023, ApJ, 950, 94

    Lower, S., Narayanan, D., Li, Q., & Davé, R. 2023, ApJ, 950, 94

  29. [37]

    C., Cormier, D., Hony, S., et al

    Madden, S. C., Cormier, D., Hony, S., et al. 2020, A&A, 643, A141

  30. [38]

    Meijerink, R., Spaans, M., & Israel, F. P. 2007, A&A, 461, 793

  31. [39]

    F., & van der Werf, P

    Meijerink, R., Spaans, M., Loenen, A. F., & van der Werf, P. P. 2011, A&A, 525, A119

  32. [40]

    D., Koda, J., Momose, R., et al

    Meyer, J. D., Koda, J., Momose, R., et al. 2011, The Astrophysical Journal, 744, 42, publisher: The American Astronomical Society

  33. [41]

    D., Koda, J., Momose, R., et al

    Meyer, J. D., Koda, J., Momose, R., et al. 2013, The Astrophysical Journal, 772, 107, publisher: The American Astronomical Society Muñoz, J. A. & Furlanetto, S. R. 2013, MNRAS, 435, 2676

  34. [42]

    2007, PASJ, 59, 43

    Muraoka, K., Kohno, K., Tosaki, T., et al. 2007, PASJ, 59, 43

  35. [43]

    2016, PASJ, 68, 18

    Muraoka, K., Takeda, M., Yanagitani, K., et al. 2016, PASJ, 68, 18

  36. [44]

    J., Hayward, C

    Narayanan, D., Cox, T. J., Hayward, C. C., Younger, J. D., & Hernquist, L. 2009, MNRAS, 400, 1919

  37. [45]

    C., & Hernquist, L

    Narayanan, D., Krumholz, M., Ostriker, E. C., & Hernquist, L. 2011, MNRAS, 418, 664

  38. [46]

    & Krumholz, M

    Narayanan, D. & Krumholz, M. R. 2014, MNRAS, 442, 1411

  39. [47]

    & Krumholz, M

    Narayanan, D. & Krumholz, M. R. 2017, MNRAS, 467, 50

  40. [48]

    R., Ostriker, E

    Narayanan, D., Krumholz, M. R., Ostriker, E. C., & Hernquist, L. 2012, MN- RAS, 421, 3127

  41. [49]

    P., Burkhart, B., Mac Low, M.-M., et al

    Olsen, K. P., Burkhart, B., Mac Low, M.-M., et al. 2021, ApJ, 922, 88

  42. [50]

    P., Greve, T

    Olsen, K. P., Greve, T. R., Brinch, C., et al. 2016, MNRAS, 457, 3306

  43. [51]

    Osterbrock, D. E. & Ferland, G. J. 2006, Astrophysics of gaseous nebulae and active galactic nuclei

  44. [52]

    P., van der Werf, P

    Papadopoulos, P. P., van der Werf, P. P., Xilouris, E. M., et al. 2012, MNRAS, 426, 2601

  45. [53]

    L., Ott, J., Klein, U., et al

    Pineda, J. L., Ott, J., Klein, U., et al. 2009, ApJ, 703, 736 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, A&A, 594, A13

  46. [54]

    S., Faisst, A

    Popping, G., Narayanan, D., Somerville, R. S., Faisst, A. L., & Krumholz, M. R. 2019, MNRAS, 482, 4906

  47. [55]

    P., Spaans, M., Trager, S

    Popping, G., Pérez-Beaupuits, J. P., Spaans, M., Trager, S. C., & Somerville, R. S. 2014, MNRAS, 444, 1301

  48. [56]

    2022, MNRAS, 510, 3321

    Popping, G., Pillepich, A., Calistro Rivera, G., et al. 2022, MNRAS, 510, 3321

  49. [57]

    2016, MNRAS, 461, 93

    Popping, G., van Kampen, E., Decarli, R., et al. 2016, MNRAS, 461, 93

  50. [58]

    2019, ApJ, 877, L23

    Puglisi, A., Daddi, E., Liu, D., et al. 2019, ApJ, 877, L23

  51. [59]

    2021, MNRAS, 508, 5217

    Puglisi, A., Daddi, E., Valentino, F., et al. 2021, MNRAS, 508, 5217

  52. [60]

    C., et al

    Ramambason, L., Lebouteiller, V ., Madden, S. C., et al. 2024, A&A, 681, A14

  53. [61]

    2012, ApJ, 757, 155 Rémy-Ruyer, A., Madden, S

    Rebolledo, D., Wong, T., Leroy, A., Koda, J., & Donovan Meyer, J. 2012, ApJ, 757, 155 Rémy-Ruyer, A., Madden, S. C., Galliano, F., et al. 2014, A&A, 563, A31

  54. [62]

    J., Faucher-Giguère, C.-A., Gurvich, A

    Richings, A. J., Faucher-Giguère, C.-A., Gurvich, A. B., Schaye, J., & Hayward, C. C. 2022, MNRAS, 517, 1557

  55. [63]

    Richings, A. J. & Schaye, J. 2016, MNRAS, 458, 270

  56. [64]

    A., Hodge, J., Walter, F., Carilli, C

    Riechers, D. A., Hodge, J., Walter, F., Carilli, C. L., & Bertoldi, F. 2011, ApJ, 739, L31

  57. [65]

    2023, A&A, 679, A129

    Rizzo, F., Roman-Oliveira, F., Fraternali, F., et al. 2023, A&A, 679, A129

  58. [66]

    Rosenberg, M. J. F., van der Werf, P. P., Aalto, S., et al. 2015, ApJ, 801, 72

  59. [67]

    A., Vegetti, S., et al

    Rybak, M., Hodge, J. A., Vegetti, S., et al. 2020, MNRAS, 494, 5542

  60. [68]

    T., Castillo, M

    Rybak, M., Jansen, J. T., Castillo, M. F., et al. 2024 [arXiv:2411.06474]

  61. [69]

    K., et al

    Saito, T., Iono, D., Xu, C. K., et al. 2017, ApJ, 835, 174

  62. [70]

    M., Leroy, A

    Sandstrom, K. M., Leroy, A. K., Walter, F., et al. 2013, ApJ, 777, 5

  63. [71]

    2024, A&A, 682, A98

    Schimek, A., Decataldo, D., Shen, S., et al. 2024, A&A, 682, A98

  64. [72]

    K., Walter, F., et al

    Schruba, A., Leroy, A. K., Walter, F., et al. 2011, AJ, 142, 37

  65. [73]

    E., Baker, A

    Sharon, C. E., Baker, A. J., Harris, A. I., & Thomson, A. P. 2013, ApJ, 765, 6

  66. [74]

    E., Tagore, A

    Sharon, C. E., Tagore, A. S., Baker, A. J., et al. 2019, ApJ, 879, 52

  67. [75]

    D., Bryan, G

    Smith, B. D., Bryan, G. L., Glover, S. C. O., et al. 2017, MNRAS, 466, 2217

  68. [76]

    S., Aravena, M., Marrone, D

    Spilker, J. S., Aravena, M., Marrone, D. P., et al. 2015, ApJ, 811, 124

  69. [77]

    M., Riechers, D

    Stanley, F., Jones, B. M., Riechers, D. A., et al. 2023, ApJ, 945, 24

  70. [78]

    M., Papadopoulos, P

    Swinbank, A. M., Papadopoulos, P. P., Cox, P., et al. 2011, ApJ, 742, 11

  71. [79]

    2023, ApJ, 957, L15

    Tadaki, K.-i., Kodama, T., Koyama, Y ., et al. 2023, ApJ, 957, L15

  72. [80]

    J., et al

    Tadaki, K.-i., Kodama, T., Nelson, E. J., et al. 2017, ApJ, 841, L25

  73. [81]

    J., Swinbank, A

    Taylor, D. J., Swinbank, A. M., Smail, I., et al. 2025, MNRAS, 536, 1149

  74. [82]

    A., Richings, A

    Thompson, O. A., Richings, A. J., Gibson, B. K., et al. 2024, MNRAS, 532, 1948

  75. [83]

    2015, SPHGR: Smoothed-Particle Hydrodynamics Galaxy Re- duction, Astrophysics Source Code Library, record ascl:1502.012

    Thompson, R. 2015, SPHGR: Smoothed-Particle Hydrodynamics Galaxy Re- duction, Astrophysics Source Code Library, record ascl:1502.012

  76. [84]

    G., Smith, R

    Tress, R. G., Smith, R. J., Sormani, M. C., et al. 2020, MNRAS, 492, 2973

  77. [85]

    2020, A&A, 641, A155

    Valentino, F., Daddi, E., Puglisi, A., et al. 2020, A&A, 641, A155

  78. [86]

    2018, MNRAS, 473, 271

    Vallini, L., Pallottini, A., Ferrara, A., et al. 2018, MNRAS, 473, 271

  79. [87]

    Vallini, L., Tielens, A. G. G. M., Pallottini, A., et al. 2019, MNRAS, 490, 4502 van der Werf, P. P., Isaak, K. G., Meijerink, R., et al. 2010, A&A, 518, L42

  80. [88]

    E., et al

    Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261

  81. [89]

    2011, ApJS, 197, 16

    Wong, T., Hughes, A., Ott, J., et al. 2011, ApJS, 197, 16

  82. [90]

    2023, PASJ, 75, 743 Article number, page 17 of 20 A&A proofs: manuscript no

    Yasuda, A., Kuno, N., Sorai, K., et al. 2023, PASJ, 75, 743 Article number, page 17 of 20 A&A proofs: manuscript no. main Appendix A: Physical properties driving changes in line excitation In this section, we show the dependence of the CO(3–2) /CO(1–0) line ratio on cloud prop...

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.