pith. sign in

arxiv: 2604.12797 · v2 · pith:QJFAQ3B5new · submitted 2026-04-14 · 🧮 math.PR · math.AP

A free boundary problem for the mean-field limit of diffusing particles with nonlinear boundary reactivity

Eliana Fausti , Andreas Sojmark This is my paper

Pith reviewed 2026-05-21 01:22 UTC · model grok-4.3

classification 🧮 math.PR math.AP
keywords mean-field limitfree boundary problemdiffusing particlesnonlinear reactivitySkorokhod topologyboundary reactionepidemic modelinert drift
0
0 comments X

The pith

A finite system of diffusing particles with reactive boundaries converges to a unique free boundary problem with nonlinear nonlocal reactivity.

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

The paper examines a system where many particles diffuse inside a domain, reflect at the boundary, and may be killed according to a mechanism that depends on the boundary's current reactivity and each particle's accumulated local time there. Every killing event moves the boundary and updates its reactivity, creating a coupled stochastic evolution. The authors establish that as the number of particles grows, the empirical measure of their positions converges to the solution of a deterministic free boundary problem whose boundary condition is both nonlinear and nonlocal. This limit supplies a rigorous macroscopic description for processes such as epidemic spreading and extends classical fixed-boundary Robin conditions to the moving-boundary setting.

Core claim

The empirical measure flows of the finite-particle system, which incorporate killing through local time at the boundary, converge weakly to the unique solution of a free boundary problem. Inside the time-dependent domain the density satisfies a diffusion equation, while the boundary position evolves according to the integrated effect of a reactivity that depends nonlinearly and nonlocally on the density history at the boundary; this condition reduces to the classical Robin condition when the boundary is fixed and reactivity is constant.

What carries the argument

Skorokhod's M1 topology on path space together with a probabilistic characterisation of near-boundary particle behaviour, used first to identify weak limit points of the empirical measures and then to pass to the limit in the nonlinear boundary condition.

If this is right

  • The microscopic stochastic model yields a well-posed macroscopic description of epidemic spreading driven by boundary reactions.
  • The analysis extends the theory of inert drift systems to the case of moving reactive boundaries.
  • A new mean-field perspective is obtained on encounter-based models for diffusion-mediated surface reactions.
  • Uniqueness of the limit follows from a combination of probabilistic decoupling and energy estimates.

Where Pith is reading between the lines

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

  • The same convergence strategy may apply to systems with multiple interacting boundaries or with killing rates that depend on additional internal particle states.
  • Large-scale particle simulations could be used to benchmark numerical solvers for the nonlinear free boundary problem.
  • The framework suggests analogous mean-field limits for other biological or chemical processes in which particle accumulation at an interface feeds back to boundary motion.

Load-bearing premise

The particles' behaviour near the boundary admits a characterisation that permits identification of weak limit points of the empirical measure flows via Skorokhod's M1 topology.

What would settle it

Numerical simulation of a large finite-particle system in which the observed boundary trajectory and reactivity evolution deviates from the solution of the proposed free boundary problem by more than a small error that vanishes with particle number.

read the original abstract

Consider a finite system of diffusing particles coupled through a reactive boundary. Each particle is reflected, but may react with the boundary according to a killing mechanism which depends on the current reactivity of the boundary and the particle's local time along it. With every such reaction, the boundary moves and its reactivity adjusts. We show that this system admits a unique mean-field limit, described by a free boundary problem with nonlinear and nonlocal reactivity. The latter generalises the classical Robin condition for the case of a fixed boundary with constant reactivity. Via Skorokhod's M1 topology and a characterisation of the particles' behaviour near the boundary, we first identify the weak limit points of the empirical measure flows with killing. Then, we combine a probabilistic decoupling technique and energy estimates to prove uniqueness and deduce convergence. Our analysis gives a rigorous mean-field description of a model of epidemic spreading. Moreover, it contributes to the literature on inert drift systems and yields a novel mean-field perspective on the recent encounter-based framework for diffusion-mediated surface reaction from [Phys. Rev. Lett. 125 (2020) 078102].

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

1 major / 2 minor

Summary. The manuscript considers a finite system of diffusing particles reflected at a boundary but subject to a killing mechanism that depends on the particle's local time and the current reactivity of the boundary. Upon each killing event the boundary moves and its reactivity is updated. The central claim is that the empirical measure of the particle system converges, as the number of particles tends to infinity, to the unique solution of a free-boundary problem whose reactivity is nonlinear and nonlocal, thereby generalizing the classical Robin condition. The proof first identifies weak limit points of the empirical-measure flows via Skorokhod's M1 topology together with a characterization of particle behavior near the boundary, then applies a probabilistic decoupling technique and energy estimates to obtain uniqueness and deduce convergence. The result is presented as a rigorous mean-field description of an epidemic-spreading model and as a contribution to the theory of inert-drift systems and encounter-based surface reactions.

Significance. If the identification and uniqueness arguments hold, the result supplies a rigorous mean-field limit for an interacting particle system with a moving reactive boundary and state-dependent killing, extending classical Robin-type conditions to a nonlinear nonlocal setting. The manuscript explicitly credits the use of Skorokhod's M1 topology to accommodate discontinuities at the boundary and the combination of decoupling with energy estimates to close uniqueness; these are appropriate and technically substantive tools for the problem.

major comments (1)
  1. [Identification of limit points (the paragraph following the statement of the boundary characterization)] The identification of weak limit points (the step that invokes Skorokhod's M1 topology and the characterization of particle behavior near the boundary) relies on a local-time-to-killing map previously established for fixed or linearly reactive boundaries. When reactivity becomes a nonlinear nonlocal functional of the current empirical measure, the effective killing rate seen by each particle depends on the global configuration. The manuscript does not supply an explicit uniform modulus of continuity or tightness argument showing that the same map remains continuous under this feedback, which is required to close the identification before the decoupling-plus-energy-estimates uniqueness argument can be applied.
minor comments (2)
  1. [Abstract] The abstract states that the limit 'generalises the classical Robin condition' but does not indicate the precise functional form of the nonlinearity; a one-sentence description would improve readability.
  2. [Notation and preliminaries] Notation for the empirical measure and the boundary position should be introduced once and used consistently; occasional re-definition of symbols in later sections slows reading.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the constructive comment on the identification of limit points. We agree that the extension of the local-time-to-killing map requires explicit justification of continuity under the nonlinear nonlocal feedback. In the revised manuscript we will add a dedicated lemma establishing the uniform modulus of continuity with respect to the M1 topology on empirical measures, using the Lipschitz property of the reactivity functional and the uniform tightness of local-time processes already available from the energy estimates. This will close the identification step before the decoupling and uniqueness arguments.

read point-by-point responses

  1. Referee: The identification of weak limit points (the step that invokes Skorokhod's M1 topology and the characterization of particle behavior near the boundary) relies on a local-time-to-killing map previously established for fixed or linearly reactive boundaries. When reactivity becomes a nonlinear nonlocal functional of the current empirical measure, the effective killing rate seen by each particle depends on the global configuration. The manuscript does not supply an explicit uniform modulus of continuity or tightness argument showing that the same map remains continuous under this feedback, which is required to close the identification before the decoupling-plus-energy-estimates uniqueness argument can be applied.

    Authors: We thank the referee for this observation. The local-time-to-killing map was indeed established earlier for fixed or linearly reactive boundaries. In the present nonlinear nonlocal setting the reactivity is a continuous functional of the empirical measure, but an explicit uniform modulus is not written out. We will add a new lemma (placed immediately after the boundary-characterization statement) that derives the required modulus of continuity for the killing map. The argument uses the assumed Lipschitz continuity of the reactivity map with respect to the M1 distance on measures together with the uniform integrability of local times that follows from the energy estimates already present in the paper. With this addition the identification of weak limit points proceeds exactly as before, and the subsequent probabilistic decoupling plus energy-estimate uniqueness argument applies without change. We believe the revision will make the extension to the nonlinear case fully rigorous. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation relies on external reference and standard tools

full rationale

The paper's core argument proceeds by first identifying weak limit points of empirical measure flows via Skorokhod M1 topology using a boundary behavior characterization, then applying probabilistic decoupling and energy estimates to obtain uniqueness and convergence to the free-boundary PDE. These steps invoke standard probabilistic machinery together with an external citation to the encounter-based framework in Phys. Rev. Lett. 125 (2020) rather than any self-referential definition, fitted parameter renamed as prediction, or load-bearing self-citation chain. No equation or limit identification reduces to its own inputs by construction, and the cited characterization is treated as an independent external input. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard background results in stochastic analysis and PDE theory together with a domain-specific characterisation of boundary behaviour; no free parameters or new postulated entities appear in the abstract.

axioms (1)
  • domain assumption Existence of weak limit points for the empirical measure flows with killing under Skorokhod's M1 topology.
    Invoked when identifying candidate limits of the particle system.

pith-pipeline@v0.9.0 · 5725 in / 1148 out tokens · 38078 ms · 2026-05-21T01:22:35.423298+00:00 · methodology

Review history (2 revisions) →

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

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

  1. [1]

    C. L. Barnes. Hydrodynamic limit and propagation of chaos for B rownian particles reflecting from a N ewtonian barrier. Ann. Appl. Probab. , 30(4):1582--1613, 2020

    work page 2020

  2. [2]

    Banerjee and B

    S. Banerjee and B. Brown. Inert drift system in a viscous fluid: steady state asymptotics and exponential ergodicity. Trans. Amer. Math. Soc. , 373(9):6369--6409, 2020

    work page 2020

  3. [3]

    Banerjee, K

    S. Banerjee, K. Burdzy, and M. Duarte. Gravitation versus B rownian motion. Ann. Inst. Henri Poincar\'e Probab. Stat. , 55(3):1531--1565, 2019

    work page 2019

  4. [4]

    Banerjee, A

    S. Banerjee, A. Budhiraja, and B. Estevez. The inert drift atlas model. Comm. Math. Phys. , 399(3):2083--2147, 2023

    work page 2083

  5. [5]

    Burdzy, Z.-Q

    K. Burdzy, Z.-Q. Chen, and J. Sylvester. The heat equation and reflected B rownian motion in time-dependent domains. Ann. Probab. , 32(1B):775--804, 2004

    work page 2004

  6. [6]

    Banerjee and D

    S. Banerjee and D. Mukherjee. Join-the-shortest queue diffusion limit in H alfin- W hitt regime: tail asymptotics and scaling of extrema. Ann. Appl. Probab. , 29(2):1262--1309, 2019

    work page 2019

  7. [7]

    Burdzy and D

    K. Burdzy and D. Nualart. Brownian motion reflected on B rownian motion. Probab. Theory Related Fields , 122(4):471--493, 2002

    work page 2002

  8. [8]

    P. C. Bressloff. Diffusion-mediated absorption by partially reactive targets: B rownian functionals and generalized propagators. J. Phys. A: Math. Theor. , 55:205001, 2022

    work page 2022

  9. [9]

    P. C. Bressloff. Spectral theory of diffusion in partially absorbing media. Proc. R. Soc. A , 478:20220319, 2022

    work page 2022

  10. [10]

    P. C. Bressloff. A generalized D ean-- K awasaki equation for an interacting B rownian gas in a partially absorbing medium. Proc. R. Soc. A , 480:20230915, 2024

    work page 2024

  11. [11]

    Baker and M

    G. Baker and M. Shkolnikov. Zero kinetic undercooling limit in the supercooled S tefan problem. Ann. Inst. Henri Poincar\' e Probab. Stat. , 58(2):861--871, 2022

    work page 2022

  12. [12]

    M. T. Barlow and M. Yor. Semimartingale inequalities via the G arsia- R odemich- R umsey lemma, and applications to local times. J. Funct. Anal. , 49(2):198--229, 1982

    work page 1982

  13. [13]

    Dai Pra and F

    P. Dai Pra and F. den Hollander. M c K ean- V lasov limit for interacting random processes in random media. J. Stat. Phys. , 84(3):735--772, 1996

    work page 1996

  14. [14]

    An interacting particle system for the front of an epidemic advancing through a susceptible population

    E. Fausti and A. S jmark. An interacting particle system for the front of an epidemic advancing through a susceptible population. arXiv preprint arXiv:2210.09286 , 2022

    work page internal anchor Pith review Pith/arXiv arXiv 2022

  15. [15]

    D. S. Grebenkov. Paradigm shift in diffusion-mediated surface phenomena. Phys. Rev. Lett. , 125:078102, 2020

    work page 2020

  16. [16]

    D. S. Grebenkov. An encounter-based approach for restricted diffusion with a gradient drift. J. Phys. A: Math. Theor. , 55:045203, 2022

    work page 2022

  17. [17]

    K. P. Hadeler. S tefan problem, traveling fronts, and epidemic spread. Discrete Contin. Dyn. Syst. Ser. B , 21(2):417--436, 2016

    work page 2016

  18. [18]

    Hambly and J

    B. Hambly and J. Meier. M c K ean- V lasov equations with positive feedback through elastic stopping times. Electron. Commun. Probab. , 27:Paper No. 41, 13, 2022

    work page 2022

  19. [19]

    Hambly and A

    B. Hambly and A. S jmark. An SPDE model for systemic risk with endogenous contagion. Finance Stoch. , 23(3):535--594, 2019

    work page 2019

  20. [20]

    K. It \^o . Distribution-valued processes arising from independent B rownian motions. Math. Z. , 182(1):17--33, 1983

    work page 1983

  21. [21]

    F. B. Knight. On the path of an inert object impinged on one side by a B rownian particle. Probab. Theory Related Fields , 121(4):577--598, 2001

    work page 2001

  22. [22]

    S. Ledger. S korokhod's M_1 topology for distribution-valued processes. Electron. Commun. Probab. , 21(34):1--11, 2016

    work page 2016

  23. [23]

    H. P. McKean. B rownian local times. Adv. Math. , 16:91--111, 1975

    work page 1975

  24. [24]

    I. Mitoma. Tightness of probabilities on C([0,1]; S ) and D([0,1]; S ) . Ann. Probab. , 11(4):989--999, 1983

    work page 1983

  25. [25]

    Pilipenko

    A. Pilipenko. An introduction to stochastic differential equations with reflection . Lectures in Pure and Applied Mathematics. Potsdam University Press, 2014

    work page 2014

  26. [26]

    Z. Qian, F. Russo, and W. Zheng. Comparison theorem and estimates for transition probability densities of diffusion processes. Probab. Theory Related Fields , 127(3):388--406, 2003

    work page 2003

  27. [27]

    Qian and X

    Z. Qian and X. Xu. Probability bounds for reflecting diffusion processes. Statist. Probab. Lett. , 199:109855, 2023

    work page 2023

  28. [28]

    Qian and W

    Z. Qian and W. Zheng. A representation formula for transition probability densities of diffusions and applications. Stochastic Process. Appl. , 111(1):57--76, 2004

    work page 2004

  29. [29]

    Revuz and M

    D. Revuz and M. Yor. Continuous martingales and B rownian motion , volume 293 of Grundlehren der mathematischen Wissenschaften . Springer-Verlag, Berlin, third edition, 1999

    work page 1999

  30. [30]

    S jmark and F

    A. S jmark and F. Wunderlich. Weak convergence of stochastic integrals on S korokhod space in S korokhod's J_1 and M_1 topologies. Probab. Theory Related Fields , 2026

    work page 2026

  31. [31]

    Wang and Y

    R. Wang and Y. Du. Long-time dynamics of a diffusive epidemic model with free boundaries. Discrete Contin. Dyn. Syst. Ser. B , 26(4):2201--2238, 2021

    work page 2021

  32. [32]

    Wang and Y

    R. Wang and Y. Du. Long-time dynamics of a nonlocal epidemic model with free boundaries: S preading-vanishing dichotomy. J. Differential Equations , 327:322--381, 2022

    work page 2022

  33. [33]

    D. White. Processes with inert drift. Electron. J. Probab. , 12:1509--1546, 2007

    work page 2007

  34. [34]

    Z. Wang, H. Nie, and Y. Du. Spreading speed for a W est N ile virus model with free boundary. J. Math. Biol. , 79(2):433--466, 2019

    work page 2019