{"id":"c528cac7-5efe-45a5-8cf0-7c54ede8bdb1","arxiv_id":"2504.14207","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Flagellar wrapping in bacteria starts as a torque-driven buckling instability of a coiled helix, with a critical speed that scales as (R/L)^-3.","lead":"This paper shows, with tabletop experiments and simulations, that a rotating soft helix buckles and wraps around a cylindrical body once the motor speed passes a critical value set by a simple scaling law. The same scaling places three wrapping-capable bacteria in the unstable region, indicating that torque-induced buckling, not a biological switch, starts flagellar wrapping.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The bacterial Fig. 5 placement may be robust, but the prefactor 0.06 in Mc = 0.06(R/L)^-3 is fit, not derived, and Appendix F multiplies omega_c by 1.115; the stability boundary itself is therefore only calibrated to order one, not validated as a quantitative prediction.","rationale":"The reader's CONDITIONAL verdict identifies the most defensible weak point: the bacterial comparison in Fig. 5 rests on treating the bundle as a single helix with Bn = nB1 and on estimated parameters (f, eta, n, L) for the three species, and errors in those inputs could move the points across the boundary. My pass agrees with that as the most load-bearing part of the biological extrapolation. However, the single most generalizable concern is that the boundary itself is a fit: the exponent (R/L)^-3 is derived by scaling, but the coefficient 0.06 is empirical, and Appendix F explicitly adjusts omega_c by 1.115 because of the ambiguity of that prefactor. This matters because the paper's strongest claim, that all wrapping bacteria sit in the unstable region, is a comparison of estimated M values against this calibrated line. It is a quantitative claim rested on a curve whose own uncertainty is not characterized. I do not think this shifts the verdict: the experimental and numerical data collapse onto the line over a wide range, the bacterial points appear to lie well inside the unstable region, and the normal-form control sits in the stable region. The mechanics of the instability and the scaling form are independently supported by the experiments and Stokesian simulations, including the demonstration that long-range hydrodynamic interactions are needed for quantitative agreement. The honest statement is that the biological placement is suggestive rather than a precise quantitative test, which the CONDITIONAL verdict already captures. I would keep the verdict UNCHANGED rather than moving it, because the concern is about the strength of the quantitative claim, not its validity.","tokens_in":22072,"tokens_out":2085,"duration_ms":16627,"concrete_test":"Re-derive the prefactor in Mc = 0.06(R/L)^-3 from the scaling argument in Section IV without using the fit to Fig. 5: compute the numeric coefficient from Eq. (2), the definition zeta = pi eta / log(0.18 P/a), the kink-torque assumption Nc ~ A/R, and the stated relation P ~ A omega / R. If the resulting coefficient differs from 0.06 by more than a factor of ~2, state the boundary with its uncertainty band and re-plot the three bacterial data points (Appendix G values) and the normal-form C. insecticola point against that band; the wrapping-motility claim survives only if all wrapping points remain in the unstable region while the normal form remains in the stable region across the band.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that wrapping onset for bacteria is set by the elastohydrodynamic stability boundary Mc = 0.06(R/L)^-3. The functional form (R/L)^-3 follows from the scaling argument in Section IV, but the prefactor 0.06 is not derived from that argument; it is obtained by fitting the experimental and numerical data in Fig. 5. Separately, Appendix F states that the omega_c values used for the critical-slowing-down analysis were multiplied by 1.115, with the justification that this adjustment is acceptable given the intrinsic ambiguity of the prefactor 0.06. Thus the boundary shown in Fig. 5 has an unquantified order-one uncertainty, and the bacterial comparison in Fig. 5 depends on comparing estimated M values against this fitted line. The reader's stated weakest assumption (bundle stiffness Bn = nB1 and estimated biophysical inputs in Appendix G) is indeed the most uncertain part of the comparison, but it compounds, rather than replaces, the less-examined issue of the empirical prefactor. If the true prefactor were, for example, 0.09 or 0.04, the margin by which the bacterial points sit in the unstable region would change, although all three wrapping-bacteria points appear deep enough in the unstable region that an order-one prefactor shift alone may not flip them. The concern is therefore not that the conclusion is wrong, but that the quantitative boundary used to place bacteria is a calibrated curve rather than a parameter-free prediction, weakening the claim that the agreement between the boundary and the bacterial points is a strong independent test.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a combined experimental, numerical, and theoretical study of torque-induced buckling of a helical filament rotating in a viscous fluid, motivated by flagellar wrapping in bacteria. A macroscale table-top model uses silicone helices in glycerin, and simulations combine the Kirchhoff elastic rod formulation with Stokesian dynamics including long-range hydrodynamic interactions. The central result is a stability diagram in terms of the dimensionless angular frequency M = ηωL^4/A versus the geometric parameter (R/L)^3, with the boundary Mc = 0.06(R/L)^−3, rationalized by a power-balance scaling argument. The authors show that long-range hydrodynamic interactions are necessary to reproduce the experimental dynamics, that available data for wrapping bacteria fall in the predicted unstable region while a normal-form flagellum falls in the stable region, and that the waiting time for buckling onset above threshold scales as ε^−1, consistent with a supercritical Hopf bifurcation. The paper concludes that wrapping onset is a generic motor-torque-induced mechanical instability rather than a dedicated biological control pathway.","tokens_in":22470,"tokens_out":5823,"duration_ms":50857,"significance":"If the result holds, it provides a quantitative physical mechanism for flagellar wrapping, with onset set by geometry, fluid viscosity, filament stiffness, and motor speed. The study's strengths are the clean combination of macroscale experiments, elastohydrodynamic simulation with HIs, and a scaling argument that yields the functional form Mc ∼ (R/L)^−3, plus falsifiable placements of bacterial species on the stability diagram. The paper is also refreshing in explicitly demonstrating, rather than assuming, the importance of long-range hydrodynamic interactions for geometrically nonlinear filament deformation. However, the quantitative predictive content is partly weakened by the empirical prefactor 0.06, the 1.115 rescaling of ωc used in the critical-slowing-down test, and the reliance on estimated biological parameters in Fig. 5; these issues limit the strength of the claim that the boundary is a parameter-free prediction.","major_comments":[{"comment":"The stability boundary is presented as \"the theoretical prediction Mc = 0.06×(R/L)^−3\" and as being \"validated\" by the data, but the prefactor 0.06 is obtained by fitting the experimental and numerical data; the scaling argument in Eq. (2) fixes only the functional form (R/L)^−3. Please state this distinction explicitly and quantify the uncertainty of the prefactor, for example from the spread of data in Fig. 5 or from a collapse analysis. Because the bacterial star symbols in Fig. 5 are compared with this fitted line, the margin by which the bacteria sit in the unstable region should be re-examined for prefactors within the fitted range.","section":"Section IV, Fig. 5"},{"comment":"The claimed ε^−1 critical slowing down in Fig. 6 is obtained after multiplying the experimentally and numerically determined ωc values by 1.115. Since ε is defined with respect to ωc, this rescaling directly affects the abscissa and therefore the fitted slope. Please show the same plot without the 1.115 factor, or with an independently estimated ωc, and report whether the slope −1 is still consistent within the scatter. As written, the quantitative support for Eq. (3) is weaker than the text suggests.","section":"Appendix F, Fig. 6"},{"comment":"The biological comparison rests on several estimated inputs: the P. putida rotation frequency f=50 Hz and viscosity η=0.2 mPa·s are assumed, the flagellar number n=3 is assumed; S. putrefaciens f=50 Hz is estimated from a movie; C. insecticola L=7.1 µm is estimated from an image; and the bundle stiffness is taken as Bn = nB1. None of these uncertainties is propagated into the star symbols in Fig. 5. Please provide a sensitivity analysis giving ranges of M for each species under plausible parameter variations, and add error bars or uncertainty regions to Fig. 5, so the claim that wrapping bacteria lie in the unstable region is quantitatively supported.","section":"Appendix G, Fig. 5"}],"minor_comments":[{"comment":"\"Recent advances in microscopy techniques has uncovered\" should be \"have uncovered\".","section":"Abstract and Section I"},{"comment":"References 9, 22, 23, and 32 contain apparent errors: in Ref. [9] \"Kesseler\" should be \"Kessler\"; in Ref. [22] the volume \"72\" should be \"12\"; Ref. [23] has \"Proc. Natl. Acad. U.S.A.\" missing a period; Ref. [32] is a preprint and should be updated with journal details if available.","section":"References"},{"comment":"In the sentence defining the spontaneous geometry, the parameters are printed as \"κ0\" twice; the second quantity is the spontaneous torsion and should be denoted τ0.","section":"Appendix D"},{"comment":"The quantity P in the power balance is not defined explicitly; if it is the motor power, define it and clarify the dimensional steps leading to P∼Aω/R and to Eq. (2).","section":"Section IV, Eq. (2)"},{"comment":"The dashed line in Fig. 6 is labeled \"0.5/ε\", but Eq. (3) states tc/T∼ε^−1; please explain the factor 0.5 or state explicitly that the line is a guide to the eye rather than a quantitative prediction.","section":"Fig. 6 and Eq. (3)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a strong combination of experiment, simulation, and scaling theory, and the central instability mechanism is credible. My main reservation is that the quantitative boundary and the critical-slowing-down test are, at least in part, calibrated rather than predictive, and the biological conclusions rely on estimated inputs. These issues are fixable with sensitivity analysis and revised wording, so I recommend major revision rather than rejection. I would also encourage the editor to ask the authors to make the simulation code or the data underlying Fig. 5 available, as that would materially strengthen the reproducibility of the biological comparison."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core physical result is real and worth referee time, but the biological comparison is shakier than the authors claim. Using the parameter values in Appendix G, the P. putida point does not sit in the unstable region. With n=3, L=8 μm, R=0.6 μm, f=50 Hz, and η=0.2 mPa·s, I get M≈24, while Mc≈0.06(R/L)^-3≈142. That is a factor of six below the boundary, not above it. So the statement that all three wrapping bacteria fall in the unstable region looks wrong as written. Either the viscosity has a typo (0.2 mPa·s is lower than water) or the boundary prefactor is very different, but in either case the biological claim needs checking.\n\nWhat is genuinely new: the unwinding-rotation geometry around a cylinder, the combination of table-top experiments, Stokesian dynamics, and scaling, and the prediction Mc~(R/L)^-3. The agreement between experiment and simulation for the helix tip trajectories is excellent and parameter-free, and the demonstration that long-range hydrodynamic interactions are essential is convincing. The stability diagram spans enough parameter space that the power-law exponent is credible. The supercritical slowing-down data is a nice addition, though the 1.115 rescaling of ωc in Appendix F is ad hoc.\n\nThe soft spots, in proportion: the prefactor 0.06 is fit, not derived. The scaling argument fixes the exponent but not the constant. That is not fatal, but the paper calls it a theoretical prediction, which overstates it. The bundle stiffness Bn=nB1 is a reasonable first guess but unvalidated, and changes in n shift M linearly. The P. putida issue above is the real problem. If the viscosity was meant to be 2 mPa·s rather than 0.2, the point moves deep into the unstable region and the story holds, but that needs to be corrected and the sensitivity to assumed parameters analyzed.\n\nWho this is for: people working on low-Reynolds-number elastohydrodynamics, flagellar mechanics, and bacterial motility. They will get a solid reference for the buckling instability and a cautionary example of how biological parameter estimation can change a conclusion. I would send it to peer review, but ask for a major revision: fix or justify the P. putida parameters, present the prefactor as calibrated, and avoid claiming a parameter-free boundary.","headline":"Solid elastohydrodynamic buckling result, but the bacterial comparison has a likely parameter error that flips one species into the stable region.","tokens_in":23005,"tokens_out":7675,"would_cite":true,"duration_ms":62643,"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":"The paper shows that the onset of flagellar wrapping is a motor-torque-induced buckling instability of a helical filament, with the stability boundary $M_c = 0.06(R/L)^{-3}$.","keywords":["flagellar wrapping","motor torque","buckling instability","helical filament","elastohydrodynamics","Stokesian dynamics","bacterial motility","rotary motor torque"],"falsifier":"Track a single wrapping flagellum or a controlled bundle of known $n$, $R$, $L$, and stiffness in a viscosity-calibrated fluid, and record the rotation frequency at which the first kink appears; if buckling starts at $M \\neq 0.06(R/L)^{-3}$, or if an independently measured bacterial point lands on the stable side of the line, the central scaling claim is false.","tokens_in":21887,"feed_emoji":"🦠","tokens_out":8556,"duration_ms":73841,"temperature":0.7,"pith_summary":"The paper claims that flagellar wrapping, in which a bacterium coils its flagellum around its body and screws through viscous or confined fluids, begins as a motor-torque-induced buckling instability of a helical filament, not as a separately regulated biological program. It builds a tabletop model of a soft elastomer helix rotating in glycerin around a rigid cylinder and matches its behavior to a Stokesian-dynamics simulation that includes long-range hydrodynamic interactions. Experiment and simulation together collapse onto the stability boundary $M_c = 0.06(R/L)^{-3}$ in terms of the dimensionless drive $M = \\eta\\omega L^4/A$. Placing published bacterial data on that diagram puts three wrapping species in the unstable region and the non-wrapping normal form in the stable region. If the claim is right, the onset of wrapping is set by geometry, fluid viscosity, filament stiffness, and motor speed, and bacteria need only a sufficiently high-torque motor to trigger it.","feed_headline":"Flagellar wrapping starts as a torque-induced buckling instability","feed_subtitle":"Tabletop helix data collapse on one stability line; three wrapping bacteria sit on the unstable side.","key_machinery":"The load-bearing object is the elastohydrodynamic buckling of a helical filament: the motor applies an unwinding torque while viscous drag resists, and when the viscous work rate overcomes the elastic cost of a kink of curvature scale $A/R$, the helix buckles. The controlling group is $M = \\eta\\omega L^4/A$, the ratio of viscous rotation to elastic relaxation, and the stability boundary is $M_c = 0.06(R/L)^{-3}$. The numerical machinery combines the Kirchhoff elastic-rod formulation with Stokesian dynamics using Rotne-Prager mobility tensors for long-range hydrodynamic interactions; switching off the long-range interactions removes the quantitative agreement with experiment. The slow dynamics near onset are described by the supercritical Hopf bifurcation scaling $t_c/T \\sim \\epsilon^{-1}$.","core_discovery":"The central discovery is that a soft helix driven by an unwinding torque in a viscous fluid loses stability through a buckling transition once the nondimensional angular velocity $M = \\eta\\omega L^4/A$ exceeds a critical value that scales as $(R/L)^{-3}$. The instability appears as a localized kink, or perversion, connecting opposite-handed sections of the helix, and when a rigid cylinder is present the buckled helix wraps around it. The paper shows that long-range hydrodynamic interactions are essential for quantitative agreement, because they reduce the effective rotational friction of the flexible helix and help the wrapping complete rapidly. Above onset, the waiting time until buckling slows as $t_c/T \\sim \\epsilon^{-1}$, consistent with a supercritical Hopf bifurcation. The paper concludes that bacteria exploit this motor-driven instability to initiate their wrapping motility, and that the coiled, large-radius polymorphic form is mechanically necessary because the normal form would require a torque beyond any bacterial motor.","pith_inferences":["A direct test not reported in the paper: attach a bead to a single flagellar motor and read the torque-speed curve during wrapping; the critical torque should sit near $0.5A/R$ regardless of species.","The $(R/L)^{-3}$ boundary is a design rule for artificial micro-swimmers: a slightly larger helical radius or a shorter filament lowers the motor speed needed for wrapping, at the cost of bulk.","The model clamps the cell body, whereas a free-swimming body counter-rotates; letting the body rotate in the model may shift the boundary slightly, a natural extension the paper leaves open.","If bundle stiffness indeed scales as $nB_1$, then bundling is itself a mechanical switch: a bacterium can move from stable to unstable by recruiting more filaments into the bundle without changing motor speed."],"forward_implications":["If the scaling claim is correct, the onset of wrapping is set by mechanical inputs alone: geometry, viscosity, filament stiffness, and motor speed.","The coiled polymorphic form is mechanically necessary; with the normal-form radius, the required torque would exceed what bacterial motors can produce.","Wrapping-capable bacteria must have high-torque motors, on the order of the 2000-4000 pN nm values reported for several species.","Raising the viscosity of the medium raises $M$, which explains why the fraction of wrapping cells increases in more viscous media.","The flexible hook at the flagellar base acts as a torque-transmitting universal joint, so real bacteria wrap without the loop seen in the clamped model, and an engineered hook-like joint could let artificial swimmers wrap in confined fluids."],"supporting_citations":[{"why":"Defines the non-dimensional parameter $M=\\eta\\omega L^4/A$ and establishes the rotating-helix experiment that the scaling analysis builds on.","marker":"[26]"},{"why":"Supplies the density-matched elastomer fabrication method and the prior observation of rotating-helix buckling that this model extends.","marker":"[29]"},{"why":"Supplies the Kirchhoff elastic-rod formulation behind the internal forces and torques in the numerical simulation.","marker":"[30]"},{"why":"Supplies the Stokesian-dynamics mobility matrix with Rotne-Prager long-range hydrodynamic interactions used in the simulation.","marker":"[31]"},{"why":"Provides the Shewanella putrefaciens wrapping observations, geometry, and frequency estimate used as one bacterial data point on the stability diagram.","marker":"[20]"},{"why":"Provides the Pseudomonas putida wrapping parameters, multiple-flagella estimate, and viscosity used as another bacterial data point.","marker":"[21]"},{"why":"Provides the Caballeronia insecticola wrapping images, filament geometry, and symbiotic context used as the third bacterial data point.","marker":"[22]"},{"why":"Supports treating a flagellar bundle as a single elastic helix by showing that the bundle stays stable even at stalled-motor conditions.","marker":"[32]"},{"why":"Supplies measured bacterial motor torques against which the predicted critical torque $0.5A/R$ is compared.","marker":"[41]"}],"fun_headline_variants":["Torque-induced buckling triggers flagellar wrapping","Bacteria use torque buckling to wrap flagella","Flagellar wrapping starts with a torque kink","Buckling under torque: the secret to flagellar wrapping","How bacteria corkscrew: torque-driven buckling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The biological comparison assumes that a bundle of $n$ flagella behaves as one elastic helix whose stiffness is $n$ times that of a single filament, and it uses motor speeds, viscosities, and lengths estimated from published images and movies; if those estimates are substantially wrong, the three bacterial data points could cross the stability boundary.","fun_headline_variants_meta":{"raw":{"variants":["Torque-induced buckling triggers flagellar wrapping","Bacteria use torque buckling to wrap flagella","Flagellar wrapping starts with a torque kink","Buckling under torque: the secret to flagellar wrapping","How bacteria corkscrew: torque-driven buckling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000537,"raw_usage":{"total_tokens":2599,"prompt_tokens":984,"completion_tokens":1615,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":600,"completion_tokens_details":{"reasoning_tokens":1543}},"tokens_in":600,"tokens_out":1615,"duration_ms":10626,"temperature":1.0,"reasoning_tokens":1543,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:54:38.643164+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track a single wrapping flagellum or a controlled bundle of known $n$, $R$, $L$, and stiffness in a viscosity-calibrated fluid, and record the rotation frequency at which the first kink appears; if buckling starts at $M \\neq 0.06(R/L)^{-3}$, or if an independently measured bacterial point lands on the stable side of the line, the central scaling claim is false.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the non-dimensional parameter $M=\\eta\\omega L^4/A$ and establishes the rotating-helix experiment that the scaling analysis builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the density-matched elastomer fabrication method and the prior observation of rotating-helix buckling that this model extends."},{"cited_title":"Audoly and Y","cited_arxiv_id":null,"evidence_quote":"Supplies the Kirchhoff elastic-rod formulation behind the internal forces and torques in the numerical simulation."},{"cited_title":"Manghi, X","cited_arxiv_id":null,"evidence_quote":"Supplies the Stokesian-dynamics mobility matrix with Rotne-Prager long-range hydrodynamic interactions used in the simulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Shewanella putrefaciens wrapping observations, geometry, and frequency estimate used as one bacterial data point on the stability diagram."},{"cited_title":"Hintsche, V","cited_arxiv_id":null,"evidence_quote":"Provides the Pseudomonas putida wrapping parameters, multiple-flagella estimate, and viscosity used as another bacterial data point."},{"cited_title":"Kinosita, Y","cited_arxiv_id":null,"evidence_quote":"Provides the Caballeronia insecticola wrapping images, filament geometry, and symbiotic context used as the third bacterial data point."},{"cited_title":"Symbiotic bac- teria break through micrometer-level narrow passages by flagella wrapping,","cited_arxiv_id":null,"evidence_quote":"Supports treating a flagellar bundle as a single elastic helix by showing that the bundle stays stable even at stalled-motor conditions."},{"cited_title":"Sowa and R","cited_arxiv_id":null,"evidence_quote":"Supplies measured bacterial motor torques against which the predicted critical torque $0.5A/R$ is compared."}],"review_version":1}