arxiv: 2605.04812 · v1 · submitted 2026-05-06 · ❄️ cond-mat.stat-mech

Recognition: unknown

Role of mass fluctuations in the diffusion of clusters of Brownian particles with activity

Antonio Suma, Daniela Moretti, Giuseppe Gonnella, Pasquale Digregorio

Authors on Pith no claims yet

Pith reviewed 2026-05-08 16:36 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech

keywords active Brownian particlescluster diffusionmass fluctuationsanomalous diffusionLangevin equationscenter of massGaussian process

0 comments

The pith

Mass fluctuations in active Brownian particle clusters produce anomalous center-of-mass diffusion scaling as N to the minus 0.63

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The authors start from the Langevin equations for individual active particles and derive two coupled stochastic equations: one tracking the cluster's center-of-mass motion and one for the fluctuating number of particles N(t) in the cluster. They model N(t) as a Gaussian process with mean N0, variance scaling as N0 to the beta, and persistence time scaling as N0 to the kappa. Solving the equations analytically shows the long-time diffusion coefficient D is the sum of a standard term scaling as N0 to the minus one from noise and active forces, plus a fluctuation term scaling as N0 to the minus delta, where delta equals 2 minus 2 over d minus beta plus kappa. When the fluctuation term dominates, this yields the observed anomalous scaling D approximately N to the minus 0.63 that matches large-scale simulations.

Core claim

The long-time diffusion coefficient of the cluster center of mass is the sum of a conventional contribution proportional to N0 inverse from thermal noise and summed active forces, and a fluctuation-driven contribution proportional to N0 to the power -delta with delta = 2 - 2/d - beta + kappa. This leads to D scaling as N to the power -alpha with alpha equal to 0.63 plus or minus 0.06 in quantitative agreement with ABP simulations of isolated clusters.

What carries the argument

Coupled stochastic equations for the center-of-mass trajectory and the cluster mass N(t), with N(t) treated as an independent Gaussian process.

If this is right

The diffusion coefficient separates into conventional and fluctuation-driven parts.
Anomalous scaling D ~ N^{-α} occurs when the fluctuation term is dominant.
The value of the exponent α is fixed by the mass fluctuation exponents β and κ together with dimension d.
The analytical predictions match the scaling observed in large-scale simulations of active Brownian particle clusters.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same fluctuation mechanism could explain anomalous diffusion in other systems with variable particle numbers, such as biological aggregates or colloidal assemblies.
Including possible correlations or back-action between the center-of-mass velocity and the attachment/detachment rates might modify the predicted exponent.
The framework might be generalized to clusters in non-stationary states or under external fields.
Experimental tests could involve tracking cluster sizes and diffusion in active colloidal systems with tunable activity.

Load-bearing premise

The cluster size N(t) can be modeled as an independent Gaussian process whose statistical properties are fixed independently of the center-of-mass dynamics.

What would settle it

Observation in simulations or experiments that the variance or persistence time of N(t) depends on the instantaneous velocity of the cluster center of mass would contradict the independence assumption.

Figures

Figures reproduced from arXiv: 2605.04812 by Antonio Suma, Daniela Moretti, Giuseppe Gonnella, Pasquale Digregorio.

**Figure 1.** Figure 1: (a) Time evolution of a stationary cluster of view at source ↗

**Figure 2.** Figure 2: (a) Schematic of fragment attachment to the main cluster. The fragment (center-of-mass view at source ↗

**Figure 3.** Figure 3: Diffusion coefficient D versus mean cluster size N0, as given by Eq. (18), using Da = 8.33, d = 2 (same value as ABPs simulations, see SM). (a) Fixed Dg = 1.4 (from fitted ABPs simulation values); and different values of δ. (b) Fixed δ = 0.63 and different values of Dg. Da = 8.33 (SM), while the fits in view at source ↗

read the original abstract

Motivated by the anomalous diffusion observed in clusters of active Brownian particles (ABPs), where the center-of-mass diffusion coefficient scales as $D\sim N^{-1/2}$ with respect to the number $N$ of particles in the cluster, we derive a minimal theoretical framework starting from the single-particle Langevin equations. The model consists of two coupled stochastic equations: one for the cluster center-of-mass trajectory and one for the mass evolution $N(t)$, explicitly accounting for stochastic displacements induced by particle attachment and detachment. We specialize and validate the framework against ABP simulations of isolated clusters in stationary conditions, where $N(t)$ follows a Gaussian process with mean $N_0$, variance $\propto N_0^\beta$, and persistence time $\propto N_0^\kappa$. Analytical solution of the coupled equations yields the long-time diffusion coefficient as the sum of two contributions: a conventional term $\propto N_0^{-1}$) due to thermal noise plus summation of active forces, and a fluctuation-driven term $\propto N_0^{-\delta}$ with $\delta=2-2/d-\beta+\kappa$, where $d$ is the spatial dimension. We demonstrate that anomalous scaling emerges whenever the second term becomes dominant. The model predicts $D\sim N^{-\alpha}$ with $\alpha=0.63\pm0.06$, in good quantitative agreement with large-scale ABP simulations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Mass fluctuations drive a derived extra term in cluster diffusion for ABPs, but the scaling exponent comes from fitting the same simulations used for validation.

read the letter

The main thing here is a derivation that ties mass fluctuations in active Brownian particle clusters to their center-of-mass diffusion. Starting from single-particle Langevin equations, they write coupled stochastic equations for the cluster position and for N(t), the instantaneous number of particles. Treating N(t) as a Gaussian process with mean N0, variance scaling as N0^β, and persistence as N0^κ, they solve for the long-time diffusion constant. It comes out as the sum of the usual thermal-plus-active term proportional to 1/N0 and a new term from the fluctuations that goes as N0 to the power of minus (2 - 2/d - β + κ). When the second term wins, you get anomalous scaling, and plugging in their measured β and κ gives α around 0.63, which lines up with the ABP runs. The derivation itself is straightforward and internally consistent. They show how the velocity autocorrelation picks up the extra contribution from the fluctuating force due to changing N. That part is new and gives a concrete expression for the fluctuation-driven diffusion. The limitation is that β and κ are extracted from the very simulations they compare against, so the match on the exponent is not an independent test. The model also assumes that N(t) evolves independently of the cluster's velocity, with no back-action from motion on attachment or detachment. If the cluster's movement affects how particles join or leave, that could introduce correlations that change the effective scaling. The paper does not check for those correlations explicitly. This work is aimed at the active matter and statistical mechanics community, particularly those studying clustering and diffusion in driven particle systems. It offers a simple framework that could be extended to other active systems where cluster size fluctuates. I think it deserves peer review because the approach is clear and the idea is worth testing further, even with the caveats on the fitted parameters and independence assumption.

Referee Report

3 major / 2 minor

Summary. The manuscript develops a minimal model for the center-of-mass diffusion of clusters of active Brownian particles (ABPs) by coupling the cluster trajectory to stochastic mass fluctuations N(t). Starting from single-particle Langevin equations, it derives two coupled SDEs, specializes to the case where N(t) is an exogenous Gaussian process with mean N0, variance ~N0^β and persistence ~N0^κ, and obtains an analytical long-time diffusion coefficient D as the sum of a conventional term ~N0^{-1} and a fluctuation-driven term ~N0^{-δ} with δ=2-2/d-β+κ. The model predicts anomalous scaling D~N^{-α} with α=0.63±0.06 when the second term dominates, in agreement with ABP simulations.

Significance. If the independence assumption for N(t) holds, the framework supplies a mechanistic account of how mass fluctuations can generate anomalous diffusion in active clusters, separating thermal/active-force contributions from fluctuation-driven ones. The analytical solution of the coupled SDEs is a technical strength that could generalize to other fluctuating-particle systems. However, because β and κ are measured from the same simulations used for validation, the quantitative match is a consistency check rather than an independent prediction, limiting the result's first-principles impact.

major comments (3)

[Abstract / model setup] Abstract and model derivation: the analytical expression for the fluctuation-driven term ∝ N0^{-δ} (with δ=2-2/d-β+κ) is obtained only by treating N(t) as an independent exogenous Gaussian process whose statistics are inserted directly into the center-of-mass equation. The manuscript provides no demonstration that cross-correlations between δN(t) and cluster velocity (arising from attachment-induced displacements or motion-dependent attachment rates) remain negligible once the equations are fully coupled; any such correlations would renormalize the velocity autocorrelation integral and change the predicted scaling.
[Abstract] Abstract and validation paragraph: the claim that the model 'predicts' α=0.63±0.06 'in good quantitative agreement' with ABP simulations is circular because the exponents β and κ that enter δ are extracted from the identical simulations. This reduces the numerical match to a post-hoc consistency test rather than an a-priori validation of the scaling form.
[Analytical solution section] Derivation of long-time D: the two-term formula for the diffusion coefficient assumes that the variance and persistence exponents of N(t) are fixed by the underlying particle dynamics and can be substituted without back-action. The paper does not quantify the regime of validity of this approximation (e.g., via a perturbative estimate of the neglected correlation terms), which is load-bearing for the central claim that the fluctuation term can dominate and produce anomalous scaling.

minor comments (2)

The manuscript would benefit from an explicit table or appendix listing the measured values of β, κ, d, and the resulting δ and α, together with their uncertainties, to make the numerical pathway from simulation inputs to predicted exponent transparent.
Notation for the two contributions to D (conventional vs. fluctuation-driven) could be introduced earlier and used consistently when discussing the crossover condition for anomalous scaling.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate the revisions we will make to clarify the model's assumptions and strengthen the validation.

read point-by-point responses

Referee: [Abstract / model setup] Abstract and model derivation: the analytical expression for the fluctuation-driven term ∝ N0^{-δ} (with δ=2-2/d-β+κ) is obtained only by treating N(t) as an independent exogenous Gaussian process whose statistics are inserted directly into the center-of-mass equation. The manuscript provides no demonstration that cross-correlations between δN(t) and cluster velocity (arising from attachment-induced displacements or motion-dependent attachment rates) remain negligible once the equations are fully coupled; any such correlations would renormalize the velocity autocorrelation integral and change the predicted scaling.

Authors: We agree that the closed-form analytical expression for the fluctuation-driven term is derived under the approximation that N(t) is an exogenous Gaussian process whose statistics are inserted into the center-of-mass equation. The manuscript begins from the single-particle Langevin equations to obtain the coupled SDEs, with attachment-induced displacements included in the center-of-mass dynamics. However, to close the equations analytically, the independence assumption is invoked. In the revised manuscript we will add an explicit check of the cross-correlation between δN(t) and cluster velocity, computed directly from the ABP simulations, and demonstrate that it remains small in the regime where the fluctuation term dominates. This discussion will be inserted in the model-setup section. revision: yes
Referee: [Abstract] Abstract and validation paragraph: the claim that the model 'predicts' α=0.63±0.06 'in good quantitative agreement' with ABP simulations is circular because the exponents β and κ that enter δ are extracted from the identical simulations. This reduces the numerical match to a post-hoc consistency test rather than an a-priori validation of the scaling form.

Authors: The referee correctly notes that β and κ are measured from the same simulations used for validation, so the numerical agreement constitutes a consistency check rather than an independent prediction. The primary contribution of the work is the derivation of the scaling relation δ = 2 − 2/d − β + κ from the coupled SDEs. We will revise the abstract and the validation paragraph to state that the model 'yields a scaling α = 0.63 ± 0.06 consistent with' the simulations, thereby accurately reflecting the nature of the comparison. revision: yes
Referee: [Analytical solution section] Derivation of long-time D: the two-term formula for the diffusion coefficient assumes that the variance and persistence exponents of N(t) are fixed by the underlying particle dynamics and can be substituted without back-action. The paper does not quantify the regime of validity of this approximation (e.g., via a perturbative estimate of the neglected correlation terms), which is load-bearing for the central claim that the fluctuation term can dominate and produce anomalous scaling.

Authors: We acknowledge that the two-term expression for D relies on fixed N(t) statistics without back-action and that the manuscript does not provide a quantitative bound on the neglected correlation terms. In the revised version we will add, in the analytical-solution section, a perturbative estimate of the leading-order cross terms together with an order-of-magnitude analysis showing under which conditions (in terms of N0 and activity parameters) these terms remain subdominant. This will delineate the regime in which the fluctuation-driven contribution can dominate and produce the reported anomalous scaling. revision: yes

Circularity Check

1 steps flagged

Numerical value of predicted scaling exponent α is computed directly from simulation-measured β and κ

specific steps

fitted input called prediction [Abstract]
"The model predicts D∼N^{-α} with α=0.63±0.06, in good quantitative agreement with large-scale ABP simulations."

α is obtained by substituting the simulation-extracted values of β and κ (which define the variance and persistence of the exogenous N(t) process) into δ=2-2/d-β+κ and then into the expression for D; the numerical agreement with ABP data is therefore fixed by those same measured inputs rather than emerging from the derivation alone.

full rationale

The analytical derivation of the two-term expression for D and the formula for δ=2-2/d-β+κ proceeds from the coupled SDEs under the stated Gaussian-process assumption for N(t) and is independent of the specific numerical values. However, the headline claim that the model 'predicts' α=0.63±0.06 reduces to inserting measured β and κ (obtained from the same ABP simulations) into that formula; the reported quantitative agreement is therefore a consistency check on the input exponents rather than an independent first-principles prediction. This matches the 'fitted input called prediction' pattern at the level of the final numerical result while leaving the core derivation non-circular.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on treating N(t) as an exogenous Gaussian process whose statistics are measured separately; this introduces two fitted exponents and assumes no back-coupling from position to attachment/detachment rates.

free parameters (2)

β
Exponent governing variance of N(t) ~ N0^β; measured from ABP simulations and inserted into δ.
κ
Exponent governing persistence time of N(t) ~ N0^κ; measured from ABP simulations and inserted into δ.

axioms (2)

domain assumption Single-particle dynamics obey standard overdamped Langevin equations with constant propulsion speed and rotational diffusion.
Invoked at the start of the derivation to obtain the center-of-mass equation.
ad hoc to paper N(t) is an independent Gaussian process with no feedback from the cluster's position or velocity.
Required to close the two-equation system and obtain the analytical solution.

pith-pipeline@v0.9.0 · 5560 in / 1662 out tokens · 48115 ms · 2026-05-08T16:36:15.023726+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

39 extracted references

[1]

Jan K. G. Dhont.An Introduction to Dynamics of Colloids, volume 2. Elsevier, 1996

1996
[2]

Oxford University Press, 1988

Masao Doi and Samuel Frederick Edwards.The Theory of Polymer Dynamics, volume 73. Oxford University Press, 1988

1988
[3]

Aug´ e, P

S. Aug´ e, P. O. Schmit, C. A. Crutchfield, et al. Nmr measure of translational diffusion and fractal dimension. application to molecular mass measurement.ACS Publications, 113:1914–1918, 2009

1914
[4]

Meakin and A

P. Meakin and A. T. Skjeltorp. Application of experimental and numerical models to the physics of multiparticle systems.Advances in Physics, 42(1):1–127, 1993

1993
[5]

Luczka, P

J. Luczka, P. H¨ anggi, and A. Gadomski. Diffusion of clusters with randomly growing masses. 51(6):5762–5769, 1995

1995
[6]

Rajesh, D

R. Rajesh, D. Das, B. Chakraborty, and M. Barma. Aggregate formation in a system of coagulating and fragmenting particles with mass-dependent diffusion rates.Phys. Rev. E, 66:056104, 2002

2002
[7]

Hidalgo-Soria and E

M. Hidalgo-Soria and E. Barkai. Hitchhiker model for laplace diffusion processes.Phys. Rev. E, 102:012109, 2020

2020
[8]

A. V. Chechkin, F. Seno, R. Metzler, and I. M. Sokolov. Brownian yet non-gaussian diffusion: From superstatistics to subordination of diffusing diffusivities.Phys. Rev. X, 7(2):021002, 2017. 8

2017
[9]

Baldovin, E

F. Baldovin, E. Orlandini, and F. Seno. Polymerization induces non-gaussian diffusion.Frontiers in Physics, 2019

2019
[10]

M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. Aditi Simha. Hydrodynamics of soft active matter.Rev. Mod. Phys., 85:1143, 2013

2013
[11]

Bechinger, R

C. Bechinger, R. Di Leonardo, H. L¨ owen, C. Reichhardt, G. Volpe, and G. Volpe. Active particles in complex and crowded environments.Rev. Mod. Phys., 88:045006, 2016

2016
[12]

Elgeti, R

J. Elgeti, R. G. Winkler, and G. Gompper. Physics of microswimmers—single particle motion and collective behavior: A review.Rep. Prog. Phys., 78:056601, 2015

2015
[13]

R. G. Winkler, J. Elgeti, and G. Gompper. Active polymers — emergent conformational and dynamical properties: A brief review.Journal of the Physical Society of Japan, 86:101014, 2017

2017
[14]

Bianco, E

V. Bianco, E. Locatelli, and P. Malgaretti. Globally and locally active polymers: shape, structure, and dynamics.Physical Review Letters, 121:217802, 2018

2018
[15]

C. P. Beatrici, R. M. C. de Almeida, and L. G. Brunnet. Mean-cluster approach indicates cell sorting time scales are determined by collective dynamics.Phys. Rev. E, 95(3):032402, 2017

2017
[16]

Anomalous mass dependency in hydra endoderm cell cluster diffusion

Aline L¨ utz, Carine Beatrici, L´ ıvio Amaral, and Leonardo Brunnet. Anomalous mass dependency in hydra endoderm cell cluster diffusion. 2024

2024
[17]

Deblais, A

A. Deblais, A. C. Maggs, D. Bonn, and S. Woutersen. Phase separation by entanglement of active polymerlike worms.Phys. Rev. Lett., 124(20):208006, 2020

2020
[18]

Tailleur and M

J. Tailleur and M. E. Cates. Statistical mechanics of interacting run-and-tumble bacteria.Phys. Rev. Lett., 100(21):218103, 2008

2008
[19]

Redner, M

G. Redner, M. Hagan, and A. Baskaran. Structure and dynamics of a phase-separating active colloidal fluid.Phys. Rev. Lett., 110:055701, 2013

2013
[20]

Microscopic theory for the phase separation of self-propelled repulsive disks.Europhys

Julian Bialk´ e, Hartmut L¨ owen, and Thomas Speck. Microscopic theory for the phase separation of self-propelled repulsive disks.Europhys. Lett., 103(3):30008, 2013

2013
[21]

Cates and J

M. Cates and J. Tailleur. Motility-induced phase separation.Annu. Rev. Condens. Matter Phys., 6:219, 2015

2015
[22]

Fily and M

Y. Fily and M. C. Marchetti. A thermal phase separation of self-propelled particles with no alignment. Phys. Rev. Lett., 108:235702, 2012

2012
[23]

Gonnella, D

G. Gonnella, D. Marenduzzo, A. Suma, and A. Tiribocchi. Motility-induced phase separation and coarsening in active matter.C. R. Phys., 16:316, 2015

2015
[24]

Cugliandolo, P

L. Cugliandolo, P. Digregorio, G. Gonnella, and A. Suma. Phase coexistence in two-dimensional passive and active dumbbell systems.Phys. Rev. Lett., 119:268002, 2017

2017
[25]

Cugliandolo, Giuseppe Gonnella, and Ignacio Pagonabarraga

Pasquale Digregorio, Demian Levis, Antonio Suma, Leticia F. Cugliandolo, Giuseppe Gonnella, and Ignacio Pagonabarraga. Full phase diagram of active brownian disks: From melting to motility-induced phase separation.Phys. Rev. Lett., 121(9):098003, 2018

2018
[26]

Digregorio, D

P. Digregorio, D. Levis, A. Suma, L. F. Cugliandolo, G. Gonnella, and I. Pagonabarraga. 2d melting and motility induced phase separation in active brownian hard disks and dumbbells.Journal of Physics: Conference Series, 1163(1):012073, 2019

2019
[27]

Caporusso, Pasquale Digregorio, Demian Levis, Leticia F

Claudio B. Caporusso, Pasquale Digregorio, Demian Levis, Leticia F. Cugliandolo, and Giuseppe Gonnella. Motility-induced microphase and macrophase separation in a two-dimensional active brownian particle system.Phys. Rev. Lett., 125:178004, Oct 2020

2020
[28]

C. B. Caporusso, L. F. Cugliandolo, P. Digregorio, G. Gonnella, D. Levis, and A. Suma. Dynamics of motility-induced clusters: Coarsening beyond ostwald ripening.Phys. Rev. Lett., 131(6):068201, 2023

2023
[29]

Ginot, I

F. Ginot, I. Theurkauff, F. Detcheverry, C. Ybert, and C. Cottin-Bizonne. Aggregation-fragmentation and individual dynamics of active clusters.Nature Communications, 9(1), 2018

2018
[30]

Self-propelled particle in an external potential: Existence of an effective temperature

Grzegorz Szamel. Self-propelled particle in an external potential: Existence of an effective temperature. Physical Review E, 90(1):012111, 2014. 9

2014
[31]

Active particles under confinement and effective force generation among surfaces.Soft matter, 14(44):9044–9054, 2018

Lorenzo Caprini and Umberto Marini Bettolo Marconi. Active particles under confinement and effective force generation among surfaces.Soft matter, 14(44):9044–9054, 2018

2018
[32]

G. H. P. Nguyen, R. Wittmann, and H. L¨ owen. Active ornstein-uhlenbeck model for self-propelled particles with inertia.Journal of Physics: Condensed Matter, 34(3):035101, 2021

2021
[33]

Martin, J

D. Martin, J. O’Byrne, M. E. Cates, et al. Statistical mechanics of active ornstein-uhlenbeck particles. Physical Review E, 103(3):032607, 2021

2021
[34]

Emerging mesoscale flows and chaotic advection in dense active matter.Physical Review Letters, 132(21):218301, 2024

Yann-Edwin Keta, Juliane U Klamser, Robert L Jack, and Ludovic Berthier. Emerging mesoscale flows and chaotic advection in dense active matter.Physical Review Letters, 132(21):218301, 2024

2024
[35]

Active ornstein-uhlenbeck particles.Physical Review E, 100(2):022601, 2019

Luis L Bonilla. Active ornstein-uhlenbeck particles.Physical Review E, 100(2):022601, 2019

2019
[36]

The parental active model: A unifying stochastic description of self-propulsion.The Journal of Chemical Physics, 156(7), 2022

Lorenzo Caprini, Alexander R Sprenger, Hartmut L¨ owen, and Ren´ e Wittmann. The parental active model: A unifying stochastic description of self-propulsion.The Journal of Chemical Physics, 156(7), 2022

2022
[37]

Dynamics of active particles with translational and rotational inertia.Journal of Physics: Condensed Matter, 35(30):305101, 2023

Alexander R Sprenger, Lorenzo Caprini, Hartmut L¨ owen, and Ren´ e Wittmann. Dynamics of active particles with translational and rotational inertia.Journal of Physics: Condensed Matter, 35(30):305101, 2023

2023
[38]

Computationally efficient recursions for top-order invariant polynomials with applications.Cambridge University Press, 1163(1):211–242, 2009

Grant Hillier, Raymond Kan, and Xiaolu Wang. Computationally efficient recursions for top-order invariant polynomials with applications.Cambridge University Press, 1163(1):211–242, 2009

2009
[39]

Vershynin.High-Dimensional Probability: An Introduction with Applications in Data Science

R. Vershynin.High-Dimensional Probability: An Introduction with Applications in Data Science. Cambridge University Press, 2018. 10

2018