Hybrid approaches for multiple-species stochastic reaction-diffusion models

Abstract

Reaction-diffusion models are used to describe systems in fields as diverse as physics, chemistry, ecology and biology. The fundamental quantities in such models are individual entities such as atoms and molecules, bacteria, cells or animals, which move and/or react in a stochastic manner. If the number of entities is large, accounting for each individual is inefficient, and often partial differential equation (PDE) models are used in which the stochastic behaviour of individuals is replaced by a description of the averaged, or mean behaviour of the system. In some situations the number of individuals is large in certain regions and small in others. In such cases, a stochastic model may be inefficient in one region, and a PDE model inaccurate in another. To overcome this problem, we develop a scheme which couples a stochastic reaction-diffusion system in one part of the domain with its mean field analogue, i.e. a discretised PDE model, in the other part of the domain. The interface in between the two domains occupies exactly one lattice site and is chosen such that the mean field description is still accurate there. This way errors due to the flux between the domains are small. Our scheme can account for multiple dynamic interfaces separating multiple stochastic and deterministic domains, and the coupling between the domains conserves the total number of particles. The method preserves stochastic features such as extinction not observable in the mean field description, and is significantly faster to simulate on a computer than the pure stochastic model.

arXiv comment38 pages, 8 figures

References in corpus (11)

• Simulation of stochastic reaction-diffusion processes on unstructured meshes
• A Convergent Reaction-Diffusion Master Equation
• Spatially hybrid computations for streamer discharges: II. Fully 3D simulations
• Spatially hybrid computations for streamer discharges with generic features of pulled fronts: I. Planar fronts
• Mesoscopic and continuum modelling of angiogenesis
• Adaptive Hybrid Simulations for Multiscale Stochastic Reaction Networks
• Hybrid method for simulating front propagation in reaction-diffusion systems
• Reaction rates for mesoscopic reaction-diffusion kinetics
• Adaptive two-regime method: application to front propagation
• Constant-complexity Stochastic Simulation Algorithm with Optimal Binning
• Convergence of methods for coupling of microscopic and mesoscopic reaction-diffusion simulations

Cited by in corpus (11)

• Variational system identification of the partial differential equations governing pattern-forming physics: Inference under varying fidelity and noise
• Variational system identification of the partial differential equations governing microstructure evolution in materials: Inference over sparse and spatially unrelated data
• Perspectives on the mathematics of biological patterning and morphogenesis
• A Perspective on Regression and Bayesian Approaches for System Identification of Pattern Formation Dynamics
• Numerical algorithms based on Galerkin methods for the modeling of reactive interfaces in photoelectrochemical (PEC) solar cells
• Stochastic multi-scale models of competition within heterogeneous cellular populations: simulation methods and mean-field analysis
• Optimisation of Simulations of Stochastic Processes by Removal of Opposing Reactions
• Patch formation driven by stochastic effects of interaction between viruses and defective interfering particles
• An equivalence framework for an age-structured multi-stage representation of the cell cycle
• Parameter Estimation in an SPDE Model for Cell Repolarisation
• A variational-level-set based partitioning method for block-structured meshes