Electronic Journal of Probability
- Electron. J. Probab.
- Volume 22 (2017), paper no. 8, 37 pp.
Uniform in time interacting particle approximations for nonlinear equations of Patlak-Keller-Segel type
We study a system of interacting diffusions that models chemotaxis of biological cells or microorganisms (referred to as particles) in a chemical field that is dynamically modified through the collective contributions from the particles. Such systems of reinforced diffusions have been widely studied and their hydrodynamic limits that are nonlinear non-local partial differential equations are usually referred to as Patlak-Keller-Segel (PKS) equations.
Solutions of the classical PKS equation may blow up in finite time and much of the PDE literature has been focused on understanding this blow-up phenomenon. In this work we study a modified form of the PKS equation that is natural for applications and for which global existence and uniqueness of solutions are easily seen to hold. Our focus here is instead on the study of the long time behavior through certain interacting particle systems.
Under the so-called “quasi-stationary hypothesis” on the chemical field, the limit PDE reduces to a parabolic-elliptic system that is closely related to granular media equations whose time asymptotic properties have been extensively studied probabilistically through certain Lyapunov functions [17, 4, 9]. The modified PKS equation studied in the current work is a parabolic-parabolic system for which analogous Lyapunov function constructions are not available. A key challenge in the analysis is that the associated interacting particle system is not a Markov process as the interaction term depends on the whole history of the empirical measure.
We establish, under suitable conditions, uniform in time convergence of the empirical measure of particle states to the solution of the PDE. We also provide uniform in time exponential concentration bounds for rate of the above convergence under additional integrability conditions. Finally, we introduce an Euler discretization scheme for the simulation of the interacting particle system and give error bounds that show that the scheme converges uniformly in time and in the size of the particle system as the discretization parameter approaches zero.
Electron. J. Probab., Volume 22 (2017), paper no. 8, 37 pp.
Received: 24 July 2016
Accepted: 8 January 2017
First available in Project Euclid: 31 January 2017
Permanent link to this document
Digital Object Identifier
Mathematical Reviews number (MathSciNet)
Zentralblatt MATH identifier
Primary: 60K35: Interacting random processes; statistical mechanics type models; percolation theory [See also 82B43, 82C43] 60H30: Applications of stochastic analysis (to PDE, etc.) 60H35: Computational methods for stochastic equations [See also 65C30] 60K40: Other physical applications of random processes 60F05: Central limit and other weak theorems
weakly interacting particle systems uniform propagation of chaos McKean-Vlasov equations kinetic equations chemotaxis reinforced diffusions Patlak-Keller-Segel equations granular media equations uniform exponential concentration bounds long time behavior uniform in time Euler approximations
Budhiraja, Amarjit; Fan, Wai-Tong Louis. Uniform in time interacting particle approximations for nonlinear equations of Patlak-Keller-Segel type. Electron. J. Probab. 22 (2017), paper no. 8, 37 pp. doi:10.1214/17-EJP25. https://projecteuclid.org/euclid.ejp/1485831704