arxiv: 2604.09353 · v1 · submitted 2026-04-10 · 🌌 astro-ph.GA

Recognition: no theorem link

Constraining the Molecular Kennicutt-Schmidt Relation with Multi-Transition CO Observations of Nearby Galaxies

Victoria G. G. Samboco , Ryan P. Keenan

Authors on Pith no claims yet

Pith reviewed 2026-05-10 17:06 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords Kennicutt-Schmidt relationmolecular gasstar formation rateCO transitionsdense gasgalaxy evolutionnearby galaxies

0 comments

The pith

Higher-J CO transitions produce increasingly linear molecular Kennicutt-Schmidt relations in nearby galaxies, with slopes falling from 1.26 to 1.07.

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

The paper measures the star formation rate surface density against molecular gas surface density using three CO lines that trace gas of increasing density. It finds the power-law index of the relation drops from 1.26 for CO(1-0) to 1.07 for CO(3-2), approaching linearity as the tracer samples denser gas. This trend implies that the underlying link between gas volume density and star formation volume density follows a power-law index near 1.5, so denser gas converts to stars more efficiently. The result refines how galaxy evolution models should connect total molecular gas to star formation by showing that tracer choice and density matter.

Core claim

Using a uniform sample of 36 nearby galaxies observed in CO(1-0), CO(2-1), and CO(3-2), the authors apply binning plus MCMC fitting to derive the slope, intercept, and scatter of the Σ_SFR–Σ_CO relation for each transition. The measured slopes are 1.26, 1.14, and 1.07 respectively. The systematic decrease toward unity is interpreted as evidence that denser molecular gas is more directly coupled to ongoing star formation, consistent with an underlying volume-density power-law index of approximately 1.5.

What carries the argument

Power-law fitting via binning and MCMC on surface-density maps of three CO transitions that progressively sample denser molecular gas.

If this is right

Denser molecular gas is more directly linked to star formation than lower-density gas.
Star formation efficiency rises in higher-density environments.
Choice of molecular tracer changes the apparent form of the Kennicutt-Schmidt relation.
Results align with earlier linear relations found using HCN or high-J CO.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Star formation prescriptions in galaxy simulations may need explicit density thresholds rather than total molecular gas.
Extending the same analysis to CO(4-3) or higher transitions could test whether the slope continues to approach 1.0.
The inferred volume-density index of 1.5 could be checked against hydrodynamic simulations that resolve individual clouds.

Load-bearing premise

The different CO lines cleanly trace increasing gas densities without major excitation biases or selection effects that would distort the measured slopes.

What would settle it

Independent maps showing no decrease in slope across the same three CO transitions, or resolved volume-density measurements yielding a power-law index far from 1.5.

Figures

Figures reproduced from arXiv: 2604.09353 by Ryan P. Keenan, Victoria G. G. Samboco.

**Figure 1.** Figure 1: The distribution of our galaxy sample in stellar mass and SFR. Our full sample of 36 galaxies is represented by blue boxes. Black outlines denote the subsample of massive, star-forming galaxies (Stellar Mass ≥ 1010 M⊙, SFR≥ 1 M⊙yr−1 ) used at points in this work. These are sources detected in CO(1–0), CO(2–1), and CO(3–2). The gray region represents the star forming galaxy main sequence (Speagle et al. 2… view at source ↗

**Figure 2.** Figure 2: Relationship between SFR surface density (ΣSFR) and molecular gas surface density (ΣCO) for the full galaxy sample, as traced by CO(1–0) (left panel), CO(2–1) (middle panel), and CO(3–2) (right panel). The Black dashed line indicates the best-fitting power law relation, and the thin red lines show 100 randomly selected fits from the MCMC chain. Individual galaxies (black circles) are used in the fit. The g… view at source ↗

**Figure 3.** Figure 3: SFR surface density as a function of molecular gas surface density for a subsample excluding low-mass and low-SFR galaxies. Data are shown for CO(1–0) (left panel), CO(2–1) (middle panel), and CO(3–2) (right panel). Individual galaxies (black circles) are fit with the black dashed line. Blue squares (binned data) are shown for illustrative purposes. These slopes are shallower than the typically observed sl… view at source ↗

**Figure 4.** Figure 4: Joint distribution of K–S slopes for the CO(1–0) and CO(3–2) lines across 50,000 bootstrap realizations of our data. In each realization the observed quantities (SFR, CO luminosities, and half-light radius) were perturbed according to their measurement uncertainties, and the K–S relation was re-fitted using orthogonal distance regression (ODR). The contours indicate the 68% (solid) and 95% (dashed) highest… view at source ↗

**Figure 5.** Figure 5: Observed molecular Kennicutt-Schmidt (K–S) index as a function of CO transition from CO(1–0) to CO(4-3). The shaded regions represent theoretical predictions from Narayanan et al. (2011) for different indices of the underlying (volumetric) SFR density-molecular gas density relation: N = 2 (green), N = 1.5 (orange), and N = 1 (blue). The blue star markers correspond to our data. The remaining data points re… view at source ↗

read the original abstract

The relationship between the star formation rate surface density and the molecular gas surface density in galaxies is key to understanding galaxy evolution. To investigate the molecular Kennicutt-Schmidt (K-S) relation and its dependence on gas density, we analyze a uniform sample of 36 nearby galaxies from the AMISS survey, focusing on the CO(1-0), CO(2-1), and CO(3-2) transitions, which trace progressively denser and warmer molecular gas. Using statistical methods that combine binning with Markov Chain Monte Carlo (MCMC) fitting, we derive the slope, scatter, and intercept of the $\Sigma_{\mathrm{SFR}}$-$\Sigma_{\mathrm{CO}}$ relation for each transition. We find power-law slopes of 1.26, 1.14, and 1.07 for CO(1-0), CO(2-1), and CO(3-2), respectively, consistent with a trend toward increasingly linear star formation relations at higher-J transitions. This behavior supports the idea that denser gas is more directly linked to ongoing star formation and is consistent with previous findings of near-linear correlations between HCN or high-J CO luminosities and global SFR. The observed trend suggests an underlying relation between gas and SFR volume densities with a power-law index of $\sim$1.5, indicating enhanced star formation efficiency in denser environments. These findings underscore the critical role of dense gas in regulating star formation and highlight the importance of tracer selection and excitation conditions when interpreting the K-S relation across different environments.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

New multi-transition slopes from a uniform 36-galaxy sample are useful empirical additions, but the jump to a ~1.5 volume-density index rests on an unstated mapping.

read the letter

The paper's core contribution is a set of consistent power-law fits to the molecular KS relation across CO(1-0), (2-1), and (3-2) in the same 36 nearby galaxies from AMISS. The slopes come out at 1.26, 1.14, and 1.07 respectively, with the expected trend toward linearity at higher-J lines. That uniformity across tracers is the real value; most prior work mixes different samples or single transitions, so having one dataset processed the same way gives a cleaner comparison. The binning-plus-MCMC approach is standard but applied evenly, which helps when comparing the relations directly. It also aligns with the broader picture from HCN and dense-gas tracers that star formation efficiency rises with density. The abstract frames this as evidence that denser gas links more directly to SFR, which is a reasonable reading of the slope trend. The soft spot is the claim that the surface-density slopes imply an underlying volume-density power-law index of ~1.5. No explicit derivation, disk-thickness model, or filling-factor assumption appears in the abstract or the stress-test note, so the step from observed slopes to that number is not transparent. Without seeing the full methods section, it is hard to judge how selection effects, excitation variations, or completeness were handled in the fits. The paper is aimed at people who model or observe the KS relation in nearby galaxies and want updated empirical anchors for dense-gas prescriptions. A reader working on simulation sub-grid recipes or multi-tracer surveys would find the slopes worth checking against their own data. It deserves peer review because the dataset and uniform fitting add a concrete constraint that is not already in the single-transition literature, even if the volume-density interpretation needs more scaffolding to stand on its own.

Referee Report

1 major / 4 minor

Summary. The paper analyzes multi-transition CO observations (CO(1-0), CO(2-1), CO(3-2)) of a uniform sample of 36 nearby galaxies from the AMISS survey to constrain the molecular Kennicutt-Schmidt relation. Using statistical binning combined with MCMC fitting, it reports power-law slopes of 1.26, 1.14, and 1.07 for the Σ_SFR–Σ_CO relations, respectively, with a trend toward linearity at higher-J transitions. The authors interpret this as evidence that denser gas is more directly linked to star formation and suggest an underlying volume-density power-law index of ∼1.5 between gas and SFR.

Significance. If the empirical slopes are robust and the physical interpretation can be placed on firmer footing, the work supplies useful observational constraints on how the KS relation varies with molecular-gas density tracer. The uniform sample and multi-transition approach, together with the MCMC fitting procedure, represent a clear methodological strength that allows direct comparison across density regimes. The results align with and extend prior findings on dense-gas tracers such as HCN, reinforcing the role of dense molecular gas in regulating star formation efficiency.

major comments (1)

[Abstract and §4] Abstract and §4 (or equivalent discussion section): The claim that the observed surface-density slopes 'suggest an underlying relation between gas and SFR volume densities with a power-law index of ∼1.5' is not accompanied by an explicit derivation, equation, or physical model. Connecting the measured surface-density power laws to a volume-density index requires a specific mapping (e.g., assumptions on how effective disk thickness, filling factor, or the width of the density PDF changes with the tracer's critical density). No such mapping is referenced or derived, rendering the headline physical conclusion unsupported even though the raw slope measurements are directly tied to the data.

minor comments (4)

[Methods] Methods section: The description of the binning procedure and MCMC implementation should include the number of bins, binning criteria (e.g., equal-number or equal-width), choice of priors, and how measurement uncertainties and upper limits are propagated into the posterior distributions.
[§2] Sample description (likely §2): Clarify the selection function for the 36 galaxies—specifically whether the sample is complete in all three CO transitions or whether non-detections in higher-J lines introduce bias—and report the detection fractions per transition.
[§3 and figures] Figure captions and §3: The plotted relations should explicitly state the fitting method, the adopted binning, and whether the reported slopes include or exclude galaxies with only upper limits; error bars on binned points should be shown.
[Discussion] Discussion: The quantitative comparison with previous HCN and high-J CO studies could be strengthened by tabulating literature slope values alongside the new measurements rather than qualitative statements of consistency.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive feedback. We address the single major comment below.

read point-by-point responses

Referee: [Abstract and §4] Abstract and §4 (or equivalent discussion section): The claim that the observed surface-density slopes 'suggest an underlying relation between gas and SFR volume densities with a power-law index of ∼1.5' is not accompanied by an explicit derivation, equation, or physical model. Connecting the measured surface-density power laws to a volume-density index requires a specific mapping (e.g., assumptions on how effective disk thickness, filling factor, or the width of the density PDF changes with the tracer's critical density). No such mapping is referenced or derived, rendering the headline physical conclusion unsupported even though the raw slope measurements are directly tied to the data.

Authors: We agree that the link between the measured surface-density slopes (1.26, 1.14, and 1.07) and an underlying volume-density power-law index of ∼1.5 is stated without an explicit derivation or referenced model in the current version. The statement is a qualitative inference drawn from the progressive approach to linearity as the tracer critical density increases, but we acknowledge that a rigorous mapping (via disk thickness, filling factor, or PDF width) is required to support it quantitatively. In the revised manuscript we will add a short derivation in §4 under the assumption of constant scale height combined with a density-dependent volume filling factor that narrows with increasing critical density; this yields an effective volume index of ∼1.5 that reproduces the observed surface-slope trend. We will also update the abstract to reference the added model and cite the relevant theoretical literature on turbulent density PDFs and the KS relation. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results are direct empirical fits

full rationale

The paper derives power-law slopes (1.26, 1.14, 1.07) via binning and MCMC fitting applied to new multi-transition CO observations of 36 galaxies. These are statistical outputs from the dataset and do not reduce to any prior fitted parameter or self-citation by the paper's own equations. The interpretive suggestion of an underlying volume-density power-law index of ~1.5 is presented as a qualitative inference from the observed trend in slopes, without any explicit mapping equation or self-referential derivation that would make it tautological. No load-bearing step matches the enumerated circularity patterns; the chain remains self-contained as analysis of fresh observational data.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on empirical power-law fits to binned data and standard astrophysical assumptions about what CO transitions trace; no new entities are postulated.

free parameters (1)

power-law slopes for each CO transition = 1.26, 1.14, 1.07
Fitted via MCMC to binned Σ_SFR-Σ_CO data for CO(1-0), (2-1), and (3-2).

axioms (1)

domain assumption Higher-J CO transitions trace progressively denser and warmer molecular gas
Invoked to interpret the observed trend toward linearity as evidence for density-dependent star formation efficiency.

pith-pipeline@v0.9.0 · 5589 in / 1259 out tokens · 48617 ms · 2026-05-10T17:06:39.396965+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

43 extracted references · 25 canonical work pages · 1 internal anchor

[1]

M., Maiolino, R., Bluck, A

Baker, W. M., Maiolino, R., Bluck, A. F., et al. 2022, Monthly Notices of the Royal Astronomical Society, 510, 3622

2022
[2]

2014, Astronomische Nachrichten, 335, 470

Bigiel, F., Cormier, D., & Schmidt, T. 2014, Astronomische Nachrichten, 335, 470

2014
[3]

2010, AJ, 140, 1194, doi: 10.1088/0004-6256/136/6/2846

Bigiel, F., Leroy, A., Walter, F., et al. 2008, The Astronomical Journal, 136, 2846–2871, doi: 10.1088/0004-6256/136/6/2846

work page doi:10.1088/0004-6256/136/6/2846 2008
[4]

and Wolfire, Mark and Leroy, Adam K

Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207, doi: 10.1146/annurev-astro-082812-140944

work page Pith review doi:10.1146/annurev-astro-082812-140944 2013
[5]

2014, A&A, 564, A65, doi: 10.1051/0004-6361/201322311

Boselli, A., Cortese, L., & Boquien, M. 2014, A&A, 564, A65, doi: 10.1051/0004-6361/201322311

work page doi:10.1051/0004-6361/201322311 2014
[6]

S., Wagg, J., Cicone, C., et al

Bothwell, M. S., Wagg, J., Cicone, C., et al. 2014, MNRAS, 445, 2599, doi: 10.1093/mnras/stu1936

work page doi:10.1093/mnras/stu1936 2014
[7]

Galactic

Chabrier, G. 2003, PASP, 115, 763, doi: 10.1086/376392 den Brok, J. S., Chatzigiannakis, D., Bigiel, F., et al. 2021, Monthly Notices of the Royal Astronomical Society, 504, 3221, doi: 10.1093/mnras/stab859

work page internal anchor Pith review doi:10.1086/376392 2003
[8]

W., Lang, D., & Goodman, J

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, Publications of the Astronomical Society of the Pacific, 125, 306

2013
[9]

J., et al

Freundlich, J., Combes, F., Tacconi, L. J., et al. 2019, Astronomy & Astrophysics, 622, A105

2019
[10]

Gao, Y., & Solomon, P. M. 2004, The Astrophysical Journal, 606, 271–290, doi: 10.1086/382999

work page doi:10.1086/382999 2004
[11]

Data analysis recipes: Fitting a model to data

Hogg, D. W., Bovy, J., & Lang, D. 2010, arXiv preprint arXiv:1008.4686

work page Pith review arXiv 2010
[12]

J., Werner, M

Hollenbach, D. J., Werner, M. W., & Salpeter, E. E. 1971, Astrophysical Journal, vol. 163, p. 165, 163, 165

1971
[13]

2017, MNRAS, 466, 4795, doi: 10.1093/mnras/stx046

Janowiecki, S., Catinella, B., Cortese, L., et al. 2017, MNRAS, 466, 4795, doi: 10.1093/mnras/stx046 Jim´ enez-Donaire, M. J., Bigiel, F., Leroy, A. K., et al. 2019, ApJ, 880, 127, doi: 10.3847/1538-4357/ab2b95

work page doi:10.1093/mnras/stx046 2017
[14]

P., Marrone, D

Keenan, R. P., Marrone, D. P., & Keating, G. K. 2025, The Astrophysical Journal, 979, 228

2025
[15]

P., Marrone, D

Keenan, R. P., Marrone, D. P., Keating, G. K., et al. 2024, ApJ, 975, 150, doi: 10.3847/1538-4357/ad7504

work page doi:10.3847/1538-4357/ad7504 2024
[16]

2020, ApJ, 889, 103, doi: 10.3847/1538-4357/ab6221

Lamperti, I., Saintonge, A., Koss, M., et al. 2020, ApJ, 889, 103, doi: 10.3847/1538-4357/ab6221

work page doi:10.3847/1538-4357/ab6221 2020
[17]

K., Walter, F., Brinks, E., et al

Leroy, A. K., Walter, F., Brinks, E., et al. 2008, The Astronomical Journal, 136, 2782–2845, doi: 10.1088/0004-6256/136/6/2782

work page doi:10.1088/0004-6256/136/6/2782 2008
[18]

K., Walter, F., Bigiel, F., et al

Leroy, A. K., Walter, F., Bigiel, F., et al. 2009, AJ, 137, 4670, doi: 10.1088/0004-6256/137/6/4670

work page doi:10.1088/0004-6256/137/6/4670 2009
[19]

K., Walter, F., Sandstrom, K., et al

Leroy, A. K., Walter, F., Sandstrom, K., et al. 2013, The Astronomical Journal, 146, 19

2013
[20]

K., Schinnerer, E., Hughes, A., et al

Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2021, ApJS, 257, 43, doi: 10.3847/1538-4365/ac17f3

work page doi:10.3847/1538-4365/ac17f3 2021
[21]

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

Leroy, A. K., Rosolowsky, E., Usero, A., et al. 2022, The Astrophysical Journal, 927, 149, doi: 10.3847/1538-4357/ac3490

work page doi:10.3847/1538-4357/ac3490 2022
[22]

L., et al

Lin, L., Pan, H.-A., Ellison, S. L., et al. 2019, ApJL, 884, L33, doi: 10.3847/2041-8213/ab4815 —. 2024, ApJ, 963, 115, doi: 10.3847/1538-4357/ad18b9

work page doi:10.3847/2041-8213/ab4815 2019
[23]

2011, A&A, 534, A102, doi: 10.1051/0004-6361/201117056

Lisenfeld, U., Espada, D., Verdes-Montenegro, L., et al. 2011, A&A, 534, A102, doi: 10.1051/0004-6361/201117056

work page doi:10.1051/0004-6361/201117056 2011
[24]

2015, The Astrophysical Journal, 810, L14, doi: 10.1088/2041-8205/810/2/l14 Molecular K–S Relation Across CO Transitions13

Liu, D., Gao, Y., Isaak, K., et al. 2015, The Astrophysical Journal, 810, L14, doi: 10.1088/2041-8205/810/2/l14 Molecular K–S Relation Across CO Transitions13

work page doi:10.1088/2041-8205/810/2/l14 2015
[25]

N., & Baars, J

Martin, R. N., & Baars, J. W. 1990, in Instrumentation in Astronomy VII, Vol. 1235, SPIE, 503–514

1990
[26]

F., & Ostriker, E

McKee, C. F., & Ostriker, E. C. 2007, Annu. Rev. Astron. Astrophys., 45, 565

2007
[27]

2017, ApJ, 837, 137, doi: 10.3847/1538-4357/aa6115

Morokuma-Matsui, K., & Muraoka, K. 2017, ApJ, 837, 137, doi: 10.3847/1538-4357/aa6115

work page doi:10.3847/1538-4357/aa6115 2017
[28]

2023, A&A, 678, A183, doi: 10.1051/0004-6361/202346951

Nagy, D., Dessauges-Zavadsky, M., Messa, M., et al. 2023, A&A, 678, A183, doi: 10.1051/0004-6361/202346951

work page doi:10.1051/0004-6361/202346951 2023
[29]

J., Hayward, C

Narayanan, D., Cox, T. J., Hayward, C. C., & Hernquist, L. 2011, Monthly Notices of the Royal Astronomical Society, 412, 287

2011
[30]

J., Bigiel, F., et al

Neumann, L., Gallagher, M. J., Bigiel, F., et al. 2023, MNRAS, 521, 3348, doi: 10.1093/mnras/stad424

work page doi:10.1093/mnras/stad424 2023
[31]

2021, A&A, 650, A134, doi: 10.1051/0004-6361/202140733

Pessa, I., Schinnerer, E., Belfiore, F., et al. 2021, A&A, 650, A134, doi: 10.1051/0004-6361/202140733

work page doi:10.1051/0004-6361/202140733 2021
[32]

E., & Kravtsov, A

Robertson, B. E., & Kravtsov, A. V. 2008, The Astrophysical Journal, 680, 1083–1111, doi: 10.1086/587796

work page doi:10.1086/587796 2008
[33]

M., et al

Roman-Duval, J., Heyer, M., Brunt, C. M., et al. 2016, The Astrophysical Journal, 818, 144

2016
[34]

2022, ARA&A, 60, 319, doi: 10.1146/annurev-astro-021022-043545

Saintonge, A., & Catinella, B. 2022, ARA&A, 60, 319, doi: 10.1146/annurev-astro-021022-043545

work page doi:10.1146/annurev-astro-021022-043545 2022
[35]

, archivePrefix = "arXiv", eprint =

Saintonge, A., Catinella, B., Tacconi, L. J., et al. 2017, ApJS, 233, 22, doi: 10.3847/1538-4365/aa97e0

work page doi:10.3847/1538-4365/aa97e0 2017
[36]

1999, The Astrophysical Journal, 523, 701 Salom´ e, Q., Krongold, Y., Longinotti, A

Sakamoto, S. 1999, The Astrophysical Journal, 523, 701 Salom´ e, Q., Krongold, Y., Longinotti, A. L., et al. 2023, Monthly Notices of the Royal Astronomical Society, 524, 3130

1999
[37]

M., Leroy, A., Walter, F., et al

Sandstrom, K. M., Leroy, A., Walter, F., et al. 2013, The Astrophysical Journal, 777, 5

2013
[38]

1959, Astrophysical Journal, vol

Schmidt, M. 1959, Astrophysical Journal, vol. 129, p. 243, 129, 243

1959
[39]

K., Walter, F., et al

Schruba, A., Leroy, A. K., Walter, F., et al. 2011, The Astronomical Journal, 142, 37

2011
[40]

C., Dullemond, C

Shetty, R., Glover, S. C., Dullemond, C. P., & Klessen, R. S. 2011, Monthly Notices of the Royal Astronomical Society, 412, 1686

2011
[41]

Silverman, J. D. 2014, The Astrophysical Journal Supplement Series, 214, 15

2014
[42]

D., Wong, T., et al

Utomo, D., Bolatto, A. D., Wong, T., et al. 2017, The Astrophysical Journal, 849, 26

2017
[43]

2021, PASJ, 73, 257, doi: 10.1093/pasj/psaa119

Yajima, Y., Sorai, K., Miyamoto, Y., et al. 2021, PASJ, 73, 257, doi: 10.1093/pasj/psaa119

work page doi:10.1093/pasj/psaa119 2021