From stochastic individual-based models to free-boundary Hamilton-Jacobi equations
Pith reviewed 2026-05-10 18:08 UTC · model grok-4.3
The pith
Stochastic branching models for trait-structured populations converge to free-boundary Hamilton-Jacobi equations with state constraints that enforce local extinctions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In the regime of large populations and small mutations analyzed in logarithmic scales, the stochastic individual-based system converges to a class of free-boundary Hamilton-Jacobi equations with state constraints. These equations go beyond classical Hamilton-Jacobi equations from deterministic models by accounting for possible extinction of the system in certain regions of the trait space. The derivation combines analysis of Hamilton-Jacobi equations with large deviation principles and tools from branching process theory.
What carries the argument
Free-boundary Hamilton-Jacobi equations with state constraints, obtained by large-deviation analysis of the branching process, which enforce zero population density outside the surviving trait regions.
If this is right
- The limiting population density must be zero in any trait region where extinction is favored, creating a sharp free boundary.
- The support of the limiting population distribution is determined by the solution of the constrained equations rather than filling the entire trait space.
- Evolutionary outcomes can be predicted with explicit inclusion of stochastic die-off in suboptimal trait regions.
- The macroscopic model distinguishes persistence from extinction more accurately than interior Hamilton-Jacobi equations alone.
Where Pith is reading between the lines
- The same scaling approach may apply to models with different mutation kernels or interaction structures, potentially yielding similar free boundaries.
- Numerical checks of the limiting support against finite-population simulations could quantify how well the state constraint approximates extinction for moderately large populations.
- This limit suggests that deterministic models without state constraints may overestimate long-term survival in marginal trait regions.
Load-bearing premise
The analysis requires a regime of large population size and small mutations, taken specifically in logarithmic scales of size and time.
What would settle it
A numerical simulation of the individual-based branching model in a parameter regime with extinction in part of the trait space, followed by comparison of the simulated population support against the zero level set of the limiting equation, would test whether the free boundary appears as predicted.
read the original abstract
We study a stochastic branching model for a population structured by a quantitative phenotypic trait and subject to births, deaths, and mutations. In a regime of large population and small mutations, and in logarithmic scales of size and time, we derive a certain class of free boundary Hamilton-Jacobi equations with state constraints from the stochastic individual-based system. This goes beyond the classical Hamilton-Jacobi equations obtained from deterministic models by taking into account the possible extinction of the system in certain regions of the trait space. The proof is obtained by combining methods for the analysis of Hamilton-Jacobi equations with probabilistic tools from the theory of large deviations and branching processes.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives a class of free-boundary Hamilton-Jacobi equations with state constraints from a stochastic individual-based branching model for a population structured by a quantitative phenotypic trait. In the regime of large population and small mutations, on logarithmic scales of size and time, the authors combine methods from Hamilton-Jacobi equations with large deviations and branching processes to obtain the limit, accounting for extinction in certain trait regions beyond classical deterministic models.
Significance. If the derivation holds, this result is significant as it provides a rigorous passage from stochastic microscopic models to macroscopic free-boundary PDEs that incorporate the possibility of local extinction, which is a key feature not captured by deterministic approximations. The probabilistic approach using large deviations and branching processes adds robustness to the analysis in the logarithmic regime.
minor comments (2)
- [Abstract] The abstract describes the proof strategy at a high level; adding one sentence on the main technical steps (e.g., how the state constraint arising from extinction is obtained via branching-process estimates) would improve accessibility.
- Notation for the trait space, the free boundary, and the logarithmic scaling should be introduced with a short table or explicit list in the introduction to aid readers following the large-deviation arguments.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our manuscript, the accurate summary of its contributions, and the recommendation for minor revision. The referee correctly identifies the key novelty in incorporating extinction effects via the stochastic branching model and the use of large-deviation and branching-process techniques in the logarithmic regime.
Circularity Check
Derivation from stochastic model to PDE limit is self-contained
full rationale
The paper starts from an explicit stochastic individual-based branching process with births, deaths, and mutations. It invokes standard large-deviation and branching-process tools to pass to the logarithmic large-population/small-mutation limit, obtaining a free-boundary Hamilton-Jacobi equation with state constraints. No equation or step is shown to be presupposed by the target PDE, no parameter is fitted to data and then relabeled a prediction, and no load-bearing uniqueness result is imported solely via self-citation. The derivation therefore remains independent of its own output.
Axiom & Free-Parameter Ledger
axioms (4)
- domain assumption The population is large
- domain assumption Mutations are small
- domain assumption Logarithmic scaling in size and time
- domain assumption The model is a branching process with births, deaths, mutations
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel (J uniqueness); branch_selection (coupling excludes additive branch) echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
ua(t,x)=sup{Ft(f):f∈AC[0,t],f(t)=x,∀s Fs(f)≥a} ... state-constrained Hamilton-Jacobi ... Ωa open ... ∂tu=p(x)H(∂xu)+R(x) on Ωa, u=a on ∂Ωa
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability; bare_distinguishability_of_absolute_floor echoes?
echoesECHOES: this paper passage has the same mathematical shape or conceptual pattern as the Recognition theorem, but is not a direct formal dependency.
only paths fs such that Fs(f)>0 for all s∈[0,t] are admissible: the population gets extinct on the way otherwise
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
V. Bansaye, J.F. Delmas, L. Marsalle and V.C. Tran. Limit theorems for Markov processes indexed by continuous time Galton-Watson trees. The Annals of Applied Probability, Vol. 21, No. 6, 2263-2314, (2011)
work page 2011
-
[2]
G. Barles. Solutions de viscosit´ edes ´ equationsde Hamilton-Jacobi. Springer-Verlag Berlin Heidelberg, 1994. [3]G. Barles,An Introduction to the Theory of Viscosity Solutions for First-Order Hamilton-Jacobi Equations and Applications,Lecture Notes in Mathematics 2074, 2013
work page 1994
- [3]
- [4]
-
[5]
J. Berestycki and E. Brunet and J.W. Harris and J. W. Harris and C. Simon and M.I. Roberts. Growth rates of the population in a branching Brownian motion with an inhomogeneous breeding potential. Stochastic Process. Appl. 125 (2015), no. 5, 2096–2145
work page 2015
-
[6]
J.D. Biggins. The growth and spread of the general branching random walk. The Annals of Applied Probability, 5(4):1008-1024, (1995)
work page 1995
-
[7]
P. Billingsley. Convergence of Probability Measures. John Wiley & Sons, New York (1968)
work page 1968
- [8]
-
[9]
V. Calvez, B. Henry, S. M´ el´ eard, V.C. Tran. Dynamics of lineages in adaptation to a gradual environmental change. Annales Henri Lebesgue, 5:729–777, 2022. [11]V. Calvez, V. and K.-Y. Lam,Uniqueness of the viscosity solution of a constrained Hamilton-Jacobi equation, Calc. Var. Partial Differ. Equ.,59(5)(2020) pp. 163
work page 2022
-
[10]
N. Champagnat. A microscopic interpretation for adaptative dynamics trait substitution sequence models. Stochastic Processes and their Applications, 116:1127–1160, 2006
work page 2006
-
[11]
N. Champagnat, R. Ferri` ere, and S. M´ el´ eard. Unifying evolutionary dynamics: from individual stochastic processes to macroscopic models via timescale separation. Theoretical Population Biology, 69:297–321, 2006
work page 2006
-
[12]
N. Champagnat and S. M´ el´ eard. Polymorphic evolution sequence and evolutionary branching. Probability Theory and Related Fields, 151(1-2):45–94, 2011
work page 2011
-
[13]
N. Champagnat, S. M´ el´ eard and V.C. Tran. Stochastic analysis of emergence of evolutionary cyclic behavior in population dynamics with transfer. Ann. Appl. Probab., 31(4), 1820–1867, 2021
work page 2021
-
[14]
N. Champagnat, S. M´ el´ eard, S. Mirrahimi and V.C. Tran. Filling the gap between individual-based evolu- tionary models and Hamilton-Jacobi equations. J. Ec. Polytechnique, 10, 1247–1275, 2023
work page 2023
-
[15]
L. Coquille and A. Kraut and C. Smadi. Stochastic individual-based models with power law mutation rate on a general finite trait space. Electron. J. Probab., 26, 1–37, 2021. [18]G. Dal Maso, H. Frankowska,Value functions for Bolza problems with discontinuous Lagrangians and Hamilton-Jacobi inequalities, ESAIM: Control, Optimisation and Calculus of Variati...
work page 2021
-
[16]
A. Dembo and O. Zeitouni. Large Deviations Techniques and Applications. Vol 38., Springer, second edition, 1998
work page 1998
-
[17]
O. Diekmann, P.-E. Jabin, S. Mischler, and B. Perthame. The dynamics of adaptation: an illuminating example and a Hamilton-Jacobi approach. Theoretical Population Biology, 67, 257–271, 2005
work page 2005
-
[18]
P. Dupuis and R.S. Ellis. A Weak Convergence Approach to the Theory of Large Deviations. Wiley Series in Probability and Statistics, 1997
work page 1997
-
[19]
R. Durrett and J. Mayberry. Travelling waves of selective sweeps. Annals of Applied Probability, 21(2), 699–744, 2011
work page 2011
-
[20]
M. Esser and A. Kraut. A general multi-scale description of metastable adaptive motion across fitness valleys. J. Math. Biol., 89(46), 2024
work page 2024
-
[21]
Stewart N. Ethier and Thomas G, Kurtz. Markov processes: characterization and convergence. John Wiley & Sons, 1986
work page 1986
-
[22]
L. C. Evans Partial differential equations Graduate Studies in Mathematics Vol. 19, American Mathematical Society, 1998
work page 1998
-
[23]
L. C. Evans and P. E. Souganidis. A PDE approach to geometric optics for certain semilinear parabolic equations. Indiana Univ. Math. J., 38(1):141–172, 1989
work page 1989
-
[24]
M. Fang and O. Zeitouni. Slowdown for Time Inhomogeneous Branching Brownian Motion. J. Statistical Physics, 149:1-9, (2012)
work page 2012
-
[25]
A. Fathi. Weak Kam Theorem in Lagrangian Dynamics. Cambridge Studies in Advanced Mathematics, Cambridge University Press, 2016
work page 2016
-
[26]
Limit theorems for large deviations and reaction-diffusion equations
M Freidlin. Limit theorems for large deviations and reaction-diffusion equations. The Annals of Probability, 13(3):639–675, 1985
work page 1985
- [27]
-
[28]
N. Fournier and S. M´ el´ eard. A Microscopic Probabilistic Description of a Locally Regulated Population and Macroscopic Approximations. Ann. Appl. Probab., 14(4):1880-1919, 2004
work page 1919
-
[29]
R. Hardy and S.C. Harris. A Spine Approach to Branching Diffusions with Applications to Lp-Convergence of Martingales. S´ eminairede Probabilit´ esXLII, Lecture Notes in Mathematics, 1979:281-330, Springer Nature, 2009
work page 1979
- [30]
-
[31]
N. Ikeda and S. Watanabe. Stochastic Differential Equations and Diffusion Processes, North-Holland Pub- lishing Company, 1989
work page 1989
-
[32]
P.-E. Jabin. Small populations corrections for selection-mutation models. Netw. Heterog. Media, 7(4), 805–836, 2012
work page 2012
- [33]
-
[34]
G. Last and M. Penrose. Lectures on the Poisson Process. IMS Textbook by Cambridge University Press (2017)
work page 2017
-
[35]
A. Lorz, S. Mirrahimi and B Perthame. Dirac mass dynamics in a multidimensional nonlocal parabolic equation. Communications in Partial Differential Equations, 36, 1071–1098 (2011)
work page 2011
-
[36]
P. Maillard and G. Raoul and J. Tourniaire. Spreading speed of locally regulated population models in macroscopically heterogeneous environments (2024). ArXiv:2105.06985
-
[37]
B. Mallein. Maximal displacement of a branching random walk in time- inhomogeneous environment. Stochastic Processes and Their Applications, 125(10), 3958–4019 (2015)
work page 2015
-
[38]
A. Marguet. Uniform sampling in a structured branching population. Bernoulli, 25, 4A, 2649-2695, 2019
work page 2019
-
[39]
S. M´ el´ eard and V.C. Tran. Nonlinear historical superprocess approximations for population models with past dependence. Electronic Journal of Probability, 17(47):1-32, 2012
work page 2012
-
[40]
S. Mirrahimi, G. Barles, B. Perthame, P.E. Souganidis. A singular Hamilton-Jacobi equation modeling the tail problem. SIAM J. Math. Anal., 44 (6), 4297–4319, 2012
work page 2012
-
[41]
B. Perthame and G. Barles. Dirac concentrations in Lotka-Volterra parabolic PDEs. Indiana Univ. Math. J. 57, 3275–3301 (2008)
work page 2008
-
[42]
B. Perthame and M. Gauduchon. Survival thresholds and mortality rates in adaptive dynamics: conciliating deterministic and stochastic simulations. Math. Med. Biol. 27(3), 195–210 (2010)
work page 2010
-
[43]
Protter.Stochastic integration and differential equations, second Edition
P.E. Protter.Stochastic integration and differential equations, second Edition. Stochastic Modelling and Applied Probability, 21, Springer, Berlin, 2004
work page 2004
-
[44]
C. Desmarais, E. Schertzer, Z. Talygi´ as. K-Branching random walk with noisy selection: large population limits and phase transitions. arXiv:2509.26254
-
[45]
D. Waxman and S. Gavrilets. 20 Questions on Adaptive Dynamics. Journal of Evolutionary Biology 18, 1139-1154 (2005)
work page 2005
- [46]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.