Pith. sign in

REVIEW 4 major objections 5 minor 29 references

Wave front propagation in the Active Coagulation Model

T0 review · 4 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read In an active coagulation model, density perturbations switch from damped traveling waves to diffusive relaxation at a critical wavenumber $k_c=|\mu-\alpha|/(2v)$, and the diffusive limit is well defined only for small reaction rates…

desk verdict A clean linear-stability crossover in a new active coagulation model, with a KPZ claim that only works in the parameter regime where the main effect is absent. read the letter →

arxiv 2502.03372 v3 pith:STRCA3GS submitted 2025-02-05 cond-mat.stat-mech cond-mat.softphysics.bio-ph

Matteo Paoluzzi This is my paper · ORCID
classification cond-mat.stat-mechcond-mat.softphysics.bio-ph MSC 82C2282C3182C27
keywords run-and-tumbleparticlescoagulationmodelabsorbingstatephasetransitionwavefrontpropagationlong-livedcurrentsdiffusivelimitKPZuniversalitylinearstabilityanalysis
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper studies a minimal population model: self-propelled run-and-tumble particles that coagulate (annihilate) at rate $\beta$ and are born at rate $\mu$. It argues that active motion does not move the absorbing-state transition point, but it changes how the system relaxes: density and current stay coupled through the coagulation rate, long-lived currents form, and perturbations of wavenumber $k$ relax as damped traveling waves when $k>|\mu-\alpha|/(2v)$ and as purely diffusive modes when $k$ is smaller. The paper further claims that in the diffusive limit the effective diffusion coefficient is $D[\rho]=v^2/(\alpha+\mu-\beta\rho)$, which can diverge or become negative, so the diffusive description is only safe at small reaction rates; in that regime a noisy version of the dynamics maps to the KPZ equation. If correct, the result identifies the tumbling rate relative to the birth rate as the control parameter that decides whether a spreading population invades as a wave or as a diffusive blur.

What carries the argument

The central object is the coupled density-current system $\dot{\rho}=-J'-\mu\rho+\beta\rho^2$ and $\dot{J}=-v^2\rho'-J(\alpha+\mu-\beta\rho)$, where $\rho=R+L$ is the particle density and $J=v(R-L)$ is the current of right- and left-moving run-and-tumble particles. The mechanism that carries the argument is the coupling term $-J(\alpha+\mu-\beta\rho)$ in the current equation, which makes the current relaxation rate depend on the local density and thereby generates long-lived currents and the wavenumber-dependent traveling-to-diffusive crossover. Linear stability analysis of this system supplies the critical wavenumber $k_c$, while the adiabatic closure $\dot{J}=0$ supplies the effective diffusion coefficient whose divergence marks the breakdown of the diffusive limit.

What would settle it

Track the relaxation of a cosinusoidal density perturbation with wavenumber $k$ in a particle simulation or exact numerical solution of the full density-current equations; the claim predicts damped oscillations with a nonzero imaginary part in the dispersion relation whenever $k>|\mu-\alpha|/(2v)$. Observing purely monotonic decay at all wavenumbers for $\alpha<\mu$ would falsify the traveling-wave regime, and observing no sign change in $D[\rho]$ at $\beta\rho=\alpha+\mu$ would falsify the predicted breakdown of the diffusive limit.

Watch

Extended reading notes

Core claim

The central discovery is that for the active coagulation model the stationary mixed state $\rho_\infty=\mu/\beta$ is unchanged by activity, but the relaxation toward it is controlled by the ratio of the tumbling rate $\alpha$ to the reaction scales. Linearizing the coupled density-current equations around the mixed state gives a critical wavenumber $k_c=|\mu-\alpha|/(2v)$: perturbations with longer wavelengths ($k<k_c$) relax diffusively, while those with shorter wavelengths ($k>k_c$) first propagate as damped waves. The same density-current coupling produces a peak in the mean-field current $J(t)$, a signature of long-lived currents, and makes the diffusive closure $J=-D[\rho]\rho'$ pathological whenever $\beta\rho$ approaches $\alpha+\mu$, because $D[\rho]=v^2/(\alpha+\mu-\beta\rho)$ diverges and changes sign. In the small-$\beta$ regime the theory is well defined, and the fluctuating equation for the perturbation reduces to a KPZ-like equation; the paper therefore concludes that pattern formation at criticality should belong to the KPZ universality class.

Load-bearing premise

The conclusions about the diffusive limit and the KPZ mapping rest on the assumption that the reaction rates $\beta$ and $\mu$ are small enough that the effective diffusion coefficient $D[\rho]$ stays positive and finite; the paper itself flags that outside this regime $D[\rho]$ diverges and then changes sign.

Editorial extensions

If this is right

  • The absorbing-state transition point and the stationary mixed-state density remain at their mean-field values even when activity is strong, so directed-percolation universality is preserved despite the changed dynamics.
  • The tumbling rate $\alpha$ tunes between traveling-wave relaxation (small $\alpha$) and purely diffusive relaxation (large $\alpha$), meaning the relaxation of density perturbations directly exposes $k_c$ and hence the active speed $v$.
  • In the large-$\alpha$ Brownian limit the model reduces to the Fisher-Kolmogorov equation, so standard reaction-diffusion wave-front results apply in that regime.
  • The diffusion coefficient becoming negative or divergent signals that the current cannot always be eliminated; a well-defined diffusive description exists only for small reaction rates.
  • Near the absorbing-state critical point ($\mu\to 0$), the noisy dynamics is formally equivalent to KPZ, predicting KPZ scaling for large-scale pattern formation.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same density-current mechanism should appear whenever active motion is coupled to any reaction with a density-dependent rate; replacing coagulation by SIS or SIR reactions and checking for the same $k_c$ crossover would be a direct test of the generality of this picture.
  • The divergence of $D[\rho]$ at $\rho_c=(\alpha+\mu)/\beta$ suggests that the diffusive limit is not merely a technical approximation but marks a regime where the current must be kept as a slow field, potentially leading to a two-field KPZ-like description instead of a single scalar equation.
  • The KPZ mapping at $\mu=0$ implies that active noise could alter the roughening exponents of an expanding front, giving a measurable signature of activity even though the phase boundary is unchanged.
  • Because the model is one-dimensional, a natural extension is to two dimensions, where the current is a vector and the traveling-to-diffusive crossover may interact with vorticity or banding instabilities; this is beyond what the paper addresses.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript presents a one-dimensional continuum model of run-and-tumble particles with coagulation and decoagulation reactions, summarized in Eqs. (4)-(5). It solves the mean-field logistic density and the current peak time, performs a linear stability analysis of the mixed stationary state, derives a diffusive closure for the current, and uses that closure to discuss Fisher-Kolmogorov front propagation and a formal mapping to the Kardar-Parisi-Zhang equation. Numerical solutions of the PDEs are used to illustrate the qualitative regimes.

Significance. The paper is clearly written and has the virtue of exact mean-field and linear-stability results with no parameter fitting; in particular, the predicted crossover wavenumber kc=|mu-alpha|/(2v) and the mean-field current peak time are explicit and falsifiable. The identification of the singular diffusive closure is honest and important. However, several algebraic errors in the linear-stability and diffusive-limit derivations, and the questionable KPZ mapping at the absorbing-state critical point, currently undermine the pattern-formation claims and require correction before the results can be fully endorsed.

major comments (4)
  1. [Linear Stability Analysis, Eqs. (24), (28), (30)] Eq. (24) is inconsistent with Eq. (28). Linearizing Eq. (5) around rho_infinity = mu/beta gives delta J_dot = -v^2 delta rho' - alpha delta J, because beta rho_infinity = mu; the printed coefficient (alpha+mu) would lead to different eigenvalues. Eq. (28), which is the correct result for the physical linearization, should be derived from the corrected equation. In addition, the text states that stability requires Re(sigma) > 0, whereas the eigenvalues in Eq. (28) have Re(sigma) < 0 for all k; the condition is Re(sigma) < 0. Finally, Eq. (30) should read kc = |mu-alpha|/(2v), not (mu-alpha)/(2v), since the threshold follows from (mu-alpha)^2 - (2vk)^2 < 0.
  2. [Diffusive limit and pattern formation, Eqs. (35)-(42)] The small-beta expansion of D[rho] is algebraically incorrect. For rho = rho_infinity + delta rho with rho_infinity = mu/beta, the denominator in Eq. (36) is alpha+mu-beta rho = alpha - beta delta rho. Hence D(rho_infinity) = v^2/alpha and the Taylor coefficient is v^2 beta/alpha^2. Equations (39) and (42) instead use tilde D = v^2/(alpha+mu) and tilde beta = beta/(beta+alpha). The resulting diffusion constant and the coefficient of (delta rho')^2 in Eq. (42) are therefore not the small-beta limit of Eq. (36). In particular, the claimed dependence of the KPZ amplitude on mu/alpha is an artifact of this expansion error; the corrected coefficient is v^2 beta/alpha^2.
  3. [Diffusive limit and pattern formation, Eqs. (18), (35)-(36)] The long-lived current effect advertised in the abstract is the mean-field peak in J(t), which by Eq. (18) exists only for beta > alpha+mu. In this regime the density rho_c = (alpha+mu)/beta lies only alpha/beta above the stationary density mu/beta, so the denominator in D[rho] is close to zero and the adiabatic closure J = -D[rho] rho' is uncontrolled. The diffusive and KPZ results are derived in the small-beta regime, where D is regular but the mean-field current decays monotonically. The manuscript acknowledges the singularity but does not supply an argument connecting the two regimes; as written, the traveling-wave/diffusive-wave crossover does not require the closure, while the FK/KPZ pattern-formation claims are confined to the opposite parameter regime. This limitation should be stated explicitly, and the phrase ranging from Fisher-Kolmogorov to Kardar-Parisi-Zhang in the abstract should be qualified.
  4. [Diffusive limit and pattern formation, Eq. (43) and Conclusions] The step from Eq. (42) to the KPZ statement is not justified. Equation (42) is an expansion around the mixed stationary state with mu > 0; its linear term is -mu delta rho. Setting mu = 0 does not describe the absorbing-state critical point in the same expansion, because at mu = 0 the stationary state is rho = 0 and the reaction term is -beta rho^2, which is not negligible at criticality. The Langevin equation normally associated with this absorbing-state transition belongs to the directed-percolation universality class, as the final paragraph itself states. The relation between the claimed KPZ universality for pattern formation at criticality and the stated DP universality class of the transition must be clarified; as it stands, the formal equivalence in Eq. (43) is not established.
minor comments (5)
  1. [Large-alpha limit, Eq. (32)] Eq. (32) has J[rho] = DA rho'; the correct closure is J[rho] = -DA rho', and the Fisher-Kolmogorov equation (33) uses the latter sign. Please fix the sign in Eq. (32).
  2. [Mean-field approximation, Eq. (18)] The statement 'with t = 0 if beta < alpha+mu' should be 'if beta <= alpha+mu', since at beta = alpha+mu the argument of the logarithm is unity.
  3. [Fig. 6 and surrounding text] The text before Fig. 6 says the two tumbling rates are alpha = 0.5 and alpha = 2, while the caption says alpha = 0.1 and alpha = 2; please reconcile the values.
  4. [Appendix A] The explicit Euler scheme with central differences is dispersive for the hyperbolic terms, and no convergence test or resolution study is reported; a brief convergence statement would strengthen the numerical evidence.
  5. [Introduction] The phrase 'birth/date processes' appears to be a typo for 'birth/death processes'.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the paper's central results are derived directly from the stated equations with free parameters, no fitted inputs, and no load-bearing self-citations.

full rationale

I walked the derivation chain from the stated run-and-tumble plus reaction equations (2)-(5). All parameters (α, β, μ, v) are free inputs; none is fitted to data or extracted from the phenomena being predicted. The stationary density ρ∞ = μ/β follows immediately from setting time and space derivatives to zero, so the activity independence of the transition point is a direct algebraic consequence of the model definition rather than a circularly assumed conclusion. The long-lived-current peak and its time scale (Eq. 18) are analytic solutions of the mean-field equations (12)-(16). The traveling-to-diffusive crossover (Eqs. 28-30) is derived from the dispersion relation of the linearized system (23)-(26). The diffusive closure (Eqs. 35-36) is obtained by the explicit adiabatic approximation d_t J = 0, and the paper openly flags the divergence and sign change of D[ρ]; this is a limitation of the approximation, not circularity. The KPZ discussion (Eqs. 40-43) adds noise explicitly and identifies the resulting equation as formally KPZ-like; it does not re-import the target universality claim as an input. Self-citations [15,17,18] appear only in contextual references and in a side remark about ∇^4 regularization; they are not load-bearing for the paper's conclusions. The acknowledged small-β restriction on the diffusive limit is a physical-consistency caveat, not a circular reduction. Therefore the paper is self-contained in its derivation chain, and no circular step is present.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claims rest on the mean-field and diffusive closures listed above. The model parameters alpha, beta, mu, and v are physical inputs scanned numerically, not fitted to data, so no free parameters are counted. No new entities are introduced.

assumptions (4)
  • domain assumption The coarse-grained PDEs (4)-(5) faithfully represent the microscopic run-and-tumble plus coagulation/decoagulation lattice model.
    The reaction rates are inserted into the hydrodynamic equations for right- and left-movers without deriving them from a lattice master equation; this assumes local density-dependent reactions and a well-mixed subgrid.
  • domain assumption The well-mixed closure drops all spatial derivatives so that rho evolves by the logistic equation independently of J.
    Used in the Well-Mixed Approximation section to derive Eqs (12)-(18); ignores spatial correlations and fluctuations.
  • domain assumption The diffusive limit is obtained by setting d_t J = 0, giving J = -D[rho] rho'.
    Used around Eqs (35)-(37); the paper itself notes this closure is pathological when D diverges or becomes negative.
  • ad hoc to paper A white-noise term with amplitude T is added to Eq (42) to represent fast degrees of freedom.
    The noise is not derived from the microscopic dynamics; it is the basis for the claimed mapping to the KPZ universality class.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Wave front propagation in the Active Coagulation Model." pith.science (2026). https://pith.science/paper/STRCA3GS

@misc{pith2026250203372,
  author       = {Pith},
  title        = {Pith review of: Wave front propagation in the Active Coagulation Model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/STRCA3GS}},
  note         = {Machine review of arXiv:2502.03372}
}
read the original abstract

Spreading processes on top of active dynamics provide a novel theoretical framework for capturing emerging collective behavior in living systems. I consider run-and-tumble dynamics coupled with coagulation/decoagulation reactions that lead to an absorbing state phase transition. While the active dynamics does not change the location of the transition point, the relaxation toward the stationary state depends on motility parameters. Because of the competition between spreading dynamics and active motion, the system can support long-living currents whose typical time scale is a nontrivial function of motility and reaction rates. Because of this interplay between time-scales, the wave front propagation qualitatively changes from traveling to diffusive waves. Moving beyond the mean-field regime, instability at finite length scales regulates a crossover from periodic to diffusive modes. Finally, it is possible to individuate different mechanisms of pattern formation on a large time scale, ranging from the Fisher-Kolmogorov to the Kardar-Parisi-Zhang equation.

Figures

Figures reproduced from arXiv: 2502.03372 by the authors.

Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Mean-Field approximation. (a) Time evolution of density [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. FIG. 4. Phase Diagram. (a) Stationary density [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figures from the paper (3 more)
Figure 5
Figure 5. Figure 5: FIG. 5. Linear stability. (a) Contour plot of ∆( [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Dynamics. The color map indicates the time evolu [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Diffusive Limit. (a) Relaxation of a Gaussian density [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

29 extracted references · 28 canonical work pages

  1. [1]

    Fruchart, R

    M. Fruchart, R. Hanai, P. B. Littlewood and V. Vitelli, Nature, 2021, 592, 363–369

    work page 2021

  2. [2]

    Sompolinsky, A

    H. Sompolinsky, A. Crisanti and H. J. Sommers, Phys. Rev. Lett., 1988, 61, 259–262

    work page 1988

  3. [3]

    Hinrichsen, Advances in physics , 2000, 49, 815–958

    H. Hinrichsen, Advances in physics , 2000, 49, 815–958

    work page 2000

  4. [4]

    Crisanti, A

    A. Crisanti, A. Puglisi and D. Villamaina, Phys. Rev. E , 2012, 85, 061127

    work page 2012

  5. [5]

    S. Saha, J. Agudo-Canalejo and R. Golestanian, Phys. Rev. X, 2020, 10, 041009

    2020

  6. [6]

    Pisegna, S

    G. Pisegna, S. Saha and R. Golestanian, Proceed- ings of the National Academy of Sciences , 2024, 121, e2407705121

    work page 2024

  7. [7]

    Dinelli, J

    A. Dinelli, J. O’Byrne, A. Curatolo, Y. Zhao, P. Sollich and J. Tailleur, Nature Communications, 2023, 14, 7035

    work page 2023

  8. [8]

    Y. Avni, M. Fruchart, D. Martin, D. Seara and V. Vitelli, arXiv preprint arXiv:2409.07481 , 2024

    work page Pith review arXiv 2024

Show all 29 references
  1. [9]

    J. D. Murray, Mathematical biology: I. An introduction , Springer Science & Business Media, 2007, vol. 17

    2007

  2. [10]

    Hallatschek, P

    O. Hallatschek, P. Hersen, S. Ramanathan and D. R. Nelson, Proceedings of the National Academy of Sciences, 2007, 104, 19926–19930

    2007

  3. [11]

    V. E. Deneke and S. Di Talia, Journal of Cell Biology , 2018, 217, 1193–1204

    2018

  4. [12]

    Peruani and G

    F. Peruani and G. J. Sibona, Phys. Rev. Lett., 2008, 100, 168103

    2008

  5. [13]

    Peruani and G

    F. Peruani and G. J. Sibona, Soft matter, 2019, 15, 497– 503

    2019

  6. [14]

    Gascuel, P

    H.-M. Gascuel, P. Rahmani, R. Bon and F. Peruani, Phys. Rev. Lett. , 2024, 133, 058301

    2024

  7. [15]

    J. P. Rodr ´ ıguez, M. Paoluzzi, D. Levis and M. Starnini, Phys. Rev. Res. , 2022, 4, 043160

    2022

  8. [16]

    Levis, A

    D. Levis, A. Diaz-Guilera, I. Pagonabarraga and M. Starnini, Phys. Rev. Res. , 2020, 2, 032056

    2020

  9. [17]

    Paoluzzi, M

    M. Paoluzzi, M. Leoni and M. C. Marchetti, Soft Matter, 2020, 16, 6317–6327

    2020

  10. [18]

    Paoluzzi, M

    M. Paoluzzi, M. Leoni and M. C. Marchetti, Phys. Rev. E, 2018, 98, 052603

    2018

  11. [19]

    Marcolongo, F

    B. Marcolongo, F. Peruani and G. J. Sibona, Physica A: Statistical Mechanics and its Applications , 2024, 648, 129916

    2024

  12. [20]

    Peruani and C

    F. Peruani and C. F. Lee, Europhysics Letters, 2013, 102, 58001

    2013

  13. [21]

    Norambuena, F

    A. Norambuena, F. J. Valencia and F. Guzm´ an-Lastra, Scientific Reports, 2020, 10, 20845

    2020

  14. [22]

    de Castro, F

    P. de Castro, F. Urbina, A. Norambuena and F. Guzm´ an- Lastra, Phys. Rev. E , 2023, 108, 044104

    2023

  15. [23]

    Lib´ al, P

    A. Lib´ al, P. Forg´ acs, A. N´ eda, C. Reichhardt, N. Hen- gartner and C. J. O. Reichhardt, Phys. Rev. E , 2023, 107, 024604. 9

    2023

  16. [24]

    Forg´ acs, A

    P. Forg´ acs, A. Lib´ al, C. Reichhardt, N. Hengartner and C. J. Reichhardt, Communications Physics, 2023, 6, 294

    2023

  17. [25]

    M. E. Cates, D. Marenduzzo, I. Pagonabarraga and J. Tailleur, Proceedings of the National Academy of Sci- ences, 2010, 107, 11715–11720

    2010

  18. [26]

    Curatolo, N

    A. Curatolo, N. Zhou, Y. Zhao, C. Liu, A. Daerr, J. Tailleur and J. Huang,Nature Physics, 2020, 16, 1152– 1157

    2020

  19. [27]

    M. J. Schnitzer, Phys. Rev. E , 1993, 48, 2553–2568

    1993

  20. [28]

    Kac, The Rocky Mountain Journal of Mathematics , 1974, 4, 497–509

    M. Kac, The Rocky Mountain Journal of Mathematics , 1974, 4, 497–509

    1974

  21. [29]

    Kardar, G

    M. Kardar, G. Parisi and Y.-C. Zhang, Phys. Rev. Lett., 1986, 56, 889–892

    1986

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.