How reaction-diffusion PDEs approximate the large-population limit of stochastic particle models
Files
Published version
Date
2021-01
Authors
Isaacson, Samuel A.
Ma, Jingwei
Spiliopoulos, Konstantinos
Version
OA Version
Published version
Citation
S.A. Isaacson, J. Ma, K. Spiliopoulos. 2021. "How Reaction-Diffusion PDEs Approximate the Large-Population Limit of Stochastic Particle Models" SIAM Journal on Applied Mathematics, Volume 81, Issue 6, pp.2622-2657. https://doi.org/10.1137/20m1365429
Abstract
Reaction-diffusion PDEs and particle-based stochastic reaction-diffusion (PBSRD) models are commonly used approaches for modeling the spatial dynamics of chemical and biological systems. Standard reaction-diffusion PDE models ignore the underlying stochasticity of spatial transport and reactions and are often described as appropriate in regimes where there are large numbers of particles in a system. Recent studies have proven the rigorous large-population limit of PBSRD models, showing the resulting mean-field models (MFMs) correspond to nonlocal systems of partial-integro differential equations. In this work we explore the rigorous relationship between standard reaction-diffusion PDE models and the derived MFM. We prove that the former can be interpreted as an asymptotic approximation to the later in the limit that bimolecular reaction kernels are short-range and averaging. As the reactive interaction length scale approaches zero, we prove the MFMs converge at second order to standard reaction-diffusion PDE models. In proving this result we also establish local well-posedness of the MFM model in time for general systems and global well-posedness for specific reaction systems and kernels. Finally, we illustrate the agreement and disagreement between the MFM, SM, and underlying particle model for several numerical examples.
Description
License
Copyright © by SIAM. Unauthorized reproduction of this article is prohibited.