Pith. sign in

REVIEW 3 major objections 5 minor 24 references

The Importance of Tidal Forces in Molecular Cloud Dynamics

T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A galaxy's tides can bind molecular clouds that look unbound, unbind clouds that look bound, and make the standard virial parameter misclassify their dynamics.

desk verdict The full-potential result is real, but the tidal-specific claim needs a control separating internal structure from external tides. read the letter →

arxiv 2507.03788 v1 pith:YWME6NOU submitted 2025-07-04 astro-ph.GA

JinWoo Lee , Alexa Saur , Mordecai-Mark Mac Low , Hui Li This is my paper · ORCID
classification astro-ph.GA
keywords molecularcloudstidalforcesvirialparametertheoremgalacticenvironmentstarformationnumericalsimulationstensor
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

This paper argues that the dynamical state of molecular clouds cannot be read off from their own gravity alone, because the tidal pull of the surrounding galaxy can either help compress a cloud or tear it apart. The authors analyze clouds identified with a dendrogram algorithm in zoom-in simulations of a Milky Way-mass galaxy, covering the galactic center, a solar-circle region, and the outskirts at three epochs separated by 2 Myr. They compute two virial parameters: the classical one that treats each cloud as an isolated homogeneous sphere, and a full one built from the total gravitational potential energy including everything outside the cloud. The two frequently disagree, with clouds switching between bound and unbound classifications in both directions. The paper concludes that the full environmental potential, not just cloud self-gravity, must be included to determine the dynamical state of molecular clouds.

What carries the argument

The load-bearing object is the full virial parameter $\alpha_{\rm full}=2K/W$, where $K$ is the cloud's kinetic energy about its center of mass and $W=-\sum_i x_i\rho_i\,\partial\Phi/\partial x_i\,\Delta V_i$ is the gravitational energy obtained by summing the full potential gradient over all cloud cells. Unlike the classical $\alpha_{\rm class}=2K/|E_g|$, which approximates the cloud as an isolated homogeneous sphere with $E_g=-(3/5)GM^2/R$, $W$ includes tidal contributions from the galactic disk, spiral structure, and neighboring gas, so its sign and magnitude separate bound ($0>\alpha_{\rm full}\ge-2$), turbulence-unbound ($|\alpha_{\rm full}|>2$), and tidally unbound ($0<\alpha_{\rm full}<2$) states. The paper also uses the volume-averaged tidal tensor, its maximum-to-minimum eigenvalue ratio, and the Toomre parameter to characterize how compressive or extensive the environmental forcing is.

What would settle it

Recompute each cloud's self-gravitational energy from the actual three-dimensional density field (dropping the homogeneous-sphere approximation) and compare $\alpha_{\rm self}=2K/|W_{\rm self}|$ with $\alpha_{\rm full}$; if the disagreement mostly disappears, the tidal attribution is wrong, while if large disagreement remains, the environmental field is confirmed as the cause.

Watch

Extended reading notes

Core claim

The central claim is that the standard virial parameter systematically misclassifies a substantial population of molecular clouds because it ignores the tide from the environment. In the simulation, clouds that look bound under the classical criterion can have positive total gravitational energy and be unbound by tidal stretching, while clouds that look unbound can be gravitationally bound because compressive external tides add to their self-gravity. Quantitatively, the paper defines the full virial parameter $\alpha_{\rm full}=2K/W$, with $W$ computed from the full potential gradient summed over the cloud, and finds significant numbers of clouds in all four dynamical states — turbulence-dominated unbound, tidally unbound, tidally bound, and kinetic-energy unbound — in every region and epoch. The authors present this as agreement with an earlier study that found tides bind and unbind clouds, and they attribute the disagreement with a contrasting study to the smaller, denser clouds of that work, whose deeper potential wells make self-gravity dominate over tides.

Load-bearing premise

The conclusion depends on attributing the gap between the classical and full virial parameters to the surrounding galaxy's tides, yet the classical measure also simplifies each cloud's own gravity to a uniform ball and the paper does not separately compute the cloud's true self-gravity from its actual lumpy density structure, so part of the gap could be internal shape rather than external tides.

Editorial extensions

If this is right

  • Surveys that use the classical virial parameter to label molecular clouds bound or unbound will mislabel a non-negligible fraction, including some clouds that are actually being torn apart and some that are actually held together by environmental compression.
  • Star formation can occur in clouds that appear unbound in the classical analysis, because compressive tidal forces can supplement self-gravity.
  • Cloud dynamical state depends on galactic location: only the central region in this simulation shows fully compressive tidal forces, while intermediate and outer regions contain clouds with at least one extensive tidal axis.
  • Simulations of cloud formation and evolution should include the galactic potential rather than only the cloud's own mass when assessing stability.

Reading between the lines

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

  • A clean way to isolate the tidal effect from internal structure would be to recompute each cloud's true self-gravitational energy from its actual density field and compare the residual gap with the external contribution; this would test whether the classical homogeneous-sphere approximation is partly responsible for the misclassification.
  • The same comparison could be made in higher-resolution simulations of individual clouds embedded in a galactic potential, where the internal density structure is better resolved, to see how the balance between self-gravity and tides shifts with scale.
  • Observationally, one could search for velocity gradients aligned with the extensive eigenvector of the galactic tidal tensor in clouds classified as tidally unbound; such streaming motions would be a direct tidal signature.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. This paper analyzes molecular clouds identified by a dendrogram algorithm in three zoom-in regions (central, intermediate, outskirts) of an AREPO Milky Way-mass galaxy simulation at three epochs. For each cloud, the authors compute a classical virial parameter based on the homogeneous-sphere self-gravity estimate and a 'full' virial parameter using the total gravitational energy W computed from the full potential field. They compare these two diagnostics, examine the sign and magnitude of W, measure tidal-tensor eigenvalue ratios, compute Toomre Q parameters, and investigate the Heyer relation. The paper's central claim is that tidal forces from the environment can both unbind clouds that appear bound by the classical virial parameter and bind clouds that appear unbound, so the full environmental potential must be included in dynamical-state assessments.

Significance. If the central tidal attribution is correct, the paper would strengthen the case that standard observational virial parameters systematically misclassify a significant fraction of molecular clouds, with implications for how cloud boundedness and star formation are inferred. The study has clear strengths: it uses a global galaxy simulation with sub-parsec zoom-in regions, covers three distinct galactic environments and three epochs, computes W directly from the full potential, and offers quantitative population fractions and comparisons to previous work by Ramírez-Galeano et al. and Ganguly et al. No parameters are fitted to the tidal conclusion, and the W>0 branch of the analysis cleanly demonstrates that external tides can unbind some clouds. However, the stronger claim that tides bind apparently unbound clouds is not isolated from the internal density structure of the clouds, because the classical comparison uses a homogeneous-sphere self-gravity estimate while W includes the true self-gravity. The paper's weaker conclusion, that the full potential matters, is supported; the specific tidal attribution is not fully established.

major comments (3)
  1. [Section 2.2, Eqs. (3)-(7), Figure 4, Table 1] The central attribution of the α_class−α_full difference to tides is not isolated. α_class uses the homogeneous-sphere estimate E_g = −3GM²/5R, while W in α_full includes the cloud's actual self-gravitational energy from its nonuniform density plus the external contribution. A cloud with no external tides but a centrally concentrated internal density profile has |W_self| > |E_g|, so α_full < α_class; such a cloud would be placed in the 'bound because of tidal forces' region of Figure 4 (α_class > 2, 0 > α_full > −2) even though W_ext = 0. The fractions in Table 1 and the abstract's claim that tides can 'bind apparently unbound clouds' are therefore under-determined. I recommend computing W_self from the actual density distribution, defining W_ext = W − W_self, and redoing Figure 4 and Table 1 with the two terms separated. The W > 0 branch already demonstrates that external tides unbind some clouds; the bound branch requires this control.
  2. [Section 2.2, Eqs. (3)-(7)] The main virial comparison uses only the bulk kinetic energy K from Eq. (5), which sums cell velocities relative to the cloud center of mass. The observational virial parameter in Eq. (1) and the αtt definition in Eq. (8) include thermal energy, but α_class and α_full do not. If thermal support is non-negligible, some clouds classified as 'bound' by 0 > α_full > −2 may actually be unbound when 2K + 2E_TE > |W|. The authors should quantify E_TE for the classified clouds or explicitly justify that thermal energy is negligible for the populations studied.
  3. [Section 3.3, Eq. (14), Table 2] Equation (14) defines κ = κ_c ≡ 2V_c(R)²/R², which has dimensions of inverse time squared, whereas the epicyclic frequency in Eqs. (11)-(12) has dimensions of inverse time. For a flat rotation curve the correct expression is κ = sqrt(2)V_c(R)/R, or equivalently κ² = 2V_c(R)²/R². As written, the formula would make the Toomre parameters in Table 2 dimensionally inconsistent, and the reported Q values should be checked against the actual implementation.
minor comments (5)
  1. [Section 2.2, text after Eq. (6)] The text says the derivative of 'the cloud's gravitational potential φ_i' was computed, but Eq. (6) and the surrounding sentences require the total gravitational potential Φ from all internal and external mass. If the implementation literally used only the cloud's own potential, W would contain no tidal term; if it used the total potential, the wording should be corrected to avoid ambiguity.
  2. [Section 2.2, Eq. (8)] The notation ETE is used for thermal energy, but the text says 'mass weighted internal energy'; please use 'mass-weighted' consistently and clarify whether this includes only thermal energy or all internal energy.
  3. [Section 3.4] The sentence 'while their magnitude indicates the its strength' contains a typo ('the its'); also, the text should state explicitly that the eigenvalue ratio uses absolute values, as the figure caption does.
  4. [Section 3.5 and final Conclusions bullet] The manuscript attributes the failure to reproduce the Heyer relation to limited numerical resolution within individual clouds. This is plausible but not tested; a resolution study or a direct statement that this is speculative would strengthen the presentation.
  5. [Section 2.2, Eqs. (1) and (3)] The relation between αobs in Eq. (1) and αclass in Eq. (3) should be stated explicitly; as written, it is not clear that the 1D line-width version and the 3D cell-velocity version are equivalent, particularly because Eq. (5) does not include thermal motions.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: the full-virial comparison is a measurement, though the tidal attribution is under-controlled.

full rationale

The paper's derivation chain is not circular. Equations (3)-(7) define two distinct virial parameters from independently computed kinetic and gravitational energies. No parameter is fitted to the α_class–α_full difference, and no conclusion is obtained by renaming an input. The statement that classification by the full potential differs from classification by the classical homogeneous-sphere approximation is a measurement result from the simulation, not a tautology, because the empirical content lies in the fractions of clouds changing state (Table 1, Figures 3-4). The main validity concern is that the paper attributes the α_class–α_full gap to environmental tides without separating the cloud's own gravitational energy computed from its true density distribution from the external tidal term. Since W in Equation (6) contains both the cloud's self-gravity and the environment, the 'bound because of tidal forces' category could be contaminated by internal density structure; the paper never computes W_self alone. This is a missing-control / inference problem rather than a circularity: the tidal attribution is not forced by the equations, and the definitions leave room for internal-structure effects. The self-citations (Li et al. 2020; Li et al. in prep.) are used as data provenance for the simulation, not as unverified theorems carrying the argument, so they are not load-bearing in a circular sense. Accordingly, the circularity score is low.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central computation rests on standard virial theorem machinery from cited prior work; the main assumption ledger is small. The key conceptual weight is carried by the interpretation that departures of alpha_full from alpha_class measure tides, which is an unverified domain assumption rather than a fitted parameter. Dendrogram thresholds and the Heyer-relation density cut are the only freely chosen numbers, and they shape the cloud sample but do not encode the conclusion.

free parameters (4)
  • dendrogram minvalue = 0.5 code density units (~150 cm^-3)
    Threshold density for identifying cloud structures in astrodendro; chosen by hand to produce Milky Way-like clouds. Changes the cloud sample and thus the virial parameter statistics.
  • dendrogram mindelta = 0.5 code density units
    Minimum contrast to surrounding density for dendrogram leaves; chosen by hand as above.
  • dendrogram minnpix = 10 voxels
    Minimum number of 1 pc voxels for a cloud; sets the minimum cloud radius and affects the smallest clouds included.
  • Heyer relation density threshold = 100 cm^-3
    Applied when plotting the Heyer relation in Figure 6; post hoc restriction that changes the slope and scatter of the relation.
assumptions (5)
  • standard math The gravitational virial term for a cloud can be computed as W = -sum x_i rho_i (dPhi/dx_i) dV_i using the full simulation potential.
    Invoked in Eq. (6) from Shu (1992), McKee & Zweibel (1992), and Ramirez-Galeano et al. (2022).
  • domain assumption The surface term in the virial theorem can be neglected for the identified clouds when interpreting alpha_full.
    The paper uses W and alpha_full directly as an energy-based boundedness criterion without evaluating surface pressure terms; this follows the prior literature but is not tested here.
  • domain assumption The AREPO simulation's gravitational potential contains an accurate representation of the total gravitational field, including dark matter, stars, and gas, down to the resolution scale of 0.1-0.2 pc.
    Used in Eq. (6) and Eq. (9); the tidal tensor eigenvalues are derived from this potential.
  • domain assumption Dendrogram structures with the chosen thresholds correspond to physically meaningful molecular clouds.
    Section 2.2 identifies clouds; the parameter values were chosen to match characteristic Milky Way cloud scales but are not derived from a physical definition.
  • domain assumption The classical virial parameter alpha_class with the homogeneous-sphere gravitational energy is the standard observational proxy for boundedness.
    Eq. (3)-(4) and the comparison in Section 3.1 rely on this as the baseline; the approximation is known to be crude for realistic clouds.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The Importance of Tidal Forces in Molecular Cloud Dynamics." pith.science (2026). https://pith.science/paper/YWME6NOU

@misc{pith2026250703788,
  author       = {Pith},
  title        = {Pith review of: The Importance of Tidal Forces in Molecular Cloud Dynamics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YWME6NOU}},
  note         = {Machine review of arXiv:2507.03788}
}
read the original abstract

We investigate the role of tidal forces in molecular cloud formation by examining how apparent boundedness, as diagnosed by the classical virial parameter, relates to the actual gravitational state of clouds subject to tidal forces from their environment. Clouds are identified by a dendrogram algorithm in zoom-in regions taken from a simulation of a Milky Way-mass galaxy with the Voronoi mesh code AREPO that resolves star-forming regions at sub-parsec resolution. To look at a range of environments, we use data from three different regions that evolve differently in the center, near the equivalent of the Solar circle, and the outskirts of the modeled galaxy, at three different times, each spaced 2 Myr apart. We compute the importance of tidal forces on all identified clouds. We then compare the boundedness of clouds including only their internal potentials to boundedness also including the external gravitational potential. This comparison shows that tidal forces can unbind apparently bound clouds and bind apparently unbound clouds. We characterize the cloud population by comparing their virial parameters to their surface densities, finding the ratio of the maximum to the minimum eigenvalues of the tidal tensor, and determining the strength of gravitational instability in each examined region. We find that it is necessary to take the total gravitational potential into account rather than just the internal self-gravity of the clouds to have an accurate understanding of cloud dynamics.

Figures

Figures reproduced from arXiv: 2507.03788 by the authors.

Figure 1
Figure 1. Gas density projection plot of the whole galaxy (left) and the three zoom-in regions: central (upper right), intermediate (mid right), and outskirt (lower right). The three zoom-in figures show gas surface density in the range of 10 − 106M⊙/pc2 . astrodendro (Rosolowsky et al. 2008) on the Carte￾sian grid. There are three parameters that control the astrodendro algorithm for determining when density perturbations sh… view at source ↗
Figure 2
Figure 2. (Top) Size distribution in meters of all clouds across three different regions at three epochs. The minimum cloud radius is derived from our choice of minnpix = 10 voxels. (Bottom) Column densities in three different angles of epoch 540 in the intermediate region only. The three panels for the bottom figure from left to right each correspond to projections along the x-, y-, and z-axis, respectively. Overlapping in g… view at source ↗
Figure 3
Figure 3. Histograms of the virial parameter αfull for clouds that are unbound, with gravitational energy W > 0 (left column), and either bound or unbound with W ≤ 0 (middle column), as well as for the full population (right column). The three different epochs in each panel across the three regions are indicated by the legend. In the middle column, the vertical line at |αfull| = 2 marks the division between bound (left) and u… view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Contour plots of the classical virial parameter αclass against the signed value of the full virial parameter αfull. The pink contours contain clouds with αfull > 0, while the dark blue contours show αfull < 0. Clouds that appear bound with αclass < 2 (below the horizon…
Figure 5
Figure 5. Figure 5: The ratio of the maximum to minimum eigenvalue of the tidal tensor plotted with αtt for the three regions and time epochs (labeled on axes). The red dots indicate fully compressive clouds (all λ < 0), which only occur in the central region, while not fully compressive …
Figure 6
Figure 6. Figure 6: The Heyer et al. (2009) relation for individual clouds at epochs 520, 540, and 560 in the three regions (given in axis labels), using a density threshold to identify clouds of n > 100 cm−3 . The solid red line is a linear regression to the data, while the black dashed …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 9 canonical work pages

  1. [1]

    Ballesteros-Paredes, J., W.Hartmann, L., Vázques-Semandeni, E., Heitsch, F., & Zamora-Avilés, M. A. 2011, Monthly Notices of the Royal Astronomical Society, 411, 65–70, 10.1111/j.1365-2966.2010.17657.x

    arXiv 2011

  2. [2]

    Bertoldi , F., & McKee , C. F. 1992, , 395, 140

    work page 1992

  3. [3]

    D., & Seifried , D

    Ganguly , S., Walch , S., Clarke , S. D., & Seifried , D. 2024, , 10.1093/mnras/stae032

    doi:10.1093/mnras/stae032 2024

  4. [4]

    1965, Monthly Notices of the Royal Astronomical Society, 130, 97–124, https://doi.org/10.1093/mnras/130.2.97

    Goldreich, P., & Lynden-Bell, D. 1965, Monthly Notices of the Royal Astronomical Society, 130, 97–124, https://doi.org/10.1093/mnras/130.2.97

    doi:10.1093/mnras/130.2.97 1965

  5. [5]

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

    doi:10.1146/annurev-astro-082214-122324 2015

  6. [6]

    Heyer, M., Krawczyk, C., Duval, J., & Jackson, J. M. 2009, The Astrophysical Journal, 699, 1092, 10.48550/arXiv.0809.1397

    arXiv doi:10.48550/arxiv.0809.1397 2009

  7. [7]

    C., Mac Low , M.-M., Klessen , R

    Ib \'a \ n ez-Mej \' a , J. C., Mac Low , M.-M., Klessen , R. S., & Baczynski , C. 2016, , 824, 41. 1511.05602

    arXiv 2016

  8. [8]

    V., & Torrey, P

    Li, H., Vogelsberger, M., Marinacci, F., Sales, L. V., & Torrey, P. 2020, Monthly Notices of the Royal Astronomical Society, 499, 5862, 10.1093/mnras/staa3122

    doi:10.1093/mnras/staa3122 2020

Show all 24 references
  1. [9]

    Li , Y., Mac Low , M.-M., & Klessen , R. S. 2006, , 639, 879

    2006

  2. [10]

    V., Vogelsberger, M., Torrey, P., & Springel, V

    Marinacci, F., Sales, L. V., Vogelsberger, M., Torrey, P., & Springel, V. 2019, Monthly Notices of the Royal Astronomical Society, 489, 4233, 10.48550/arXiv.1905.08806

    2019 doi arXiv

  3. [11]

    F., & Zweibel , E

    McKee , C. F., & Zweibel , E. G. 1992, , 399, 551

    1992

  4. [12]

    Y., & Li , H

    Meng , X., Gnedin , O. Y., & Li , H. 2019, , 486, 1574, 10.1093/mnras/stz925

    2019 doi

  5. [13]

    F., Frenk , C

    Navarro , J. F., Frenk , C. S., & White , S. D. M. 1997, , 490, 493, 10.1086/304888

    1997 doi

  6. [14]

    J., Camacho, V., & Zamora-Avilés, M

    Ramírez-Galeano, L., Ballesteros-Paredes, J., Smith, R. J., Camacho, V., & Zamora-Avilés, M. 2022, Monthly Notices of the Royal Astronomical Society, 515, 2822, 10.48550/arXiv.2206.09187

    2022 doi arXiv

  7. [15]

    J., & Christian M., B

    Renaud, F., Gieles, M., Smith, R. J., & Christian M., B. 2011, Monthly Notices of the Royal Astronomical Society, 418, 759–769, 10.1111/j.1365-2966.2011.19531.x

    2011

  8. [16]

    B., & Wiegert , J

    Romeo , A. B., & Wiegert , J. 2011, , 416, 1191

    2011

  9. [17]

    W., Pineda, J

    Rosolowsky, E. W., Pineda, J. E., Kauffmann, J., & Goodman, A. A. 2008, The Astrophysical Journal, 679, 1338, 10.1086/587685

    2008 doi

  10. [18]

    2017, Monthly Notices of the Royal Astronomical Society, 472, 4797–4818, https://doi.org/10.1093/mnras/stx2343

    Seifried, D., Walch, S., Girichidis, P., et al. 2017, Monthly Notices of the Royal Astronomical Society, 472, 4797–4818, https://doi.org/10.1093/mnras/stx2343

    2017 doi

  11. [19]

    Shu, F. H. 1992, The Physics of Astrophysics, Vol. II (Sausalito, California: University Science Books)

    1992

  12. [20]

    H., Adams , F

    Shu , F. H., Adams , F. C., & Lizano , S. 1987, , 25, 23, 10.1146/annurev.aa.25.090187.000323

    1987

  13. [21]

    2010, , 401, 791

    Springel , V. 2010, , 401, 791

    2010

  14. [22]

    1964, , 139, 1217

    Toomre , A. 1964, , 139, 1217

    1964

  15. [23]

    E., et al

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

    2020 doi

  16. [24]

    1994, , 427, 759

    Wang , B., & Silk , J. 1994, , 427, 759

    1994

Pith tools

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