{"id":"c63aba01-2bb6-4a30-9922-56ab95feea97","arxiv_id":"2501.03027","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A two-component model combining self-gravity and galactic shear reproduces the steep velocity dispersion-size relation of molecular clouds and its radial variation across the Milky Way.","lead":"Molecular clouds show a puzzlingly steep relation between their size and internal motion, and this paper proposes that the cause is a shift from self-gravity at small scales to galactic shear at large scales. The model reproduces the steep slope and predicts how the relation changes with distance from the Milky Way's center, so a generalist can see one way that cloud-scale and galaxy-scale physics connect.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The solar-vicinity calibration fits the model to a synthetic velocity-dispersion target drawn from the same steep relation the model claims to explain; the independent samples then fix f=1 rather than measure it, so the paper does not yet establish the gravity-shear transition.","rationale":"The paper presents a physically plausible two-component model and shows some encouraging agreement in the independent predictions for Miville-Deschenes et al. (2017) and Sun et al. (2024). However, the load-bearing parameter f is only calibrated on a synthetic target that already encodes the steep slope being explained. The independent samples are used with f fixed, so they do not provide a genuine test of the mechanism. A two-parameter fit to those samples would settle the issue, as described in the concrete test. This concern matches the Reader's weakest-assumption and rationale, so the verdict remains REJECT.","tokens_in":13580,"tokens_out":8300,"duration_ms":80464,"concrete_test":"Refit Eq. (1) to the Sun et al. (2024) resolved clouds (R>2 pc, Rgal<20 kpc) and the Miville-Deschenes et al. (2017) sample with both A and f free, using the same shear timescale (Eq. 4) and XCO calibration as the paper, and report best-fit f with bootstrap uncertainties. If f is consistent with 0.96 and the equality scale in Eq. (8) is ~100 pc in these independent samples, the solar calibration is transportable and the central claim is supported; if f is unconstrained or differs significantly from 1, the paper's fixed-f application to Miville and Sun is invalid and the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 3.2.2, the model parameters A and f are determined by assigning every Xie et al. (2024) cloud a velocity dispersion from the Zhou et al. (2022) power-law relation via Eq. (6), sigma_v(1D)=0.52(R/pc)^0.67 km/s, plus 20% Gaussian noise. The fit target is therefore the very steep relation the paper claims to explain. Because Xie et al. clouds obey M proportional to R^1.94 (Eq. 11), the gravity term in Eq. (1) scales as R^0.47 and the shear term as R; over a limited size range a linear combination can approximate an R^0.67 power law, so the good fit and the inferred f ~ 0.96 are largely built into the input, not discovered from it. The ~100 pc transition scale is then a derived property of the assumed Zhou slope and the mass-size relation. The applications to Miville-Deschenes et al. (2017) and Sun et al. (2024) fit only the normalization A while fixing f=1 based on the solar fit; hence the key parameter that controls the gravity-shear transition is never tested against independent data. Since Sec. 3.1 states that A and f are expected to vary with tracer and cloud definition, propagating a single f without uncertainty or a free fit is an unjustified leap. The virial-parameter trend (Eq. 10) is also a direct consequence of the model definition (Eq. 9), not an independent check.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes that the observed steep velocity dispersion-size relation (β > 0.6) in molecular clouds above a few parsecs results from a gradual transition between self-gravity-dominated and galactic-shear-dominated regimes. The model is σ_v,total = A[(GM/R)^{1/2} + f(R/t_shear)]. The authors calibrate A and f using the Zhou et al. (2022) YSO-association scaling relation imposed on the Xie et al. (2024) dust clouds, infer a transition scale of about 100 pc, and then apply the model with f = 1 to the Miville-Deschênes et al. (2017) and Sun et al. (2024) samples to predict the radial variation of the scaling relation. The paper concludes that the steep slope and its environmental variation are explained by the interplay of internal gravity and external shear.","tokens_in":13941,"tokens_out":7172,"duration_ms":132200,"significance":"If validated, the model would provide a simple physical explanation for a currently puzzling observational result and would make concrete, testable predictions for resolved cloud surveys and simulations. The analytical formula is transparent, the manuscript uses published catalogues, and the radial predictions are specific and falsifiable. However, the current significance is limited because the solar-vicinity calibration uses a synthetic velocity-dispersion target, and the external-sample comparisons fix the crucial parameter f rather than measure it. With an independent calibration or a held-out test, the paper could be a valuable contribution to the interpretation of molecular cloud scaling relations.","major_comments":[{"comment":"The calibration target is synthetic, not measured. The velocity dispersion of every Xie et al. (2024) cloud is generated from the Zhou et al. (2022) relation through Eq. (6), σ_v(1D)=0.52(R/pc)^0.67 km/s, with 20% Gaussian noise, rather than being a directly measured quantity for those clouds. Fitting Eq. (1) to this target therefore cannot validate the model: the fitted slope (Eq. 7) and f ≈ 0.96 largely recover the input relation, and the ~100 pc transition scale in Sec. 3.2.3 is a derived property of the assumed Zhou slope and the Xie et al. mass-size relation (Eq. 11). The central claim that the model explains the steep observed relation is not supported by this exercise. Please re-fit using directly measured cloud velocity dispersions, or explicitly present the solar-vicinity result as a demonstration of the model's functional form rather than as an observational test.","section":"Sec. 3.2.2, Eq. (6)"},{"comment":"The external applications do not test the key parameter of the model. For both Miville-Deschênes et al. (2017) and Sun et al. (2024), f is fixed to 1 based on the solar calibration, and A is fitted on the same dataset whose relation is then compared with the prediction. Since Eq. (1) is linear in A, the normalization agreement is guaranteed by construction. The gravity-shear transition, controlled by f, is therefore never measured independently. The authors should fit A and f freely on one sample and predict a held-out sample, or at least test whether f=1 is consistent with the Miville-Deschênes and Sun data, and propagate uncertainties into the predicted slopes.","section":"Sec. 3.3 and Appendix A"},{"comment":"The size calibration assumes a fixed 45° angle between the cloud major axis and the line of sight. The projected size enters all subsequent fits, so the inferred A, f, the ~100 pc transition scale, and the virial-parameter slope all depend on this assumption, but no sensitivity study or uncertainty is given. Please vary the projection angle over a plausible range or marginalize over it.","section":"Sec. 3.2.1, Eq. (5)"},{"comment":"The virial-parameter trend is not an independent check. As the authors note in Eq. (12), the positive slope β3 = 2β2 - β1 + 1 is forced by the fitted mass-size relation (Eq. 11) and velocity dispersion-size relation (Eq. 7) once the model form is adopted. The text should state clearly that this is a derived consistency relation, not empirical support for the gravity-shear mechanism.","section":"Sec. 3.2.4, Eq. (12)"},{"comment":"The assumption that f remains equal to 1 at all Galactocentric distances is not justified. Section 3.1 explicitly states that A and f are expected to vary with cloud definition and observational tracer, and the shear efficiency relative to gravity could plausibly vary with radius, yet the radial predictions fix f=1 and only vary A and the shear timescale. Without a physical argument or a fit of f in the radial bins, the predicted radial variation is conditional on an untested assumption.","section":"Sec. 3.3, Eq. (17); Appendix A, Eq. (A4)"}],"minor_comments":[{"comment":"The conversion σ_v,2D = √2 σ_v,1D assumes isotropic velocity dispersions; please state this assumption and discuss its effect on the calibration.","section":"Sec. 3.2.1"},{"comment":"The caption contains 'km s□1' in several places instead of 'km s^{-1}', and Eq. (A1) quotes σ_v(Sun)=0.19(R/pc)^0.78 while the text mentions β=0.65±0.004 for the whole resolved sample without clarifying which sample that slope refers to; please correct and clarify.","section":"Fig. 5 caption, Eq. (A1)"},{"comment":"The Facilities line contains 'FL WO:2MASS, CTIO:2MASS', which appears to be a formatting artifact; please correct.","section":"Facilities line"},{"comment":"Best-fit values A and f are quoted without uncertainties; please provide error bars and, for the transition scale, a confidence interval.","section":"Secs. 3.2.2, 3.3, Appendix A"},{"comment":"The identification of the ~100 pc transition scale with the molecular gas disk scale height is only qualitative; please quantify the comparison (e.g., with a reference value and its radial variation) or label it as a suggestion.","section":"Sec. 3.2.3"},{"comment":"The manuscript uses 'explain' for what is currently a reproduction of an assumed relation; please use more cautious wording such as 'is consistent with' or 'can reproduce' unless an independent test is added.","section":"Abstract and Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The main issue is the circular calibration: the solar-vicinity fit uses a synthetic velocity-dispersion target, and the external samples fix f rather than measure it. The model is simple and potentially useful, but the current manuscript does not provide an independent test. If the authors can obtain direct velocity-dispersion measurements for the Xie et al. clouds, or fit f freely on one of the large samples and test on a held-out sample, the paper could become publishable. I recommend major revision rather than rejection because the required additional analysis is plausible within the scope of the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read the paper. The central calibration is circular, and that kills the main claim as presented. In Sec. 3.2.2 they take the Xie et al. clouds, compute each cloud's projected radius with an assumed 45-degree angle, and then assign each cloud a velocity dispersion from the Zhou et al. relation via Eq. (6), plus 20% noise. That synthetic sigma_v is exactly the power law they say they are explaining. Fitting Eq. (1) to it and getting f ~ 0.96 tells you about the relation you put in, not about the physics. The derived 100 pc transition scale is then a consequence of the assumed Zhou slope and the mass-size relation, not an independent measurement.\n\nThe later applications do not rescue it. For Miville-Deschenes et al. and Sun et al. they fix f = 1 (based on the solar fit) and fit only the normalization A. So the key parameter that sets the gravity-shear balance is never actually measured on independent data. The close agreement with the Sun et al. radial slopes (within ~0.07) is a point in the model's favor, but a weak one: A is re-fit per sample, f is fixed, so the slope variation is mostly coming from the density and shear inputs, not from a free f.\n\nCredit where it's due: the paper lays out the observational puzzle well, the two-component additive formula is a reasonable ansatz, and the idea that the transition scale might be tied to the disk scale height is worth thinking about. The authors are honest about the synthetic calibration in the text; they do not hide it.\n\nBut the virial-parameter trend (Eq. 10) is not a check; it is built into the model definition and the fitted sigma_v-R and M-R relations. The 45-degree projection angle (Eq. 5) is an arbitrary assumption with no sensitivity analysis. No error bars are given for A or f. And the paper's own statement that A and f are expected to vary with tracer and cloud definition sits awkwardly with fixing f = 1 in the later sections.\n\nNet: the mechanism is plausible and worth testing, but this paper does not establish it. The fix is straightforward: take a sample with measured sizes, masses, and velocity dispersions, fit both A and f with uncertainties, and check whether f ~ 1 holds out of sample. As presented, I would recommend rejection, not because the idea is bad but because the evidence is circular.","headline":"The steep sigma_v-R slope is explained by fitting a model to a synthetic target drawn from the same relation, so the gravity-shear transition is not yet established.","tokens_in":14439,"tokens_out":4377,"would_cite":false,"duration_ms":36243,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A two-mechanism sum explains the steep cloud velocity-size law","keywords":["molecular clouds","velocity dispersion-size relation","galactic shear","self-gravity","Larson's relations","interstellar dynamics","Galactic rotation","Milky Way disk"],"falsifier":"Measure, for a sample of molecular clouds spanning 10–300 pc in the solar vicinity, the internal velocity dispersion and the velocity gradient expected from Galactic rotation across each cloud's projected extent; the model predicts that the shear term rises linearly with size and becomes comparable to the self-gravity term near 100 pc. A cleaner test is the predicted shape of the $\\sigma_{v}$–$R$ relation itself: it should curve from slope about 0.5 at small sizes to about 1 at large sizes, with the transition at roughly the local disk scale height; if the relation is a single power law with no detectable steepening over 10–200 pc, or if the transition scale differs by more than a factor of two from the disk scale height, the gravity–shear transition described here would be ruled out.","tokens_in":13352,"feed_emoji":"🌌","tokens_out":6329,"duration_ms":52604,"temperature":0.7,"pith_summary":"Molecular clouds larger than a few parsecs show a velocity dispersion–size relation steeper than the classic Larson law with slope about 0.5, and no single mechanism explains the observed slopes above 0.6. This paper argues that the steep slope arises from a gradual transition between two regimes: small clouds whose internal motions are set by self-gravity, and large clouds whose motions are set by the shear of Galactic rotation. The two are combined in one analytical formula, $\\sigma_{v,\\rm total}=A[(GM/R)^{1/2}+f(R/t_{\\rm shear})]$, and the transition scale comes out near 100 pc, about the scale height of the Milky Way's molecular gas disk. If the model is right, the steep slope and its variation with Galactocentric distance are not separate puzzles but the same interplay of internal gravity and external shear. The paper confronts the formula with observations in the solar vicinity and across the Galactic disk, and reports agreement with the measured slopes.","feed_headline":"Steep cloud velocity-size slopes come from a gravity-to-shear switch","feed_subtitle":"One formula with internal gravity and galactic shear reproduces slopes above 0.6 and their drift across the Milky Way.","key_machinery":"The carrying object is Equation (1), a two-component velocity dispersion formula $\\sigma_{v,\\rm total}=A[(GM/R)^{1/2}+f(R/t_{\\rm shear})]$, where $M$ and $R$ are the cloud mass and size, $t_{\\rm shear}=\\kappa^{-1}=(2A_{\\rm Oort})^{-1}$ is the shear timescale from Galactic rotation, and $A$ and $f$ are fitted parameters that set the normalization and the relative efficiency of the two channels. The self-gravity term scales as $R^{1/2}$ under the Larson mass–size relation and the shear term scales as $R^{1}$, so their relative weight changes with cloud size, making the slope of the combined relation interpolate continuously between 0.5 and 1. The paper uses the ratio $\\lambda_{1}\\sigma_{v,g}/\\lambda_{2}\\sigma_{v,\\rm shear}=t_{\\rm shear}(GM/R^{3})^{1/2}/f$ to locate the transition scale and to map gravity-dominated versus shear-dominated regions of the Galaxy; this ratio is the diagnostic that turns the formula from a curve fit into a physical classification of cloud dynamics.","core_discovery":"The paper's central claim is that the steep velocity dispersion–size relation $\\sigma_{v}\\sim R^{\\beta}$ with $\\beta\\sim0.6$–$0.8$, observed for molecular clouds above several parsecs, is produced by the combined action of self-gravity and Galactic shear. Gravity alone gives $\\sigma_{v}\\sim R^{1/2}$ and shear alone gives $\\sigma_{v}\\sim R^{1}$; neither reproduces the observed slope. The authors posit a two-component sum $\\sigma_{v,\\rm total}=A[(GM/R)^{1/2}+f(R/t_{\\rm shear})]$ with $\\sigma_{v,\\rm g}=(GM/R)^{1/2}$ and $\\sigma_{v,\\rm shear}=R/t_{\\rm shear}$, where $t_{\\rm shear}$ is the local shear timescale. Under the Larson mass–size relation $M\\sim R^{2}$, each term retains its own slope, so the combined relation spans the observed intermediate slopes through a gradual transition. Fitting the formula to solar-vicinity data yields $A=1.94$, $f=0.96$, and a transition scale near 100 pc at which gravity and shear contribute equally; small clouds are gravity-dominated and roughly virialized, while large clouds are shear-dominated and supervirial. Applied to Galactic disk samples at different Galactocentric distances, the same formula predicts the observed normalization and slope changes, tracing them to variations in cloud density structure and shear rate.","pith_inferences":["The same two-component formula, with $t_{\\rm shear}$ set by the local rotation curve, could be applied to other galaxies: their disk scale heights and rotation curve shapes would set their own gravity–shear transition scales, turning the transition scale into a diagnostic of galactic environment.","Because the gravity term depends on cloud mass, the model implies that the observed slope and normalization depend on how clouds are defined and on the tracer used; comparing different tracers with the same fitted parameters would test whether $A$ and $f$ absorb real physics or just sample-selection effects.","A direct test would measure velocity gradients across large clouds that are attributable to Galactic rotation and compare them with the internal velocity dispersion; if shear contributions are absent in clouds near 100 pc, the claimed transition scale would fail.","The predicted rise of virial parameter with cloud size, if confirmed with independent samples, would contradict the frequently reported declining virial parameter–size trend, suggesting that boundary definitions rather than physics drive much of that discrepancy."],"forward_implications":["Small clouds should be roughly in virial equilibrium while large clouds are supervirial, with a virial parameter that rises with size as $\\alpha_{\\rm vir}\\sim R^{0.5}$ for the calibrated solar-vicinity sample.","The observed normalization of the $\\sigma_{v}$–$R$ relation should decline with Galactocentric distance as both the shear rate and cloud surface density fall; the model predicts $\\sigma_{v}=0.52\\,R^{0.61}$ at 5–6 kpc and $0.24\\,R^{0.67}$ at 10–11 kpc for the Miville-Deschênes et al. (2017) sample.","The slope of the relation should steepen slightly with Galactocentric distance in the outer Galaxy, from 0.81 to 0.84 across the Sun et al. (2024) radial bins, driven by the declining shear rate.","Below a few parsecs the relation is expected to flatten or decouple from cloud size, consistent with observed breaks, because the shear–gravity combination no longer controls the dynamics on those scales."],"supporting_citations":[{"why":"Supplies the steep two-dimensional velocity dispersion–size relation for YSO associations in the solar vicinity that anchors the calibration of A and f.","marker":"Zhou et al. (2022)"},{"why":"Provides the cloud masses and sizes in the solar vicinity onto which the Zhou et al. relation is imposed to build the fit target.","marker":"Xie et al. (2024)"},{"why":"Galactic-disk cloud sample with size, mass, and velocity dispersion used to predict the radial variation of the relation.","marker":"Miville-Deschênes et al. (2017)"},{"why":"Outer-disk cloud sample used to test the predicted slopes at different Galactocentric distances.","marker":"Sun et al. (2024)"},{"why":"Establishes the velocity dispersion–size and mass–size relations that the two-component model extends.","marker":"Larson (1981)"},{"why":"Provides the Oort constant A_Oort = 15.1 km/s/kpc that sets the solar shear timescale.","marker":"Li et al. (2019)"},{"why":"Supplies the Galactic rotation curve used to compute shear timescales at different Galactocentric distances.","marker":"Mróz et al. (2019)"},{"why":"Provides the radial X_CO calibration used to convert the Miville-Deschênes et al. cloud masses.","marker":"Lada & Dame (2020)"}],"fun_headline_variants":["Cloud turbulence steep slopes: gravity plus shear formula fits","Why molecular cloud size-velocity slopes exceed 0.5: gravity to shear","Gravity and shear together set steep cloud size-velocity relation","Two forces, one formula: steep cloud scaling explained","Shear not gravity alone steepens cloud velocity-size slope"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calibration assumes that the velocity dispersion–size relation measured for YSO associations by Zhou et al. (2022) can be assigned, point by point, to the dust clouds of Xie et al. (2024) as their individual velocity dispersions, making the fit target synthetic rather than directly measured; if this cross-sample equivalence or the assumed 45-degree projection fails, the fitted parameters, the 100 pc transition scale, and the radial predictions lose their anchor.","fun_headline_variants_meta":{"raw":{"variants":["Cloud turbulence steep slopes: gravity plus shear formula fits","Why molecular cloud size-velocity slopes exceed 0.5: gravity to shear","Gravity and shear together set steep cloud size-velocity relation","Two forces, one formula: steep cloud scaling explained","Shear not gravity alone steepens cloud velocity-size slope"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000754,"raw_usage":{"total_tokens":3469,"prompt_tokens":1177,"completion_tokens":2292,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":793,"completion_tokens_details":{"reasoning_tokens":2207}},"tokens_in":793,"tokens_out":2292,"duration_ms":15178,"temperature":1.0,"reasoning_tokens":2207,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:58:40.654366+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, for a sample of molecular clouds spanning 10–300 pc in the solar vicinity, the internal velocity dispersion and the velocity gradient expected from Galactic rotation across each cloud's projected extent; the model predicts that the shear term rises linearly with size and becomes comparable to the self-gravity term near 100 pc. A cleaner test is the predicted shape of the $\\sigma_{v}$–$R$ relation itself: it should curve from slope about 0.5 at small sizes to about 1 at large sizes, with the transition at roughly the local disk scale height; if the relation is a single power law with no detectable steepening over 10–200 pc, or if the transition scale differs by more than a factor of two from the disk scale height, the gravity–shear transition described here would be ruled out.","supporting_citations":[],"review_version":1}