Hybrid stochastic simulation of reaction-diffusion systems with slow and fast dynamics
- Department of Mathematics, Ryerson University, Toronto, Ontario M5B 2K3 (Canada)
In this paper, we present a novel hybrid method to simulate discrete stochastic reaction-diffusion models arising in biochemical signaling pathways. We study moderately stiff systems, for which we can partition each reaction or diffusion channel into either a slow or fast subset, based on its propensity. Numerical approaches missing this distinction are often limited with respect to computational run time or approximation quality. We design an approximate scheme that remedies these pitfalls by using a new blending strategy of the well-established inhomogeneous stochastic simulation algorithm and the tau-leaping simulation method. The advantages of our hybrid simulation algorithm are demonstrated on three benchmarking systems, with special focus on approximation accuracy and efficiency.
- OSTI ID:
- 22493341
- Journal Information:
- Journal of Chemical Physics, Vol. 143, Issue 23; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
HRSSA – Efficient hybrid stochastic simulation for spatially homogeneous biochemical reaction networks
Slow-scale split-step tau-leap method for stiff stochastic chemical systems