Abstract:Stochastic reaction-diffusion models are employed to represent many complex physical, biological, societal, and ecological systems. The macroscopic reaction rates describing the large-scale kinetics in such systems are effective, scale-dependent parameters that need to be either measured experimentally or computed using a microscopic model. In a Monte Carlo simulation of stochastic reaction-diffusion systems, microscopic probabilities for specific events to happen serve as the input control parameters. Finding the functional dependence of emergent macroscopic rates on the microscopic probabilities is a difficult problem, and there is no systematic analytical way to achieve this goal. Therefore, we introduce a straightforward numerical method of using Monte Carlo simulations to evaluate the macroscopic reaction rates by directly obtaining the count statistics of how many events occur per simulation time step. Our technique is first tested on well-understood fundamental examples, namely restricted birth processes, diffusion-limited two-particle coagulation, and two-species pair annihilation kinetics. Next we utilize the thus gained experience to investigate how the microscopic algorithmic probabilities become coarse-grained into effective macroscopic rates in more complex model systems such as the Lotka--Volterra model for predator-prey competition and coexistence, as well as the rock-paper-scissors or cyclic Lotka--Volterra model as well as its May--Leonard variant that capture population dynamics with cyclic dominance motifs. Thereby we achieve a deeper understanding of coarse-graining in spatially extended stochastic reaction systems and the nontrivial relationships between the associated microscopic and macroscopic model parameters. The proposed technique should generally provide a useful means to better fit Monte Carlo simulation results to experimental or observational data.

Monte Carlo simulation with fixed steplength for diffusion processes in nonhomogeneous media

Monte-Carlo Simulation of a Stochastic Differential Equation

Grand canonical diffusion-influenced reactions: a stochastic theory with applications to multiscale reaction-diffusion simulations

Reaction-controlled diffusion: Monte Carlo simulations

On the Mathematical foundations of Diffusion Monte Carlo

Monte Carlo simulation of the reaction–diffusion process of a complex chemical scheme on a fractal lattice

Review of Monte Carlo methods for particle transport in continuously-varying media

A Monte Carlo simulation for kinetic chemotaxis models: an application to the traveling population wave

Stochastic method for accommodation of equilibrating basins in kinetic Monte Carlo simulations

Computing macroscopic reaction rates in reaction-diffusion systems using Monte Carlo simulations

Asymptotic-Preserving Monte Carlo Methods for Transport Equations in the Diffusive Limit

Monte-Carlo Simulation of Crystal Growth with Diffusion Employing an Anisotropic Model

Stochastic Differential Equations for Practical Simulation of Gene Circuits

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

Discrete stochastic modeling of calcium channel dynamics

Finite Volume Model to Study Calcium Diffusion in Neuron Involving JRYR, JSERCA and JLEAK

Monte Carlo simulation of electrostatic interactions in inhomogeneous dielectric media: Correct sampling for the local lattice simulation algorithm

Equation-free, multiscale computation for unsteady random diffusion

First-Passage Kinetic Monte Carlo method

Effects of different discretisations of the Laplacian upon stochastic simulations of reaction-diffusion systems on both static and growing domains

Monte Carlo and kinetic Monte Carlo methods