pith. machine review for the scientific record. sign in

arxiv: 2604.06354 · v3 · submitted 2026-04-07 · 🌌 astro-ph.GA

Recognition: 2 theorem links

· Lean Theorem

Multi-scale Gas Structure and Dynamics in an Extragalactic Central Molecular Zone

Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:18 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords molecular gascentral molecular zoneNGC 3351ALMAvelocity dispersionfree-fall timedendrogramgas dynamics
0
0 comments X

The pith

Gravitational free-fall times match ordered gas crossing times on all scales in an extragalactic central molecular zone, while random motions cross faster below 10 parsecs.

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

This paper maps the hierarchical structure of molecular gas in the central molecular zone of NGC 3351 using high-resolution ALMA observations. It separates the motions inside each gas structure into ordered flows versus random components and compares both to the time gas would take to collapse under its own gravity. Ordered motions appear stronger on scales above 30 parsecs while random motions dominate below that transition. The free-fall time stays comparable to the ordered crossing time across the full range of sizes examined, yet both exceed the random crossing time at the smallest scales. These timescale comparisons show how different physical drivers operate at different sizes in a Milky Way-like galactic center.

Core claim

Using ALMA CO(3-2) data at 5 parsec resolution and the dendrogram method to identify structures across two decades in size, the observed velocity dispersion is decomposed into ordered and random parts. Ordered motions dominate above 30 parsecs while random motions become dominant below that scale. Modulo conversion-factor uncertainties, the gravitational free-fall time remains comparable to the ordered-motion crossing time at every scale, and both become longer than the random-motion crossing time at scales below 10 parsecs.

What carries the argument

Decomposition of velocity dispersion inside dendrogram-identified molecular structures into ordered versus random components, followed by direct comparison of those crossing times against the gravitational free-fall time derived from structure mass and size.

Load-bearing premise

The factor that converts observed CO intensity into molecular gas mass is assumed accurate enough that its uncertainties do not reverse the relative ordering of free-fall times and crossing times.

What would settle it

An independent mass measurement, for example from dust emission calibrated without CO, that yields free-fall times substantially shorter than ordered crossing times at scales above 30 parsecs.

Figures

Figures reproduced from arXiv: 2604.06354 by Elias Oakes, Eric W. Koch, Erik Rosolowsky, Jiayi Sun, Liam M. Wang, Sanghyuk Moon, Yu-Hsuan Teng.

Figure 1
Figure 1. Figure 1: (Top) CO (3–2) moment-0 and moment-1 maps for the CMZ of NGC 3351. Our following structural analysis focuses on the outer star-forming ring (outside the hatched region; see section 2), where one sees the richest gas structure hierarchy. (Bottom) A dendrogram of the hierarchical gas structures across NGC 3351’s star-forming ring. Vertical lines near the top of each structure tree represents local maxima in … view at source ↗
Figure 2
Figure 2. Figure 2: CO (3–2) moment-1 maps (left), model velocity gradient maps (middle), and velocity residual maps (right) for two identified structures in the dendrogram. The velocity gradient model and its residual provides one way to distinguish ordered versus random motions within each structure (see subsection 2.4). The top row shows a large trunk (structure #208, rphy = 47 pc) located on the southwest side of NGC 3351… view at source ↗
Figure 3
Figure 3. Figure 3: (Left) The size–mass relation among all identified gas structures, with the best-fit power law relation Mmol ∝ r 2.54 phy shown by a black line. The vertical shaded region shows where the measured size becomes unreliable due to spatial resolution limit. Besides, the gas mass measurements for smaller structures may also be inaccurate (gray dotted line) due to the coarser spatial resolution of the αCO(3−2) c… view at source ↗
Figure 4
Figure 4. Figure 4: (Top) Estimated velocity dispersion of random motions as a function of structure size. The extent of each vertical bar represents the upper and lower limits on σrand (see subsection 2.4). (Middle) Similar to the top panel, but shows the effective velocity dispersion of ordered motion, σord. (Bottom) Ratio between σrand and σord (in logarithmic scale) as a function of structure size. The gray data points ma… view at source ↗
Figure 5
Figure 5. Figure 5: (Left) Random motion crossing time (orange), ordered motion crossing time (red), and gravitational free-fall time (blue) as functions of structure size. As in [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: A comparison of our CO data to a pure galactic rotation model. The top row displays the CO moment 1 maps at different scales while the bottom row displays the galactic rotation model. We find that for larger structures, the galactic rotation model seems to describe the observed velocity gradients well, whereas in smaller structures, the model fails to capture the observed velocity gradient. Our interpretat… view at source ↗
Figure 7
Figure 7. Figure 7: (Left) Velocity gradient amplitude between the best-fit model and the galactic rotation model as a function of radius. We find that for larger structures (> 10pc) the ratio approaches unity, which indicates that galactic rotation alone dominates the modeled ordered motion. For smaller structures, the galactic rotation model tends to underestimate the velocity gradient amplitude, suggesting substantial cont… view at source ↗
read the original abstract

The structures and dynamics of the interstellar medium are governed by a combination of self-gravity, external gravity, and various sources of ordered and random motions on different spatial scales. This paper uses ALMA CO (3-2) observations at 0.1" $\approx$ 5 pc resolution to examine the scale dependence of molecular gas structure and dynamics in the central molecular zone (CMZ) of a nearby galaxy, NGC 3351. We use the dendrogram technique to characterize hierarchical molecular gas structures spanning two decades in spatial scales and measure their size, gas mass, and velocity dispersion. Their size-linewidth relation shows a power-law slope of 0.58, comparable to measurements for CMZs in other galaxies and suggestive of significant contribution from ordered motion on large scales. We further decompose the observed velocity dispersion in each gas structure into ordered versus random motions. The former appears stronger in gas structures at $\gtrsim$ 30 pc while the latter becomes more dominant at $\lesssim$ 30 pc. Modulo uncertainties with the CO-to-H$_2$ conversion factor, the estimated gravitational free-fall time is comparable to the crossing time of ordered motions for structures on all spatial scales, and both becomes longer than the crossing time of random motions at small, $\lesssim$ 10 pc scales. Our results highlight the varying sources and drivers of gas motions on different spatial scales in the CMZ of a Milky Way-like galaxy.

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.

Referee Report

2 major / 2 minor

Summary. The paper analyzes ALMA CO(3-2) observations at ~5 pc resolution of the CMZ in NGC 3351. Using dendrograms, it identifies hierarchical molecular structures spanning two decades in scale, measures their sizes, masses, and velocity dispersions, and reports a size-linewidth power-law slope of 0.58. Velocity dispersions are decomposed into ordered and random components, with ordered motions dominating at ≳30 pc and random motions at ≲30 pc. Modulo X_CO uncertainties, gravitational free-fall times are found comparable to ordered-motion crossing times on all scales and longer than random-motion crossing times at ≲10 pc scales.

Significance. If robust, the results provide direct observational constraints on the scale-dependent balance between self-gravity, ordered motions, and random motions in an extragalactic CMZ analogous to the Milky Way. The size-linewidth slope and ordered/random decomposition are directly supported by the ALMA data and standard methods, extending similar analyses from other galaxies. The time-scale comparisons, if confirmed against X_CO variations, would strengthen understanding of what drives gas dynamics across scales in such environments.

major comments (2)
  1. [§4] §4 (time-scale comparisons): the central claim that t_ff is comparable to t_cross,ordered on all scales (and exceeds t_cross,random at ≲10 pc) rests on masses derived from a single X_CO conversion factor applied to integrated CO intensity. Because t_ff ∝ 1/sqrt(ρ) and ρ ∝ M, a uniform factor-of-2 shift in X_CO scales all t_ff values by ~1.4 while leaving both crossing times unchanged. No propagated uncertainties, error bands on the time ratios, or sensitivity runs with alternative X_CO values are presented, so the qualitative statements of 'comparable' and 'longer' remain untested against the uncertainties explicitly flagged in the abstract.
  2. [§3.2] §3.2 (velocity dispersion decomposition): while the method for separating ordered and random motions is described as standard, the paper does not quantify how the decomposition threshold or spatial filtering choices affect the reported transition scale (~30 pc) or the subsequent time-scale ordering. This choice is load-bearing for the claim that random motions dominate only at small scales.
minor comments (2)
  1. [Abstract] Abstract: the clause 'both becomes longer than the crossing time of random motions' contains a subject-verb agreement error and should read 'both become longer'.
  2. [Throughout] Figure captions and text: several instances of 'CO-to-H2' use inconsistent subscript formatting; standardize to CO-to-H₂ throughout.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments, which have helped us identify areas where the manuscript can be strengthened. We address each major comment below and describe the revisions we will implement.

read point-by-point responses
  1. Referee: [§4] §4 (time-scale comparisons): the central claim that t_ff is comparable to t_cross,ordered on all scales (and exceeds t_cross,random at ≲10 pc) rests on masses derived from a single X_CO conversion factor applied to integrated CO intensity. Because t_ff ∝ 1/sqrt(ρ) and ρ ∝ M, a uniform factor-of-2 shift in X_CO scales all t_ff values by ~1.4 while leaving both crossing times unchanged. No propagated uncertainties, error bands on the time ratios, or sensitivity runs with alternative X_CO values are presented, so the qualitative statements of 'comparable' and 'longer' remain untested against the uncertainties explicitly flagged in the abstract.

    Authors: We agree that a quantitative assessment of X_CO uncertainties is needed to support the time-scale claims. In the revised manuscript, we will add a sensitivity analysis varying X_CO by a factor of 2 (consistent with typical extragalactic CMZ uncertainties), propagate these into the free-fall and crossing time ratios, and include error bands on the relevant figures and text. This will test and confirm that the qualitative statements of comparability and ordering hold within the expected range. revision: yes

  2. Referee: [§3.2] §3.2 (velocity dispersion decomposition): while the method for separating ordered and random motions is described as standard, the paper does not quantify how the decomposition threshold or spatial filtering choices affect the reported transition scale (~30 pc) or the subsequent time-scale ordering. This choice is load-bearing for the claim that random motions dominate only at small scales.

    Authors: We will strengthen the robustness analysis of the velocity dispersion decomposition. The revised manuscript will include additional tests varying the spatial filtering scale and the threshold separating ordered versus random components. These tests will demonstrate that the ~30 pc transition scale remains stable, with the results presented in a new subsection or appendix to support the scale-dependent dominance of motions. revision: yes

Circularity Check

0 steps flagged

No circularity: all quantities derived from direct ALMA measurements using standard formulas

full rationale

The paper identifies hierarchical structures with dendrograms on ALMA CO(3-2) data, measures sizes, masses (via integrated intensity and a fixed X_CO), and velocity dispersions, then decomposes the latter into ordered/random components and computes t_ff and crossing times with the usual expressions t_ff ∝ 1/sqrt(G rho) and t_cross = R/v. None of these steps defines a quantity in terms of itself, fits a parameter on a subset then re-predicts a related quantity from the same data, or relies on a self-citation chain for the central time-scale comparison. The X_CO uncertainty is explicitly flagged and scales all masses uniformly, but this is an external systematic, not a circular reduction within the paper's equations.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The analysis rests on the standard dendrogram algorithm for hierarchical structure identification and the CO-to-H2 conversion factor for mass estimates; no new entities are postulated.

free parameters (1)
  • CO-to-H2 conversion factor
    Converts observed CO(3-2) intensity to H2 mass; uncertainties directly affect mass and free-fall time estimates.
axioms (1)
  • domain assumption Dendrogram algorithm correctly identifies physically meaningful hierarchical molecular structures
    Invoked when applying the technique to the ALMA cube to extract size, mass, and velocity dispersion across scales.

pith-pipeline@v0.9.0 · 5584 in / 1441 out tokens · 43092 ms · 2026-05-10T19:18:39.949198+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

38 extracted references · 38 canonical work pages

  1. [1]

    S., Lee, J

    Anand, G. S., Lee, J. C., Van Dyk, S. D., et al. 2021, MNRAS, 501, 3621, doi: 10.1093/mnras/staa3668 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, ˚ a, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, The Astronomical Journal, 156, 123, doi: 10.3847/1538-3881/aab...

  2. [2]

    2008, ApJ, 686, 948, doi: 10.1086/591513

    Blitz, L. 2008, ApJ, 686, 948, doi: 10.1086/591513

  3. [3]

    2000, , 543, 822, 10.1086/317122

    Burkert, A., & Bodenheimer, P. 2000, ApJ, 543, 822, doi: 10.1086/317122

  4. [4]

    J., Shivaei, I., et al

    Calzetti, D., Battisti, A. J., Shivaei, I., et al. 2021, ApJ, 913, 37, doi: 10.3847/1538-4357/abf118

  5. [5]

    2023, MNRAS, 522, 4078, doi: 10.1093/mnras/stad1211

    Choi, W., Liu, L., Bureau, M., et al. 2023, MNRAS, 522, 4078, doi: 10.1093/mnras/stad1211

  6. [6]

    2024, MNRAS, 531, 4045, doi: 10.1093/mnras/stae1394

    Choi, W., Bureau, M., Liu, L., et al. 2024, MNRAS, 531, 4045, doi: 10.1093/mnras/stae1394

  7. [7]

    2021,, 2.0.0 Zenodo, doi: 10.5281/zenodo.4905459

    Comrie, A., et al. 2021, CARTA: The Cube Analysis and Rendering Tool for Astronomy (v2.0.0), doi: 10.5281/zenodo.4905459

  8. [8]

    2010, ARA&A, 48, 547, doi: 10.1146/annurev-astro-081309-130854

    Fukui, Y., & Kawamura, A. 2010, ARA&A, 48, 547, doi: 10.1146/annurev-astro-081309-130854

  9. [9]

    2019, radio-astro-tools/spectral-cube: Release v0.4.5, v0.4.5, Zenodo, doi: 10.5281/zenodo.3558614

    Ginsburg, A., Koch, E., Robitaille, T., et al. 2019, radio-astro-tools/spectral-cube: Release v0.4.5, v0.4.5, Zenodo, doi: 10.5281/zenodo.3558614

  10. [10]

    A., Benson, P

    Goodman, A. A., Benson, P. J., Fuller, G. A., & Myers, P. C. 1993, ApJ, 406, 528, doi: 10.1086/172465

  11. [11]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2

  12. [12]

    Heyer, M., & Dame, T. M. 2015, ARA&A, 53, 583, doi: 10.1146/annurev-astro-082214-122324

  13. [13]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi: 10.1109/MCSE.2007.55

  14. [14]

    2011, The Astrophysical Journal, 732, 78, doi: 10.1088/0004-637X/732/2/78

    Imara, N., & Blitz, L. 2011, The Astrophysical Journal, 732, 78, doi: 10.1088/0004-637X/732/2/78

  15. [15]

    D., Koch, E

    Krieger, N., Bolatto, A. D., Koch, E. W., et al. 2020, ApJ, 899, 158, doi: 10.3847/1538-4357/aba903

  16. [16]

    E., Rosolowsky, E., et al

    Lang, P., Meidt, S. E., Rosolowsky, E., et al. 2020, ApJ, 897, 122, doi: 10.3847/1538-4357/ab9953

  17. [17]

    Larson, R. B. 1981, MNRAS, 194, 809, doi: 10.1093/mnras/194.4.809

  18. [18]

    , keywords =

    Leroy, A. K., Bolatto, A. D., Ostriker, E. C., et al. 2015, ApJ, 801, 25, doi: 10.1088/0004-637X/801/1/25

  19. [19]

    K., Hughes, A., Liu, D., et al

    Leroy, A. K., Schinnerer, E., Hughes, A., et al. 2021, The Astrophysical Journal Supplement Series, 257, 43, doi: 10.3847/1538-4365/abec80

  20. [20]

    2021, MNRAS, 505, 4048, doi: 10.1093/mnras/stab1537

    Liu, L., Bureau, M., Blitz, L., et al. 2021, MNRAS, 505, 4048, doi: 10.1093/mnras/stab1537

  21. [21]

    F., & Ostriker, E

    McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565, doi: 10.1146/annurev.astro.45.051806.110602

  22. [22]

    E., Leroy, A

    Meidt, S. E., Leroy, A. K., Rosolowsky, E., et al. 2018, ApJ, 854, 100, doi: 10.3847/1538-4357/aaa290

  23. [23]

    K., Faesi, C

    Oakes, E. K., Faesi, C. M., Rosolowsky, E., et al. 2025, The Astrophysical Journal, 993, 193, doi: 10.3847/1538-4357/ae0642

  24. [24]

    Robitaille, T. P. 2019, Astrodendro: Dendrograms for Astronomical Data

  25. [25]

    P., & Bressert, E

    Robitaille, T. P., & Bressert, E. 2012, APLpy: Astronomical Plotting Library in Python

  26. [26]

    Goodman, A. A. 2008, ApJ, 679, 1338, doi: 10.1086/587685

  27. [27]

    A., Prandoni, I., et al

    Ruffa, I., Davis, T. A., Prandoni, I., et al. 2019, MNRAS, 489, 3739, doi: 10.1093/mnras/stz2368

  28. [28]

    Schinnerer, E., & Leroy, A. K. 2024, ARA&A, 62, 369, doi: 10.1146/annurev-astro-071221-052651

  29. [29]

    , archivePrefix = "arXiv", eprint =

    Shetty, R., Beaumont, C. N., Burton, M. G., Kelly, B. C., & Klessen, R. S. 2012, MNRAS, 425, 720, doi: 10.1111/j.1365-2966.2012.21588.x

  30. [30]

    , keywords =

    Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730, doi: 10.1086/165493

  31. [31]

    and Schruba, Andreas and Rosolowsky, Erik and Hughes, Annie and Kruijssen, J

    Sun, J., Leroy, A. K., Schruba, A., et al. 2018, ApJ, 860, 172, doi: 10.3847/1538-4357/aac326

  32. [32]

    K., Schinnerer, E., et al

    Sun, J., Leroy, A. K., Schinnerer, E., et al. 2020, ApJL, 901, L8, doi: 10.3847/2041-8213/abb3be 13

  33. [33]

    K., Rosolowsky, E., et al

    Sun, J., Leroy, A. K., Rosolowsky, E., et al. 2022, The Astronomical Journal, 164, 43, doi: 10.3847/1538-3881/ac74bd

  34. [34]

    ApJ967(2), 133 (2024) https://doi.org/10.3847/1538-4357/ad3de6 arXiv:2401.14453 [astro-ph.GA]

    Sun, J., He, H., Batschkun, K., et al. 2024, The Astrophysical Journal, 967, 133, doi: 10.3847/1538-4357/ad3de6

  35. [35]

    M., Sun, J., et al

    Teng, Y.-H., Sandstrom, K. M., Sun, J., et al. 2022, ApJ, 925, 72, doi: 10.3847/1538-4357/ac382f

  36. [36]

    , keywords =

    Teng, Y.-H., Sandstrom, K. M., Sun, J., et al. 2023, The Astrophysical Journal, 950, 119, doi: 10.3847/1538-4357/accb86 V´ azquez-Semadeni, E., Palau, A., Ballesteros-Paredes, J., G´ omez, G. C., & Zamora-Avil´ es, M. 2019, MNRAS, 490, 3061, doi: 10.1093/mnras/stz2736

  37. [37]

    Virtanen, R

    Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261, doi: 10.1038/s41592-019-0686-2

  38. [38]

    2025, The Astrophysical Journal Letters, 983, L21, doi: 10.3847/2041-8213/adc095

    Xie, Y.-H., & Li, G.-X. 2025, The Astrophysical Journal Letters, 983, L21, doi: 10.3847/2041-8213/adc095