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, Journal Name: Journal of Chemical Physics Journal Issue: 23 Vol. 143; ISSN JCPSA6; ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
HRSSA – Efficient hybrid stochastic simulation for spatially homogeneous biochemical reaction networks
An adaptive tau-leaping method for stochastic simulations of reaction-diffusion systems