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

arXiv:2605.29958 · q-bio.PE, cond-mat.stat-mech, math.PR, nlin.PS · Submitted 2026-05-28 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Genomics Radio. Generated commentary on the latest computational biology and genomics papers.

Ines: We're looking at the paper: "Lattice Brownian bees with cooperative reproduction".

Marcus: Here's the gist of the paper. The gist The authors extend a Brownian bees model to cooperative reproduction, determining steady states and dynamical regimes based on the reproduction rate parameter k, finding stability transitions at k=3 and distinct behaviors for k<3, k=3, and k≥4.

Ines: Model Formulation The study extends the “Brownian bees” model of Berestycki et al.

Marcus: (2021, 2022) to cooperative reproduction where each birth event involves the production of one offspring from a group of k individuals in the HD limit N → ∞ The authors formulate a hydrodynamic free-boundary problem for this model using this formalism.

Ines: In the microscopic model, N ≫ 1 random walkers on a one-dimensional lattice are considered, where k co-located particles can undergo the cooperative branching reproduction kA → (k+1)A with rate Λn k/k!

Title and authors: Ines: We just talked about the title and authors of "Lattice Brownian bees with cooperative reproduction: steady states, collapse, and spreading," focusing on how cooperation affects population dynamics based on the parameter k.

Marcus: Yeah, I think it’s important to see this as taking a foundational model—the Brownian bees—and adding the complexity of cooperative reproduction, which is represented by that k value.

Yuki: From my side, I’m just thinking about how these models connect to the broader history of species. It seems like we can use this structure to predict if a population will sustain itself or fragment based on how strongly they cooperate.

Ines: Exactly. The paper shows that the resulting dynamics are highly sensitive to k, going from stable configurations for small k up to distinct collapse or spreading regimes when k is large.

Marcus: And the authors use this framework to find specific transition points, like at k=three where we see a whole family of steady states emerge, which is a key feature they identified <ref:2605.29958#pg1>.

Yuki: So instead of just looking at one outcome, we can map out an entire set of possible macroscopic states for that system when k hits three. That’s pretty powerful for hypothesis testing.

Ines: It’s about understanding the full range of possibilities, not just a single equilibrium point. They also establish the linear stability criteria based on that reproduction and diffusion ratio to tell us which regimes are viable.

Marcus: And they do tie this to the underlying mathematical structure by showing how that mass constraint exponent, (three-k)/two controls the qualitative behavior of those steady states <ref:2605.29958#pg1>.

Yuki: That exponent is a big piece of information because it directly tells us how much weight the reproduction term carries compared to the diffusion term in determining stability.

Ines: So, basically, they’re giving us a mathematical tool to predict if our biological setup, defined by k, will lead to stable coexistence or population collapse.

The paper's summary: Marcus: Now we look at the actual summary of "Lattice Brownian bees with cooperative reproduction: steady states, collapse, and spreading." Essentially, they take a simple branching process on a lattice and analyze it in the hydrodynamic limit.

Ines: They formulate it as a free-boundary reaction-diffusion problem where the term lambda u k introduces that nonlinear local density dependence related to k. This nonlinearity is what fundamentally changes the spatial structure of the population.

Yuki: That nonlinearity is exactly what we see in quorum sensing or cooperative breeding, right? It’s how local density dictates reproduction, which makes this model relevant to those biological contexts.

Marcus: The paper shows that for k two they get robust invasion fronts, meaning traveling-wave solutions, which implies stable population invasion and a pretty predictable spread <ref:2605.29958#pg1>.

Ines: But when k is larger than two, the situation changes dramatically. They find that at k four the steady state is unstable and splits into regimes of finite-time collapse or diffusive spreading <ref:2605.29958#pg1>.

Yuki: So we move from simple invasion to a choice between being squeezed into a singularity or dispersing over time, depending on whether k is small, critical, or large.

Marcus: And for k=four they characterize the collapse dynamics by showing that reproduction dominates in the bulk while diffusion happens in narrow boundary layers near the edges <ref:2605.29958#pg1>.

Ines: That scale separation requires a matched-asymptotic description because you need two different ways to describe what’s happening simultaneously across those vastly different spatial scales.

Yuki: It sounds like this paper provides a roadmap for understanding how cooperative mechanisms translate into macroscopic spatial patterns in populations, whether that's invasion or collapse.

The paper's improvements: Ines: The paper suggests that the real improvement lies in moving beyond just finding these static steady states and actually analyzing the dynamic regimes.

Marcus: Right. They go into detailing the linear stability analysis, which uses eigenvalues of a Schrödinger-type equation to definitively tell us if a steady state is stable or unstable before we even look at time evolution.

Yuki: That’s important because it moves the discussion from just finding where things *can* be stable to predicting which configurations will actually persist long-term in nature.

Ines: And for that marginal case, k=three they suggest a specific analysis of the critical ratio alpha = lambda/D <ref:2605.29958#pg1>. Above that ratio, you get collapse; below it, you get diffusive spreading.

Marcus: So the improvement is moving from just finding a stability boundary to mapping out the precise dynamical regimes based on where those parameters land relative to each other.

Yuki: That gives us a clearer rule: if k is above three, we are likely looking at one of those two extreme behaviors—collapse or spreading—and k=three is the dividing line <ref:2605.29958#pg1>.

Ines: And for k=four they specifically point out that the collapse dynamics exhibit this scale separation, which necessitates that matched-asymptotic description we mentioned earlier <ref:2605.29958#pg1>.

Marcus: So the suggested improvement is to provide a complete dynamical picture: stability boundaries, critical ratios, and specific descriptions of the collapse or spreading dynamics for different values of k.

Conclusion: Ines: So to wrap up on "Lattice Brownian bees with cooperative reproduction: steady states, collapse, and spreading," the main finding is that the value of k completely dictates whether a population self-organizes stably, collapses to a point, or spreads out diffusively.

Marcus: For k less than or equal to two, you get stable macroscopic collective behavior. For k four you don't get that stable configuration at all; it either collapses or spreads diffusively <ref:2605.29958#pg1>.

Yuki: And the critical value k=three sits right at the boundary between those two major dynamical regimes, which is a really neat way to visualize the transition in terms of population structure <ref:2605.29958#pg1>.

Ines: It confirms that cooperative reproduction isn't just adding a small tweak; it introduces entirely new dynamical possibilities that depend heavily on the interaction strength k.

Marcus: This paper gives us a concrete way to predict what happens in systems where group-level reproduction matters, whether we’re looking at biology or something else.

Yuki: It provides a powerful tool for understanding how population structure evolves under different levels of cooperative pressure.

Ines: We’ve seen how this work uses the hydrodynamic limit and the Lane-Emden equation to map out these complex dynamics based on k. That’s what we have today.

Racah Institute of Physics · Hebrew University of Jerusalem

q-bio.PE, cond-mat.stat-mech, math.PR, nlin.PS

Submitted: 2026-05-28

Updated: 2026-10-21

Comments: 26 one-column pages, 6 figures

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 78/100

The gist: The gist The authors extend a Brownian bees model to cooperative reproduction, determining steady states and dynamical regimes based on the reproduction rate parameter k, finding stability

Key concepts

Brownian Bees Model: This is a mathematical model describing random walkers on a lattice where individuals reproduce cooperatively. The key feature is that reproduction depends on the presence of multiple co-located particles, leading to collective growth rather than simple independent movement.
Hydrodynamic Limit (HD Limit): This approximation simplifies the microscopic model by considering a very large number of walkers ($N o ext{infinity}$). It transforms the complex microscopic process into a simpler reaction-diffusion equation that describes the macroscopic behavior of the population density over time.
Linear Stability Analysis: This technique determines if a steady state (a stable configuration where the population size is constant) will persist. By examining how small perturbations grow, researchers can identify critical points (like k=3) where stability changes, predicting whether the system settles into a fixed pattern or moves toward instability.
Finite-Time Collapse vs. Diffusive Spreading: These describe the two possible outcomes when the population becomes unstable (for k>=4). Finite-time collapse means the population rapidly shrinks to a single point, while diffusive spreading means the population disperses outwards over time, governed by diffusion.

Summary

The gist The authors extend a Brownian bees model to cooperative reproduction, determining steady states and dynamical regimes based on the reproduction rate parameter k, finding stability transitions at k=3 and distinct behaviors for k<3, k=3, and k≥4.

Model Formulation

The study extends the “Brownian bees” model of Berestycki et al. (2021, 2022) to cooperative reproduction where each birth event involves the production of one offspring from a group of k individuals in the HD limit N → ∞ The authors formulate a hydrodynamic free-boundary problem for this model using this formalism. In the microscopic model, N ≫ 1 random walkers on a one-dimensional lattice are considered, where k co-located particles can undergo the cooperative branching reproduction kA → (k+1)A with rate Λn k/k! The HD approximation leads to the reaction-diffusion equation ∂tu = D u'' + λ uk.

Steady State Analysis

The steady-state problem is described by a one-dimensional Lane-Emden equation with specific boundary conditions and normalization. For k ≤ 2, the steady state is linearly stable, while at k = 3 there is a whole continuous family of steady states at a single critical ratio of the reproduction and diffusion rates. For k ≥ 4, the steady state is linearly unstable and separates regimes of a finite-time collapse and a diffusive spreading. The exponent (3 − k)/2 in the mass constraint controls the qualitative behavior of the steady states.

Linear Stability Transitions

The linear stability analysis reveals a transition at k = 3, where for k ≤ 2 the steady state is linearly stable. For k > 3, the steady state is linearly unstable and separates regimes of a finite-time collapse and a diffusive spreading. Specifically, for k ≥ 4 the unstable steady state separates regimes of a finite-time collapse and a diffusive spreading. The instability for k > 3 is analytically established via a mass-slope criterion.

Dynamical Regimes

The dynamics are qualitatively different based on the value of k. For k = 3, the system exhibits marginal stability at a single critical value α = λ/D. Above criticality, the population undergoes an asymptotically self-similar finite-time collapse to the origin. Below criticality, the population spreads diffusively, but reproduction remains quantitatively relevant. For k ≥ 4, the unstable steady state separates regimes of a finite-time collapse and a diffusive spreading. The collapse dynamics for k = 4 exhibit a scale separation requiring a matched-asymptotic description.

Collapse and Spreading Dynamics

For k = 4, the collapse dynamics are asymptotically self-similar, and the population density exhibits a scale separation requiring a matched-asymptotic description. In this regime, the reproduction dominates in the bulk while diffusion dominates in narrow boundary layers near x = L(t). For k = 4, the long-time spreading is governed by free diffusion, but with a logarithmic correction to diffusive scaling where L(t) grows as √2Dtln t.

Conclusion

The results confirm that for k ≤ 2 the population self-organizes into a stable macroscopic collective. For k ≥ 4, no such stable configuration exists, as the population either collapses to a point or spreads diffusively. The marginal case k = 3 sits precisely at the boundary between these two regimes.

--- Page 26 ---

References

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

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

[3] A. De Masi, P. A. Ferrari, E. Presutti, N. Soprano-Loto, Hydrodynamics of the N-BBM process, in: Stochastic Dynamics Out of Equilibrium, Vol. 282 of Springer Proc. Math. and Statistics, Springer, Cham

[4] J. Berestycki, É. Brunet, B. Derrida, Exact solution and precise asymptotics of a Fisher–KPP type front, J. Phys. A: Math. Theor. 51 (2018) 035204

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

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

[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

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

[9] 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

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

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

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

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

[14] P. Sasorov, A. Vilenkin, N. R. Smith, Probabilities of moderately atypical fluctuations of the size of a swarm of Brownian bees, Phys. Rev. E 107 (2023) 014140

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

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

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

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

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

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

[21] C. Gardiner, Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences, Springer, Berlin, 2004

[23] R. Courant, D.

Improvements for AI systems

Improved population dynamics modeling for cooperative reproduction systems can predict whether a population will undergo finite-time collapse to a point or diffusive spreading based on the parameter k.
The system can be analyzed to determine if a steady state is linearly stable or unstable by examining the eigenvalues of the linearized Schrödinger-type equation, which is crucial for predicting long-term survival of the population configuration.
For the critical case k equals 3, an AI system can identify if a whole continuous family of steady states exists at a single critical ratio of reproduction and diffusion rates, allowing it to map out all possible macroscopic states.
When k greater than 3, the AI can predict that the steady state is linearly unstable and separate the dynamics into regimes of finite-time collapse or diffusive spreading based on initial conditions.
For k equal to 4, a matched-asymptotic description can be used to characterize the collapse dynamics by separating the bulk region where reproduction dominates from narrow boundary layers where diffusion is dominant.

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 infinity, 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 2 and instability for k 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 4, the unstable steady state separates regimes of a finite-time collapse and 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.

Sources

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