A moment closure method for stochastic reaction networks

Chang Hyeong Lee, Kyeong Hun Kim, Pilwon Kim

Research output: Contribution to journalArticle

72 Citations (Scopus)

Abstract

In this paper we present a moment closure method for stochastically modeled chemical or biochemical reaction networks. We derive a system of differential equations which describes the dynamics of means and all central moments from a chemical master equation. Truncating the system for the central moments at a certain moment term and using Taylor approximation, we obtain explicit representations of means and covariances and even higher central moments in recursive forms. This enables us to deal with the moments in successive differential equations and use conventional numerical methods for their approximations. Furthermore, we estimate the errors in the means and central moments generated by the approximation method. We also find the moments at equilibrium by solving truncated algebraic equations. We show in examples that numerical solutions based on the moment closure method are accurate and efficient by comparing the results to those of stochastic simulation algorithms.

Original languageEnglish
Article number134107
JournalJournal of Chemical Physics
Volume130
Issue number13
DOIs
Publication statusPublished - 2009 Apr 20

Fingerprint

closures
Differential equations
moments
Numerical methods
differential equations
approximation
chemical reactions
estimates

ASJC Scopus subject areas

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

Cite this

A moment closure method for stochastic reaction networks. / Lee, Chang Hyeong; Kim, Kyeong Hun; Kim, Pilwon.

In: Journal of Chemical Physics, Vol. 130, No. 13, 134107, 20.04.2009.

Research output: Contribution to journalArticle

@article{c13e881884704ff4a47a46b0eb4bfac7,
title = "A moment closure method for stochastic reaction networks",
abstract = "In this paper we present a moment closure method for stochastically modeled chemical or biochemical reaction networks. We derive a system of differential equations which describes the dynamics of means and all central moments from a chemical master equation. Truncating the system for the central moments at a certain moment term and using Taylor approximation, we obtain explicit representations of means and covariances and even higher central moments in recursive forms. This enables us to deal with the moments in successive differential equations and use conventional numerical methods for their approximations. Furthermore, we estimate the errors in the means and central moments generated by the approximation method. We also find the moments at equilibrium by solving truncated algebraic equations. We show in examples that numerical solutions based on the moment closure method are accurate and efficient by comparing the results to those of stochastic simulation algorithms.",
author = "Lee, {Chang Hyeong} and Kim, {Kyeong Hun} and Pilwon Kim",
year = "2009",
month = "4",
day = "20",
doi = "10.1063/1.3103264",
language = "English",
volume = "130",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "13",

}

TY - JOUR

T1 - A moment closure method for stochastic reaction networks

AU - Lee, Chang Hyeong

AU - Kim, Kyeong Hun

AU - Kim, Pilwon

PY - 2009/4/20

Y1 - 2009/4/20

N2 - In this paper we present a moment closure method for stochastically modeled chemical or biochemical reaction networks. We derive a system of differential equations which describes the dynamics of means and all central moments from a chemical master equation. Truncating the system for the central moments at a certain moment term and using Taylor approximation, we obtain explicit representations of means and covariances and even higher central moments in recursive forms. This enables us to deal with the moments in successive differential equations and use conventional numerical methods for their approximations. Furthermore, we estimate the errors in the means and central moments generated by the approximation method. We also find the moments at equilibrium by solving truncated algebraic equations. We show in examples that numerical solutions based on the moment closure method are accurate and efficient by comparing the results to those of stochastic simulation algorithms.

AB - In this paper we present a moment closure method for stochastically modeled chemical or biochemical reaction networks. We derive a system of differential equations which describes the dynamics of means and all central moments from a chemical master equation. Truncating the system for the central moments at a certain moment term and using Taylor approximation, we obtain explicit representations of means and covariances and even higher central moments in recursive forms. This enables us to deal with the moments in successive differential equations and use conventional numerical methods for their approximations. Furthermore, we estimate the errors in the means and central moments generated by the approximation method. We also find the moments at equilibrium by solving truncated algebraic equations. We show in examples that numerical solutions based on the moment closure method are accurate and efficient by comparing the results to those of stochastic simulation algorithms.

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

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

U2 - 10.1063/1.3103264

DO - 10.1063/1.3103264

M3 - Article

C2 - 19355717

AN - SCOPUS:64549159754

VL - 130

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 13

M1 - 134107

ER -