{"id":"51151d18-3b91-4635-9c99-9546cf0969b4","arxiv_id":"2412.14419","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Cluster stability in photoactive granular particles is set by the balance of adsorption and desorption, captured by κ = const·(D/v)/N_T^{3/2}, with stable clusters only when κ < 0.3849 and initial size exceeds a critical value.","lead":"Experiments with light-powered centimeter-scale robots show that clusters form and persist at low light and high robot numbers, while high light and low numbers dissolve them. A simple adsorption-desorption model with one control parameter reproduces the phase diagram and predicts a critical cluster size for survival.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Zero-order desorption in Eq. (S2) is the load-bearing unmeasured assumption; the predicted ncrit=1/3 and the bistable phase structure depend on it, although the dissolution experiment partially corroborates it.","rationale":"The reader identified the zero-order desorption assumption as the weakest point, and I agree it is the most load-bearing condition for the central claim: without a constant desorption term, the model loses the lower unstable fixed point that produces the predicted critical initial cluster size ncrit and the specific value ncrit=1/3 at κ=κcrit. However, I also recognize that the dissolution experiment in Fig. 5(e) provides nontrivial, parameter-free support for this assumption: the observed threshold near n0≈0.33 is exactly what the zero-order model predicts at the stability boundary, and it would be absent under pure perimeter or first-order desorption. This independent check is the paper's strongest evidence and prevents the concern from being fatal. The residual doubt is that the evidence is indirect—no direct measurement of kd(N) is reported—and a mixed-order desorption term could shift ncrit while preserving a threshold-like feature. This is an addressable, empirical gap rather than an internal inconsistency. Therefore the conditional verdict is appropriate, and the concrete test described would either retire the concern or force a re-derivation of the model's bistable phase diagram. I did not identify a stronger objection: the growth-fit halving of NT is a post-hoc rationalization but is not essential to the phase diagram or the dissolution test, and the calibrated P↔D/v mapping is partly validated by the predicted D/v value at the NT=80 boundary. No ad hominem is intended; the critique is purely about the empirical grounding of the desorption term.","tokens_in":13342,"tokens_out":6949,"duration_ms":61026,"concrete_test":"Repeat the dissolution protocol of Fig. 5 at the same highest illumination, but prepare clusters over a range of initial sizes, e.g., n0 ≈ 0.1, 0.2, 0.3, 0.4, 0.5, by switching the illumination at different growth times. At the switch time t0, measure the initial dissolution slope dn/dt at t0+ for each n0. Under the zero-order model, after subtracting the independently estimated adsorption contribution using the measured free-particle velocity and cluster geometry, the inferred kd should be independent of n0. Then fit the generalized desorption kd(N) = a + b N^{1/2} + c N to the same data. If the b or c terms are statistically significant, or if a varies systematically with illumination intensity, the zero-order assumption is falsified and the predicted ncrit must be re-derived.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative predictions of the kinetic model—the existence of a critical initial cluster size ncrit and the value ncrit=1/3 at the stability boundary—follow directly from the zero-order desorption term kd = DΦc/Ap in Eq. (S2). If desorption instead scales with the number of surface particles (perimeter, proportional to N^{1/2}) or with cluster size N, the balance equation dn/dτ = (1−n)n^{1/2} − κ is replaced by a different form. Pure perimeter desorption gives dn/dτ = n^{1/2}[(1−n)−κ'], which has no ncrit threshold: any positive cluster grows whenever κ'<1. Pure first-order desorption gives dn/dτ = (1−n)n^{1/2} − κ''n, again without the lower unstable fixed point that creates ncrit. The measured dissolution threshold near n0≈0.33 in Fig. 5(e) rules out these pure alternatives and is therefore meaningful evidence for the zero-order assumption. However, the support is indirect: a mixed desorption rate kd(N) = a + b N^{1/2} + c N with small b or c would shift ncrit continuously, and finite-size or cluster-geometry effects could mimic the threshold. The paper does not directly measure the detachment rate as a function of N, so the exact functional form of kd(N) remains an assumption. Because this assumption is what generates the bistable structure and the ncrit=1/3 prediction, it is the most load-bearing element of the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies cluster formation and dissolution in a macroscopic photoactive granular system: modified HEXBUG particles powered by photovoltaic cells under programmable illumination. For low light intensity (lower particle activity) and larger population sizes, clusters are stable; for high intensity and small populations, clusters dissolve. The authors measure cluster-duration survival functions, extract a power-law exponent alpha, and build a phase diagram in the (intensity, population-size) plane. They introduce a mean-field kinetic model in which a cluster grows by adsorption of free particles and shrinks by a zero-order desorption rate, giving dn/dtau = (1-n)n^{1/2} - kappa. The model predicts a critical kappa_crit and a critical initial cluster size ncrit. A dedicated experiment with a light-intensity switch is used to test the predicted growth law and the existence of a threshold initial size for dissolution around n0 ~ 0.33.","tokens_in":13745,"tokens_out":6746,"duration_ms":53830,"significance":"If the claims hold, the paper provides a simple, analytically tractable model for cluster stability in inertial dry active matter, together with a versatile experimental platform (light-controlled macroscopic particles) and a dedicated test of a nontrivial prediction (critical initial cluster size). The authors openly provide data and code, repeat experiments five-fold, and include robustness checks with different wall geometries. The model is not derived from the target data, and the dissolution experiment partially corroborates the central zero-order desorption assumption. The critical size prediction and phase diagram are meaningful qualitative results, but the quantitative comparison is weakened by a calibrated P-to-D/v mapping and by the absence of error bars on the experimental phase boundary.","major_comments":[{"comment":"The zero-order desorption rate kd = D*Phi_c/A_p is the load-bearing assumption that generates the lower unstable fixed point and the critical initial cluster size ncrit. The paper does not directly measure kd as a function of cluster size N. The dissolution experiment in Fig. 5(e) provides indirect support, but alternative desorption forms kd(N) = a + b*N^{1/2} + c*N would shift or eliminate ncrit. The authors should either measure the detachment rate for clusters of controlled size or provide a robustness analysis showing that the predicted threshold is stable to plausible perimeter- and size-dependent desorption terms.","section":"Kinetic model, SM Eq. (S2)"},{"comment":"The experimental phase diagram rests on the fitted power-law exponent alpha, but the paper reports no error bars, no fitting range, and no goodness-of-fit for the alpha values. The dashed boundary in Fig. 3(c) is drawn by eye, and the transition band alpha = 2 +/- 0.1 is arbitrary. Since the phase boundary is the central quantity the model is claimed to reproduce, the authors should provide confidence intervals for alpha (e.g., via bootstrap or by fitting over several t-ranges) and state a quantitative criterion for assigning stable/unstable phases.","section":"Results, Fig. 3(a-c)"},{"comment":"The mapping between the experimental light power P and the model parameter D/v is calibrated from the phase boundaries; the text states that 'based on the phase boundaries, one can establish a relation between these variables.' The agreement between Fig. 3(c) and Fig. 4(b) is therefore partly a fit rather than a parameter-free prediction. The model still makes a nontrivial qualitative prediction through the NT-dependence, but the claim of reproducing the phase space should be qualified, or the calibration should be anchored by independent measurements of D and v from the individual-particle characterization in SM Sec. V.","section":"Kinetic model, Fig. 4(b)"},{"comment":"The dashed line in Fig. 5(e) is called the prediction of the kinetic model for the critical cluster size, evaluated at 'the highest value of D/v that allows the development of stable clusters.' This choice sets kappa = kappa_crit, for which the model's lower stable size is ncrit = 1/3 by construction. The data showing a threshold near 0.33 are consistent, but this is not an independent test of the predicted numeric value. An independent estimate of D/v at the highest intensity, or a distribution of predicted thresholds from particle-level measurements, would make the comparison meaningful.","section":"Dissolution experiment, Fig. 5(e)"}],"minor_comments":[{"comment":"The expression for the association phase, nassoc(t) = (NT/N) tanh^2(1/2 sqrt(NT) (Ct + c)), appears to be a typo: the solution of dn/dtau = (1-n)n^{1/2} in rescaled time is n(t) = tanh^2(1/2 sqrt(NT) (Ct + c)) up to a constant determined by initial conditions. If the factor NT/N is intentional, the derivation should be shown.","section":"Kinetic model, near Eq. (3)"},{"comment":"The caption refers to 'the phase diagram obtained with the model [Fig. 4(c)]', but Fig. 4 contains only panels (a) and (b); the intended reference appears to be Fig. 4(b).","section":"Caption of Fig. 3(c)"},{"comment":"The text uses 'absorption rate' for ka, while 'adsorption' is the standard term for particle attachment to a cluster; the terminology should be made consistent.","section":"SM Sec. III"},{"comment":"The fitted power-law exponent alpha is not shown on the survival-function plots; adding fitted lines and the obtained alpha values for each curve would help the reader judge the quality of the power-law description.","section":"Fig. 3(a,b)"},{"comment":"The wall-geometry comparison changes both the wall structure and the population size NT between panels; state explicitly that this is a qualitative robustness check rather than a controlled comparison of wall geometry alone.","section":"SM Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The manuscript describes well-executed experiments with open data and code, and the model is transparent. The main risk is overstatement of predictive power: the P-to-D/v mapping is calibrated to the phase boundary, and the zero-order desorption rate that produces the ncrit prediction is not directly measured. The dissolution experiment is a meaningful partial test, but the threshold comparison is partly circular because the parameters are chosen at kappa_crit. With a robustness analysis of the desorption form and quantitative error bars on the experimental phase diagram, the central claims would be substantially stronger. The topic fits the scope of the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. This is a genuinely useful experimental paper, not a revolution. The new thing is the centimeter-scale photoactive particle platform with programmable light, and the phase diagram of cluster stability as a function of population size and activity. The kinetic model is a close cousin of earlier adsorption-desorption treatments (Peruani et al., Ginot et al.), but the paper adds a nontrivial prediction: a critical initial cluster size ncrit below which clusters dissolve, then tests it directly with a controlled experiment. That dissolution experiment is the best part of the paper. The threshold near n0≈0.33 matches the model's prediction and, as the stress-test note says, it rules out pure perimeter or first-order desorption alternatives. That is real evidence, not just curve fitting.\n\nThe soft spots are real but not fatal. The α exponents in the phase diagram are reported without error bars, and the phase boundary in Fig. 3(c) is a guide to the eye. More importantly, the P↔D/v mapping is calibrated from the phase boundaries themselves, so the model's 'reproduction' of the experimental phase diagram is partly circular. The growth-curve fit also needs to halve NT to work, which the authors explain by two-cluster formation but is still a post-hoc parameter. And the zero-order desorption rate, Eq. S2, is load-bearing: if desorption scaled with perimeter or cluster size, the ncrit = 1/3 prediction would shift. The dissolution experiment gives indirect support, but the paper never directly measures kd as a function of N. These are the questions a referee should push on.\n\nNet: this paper is worth publishing after revision. The experimental platform is new, the data are open, and the critical-size prediction is falsifiable and partially confirmed. I would send it to peer review and ask for error bars on α, a clearer statement of which aspects of the model comparison are calibrated, and a direct or at least more systematic test of the desorption scaling. For a soft matter / active granular audience, it's a solid contribution.","headline":"Solid experimental Letter with a genuine new prediction—critical initial cluster size—supported by a clean dissolution experiment; the main caveats are calibrated model mapping and an unmeasured zero-order desorption assumption, both addressable in revision.","tokens_in":14199,"tokens_out":1808,"would_cite":true,"duration_ms":15246,"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":"This paper claims that cluster stability in macroscopic photoactive particles is set by a single kinetic control parameter κ, with clusters dissolving above κcrit = 2√3/9 ≈ 0.385 or below a critical initial size.","keywords":["photoactive granular matter","active matter","cluster dynamics","adsorption-desorption model","cluster stability phase diagram","light-controlled activity","self-propelled particles"],"falsifier":"Measure the per-particle detachment rate for clusters of controlled size at fixed illumination: the model requires that this rate be independent of $N$, so a clear dependence on cluster size, such as scaling like $N^{-1/2}$, would refute the zero-order desorption assumption. Alternatively, prepare initial clusters with sizes straddling the predicted $n_0 \\approx 0.33$ under the highest activity that still permits stable clusters and check that the dissolution probability drops sharply at that size.","tokens_in":13192,"feed_emoji":"🔆","tokens_out":8739,"duration_ms":72396,"temperature":0.7,"pith_summary":"This paper reports experiments on centimeter-scale photoactive particles that move when illuminated, and proposes a minimal kinetic explanation for when their clusters survive. The authors find that stable clustering occurs only at low particle activity and large population sizes, while high activity or small populations produce clusters that form and dissolve rapidly. They reduce the cluster dynamics to a single adsorption–desorption equation, $dn/d\\tau = (1-n)\\sqrt{n} - \\kappa$, with $\\kappa$ combining geometry, population, and the ratio of diffusion to velocity. The model predicts a sharp threshold $\\kappa_{\\mathrm{crit}} = 2\\sqrt{3}/9 \\approx 0.385$ above which no cluster is stable, plus a critical initial cluster size below which even stable-phase clusters dissolve; controlled experiments with light-switched activity corroborate both features. If correct, this gives a simple design rule for programmable assembly and disassembly in active granular matter.","feed_headline":"A one-line kinetic rule predicts which clusters survive","feed_subtitle":"Cluster stability in light-driven particles hinges on one parameter κ with a critical value and a minimum cluster size.","key_machinery":"The central object is the dimensionless adsorption–desorption balance $dn/d\\tau = (1-n)\\sqrt{n} - \\kappa$. The adsorption term $(1-n)\\sqrt{n}$ encodes the collision rate of free particles with a cluster: the gas density contributes $(1-n)$ and the cluster perimeter contributes $\\sqrt{n}$, so larger clusters capture particles faster. The desorption term $\\kappa$ is zero-order, treated as independent of cluster size. The balance curve $\\kappa = (1-n)\\sqrt{n}$ has a single maximum at $\\kappa_{\\mathrm{crit}} = 2\\sqrt{3}/9$, which separates a regime where clusters can stabilize from one where every cluster decays, and the lower branch of the same curve fixes the minimum seed size needed for stable growth. This one-equation framework carries the entire phase diagram.","core_discovery":"The authors establish a transition between unstable and stable clustering in a macroscopic active-matter system. By analyzing the survival functions of cluster durations, they find a power-law exponent $\\alpha > 2$ (finite mean duration, unstable clusters) at high activity and small populations, and $\\alpha \\le 2$ (diverging mean duration, stable clusters) at low activity and large populations. They then derive a deterministic kinetic model, $dN/dt = v(N_T-N)/A_T\\sqrt{NA_p/\\Phi_c} - D\\Phi_c/A_p$, which rescales to $dn/d\\tau = (1-n)\\sqrt{n} - \\kappa$ with $\\kappa = (\\Phi_c^{3/2}A_T)/(A_p^{3/2}N_T^{3/2}) \\cdot (D/v)$. The model predicts that above $\\kappa_{\\mathrm{crit}} = 2\\sqrt{3}/9 \\approx 0.3849$ all clusters disassociate, while below it a cluster must start above a size-dependent critical value $n_{\\mathrm{crit}}(\\kappa)$ to grow and stabilize. Controlled experiments, in which half of the arena is dark to grow a cluster and then the whole arena is illuminated to dissolve it, recover the predicted $\\tanh^2$ growth curve and show the dissociation probability starting to fall near the predicted initial size $n_0 \\approx 0.33$.","pith_inferences":["If the zero-order desorption term is the right caricature, the same model should transfer to other contact-driven active agents, such as vibrated granular matter or robot swarms, by reinterpreting $D$, $v$, and $\\Phi_c$; the predicted critical value $\\kappa_{\\mathrm{crit}}$ would then be a universal number for this class of systems.","Because $\\kappa$ factors into geometry, population, and activity, an engineering prescription follows: choose the arena size and particle number so the system sits below $\\kappa_{\\mathrm{crit}}$, then use light gradients to seed clusters larger than $n_{\\mathrm{crit}}$; this turns cluster stability into a programmable switch.","The observation that boundary clustering persists at low activity even with inward-guiding walls points to a quantity the model does not yet contain: a wall-trapping time that could be measured by tracking single-particle residence times near the boundary, potentially extending the kinetic equation to spatially structured light fields."],"forward_implications":["With the same arena and particles, lowering $D/v$ (lower light intensity) or raising the population $N_T$ moves the system below the critical $\\kappa$, turning transient clustering into persistent clustering.","A cluster that starts smaller than $n_{\\mathrm{crit}}(\\kappa)$ will dissolve even under conditions that otherwise support stable clustering, so assembly protocols must seed clusters above this threshold.","For $\\kappa < \\kappa_{\\mathrm{crit}}$, the model's nullcline $\\kappa=(1-n)\\sqrt{n}$ determines the final stable cluster size reached as $\\tau \\to \\infty$, giving a quantitative prediction for the long-time cluster population.","Switching illumination from low to high should reverse cluster growth along the dissociation branch of the kinetic equation, a prediction the authors confirm in controlled $N_T=100$ experiments.","The cluster-duration statistics should shift from a power law with exponent $\\alpha>2$ to $\\alpha\\le 2$ exactly where the control parameter crosses $\\kappa_{\\mathrm{crit}}$, linking the measured survival functions to the phase boundary."],"supporting_citations":[{"why":"Supplies the cluster-dynamics kinetic framework of adsorption–desorption balances that the model adapts to photoactive particles.","marker":"[51]"},{"why":"Provides the aggregation–fragmentation and individual-dynamics treatment of active clusters that motivates the desorption term.","marker":"[52]"},{"why":"Documents wall accumulation of self-propelled particles, the seed mechanism for boundary clustering that the flower-shaped walls are designed to counteract.","marker":"[43]"},{"why":"Gives the rationale for the inward-guiding wall design used to prevent boundary clustering in the main arena.","marker":"[49]"}],"fun_headline_variants":["Light controls cluster fate via one critical parameter","Kinetic rule predicts which clusters survive in light","One parameter decides cluster stability in active matter","Critical kappa: the switch for cluster survival in light"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that particles leave a cluster at a constant rate independent of cluster size, so if small clusters actually lose surface particles faster than large ones the predicted stable phase and critical initial size would shift.","fun_headline_variants_meta":{"raw":{"variants":["Light controls cluster fate via one critical parameter","Kinetic rule predicts which clusters survive in light","One parameter decides cluster stability in active matter","Critical kappa: the switch for cluster survival in light"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000479,"raw_usage":{"total_tokens":2372,"prompt_tokens":949,"completion_tokens":1423,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":1364}},"tokens_in":565,"tokens_out":1423,"duration_ms":9905,"temperature":1.0,"reasoning_tokens":1364,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:15:40.863355+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the per-particle detachment rate for clusters of controlled size at fixed illumination: the model requires that this rate be independent of $N$, so a clear dependence on cluster size, such as scaling like $N^{-1/2}$, would refute the zero-order desorption assumption. Alternatively, prepare initial clusters with sizes straddling the predicted $n_0 \\approx 0.33$ under the highest activity that still permits stable clusters and check that the dissolution probability drops sharply at that size.","supporting_citations":[{"cited_title":"Peruani, L","cited_arxiv_id":null,"evidence_quote":"Supplies the cluster-dynamics kinetic framework of adsorption–desorption balances that the model adapts to photoactive particles."},{"cited_title":"Ginot, I","cited_arxiv_id":null,"evidence_quote":"Provides the aggregation–fragmentation and individual-dynamics treatment of active clusters that motivates the desorption term."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents wall accumulation of self-propelled particles, the seed mechanism for boundary clustering that the flower-shaped walls are designed to counteract."},{"cited_title":"Kumar, H","cited_arxiv_id":null,"evidence_quote":"Gives the rationale for the inward-guiding wall design used to prevent boundary clustering in the main arena."}],"review_version":1}