Open Access
2015 Perturbation and Truncation of Probability Generating Function Methods for Stiff Chemical Reactions
Soyeong Jeong, Pilwon Kim, Chang Hyeong Lee
J. Appl. Math. 2015: 1-8 (2015). DOI: 10.1155/2015/864238

Abstract

One can reformulate chemical master equations of the stochastic reaction network into a partial differential equation (PDE) of a probability generating function (PGF). In this paper, we present two improvements in such PGF-PDE approach, based on perturbation and double-truncation, respectively. The stiff system that involves fast and slow reactions together often requires high computational cost. By applying the perturbation method to PGF-PDEs, we expand the equation in terms of a small reaction rate which is often responsible for such stiffness of the system. Also by doubly truncating, we dump relatively small terms and reduce the computational load significantly at each time step. The terms corresponding to rare events are sieved out through truncations of Taylor expansion. It is shown through numerical examples of enzyme kinetics, transition model, and Brusselator model that the suggested method is accurate and efficient for approximation of the state probabilities.

Citation

Download Citation

Soyeong Jeong. Pilwon Kim. Chang Hyeong Lee. "Perturbation and Truncation of Probability Generating Function Methods for Stiff Chemical Reactions." J. Appl. Math. 2015 1 - 8, 2015. https://doi.org/10.1155/2015/864238

Information

Published: 2015
First available in Project Euclid: 11 June 2015

zbMATH: 07132086
MathSciNet: MR3354363
Digital Object Identifier: 10.1155/2015/864238

Rights: Copyright © 2015 Hindawi

Vol.2015 • 2015
Back to Top