{"id":"8d2e69e2-bc19-43b2-b104-52932b1341bd","arxiv_id":"2501.08065","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Hydrodynamically coupled, non-motile oscillating dumbbells show a density-dependent transition to collective motion, hydrodynamic phase separation, and phase-locking of shape oscillations.","lead":"A simulation study shows that dumbbell-shaped cells that only wobble in place can start moving together when crowded, purely through forces they exert on the surrounding water. The result suggests a route to collective cell migration that does not require individual cells to be motile or polarized.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Missing control with hydrodynamic coupling switched off leaves open that steric repulsion, not hydrodynamics, drives the observed phase-locking at high density.","rationale":"In good faith, the paper presents a coherent minimal model: a bead-spring dumbbell with myosin kinetics exhibits a Hopf limit cycle; single dumbbells are non-motile by the scallop theorem; hydrodynamic coupling between such dumbbells can produce pair-level drift (Eqs. 15-16) and, in the full 1D periodic simulations, density-dependent collective motion, phase separation, and synchronization. The qualitative phenomena are internally consistent, and the connection to observed shape oscillations of suspended cells gives biological motivation. The pair-level analytical argument is stated without derivation, and the vs ~ Φ^4 scaling is fit from few points, but these are secondary to the main simulation evidence. The most load-bearing gap is the absence of a hydrodynamic-off control. The model includes a short-range repulsive force that is active precisely in the high-density regime where the claimed phase-locking and clustering appear; this force mechanically couples the oscillators and could produce synchronization independently of hydrodynamics. Because the title and central claim are specifically that hydrodynamics drives phase-locking and collective behavior, failing to isolate hydrodynamics from steric coupling leaves the causality unresolved. The reader's weakest assumption concerned 1D geometry and unreported parameters; my concern is related but distinct, focusing on the confound between hydrodynamic and steric inter-dumbbell couplings. A simple G=0 control would settle this. The verdict should remain CONDITIONAL, with the added condition that the authors either run the no-hydrodynamics control and show the phenomena vanish, or explicitly characterize the repulsive interactions and reframe the title/claims accordingly. This is not a rejection because the observations may well be hydrodynamic in origin, and the missing control is a standard, easy-to-run test.","tokens_in":13849,"tokens_out":11429,"duration_ms":123791,"concrete_test":"Rerun the pair simulations behind Fig. 8 and the N=64 suspensions at Φ=0.7 and Φ=0.9 with the hydrodynamic coupling set to zero (G_{αβ}=0 in Eq. 7), keeping the repulsion of Eq. (5), myosin kinetics, and all parameters otherwise identical. Compute mutual information I(c1;c2), correlation R, kymographs of bound myosin, cluster counts, and center-of-mass MSD. If phase-locking and cluster formation still appear, the 'hydrodynamics-driven' phase-locking claim is unsupported; if they disappear, hydrodynamics is the essential driver. As a sanity check, the G=0 center-of-mass MSD should remain identically zero (no net external force), confirming that any observed drift requires hydrodynamics.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that hydrodynamic interactions are the driver of all three phenomena. However, Eq. (7) also contains a second inter-dumbbell coupling: the soft repulsion of Eq. (5), which turns on exactly in the high-density regime where phase-locking, clustering, and collective motion emerge. At Φ ≥ 0.4, beads of neighboring dumbbells are frequently within 2a, so repulsive forces directly couple length oscillations of different dumbbells through dl_i/dt in Eq. (2), providing a purely mechanical (non-hydrodynamic) synchronization path. The paper never runs a control with G_{αβ}=0 in Eq. (7); without it, the claim that hydrodynamic interactions drive the phase-locking transition is not established. If phase-locking persists without hydrodynamics, the title and abstract overstate the mechanism: collective drift would still need hydrodynamics, but synchronization and the resulting cluster organization could be steric in origin. This is load-bearing because the paper's novelty is precisely the hydrodynamic origin of synchronization and phase separation.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript examines a one-dimensional suspension of bead-spring dumbbells whose internal myosin-like activity produces limit-cycle oscillations in dumbbell length and bound-myosin concentration. In isolation a dumbbell is non-motile because its shape changes are reciprocal; the authors couple the beads to a Stokes fluid via Oseen hydrodynamic interactions and include a soft repulsion between beads. From deterministic simulations of 64 dumbbells with periodic boundary conditions, they report three density-dependent phenomena: a transition from diffusive to ballistic center-of-mass motion near Φ≈0.4 with collective velocity scaling as Φ^4, hydrodynamic phase separation into clusters for intermediate densities, and phase-locking of the myosin and length oscillations at high density. A pair-level analytical argument (Eqs. (11)–(16)) is presented to show that phase differences between dumbbells can produce net translation and phase-dependent attraction/repulsion, and pair simulations quantify phase-locking through mutual information and correlation.","tokens_in":14093,"tokens_out":8026,"duration_ms":79167,"significance":"The reported phenomena, if fully supported, would be a valuable conceptual contribution: they demonstrate a route to collective motility in suspensions of agents that have no intrinsic polarity, no self-propulsion, and no adhesive interactions, with a density threshold and cooperative cluster formation arising from hydrodynamic coupling alone. The paper benefits from systematic density sweeps, complementary order parameters (mutual information and correlation), cluster counting, and a candid discussion of the model's simplifications. However, the manuscript does not yet provide a derivation of the central pair formulas, a parameter table, integration details, code, or the necessary control simulation with hydrodynamics switched off. These omissions leave the quantitative threshold and the attribution of the mechanism to hydrodynamics provisional. I still regard the qualitative scenario as credible and worth publishing after a substantive revision.","major_comments":[{"comment":"The central attribution of phase-locking and clustering to hydrodynamic interactions is not established because no control simulation is reported with the hydrodynamic coupling switched off (G_αβ=0 in Eq. (7)). At Φ≥0.4 the soft repulsion of Eq. (5) directly couples beads of neighboring dumbbells and can influence dl_i/dt in Eq. (2), providing a steric synchronization mechanism that may reproduce the high-density MI/correlation and cluster signals. Please add simulations without the Oseen coupling (and, if feasible, without repulsion) and show that the phase-locking transition, cluster statistics, and diffusive-to-ballistic transition are absent or strongly modified. Without such a control, Eqs. (2)–(8) do not justify the phrase 'hydrodynamics-driven' in the title and abstract.","section":"§2.3, Eq. (7); §3.3"},{"comment":"The analytical pair result is stated without derivation, and the derivation assumes ω1=ω2 and neglects the effect of inter-dumbbell interactions on myosin binding. Because equal oscillation frequencies are the definition of phase-locking, the argument assumes one of the phenomena it is supposed to explain; it can at most show that, once two dumbbells oscillate at the same frequency with a phase difference θ, their average velocities are nonzero. Please provide the full expansion in an appendix, list all small parameters (including the ordering of ẽ/l̄, a/d, and repulsion range), and reframe the text so the numerical simulations, not Eqs. (15)–(16), carry the burden of demonstrating the emergence of phase-locking.","section":"§3.2.1, Eqs. (15)–(16)"},{"comment":"No simulation parameters are given: k1, k3, l0, a, η, τ, c0, μ, ε, n, N, L, and the initial phase and position distributions are all absent, as are the integration scheme, time step, number of realizations per Φ, and equilibration time. Figure 2C even omits the reference μ used to define T0. Without these values the density threshold Φ≈0.4, the exponent α, and the vs∼Φ^4 scaling cannot be checked or reproduced. Please add a parameter table, a numerical-methods subsection, and a code/data availability statement.","section":"§2.3–§2.4, §3.2.2"},{"comment":"The extraction of α and v_s is under-specified. The text notes that individual realizations at low Φ can be sub- or super-diffusive and that only the average over several simulations gives α≈1, but the fitting window, the number of realizations, and the averaging procedure are not stated. The inset in Fig. 4B spans roughly one decade of Φ; please give the fitted range, the fit uncertainty, and a power-law exponent with error bars. This matters because the diffusive-to-ballistic transition is one of the paper's three headline claims.","section":"§3.2.2, Fig. 4"}],"minor_comments":[{"comment":"The reference value of μ used to define T0 is missing in the text ('for c0=1.0 and µ = is taken'); please complete and state all parameter values in the figure caption.","section":"§2.4, Fig. 2C"},{"comment":"The symbol k in Eqs. (15)–(16) is not defined; presumably it denotes a combination of elastic parameters, and r should be related to the center-to-center separation d introduced in §3.2.1.","section":"Eqs. (1), (15)–(16)"},{"comment":"The sentence 'an occurrence of clusters of size Nc > nrequires at least n consecutive dumbbell pairs...' contains a typo and does not explain how n is chosen relative to the threshold Δd_c; please rewrite.","section":"§3.2.3"},{"comment":"The panel labels such as 'A2' and 'B3' are hard to follow; please label subpanels explicitly with (a), (b), etc., in the figure and refer to them accordingly.","section":"§3.3, Fig. 8"},{"comment":"The claim that the qualitative features will hold in two or three dimensions is explicitly an expectation, not a demonstrated result; please mark it as a conjecture.","section":"§4"},{"comment":"Please unify the spelling ('dumbell' vs 'dumbbell') and avoid undefined terms such as 'sessile' when 'non-motile' is meant.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and the idea is appealing, but the revision needs to be judged on whether the authors supply the full numerical details and the hydrodynamic-off control; if the control shows steric synchronization, the paper's framing must be revised substantially. I would not recommend rejection because the missing elements are, in principle, available to the authors within the existing model."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a plausible simulation study showing three density-dependent phenomena in a suspension of oscillating non-motile dumbbells. The new piece is the suspension-level behavior, not the single-agent dynamics. The mechanism claim, however, is under-supported because there is no control run with hydrodynamic coupling switched off.\n\nThe good: The model is well grounded in existing bead-spring and Dierkes-type oscillator work. The three phenomena — diffusive-to-ballistic transition, phase separation, and phase-locking — are clearly demonstrated with reasonable diagnostics. The paper is honest about the 1D setup and the need for extensions. The citation pattern looks fair; they cite earlier work on bead-spring swimmers and synchronization.\n\nThe soft spots, in order. First, the missing control. Equation (7) includes soft repulsion between beads, which becomes active at Φ≥0.4, exactly where phase-locking emerges. Without a simulation with G=0 (no hydrodynamics), you cannot rule out that steric repulsion alone synchronizes the dumbbells. That would still be a result, but it would not be hydrodynamic phase-locking. The paper never addresses this. Second, the analytical pair argument in Eqs. (15)-(16) is stated without derivation, assumes equal frequencies, and ignores the myosin coupling that the numerics show matters. So it is circular as evidence for phase-locking; the numerics carry that weight. Third, the paper does not give parameter values, initial-condition details, integration scheme, or code, making reproduction difficult. Fourth, the vs~Φ^4 scaling is fit from a few points, and the dimensional translation in the Discussion looks off by an order of magnitude: 10^-2 × 20 µm / 10 s should be ~1 µm/min, not 0.1 µm/min. Minor but sloppy.\n\nNone of these are fatal. The central idea is plausible and the qualitative results are probably robust. But the paper needs a revision that adds the hydrodynamic-off control, discloses numerical details, and fixes the dimension error. With that, it could be a solid contribution.\n\nI'd send it to review. The referee should focus on the control and the missing methods. It's not a paradigm shift, but it's a legitimate mechanism worth putting on the table. I wouldn't cite it yet, but I'd bring it to the reading group as a case study in controls for active-matter simulations.","headline":"Plausible simulation-level mechanism for collective motion in non-motile oscillating dumbbells, but the hydrodynamics-centric claim needs a control without hydrodynamic coupling.","tokens_in":14590,"tokens_out":3724,"would_cite":false,"duration_ms":34836,"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":"Hydrodynamic interactions alone can make a suspension of non-motile, oscillating dumbbells move collectively, with a sharp transition near packing fraction 0.4.","keywords":["collective motility","active matter","hydrodynamic interactions","phase-locking","sessile active dumbbells","limit-cycle oscillations","Stokes flow","phase separation"],"falsifier":"Run the same dumbbell dynamics in two or three dimensions, with larger $N$ and varied spring or myosin parameters, and check whether the center-of-mass exponent $\\alpha$ still switches from $\\approx1$ to $\\approx2$ near $\\Phi\\approx0.4$ and whether phase-locking still appears; if the ballistic transition disappears or the velocity scaling departs strongly from $\\Phi^4$, the central claim fails. An experimental equivalent would be to suspend non-adherent oscillating cells (or synthetic oscillatory colloids) at increasing density in a viscous fluid and look for synchronized shape oscillations and a mean drift.","tokens_in":13625,"feed_emoji":"🦠","tokens_out":7827,"duration_ms":69567,"temperature":0.7,"pith_summary":"The paper tries to establish that hydrodynamic interactions alone can turn a collection of non-motile, oscillating objects into a collectively moving suspension. It models each cell as a two-bead dumbbell whose myosin-driven contractions produce limit-cycle oscillations in length; by the scallop theorem such a dumbbell cannot swim on its own. In Stokes-flow simulations of 64 dumbbells arranged in a one-dimensional periodic array, the paper finds that above a suspension density of about $\\Phi=0.4$ the center of mass switches from diffusive to ballistic motion, the suspension phase-separates into clusters, and the oscillations phase-lock. If true, this provides a mechanism for collective cell motility that requires no intrinsic polarity, adhesion, or alignment rule.","feed_headline":"Non-motile dumbbells gain collective motion above packing fraction 0.4","feed_subtitle":"Above packing fraction 0.4, hydrodynamic coupling makes oscillating two-bead dumbbells phase-lock and move ballistically.","key_machinery":"The central object is the active bead-spring dumbbell: two identical beads connected by a nonlinear spring, with myosin-driven contractility represented by a bound-myosin concentration whose dynamics includes turnover and a conservation term; this system undergoes a Hopf bifurcation to limit-cycle oscillations in length and myosin concentration. The second ingredient is the Oseen tensor, the Green's function of Stokes flow that gives the velocity at one bead caused by a force on another bead. The argument is carried by the phase variable $\\psi$ describing each dumbbell's oscillation: pair-level expressions for the cycle-averaged translational velocities show that collective and relative velocities depend on the phase difference $\\theta$, giving phase-dependent attraction and repulsion as long as the phases differ. The simulations then show that hydrodynamic interactions feed back into myosin dynamics and produce phase-locking at high density, with finite phase differences among locked dumbbells.","core_discovery":"The central claim is that hydrodynamic coupling between sessile active dumbbells produces density-dependent collective dynamics: at low densities the dumbbells oscillate independently and the system's center of mass diffuses, while above $\\Phi\\approx0.4$ the dumbbells lock their oscillation phases, form hydrodynamic clusters, and the center of mass moves ballistically with velocity scaling roughly as $v_s\\sim\\Phi^4$. A pair-level analytical argument shows that two dumbbells with a phase difference $\\theta$ acquire a cycle-averaged collective velocity and a phase-dependent effective attraction or repulsion, so the same coupling that creates clusters also drives migration once the phases are coordinated. The simulations show that this coordination is absent below $\\Phi\\approx0.4$, becomes bimodal at intermediate densities (locked or unlocked states), and approaches near-perfect synchronization at $\\Phi\\approx0.8$–$0.9$. The paper concludes that collective motion of a purely sessile population is possible without polarity or adhesion.","pith_inferences":["The pair-level argument implies a testable design rule the paper does not state: a collection of reciprocal oscillators will migrate whenever their phases lock at a nonzero average difference, so tuning density (and thus coupling strength) is a control knob for collective transport without adding polarity.","Because the paper simulates only one-dimensional arrays, whether the transition survives in two or three dimensions is open; if it does, the dumbbell axes should also develop orientational order through hydrodynamic torques, making a 2D or 3D simulation a natural test of the mechanism.","The bimodal mutual-information distributions at intermediate densities suggest bistability between locked and unlocked states; a direct experimental signature would be hysteresis in the phase-correlation as density is swept up and down.","The $\\Phi^4$ scaling is measured only for the specific model parameters and the 64-particle system; checking whether the exponent is robust to changes in $N$, $\\tau$, and bead size would separate a generic hydrodynamic mechanism from a parameter-dependent coincidence."],"forward_implications":["Above $\\Phi\\approx0.4$, the suspension's center of mass crosses from diffusive ($\\alpha\\approx1$) to ballistic ($\\alpha\\approx2$) motion, and the collective velocity grows roughly as $\\Phi^4$.","Hydrodynamic interactions drive phase separation: at $\\Phi\\approx0.5$–$0.8$ the homogeneous starting state breaks into several clusters, and near $\\Phi\\approx0.9$ the dumbbells form a single spanning cluster.","Phase-locking is density-dependent and incomplete: at intermediate densities dumbbells lock at finite phase differences, and at high densities both in-phase and anti-phase locked states are reached depending on initial conditions; the finite phase differences are what generate collective motion.","Translated to cell scales, the collective speed is about $10^{-1}\\ \\mu$m/min, comparable to slow mesenchymal cell motility, but without polarity or adhesion."],"supporting_citations":[{"why":"Supplies the active contractile-element model with Hopf bifurcation and limit-cycle oscillations that each dumbbell inherits.","marker":"(Dierkes et al., 2014)"},{"why":"Establishes the scallop theorem, the reason a single two-bead dumbbell with reciprocal shape changes is sessile.","marker":"(Purcell, 1977)"},{"why":"Provides the bead-spring microswimmer framework that the dumbbell-pair collective-motion argument extends.","marker":"(Pande and Smith, 2015)"},{"why":"Demonstrates hydrodynamic synchronization of flagellar beating, the phenomenon the paper generalizes to phase-locking of dumbbells.","marker":"(Brumley et al., 2014)"},{"why":"Gives generic conditions for hydrodynamic synchronization, supporting the claim that the phase-locking is hydrodynamic in origin.","marker":"(Uchida and Golestanian, 2011)"},{"why":"Defines motility-induced phase separation, which the paper distinguishes from its own hydrodynamics-driven phase separation.","marker":"(Cates and Tailleur, 2015)"},{"why":"Reports shape oscillations in non-adhering fibroblasts, one biological motivation for modeling cells as oscillating dumbbells.","marker":"(Salbreux et al., 2007)"},{"why":"Reports shape oscillations in suspended neutrophils, the other biological inspiration for sessile oscillatory agents.","marker":"(Ehrengruber et al., 1996)"}],"fun_headline_variants":["Hydrodynamic coupling turns oscillating dumbbells into a moving swarm","Sessile dumbbells synchronize and swim when densely packed","Phase-locking gives sessile dumbbells collective motion above Phi=0.4","Oscillating dumbbells phase-lock and move ballistically via hydrodynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a strictly one-dimensional periodic array of 64 identical dumbbells with a single set of model parameters represents real suspensions of oscillating cells; if the phase-locking and ballistic motion require this 1D ordering or those specific parameters, the general claim fails.","fun_headline_variants_meta":{"raw":{"variants":["Hydrodynamic coupling turns oscillating dumbbells into a moving swarm","Sessile dumbbells synchronize and swim when densely packed","Phase-locking gives sessile dumbbells collective motion above Phi=0.4","Oscillating dumbbells phase-lock and move ballistically via hydrodynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000773,"raw_usage":{"total_tokens":3442,"prompt_tokens":983,"completion_tokens":2459,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":2376}},"tokens_in":599,"tokens_out":2459,"duration_ms":17642,"temperature":1.0,"reasoning_tokens":2376,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:30:22.742729+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same dumbbell dynamics in two or three dimensions, with larger $N$ and varied spring or myosin parameters, and check whether the center-of-mass exponent $\\alpha$ still switches from $\\approx1$ to $\\approx2$ near $\\Phi\\approx0.4$ and whether phase-locking still appears; if the ballistic transition disappears or the velocity scaling departs strongly from $\\Phi^4$, the central claim fails. An experimental equivalent would be to suspend non-adherent oscillating cells (or synthetic oscillatory colloids) at increasing density in a viscous fluid and look for synchronized shape oscillations and a mean drift.","supporting_citations":[],"review_version":1}