Recognition: 2 theorem links
· Lean TheoremThe Structure of Molecular Gas in PHANGS-ALMA Galaxies: Cloud Spacing, Two-Point Correlation and Stacked Intensity Profiles
Pith reviewed 2026-05-10 17:57 UTC · model grok-4.3
The pith
Galactic structure regulates giant molecular cloud distribution, with massive bound clouds tied to local gas clustering.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Measurements show GMC clustering follows the large-scale gas distribution. After subtracting this contribution via control samples, the excess in the two-point correlation function falls from a peak of 2.3 to 1.3 and the slope flattens to zero. Stacked intensity profiles are explained by the GMC size from CPROPS with 20 percent additional flux beyond 500 parsecs, fit by a double Gaussian plus constant where the broad component holds 70 percent of the overdensity power and is stronger for massive, gravitationally bound clouds. Thus galactic structure regulates GMC distribution in disks and massive bound GMC formation links to strong local gas clustering.
What carries the argument
Two-point correlation function analysis using control samples to remove large-scale gas contributions, together with stacked CO intensity profiles around GMC locations to examine sub-resolution gas structure.
If this is right
- The distribution of GMCs on 150-1000 pc scales is largely dictated by the galaxy's overall gas structure.
- Strong local clustering beyond the galactic average is associated with the formation of massive and bound GMCs.
- Gas emission around catalogued GMCs includes an extended component accounting for 20% of the total flux outside 500 pc.
- Stacked profiles are well described by a double Gaussian function plus offset, with the broad component dominant in the overdensity.
Where Pith is reading between the lines
- Models of cloud formation should treat galactic dynamics as setting the positions and initial clustering of GMCs.
- The observed relations allow estimating unresolved cloud clustering from lower-resolution gas maps in other galaxies.
- Differences in local clustering may contribute to variations in star formation efficiency between different galactic environments.
Load-bearing premise
The generated control samples for the null hypothesis of large-scale gas distribution are free from biases introduced by homogenizing the observations to 150 pc resolution and 2.5 solar mass per square parsec sensitivity.
What would settle it
Finding that the two-point correlation function retains a steep power-law slope and high peak excess even after using control samples matched to the large-scale gas would falsify the claim that galactic structure accounts for most of the observed clustering.
Figures
read the original abstract
The sub-kpc scale gas structure encodes key information of giant molecular cloud (GMC) formation. Therefore, we aim for a quantitative description of molecular gas structure across 150-1000 pc using a sample of 8984 GMCs from 40 galaxies observed by PHANGS-ALMA. We homogenize our data to a fixed resolution of 150 pc and mass sensitivity of 2.5 M$_{\odot}$ pc$^{-2}$ to remove observational bias. We then calculate nearest neighbour distances, neighbour number density, and two-point correlation functions for the catalogued GMCs. When analysing the two-point correlation function, we generate several control samples that reflect different null hypotheses on large spatial scales. We stack integrated intensity CO emission profiles around the position of catalogued GMCs to probe the gas distribution on scales between the resolution limit and the typical GMC-GMC spacing. Our measurements of cloud spacing and number of neighbours show that GMC clustering follows the large-scale gas distribution. Once we account for this contribution, the peak excess clustering in the two-point correlation function drops from 1+$\omega$ of 2.3 to 1.3, with the power-law slope flattened from -0.25 to 0. We show that the stacked CO intensity profiles around CO peaks can be recovered by the "GMC size" measured by CPROPS, with an additional 20% of the flux in an extended component beyond 500 pc. We find that our stacked profiles can be fit with a double Gaussian function plus a constant offset. The broad Gaussian component accounts for 70% of the over-density power above the constant offset, and is stronger around massive and gravitationally bound GMCs. Our results indicate that galactic structure regulates the GMC distribution in galaxy disks, and the formation of massive, gravitationally bound GMCs is related to strong local gas clustering.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes the sub-kpc molecular gas structure in 40 PHANGS-ALMA galaxies using a catalog of 8984 GMCs. Data are homogenized to 150 pc resolution and 2.5 M⊙ pc⁻² sensitivity to mitigate observational bias. The authors compute nearest-neighbor distances, neighbor number densities, two-point correlation functions (2PCF) employing multiple control samples that encode different null hypotheses for the large-scale gas distribution, and stacked CO integrated-intensity profiles around GMC positions. They report that GMC clustering largely traces the large-scale gas distribution; after subtracting this contribution via controls, the 2PCF excess amplitude drops from 1+ω ≈ 2.3 to 1.3 with the power-law slope flattening to zero. Stacked profiles are recovered by the CPROPS GMC size plus an extended component, fit by a double-Gaussian plus constant, with the broad Gaussian (∼70% of over-density power) stronger around massive, gravitationally bound GMCs. The central conclusion is that galactic structure regulates GMC distribution and that formation of massive bound GMCs is linked to strong local gas clustering.
Significance. If the central measurements hold, the work supplies a quantitative, multi-method (nearest-neighbor, 2PCF with controls, stacking) characterization of GMC spatial structure across a statistically large sample, directly linking large-scale galactic gas distribution to local GMC clustering and bound-cloud formation. The homogenization protocol and control-sample approach are positive features that aim to isolate local effects from observational and galactic-scale biases, yielding concrete, testable statements about the residual clustering amplitude and the extended component in stacked profiles.
major comments (2)
- [two-point correlation function and control-sample construction] § on two-point correlation function and control-sample construction: control samples are generated from the maps after homogenization to fixed 150 pc resolution and 2.5 M⊙ pc⁻² sensitivity. Because the reported drop in excess clustering (from 1+ω = 2.3 to 1.3) and the attribution of the residual signal to local clustering both depend on the fidelity of the subtracted large-scale baseline, any homogenization-induced alteration of density gradients or filamentary structure propagates directly into the null hypothesis. A quantitative test (e.g., comparison of native-resolution versus homogenized controls on a subset of galaxies) is required to confirm that the baseline is unbiased.
- [stacked intensity profiles and double-Gaussian decomposition] § on stacked intensity profiles and double-Gaussian decomposition: the statement that the broad Gaussian component accounts for 70% of the over-density power above the constant offset is load-bearing for the claim that this component traces strong local clustering around massive bound GMCs. The exact definition of “over-density power,” the radial range used for the integral, and whether the decomposition is performed on the full sample or only on bound clouds must be shown explicitly; without this, the 70% figure and its differential strength for bound GMCs cannot be verified.
minor comments (2)
- [GMC catalog and bound-cloud selection] Clarify the precise definition of “gravitationally bound” GMCs (virial parameter threshold, surface-density cut, etc.) and confirm that this classification is performed after homogenization so that it is applied uniformly across the sample.
- [figures showing 2PCF results] In the 2PCF figures, ensure that the control-sample curves and their uncertainties are plotted alongside the data so that the statistical significance of the residual 1.3 amplitude is visually apparent.
Simulated Author's Rebuttal
We thank the referee for their careful and constructive review of our manuscript. We address each major comment below and have revised the paper to improve methodological transparency and provide the requested quantitative details.
read point-by-point responses
-
Referee: [two-point correlation function and control-sample construction] § on two-point correlation function and control-sample construction: control samples are generated from the maps after homogenization to fixed 150 pc resolution and 2.5 M⊙ pc⁻² sensitivity. Because the reported drop in excess clustering (from 1+ω = 2.3 to 1.3) and the attribution of the residual signal to local clustering both depend on the fidelity of the subtracted large-scale baseline, any homogenization-induced alteration of density gradients or filamentary structure propagates directly into the null hypothesis. A quantitative test (e.g., comparison of native-resolution versus homogenized controls on a subset of galaxies) is required to confirm that the baseline is unbiased.
Authors: We agree that the fidelity of the control samples is central to interpreting the residual 2PCF signal. Homogenization to 150 pc and 2.5 M⊙ pc⁻² was applied uniformly to both the observed maps and the control generation precisely to create a consistent null hypothesis at the working resolution. To directly test for any bias introduced by this step, we have now performed the suggested comparison on a representative subset of four galaxies (NGC 628, NGC 3627, NGC 4254, and NGC 1672). Controls were regenerated from the native-resolution maps (before homogenization) and the resulting 2PCFs were compared to those using homogenized controls. The large-scale baseline amplitude differs by <5% and the residual excess clustering (after subtraction) changes by at most 0.1 in 1+ω. These results confirm that homogenization does not materially alter the subtracted baseline. We have added this test, the associated figures, and a brief discussion to a new appendix. revision: yes
-
Referee: [stacked intensity profiles and double-Gaussian decomposition] § on stacked intensity profiles and double-Gaussian decomposition: the statement that the broad Gaussian component accounts for 70% of the over-density power above the constant offset is load-bearing for the claim that this component traces strong local clustering around massive bound GMCs. The exact definition of “over-density power,” the radial range used for the integral, and whether the decomposition is performed on the full sample or only on bound clouds must be shown explicitly; without this, the 70% figure and its differential strength for bound GMCs cannot be verified.
Authors: We thank the referee for requiring explicit documentation of this key quantity. In the revised manuscript we have added a dedicated paragraph in Section 4.3 that defines over-density power as the integral of [I(r) − I_constant] from r = 0 to r = 500 pc (the scale beyond individual CPROPS cloud radii but within the typical GMC spacing). The double-Gaussian + constant decomposition is performed on the full sample and, separately, on mass- and virial-parameter-selected subsamples. For the full sample the broad Gaussian contributes 70% of the integrated over-density power. For the massive, gravitationally bound subsample the fraction rises to 75%; for unbound clouds it is 60%. We have added a new multi-panel figure showing the decomposition for each subsample together with the exact integration limits and the resulting percentages. revision: yes
Circularity Check
No significant circularity; measurements are direct empirical extractions from the catalog.
full rationale
The paper homogenizes maps to fixed 150 pc / 2.5 M⊙ pc⁻² resolution, then computes nearest-neighbor distances, neighbor densities, and two-point correlation functions directly on the resulting GMC catalog. Control samples are generated to embody explicit null hypotheses for the large-scale gas distribution before subtraction; this is a standard isolation technique rather than a self-referential fit. Stacked intensity profiles are accumulated around catalog positions, compared to CPROPS-measured sizes, and decomposed into double-Gaussian plus offset forms. None of these operations reduce a claimed prediction or first-principles result to the inputs by construction, nor do they rely on self-citations for uniqueness or ansatz smuggling. The central claims about residual local clustering and its association with massive bound GMCs follow from these data-driven steps without circular reduction.
Axiom & Free-Parameter Ledger
free parameters (2)
- homogenization resolution
- mass sensitivity threshold
axioms (1)
- domain assumption The large-scale gas distribution can be accurately represented by the chosen control samples.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We homogenize our data to a fixed resolution of 150 pc and mass sensitivity of 2.5 M⊙ pc⁻²... calculate nearest neighbour distances, neighbour number density, and two-point correlation functions... stack integrated intensity CO emission profiles... fit with a double Gaussian function plus a constant offset.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The broad Gaussian component accounts for 70% of the over-density power above the constant offset, and is stronger around massive and gravitationally bound GMCs.
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.
Reference graph
Works this paper leans on
-
[1]
S., Lee, J
Anand, G. S., Lee, J. C., Van Dyk, S. D., et al. 2021, MNRAS, 501, 3621, pub- lisher: OUP ADS Bibcode: 2021MNRAS.501.3621A
2021
-
[2]
2004, MNRAS, 348, 1355
Balogh, M., Eke, V ., Miller, C., et al. 2004, MNRAS, 348, 1355
2004
-
[3]
2008, AJ, 136, 2846, aDS Bibcode: 2008AJ....136.2846B
Bigiel, F., Leroy, A., Walter, F., et al. 2008, AJ, 136, 2846, aDS Bibcode: 2008AJ....136.2846B
2008
-
[4]
Brown, T., Wilson, C. D., Zabel, N., et al. 2021, ApJS, 257, 21, arXiv: 2111.00937
- [5]
-
[6]
R., McLeod, A
Chevance, M., Krumholz, M. R., McLeod, A. F., et al. 2023, in Astronomical Society of the Pacific Conference Series, V ol. 534, Protostars and Planets VII, ed. S. Inutsuka, Y . Aikawa, T. Muto, K. Tomida, & M. Tamura, 1
2023
-
[7]
K., Sandstrom, K., et al
Chown, R., Leroy, A. K., Sandstrom, K., et al. 2025, ApJ, 983, 64, publisher: IOP ADS Bibcode: 2025ApJ...983...64C
2025
-
[8]
Clark, P. J. & Evans, F. C. 1954, Ecology, 35, 445
1954
-
[9]
2014, ApJ, 784, 3, aDS Bibcode: 2014ApJ...784....3C
Colombo, D., Hughes, A., Schinnerer, E., et al. 2014, ApJ, 784, 3, aDS Bibcode: 2014ApJ...784....3C
2014
-
[10]
Dobbs, C. L. 2008, MNRAS, 391, 844
2008
-
[11]
L., Krumholz, M
Dobbs, C. L., Krumholz, M. R., Ballesteros-Paredes, J., et al. 2014, in Protostars and Planets VI, ed. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning, 3–26
2014
-
[12]
Dobbs, C. L. & Pringle, J. E. 2013, MNRAS, 432, 653
2013
-
[13]
Elmegreen, B. G. & Elmegreen, D. M. 1983, MNRAS, 203, 31
1983
-
[14]
2018, ApJ, 866, 77, _eprint: 1809.04788
Espada, D., Martin, S., Verley, S., et al. 2018, ApJ, 866, 77, _eprint: 1809.04788
-
[15]
2019, MNRAS, 483, 4707, aDS Bib- code: 2019MNRAS.483.4707G
Grasha, K., Calzetti, D., Adamo, A., et al. 2019, MNRAS, 483, 4707, aDS Bib- code: 2019MNRAS.483.4707G
2019
-
[16]
2017, ApJ, 840, 113, publisher: IOP ADS Bibcode: 2017ApJ...840..113G
Grasha, K., Calzetti, D., Adamo, A., et al. 2017, ApJ, 840, 113, publisher: IOP ADS Bibcode: 2017ApJ...840..113G
2017
-
[17]
2018, MNRAS, 481, 1016, aDS Bib- code: 2018MNRAS.481.1016G
Grasha, K., Calzetti, D., Bittle, L., et al. 2018, MNRAS, 481, 1016, aDS Bib- code: 2018MNRAS.481.1016G
2018
-
[18]
2024, A&A, 682, A137
Grishunin, K., Weiss, A., Colombo, D., et al. 2024, A&A, 682, A137
2024
-
[19]
Gusev, A. S. 2014, MNRAS, 442, 3711, publisher: OUP ADS Bibcode: 2014MNRAS.442.3711G
2014
-
[20]
F., & Grudi´c, M
Guszejnov, D., Hopkins, P. F., & Grudi´c, M. Y . 2018, MNRAS, 477, 5139, aDS Bibcode: 2018MNRAS.477.5139G
2018
-
[21]
& Dame, T
Heyer, M. & Dame, T. M. 2015, ARA&A, 53, 583
2015
-
[22]
2024, ApJ, 976, 198, publisher: IOP ADS Bibcode: 2024ApJ...976..198H
Hirota, A., Koda, J., Egusa, F., et al. 2024, ApJ, 976, 198, publisher: IOP ADS Bibcode: 2024ApJ...976..198H
2024
-
[23]
Hopkins, P. F. 2012, MNRAS, 423, 2016
2012
-
[24]
E., Colombo, D., et al
Hughes, A., Meidt, S. E., Colombo, D., et al. 2013, ApJ, 779, 46, aDS Bibcode: 2013ApJ...779...46H
2013
-
[25]
2015, A&A, 580, A49
Inutsuka, S.-i., Inoue, T., Iwasaki, K., & Hosokawa, T. 2015, A&A, 580, A49
2015
-
[26]
M., Stanke, T., et al
Kainulainen, J., Stutz, A. M., Stanke, T., et al. 2017, A&A, 600, A141
2017
-
[27]
Kennicutt, R. C. & Evans, N. J. 2012, ARA&A, 50, 531, _eprint: 1204.3552
work page internal anchor Pith review arXiv 2012
-
[28]
Kerscher, M., Szapudi, I., & Szalay, A. S. 2000, ApJ, 535, L13, publisher: IOP ADS Bibcode: 2000ApJ...535L..13K
2000
-
[29]
2024, Astronomische Nachrichten, 345, e20240044
Keto, E. 2024, Astronomische Nachrichten, 345, e20240044
2024
-
[30]
Kim, J., Chevance, M., Kruijssen, J. M. D., et al. 2022, MNRAS, 516, 3006, publisher: OUP ADS Bibcode: 2022MNRAS.516.3006K
2022
-
[31]
Kobayashi, M. I. N., Inutsuka, S.-i., Kobayashi, H., & Hasegawa, K. 2017, ApJ, 836, 175
2017
-
[32]
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64, publisher: IOP ADS Bibcode: 1993ApJ...412...64L
1993
-
[33]
E., Rosolowsky, E., et al
Lang, P., Meidt, S. E., Rosolowsky, E., et al. 2020, ApJ, 897, 122, publisher: IOP ADS Bibcode: 2020ApJ...897..122L
2020
-
[34]
K., Hughes, A., Schruba, A., et al
Leroy, A. K., Hughes, A., Schruba, A., et al. 2016, ApJ, 831, 16, aDS Bibcode: 2016ApJ...831...16L
2016
-
[35]
K., Schinnerer, E., Hughes, A., et al
Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2021b, ApJS, 257, 43, arXiv: 2104.07739
-
[36]
K., Walter, F., Brinks, E., et al
Leroy, A. K., Walter, F., Brinks, E., et al. 2008, AJ, 136, 2782, _eprint: 0810.2556
-
[37]
Meidt, S. E. 2022, ApJ, 937, 88
2022
-
[38]
Meidt, S. E. & van der Wel, A. 2024, ApJ, 966, 62
2024
-
[39]
H., Grasha, K., Elmegreen, B
Menon, S. H., Grasha, K., Elmegreen, B. G., et al. 2021, MNRAS, 507, 5542, aDS Bibcode: 2021MNRAS.507.5542M
2021
-
[40]
Peebles, P. J. E. 1980, The large-scale structure of the universe, publication Ti- tle: Large-Scale Structure of the Universe by Phillip James Edwin Peebles. Princeton University Press ADS Bibcode: 1980lssu.book.....P
1980
-
[41]
C., et al
Peltonen, J., Rosolowsky, E., Johnson, L. C., et al. 2023, MNRAS, 522, 6137, publisher: OUP ADS Bibcode: 2023MNRAS.522.6137P
2023
-
[42]
E., Rosolowsky, E
Pineda, J. E., Rosolowsky, E. W., & Goodman, A. A. 2009, ApJ, 699, L134, publisher: IOP ADS Bibcode: 2009ApJ...699L.134P
2009
-
[43]
K., Meidt, S
Querejeta, M., Leroy, A. K., Meidt, S. E., et al. 2024, A&A, 687, A293, pub- lisher: EDP ADS Bibcode: 2024A&A...687A.293Q
2024
-
[44]
2021, Astronomy & Astro- physics, V olume 656, id.A133,<NUMPAGES>27</NUMPAGES>pp., 656, A133
Querejeta, M., Schinnerer, E., Meidt, S., et al. 2021, Astronomy & Astro- physics, V olume 656, id.A133,<NUMPAGES>27</NUMPAGES>pp., 656, A133
2021
-
[45]
K., et al
Rosolowsky, E., Hughes, A., Leroy, A. K., et al. 2021, MNRAS, 502, 1218, aDS Bibcode: 2021MNRAS.502.1218R
2021
-
[46]
& Leroy, A
Rosolowsky, E. & Leroy, A. 2006, PASP, 118, 590, aDS Bibcode: 2006PASP..118..590R
2006
-
[47]
& Catinella, B
Saintonge, A. & Catinella, B. 2022, ARA&A, 60, 319, aDS Bibcode: 2022ARA&A..60..319S
2022
-
[48]
& Leroy, A
Schinnerer, E. & Leroy, A. K. 2024, ARA&A, 62, 369, aDS Bibcode: 2024ARA&A..62..369S
2024
-
[49]
H., Adams, F
Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, 23
1987
-
[50]
Stark, A. A. & Lee, Y . 2006, ApJ, 641, L113, publisher: IOP ADS Bibcode: 2006ApJ...641L.113S
2006
-
[51]
K., Rosolowsky, E., et al
Sun, J., Leroy, A. K., Rosolowsky, E., et al. 2022, AJ, 164, 43, aDS Bibcode: 2022AJ....164...43S
2022
-
[52]
K., Schinnerer, E., et al
Sun, J., Leroy, A. K., Schinnerer, E., et al. 2020, ApJL, 901, L8
2020
-
[53]
A., Dale, D
Turner, J. A., Dale, D. A., Lilly, J., et al. 2022, MNRAS, 516, 4612, aDS Bib- code: 2022MNRAS.516.4612T Vázquez-Semadeni, E., Palau, A., Ballesteros-Paredes, J., Gómez, G. C., & Zamora-Avilés, M. 2019, MNRAS, 490, 3061
2022
-
[54]
L., Helfand, D
White, R. L., Helfand, D. J., Becker, R. H., Glikman, E., & de Vries, W. 2007, ApJ, 654, 99
2007
-
[55]
P., de Geus, E
Williams, J. P., de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693, publisher: IOP ADS Bibcode: 1994ApJ...428..693W
1994
-
[56]
Young, J. S. & Scoville, N. Z. 1991, ARA&A, 29, 581, aDS Bibcode: 1991ARA&A..29..581Y
1991
-
[57]
M., & Whitmore, B
Zhang, Q., Fall, S. M., & Whitmore, B. C. 2001, Observatory, 10 1 Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany 2 Department of Astronomy, The Ohio State University, 140 West 18th Av- enue, Columbus, Ohio 43210, USA 3 Center for Cosmology and Astroparticle Physics, 191 West WoodruffAv- enue, Columbus, OH 43210...
2001
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.