A modified next reaction method for simulating chemical systems with time dependent propensities and delays
J. Chem. Phys. 127, 214107 (2007); doi:10.1063/1.2799998
Published 6 December 2007 | See: Publisher's Note
You are not logged in to this journal. Log in
Chemical reaction systems with a low to moderate number of molecules are typically modeled as discrete jump Markov processes. These systems are oftentimes simulated with methods that produce statistically exact sample paths such as the Gillespie algorithm or the next reaction method. In this paper we make explicit use of the fact that the initiation times of the reactions can be represented as the firing times of independent, unit rate Poisson processes with internal times given by integrated propensity functions. Using this representation we derive a modified next reaction method and, in a way that achieves efficiency over existing approaches for exact simulation, extend it to systems with time dependent propensities as well as to systems with delays.
©2007 American Institute of Physics
| History: | Received 20 June 2007; accepted 26 September 2007; published 6 December 2007; publisher error corrected 28 January 2008 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/127/214107/1 |
ERRATUM
- Publisher's Note: “A modified next reaction method for simulating chemical systems with time dependent propensities and delays” [J. Chem. Phys. 127, 214107 (2007)]
David F. Anderson
J. Chem. Phys. 128, 109903 (2008)
KEYWORDS and PACS
RELATED DATABASES
PUBLICATION DATA
0021-9606 (print)
1089-7690 (online)
REFERENCES (20)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- A. Arkin, J. Ross, and H. H. McAdams,
Genetics 149, 1633 (1998) . - H. H. McAdams and A. Arkin,
Proc. Natl. Acad. Sci. U.S.A. 94, 814 (1997) . - E. M. Ozbudak, M. Thattai, I. Kurtser, A. D. Grossman, and A. van Oudenaarden,
Nat. Genet. 31, 69 (2002) . - H. E. Samad, M. Khammash, L. Petzold, and D. Gillespie,
Int. J. Robust Nonlinear Control 15, 691 (2005) . - D. T. Gillespie,
J. Comput. Phys. 22, 403 (1976) . - D. T. Gillespie,
J. Phys. Chem. 81, 2340 (1977) . - M. Gibson and J. Bruck,
J. Phys. Chem. A 105, 1876 (2000) . - T. G. Kurtz,
Ann. Probab. 8, 682 (1980) . - S. N. Ethier and T. G. Kurtz, Markov Processes: Characterization and Convergence (Wiley, New York, 1986).
- K. Ball, T. G. Kurtz, L. Popovic, and G. Rempala,
Ann. Appl. Probab. 16, 1925 (2006) . - D. Bratsun, D. Volfson, L. S. Tsimring, and J. Hasty,
Proc. Natl. Acad. Sci. U.S.A. 102, 14593 (2005) . - M. Barrio, K. Burrage, A. Leier, and T. Tian, PLOS Comput. Biol. 2, 1017 (2006).
- X. Cai, J. Chem. Phys. 126, 124108 (2007).
- D. Y. Burman,
Adv. Appl. Probab. 13, 846 (1981) . - R. Schassberger,
Adv. Appl. Probab. 10, 836 (1978) . - P. W. Glynn,
Proc. IEEE 77, 14 (1989) . - P. J. Haas, Stochastic Petri Nets: Modelling Stability, Simulation, 1st ed. (Springer, New York, 2002).
- D. F. Anderson, “Incorporating postleap checks in tau-leaping,” J. Chem. Phys. (to be published).
- T. Lu, D. Volfson, L. Tsimring, and J. Hasty, IEE Syst. Biol. 1, 1 (2004).
- G. A. Rempala, K. S. Ramos, and T. Kalbfleisch,
J. Theor. Biol. 242, 101 (2006) .








