Nested stochastic simulation algorithms for chemical kinetic systems with multiple time scales

We present an efficient numerical algorithm for simulating chemical kinetic systems with multiple time scales. This algorithm is an improvement of the traditional stochastic simulation algorithm (SSA), also known as Gillespie's algorithm. It is in the form of a nested SSA and uses an outer SSA to simulate the slow reactions with rates computed from realizations of inner SSAs that simulate the fast reactions. The algorithm itself is quite general and seamless, and it amounts to a small modification of the original SSA. Our analysis of such multi-scale chemical kinetic systems allows us to identify the slow variables in the system, derive effective dynamics on the slow time scale, and provide error estimates for the nested SSA. Efficiency of the nested SSA is discussed using these error estimates, and illustrated through several numerical examples.

[1]  Eric Vanden Eijnden.  Numerical techniques for multi-scale dynamical systems with stochastic effects , 2003 .

[2]  L. You,et al.  Stochastic vs. deterministic modeling of intracellular viral kinetics. , 2002, Journal of theoretical biology.

[3]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[4]  R. Z. Khasminskii,et al.  Constructing asymptotic series for probability distributions of Markov chains with weak and strong interactions , 1997 .

[5]  W. Fontana,et al.  Small Numbers of Big Molecules , 2002, Science.

[6]  Eric Vanden-Eijnden,et al.  Nested stochastic simulation algorithm for chemical kinetic systems with disparate rates. , 2005, The Journal of chemical physics.

[7]  J. Rawlings,et al.  Approximate simulation of coupled fast and slow reactions for stochastic chemical kinetics , 2002 .

[8]  A. Bensoussan,et al.  Asymptotic analysis for periodic structures , 1979 .

[9]  R. Has’minskiĭ.  On Stochastic Processes Defined by Differential Equations with a Small Parameter , 1966 .

[10]  A. B. Bortz,et al.  A new algorithm for Monte Carlo simulation of Ising spin systems , 1975 .

[11]  S. Meyn,et al.  Stability of Markovian processes I: criteria for discrete-time Chains , 1992, Advances in Applied Probability.

[12]  R. Khas'minskii.  A Limit Theorem for the Solutions of Differential Equations with Random Right-Hand Sides , 1966 .

[13]  C. Rao,et al.  Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm , 2003 .

[14]  Thomas G. Kurtz,et al.  A limit theorem for perturbed operator semigroups with applications to random evolutions , 1973 .

[15]  Björn Engquist,et al.  Heterogeneous multiscale methods for stiff ordinary differential equations , 2005, Math. Comput..

[16]  E. Vanden-Eijnden,et al.  The Heterogeneous Multiscale Method: A Review , 2007 .

[17]  Yang Cao,et al.  Multiscale stochastic simulation algorithm with stochastic partial equilibrium assumption for chemically reacting systems , 2005 .

[18]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[19]  Eric Vanden-Eijnden,et al.  NUMERICAL TECHNIQUES FOR MULTI-SCALE DYNAMICAL SYSTEMS WITH STOCHASTIC EFFECTS ⁄ , 2003 .

[20]  C. Petri.  Kommunikation mit Automaten , 1962 .

[21]  Masaru Tomita,et al.  A multi-algorithm, multi-timescale method for cell simulation , 2004, Bioinform..

[22]  D G Vlachos,et al.  Overcoming stiffness in stochastic simulation stemming from partial equilibrium: a multiscale Monte Carlo algorithm. , 2005, The Journal of chemical physics.

[23]  Linda R Petzold,et al.  The slow-scale stochastic simulation algorithm. , 2005, The Journal of chemical physics.

[24]  I. Amdur,et al.  Kinetic Theory of Gases , 1959 .

[25]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[26]  E. Vanden-Eijnden,et al.  Analysis of multiscale methods for stochastic differential equations , 2005 .