Efficient exact and K -skip methods for stochastic simulation of coupled chemical reactions

Xiaodong Cai, Ji Wen

Research output: Contribution to journalArticle

7 Citations (Scopus)

Abstract

Gillespie's direct method (DM) [D. Gillespie, J. Chem. Phys. 81, 2340 (1977)] for exact stochastic simulation of chemical reaction systems has been widely adopted. It is easy to implement but requires large computation for relatively large systems. Recently, two more efficient methods, next reaction method (NRM) [M. A. Gibson and J. Bruck, J. Phys. Chem. A 105, 1876 (2000)] and optimized DM (ODM) [Y. Cao, J. Chem. Phys. 121, 4059 (2004)], have been developed to improve simulation speed. It has been demonstrated that the ODM is the state-of-the-art most efficient method for exact stochastic simulation of most practical reaction systems. In this paper, we first develop an exact stochastic simulation algorithm named ODMK that is more efficient than the ODM. We then develop an approximate method named K -skip method to further accelerate simulation. Using two chemical reaction systems, we demonstrate that our ODMK and K -skip method can save 20%-30% and 70%-80% simulation time, respectively, comparing to the ODM. We also show that our ODMK and K -skip method provide almost the same simulation accuracy as the ODM.

Original languageEnglish
Article number064108
JournalJournal of Chemical Physics
Volume131
Issue number6
DOIs
StatePublished - Aug 28 2009

Fingerprint

Chemical reactions
chemical reactions
simulation

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry
  • Medicine(all)

Cite this

Efficient exact and K -skip methods for stochastic simulation of coupled chemical reactions. / Cai, Xiaodong; Wen, Ji.

In: Journal of Chemical Physics, Vol. 131, No. 6, 064108, 28.08.2009.

Research output: Contribution to journalArticle

@article{3a2b733cbd424541b4d244c7ff0b865e,
title = "Efficient exact and K -skip methods for stochastic simulation of coupled chemical reactions",
abstract = "Gillespie's direct method (DM) [D. Gillespie, J. Chem. Phys. 81, 2340 (1977)] for exact stochastic simulation of chemical reaction systems has been widely adopted. It is easy to implement but requires large computation for relatively large systems. Recently, two more efficient methods, next reaction method (NRM) [M. A. Gibson and J. Bruck, J. Phys. Chem. A 105, 1876 (2000)] and optimized DM (ODM) [Y. Cao, J. Chem. Phys. 121, 4059 (2004)], have been developed to improve simulation speed. It has been demonstrated that the ODM is the state-of-the-art most efficient method for exact stochastic simulation of most practical reaction systems. In this paper, we first develop an exact stochastic simulation algorithm named ODMK that is more efficient than the ODM. We then develop an approximate method named K -skip method to further accelerate simulation. Using two chemical reaction systems, we demonstrate that our ODMK and K -skip method can save 20{\%}-30{\%} and 70{\%}-80{\%} simulation time, respectively, comparing to the ODM. We also show that our ODMK and K -skip method provide almost the same simulation accuracy as the ODM.",
author = "Xiaodong Cai and Ji Wen",
year = "2009",
month = "8",
day = "28",
doi = "10.1063/1.3204422",
language = "English",
volume = "131",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "6",

}

TY - JOUR

T1 - Efficient exact and K -skip methods for stochastic simulation of coupled chemical reactions

AU - Cai, Xiaodong

AU - Wen, Ji

PY - 2009/8/28

Y1 - 2009/8/28

N2 - Gillespie's direct method (DM) [D. Gillespie, J. Chem. Phys. 81, 2340 (1977)] for exact stochastic simulation of chemical reaction systems has been widely adopted. It is easy to implement but requires large computation for relatively large systems. Recently, two more efficient methods, next reaction method (NRM) [M. A. Gibson and J. Bruck, J. Phys. Chem. A 105, 1876 (2000)] and optimized DM (ODM) [Y. Cao, J. Chem. Phys. 121, 4059 (2004)], have been developed to improve simulation speed. It has been demonstrated that the ODM is the state-of-the-art most efficient method for exact stochastic simulation of most practical reaction systems. In this paper, we first develop an exact stochastic simulation algorithm named ODMK that is more efficient than the ODM. We then develop an approximate method named K -skip method to further accelerate simulation. Using two chemical reaction systems, we demonstrate that our ODMK and K -skip method can save 20%-30% and 70%-80% simulation time, respectively, comparing to the ODM. We also show that our ODMK and K -skip method provide almost the same simulation accuracy as the ODM.

AB - Gillespie's direct method (DM) [D. Gillespie, J. Chem. Phys. 81, 2340 (1977)] for exact stochastic simulation of chemical reaction systems has been widely adopted. It is easy to implement but requires large computation for relatively large systems. Recently, two more efficient methods, next reaction method (NRM) [M. A. Gibson and J. Bruck, J. Phys. Chem. A 105, 1876 (2000)] and optimized DM (ODM) [Y. Cao, J. Chem. Phys. 121, 4059 (2004)], have been developed to improve simulation speed. It has been demonstrated that the ODM is the state-of-the-art most efficient method for exact stochastic simulation of most practical reaction systems. In this paper, we first develop an exact stochastic simulation algorithm named ODMK that is more efficient than the ODM. We then develop an approximate method named K -skip method to further accelerate simulation. Using two chemical reaction systems, we demonstrate that our ODMK and K -skip method can save 20%-30% and 70%-80% simulation time, respectively, comparing to the ODM. We also show that our ODMK and K -skip method provide almost the same simulation accuracy as the ODM.

UR - http://www.scopus.com/inward/record.url?scp=69149100616&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=69149100616&partnerID=8YFLogxK

U2 - 10.1063/1.3204422

DO - 10.1063/1.3204422

M3 - Article

VL - 131

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 6

M1 - 064108

ER -