An adaptive tau-leaping method for stochastic simulations of reaction-diffusion systems
- Department of Mathematics, Ryerson University, Toronto, ON, M5B 2K3 (Canada)
Stochastic modelling is critical for studying many biochemical processes in a cell, in particular when some reacting species have low population numbers. For many such cellular processes the spatial distribution of the molecular species plays a key role. The evolution of spatially heterogeneous biochemical systems with some species in low amounts is accurately described by the mesoscopic model of the Reaction-Diffusion Master Equation. The Inhomogeneous Stochastic Simulation Algorithm provides an exact strategy to numerically solve this model, but it is computationally very expensive on realistic applications. We propose a novel adaptive time-stepping scheme for the tau-leaping method for approximating the solution of the Reaction-Diffusion Master Equation. This technique combines effective strategies for variable time-stepping with path preservation to reduce the computational cost, while maintaining the desired accuracy. The numerical tests on various examples arising in applications show the improved efficiency achieved by the new adaptive method.
- OSTI ID:
- 22611622
- Journal Information:
- AIP Advances, Vol. 6, Issue 3; Other Information: (c) 2016 Author(s); Country of input: International Atomic Energy Agency (IAEA); ISSN 2158-3226
- Country of Publication:
- United States
- Language:
- English
Similar Records
Stochastic simulation of reaction-diffusion systems: A fluctuating-hydrodynamics approach
An adaptive multi-level simulation algorithm for stochastic biological systems