Abstract: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.

A particle method for non-local advection–selection–mutation equations

A particle method for non-local advection-selection-mutation equations

A Multilevel Monte Carlo Asymptotic-Preserving Particle Method for Kinetic Equations in the Diffusion Limit

Particle method for the numerical simulation of the path-dependent McKean-Vlasov equation

Regularity and trend to equilibrium for a non-local advection-diffusion model of active particles

Compound Poisson particle approximation for McKean-Vlasov SDEs

A Particle Method for the Semiclassical Limit of the Schrödinger Equation and the Vlasov--Poisson Equations

Nonlocal particle approximation for linear and fast diffusion equations

The approximation of the quadratic porous medium equation via nonlocal interacting particles subject to repulsive Morse potential

On a practical implementation of particle methods

Particle Method for the McKean-Vlasov equation with common noise

Uniform in time interacting particle approximations for nonlinear equations of Patlak-Keller-Segel type

A meshless stochastic method for Poisson-Nernst-Planck equations

A Particle Method for the KdV Equation

Convergence Of A Particle Method And Global Weak Solutions Of A Family Of Evolutionary Pdes

A statistical approach for simulating the density solution of a McKean-Vlasov equation

Particle Methods for Dispersive Equations

Particle method and quantization-based schemes for the simulation of the McKean-Vlasov equation

A High Order Stochastic Asymptotic Preserving Scheme for Chemotaxis Kinetic Models with Random Inputs

A Practical Guide to Deterministic Particle Methods