REVIEW 4 major objections 4 minor 66 references
Kinetically arrested clusters in active filament arrays
T0 review · 4 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read A pinned-free array of active elastic filaments self-assembles into regularly spaced, kinetically arrested clusters whose size and spacing are set by activity, elasticity, and grafting density.
desk verdict Useful simulation study of clustered active filament carpets; the clustering phenomenology is solid, but the headline cluster-size scaling law is not yet supported by independent parameter sweeps. 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 central object is the dimensionless activity parameter $\beta = f \sigma^2 (N_m-1)^3/\kappa$, the ratio of active compressive forcing to bending stiffness, together with the dimensionless anchoring spacing $\Delta = \delta/\sigma$. The argument also depends on a second structural feature: setting the bond length equal to the bead diameter, $\ell_0=\sigma$, makes each filament a corrugated chain, so bead-bead contacts supply tangential resistance that acts as friction. These two ingredients do the work: $\beta$ drives the global rotation and buckling instability of individual filaments, while $\Delta$ controls how many rotated filaments can be packed before contact torques balance the active torques. The paper's scaling theory combines the single-filament coil radius $\sim \ell/\beta^{1/3}$ with the grafting spacing to arrive at the cluster size.
What would settle it
Repeat the same parameter sweep with continuous smooth filaments or with overlapping beads ($\ell_0<\sigma$) so that tangential bead-locking is absent, holding $\beta$ and $\Delta$ fixed; if clusters no longer appear or the mean cluster size no longer follows the predicted scaling, the corrugation-friction mechanism is the load-bearing part of the claim.
Extended reading notes
Core claim
The paper's central claim is that activity alone, with no attractive interactions, organizes an array of pinned-free elastic filaments into a periodic pattern of compact towers. For activity parameter $\beta$ above the single-filament buckling threshold, each filament buckles and rotates; normal steric repulsion prevents free rotation, and tangential contact forces produced by the bead-scale corrugation lock neighboring filaments into arrested bundles. The authors show that cluster outlines at different grafting densities collapse onto a common shape after rescaling lateral widths by $\Delta^{a(\beta)}$, and they derive a scaling expression in which the mean number of filaments per cluster, $\langle N_f\rangle$, decreases with $\beta$ and with the anchoring spacing $\Delta = \delta/\sigma$.
Load-bearing premise
The arrest mechanism requires the filament surfaces to be corrugated: only because the bead spacing equals the bead diameter do tangential contact forces mimic friction and lock filaments in place; with smoother filaments the clusters might glide past each other or oscillate instead of jamming.
Editorial extensions
If this is right
- Above a threshold activity (around $\beta \approx 29.7$ in the simulations), sufficiently large arrays form stable kinetically arrested clusters, while arrays below the threshold tilt or fluctuate without clustering.
- For fixed activity, decreasing the anchoring spacing $\Delta$ produces wider clusters containing more filaments, whereas increasing $\beta$ produces more compact tower-like clusters.
- Cluster outlines for different $\Delta$ can be collapsed onto a self-similar shape by rescaling lateral widths by $\Delta^{a(\beta)}$, with $a(\beta)$ decreasing as activity increases.
- The average number of filaments in a cluster decreases monotonically with $\Delta$ and follows the scaling-theory prediction, while the number of clusters grows approximately linearly with $\Delta$ at a slope set by $\beta$.
- Filaments that oscillate between adjacent clusters do so with a well-defined frequency, and the geometric constraint of arrest makes pinned filaments behave like effectively clamped, shorter filaments.
Reading between the lines
- The geometric friction mechanism suggests a tunable design axis: changing bead shape, size, or bond length relative to bead diameter while holding $\beta$ and $\Delta$ fixed should continuously vary how strongly clusters are arrested.
- Because the simulations use the freely draining approximation, full hydrodynamic coupling between clustered filaments could alter the arrest threshold; comparing cluster statistics with and without pair-mobility hydrodynamic interactions would test that sensitivity.
- If the scaling law holds beyond the simulated range, it implies that cluster size can be actively toggled by modulating activity over time, allowing reversible assembly and disassembly of filament-based materials.
- The effectively clamped behavior of arrested filaments suggests that inter-cluster oscillators could be engineered as localized, synchronized actuators embedded in an otherwise static filament carpet.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript reports Brownian dynamics simulations of a two-dimensional array of 300 active elastic bead-spring filaments, pinned at one end and free at the other, with compressive follower forces acting along each filament and WCA steric interactions between beads. For inter-filament spacing Δ = δ/σ between 2 and 5 and activity β = f σ^2 (Nm−1)^3/κ above a threshold, the authors observe that the array self-assembles into regularly spaced, compact clusters that appear kinetically arrested. They characterize cluster outlines by rescaling them with an activity-dependent exponent a(β), and they propose a scaling theory, summarized as a dashed line in Fig. 4(a), that relates the mean number of filaments per cluster to β and Δ. They also compute local pressure fields and active, elastic, and contact torques in an effort to identify the mechanism of cluster stabilization.
Significance. If the central quantitative claims hold, the paper provides a useful design principle: cluster size and spacing in pinned active filament arrays are controlled by β and Δ. The simulations are presented with some care: cluster-size data are averaged over time and initial conditions with error bars, and clustering appears robust across Δ = 2–5 and β in the range 23.97–59.32. However, the quantitative theory is not self-contained in the main text: the scaling law for ⟨Nf⟩ is not displayed, and the dependences on elasticity and filament length are not independently tested. These gaps currently make the predictive claim hard to assess rather than directly contradicted.
major comments (4)
- [§4.3, Fig. 4(a)] The abstract and Introduction promise 'theoretical expressions' for the number of filaments in a cluster and the spacing between clusters, but the main text never writes such an expression. The only quantitative comparison is the dashed line in Fig. 4(a), described as the prediction from the scaling theory, with details relegated to the ESM. Without the equation, including prefactors and any fitted constants, a reader cannot reproduce the prediction or see how β and Δ enter. The inter-cluster spacing prediction is similarly not shown. This is load-bearing because the claimed derivation of a scaling law is a central advertised result.
- [Table 1 and §4.3] The simulation campaign varies f and δ, while κ = 2×10^4, Nm = 40, and ℓ = 39 are fixed in the main runs. The dependence of cluster size on elasticity and filament length is therefore only implied through the combination β. If, for example, κ were doubled and f were doubled so that β remained fixed, the current data cannot rule out a different cluster size. Given the abstract's claim that elasticity crucially influences cluster size, at least a few runs at different κ and Nm (or an explicit dimensional argument showing why such sweeps are unnecessary) are needed.
- [§3.4, Fig. 3(d)–(f)] The self-similar exponents a(β) are reported as best-fit values 1.6, 1.3, and 1.05 with no uncertainties, and the collapse is visibly imperfect near the pinned ends, as the authors acknowledge. Since these exponents are one of the paper's quantitative outputs, the authors should report fit uncertainties, define a collapse metric, and state how many clusters and Δ values contribute to each fit.
- [§3.5 and Fig. 4] The notion of 'kinetic arrest' is not quantitatively established. Cluster size counts only filaments that are always part of a cluster, explicitly excluding inter-cluster oscillating filaments, but no operational criterion or time scale for arrest (for example, a plateau in cluster-member mobility or an upper bound on detachment rate) is reported. Because the title and central claim depend on the arrest mechanism, a quantitative definition of arrest should be added.
minor comments (4)
- [Eq. (4)] Equation (4) writes a gradient of Φα_κ for the bending term, but the bending potential is defined as Φα_B in Eq. (2); the notation should be made consistent.
- [Fig. 5 caption] The caption lists β = 44.94 while the text and other figures use β = 44.49; please check which value is correct.
- [Fig. 4(b)] The text quotes linear slopes for ⟨Nc⟩ versus Δ as 3.8, 5.5, and 6.8 without uncertainties; since these slopes are presented as quantitative trends, giving confidence intervals would be useful.
- [§3.2] The discussion of 'short system sizes' suppressing clustering does not report the system sizes tested or show the corresponding data; a supplementary figure or table would make this claim checkable.
Circularity Check
No circularity found; the central clustering observation is a direct simulation result and the scaling theory is not shown to reduce to its own inputs.
full rationale
The paper's load-bearing claims are produced by Brownian dynamics simulations of a defined model, and the scaling theory is invoked only after those simulations produce the clustered states. The single-filament instability threshold (βc ≈ 30.6) is imported from earlier work by the same group, but that result is a parameter-free prediction from the stated follower-force model, is re-examined in Section 4.1 via the solvability condition, and is consistent with the current f = 10 simulations (β = 29.66, near threshold), so it counts as independent support rather than circular input. The tangential 'effective friction' mechanism is an explicit modeling choice (ℓ0 = σ making the filament surface corrugated), not a fitted parameter renamed as a prediction. The self-similar collapse exponents a(β) are reported as best-fit descriptions of observed cluster outlines, not as predictions derived from themselves. The central quantitative claim, ⟨Nf⟩ as a function of β and Δ, is stated as a scaling theory, but the provided text does not display a closed-form relation or its fitted constants; consequently no specific equation can be exhibited that reduces the prediction to the simulation data. The skeptic's concern that κ and Nm are not independently swept is a correctness and generality risk, not a circularity. Under the hard requirement of quoting a specific reduction, no circular step is established.
Assumptions & free parameters
free parameters (2)
- Cluster shape scaling exponent a(β) =
a ≈ 1.6 for β = 29.66; 1.3 for β = 44.49; 1.05 for β = 59.32
- Empirical slope of ⟨Nc⟩ versus Δ =
3.8, 5.5, and 6.8 for β = 29.66, 44.49, and 59.32
assumptions (4)
- domain assumption Freely draining limit: hydrodynamic interactions between beads and filaments are neglected; drag is local and isotropic.
- domain assumption Filament inextensibility in the analytical stability and scaling theory.
- domain assumption Equating bead spacing and bead diameter (ℓ0 = σ) to generate corrugation-induced tangential friction.
- ad hoc to paper Replacement of the distributed follower force by a concentrated tip force in the simplified linear stability analysis.
Cite this review
Pith. "Pith review of Kinetically arrested clusters in active filament arrays." pith.science (2026). https://pith.science/paper/VAS5BDY7
@misc{pith2026241220536,
author = {Pith},
title = {Pith review of: Kinetically arrested clusters in active filament arrays},
year = {2026},
howpublished = {\url{https://pith.science/paper/VAS5BDY7}},
note = {Machine review of arXiv:2412.20536}
}
read the original abstract
We use Brownian dynamics simulations and theory to study the over-damped spatiotemporal dynamics and pattern formation in a fluid-permeated array of equally spaced, active, elastic filaments that are pinned at one end and free at the other. The filaments are modeled as connected colloidal chains with activity incorporated via compressive follower forces acting along the filament backbone. The length of the chains is smaller than the thermal persistence length. For a range of filament separation and activity values, we find that the filament array eventually self-assembles into a series of regularly spaced, kinetically arrested, compact clusters. Filament activity, geometry, elasticity, and grafting density are each seen to crucially influence the size, shape, and spacing of emergent clusters. Furthermore, cluster shapes for different grafting densities can be rescaled into self-similar forms with activity-dependent scaling exponents. We derive theoretical expressions that relate the number of filaments in a cluster and the spacing between clusters, to filament activity, filament elasticity, and grafting density. Our results provide insight into the physical mechanisms involved in the initiation of clustering and suggest that steric contact forces and friction balance active forces and filament elasticity to stabilize the clusters. Our simulations suggest design principles to realize filament-based clusters and similar self-assembling biomimetic materials using active colloids or synthetic microtubule-motor systems.
Figures
Reference graph
Works this paper leans on
-
[1]
G. M. Whitesides and B. Grzybowski, Science, 2002, 295, 2418--2421
work page 2002
-
[2]
V. N. Manoharan, Science, 2015, 349, 1253751
work page 2015
-
[3]
L. L. Ong, N. Hanikel, O. K. Yaghi, C. Grun, M. T. Strauss, P. Bron, J. Lai-Kee-Him, F. Schueder, B. Wang, P. Wang et al., Nature, 2017, 552, 72--77
work page 2017
-
[4]
W. B. Rogers, W. M. Shih and V. N. Manoharan, Nature Reviews Materials, 2016, 1, 1--14
work page 2016
-
[5]
W. Zhou, Y. Lim, H. Lin, S. Lee, Y. Li, Z. Huang, J. S. Du, B. Lee, S. Wang, A. S \'a nchez-Iglesias et al., Nature materials, 2024, 23, 424--428
work page 2024
-
[6]
A. McMullen, M. Mu \ n oz Basagoiti, Z. Zeravcic and J. Brujic, Nature, 2022, 610, 502--506
work page 2022
-
[7]
Y. Cui, H. Zhu, J. Cai and H. Qiu, Nature Communications, 2021, 12, 5682
work page 2021
-
[8]
M. Liu, X. Zheng, V. Grebe, D. J. Pine and M. Weck, Nature Materials, 2020, 19, 1354--1361
work page 2020
Show all 66 references
-
[9]
Nishiguchi, J
D. Nishiguchi, J. Iwasawa, H.-R. Jiang and M. Sano, New Journal of Physics, 2018, 20, 015002
2018
-
[10]
Huang, C
Y. Huang, C. Wu, J. Chen and J. Tang, Angewandte Chemie International Edition, 2024, 63, e202313885
2024
-
[11]
Bricard, J.-B
A. Bricard, J.-B. Caussin, N. Desreumaux, O. Dauchot and D. Bartolo, Nature, 2013, 503, 95
2013
-
[12]
Palacci, S
J. Palacci, S. Sacanna, A. P. Steinberg, D. J. Pine and P. M. Chaikin, Science, 2013, 339, 936--940
2013
-
[13]
A. E. Patteson, A. Gopinath and P. E. Arratia, Current Opinion in Colloid & Interface Science, 2016, 21, 86--96
2016
-
[14]
S. A. Mallory, C. Valeriani and A. Cacciuto, Annual review of physical chemistry, 2018, 69, 59--79
2018
-
[15]
Z. Wang, Z. Wang, J. Li, C. Tian and Y. Wang, Nature communications, 2020, 11, 2670
2020
-
[16]
Ramaswamy, Annual Review of Condensed Matter Physics, 2010, 1, 323--345
S. Ramaswamy, Annual Review of Condensed Matter Physics, 2010, 1, 323--345
2010
-
[17]
M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao and R. A. Simha, Rev. Mod. Phys., 2013, 85, 1143--1189
2013
-
[18]
Y. Fily, A. Baskaran and M. F. Hagan, Physical Review E, 2015, 91, 012125
2015
-
[19]
G. S. Redner, M. F. Hagan and A. Baskaran, Phys. Rev. Lett., 2013, 110, 055701
2013
-
[20]
M. E. Cates and J. Tailleur, Annual Review of Condensed Matter Physics, 2015, 6, 219--244
2015
-
[21]
Digregorio, D
P. Digregorio, D. Levis, A. Suma, L. F. Cugliandolo, G. Gonnella and I. Pagonabarraga, Physical review letters, 2018, 121, 098003
2018
-
[22]
C. B. Caporusso, P. Digregorio, D. Levis, L. F. Cugliandolo and G. Gonnella, Physical Review Letters, 2020, 125, 178004
2020
-
[23]
Sanoria, R
M. Sanoria, R. Chelakkot and A. Nandi, Soft Matter, 2024
2024
-
[24]
Das and R
S. Das and R. Chelakkot, Soft Matter, 2020, 16, 7250--7255
2020
-
[25]
S. Das, S. Ghosh and R. Chelakkot, Physical Review E, 2020, 102, 032619
2020
-
[26]
Kushwaha, S
P. Kushwaha, S. Maity, A. Menon, R. Chelakkot and V. Chikkadi, Soft Matter, 2024
2024
-
[27]
Sanoria, R
M. Sanoria, R. Chelakkot and A. Nandi, Phys. Rev. E, 2022, 106, 034605
2022
-
[28]
B \"a r, R
M. B \"a r, R. Gro mann, S. Heidenreich and F. Peruani, Annual Review of Condensed Matter Physics, 2020, 11, 441--466
2020
-
[29]
Sanchez, D
T. Sanchez, D. Welch, D. Nicastro and Z. Dogic, Science, 2011, 333, 456--459
2011
-
[30]
S. A. Yadav, D. Khatri, A. Soni, N. Khetan and C. A. Athale, Biophysical Journal, 2024, 123, 509--524
2024
-
[31]
Nitta, Y
T. Nitta, Y. Wang, Z. Du, K. Morishima and Y. Hiratsuka, Nature materials, 2021, 20, 1149--1155
2021
-
[32]
Xu, C.-r
T.-l. Xu, C.-r. Qin, B. Tang, J.-c. Gao, J. Zhou, K. Chen, T. H. Zhang and W.-d. Tian, The Journal of Chemical Physics, 2024, 161, 064905
2024
-
[33]
Zheng, M
E. Zheng, M. Brandenbourger, L. Robinet, P. Schall, E. Lerner and C. Coulais, Physical Review Letters, 2023, 130, 178202
2023
-
[34]
Tiwari, P
I. Tiwari, P. Parmananda and R. Chelakkot, Soft Matter, 2020, 16, 10334--10344
2020
-
[35]
De Canio, E
G. De Canio, E. Lauga and R. E. Goldstein, Journal of The Royal Society Interface, 2017, 14, 20170491
2017
-
[36]
Y. Fily, P. Subramanian, T. M. Schneider, R. Chelakkot and A. Gopinath, Journal of the Royal Society Interface, 2020, 17, 20190794
2020
-
[37]
Chelakkot, A
R. Chelakkot, A. Gopinath, L. Mahadevan and M. F. Hagan, Journal of The Royal Society Interface, 2014, 11, 20130884
2014
-
[38]
Elgeti, R
J. Elgeti, R. G. Winkler and G. Gompper, Reports on progress in physics, 2015, 78, 056601
2015
-
[39]
A. S. Sangani and A. Gopinath, Physical Review Fluids, 2020, 5, 083101
2020
-
[40]
Fatehiboroujeni, A
S. Fatehiboroujeni, A. Gopinath and S. Goyal, Journal of Computational and Nonlinear Dynamics, 2018, 13, 121005--121005--8
2018
-
[41]
Fatehiboroujeni, A
S. Fatehiboroujeni, A. Gopinath and S. Goyal, Physical Review E, 2021, 103, 013005
2021
-
[42]
X. Liao, P. K. Purohit and A. Gopinath, The Journal of Chemical Physics, 2020, 153, 194901
2020
-
[43]
R. G. Winkler, J. Elgeti and G. Gompper, Journal of the Physical Society of Japan, 2017, 86, 101014
2017
-
[44]
Eisenstecken, G
T. Eisenstecken, G. Gompper and R. Winkler, Polymers, 2016, 8, 304
2016
-
[45]
S. K. Anand and S. P. Singh, Physical Review E, 2018, 98, 042501
2018
-
[46]
S. K. Anand and l. P. Singh, Soft Matter, 2019, 15, 4008--4018
2019
-
[47]
Jayaraman, S
G. Jayaraman, S. Ramachandran, S. Ghose, A. Laskar, M. S. Bhamla, P. B. S. Kumar and R. Adhikari, Phys. Rev. Lett., 2012, 109, 158302
2012
-
[48]
Laskar, R
A. Laskar, R. Singh, S. Ghose, G. Jayaraman, P. B. S. Kumar and R. Adhikari, Scientific Reports, 2013, 3, 1964
2013
-
[49]
Chakrabarti and D
B. Chakrabarti and D. Saintillan, Physical Review Fluids, 2019, 4, 043102
2019
-
[50]
Elgeti and G
J. Elgeti and G. Gompper, Proceedings of the National Academy of Sciences, 2013, 110, 4470--4475
2013
-
[51]
Chakrabarti, S
B. Chakrabarti, S. F \"u rthauer and M. J. Shelley, Proceedings of the National Academy of Sciences, 2022, 119, e2113539119
2022
-
[52]
D. J. Hickey, R. Golestanian and A. Vilfan, Proceedings of the National Academy of Sciences, 2023, 120, e2307279120
2023
-
[53]
von Kenne, M
A. von Kenne, M. B \"a r and T. Niedermayer, Physical Review E, 2024, 109, 054407
2024
-
[54]
Chakrabarti and D
B. Chakrabarti and D. Saintillan, Physical review letters, 2019, 123, 208101
2019
-
[55]
W. Zhou, Z. Hao and N. Gravish, Physical Review X, 2021, 11, 031051
2021
-
[56]
W. Zhou, J. D. Peralta, Z. Hao and N. Gravish, Physical Review E, 2022, 105, 054604
2022
-
[57]
Quillen, A
A. Quillen, A. Peshkov, E. Wright and S. McGaffigan, Physical Review E, 2021, 104, 014412
2021
-
[58]
Chelakkot, M
R. Chelakkot, M. F. Hagan and A. Gopinath, Soft matter, 2021, 17, 1091--1104
2021
-
[59]
Gopinath and L
A. Gopinath and L. Mahadevan, Proceedings of the Royal Society A, 2011, 467, 1665--1685
2011
-
[60]
R. E. Isele-Holder, J. J \"a ger, G. Saggiorato, J. Elgeti and G. Gompper, Soft Matter, 2016, 12, 8495--8505
2016
-
[61]
F. Ling, H. Guo and E. Kanso, Journal of the Royal Society Interface, 2018, 15, 20180594
2018
-
[62]
Vilfan, S
A. Vilfan, S. Subramani, E. Bodenschatz, R. Golestanian and I. Guido, Nano Letters, 2019, 19, 3359--3363
2019
-
[63]
Laskar and R
A. Laskar and R. Adhikari, New Journal of Physics, 2017, 19, 033021
2017
-
[64]
A. L. Chau, C. D. Pugsley, M. E. Miyamoto, Y. Tang, C. D. Eisenbach, T. E. Mates, C. J. Hawker, M. T. Valentine and A. A. Pitenis, Tribology Letters, 2023, 71, 108
2023
-
[65]
deviated-bending
P. G. Gillespie and U. M \"u ller, Cell, 2009, 139, 33--44 mcitethebibliography Active-Clusters.out0000664000000000000000000000212214734304355013331 0ustar rootroot [1][-] section.1 Introduction [1][-] section.2 Model and simulation details [1][-] section.3 Simulation results ...
2009
-
[66]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.