pith. sign in

arxiv: 2605.29958 · v1 · pith:DZ3YKKWOnew · submitted 2026-05-28 · 🧬 q-bio.PE · cond-mat.stat-mech· math.PR· nlin.PS

Lattice Brownian bees with cooperative reproduction: steady states, collapse, and spreading

Ohad Vilk , Baruch Meerson This is my paper

Pith reviewed 2026-06-28 23:38 UTC · model grok-4.3

classification 🧬 q-bio.PE cond-mat.stat-mechmath.PRnlin.PS
keywords brownian beescooperative reproductionhydrodynamic limitfree-boundary problempopulation collapsediffusive spreadingrandom walkerslinear stability
0
0 comments X

The pith

Cooperative reproduction makes steady states of Brownian bees stable only for k≤2, with a critical family at k=3 separating collapse from spreading and instability for k≥4.

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

The paper extends the Brownian bees model of random walkers with removal of the farthest particle to include cooperative reproduction of order k. In the infinite-population limit it replaces the microscopic rules with a hydrodynamic free-boundary problem whose steady-state densities can be found for every k. Linear stability analysis shows these states are stable when k is at most 2 and unstable when k is at least 4. At the marginal value k=3 a continuous family of steady states exists only at one critical ratio of reproduction to diffusion rates; above that ratio the population collapses in finite time in a self-similar way, while below it the population spreads diffusively with reproduction still quantitatively important. For k greater than or equal to 4 the unstable steady state separates a collapse regime from a spreading regime, the collapse requiring matched asymptotics because of an inner scale separation.

Core claim

In the N to infinity limit the model is described by a hydrodynamic free-boundary problem whose steady-state solutions are linearly stable for k less than or equal to 2 and unstable for k greater than or equal to 4. For k equals 3 there exists a continuous family of steady states at one critical ratio of reproduction to diffusion rates, with asymptotically self-similar finite-time collapse above criticality and diffusive spreading below it. For k greater than or equal to 4 the unstable steady state separates regimes of collapse and spreading, the collapse dynamics being asymptotically self-similar with the population density exhibiting scale separation that requires a matched-asymptotic desc

What carries the argument

The hydrodynamic free-boundary problem obtained in the N to infinity limit, which governs the evolution of the population density inside a moving boundary set by the farthest particle.

If this is right

  • Explicit steady-state densities exist and can be computed for every integer k.
  • Linear stability holds exactly for k less than or equal to 2 and fails for k greater than or equal to 4.
  • At k=3 the sign of the deviation from the critical reproduction-to-diffusion ratio decides between finite-time collapse and diffusive spreading.
  • For k greater than or equal to 4 the unstable steady state forms the boundary between the two regimes, and collapse inside that regime needs matched asymptotics.
  • Both the hydrodynamic solutions and the microscopic Monte Carlo runs confirm the predicted thresholds and scaling laws.

Where Pith is reading between the lines

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

  • The same free-boundary structure may appear in other lattice models where birth removes a distant particle, allowing the stability thresholds found here to be tested without solving the full microscopic dynamics.
  • If the hydrodynamic description remains valid at moderate N, the critical ratio at k=3 supplies a concrete, parameter-free prediction for the onset of collapse that can be checked in finite-population simulations.
  • The scale separation that appears for k greater than or equal to 4 suggests that similar matched-asymptotic constructions could be needed in other free-boundary problems with strong nonlinearity.
  • The diffusive spreading regime below criticality still carries a quantitative imprint of reproduction, so long-time density profiles in that regime are not purely Gaussian even though the front moves diffusively.

Load-bearing premise

The hydrodynamic free-boundary problem written in the N to infinity limit accurately captures the long-time behavior of the underlying microscopic lattice model with cooperative reproduction kA to (k+1)A.

What would settle it

Monte Carlo simulation of the original lattice model for k=3 at the analytically predicted critical ratio, checking whether the population collapses in finite time or spreads diffusively as the ratio is varied across that value.

Figures

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

Figure 1
Figure 1. Figure 1: Lattice Brownian bees with cooperative reproduction. Particles (circles) on a [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Rescaled steady-state density profile for [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Population dynamics for k = 3 in the three regimes: expansion (α = 4, pan￾els a,b), steady state (α = αc = π 2/2, panels c,d), and collapse (α = 6, panels e,f). (a,c,e) Density profiles u(x, t) in physical coordinates as functions of x at successive times. Dashed gray line (clearly visible on panels a and c) shows the initial condition: the Jacobi cn-steady state (20) at half-width L0 = 100. Colored solid … view at source ↗
Figure 4
Figure 4. Figure 4: Population dynamics for k = 4, starting from rescaled copies of unstable steady￾state profile (14) with different initial half-widths L(0). Dashed gray curve (clearly visible only on panel c): steady-state profile. Insets show L(t): blue band = MC interquartile range (25–75 percentile), black curve = HD, dashed gray = steady-state value LSS ≃ 15.7. (a) Collapse (narrow IC, L(0) = 15.3): HD density profiles… view at source ↗
Figure 5
Figure 5. Figure 5: Self-similar collapse for k = 4 near the blowup time. Solid curves: HD density profiles u(x) at six logarithmically spaced times approaching the blowup time T (from dark to light); dashed curves: matched-asymptotic theory predictions (38) and (39) at the same times. Insets: the half-width of support L (top) and max(u) (bottom) versus τ = T − t on log-log scale; solid = HD, dashed = leading-order scaling L … view at source ↗
Figure 6
Figure 6. Figure 6: Diffusive spreading for k = 4. Shown are (a) Rescaled density profiles, u √ 4πDt vs. x/√ 2Dt: markers = MC (t = 5×103 , 5×104 , 5×105 ; 10 realizations), solid lines = HD numerics (t = 5×104 and 5×105 ); dashed line: the Gaussian (40). Inset: variance σ 2 (t) (circles = MC, solid = HD) vs. the purely diffusive scaling 2Dt (dashed). (b) Boundary position L(t) from HD numerics (solid) compared with theoretic… view at source ↗
read the original abstract

We extend the ``Brownian bees'' model of Berestycki et al. (2021, 2022) to cooperative reproduction, $kA\to(k{+}1)A$, of a population of $N$ symmetric random walkers with removal, at each birth event, of the particle farthest from the origin. Working in the limit $N\to\infty$, we formulate a hydrodynamic free-boundary problem for this model. Using this formalism, we determine steady state population densities for all~$k$ and prove their linear stability for $k\le 2$ and instability for $k\ge 4$. In the marginal case $k=3$, there is a whole continuous family of steady states at a single, critical ratio of the reproduction and diffusion rates. Above criticality the population undergoes an asymptotically self-similar finite-time collapse to the origin. Below the criticality the population spreads diffusively, but the reproduction remains quantitatively relevant. For $k\ge 4$, the unstable steady state separates regimes of a finite-time collapse and a diffusive spreading. Here the collapse dynamics is asymptotically self-similar, and the population density exhibits a scale separation requiring a matched-asymptotic description. Our analytical predictions are confirmed by numerical solutions of the hydrodynamic free-boundary problem and by Monte Carlo simulations of the original microscopic model.

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.

Referee Report

2 major / 2 minor

Summary. The paper extends the Brownian bees model of Berestycki et al. to cooperative reproduction kA → (k+1)A for N symmetric random walkers on a lattice, with removal of the farthest particle at each birth. In the N→∞ limit the authors formulate a hydrodynamic free-boundary problem, derive explicit steady-state densities for every k, prove linear stability for k≤2 and instability for k≥4, and analyze the marginal case k=3 where a continuous family of steady states exists at a critical reproduction-to-diffusion ratio. Above criticality they obtain asymptotically self-similar finite-time collapse; below it they obtain diffusive spreading with quantitatively relevant reproduction. For k≥4 the unstable steady state separates collapse and spreading regimes, the former requiring matched asymptotics. All predictions are checked against numerical solutions of the free-boundary PDE and Monte-Carlo simulations of the microscopic process.

Significance. If the hydrodynamic free-boundary problem is the correct continuum limit, the work supplies a complete classification of steady states and their stability thresholds together with explicit self-similar collapse and spreading solutions. The explicit densities for all k, the sharp stability transition at k=3, and the matched-asymptotic description for k≥4 constitute a substantial analytic advance over the original Brownian-bees results. The dual confirmation by PDE numerics and direct Monte-Carlo sampling of the lattice model is a clear strength.

major comments (2)
  1. [§2] §2 (formulation of the hydrodynamic problem): the central claims—explicit steady-state densities, linear stability for k≤2, instability for k≥4, the critical family at k=3, and the self-similar collapse/spreading regimes—all rest on the free-boundary PDE being the accurate N→∞ limit of the microscopic generator. The manuscript states that the problem is “formulated” in this limit but supplies neither a derivation, tightness argument, nor quantitative error bound between the lattice process and the continuum evolution. This is load-bearing for every subsequent analytic result.
  2. [Stability analysis] Stability analysis (linearized operator around the steady states): the proofs of stability for k≤2 and instability for k≥4 are asserted, yet the explicit form of the linearized free-boundary problem, the function space, and the spectral estimates are not displayed. Without these details it is impossible to verify that the free-boundary condition and the cooperative birth term have been correctly incorporated into the eigenvalue problem.
minor comments (2)
  1. [§3] Notation for the critical ratio at k=3 should be introduced once and used consistently; the current text alternates between “critical ratio” and an unnamed parameter.
  2. [Numerical results] Figure captions for the Monte-Carlo panels should state the number of independent realizations and the lattice size used, to allow direct comparison with the hydrodynamic numerics.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of the analytic contributions and for highlighting the two points that require clarification. Both concerns are valid and will be addressed by expanding the relevant sections in the revised manuscript.

read point-by-point responses

  1. Referee: [§2] §2 (formulation of the hydrodynamic problem): the central claims—explicit steady-state densities, linear stability for k≤2, instability for k≥4, the critical family at k=3, and the self-similar collapse/spreading regimes—all rest on the free-boundary PDE being the accurate N→∞ limit of the microscopic generator. The manuscript states that the problem is “formulated” in this limit but supplies neither a derivation, tightness argument, nor quantitative error bound between the lattice process and the continuum evolution. This is load-bearing for every subsequent analytic result.

    Authors: We agree that an explicit derivation of the hydrodynamic limit strengthens the foundation. In the revision we will add a dedicated subsection to §2 that derives the free-boundary PDE from the microscopic generator, adapting the tightness and convergence arguments of Berestycki et al. (2021) to the cooperative case kA→(k+1)A and indicating the expected rate of convergence in the N→∞ limit. revision: yes

  2. Referee: [Stability analysis] Stability analysis (linearized operator around the steady states): the proofs of stability for k≤2 and instability for k≥4 are asserted, yet the explicit form of the linearized free-boundary problem, the function space, and the spectral estimates are not displayed. Without these details it is impossible to verify that the free-boundary condition and the cooperative birth term have been correctly incorporated into the eigenvalue problem.

    Authors: The linearized operator appears in Section 4, but the function-space setting and the key spectral estimates are only sketched. We will expand this material by adding an appendix that states the precise linearized free-boundary eigenvalue problem (including the cooperative birth term and the moving-boundary condition), specifies the function space, and outlines the spectral estimates used to establish stability for k≤2 and instability for k≥4. revision: yes

Circularity Check

0 steps flagged

No circularity: hydrodynamic formulation used as starting point for independent derivations

full rationale

The paper explicitly states it formulates the hydrodynamic free-boundary problem in the N→∞ limit and then derives steady-state densities, linear stability proofs, and self-similar collapse/spreading regimes from that PDE system. No quoted step reduces a claimed output (e.g., stability thresholds or critical ratio) to a fitted parameter or self-citation by construction. The cited prior Brownian bees work is external and not load-bearing for the new cooperative-reproduction results. Monte Carlo confirmation is reported separately and does not create a closed loop. This is a standard non-circular modeling-plus-analysis structure.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the hydrodynamic limit being valid; no additional free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption The hydrodynamic free-boundary problem in the N→∞ limit accurately captures the microscopic lattice dynamics with cooperative reproduction.
    Invoked immediately when the authors formulate the continuous description from the particle rules.

pith-pipeline@v0.9.1-grok · 5775 in / 1370 out tokens · 24488 ms · 2026-06-28T23:38:47.522445+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

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

  1. [1]

    Brunet, B

    É. Brunet, B. Derrida, Shift in the velocity of a front due to a cutoff, Phys. Rev. E 56 (1997) 2597–2604

    1997

  2. [2]

    Maillard, Speed and fluctuations of N-particle branching Brownian motion with spatial selection, Probab

    P. Maillard, Speed and fluctuations of N-particle branching Brownian motion with spatial selection, Probab. Theory Relat. Fields 166 (2016) 1061–1173

    2016

  3. [3]

    De Masi, P

    A. De Masi, P. A. Ferrari, E. Presutti, N. Soprano-Loto, Hydrodynamics oftheN-BBMprocess, in: StochasticDynamicsOutofEquilibrium, Vol. 282 of Springer Proc. Math. and Statistics, Springer, Cham, 2019, pp. 523–549

    2019

  4. [4]

    Berestycki, É

    J. Berestycki, É. Brunet, B. Derrida, Exact solution and precise asymp- totics of a Fisher–KPP type front, J. Phys. A: Math. Theor. 51 (2018) 035204. 23

    2018

  5. [5]

    Berestycki, É

    J. Berestycki, É. Brunet, J. Nolen, S. Penington, Brownian bees in the infinite swarm limit, Ann. Probab. 50 (2022) 2133–2177

    2022

  6. [6]

    Berestycki, É

    J. Berestycki, É. Brunet, J. Nolen, S. Penington, A free boundary prob- lem arising from branching Brownian motion with selection, Trans. Am. Math. Soc. 374 (2021) 6269–6329

    2021

  7. [7]

    O. Vilk, M. Assaf, B. Meerson, Fluctuations and first-passage properties of systems of Brownian particles with reset, Phys. Rev. E 106 (2022) 024117

    2022

  8. [8]

    R. A. Fisher, The wave of advance of advantageous genes, Ann. Eugenics 7 (1937) 355–369

    1937

  9. [9]

    Kolmogorov, I

    A. Kolmogorov, I. Petrovsky, N. Piskunov, Étude de l’équation de la diffusion avec croissance de la quantité de matière et son application à un problème biologique, Moscow Univ. Math. Bull. 1 (1937) 1–16

    1937

  10. [10]

    van Saarloos, Front propagation into unstable states, Phys

    W. van Saarloos, Front propagation into unstable states, Phys. Rep. 386 (2003) 29–222

    2003

  11. [11]

    Brunet, B

    É. Brunet, B. Derrida, A. H. Mueller, S. Munier, Effect of selection on ancestry: An exactly soluble case and its phenomenological generaliza- tion, Phys. Rev. E 73 (2006) 056126

    2006

  12. [12]

    Meerson, P

    B. Meerson, P. Sasorov, Persistent fluctuations of the swarm size of Brownian bees, Phys. Rev. E 103 (2021) 032140

    2021

  13. [13]

    Siboni, P

    M. Siboni, P. Sasorov, B. Meerson, Fluctuations of a swarm of Brownian bees, Phys. Rev. E 104 (2021) 054131

    2021

  14. [14]

    Sasorov, A

    P. Sasorov, A. Vilenkin, N. R. Smith, Probabilities of moderately atyp- ical fluctuations of the size of a swarm of Brownian bees, Phys. Rev. E 107 (2023) 014140

    2023

  15. [15]

    Benaïm, B

    M. Benaïm, B. Cloez, F. Panloup, Barycentric Brownian bees, Ann. Appl. Probab. 32 (2022) 2504–2539

    2022

  16. [16]

    Meerson, O

    B. Meerson, O. Vilk, Age-structured hydrodynamics of ensembles of anomalously diffusing particles with renewal resetting, Phys. Rev. Res. 8 (2026) 023103. 24

    2026

  17. [17]

    M. B. Miller, B. L. Bassler, Quorum sensing in bacteria, Annu. Rev. Microbiol. 55 (2001) 165–199

    2001

  18. [18]

    D. R. Levitan, Influence of body size and population density on fertil- ization success and reproductive output in a free-spawning invertebrate, Biol. Bull. 181 (1991) 261–268

    1991

  19. [19]

    Clutton-Brock, Breeding together: kin selection and mutualism in cooperative vertebrates, Science 296 (2002) 69–72

    T. Clutton-Brock, Breeding together: kin selection and mutualism in cooperative vertebrates, Science 296 (2002) 69–72

    2002

  20. [20]

    O. Vilk, B. Meerson, Emergent population dynamics of random walkers with cooperative reproduction and spatial selection, arXiv:2605.01770 (2026)

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  21. [21]

    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

    2004

  22. [22]

    F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, M. A. McClain, NIST Digital Library of Mathematical Functions, 2024, release 1.2.1

    2024

  23. [23]

    Courant, D

    R. Courant, D. Hilbert, Methods of Mathematical Physics, Vol. 1, In- terscience, New York, 1953

    1953

  24. [24]

    E. A. Coddington, N. Levinson, Theory of Ordinary Differential Equa- tions, McGraw-Hill, New York, 1955

    1955

  25. [25]

    E. T. Whittaker, G. N. Watson, A Course of Modern Analysis, 4th Edition, Cambridge University Press, Cambridge, 2013

    2013

  26. [26]

    C. M. Bender, S. A. Orszag, Advanced Mathematical Methods for Sci- entists and Engineers: Asymptotic Methods and Perturbation Theory, Springer, New York, 1999

    1999

  27. [27]

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

    1992

  28. [28]

    Virtanen, et al., SciPy 1.0: Fundamental algorithms for scientific computing in Python, Nature Methods 17 (2020) 261–272

    P. Virtanen, et al., SciPy 1.0: Fundamental algorithms for scientific computing in Python, Nature Methods 17 (2020) 261–272

    2020

  29. [29]

    D. T. Gillespie, Exact stochastic simulation of coupled chemical reac- tions, J. Phys. Chem. 81 (1977) 2340–2361. 25

    1977