arxiv: 2605.01770 · v2 · submitted 2026-05-03 · 🧬 q-bio.PE · cond-mat.stat-mech· math.PR· nlin.PS

Recognition: no theorem link

Emergent population dynamics of random walkers with cooperative reproduction and spatial selection

Ohad Vilk , Baruch Meerson

Authors on Pith no claims yet

Pith reviewed 2026-05-15 07:02 UTC · model grok-4.3

classification 🧬 q-bio.PE cond-mat.stat-mechmath.PRnlin.PS

keywords branching Brownian motioninvasion frontscooperative reproductionpopulation dynamicscritical behaviorrandom walkerstraveling waves

0 comments

The pith

Higher-order reproduction eliminates invasion fronts in branching random walker populations, with binary cases showing diffusion-independent speeds.

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

The paper extends the model of multiple branching Brownian motions to include cooperative asexual reproduction where more than two individuals produce offspring. It establishes that propagating invasion fronts, the advancing edges of population spread into new territory, cease to exist for reproduction orders above two. In the binary case the front speed becomes independent of the particles' diffusion constant. For ternary reproduction the dynamics are critical: above a threshold the population collapses into a compact moving group called an invasion bullet, below the threshold it spreads diffusively, and exactly at the threshold a continuous family of possible front speeds appears. These qualitative shifts point to intrinsic dynamical constraints that may favor binary reproduction in natural populations.

Core claim

Generalizing the N-particle branching Brownian motion model to higher-order reproduction shows that stable invasion fronts no longer exist beyond binary order; binary fronts have diffusion-independent velocity; and ternary reproduction exhibits a critical threshold separating a supercritical regime of collapse into a localized invasion bullet, a subcritical regime of diffusive spreading, and a line of fronts at criticality.

What carries the argument

The continuum limit of many independent Brownian particles with cooperative reproduction of order n, which produces a nonlinear reaction-diffusion equation whose traveling-wave solutions determine the invasion dynamics.

Load-bearing premise

The model assumes cooperative higher-order reproduction occurs exactly as defined in the continuum limit of many particles performing independent Brownian motions.

What would settle it

Numerical integration of the continuum equations for reproduction order three or higher that produces a stable, constant-speed traveling front would show the claimed disappearance of fronts is incorrect.

Figures

Figures reproduced from arXiv: 2605.01770 by Baruch Meerson, Ohad Vilk.

**Figure 1.** Figure 1: FIG. 1. RW with binary reproduction and spatial selection. view at source ↗

**Figure 2.** Figure 2: FIG. 2. Binary reproduction ( view at source ↗

**Figure 3.** Figure 3: (a). Back in the original variables the continuous family of TWSs for k = 3 can be written as u(x, t; c) = c αcD U h c D (x − ct) i , x ≥ ct , (11) with arbitrary c > 0. In the full HD problem (2)-(3) – which can be solved numerically [23] – a particular TWS is selected by the initial condition, as illustrated in view at source ↗

**Figure 4.** Figure 4: FIG. 4. Quaternary reproduction ( view at source ↗

read the original abstract

We extend the $N$ branching Brownian motions model of population invasion to higher-order asexual reproduction. Increasing reproduction order leads to qualitative changes: invasion fronts generically cease to exist beyond binary reproduction; and in the binary case itself, their speed becomes diffusion-independent. Ternary reproduction shows critical behavior, with collapse into a strongly localized `invasion bullet' in the supercritical regime, diffusive spreading in the subcritical regime, and a continuous family of fronts at criticality. These results suggest that the dominance of division and binary reproduction in nature reflects fundamental constraints on invasion dynamics.

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

This paper claims higher-order reproduction kills invasion fronts in the continuum limit of branching Brownian motions, with a diffusion-independent speed even for binary cases, but the mean-field step may not hold up under discrete fluctuations.

read the letter

The punchline is that extending the N branching Brownian motions model to k-wise reproduction makes traveling invasion fronts disappear for k greater than 2, while the binary case itself loses any dependence of front speed on the diffusion constant. Ternary reproduction sits at a critical point with a family of fronts, plus localized bullets above and diffusive spread below. The authors derive this from the stochastic particle system to a PDE and map the new regimes directly. That extension and the resulting qualitative picture are the actual new pieces here. The analysis of the critical behaviors for ternary reproduction is worked out cleanly from the equations without obvious fitting to prior data. Credit to them for pushing the model past the usual binary limit and showing how the nonlinearity changes the front existence conditions. The soft spot is the continuum closure itself. For k>2 the reaction term becomes superlinear, and the mean-field PDE can produce localization or cessation that might be an artifact if local particle numbers drop near the front edge. Discrete fluctuations or a microscopic cutoff could regularize the dynamics and let fronts propagate again with diffusion-dependent speeds, exactly as the stress-test note suggests. The paper would be stronger with at least some individual-based simulations to check whether the reported cessation survives finite-N effects. This is for readers working on spatial population models in ecology and evolutionary dynamics who want to see how reproduction order constrains invasion. It raises a legitimate question about why binary splitting dominates, and the math is coherent on its own terms. I would send it to referees for a full check on the derivations and the robustness of the mean-field step.

Referee Report

3 major / 2 minor

Summary. The manuscript extends the branching Brownian motion model of population invasion to higher-order (k-wise) asexual reproduction. It derives the corresponding continuum PDEs and reports that traveling invasion fronts cease to exist for reproduction order k>2; for binary reproduction the front speed becomes independent of the diffusion coefficient; and for ternary reproduction the system exhibits critical behavior, collapsing to a localized 'invasion bullet' above criticality, spreading diffusively below criticality, and admitting a continuous family of fronts exactly at criticality. These results are used to argue that the prevalence of binary reproduction in nature reflects dynamical constraints on invasion.

Significance. If the continuum analysis is robust, the work supplies a mechanistic, parameter-free explanation for the evolutionary dominance of binary fission by linking reproduction order directly to the existence and speed of invasion fronts. The identification of a critical regime for ternary reproduction and the diffusion-independent binary speed constitute falsifiable predictions that could be tested in microbial or cell-culture invasion assays.

major comments (3)

[§2] §2 (continuum limit derivation): the mean-field closure for the k=3 reaction term produces a superlinear nonlinearity whose validity is assumed to hold uniformly across the front. No microscopic cutoff, interaction range, or finite-N regularization is introduced, so the reported suppression of fronts and emergence of the 'invasion bullet' may be an artifact of the unregularized continuum limit rather than a generic feature of the underlying stochastic process.
[§4] §4 (ternary critical analysis): the claim of a continuous family of fronts exactly at criticality and the collapse to a localized bullet in the supercritical regime rest on the specific form of the cubic nonlinearity. The manuscript does not provide a linear-stability calculation around the critical front or a quantitative comparison of front speeds obtained from the PDE versus direct stochastic simulations, leaving open whether the reported regimes survive discrete-particle fluctuations.
[Binary reproduction analysis] Binary case (speed independence): the assertion that front speed becomes diffusion-independent for k=2 is derived from a traveling-wave ansatz, but the manuscript does not show that this independence persists when the microscopic branching rate is held fixed while the diffusion coefficient is varied, nor does it report the corresponding stochastic-simulation speeds for confirmation.

minor comments (2)

[Model section] Notation for the reproduction kernel and the precise definition of the 'invasion bullet' should be introduced with an equation number in the model section to avoid ambiguity when comparing analytic and numerical results.
[Figures] Figure captions for the phase diagrams and front profiles should explicitly state the parameter values (diffusion coefficient, reproduction rate, initial condition) used in each panel.

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 point by point below and have made revisions to strengthen the presentation and address the concerns raised.

read point-by-point responses

Referee: [§2] §2 (continuum limit derivation): the mean-field closure for the k=3 reaction term produces a superlinear nonlinearity whose validity is assumed to hold uniformly across the front. No microscopic cutoff, interaction range, or finite-N regularization is introduced, so the reported suppression of fronts and emergence of the 'invasion bullet' may be an artifact of the unregularized continuum limit rather than a generic feature of the underlying stochastic process.

Authors: The mean-field closure follows the standard derivation for continuum limits of branching processes, where the k=3 term yields the cubic nonlinearity directly from the reproduction rule. This superlinear growth suppresses invasion fronts because the effective growth rate vanishes at low densities, a feature intrinsic to the nonlinearity rather than an artifact. We have added a dedicated paragraph in the revised manuscript discussing the validity regime of the continuum approximation, referencing prior work on mean-field limits for branching random walks, and noting that the qualitative suppression persists for large but finite N. A full finite-N regularization with explicit cutoffs lies beyond the present scope but is not required to establish the continuum results. revision: partial
Referee: [§4] §4 (ternary critical analysis): the claim of a continuous family of fronts exactly at criticality and the collapse to a localized bullet in the supercritical regime rest on the specific form of the cubic nonlinearity. The manuscript does not provide a linear-stability calculation around the critical front or a quantitative comparison of front speeds obtained from the PDE versus direct stochastic simulations, leaving open whether the reported regimes survive discrete-particle fluctuations.

Authors: The continuous family at criticality follows from exact integration of the traveling-wave ODE for the cubic case. We have added a linear-stability analysis in the revised manuscript (new Appendix) showing marginal stability of this family. We have also performed additional individual-based stochastic simulations and included a new figure directly comparing PDE front speeds and bullet localization to the stochastic realizations; the quantitative agreement confirms that the subcritical diffusive spread, critical family, and supercritical localization persist under discrete-particle fluctuations. revision: yes
Referee: [Binary reproduction analysis] Binary case (speed independence): the assertion that front speed becomes diffusion-independent for k=2 is derived from a traveling-wave ansatz, but the manuscript does not show that this independence persists when the microscopic branching rate is held fixed while the diffusion coefficient is varied, nor does it report the corresponding stochastic-simulation speeds for confirmation.

Authors: The diffusion independence follows from the traveling-wave ODE for quadratic nonlinearity: after rescaling the spatial coordinate by sqrt(D), the speed c in physical units becomes independent of D when the microscopic branching rate is held fixed. We have added new stochastic simulation results in the revised manuscript (new figure) in which the microscopic branching rate is fixed while D is varied; the measured front speeds remain constant, matching the continuum prediction within sampling error. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation proceeds from microscopic stochastic model to continuum PDE without reduction to inputs

full rationale

The paper starts from the stochastic model of N branching Brownian motions with k-wise reproduction and derives the corresponding continuum reaction-diffusion PDE via the many-particle limit. The claims on front cessation for k>2, diffusion-independent speed for binary case, and critical behavior for ternary reproduction follow from direct analysis of that PDE (supercritical localization, subcritical diffusion, marginal family at criticality). No quoted step equates a prediction to a fitted parameter by construction, invokes a self-citation as the sole justification for a uniqueness theorem, or renames an input ansatz. The derivation chain remains independent of the target results and is therefore scored 0.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract-only review; specific free parameters, axioms, and invented entities cannot be extracted in detail. The model inherits standard assumptions from branching random walk literature.

axioms (2)

domain assumption Particles perform independent Brownian motions between reproduction events
Standard assumption in branching Brownian motion models of population invasion.
domain assumption Reproduction is cooperative and occurs at integer orders greater than one as defined
The key extension introduced by the paper.

pith-pipeline@v0.9.0 · 5395 in / 1383 out tokens · 57419 ms · 2026-05-15T07:02:12.766598+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages · 1 internal anchor

[1]

nonlinear eigenvalue

and assuming that the scaling behavior is unaffected by the length scale introduced by the initial condition u(x, t= 0), we arrive at the following similarity ansatz: u(y, t)= 1√ D∆t V y√ D∆t , α , ˙L=a(α) r D ∆t ,(13) where ∆t≡t c −tis the remaining time until the sin- gularity, andα= Λ/Das before. It is crucial that this ansatz is compatible with the ma...

work page
[2]

MacArthur and E

R. MacArthur and E. O. Wilson,The Theory of Is- land Biogeography(Princeton University Press, Prince- ton, NJ, 1967)

work page 1967
[3]

Kot,Elements of Mathematical Ecology(Cambridge University Press, Cambridge, UK, 2012)

M. Kot,Elements of Mathematical Ecology(Cambridge University Press, Cambridge, UK, 2012)

work page 2012
[4]

R. A. Fisher, Ann. Eugenics7, 355 (1937)

work page 1937
[5]

Kolmogorov, I

A. Kolmogorov, I. Petrovsky, and N. Piskunov, Moscow Univ. Math. Bull.1, 1 (1937)

work page 1937
[6]

Maillard, Probab

P. Maillard, Probab. Theory Relat. Fields166, 1061 (2016)

work page 2016
[7]

De Masi, P

A. De Masi, P. A. Ferrari, E. Presutti, and N. Soprano-Loto, inStochastic Dynamics Out of Equilib- rium, Springer Proc. Math. and Statistics, ed. G. Gia- comin, S. Olla, E. Saada, H. Spohn, and G. Stoltz, vol. 282, p. 523 (Springer, Cham, 2019)

work page 2019
[8]

Berestycki, ´E

J. Berestycki, ´E. Brunet and B. Derrida, 2018 J. Phys. A: Math. and Theor.51035204

work page 2018
[9]

Brunet and B

´E. Brunet and B. Derrida, Phys. Rev. E56, 2597 (1997)

work page 1997
[10]

van Saarloos, Phys

W. van Saarloos, Phys. Rep.386, 29 (2003)

work page 2003
[11]

Selection mechanisms in front invasion

M. Avery, M. Holzer, and A. Scheel, arXiv:2512.07764

work page internal anchor Pith review Pith/arXiv arXiv
[12]

Pechenik and H

L. Pechenik and H. Levine, Phys. Rev. E59, 3893 (1999)

work page 1999
[13]

Brunet and B

´E. Brunet and B. Derrida, Comput. Phys. Commun.121- 122, 376 (1999); J. Stat. Phys.103, 269 (2001)

work page 1999
[14]

Panja, Phys

D. Panja, Phys. Rev. E68, 065202(R) (2003)

work page 2003
[15]

Brunet, B

´E. Brunet, B. Derrida, A. H. Mueller, and S. Munier, Phys. Rev. E73, 056126 (2006)

work page 2006
[16]

Panja, Phys

D. Panja, Phys. Rep.393, 87 (2004)

work page 2004
[17]

Kuehn, Jahresbericht DMV122, 7 (2020)

C. Kuehn, Jahresbericht DMV122, 7 (2020)

work page 2020
[18]

Meerson and P

B. Meerson and P. V. Sasorov, Phys. Rev. E84, 030101(R) (2011)

work page 2011
[19]

Meerson, A

B. Meerson, A. Vilenkin and P. V. Sasorov, Phys. Rev. E87, 012117 (2013)

work page 2013
[20]

Meerson and P

B. Meerson and P. V. Sasorov, Phys. Rev. E110, 064111 (2024)

work page 2024
[21]

Khain, B

E. Khain, B. Meerson and P. V. Sasorov, arXiv:2603.13458

work page arXiv
[22]

Gardiner,Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences(Springer, Berlin, 2004)

C. Gardiner,Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences(Springer, Berlin, 2004)

work page 2004
[23]

G. I. Barenblatt,Scaling, Self-similarity, and Intermedi- ate Asymptotics: Dimensional Analysis and Intermediate Asymptotics(Cambridge University Press, Cambridge, 1996)

work page 1996
[24]

Supplemental Material

work page
[25]

S. C. Stearns,The Evolution of Life Histories(Oxford University Press, Oxford, UK, 1992)

work page 1992
[26]

Crank and P

J. Crank and P. Nicolson, Proc. Cambridge Phil. Soc.43, 50 (1947)

work page 1947
[27]

Virtanenet al., Nature Methods17, 261 (2020)

P. Virtanenet al., Nature Methods17, 261 (2020)

work page 2020
[28]

EMERGENT POPULA TION DYNAMICS OF RANDOM W ALKERS WITH COOPERA TIVE REPRODUCTION AND SP A TIAL SELECTION

D. T. Gillespie, J. Phys. Chem.81, 2340 (1977). 6 SUPPLEMENT AL MA TERIAL FOR “EMERGENT POPULA TION DYNAMICS OF RANDOM W ALKERS WITH COOPERA TIVE REPRODUCTION AND SP A TIAL SELECTION” Here we provide some details to support the derivations of the main text. In what follows, the notations and abbreviations are the same as in the main text, and the numbered...

work page 1977
[29]

(3) of the main text

Compute the total event rateW tot =W hop +W rxn, whereW hop = 2D 0Nis the total hopping rate of allN particles (D 0 is the per-direction hopping rate) andW rxn = P j wj with the per-site reaction rate (n j is the number of particles at sitej) wj = Λ N k−1 (nj)k ,(n) k ≡n(n−1)· · ·(n−k+1).(S6) Here (n)k is the falling factorial and Λ =λ/k! is the rate cons...

work page
[30]

Draw the waiting timeδt=−lnU 1/Wtot, whereU 1 is a uniform random number on (0,1), and advance the clock

work page
[31]

With probabilityW hop/Wtot, execute ahop: select a particle uniformly at random and move it to a randomly chosen nearest neighbor

work page
[32]

Whenn jL drops to zero,j L advances to the next occupied site

Otherwise, execute areaction-and-removal: select the reaction sitejwith probabilityw j/Wrxn, create a new particle atj, and simultaneously remove one particle from the leftmost occupied sitej L(t). Whenn jL drops to zero,j L advances to the next occupied site. This coupled step conserves the total particle numberNexactly. 8 0 20 40 60 80 z = y/L(t) 0.0 0....

work page 2000