arxiv: 2605.04860 · v1 · submitted 2026-05-06 · 🧮 math.PR

Recognition: unknown

Branching Brownian motion with rank-based selection and reaction-diffusion equations

Jacob Mercer

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

classification 🧮 math.PR

keywords branching Brownian motionrank-based selectionhydrodynamic limitreaction-diffusion equationasymptotic velocitytravelling wavesweak selection principle

0 comments

The pith

Branching Brownian motion with rank-based selection converges to a reaction-diffusion equation with nonlinearity given by the killing function.

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

This paper studies branching Brownian motion where particles branch at time-dependent rate r and are killed based on their rank according to a function ψ. It proves that under fairly minimal conditions the empirical measure of the particles converges to the solution of the reaction-diffusion PDE U_t = 1/2 U_xx + r(t) G(U), where G is derived from ψ. This generalizes the N-BBM process and similar models, enabling the transfer of asymptotic results between the particle system and the PDE. A sympathetic reader would care because it provides a microscopic justification for certain reaction-diffusion models and a way to analyze selection in large populations via PDE techniques.

Core claim

We show that the hydrodynamic limit of the branching-selection particle system is the reaction-diffusion equation U_t = ½ U_xx + r(t)G(U) with G(U) a function of ψ. General conditions are given under which the system has an asymptotic velocity described up to order (log N)^{-2}, and this velocity connects to the spreading speeds and travelling waves of the PDE, providing a partial weak selection principle.

What carries the argument

The killing function ψ that depends on particle rank and determines the form of the nonlinearity G in the limiting reaction-diffusion equation.

Load-bearing premise

The result relies on the killing function ψ and branching rate r(t) satisfying fairly minimal but unspecified conditions that guarantee the hydrodynamic limit holds.

What would settle it

Numerical simulation of the particle system for large N with a chosen ψ showing that the front position does not match the predicted speed from the PDE travelling wave, or a mathematical counterexample where the limit PDE differs from the claimed form.

Figures

Figures reproduced from arXiv: 2605.04860 by Jacob Mercer.

**Figure 1.** Figure 1: Functions G and the corresponding selection functions ψ for Examples 1, 2, & 3. Example 4: Let us consider another way in which the PDE connection may be insightful. Consider a branching-selection particle system with rank-dependent selection such that particles of intermediate rank are deleted, whereas particles of extreme ranks (ie. close to 0 or N) are never deleted. What should we expect the behaviour … view at source ↗

**Figure 2.** Figure 2: Selection function ψ(x) = max{0, K(x − 0.8)(0.3 − x)} and corresponding G(x). 18 view at source ↗

read the original abstract

We consider a family of branching-selection particle systems in which particles branch at time dependent rate $r$ and are killed with a probability which is dependent on their rank via some function $\psi$. We show that, under fairly minimal conditions, the hydrodynamic limit of such a system is given by the reaction-diffusion equation $U_t = \frac12 U_{xx} + r(t)G(U)$ with nonlinearity $G(U)$ which is a function of $\psi$. This is a significant generalisation of the well-studied $N$-BBM process, and is similar to the family of `$(b,D)$-BBM' processes described by Groisman \& Soprano-Loto (arXiv:2008.09460). On the one hand, this allows us to understand common reaction-diffusion equations as limits of interacting particle systems with simple descriptions. On the other hand, the asymptotic behaviour of solutions of the reaction-diffusion PDEs can help us predict the asymptotic properties of the associated particle systems. We give general conditions under which the branching-selection particle system has an asymptotic velocity, and describe the velocity up to order $(\log N)^{-2}$; furthermore, we describe the connection between this velocity and the spreading speeds and travelling waves of the corresponding reaction-diffusion equation. This provides a partial weak selection principle.

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 generalizes N-BBM to arbitrary rank-based killing via ψ and links the resulting particle velocity to the PDE traveling-wave speed under explicit conditions.

read the letter

The paper extends the fixed-N selection in N-BBM to a broader family where particles are killed according to their rank through a function ψ. Under stated conditions on ψ and the branching rate r(t), the empirical measure converges to the solution of the reaction-diffusion equation U_t = ½ U_xx + r(t) G(U), with G determined by ψ. It also gives an asymptotic velocity for the particle system up to order (log N)^{-2} and relates that velocity to the spreading speed of the corresponding PDE traveling wave. This supplies a partial weak selection principle and shows how simple rank rules can produce a range of nonlinearities in the limit. The stress-test note confirms that monotonicity, bounded variation, and integrability requirements on ψ are written into the theorems, and the proof uses standard tightness plus generator identification without hidden steps. The velocity asymptotics rest on applying known PDE results to the limit equation rather than deriving new PDE estimates, which keeps the contribution focused on the particle side. The argument structure is clean and the mapping from ψ to G is the main new piece relative to the cited N-BBM and (b,D)-BBM work. This is useful for people working on branching processes or hydrodynamic limits who need flexible selection mechanisms. It is worth sending to peer review because the claims are grounded and the generalization is concrete.

Referee Report

2 major / 3 minor

Summary. The manuscript studies branching Brownian motion where particles branch at time-dependent rate r(t) and are killed according to a rank-dependent function ψ. It proves that, under explicit conditions on ψ (monotonicity, bounded variation, integrability) and r (continuous and positive), the empirical measure converges in probability to the solution of the reaction-diffusion equation U_t = ½ U_xx + r(t) G(U), with nonlinearity G derived from ψ. The paper further establishes an asymptotic velocity for the particle system up to order (log N)^{-2} and connects this velocity to the spreading speed and traveling waves of the limiting PDE, yielding a partial weak selection principle. This generalizes the N-BBM process and the (b,D)-BBM family.

Significance. If the hydrodynamic limit and velocity results hold, the work supplies a microscopic particle-system foundation for a broad class of reaction-diffusion equations with rank-based nonlinearities, enabling transfer of techniques between the two settings. The explicit velocity expansion and link to PDE traveling-wave speeds constitute a concrete advance beyond existing N-BBM analyses, while the proof via tightness and generator identification under minimal assumptions is a methodological strength.

major comments (2)

[§2.2, Theorem 2.3] §2.2, Theorem 2.3: the identification of limit points in the mild form of the PDE relies on the generator of the rank-based killing; the argument that the test-function integral converges to the integral against G(U) appears to require an additional uniform integrability step that is not explicitly verified under the stated bounded-variation assumption on ψ.
[§4.1, Proposition 4.2] §4.1, Proposition 4.2: the velocity expansion up to (log N)^{-2} is obtained by combining the hydrodynamic limit with standard traveling-wave analysis; however, the error term inherited from the hydrodynamic approximation is only controlled to o(1), which is insufficient to justify the claimed second-order term without a quantitative rate in the convergence of the empirical measure.

minor comments (3)

[Abstract] The abstract states 'fairly minimal conditions' on ψ and r(t); these should be restated verbatim in the introduction or in the statement of the main theorems for immediate readability.
[§1 and §3] Notation for the empirical measure and the rank function is introduced in §1 but reused with slight variations in §3; a single consolidated definition would reduce ambiguity.
[Introduction] The comparison with Groisman & Soprano-Loto (arXiv:2008.09460) is mentioned but lacks a precise statement of which features are new versus recovered.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable comments on our manuscript. We address the two major comments point by point below, indicating the revisions we plan to make.

read point-by-point responses

Referee: [§2.2, Theorem 2.3] §2.2, Theorem 2.3: the identification of limit points in the mild form of the PDE relies on the generator of the rank-based killing; the argument that the test-function integral converges to the integral against G(U) appears to require an additional uniform integrability step that is not explicitly verified under the stated bounded-variation assumption on ψ.

Authors: We appreciate this observation. The proof of Theorem 2.3 does rely on identifying the limit via the generator, and while the bounded variation of ψ ensures the necessary regularity for G, an explicit uniform integrability step is indeed helpful to rigorously pass to the limit in the integral term. We will revise the manuscript to include this verification, using the integrability condition on ψ to bound the relevant terms and apply Vitali's convergence theorem or an equivalent argument in the space of measures. revision: yes
Referee: [§4.1, Proposition 4.2] §4.1, Proposition 4.2: the velocity expansion up to (log N)^{-2} is obtained by combining the hydrodynamic limit with standard traveling-wave analysis; however, the error term inherited from the hydrodynamic approximation is only controlled to o(1), which is insufficient to justify the claimed second-order term without a quantitative rate in the convergence of the empirical measure.

Authors: This is a fair point regarding the precision needed for the second-order term. Our current proof combines the hydrodynamic limit with traveling-wave analysis, and the o(1) error is absorbed in the logarithmic correction due to the specific structure of the front propagation. Nevertheless, to strengthen the argument, we will include a more detailed error analysis in the revised manuscript, possibly deriving a rate under additional assumptions or clarifying why the o(1) suffices for the claimed order. revision: partial

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper establishes the hydrodynamic limit from the rank-based branching-selection particle system to the reaction-diffusion PDE via tightness of the empirical measure, generator identification of limit points, and passage to the mild form of the PDE under explicitly stated conditions on the killing function ψ and rate r(t). This constitutes a forward derivation from microscopic dynamics to macroscopic equation rather than any reduction of outputs to inputs by construction. The asymptotic velocity description up to order (log N)^{-2} and its link to PDE spreading speeds follow by applying standard, independent traveling-wave analysis to the limiting PDE; no parameters are fitted from the particle system and then renamed as predictions, no self-definitional steps appear, and no load-bearing self-citations are invoked. The argument is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper rests on standard tools of stochastic analysis and PDE theory with no new free parameters or invented entities.

axioms (2)

standard math Tightness and convergence of empirical measures for branching particle systems
Required to pass from the microscopic process to the hydrodynamic PDE limit.
domain assumption Existence of spreading speeds and traveling waves for the reaction-diffusion equation with the derived nonlinearity G
Used to identify the asymptotic velocity of the particle system.

pith-pipeline@v0.9.0 · 5521 in / 1274 out tokens · 55514 ms · 2026-05-08T16:51:37.727227+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

22 extracted references · 3 canonical work pages

[1]

A weak formulation of free boundary problems and its application to hydrodynamic limits of particle systems with selection.Ann

Rami Atar. A weak formulation of free boundary problems and its application to hydrodynamic limits of particle systems with selection.Ann. Probab., 53(5):1825–1856, 2025. 26

2025
[2]

PhD thesis, Duke University, 2019

Erin M Beckman.Asymptotic Behaviour of Certain Branching Processes. PhD thesis, Duke University, 2019

2019
[3]

Selection principle for theN-BBM.https://arxiv.org/abs/ 2407.05792, 2024.arXiv:2407.05792

Julien Berestycki and Oliver Tough. Selection principle for theN-BBM.https://arxiv.org/abs/ 2407.05792, 2024.arXiv:2407.05792

work page arXiv 2024
[4]

Billingsley.Convergence of probability measures

Patrick. Billingsley.Convergence of probability measures. Wiley series in probability and mathematical statistics. Wiley, New York, 1968

1968
[5]

V. I. Bogachev.Measure theory. Springer, Berlin, 2007

2007
[6]

Ferrari, Errico Pesutti, and Nahuel Soprano-Loto

Anna De Masi, Pablo A. Ferrari, Errico Pesutti, and Nahuel Soprano-Loto. Hydrodynamics of the N-BBM Process. InStochastic Dynamics Out of Equilibrium, volume 282 ofSpringer Proceedings in Mathematics & Statistics, pages 523–549. Springer International Publishing AG, Switzerland, 2019

2019
[7]

Convergence and wave propagation for a system of branching rank-based interacting brownian particles.https://arxiv.org/abs/2505.08563, 2025

Mete Demircigil and Milica Tomasevic. Convergence and wave propagation for a system of branching rank-based interacting brownian particles.https://arxiv.org/abs/2505.08563, 2025. arXiv:2505.08563

work page arXiv 2025
[8]

Existence and convergence to a propagating terrace in one-dimensional reaction-diffusion equations.Transactions of the American Mathematical Society, 366(10):5541–5566, 2014

Arnaud Ducrot, Thomas Giletti, and Hiroshi Matano. Existence and convergence to a propagating terrace in one-dimensional reaction-diffusion equations.Transactions of the American Mathematical Society, 366(10):5541–5566, 2014

2014
[9]

Etheridge.An introduction to superprocesses

Alison. Etheridge.An introduction to superprocesses. University lecture series; 20. American Mathe- matical Society, Providence, R.I, 2000

2000
[10]

Ethier and Thomas G

Stewart N. Ethier and Thomas G. Kurtz.Markov processes : characterization and convergence. Wiley series in probability and mathematical statistics. Wiley, New York, first edition, 1986

1986
[11]

R. A. Fisher. The wave of advantageous genes.Annals of Eugenics, 7(4):355–369, 1937

1937
[12]

Existence and uniqueness of propagating terraces.Communications in contemporary mathematics, 22(6):1950055–, 2020

Thomas Giletti and Hiroshi Matano. Existence and uniqueness of propagating terraces.Communications in contemporary mathematics, 22(6):1950055–, 2020

2020
[13]

Front propagation and quasi-stationary distributions: Two faces of the same coin

Pablo Groisman and Matthieu Jonckheere. Front propagation and quasi-stationary distributions: Two faces of the same coin. InSojourns in Probability Theory and Statistical Physics - III, volume 300 of Springer Proceedings in Mathematics & Statistics, pages 242–254. Springer, 2019

2019
[14]

F-KPP scaling limit and selection principle for a Brunet-Derrida type particle system.ALEA Lat

Pablo Groisman, Matthieu Jonckheere, and Juli´ an Mart´ ınez. F-KPP scaling limit and selection principle for a Brunet-Derrida type particle system.ALEA Lat. Am. J. Probab. Math. Stat., 17(1):589–607, 2020

2020
[15]

Rank dependent branching-selection particle systems.Elec- tronic Journal of Probability, 26, 2021

Pablo Groisman and Nahuel Soprano-Loto. Rank dependent branching-selection particle systems.Elec- tronic Journal of Probability, 26, 2021

2021
[16]

Spreading properties and complex dynamics for monostable reaction-diffusion equations.Communications in partial differential equations, 37(3):511–537, 2012

Fran¸ cois Hamel and Gr´ egoire Nadin. Spreading properties and complex dynamics for monostable reaction-diffusion equations.Communications in partial differential equations, 37(3):511–537, 2012

2012
[17]

Kolmogorov, I

A. Kolmogorov, I. Petrovskii, and N. Piskunov. A study of the diffusion equation with increase in the amount of substance, and its application to a biological problem. In V. M. Tokhomirov, editor,Selected Works of A. N. Kolmogorov I, chapter 38, pages 224–270. Kluwer, 1991. 27

1991
[18]

Speed and fluctuations of n-particle branching brownian motion with spatial selection

Pascal Maillard. Speed and fluctuations of n-particle branching brownian motion with spatial selection. Probability theory and related fields, 166(3-4):1061–1173, 2016

2016
[19]

Critical Drift for Brownian Bees and a Reflected Brownian Motion Invariance Principle

Jacob Mercer. Critical Drift for Brownian Bees and a Reflected Brownian Motion Invariance Principle. https://arxiv.org/abs/2412.04527, 2024.arXiv:2412.04527

work page arXiv 2024
[20]

A criterion of convergence of measure-valued processes: application to mea- sure branching processes.Stochastics, 17(1-2):43–65, 1986

Sylvie Roelly-Coppoletta. A criterion of convergence of measure-valued processes: application to mea- sure branching processes.Stochastics, 17(1-2):43–65, 1986

1986
[21]

Schilling

Ren´ e L. Schilling. When does a c` adl` ag process have continuous sample paths?Exposition. Math., 12(3):255–261, 1994

1994
[22]

A. I. Volpert, Vitaly A. Volpert, Vladimir A. Volpert, and American Mathematical Society.Travel- ing wave solutions of parabolic systems. Translations of mathematical monographs, v. 140. American Mathematical Society, Providence, R.I, 1994. 28

1994