The numerical stability of leaping methods for stochastic simulation of chemically reacting systems.

Tau-leaping methods have recently been proposed for the acceleration of discrete stochastic simulation of chemically reacting systems. This paper considers the numerical stability of these methods. The concept of stochastic absolute stability is defined, discussed, and applied to the following leaping methods: the explicit tau, implicit tau, and trapezoidal tau.

[1]  D. Gillespie.  The chemical Langevin equation , 2000 .

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

[3]  M. A. Akanbi,et al.  Numerical solution of initial value problems in differential - algebraic equations , 2005 .

[4]  A. Arkin,et al.  Stochastic mechanisms in gene expression. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

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

[6]  A. Arkin,et al.  Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. , 1998, Genetics.

[7]  Sheldon M. Ross.  Introduction to Probability Models. , 1995 .

[8]  Linda R. Petzold,et al.  Numerical solution of initial-value problems in differential-algebraic equations , 1996, Classics in applied mathematics.

[9]  Hong Li,et al.  Efficient formulation of the stochastic simulation algorithm for chemically reacting systems. , 2004, The Journal of chemical physics.

[10]  Sheldon M. Ross,et al.  Introduction to probability models , 1975 .

[11]  D. Gillespie.  Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .

[12]  Melvin K. Simmons,et al.  Hybrid simulation of cellular behavior , 2004, Bioinform..

[13]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[14]  D. Gillespie.  The Chemical Langevin and Fokker−Planck Equations for the Reversible Isomerization Reaction† , 2002 .

[15]  宁北芳,et al.  疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A , 2005 .

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

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

[18]  Muruhan Rathinam,et al.  Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method , 2003 .

[19]  A. Arkin,et al.  It's a noisy business! Genetic regulation at the nanomolar scale. , 1999, Trends in genetics : TIG.