Reversible-equivalent-monomolecular tau: A leaping method for 'small number and stiff' stochastic chemical systems
- Department of Mathematics and Statistics, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250 (United States)
- Department of Biochemistry and Biophysics, California Institute for Quantitative Biomedical Research, University of California, San Francisco, 1700, 4th Street, San Francisco, CA 94143-2542 (United States)
Leaping methods provide for efficient and approximate time stepping of chemical reaction systems modeled by continuous time discrete state stochastic dynamics. We investigate the application of leaping methods for 'small number and stiff' systems, i.e. systems whose dynamics involve different time scales and have some molecular species present in very small numbers, specifically in the range 0 to 10. We propose a new explicit leaping scheme, reversible-equivalent-monomolecular tau (REMM-{tau}), which shows considerable promise in the simulation of such systems. The REMM-{tau} scheme is based on the fact that the exact solution of the two prototypical monomolecular reversible reactions S {sub 1} {r_reversible} S {sub 2} and S {r_reversible} 0 as a function of time takes a simple form involving binomial and/or Poisson random variables. The REMM-{tau} method involves approximating bimolecular reversible reactions by suitable monomolecular reversible reactions as well as considering each reversible pair of reactions in the system to be operating in isolation during the time step {tau}. We illustrate the use of the REMM-{tau} method through a number of biologically motivated examples and compare its performance to those of the implicit-{tau} and trapezoidal implicit-{tau} algorithms. In most cases considered, REMM-{tau} appears to perform better than these two methods while having the important advantage of being computationally faster due to the explicit nature of the method. Furthermore when stepsize {tau} is increased the REMM-{tau} exhibits a more robust performance than the implicit-{tau} or the trapezoidal implicit-{tau} for small number stiff problems.
- OSTI ID:
- 20991590
- Journal Information:
- Journal of Computational Physics, Vol. 224, Issue 2; Other Information: DOI: 10.1016/j.jcp.2006.10.034; PII: S0021-9991(06)00554-7; Copyright (c) 2006 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
An investigation into the accuracy, stability and parallel performance of a highly stable explicit technique for stiff reaction-transport PDEs
Highly stable explicit technique for stiff reaction-transport PDEs