{"id":"6268548a-07d3-45a7-8e15-b65817002960","arxiv_id":"2505.17394","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Actin-propelled beads flock through self-generated actin monomer gradients, with the transition controlled by the balance between monomer diffusion and polymerization.","lead":"Beads that swim by growing actin comet tails form moving flocks because they compete for the same building blocks, the actin monomers. The paper combines experiments with a reaction-diffusion model to show that this resource competition creates short-range attractions, a mechanism that may apply to other active materials and to cells.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Fig. 4I PEO experiment is the only experimental evidence that D_m controls flocking, but 0.8% PEO 400k may alter actin kinetics, induce depletion attraction, or change motility, so the D_m dependence is not isolated.","rationale":"The reader identified the PEO assumption as the weakest point, and I agree this is the most load-bearing concern. The theoretical model predicts a D_m-dependence (Fig. 3H), and the ξ/W=1 collapse is a simulation result, but the abstract and conclusions assert a verified dependence of the experimental flocking transition on D_m. The only experimental evidence for that assertion is Fig. 4I, which is confounded. Without an isolated D_m control, the experimental case for resource competition as the mechanism is weakened, but not destroyed: the pVCA and bead-size experiments are consistent with the model, and the wall accumulation experiment (Fig. 5) is a distinct, non-PEO piece of support. Therefore the concern does not warrant rejection, but it does warrant the CONDITIONAL verdict already assigned. I also considered whether the ξ/W=1 criterion depends on the unphysical phase-field interface width W; while this is a potential issue, it is secondary because the qualitative mechanism does not rest on the precise value of W, and the PEO confound directly threatens the experimental confirmation of the central claim. The concrete test would settle the PEO issue by reproducing the viscosity experiment with a different viscosity agent and measuring D_m and k_on directly.","tokens_in":24354,"tokens_out":7443,"duration_ms":65002,"concrete_test":"Repeat the Fig. 4I viscosity experiment using a viscosity-increasing agent without depletion/crowding chemistry (e.g., glycerol or sucrose) at a concentration matched to the viscosity of 0.8% PEO 400k, while measuring D_m directly by FCS on labeled actin monomers and polymerization kinetics by a pyrene-actin bulk assay in both buffers. If the non-flocking-to-flocking transition does not occur when D_m is reduced without PEO, then the original result does not isolate D_m and the central claim loses its key experimental support.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that flocking emerges from local competition for actin monomers and is controlled by the monomer diffusivity D_m. The only experimental manipulation that targets D_m is the addition of 0.8% (w/v) PEO MW 400k in the viscosity experiment (Materials and Methods II.D; Fig. 4I). The paper assumes this changes only D_m and thereby converts a non-flocking population into a flocking one. However, PEO at this concentration is a non-adsorbing crowder that can (i) alter actin nucleation/elongation rates through excluded-volume effects, (ii) produce weak depletion attraction between colloids with range set by the polymer radius of gyration, and (iii) modify the bead speed and collision dynamics through increased viscous drag. Any of these effects could independently shift the flocking parameter. No control is provided for these alternatives: there is no direct measurement of D_m in the PEO buffer, no pyrene-actin assay of polymerization kinetics in the same buffer, and no test of passive bead pairing under PEO. Thus the experiment does not isolate D_m as the control parameter, and the experimental support for the central claim's diffusivity dependence is not uniquely established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports experiments and a phase-field reaction-diffusion model for micron-sized beads propelled by actin comet tails. The authors show that beads generate asymmetric actin monomer gradients, that they are attracted to regions of monomer depletion (anti-chemotaxis), and that this self-generated gradient leads to short-range attractive interactions and flocking. The model, in which actin polymerization generates a stress proportional to the on-rate, reproduces spontaneous symmetry breaking, speed dependencies, monomer/polymer asymmetry, flocking, and a phase boundary in the (Dm,R) plane that collapses to xi/W=1. Experimental variations of pVCA coverage, viscosity, and bead radius are consistent with the model, and a 3D wall experiment shows accumulation near non-porous but not porous boundaries.","tokens_in":24691,"tokens_out":5702,"duration_ms":64779,"significance":"If confirmed, this work establishes a new, generic mechanism for collective behavior in active matter: chemo-mechanical feedback via self-generated resource depletion, rather than steric or hydrodynamic interactions alone. The paper's strengths include the convergence of several independent experimental trends (bead size, pVCA coverage, wall accumulation) and the public availability of the simulation code and stated data availability. The model is simple and yields a falsifiable collapse criterion (xi/W=1). However, the quantitative link between model and experiment is not yet established, and the only experiment targeting Dm uses a viscosity agent with multiple potential side effects.","major_comments":[{"comment":"The viscosity experiment is the only experimental test of the claim that the flocking transition is controlled by monomer diffusivity Dm. Adding 0.8% (w/v) PEO (MW 400,000) is assumed to change only Dm, but at this concentration PEO can also alter actin nucleation/elongation kinetics through excluded-volume effects, generate weak depletion attraction between colloids with range of the polymer radius of gyration, and change bead speed and collision dynamics through increased drag. No control measurements are provided: there is no direct measurement of Dm in the PEO buffer, no pyrene-actin polymerization assay in the same buffer, and no test of passive bead pairing under PEO. Consequently the non-flocking-to-flocking transition in Fig. 4I does not uniquely establish Dm as the control parameter. Please provide at least one of these controls, or explicitly restrict the claim to 'increased viscosity' rather than 'decreased Dm'.","section":"Materials and Methods II.D and Fig. 4I"},{"comment":"The model is never quantitatively calibrated to the experiments: the dimensionless parameters in Table S2 are given as ranges, and the comparisons in Figs. 2B-C and 4H-J are qualitative, with no model curves overlaid on the experimental data and no error bars from the simulations. Because the central prediction is the collapse of the phase boundary at xi/W=1 (Fig. S7), and xi is measured only in simulations, the paper does not demonstrate that this criterion has quantitative experimental support. Please provide at least one quantitative comparison, for example a model curve for the flocking parameter as a function of Dm, R, or k2 overlaid on the experimental points of Fig. 4H-J, or a mapping from experimental parameters to the dimensionless model parameters.","section":"Table S2 and Figs. 2B-C, 4H-J"},{"comment":"The model builds in the hypothesized mechanism: polymerization consumes monomers in Eqs. S14-S15, and the active stress is proportional to Ron and oriented by the phase-field gradient in Eq. S11. The observation of attraction and flocking in the simulations (Fig. 3) is therefore partly a restatement of the model inputs. The agreement between simulation and experiment should not be presented as independent validation of the monomer-depletion mechanism; the independent evidence comes from the experimental trends (e.g., density-dependent speed in Fig. S5, wall accumulation in Fig. 5D). Please add an explicit acknowledgment that the model assumes the chemo-mechanical coupling, and clearly distinguish which claims rest on the model alone versus on the experiments.","section":"Eqs. S11, S14, S15, and Fig. 3"}],"minor_comments":[{"comment":"In the sentence 'Finally, decreasing the size of the beads led to a transition from a flocking to a non-flocking phase (Fig. 5J)', the reference should be to Fig. 4J, not Fig. 5J.","section":"Main text, 'Flocking transition' paragraph"},{"comment":"The term 'anti-chemotactic' is used repeatedly but never formally defined; please add a one-sentence definition (e.g., motion toward regions of lower chemoattractant concentration, or self-generated gradient sensing that attracts beads to depleted zones).","section":"Abstract and Introduction"},{"comment":"The text states that the Peclet number for actin monomers Pe_m << 1 and therefore convection is turned off in the chemical transport equations; please give an estimate of Pe_m from the experimental parameters (bead speed, monomer diffusivity, bead radius) to justify this approximation.","section":"Model section, Eqs. S14-S15"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a strong fit for cond-mat.soft. The main scientific risk is the PEO experiment, which does not isolate Dm as the control parameter; the authors should be encouraged to add controls or temper the claim. The model's dimensionless parameters are not calibrated, which is acceptable for a mechanism paper but should be stated explicitly. No concerns about novelty disclosure or citation patterns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a genuine new result — actin-comet beads flock through self-generated monomer depletion, not just steric or hydrodynamic alignment — and the evidence is decent. The main soft spot is the PEO experiment, which does not isolate monomer diffusivity as cleanly as the paper claims.\n\nWhat is actually new: the flocking phase for actin-propelled beads, the anti-chemotactic response, and the xi/W=1 criterion that organizes the phase boundary. Prior work on autochemotactic droplets and actin network competition did not report this collective phase. The model adapts existing phase-field approaches (refs 54-56), but the coupling of reaction-diffusion to active stress for this system is sensible and reproduces the main experimental trends: speed versus bead size, pVCA coverage, monomer asymmetry at the trailing edge, and wall accumulation. The SI also explicitly rules out hydrodynamic interactions and MIPS, which is honest and helpful.\n\nWhat the paper does well: multiple independent experiments support the mechanism — bead size, pVCA surface coverage, density-dependent speed and tail length, quasi-2D confinement necessity, and the wall-versus-porous-boundary accumulation. The model is clearly laid out, dimensionless parameters are listed, and the xi/W collapse across about 100 simulations is a clean way to organize the phase diagram. The experiments and simulations are presented side by side, and the qualitative agreement is real.\n\nSoft spots, in order of severity:\n\n1. The PEO viscosity experiment (Fig. 4I) is the only experimental handle on D_m, and 0.8% PEO 400k is not a clean knob. PEO can alter actin nucleation and elongation through excluded-volume effects, induce depletion attraction between beads, and change bead drag. The paper provides no control — no direct D_m measurement, no pyrene-actin assay, no passive bead pairing test. The stress-test note is correct. This does not kill the paper, because the pVCA and bead-size experiments support the mechanism independently, but the specific claim that \"flocking is suppressed when diffusivity dominates\" leans heavily on that one assay.\n\n2. Model parameters are not calibrated to experiments. The dimensionless parameters are plausible, but there is no mapping to physical units, so the xi/W=1 boundary is an empirical simulation collapse rather than a quantitative prediction. This is a limitation, not a fatal flaw.\n\n3. Data and code are promised on Dryad but not yet deposited. For a preprint that is common, but for a serious review it should be available.\n\nThe circularity concern — the model builds in the hypothesized mechanism — is real but mitigated: the experiments stand independently, and the simulations demonstrate plausibility rather than prove the mechanism.\n\nWho this is for: active matter experimentalists and theorists, soft matter physicists, biophysicists. It deserves a serious referee. The central claim is likely right, and the gaps are addressable. I would engage with it, but I would want the PEO control and a quantitative parameter mapping before citing it as established.","headline":"Solid new flocking mechanism in actin-comet colloids, but the PEO experiment doesn't isolate D_m and the model isn't quantitatively calibrated.","tokens_in":25238,"tokens_out":2059,"would_cite":true,"duration_ms":17935,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Actin-comet-propelled colloids are drawn together by the monomer-depletion gradients they generate themselves, with the free-swimming-to-flocking switch governed by a single length-scale ratio, ξ/W = 1.","keywords":["actin treadmilling","anti-chemotaxis","flocking","reaction-diffusion","phase-field model","active colloids","monomer depletion","biomimetic motility"],"falsifier":"Image the actin monomer field directly between two approaching beads using fluorescently labeled G-actin: the model requires a measurable depletion well between them whose decay length $\\xi$ crosses the flocking boundary at $\\xi/W = 1$. Observing stable flocks while the inter-bead monomer profile stays flat and symmetric, or finding that the flocking transition stops tracking $\\xi/W = 1$ when monomer diffusivity is varied by a chemically inert route such as temperature or solvent isotope, would falsify the central mechanism.","tokens_in":24191,"feed_emoji":"🦠","tokens_out":22311,"duration_ms":155302,"temperature":0.7,"pith_summary":"Beads propelled by actin treadmilling — the same comet-tail mechanism that drives pathogens such as Listeria — interact through the fuel they consume. Each bead depletes actin monomers from the solution to build its tail, and because the local monomer concentration sets the polymerization rate and hence the propulsive stress, a bead is steered by the very gradient it creates: it is anti-chemotactic, drawn toward the depleted region, which for two neighbors is the gap between them. In quasi-2D confinement this attraction turns collisions into finite-size 'U'-shaped flocks that constantly merge and split. The paper captures this in a phase-field reaction-diffusion model with polymerization-driven active stress, and shows that about 100 simulations collapse onto one phase diagram whose flocking boundary is $\\xi/W = 1$, where $\\xi$ is the decay length of the monomer gradient flanking a bead and $W$ is the bead's effective surface thickness. Experiments varying bead size, nucleation density, and solution viscosity cross that boundary in the predicted directions, and in 3D the same depletion logic accumulates beads on impermeable walls but not on porous boundaries.","feed_headline":"Actin-propelled beads flock by stealing each other's fuel","feed_subtitle":"The switch between solo swimming and flocking collapses onto one length-scale ratio: ξ/W = 1.","key_machinery":"The load-bearing object is a chemomechanical reaction-diffusion phase-field model. Each bead is a phase field $\\phi_i$, and monomeric and polymerized actin ($c_m$ and $c_p$) diffuse with diffusivities $D_m$ and $D_p$, reacting through the surface-localized polymerization rate $R_{\\mathrm{on}} = c_m(k_1 + k_2 c_p^2/(k_d^2 + c_p^2))$ — a two-stage nucleation-and-growth form with Hill-type saturation — balanced by depolymerization at rate $k_{\\mathrm{off}}$. Bead motion obeys the overdamped Stokes equation with substrate friction and an active stress $-(\\alpha R_{\\mathrm{on}} + \\Gamma)\\nabla\\phi_i \\otimes \\nabla\\phi_i$, which coarse-grains the actin polymerization ratchet into a surface tension whose gradients drive motion like a Marangoni flow; the sign makes a bead steer toward the depleted side where polymerization is weaker. The decisive derived quantity is $\\xi$, the shorter exponential decay length of the monomer profile along the bead's flank, measured in single-bead simulations; the central result is that multi-bead simulations flock precisely when $\\xi/W > 1$, with $W$ the phase-field interface width.","core_discovery":"The central claim is that the collective motion of these biomimetic colloids is chemomechanical: a bead's actin comet consumes monomers, the resulting asymmetric monomer gradient modulates the local polymerization rate $R_{\\mathrm{on}}$, and because polymerization exerts an active stress normal to the bead surface, gradients in $R_{\\mathrm{on}}$ become forces that steer the bead toward the depleted region. This anti-chemotactic coupling is what breaks symmetry and sustains directed motion for an isolated bead, and it is what draws two beads together: the depletion well between them weakens their facing sides and reorients them toward each other, producing a short-range effective attraction. The quantitative claim is a collapse: in the $(D_m, R, k_2)$ parameter space, motile beads flock exactly when the exponential decay length $\\xi$ of the monomer gradient perpendicular to the direction of motion exceeds the phase-field interface thickness $W$, with flocking, non-flocking, and non-motile regions separated by $\\xi/W = 1$. The experiments match the model's three dials: lowering the nucleation factor pVCA, raising the solution viscosity (lowering monomer diffusivity $D_m$), and shrinking the beads all suppress flocking as predicted. In 3D, the same mechanism predicts — and the experiments confirm — that beads accumulate along a solid wall, which blocks monomer replenishment and creates the asymmetric gradient, but not along a porous membrane connected to a monomer reservoir. The paper concludes that active stress generation coupled to reaction-diffusion is a generic route from single-agent motility to collective pattern formation when active agents remodel their environment.","pith_inferences":["The ξ/W = 1 criterion could be pinned down more cleanly by varying monomer diffusivity through a chemically inert route such as temperature or solvent isotope, so that the diffusion length changes without any possibility of altering the actin chemistry or the bead's propulsion.","The ξ/W criterion is a transferable design rule: in any suspension of fuel-consuming motile particles, the flocking threshold should be set by the ratio of the fuel's depletion length to the particle's effective surface thickness — testable in catalytic colloids, enzyme-coated swimmers, or bacterial suspensions without actin.","Because a bead's depletion field depends on its own consumption rate and size, the attraction between two unequal beads need not be mutual; the model's logic implies non-reciprocal interactions between beads of different radii or activities, which could produce chasing, sorting, or asymmetric flock morphologies that the paper does not report.","The wall-accumulation result points to a generic physical mechanism for boundary localization of resource-consuming motile agents wherever their fuel is not replenished — a mechanism the paper's framing invites one to look for in intracellular organization and in ecological or microbial pattern formation, though it does not test those settings."],"forward_implications":["The flocking transition is governed by a single dimensionless ratio: motile beads flock when their self-generated monomer gradient reaches beyond the bead's effective surface (ξ/W > 1), and swim past one another when it does not, even though collisions still reorient them.","Three independent experimental dials — bead radius, pVCA nucleation density, and fluid viscosity (monomer diffusivity) — move the system across the ξ/W = 1 boundary in the predicted directions, so a single control parameter organizes the outcomes of three separate experiments.","Quasi-2D confinement is necessary but not sufficient for flocking: steric realignment only enables the attraction to act, and without a long-enough monomer gradient small beads collide and reorient but never flock.","The same depletion mechanism produces a distinct 3D signature: beads accumulate on impermeable walls, which block monomer replenishment, but not on porous boundaries open to a monomer reservoir.","Because the mechanism only requires that motile agents consume a diffusing resource, it should generalize to other active systems: any self-propelled object that depletes its own fuel should attract neighbors through depletion gradients and flock when the depletion length exceeds its interaction surface."],"supporting_citations":[{"why":"Reconstitutes actin-based motility of Listeria and Shigella from purified proteins; supplies the experimental assay the bead-comet system is built on.","marker":"[51]"},{"why":"Establishes the barbed-end polymerization mechanism behind comet-tail propulsion, which the model's surface-localized reaction rate coarse-grains.","marker":"[42]"},{"why":"Supplies the phase-field chemomechanical modeling strategy — coupling actin reaction-diffusion to morphology and motion — that the paper's model adapts.","marker":"[54]"},{"why":"Provides the elastic-ratchet theory of polymerization-driven force generation that the active stress term represents at coarse-grained level.","marker":"[47]"},{"why":"Shows that local competition for actin monomers can slow and steer actin network growth; the flocking mechanism extends this to bead-bead interactions.","marker":"[29]"},{"why":"Reports that global competition for a limited actin pool affects network growth and structure; a direct precursor for the resource-competition premise.","marker":"[30]"},{"why":"Demonstrates chemotaxis and autochemotaxis in self-propelled droplets, the closest field-mediated-interaction system this work distinguishes from monomer-depletion flocking.","marker":"[34]"},{"why":"Models actin-based propulsion of soft droplets through surface stress, motivating the paper's representation of polymerization as an effective active surface tension.","marker":"[58]"},{"why":"Provides the pseudo-spectral PDE solver used to run the phase-field simulations from which the decay length and the phase diagram are computed.","marker":"[57]"}],"fun_headline_variants":["Flocking colloids chase depleted actin fuel","One length ratio predicts when actin-propelled beads flock","Anti-chemotactic beads flock by chasing depleted actin","Actin depletion drives flocking in biomimetic colloids","Beads flock when monomer depletion exceeds a critical length"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The experiments that attribute the flocking transition to monomer diffusivity assume that adding 0.8% polyethylene oxide changes only how fast actin monomers diffuse through the fluid; if the polymer also alters the actin polymerization chemistry, the strength of the polymerization-generated force, or the effective interactions between beads, then the transition may not be isolating the theory's chosen control parameter $D_m$.","fun_headline_variants_meta":{"raw":{"variants":["Flocking colloids chase depleted actin fuel","One length ratio predicts when actin-propelled beads flock","Anti-chemotactic beads flock by chasing depleted actin","Actin depletion drives flocking in biomimetic colloids","Beads flock when monomer depletion exceeds a critical length"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001049,"raw_usage":{"total_tokens":4494,"prompt_tokens":1121,"completion_tokens":3373,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":737,"completion_tokens_details":{"reasoning_tokens":3295}},"tokens_in":737,"tokens_out":3373,"duration_ms":19727,"temperature":1.0,"reasoning_tokens":3295,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:47:59.089865+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Image the actin monomer field directly between two approaching beads using fluorescently labeled G-actin: the model requires a measurable depletion well between them whose decay length $\\xi$ crosses the flocking boundary at $\\xi/W = 1$. Observing stable flocks while the inter-bead monomer profile stays flat and symmetric, or finding that the flocking transition stops tracking $\\xi/W = 1$ when monomer diffusivity is varied by a chemically inert route such as temperature or solvent isotope, would falsify the central mechanism.","supporting_citations":[{"cited_title":"cuPSS: a package for pseudo-spectral integration of stochastic PDEs","cited_arxiv_id":"2405.02410","evidence_quote":"Shows that local competition for actin monomers can slow and steer actin network growth; the flocking mechanism extends this to bead-bead interactions."},{"cited_title":"Toner, Y","cited_arxiv_id":null,"evidence_quote":"Reports that global competition for a limited actin pool affects network growth and structure; a direct precursor for the resource-competition premise."}],"review_version":1}