Quantitative approximation of a Keller--Segel PDE by a branching moderately interacting particle system and suppression of blow-up
Thomas Cavallazzi, Alexandre Richard, Milica Tomasevic
TL;DR
The paper analyzes the Keller–Segel chemotaxis model with logistic damping on $\mathbb{R}^d$ and proves global well-posedness under a damping threshold, with special suppression of blow-up in $d=2$. It then constructs a microscopic, moderately interacting particle system with a singular Coulomb-type kernel and a branching mechanism that encodes logistic damping, establishing a quantitative propagation of chaos-type result. A main contribution is a rigorous rate of convergence for the mollified empirical measure $u^N$ to the PDE solution $u$, with $\sup_{t\in[0,T]} \|u^N_t-u_t\|_{L^1\cap L^r} \lesssim N^{-\varrho}$ where $\varrho$ depends on the regularization and integrability parameters, and a corollary giving convergence in Kantorovich–Rubinstein distance for the un-mollified empirical measure. The results bridge macroscopic PDE dynamics and microscopic stochastic particle models in a setting with singular interactions and branching, providing insight into blow-up suppression through demographic effects and offering explicit convergence rates useful for numerical schemes and simulations.
Abstract
The Keller--Segel PDE is a model for chemotaxis known to exhibit possible finite-time blow-up. Following a seminal work by Tello and Winkler, a logistic damping term is added in this PDE and local well-posedness of mild solutions is proven. When the space dimension is $2$ or when the damping is strong enough, the solution is global in time. In the second part of this work, a microscopic description of this model is introduced in terms of a system of stochastic moderately interacting particles. This system features two main characteristics: the interaction between particles happens through a singular (Coulomb-type) kernel which is attractive; and the particles are subject to demographic events, birth and death due to local competition with other particles. The latter induces a branching structure of the particle system. Then the main result of this work is the convergence of the empirical measure of the particle system towards the Keller--Segel PDE with logistic damping, with a rate of order $N^{-\frac{1}{2(d+1)}}$.
