{"id":"3a1ed3d1-e6e3-4ae5-b5df-be8cc303d905","arxiv_id":"2412.20536","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Active pinned-free filaments in a dense carpet assemble into kinetically arrested clusters whose size follows a scaling law in activity and spacing.","lead":"Simulations show that a row of wiggling, one-end-pinned elastic filaments can lock into regularly spaced, frozen clusters when they push hard enough and are spaced closely enough. The results suggest simple rules for how many filaments join a cluster and how the tower-like shape scales, which could guide designs for active colloid and microtubule materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The cluster-size scaling law is not tested against independent variation of elasticity or filament length, so the central quantitative claim rests on a single-parameter fit.","rationale":"The reader's verdict of CONDITIONAL is sound: the simulations plausibly demonstrate that pinned-free active filaments can form long-lived, clustered structures, and the corrugation mechanism (ℓ0 = σ) is explicitly identified as the source of effective tangential resistance. I do not see an internal contradiction that would justify rejection, and the absence of machine-checked proofs or released code is appropriately reflected in the conditional verdict. My stress-test concern is different in emphasis from the reader's weakest assumption: rather than questioning the arrest mechanism itself, I question whether the quantitative scaling law is actually supported by the data. The paper fixes κ, Nm, and ℓ in the main simulations, so the only varied inputs are f and δ. The dimensionless parameter β is a composite of f, κ, and Nm; claiming that cluster size is controlled by elasticity and filament length therefore requires an independent variation of κ or Nm, not just a dimensional rescaling. The missing derivation in §4.3 and the lack of error bars on the fitted exponents make this harder to check. A direct κ-sweep at fixed β would settle whether the scaling variable is sufficient; if it fails, the central predictive claim would need to be revised to a statement about f and δ only, while the phenomenological clustering results would remain. Since this concern reinforces the existing CONDITIONAL verdict rather than moving it, I recommend UNCHANGED.","tokens_in":112067,"tokens_out":6271,"duration_ms":78976,"concrete_test":"Run the identical NA = 300 pinned-free array at κ = 4×10^4 and κ = 1×10^4, adjusting f so that β = f σ^2 (Nm−1)^3/κ takes the same values used in Figure 4 (29.66, 44.49, 59.32), at fixed Δ ∈ {2, 3, 4}. If the measured ⟨Nf⟩ and cluster widths do not agree with the existing curve within the reported vertical-line fluctuations, then the scaling law does not collapse onto β, and the claimed dependence on elasticity fails. Independently, extract the full expression for ⟨Nf⟩ from ESM §A2 and verify that it contains no fitted exponent or undetermined prefactor.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative assertion is that the mean number of filaments in a cluster follows from a scaling theory in the activity parameter β = f σ^2 (Nm−1)^3/κ and the spacing Δ. In the provided text, however, the derivation of this scaling law is not displayed: §4.3 references Figure 8 and an ESM analysis, but no closed-form expression for ⟨Nf⟩ is given, and the 'prediction' in Figure 4(a) is not accompanied by an equation, prefactors, or error analysis. More importantly, the simulation campaign varies only f and δ, with κ = 2×10^4, Nm = 40, and ℓ = 39 held fixed throughout the main runs. The scaling variable β therefore encodes the elasticity and length dependence only through a dimensional combination; that dependence is never checked by an independent sweep over κ or Nm. If a simulation at doubled κ with correspondingly doubled f (keeping β fixed) produced a different cluster size, the claim that cluster size is controlled by elasticity—rather than merely by the force magnitude f—would fail, even though the clustering phenomenon itself would still occur. The fitted self-similar exponents a(β) ≈ 1.6, 1.3, 1.05 in §3.4 are reported without uncertainties, adding to the risk that the quantitative collapse is in-sample. This is a correctness risk for the predictive part of the central claim, not merely a presentation gap.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":112332,"tokens_out":4766,"duration_ms":55552,"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":[{"comment":"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.","section":"§4.3, Fig. 4(a)"},{"comment":"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.","section":"Table 1 and §4.3"},{"comment":"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.","section":"§3.4, Fig. 3(d)–(f)"},{"comment":"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.","section":"§3.5 and Fig. 4"}],"minor_comments":[{"comment":"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.","section":"Eq. (4)"},{"comment":"The caption lists β = 44.94 while the text and other figures use β = 44.49; please check which value is correct.","section":"Fig. 5 caption"},{"comment":"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.","section":"Fig. 4(b)"},{"comment":"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.","section":"§3.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript presents a visually robust simulation phenomenon and a plausible mechanistic picture. The main obstacle is that the central scaling law is not actually displayed in the main text and is not tested against independent variation of κ or Nm. If the authors can supply the explicit prediction, its derivation, and at least a few such validation runs, the paper would be a solid candidate for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth knowing: this is a competent Brownian dynamics study of pinned-free active filament carpets, and the qualitative result—that corrugation-induced tangential contact forces produce kinetically arrested clusters—is credible and visually robust. The novelty is real but incremental: their earlier work flagged jammed clusters for pinned-free arrays; this manuscript adds systematic sweeps over f and δ, the self-similar collapse of cluster shapes with activity-dependent exponents, and a scaling expression for mean cluster size. The pressure-field analysis is a nice extra and strengthens the mechanistic story.\n\nWhat is new: cluster size decreases monotonically with Δ and β; the width collapse xΔ^a(β) works reasonably; pressure concentrates at the cluster waist. These are concrete, potentially design-relevant results.\n\nSoft spots: the scaling theory for ⟨Nf⟩ is referenced but not actually derived in the main text; Fig 4(a) compares to a dashed line without a closed-form expression, prefactor, or error estimate. The self-similar exponents a ≈ 1.6, 1.3, 1.05 are best fits with no uncertainties and visibly imperfect collapse near the pinned ends. More importantly, all main runs use κ = 2×10^4, Nm = 40, ℓ = 39 fixed; only f and δ are varied. So β encodes elasticity and length through a dimensional combination, but that dependence is never tested independently. The claim that cluster size is controlled by elasticity rather than just force magnitude would fail if a κ sweep at fixed β changed ⟨Nf⟩. Also, the arrest mechanism relies on ℓ0 = σ making filament surfaces corrugated; the paper states this but does not explore smoother fibers, so the design principle is conditional on one microscopic choice.\n\nNone of this breaks the central qualitative claim. The clustering phenomenon is robust and the simulations appear honestly described. But the quantitative centerpiece needs a displayed derivation, error bars on exponents, and preferably a κ or filament-length sweep at fixed β. This is a solid candidate for peer review, not a desk reject, and not accept yet.\n\nWho this is for: active-matter and soft-matter readers, especially synthetic cilia and microtubule-motor communities. I would bring it to a reading group and would cite the clustering phenomenology, though not the scaling law until it is tested.","headline":"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.","tokens_in":112876,"tokens_out":1688,"would_cite":true,"duration_ms":22562,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["active filaments","kinetically arrested clusters","follower forces","Brownian dynamics","filament arrays","self-assembly","steric friction","activity-elasticity scaling"],"falsifier":"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.","tokens_in":111852,"feed_emoji":"🧵","tokens_out":7079,"duration_ms":68483,"temperature":0.7,"pith_summary":"The paper claims that a fluid-permeated array of elastic filaments, pinned at one end and free at the other, spontaneously organizes into regularly spaced, compact clusters that are kinetically arrested rather than equilibrium structures. The clustering is driven by a global buckling instability of each active filament, and it stops when steric contact forces plus effective surface friction balance the active compressive forces and bending elasticity. The paper derives a scaling theory in which the mean number of filaments in a cluster and the spacing between clusters are set by two dimensionless parameters: activity $\\beta = f\\sigma^2(N_m-1)^3/\\kappa$ and anchoring gap $\\Delta = \\delta/\\sigma$. If correct, cluster size is not a materials accident but a predictable function of activity, elasticity, and grafting density, which gives a design rule for self-assembled filament fabrics.","feed_headline":"Active filaments self-assemble into regular, jammed clusters","feed_subtitle":"Cluster size and spacing follow a simple activity-elasticity scaling, giving design rules for self-assembling filament materials.","key_machinery":"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.","core_discovery":"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$.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Shows that a straight pinned-free filament under compressive follower forces becomes unstable via a global bifurcation to coiled-buckled, rotating shapes; this is the single-filament instability the clusters are built from.","marker":"36"},{"why":"Provides the bead-chain Brownian dynamics model and its stochastic extension that the present simulations adapt, including the nonlinear coiled states used to interpret clustering.","marker":"37"},{"why":"Earlier active-filament-array study identifying the ratio $\\ell_0/\\sigma$ as an effective tangential friction parameter and showing that clamped versus pinned boundary conditions switch between waves and jammed clusters.","marker":"58"},{"why":"Continuum modeling showing that local hydrodynamic drag captures collective single- and multi-filament dynamics, used here to justify the freely draining approximation.","marker":"39"},{"why":"Analytical treatment of distributed follower-force filaments that supplies the instability threshold and rotating solution branch used to interpret the onset of clustering.","marker":"40"},{"why":"Analytical study of follower-force elastic filaments whose critical $\\beta$ and nonlinear beating or rotation results set the activity baseline for the clustering regime.","marker":"61"},{"why":"Experiments on electrically actuated colloidal chains that display the same buckling and rotation as the model, supporting the follower-force description of active filaments.","marker":"9"},{"why":"Experiments on motor-driven filament assemblies showing flagella-like beating and coiling, providing the synthetic and biological context for the arrested clusters.","marker":"30"}],"fun_headline_variants":["Filament arrays jam into periodic clusters","Active filaments buckle into arrested towers","No attractions needed: active filaments cluster","Filament activity drives self-assembled clusters","Buckling filaments form regular jammed clusters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Filament arrays jam into periodic clusters","Active filaments buckle into arrested towers","No attractions needed: active filaments cluster","Filament activity drives self-assembled clusters","Buckling filaments form regular jammed clusters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00016,"raw_usage":{"total_tokens":1221,"prompt_tokens":922,"completion_tokens":299,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":538,"completion_tokens_details":{"reasoning_tokens":236}},"tokens_in":538,"tokens_out":299,"duration_ms":3552,"temperature":1.0,"reasoning_tokens":236,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T23:18:44.046239+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that a straight pinned-free filament under compressive follower forces becomes unstable via a global bifurcation to coiled-buckled, rotating shapes; this is the single-filament instability the clusters are built from."},{"cited_title":"Chelakkot, A","cited_arxiv_id":null,"evidence_quote":"Provides the bead-chain Brownian dynamics model and its stochastic extension that the present simulations adapt, including the nonlinear coiled states used to interpret clustering."},{"cited_title":"Chelakkot, M","cited_arxiv_id":null,"evidence_quote":"Earlier active-filament-array study identifying the ratio $\\ell_0/\\sigma$ as an effective tangential friction parameter and showing that clamped versus pinned boundary conditions switch between waves and jammed clusters."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Continuum modeling showing that local hydrodynamic drag captures collective single- and multi-filament dynamics, used here to justify the freely draining approximation."},{"cited_title":"Fatehiboroujeni, A","cited_arxiv_id":null,"evidence_quote":"Analytical treatment of distributed follower-force filaments that supplies the instability threshold and rotating solution branch used to interpret the onset of clustering."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Analytical study of follower-force elastic filaments whose critical $\\beta$ and nonlinear beating or rotation results set the activity baseline for the clustering regime."},{"cited_title":"Nishiguchi, J","cited_arxiv_id":null,"evidence_quote":"Experiments on electrically actuated colloidal chains that display the same buckling and rotation as the model, supporting the follower-force description of active filaments."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experiments on motor-driven filament assemblies showing flagella-like beating and coiling, providing the synthetic and biological context for the arrested clusters."}],"review_version":1}