REVIEW 2 major objections 5 minor 44 references
Collective ballistic motion explains fast aggregation in adhesive active matter
T0 review · 2 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Strong self-alignment in adhesive active matter triggers collective ballistic aggregation, in which the average cluster mass grows as M(t) ~ t², faster than diffusion- or persistence-limited aggregation.
desk verdict Nice simulation study with a plausible persistence-length mechanism for z≈2 in adhesive active matter, but the analytic derivation of the headline exponent has a real gap at γ=0. 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 comparison is between two lengths: the cluster persistence length lp = Vcτp, the distance a cluster travels before its direction decorrelates, and the mean intercluster distance lc = (M/ρ)^{1/2} set by the global density. Self-alignment changes the mass scalings of single-cluster dynamics—strong alignment drives the polar order parameter toward 1, saturates Vc at v0 (γ = 0), and makes τp grow linearly with M (ν = 1)—so lp ∝ M while lc ∝ $M^{{1/2}}$, and the criterion 2lp/lc > 1 marks entry into collective ballistic aggregation, with lp/lc ∝ $M^{{1/2}}$ making that regime self-sustaining. The companion machinery is a generalized Smoluchowski coagulation theory whose kernels encode these scalings, giving z = 1/(1−α) for diffusion-limited aggregation and z = 2/(1−2γ) for ballistic aggregation in two dimensions.
What would settle it
Measure the coalescence rate of two equal-mass clusters moving at the same speed in random directions in the strong-alignment regime; if that rate tends to zero as the speed difference vanishes, the t² law lacks the required source of relative motion and the collective-ballistic explanation, as written, fails.
Extended reading notes
Core claim
The paper's central claim is that the anomalous aggregation exponent z ≈ 2 (mean cluster mass growing as M(t) ~ t²) is produced by a flocking transition inside clusters, not by the persistence of individual particles. Single-cluster simulations in the strong-alignment limit show that the cluster speed becomes independent of mass, Vc = v0 (exponent γ = 0), while the persistence time grows linearly with mass, τp ∝ M (exponent ν = 1). The persistence length lp = Vcτp therefore scales as M, while the average intercluster distance lc = (M/ρ)^{1/2} scales as $M^{{1/2}}$, so the ratio lp/lc grows as $M^{{1/2}}$; once the criterion 2lp/lc > 1 is met, clusters collide ballistically and keep doing so as they grow, making the t² regime asymptotic rather than transient. A generalized Smoluchowski coagulation equation with diffusion kernel ∝ (D(M)+D(M′))(R+R′)^{d−2} and ballistic kernel ∝ |Vc(M)−Vc(M′)|(R+R′)^{d−1} yields the exponent formulas z = 1/(1−α) and z = 2/(1−2γ) in d = 2, reproducing the simulations and uniting diffusion-limited (z ≈ 1/2–1), non-collective ballistic (z ≈ 1), and collective ballistic (z ≈ 2) aggregation in a single picture.
Load-bearing premise
The t² law rests on the unstated premise that clusters of equal mass and equal speed still meet at a nonzero rate because their motion directions differ; if collisions required a difference in speed, equal-speed clusters would never coalesce and the predicted fast regime would not follow.
Editorial extensions
If this is right
- Cell-sorting and tissue experiments that report aggregation exponents between 1 and 2 can be interpreted as sitting in the crossover between diffusion-limited and collective ballistic aggregation, without invoking new physics.
- The criterion 2lp/lc > 1 is directly measurable: tracking cluster centroids to obtain Vc(M) and τp(M) predicts the aggregation regime of a given adhesive active system.
- Because lp/lc keeps growing as M^{1/2} in the collective regime, the t² law is self-stabilizing: larger clusters are even further into the ballistic regime, so no parameter tuning is needed to sustain fast aggregation.
- Systems in collective ballistic aggregation develop a power-law cluster-mass distribution with dynamical scaling (exponent λ ≈ 1.1), whereas diffusion-limited aggregation retains a characteristic cluster size, giving experiments a statistical fingerprint that distinguishes the regimes.
- High single-particle persistence without alignment produces a long-lived transient with z ≈ 1 before the crossover to diffusion-limited behavior, so early-time exponents in biological assays need not reflect the asymptotic regime.
Reading between the lines
- The paper leaves implicit that the ballistic kernel written in the Supplement, |Vc(M) − Vc(M′)|, vanishes when the cluster speed is mass-independent (γ = 0); obtaining z = 2 from the written equations requires the added premise that equal-speed clusters still approach each other with nonzero relative velocity because their directions differ.
- With that premise made explicit, the Smoluchowski formalism extends naturally to an angle-averaged kernel proportional to Vc, which would give an explicit d-dimensional prediction for z(γ, d) instead of the collinear-kernel formula z = d/(1 − dγ).
- A testable extension suggested by the mechanism: reducing adhesion so that clusters can fragment should cut off the linear growth of lp with M and lower the measured exponent—the paper lists fragmentation as future work.
- The persistence-length criterion maps onto an experimental protocol: measuring Vc and τp for cell aggregates of increasing size should show lp(M) turning linear at the same alignment strength at which the aggregation exponent jumps toward 2.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies aggregation of adhesive active Brownian disks with self-alignment, varying the alignment strength J and the Péclet number Pe. It identifies three kinetic regimes: diffusion-limited cluster aggregation (DLCA, z between 1/2 and 1), non-collective ballistic aggregation (z ≈ 1), and a collective ballistic aggregation (CBA) regime with z ≈ 2. Single-cluster simulations are used to extract mass-dependent cluster speed, persistence time, and persistence length, and a criterion comparing persistence length with intercluster distance is proposed to locate the crossover to ballistic aggregation. A generalized Smoluchowski coagulation theory with diffusion- or velocity-based kernels is then used to express the aggregation exponent z in terms of measured single-cluster exponents, and the results are summarized in a Pe–J regime diagram.
Significance. If the CBA mechanism is correct, the paper provides a simple and appealing explanation for anomalously fast aggregation exponents observed in adhesive active matter: flocking makes the cluster persistence length grow linearly with mass while the intercluster distance grows only as the square root of mass, so large clusters move in a collective ballistic regime with M(t) ~ t^2. The paper's strengths are its systematic simulations, the clean single-cluster observable lp(M), the dynamical-scaling collapse of P(M,t), and the unified regime diagram. The main theoretical derivation of the headline exponent z=2, however, has a gap at the precise point γ=0, because the ballistic kernel written in the Supplement vanishes identically for mass-independent cluster speed. The result is likely repairable by adding an explicit orientational-disorder premise, but as written the central claim is not derived from the stated equations.
major comments (2)
- [Supplement, Eqs. (S24)–(S31); main text after Eq. (S31)] The ballistic kernel is written as K(M,M′) = B |M^γ − M′^γ| (M^{1/df} + M′^{1/df})^{d−1}. For the CBA regime the paper uses γ=0, corresponding to mass-independent cluster speed. At γ=0 the factor |M^γ − M′^γ| is identically zero, so I2 in Eq. (S27) vanishes, C2=0, and Eq. (S26) gives d⟨M⟩/dt = 0 rather than ⟨M⟩ ~ t^2. Consequently z = 2/(1−2γ) evaluated at γ=0 is not derivable from the kernel as stated. The intended physical picture must be that clusters have equal scalar speeds but randomly oriented velocities, so the mean relative speed is a nonzero, mass-independent constant. With that premise the kernel becomes K ~ (R+R′)^{d−1}, which for d=df=2 is homogeneous of degree 1/2 and indeed gives z=2. The authors should state this orientational-disorder premise explicitly and revise Eq. (S24) accordingly, since the current text presents γ=0 as a sufficient condition.
- [Main text, 'To explain z...' paragraph and Fig. 3] The agreement between measured and predicted aggregation exponents is a consistency check between two observables of the same simulation, not an independent prediction. The exponents γ, ν, and α are extracted from single-cluster simulations of the same model and then inserted into the Smoluchowski scaling formulas to obtain z. This is acceptable as a scaling analysis, but the text should be more precise: it says 'Smoluchowski theory predicts z=2', whereas in fact the theory uses the measured mass-dependence of cluster motion as input. Please clarify which ingredients are measured and which are derived, so that readers can judge the explanatory content of the framework.
minor comments (5)
- [Main text, after 'We obtain z = 2/(1−2γ)'] The word 'yelding' should be 'yielding'.
- [Supplement, Figs. S1 and S2] The captions of Figs. S1 and S2 appear inconsistent with the surrounding text: the text says Fig. S1 is for low Pe and Fig. S2 for high Pe, but the Fig. S1 caption states 'high Péclet number regime (Pe=10^3)'. Please correct the labeling.
- [Fig. 3] The measured exponents γ, ν, and α are central inputs to the theory, yet the panels in Fig. 3 show no error bars or fit ranges. Adding uncertainty estimates or at least stating how the power-law fits were performed would strengthen the quantitative claims.
- [High persistence, J=0 paragraph] The statement that masses M < 225.73 are in the ballistic regime and larger masses follow DLCA is based on the 2lp/lc criterion, but the M(t) data in Fig. 4a do not show a visible crossover in the displayed time window. A quantitative comparison of the predicted crossover mass with the simulation time scale would make this argument more convincing.
- [General] There is no data- or code-availability statement; providing simulation details or a reproducibility statement would be helpful, especially because several claims (robustness to packing fraction, robustness to alignment rule) are mentioned without a figure.
Circularity Check
The z=2 derivation evaluates a ballistic kernel that vanishes identically at γ=0, and the analytic theory inserts measured single-cluster exponents into standard Smoluchowski scaling; the central analytic claim is not derived from the equations as written.
-
other
[Supplemental Material, 'Ballistic aggregation', Eqs. (S24)-(S29); main text paragraph 'For J ≥ 0.8...']
"K(M, M′) = B |V (M ) − V (M ′)| (R(M ) + R(M ′))d−1 = B |M γ − M ′γ| (M 1/df + M ′1/df )d−1. ... For clusters composed of particles with uncorrelated self-propulsions, γ = −1/2, whereas for fully aligned clusters, γ = 0. ... For J ≥ 0.8, ballistic aggregation occurs with mass-independent cluster speed ( γ = 0), yelding z = 2, consistent with our simulations (Fig. 4a)."
At γ = 0, the relative-speed factor |M^γ − M′^γ| in Eq. (S24) is identically zero for every pair of masses, so the integral I2 in Eq. (S27) vanishes, giving C2 = 0 and d⟨M⟩/dt = 0 in Eq. (S26), not a growing t^2 law. The paper nevertheless evaluates the asymptotic exponent z = 2/(1 − 2γ) at γ = 0, silently replacing the zero prefactor with a nonzero constant. That constant can only come from unstated orientational disorder (equal-speed clusters moving in different directions), which appears neither in Eq. (S24) nor in the stated premise of mass-independent cluster speed. The headline z = 2 is therefore not a consequence of the written equations; the t^2 growth is imported by hand.
-
fitted input called prediction
[Main text, 'Model' section and Conclusions; Fig. 3 caption]
"we will assume power-law dependencies in our theory below in order to capture the limiting values of z, as well as its trend, as J is varied. ... For J ≥ 0.8, ballistic aggregation occurs with mass-independent cluster speed ( γ = 0), yelding z = 2, consistent with our simulations (Fig. 4a). ... With mass-independent cluster speed ( γ = 0, Fig. 3a), Smoluchowski theory predicts z = 2."
The exponent γ = 0 is measured from single-cluster simulations of the same model (Fig. 3a) and then inserted into the standard Smoluchowski scaling formula z = 2/(1 − 2γ) to obtain z = 2, which is compared with aggregation simulations of the same model. This is a consistency check between two simulation analyses of the same particle model, not an independent first-principles prediction of the aggregation exponent. The paper is transparent about assuming the power-law scalings, so the circularity is partial: the Smoluchowski scaling relation is standard and the persistence-length mechanism is separately motivated, but the advertised theory 'predicts' an exponent whose key input (γ = 0) is itself fitted from the model.
full rationale
The numerical study is essentially self-contained: the model is simulated, exponents are measured, and the regime diagram is constructed from those measurements. The self-citations to prior work by the same group (refs. [10] and [23]) are contextual and not load-bearing: the l_c = sqrt(M/ρ) estimate and the observation that alignment modifies the diffusivity-mass relation are supporting background rather than the argument's foundation. No uniqueness theorem or machine-checked external result is invoked to force the paper's choice. The serious problem is internal to the analytic derivation: the ballistic kernel written in Eq. (S24) is exactly zero when γ = 0, the value the main text identifies with the collective ballistic regime, so Eq. (S26) yields d⟨M⟩/dt = 0 unless an unstated nonzero relative-speed term from orientational disorder is added. Evaluating the homogeneous scaling exponent z = 2/(1 − 2γ) at γ = 0 hides this vanishing prefactor. In addition, the analytic theory's inputs (γ, ν, α) are measured from single-cluster simulations of the same model, making the subsequent 'Smoluchowski theory predicts z = 2' a consistency check rather than an independent derivation. The crossover criterion 2l_p/l_c > 1 is a heuristic comparison, not circular by construction, though it is calibrated on the same data. Overall, the empirical aggregation phenomenology is credible and self-contained, but the central analytic claim that collective ballistic aggregation with γ = 0 yields z = 2 is not derivable from the paper's own equations as written; the derivation requires an additional physical premise that is never stated, and the exponents fed into the scaling theory are fitted inputs. This warrants a score of 6 rather than a finding of no circularity.
Assumptions & free parameters
free parameters (3)
- cluster speed exponent γ =
from -1/2 (low J) to 0 (high J)
- cluster persistence time exponent ν =
from 0 (low J) to ≈1 (high J)
- cluster diffusion exponent α =
from -1 (low J) to 1 (high J)
assumptions (5)
- domain assumption Smoluchowski mean-field coagulation equation with spatially homogeneous cluster concentrations.
- domain assumption Dynamical scaling ansatz P(M,t) = M(t)^{-2} f(M/M(t)).
- domain assumption Cluster radius R ∼ M^{1/d_f} with d_f = d = 2; clusters are compact.
- ad hoc to paper Power-law forms Vc ∼ M^γ and τp ∼ M^ν hold over the aggregation range.
- ad hoc to paper In the collective regime, clusters have equal scalar speeds but randomly oriented velocities, supplying a non-zero mass-independent relative speed in the ballistic kernel.
Cite this review
Pith. "Pith review of Collective ballistic motion explains fast aggregation in adhesive active matter." pith.science (2026). https://pith.science/paper/XCEGH6IY
@misc{pith2026250811793,
author = {Pith},
title = {Pith review of: Collective ballistic motion explains fast aggregation in adhesive active matter},
year = {2026},
howpublished = {\url{https://pith.science/paper/XCEGH6IY}},
note = {Machine review of arXiv:2508.11793}
}
read the original abstract
Inspired by motile cells in tissue formation, we find that active systems of self-aligning adhesive particles undergo ballistic aggregation through a flocking transition. This kinetic regime emerges when the cluster persistence length grows faster with cluster mass than the intercluster distance does. We also identify and explain distinct non-collective kinetic regimes, including biologically relevant long-lived transients. Our analytical and numerical results offer a unified framework explaining the broad range of experimentally observed aggregation exponents in cellular systems and reveal physical principles potentially critical for timely tissue organization.
Figures
Reference graph
Works this paper leans on
-
[1]
C. Xiong and S. Friedlander, Proceedings of the National Academy of Sciences 98, 11851 (2001)
work page 2001
-
[2]
A. A. Hyman, C. A. Weber, and F. Jülicher, Annual re- view of cell and developmental biology 30, 39 (2014)
work page 2014
-
[3]
S. V. Franklin and M. D. Shattuck, Handbook of granular materials (CRC Press, 2016)
work page 2016
-
[4]
P. de Castro, F. Urbina, A. Norambuena, and F. Guzmán-Lastra, Physical Review E 108, 044104 (2023)
work page 2023
-
[5]
C. B. Caporusso, L. F. Cugliandolo, P. Digregorio, G. Gonnella, D. Levis, and A. Suma, Physical Review Letters 131, 068201 (2023)
work page 2023
- [6]
-
[7]
G. O’Toole, H. B. Kaplan, and R. Kolter, Annual Re- views in Microbiology 54, 49 (2000)
work page 2000
-
[8]
N. C. Heer and A. C. Martin, Development 144, 4249 (2017)
work page 2017
Show all 44 references
-
[9]
S. G. Krens and C.-P. Heisenberg, Current topics in de- velopmental biology 95, 189 (2011)
2011
-
[10]
E. F. Teixeira, C. P. Beatrici, H. C. Fernandes, and L. G. Brunnet, Physical Review Letters 134, 138401 (2025)
2025
-
[11]
Bothe, E
M. Bothe, E. Lardet, A. Poliakov, G. Pruess- ner, T. Bertrand, and I. Bordeu, arXiv preprint arXiv:2503.13714 (2025)
2025 arXiv
-
[12]
Méhes, E
E. Méhes, E. Mones, V. Nemeth, and T. Vicsek, PloS one 7, e31711 (2012)
2012
-
[13]
Leyvraz, Physics Reports 383, 95 (2003)
F. Leyvraz, Physics Reports 383, 95 (2003)
2003
-
[14]
Nakajima and S
A. Nakajima and S. Ishihara, New Journal of Physics 13, 033035 (2011)
2011
-
[15]
Cremer and H
P. Cremer and H. Löwen, Physical Review E 89, 022307 (2014)
2014
-
[16]
Durand, Plos Computational Biology 17, e1008576 (2021)
M. Durand, Plos Computational Biology 17, e1008576 (2021)
2021
-
[17]
Jiang and F
Y. Jiang and F. Leyvraz, Journal of Physics A: Mathe- matical and General 26, L179 (1993)
1993
-
[18]
P. L. Krapivsky, S. Redner, and E. Ben-Naim, A kinetic view of statistical physics (Cambridge University Press, 2010)
2010
-
[19]
Mones, A
E. Mones, A. Czirók, and T. Vicsek, New journal of physics 17, 063013 (2015)
2015
-
[20]
Baconnier, O
P. Baconnier, O. Dauchot, V. Démery, G. Düring, S. Henkes, C. Huepe, and A. Shee, Reviews of Modern Physics 97, 015007 (2025)
2025
-
[21]
J. M. Belmonte, G. L. Thomas, L. G. Brunnet, R. M. C. de Almeida, and H. Chaté, Phys. Rev. Lett. 100, 248702 (2008)
2008
-
[22]
T. Paul, N. Vadakkayil, and S. K. Das, Physical Review E 109, 064607 (2024)
2024
-
[23]
C. P. Beatrici, R. M. de Almeida, and L. G. Brunnet, Physical Review E 95, 032402 (2017)
2017
-
[24]
M. v. Smoluchowski, Z. Physik. 17, 585 (1916)
1916
-
[25]
Szabo, G
B. Szabo, G. Szöllösi, B. Gönci, Z. Jurányi, D. Selmeczi, and T. Vicsek, Physical Review E 74, 061908 (2006)
2006
-
[26]
Vicsek, A
T. Vicsek, A. Czirók, E. Ben-Jacob, I. Cohen, and O. Shochet, Physical review letters 75, 1226 (1995)
1995
-
[27]
Henkes, Y
S. Henkes, Y. Fily, and M. C. Marchetti, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics 84, 040301 (2011)
2011
-
[28]
Sarkar, G
D. Sarkar, G. Gompper, and J. Elgeti, Communications Physics 4, 36 (2021)
2021
-
[29]
D. L. Barton, S. Henkes, C. J. Weijer, and R. Sknepnek, PLoS computational biology 13, e1005569 (2017)
2017
-
[30]
Martín-Gómez, D
A. Martín-Gómez, D. Levis, A. Díaz-Guilera, and I. Pag- onabarraga, Soft matter 14, 2610 (2018)
2018
-
[31]
Y. Zhao, T. Ihle, Z. Han, C. Huepe, and P. Romanczuk, Physical Review E 104, 044605 (2021)
2021
-
[32]
[ 13, 17, 18, 23, 24, 40–42]
See Supplemental Material at [URL will be inserted by publisher] for videos and analytical details, which in- cludes Refs. [ 13, 17, 18, 23, 24, 40–42]
-
[33]
Peruani, J
F. Peruani, J. Starruß, V. Jakovljevic, L. Søgaard- Andersen, A. Deutsch, and M. Bär, Physical review let- ters 108, 098102 (2012)
2012
-
[34]
Fily and M
Y. Fily and M. C. Marchetti, Physical review letters 108, 235702 (2012)
2012
-
[35]
Y. Fily, S. Henkes, and M. C. Marchetti, Soft matter 10, 2132 (2014)
2014
-
[36]
M. C. Marchetti, Y. Fily, S. Henkes, A. Patch, and D. Yl- lanes, Current Opinion in Colloid & Interface Science 21, 34 (2016)
2016
-
[37]
Villa-Torrealba, C
A. Villa-Torrealba, C. Chávez-Raby, P. de Castro, and R. Soto, Physical Review E 101, 062607 (2020)
2020
-
[38]
de Castro, S
P. de Castro, S. Diles, R. Soto, and P. Sollich, Soft Matter 17, 2050 (2021)
2021
-
[39]
Kolb, Physical review letters 53, 1653 (1984)
M. Kolb, Physical review letters 53, 1653 (1984)
1984
-
[40]
Moncho-Jordá, F
A. Moncho-Jordá, F. Martínez-López, M. Quesada- Pérez, M. Cabrerizo-Vílchez, and R. Hidalgo-Álvarez, Surface and Colloid Science , 113 (2004)
2004
-
[41]
Moncho-Jordá, F
A. Moncho-Jordá, F. Martınez-López, and R. Hidalgo- Álvarez, Physica A: Statistical Mechanics and its Appli- cations 282, 50 (2000)
2000
-
[42]
Moncho-Jordá, G
A. Moncho-Jordá, G. Odriozola, F. Martınez-López, A. Schmitt, and R. Hidalgo-Álvarez, The European Phys- ical Journal E 5, 471 (2001)
2001
-
[43]
Cochet-Escartin, T
O. Cochet-Escartin, T. T. Locke, W. H. Shi, R. E. Steele, and E.-M. S. Collins, Biophysical journal 113, 2827 (2017)
2017
-
[44]
D. A. Beysens, G. Forgacs, and J. A. Glazier, Proceedings of the National Academy of Sciences 97, 9467 (2000). 1 Supplementary Material: Collective ballistic motion explains fast aggregation in adhesive active matter ADDITIONAL SIMULA TION DET AILS We use kc = 200 v0/(σµ), kad...
2000
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.