{"id":"36e6cec0-7fc0-4631-9a5f-cc3480b93e9f","arxiv_id":"2507.03788","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In AREPO zoom-in simulations of a Milky Way-like galaxy, accounting for the full gravitational potential changes the bound state of a substantial fraction of molecular clouds relative to the self-gravity-only virial parameter.","lead":"This paper tests whether molecular clouds are bound by their own gravity or by the gravity of the surrounding galaxy, using a high-resolution simulation of a Milky Way-like galaxy. It finds that the standard virial parameter, which ignores the galaxy's tidal pull, often labels clouds as unbound when the full gravitational field says they are bound.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No control separating W_self from W_ext: α_class vs α_full conflates internal cloud structure with environmental tides, so the tidal attribution is unproven.","rationale":"The reader's weakest assumption correctly identifies the critical gap: the difference between α_class and α_full is attributed to tides without a self-gravity-only baseline. I considered whether a more severe issue exists—e.g., whether W in Eq. 6 is computed with coordinates relative to the cloud's center of mass; the text only states 'the position vector x_i' and thus leaves this ambiguous. If absolute coordinates were used, W would be origin-dependent and the entire analysis invalid. However, the references cited (Shu 1992; McKee & Zweibel 1992) define W in the center-of-mass frame, and the authors explicitly specify center-of-mass coordinates for the kinetic energy in Eq. 5, so it is plausible they did the same for W. Without access to the code I cannot confirm, but this ambiguity is secondary because even with proper coordinates the self-gravity contamination remains. A second candidate—the tidal tensor sign convention (trace vs traceless)—does not affect the virial energy computation. Thus the single most load-bearing concern is the missing W_self. The paper's broad conclusion that the full potential matters is supported by the W>0 population (α_class<2, α_full>0), which cannot arise from self-gravity alone, but the specific and novel claim that tides bind clouds is not demonstrated. A W_self control directly resolves the attribution and should be reported with the W_ext/W_self ratio.","tokens_in":10405,"tokens_out":9732,"duration_ms":106478,"concrete_test":"Compute for every cloud the self-gravitational energy W_self using the same Eq. 6 but with the potential generated only by cells belonging to that cloud, keeping the identical grid and finite differencing. Form W_ext = W − W_self and recalculate α_self = 2K/W_self. If, among clouds with α_class > 2 and 0 > α_full > −2, the majority have α_self < 2, then their apparent binding is internal structure, not tides, and the tidal-binding claims in §3.1 and the abstract would need revision. Also report the distribution of W_ext/W_self to quantify the tidal contribution directly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that tidal forces bind and unbind clouds rests on the gap between α_class (Eq. 3, with E_g = −(3/5)GM^2/R from a homogeneous sphere, Eq. 4) and α_full (Eq. 7, computed from W in Eq. 6, which is the total gravitational energy from the full potential including both the cloud's own inhomogeneous mass distribution and all external matter). In Figure 4, clouds with α_class > 2 but 0 > α_full > −2 are labeled 'bound because of tidal forces,' but W contains W_self from the actual (nonuniform) cloud density. A centrally concentrated cloud can have |W_self| > (3/5)GM^2/R, driving α_full below 2 even with zero external tides. Since W_self is never computed and compared to α_class, the sign and magnitude of W_ext are not isolated. The fractions in Table 1 and the abstract's statement that 'tidal forces can ... bind apparently unbound clouds' are therefore under-determined; they could reflect the internal density structure of the clouds rather than the environment. The weaker conclusion that the full potential matters is supported by the W>0 branch (which requires external tides), but the paper's specific tidal attribution, its novelty over the classical virial picture, is not.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":10643,"tokens_out":13084,"duration_ms":163795,"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":[{"comment":"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.","section":"Section 2.2, Eqs. (3)-(7), Figure 4, Table 1"},{"comment":"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.","section":"Section 2.2, Eqs. (3)-(7)"},{"comment":"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.","section":"Section 3.3, Eq. (14), Table 2"}],"minor_comments":[{"comment":"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.","section":"Section 2.2, text after Eq. (6)"},{"comment":"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.","section":"Section 2.2, Eq. (8)"},{"comment":"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.","section":"Section 3.4"},{"comment":"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.","section":"Section 3.5 and final Conclusions bullet"},{"comment":"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.","section":"Section 2.2, Eqs. (1) and (3)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for this journal and the required control is feasible with the existing simulation data: computing W_self from the actual density field and isolating W_ext would directly test the tidal attribution. I do not recommend rejection, because the weaker claim that the full potential matters is supported and the data set is valuable. The authors should also clarify the relation to the simulation paper by Li et al. (in prep.) and consider whether the thermal-energy omission changes the bound/unbound classification."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper has a new and useful simulation result, but the headline claim about tidal forces is not yet proven. The authors compare the classical virial parameter (which assumes a homogeneous sphere) with the full virial parameter computed from the actual potential of the cloud plus everything around it. They find that many clouds flip between bound and unbound, and they attribute the flips to tides. That attribution is under-supported because the full potential includes the cloud's own nonuniform internal mass distribution, not just external tides. A centrally concentrated cloud would have stronger self-gravity than the homogeneous estimate, which could cause the same flips with zero external tides. The paper never isolates W_self from W_ext.\n\nWhat's genuinely new is applying this analysis to a global AREPO zoom-in run with three environments at three epochs, and showing that only the central region has fully compressive tides. The W > 0 branch is a solid result: if the total potential energy is positive, external tides must be net extensive, so that part does isolate an environmental effect. The paper also honestly reports that its clouds do not reproduce the Heyer relation, likely due to resolution, rather than over-claiming.\n\nThe main soft spot is the missing control. Clouds with α_class > 2 but 0 > α_full > -2 are labeled 'bound because of tidal forces,' but W_self is never computed from the actual density distribution. The fractions in Table 1 and the abstract's statement that tidal forces can bind apparently unbound clouds are therefore not cleanly attributed. This is not a minor point—it is the core of the paper. The fix is straightforward: compute W_self for each cloud from its own density field, compare it with the homogeneous-sphere estimate, and either subtract it from W to obtain W_ext or present the distribution of central-concentration corrections. I suspect the qualitative result will survive for at least some clouds, but the quantitative fractions will shift.\n\nMinor concerns: the epicyclic approximation for κ is crude but adequate here; the discussion of the Ramirez-Galeano vs Ganguly discrepancy is reasonable, though necessarily speculative about which simulation property drives the difference.\n\nWho this is for: people working on cloud virial parameters, galaxy-scale star formation simulations, and observers trying to interpret the observational virial parameter. It deserves a serious referee—the full-potential comparison is worth publishing—but the referee should require the W_self/W_ext separation before the tidal conclusion goes in the abstract. As is, it is a strong conference talk with a title that runs ahead of the evidence.","headline":"The full-potential result is real, but the tidal-specific claim needs a control separating internal structure from external tides.","tokens_in":11187,"tokens_out":4896,"would_cite":true,"duration_ms":49781,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["molecular clouds","tidal forces","virial parameter","virial theorem","galactic environment","star formation","numerical simulations","tidal tensor"],"falsifier":"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.","tokens_in":10200,"feed_emoji":"🌌","tokens_out":11993,"duration_ms":124998,"temperature":0.7,"pith_summary":"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.","feed_headline":"Tidal forces can bind or unbind molecular clouds","feed_subtitle":"Simulations show the standard virial check mislabels many clouds; the galaxy's pull must be counted.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Defines the classical observational virial parameter that the paper uses as the baseline for apparent boundedness.","marker":"Bertoldi & McKee 1992"},{"why":"Provides the time-dependent virial theorem underlying both the classical and full virial parameters.","marker":"McKee & Zweibel 1992"},{"why":"Supplies the formula for total gravitational energy W used in Equation (6).","marker":"Shu 1992"},{"why":"Introduced the full virial parameter including the environmental potential and reported that tides bind and unbind clouds, the result this paper reproduces.","marker":"Ramírez-Galeano et al. 2022"},{"why":"The contrasting study finding tides only deform existing clouds; the paper's explanation of the disagreement rests on their denser, smaller clouds.","marker":"Ganguly et al. 2024"},{"why":"Provided the dendrogram algorithm used to identify clouds in the density field.","marker":"Rosolowsky et al. 2008"},{"why":"The moving-mesh hydrodynamics code used to run the galaxy and zoom-in simulations.","marker":"Springel 2010"},{"why":"The parent Milky Way-mass galaxy simulation that the zoom-in technique refines to sub-parsec resolution.","marker":"Li et al. 2020"}],"fun_headline_variants":["Tides bind and unbind molecular clouds","Virial check mislabels clouds, tides matter","Galactic tides flip cloud fate","Clouds bound? Tides decide","Molecular clouds: tides tip the balance"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Tides bind and unbind molecular clouds","Virial check mislabels clouds, tides matter","Galactic tides flip cloud fate","Clouds bound? Tides decide","Molecular clouds: tides tip the balance"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000176,"raw_usage":{"total_tokens":1304,"prompt_tokens":975,"completion_tokens":329,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":266}},"tokens_in":591,"tokens_out":329,"duration_ms":3684,"temperature":1.0,"reasoning_tokens":266,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:02:54.061174+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the classical observational virial parameter that the paper uses as the baseline for apparent boundedness."},{"cited_title":"F., & Zweibel , E","cited_arxiv_id":null,"evidence_quote":"Provides the time-dependent virial theorem underlying both the classical and full virial parameters."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the formula for total gravitational energy W used in Equation (6)."},{"cited_title":"Why most molecular clouds are gravitationally dominated","cited_arxiv_id":"2206.09187","evidence_quote":"Introduced the full virial parameter including the environmental potential and reported that tides bind and unbind clouds, the result this paper reproduces."}],"review_version":1}