REVIEW 3 major objections 3 minor 1 cited by
Emergence of Anti-chemotactic Flocking in Active Biomimetic Colloids
T0 review · 3 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict Solid new flocking mechanism in actin-comet colloids, but the PEO experiment doesn't isolate D_m and the model isn't quantitatively calibrated. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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$.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Materials and Methods II.D and Fig. 4I] 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'.
- [Table S2 and Figs. 2B-C, 4H-J] 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.
- [Eqs. S11, S14, S15, and Fig. 3] 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.
minor comments (3)
- [Main text, 'Flocking transition' paragraph] 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.
- [Abstract and Introduction] 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).
- [Model section, Eqs. S14-S15] 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.
Circularity Check
No significant circularity: the model encodes the proposed depletion mechanism, but its predictions (including the xi/W collapse) are nontrivial outputs, not fits, and the central claims are independently grounded in experiments.
full rationale
The paper proposes a reaction-diffusion-phase-field model in which the polymerization rate depends on the local monomer concentration and the active stress is proportional to that rate. This encodes the hypothesized competition-for-monomers mechanism, so the simulations re-deriving attraction are consistency checks rather than independent evidence for the mechanism. However, the manuscript does not stop there: the experimentally measured trends in bead speed versus density, pVCA coverage, bead size, and viscosity are compared with model trends without fitting model parameters to data, and the xi/W=1 phase boundary is an emergent collapse of simulation outputs rather than a fitted prediction. The monomer decay length xi is measured from single-bead simulations and then compared with multi-bead phase behavior, so the boundary is not defined by construction. The viscosity experiment assumes PEO changes only D_m, which is a potential confound and a correctness risk, not a circularity: no parameter is extracted from that experiment and then renamed as a prediction. Self-citations appear only as software use (cuPSS, Ref. [57], by a co-author) and as a background reference on active colloids (Ref. [59]); neither carries the load of the central claim. No uniqueness theorem, ansatz-by-citation, or renaming of a known result was found. Overall, the derivation chain is self-consistent and the experimental grounding is independent, so the observed circularity is minor.
Assumptions & free parameters
free parameters (9)
- actin monomer diffusivity Dm (tilde Dm) =
1-100 (varied)
- nonlinear growth rate k2 (tilde k2) =
1-20 (varied)
- bead radius R (tilde R) =
10-40 (varied)
- active stress factor alpha (tilde alpha) =
5-30 (varied)
- nucleation rate k1 (tilde k1) =
0.1-2
- Hill dissociation constant kd (tilde kd) =
0.3
- depolymerization rate koff (tilde koff) =
0.006
- polymer diffusivity Dp (tilde Dp) =
0.1
- ratio epsilon =
0.2
assumptions (5)
- domain assumption Polymerization rate obeys a Hill function in polymer concentration (Eq. S5).
- ad hoc to paper Actin polymerization generates an active stress proportional to Ron and oriented by the phase-field gradient (Eq. S11).
- domain assumption Beads are represented by Cahn-Hilliard phase fields with steric repulsion and volume conservation (Eqs. S1-S4).
- domain assumption The flow field obeys the overdamped Stokes equation with substrate friction (Eq. S9).
- domain assumption PEO addition changes only Dm, not chemistry or interactions.
Cite this review
Pith. "Pith review of Emergence of Anti-chemotactic Flocking in Active Biomimetic Colloids." pith.science (2026). https://pith.science/paper/62GMYL2X
@misc{pith2026250517394,
author = {Pith},
title = {Pith review of: Emergence of Anti-chemotactic Flocking in Active Biomimetic Colloids},
year = {2026},
howpublished = {\url{https://pith.science/paper/62GMYL2X}},
note = {Machine review of arXiv:2505.17394}
}
read the original abstract
Competition for resources is a fundamental constraint that guides the self-organization of natural, biological, and human systems, ranging from urban planning and ecosystem development to intracellular pattern formation. Here, we reveal that competition for resources is at the origin of the collective dynamics that emerge in a population of colloids propelled by actin treadmilling, an out-of-equilibrium process where filaments grow from one end while shrinking from the other. Using a combination of experiments and theory, we show that symmetry-breaking, self-propulsion, and flocking emerge from the local competition for actin monomers. We demonstrate that beads propelled by actin treadmilling are anti-chemotactic and spontaneously generate asymmetric actin gradients that trigger and sustain directed motility. Flocking emerges from the combined effects of anti-chemotaxis and local competition for monomers. The flocking transition depends on the actin polymerization rate, actin monomer diffusivity, and the bead's motility, whose interplay controls the emergence of short-range attractive interactions between the colloids. Our findings demonstrate that active stress generation coupled to reaction-diffusion is a generic mechanism that can lead to a multiscale cascade of behaviors when active agents remodel their environment. Actin treadmilling offers a platform to study how motile agents that interact through a field self-organize in novel dynamical phases, with potential applications in non-reciprocal and trainable active matter.
Forward citations
Cited by 1 Pith paper
-
Growing, Buckling, and Swirling: motility from polymerization
Polymerization-driven compressive forces trigger long-wavelength buckling of surface filament carpets, and nucleation patterning selects spinning, directed, or chiral swimming of force-free particles.
Reference graph
Works this paper leans on
-
[1]
Skeletal muscle actin from rabbit muscle acetone powder Rabbit skeletal muscle actin is purified following a protocol from Bruce Goode’s lab at Brandeis University based on Spudich and Watt [1] with the following modifications. After centrifugation of the polymerized actin filaments, we resuspend the pellets in 30 mL of G-buffer (5 mM Tris HCl pH 8.0, 0.5...
-
[2]
Capping protein Mus MusculusCapZ is purified using a DNA plasmid supplied by Antoine J´ egou and Guillaume Romet-Lemonne (Addegene plasmid # 89950)[2]. The original protocol [3] utilized Ni-NTA beads which are tranferred to a disposable column, but here we use an ¨AKTA FPLC. After, the cells are lysed in buffer (40 mM NaH 2PO4 pH 7.4, 500 mM NaCl, 30 mM I...
-
[3]
Srv2 His-tagged Srv2 is purified with nickel column affinity with the following modifications to the original protocol [4]. After applying the cell lysate to the nickel column, 5 column volumes of wash buffer containing 50 mM NaH 2PO4 pH 8.0, 300 mM NaCl, 140 mM Imidazole, and 1 mM DTT is flowed over the column. Srv2 is eluted using a linear gradient from...
-
[4]
Human Profilin-1 Human profilin-1 is purified by affinity to poly-L-proline agarose [5]. The protocol followed is adapted from the original source above by the Mullins lab at the University of California San Francisco [6]. Cell lysate is applied to a homemade poly-L-proline sepharose column (following an online protocol written by the previously mentioned...
-
[5]
Human Cofilin-1 Recombinant untagged human cofilin-1 is purified by closely following [7] via cation and anion exchange chromatog- raphy
-
[6]
The plasmid also includes a His-tag
GFP-pVCA The VCA domain of the Whiskott aldrich syndrome protein is expressed tagged with a GFP molecule in BL21 DE3 cells. The plasmid also includes a His-tag. Expressed cells are lysed in lysis buffer (50 mM NaH 2PO4, 300 mM NaCl, 10 mM Imidazole, and 5% glycerol) supplemented with 1 mM PMSF, 1 4 PIC tab, 15 uL DNase, 5mg/ml lysozyme, 2 mM DTT, 0.5% tri...
-
[7]
S. Jansen, A. Collins, S. M. Chin, C. A. Ydenberg, J. Gelles, and B. L. Goode, Nature Communications6(2015), 10.1038/ncomms8202
-
[8]
Sanchez, D
T. Sanchez, D. T. N. Chen, S. J. DeCamp, M. Heymann, and Z. Dogic, Nature491, 431–434 (2012)
2012
Show all 30 references
-
[9]
M. R. Stehnach, N. Waisbord, D. M. Walkama, and J. S. Guasto, Nature Physics17, 926 (2021)
2021
-
[10]
Particle tracking The colloids were segmented and tracked with a subpixel accuracy [11–13]. We define the individual colloid velocities from their displacements over two subsequent frames with time intervalδt:v i(t) =r i(t+δt)−r i(t), where ri(t) andv i(t) are respectively the...
-
[11]
Criteria to detect flocks We define a flock as a group of at least two colloids that are close to each other and move in the same direction. In practice, we apply the following criteria to determine if two colloidsiandjare flocking at timet: •The distance between their center ...
-
[12]
The flocking parameter for a given experiment is the time average of the instantaneous flocking parameter
Flocking parameter The instantaneous flocking parameter is defined as the ratio between the number of colloids belonging to a flock (ni(t)>1) and the total number of colloids in the field of view, at timet. The flocking parameter for a given experiment is the time average of t...
-
[13]
A threshold is applied to each confocal image, and detected objects are filtered by size such that small objects of only a few pixels are removed
Analysis of confocal data pVCA-covered Beads in confocal image stacks are segmented using MATLAB. A threshold is applied to each confocal image, and detected objects are filtered by size such that small objects of only a few pixels are removed. Then, the circularity of each ob...
-
[14]
Such profiles are shown in Fig
Measurement of exponential decay lengthξfrom single-bead simulations Simulations of single motile beads reveals the spatial distribution of actin monomers around a motile bead. Such profiles are shown in Fig. 2F,G of the main text. We measured the decay length of the actin mon...
-
[15]
2B and 2C)
Statistical analysis and number of independent replicates When appropriate, mean and standard deviations of the experimental data are shown as a disk and an errorbar (Insert of Fig. 2B and 2C). In Fig. 4H,I,J, each dot represents the time average flocking parameter for an inde...
-
[16]
Kaschube, and E.F
Martin, A.C., M. Kaschube, and E.F. Wieschaus, Pulsed contractions of an actin-myosin network drive apical constriction. Nature, 2009. 457 (7228): p. 495-9. 17. Fernandez-Gonzalez, R. and J.A. Zallen, Oscillatory behaviors and hierarchical assembly of contractile structures in...
2009
-
[17]
J. A. Spudich and S. Watt, Journal of Biological Chemistry246, 4866–4871 (1971)
1971
-
[18]
Shekhar, M
S. Shekhar, M. Kerleau, S. K¨ uhn, J. Pernier, G. Romet-Lemonne, A. J´ egou, and M.-F. Carlier, Nature Communications 6(2015), 10.1038/ncomms9730
2015 doi
-
[19]
Ulrichs, I
H. Ulrichs, I. Gaska, and S. Shekhar, Nature Communications14(2023), 10.1038/s41467-023-39655-3
2023 doi
-
[20]
Quintero-Monzon, E
O. Quintero-Monzon, E. M. Jonasson, E. Bertling, L. Talarico, F. Chaudhry, M. Sihvo, P. Lappalainen, and B. L. Goode, Journal of Biological Chemistry284, 10923–10934 (2009)
2009
-
[21]
D. A. Kaiser, P. J. Goldschmidt-Clermont, B. A. Levine, and T. D. Pollard, Cell Motility14, 251–262 (1989)
1989
-
[22]
Profilin Purification — Mullins Lab — mullinslab.ucsf.edu,
“Profilin Purification — Mullins Lab — mullinslab.ucsf.edu,”https://mullinslab.ucsf.edu/profilin-purification/ (2003)
2003
-
[25]
K. W. Ebagninin, A. Benchabane, and K. Bekkour, Journal of colloid and interface science336, 360 (2009)
2009
-
[26]
P. J. Lu, P. A. Sims, H. Oki, J. B. Macarthur, and D. A. Weitz, Optics Express15, 8702 (2007)
2007
-
[27]
J. C. Crocker and D. G. Grier, Journal of Colloid and Interface Science179, 298–310 (1996)
1996
-
[28]
The Matlab Particle Tracking Code Respository,
E. D. D Blair, “The Matlab Particle Tracking Code Respository,”https://site.physics.georgetown.edu/matlab/index. html(2008)
2008
-
[29]
Caballero, arXiv preprint arXiv:2405.02410 (2024)
F. Caballero, arXiv preprint arXiv:2405.02410 (2024)
2024 arXiv
-
[30]
Toner, Y
J. Toner, Y. Tu, and S. Ramaswamy, Annals of Physics318, 170 (2005), special Issue
2005
-
[31]
M. E. Cates and J. Tailleur, Annual Review of Condensed Matter Physics6, 219 (2015)
2015
-
[38]
Bugyi, B. and M.F. Carlier, Control of actin filament treadmilling in cell motility. Annu Rev Biophys, 2010. 39 : p. 449-70. 39. McQuillen, R. and J. Xiao, Insights into the Structure, Function, and Dynamics of the Bacterial Cytokinetic FtsZ-Ring. Annual Review of Biophysics, ...
2010
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.