Approximate method for stochastic chemical kinetics with two-time scales by chemical Langevin equations
- School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan, Hubei 430074 (China)
- School of Mathematical Sciences, Monash University, Melbourne, Vic 3800 (Australia)
- Department of Mathematics, Wayne State University, Detroit, Michigan 48202 (United States)
The frequently used reduction technique is based on the chemical master equation for stochastic chemical kinetics with two-time scales, which yields the modified stochastic simulation algorithm (SSA). For the chemical reaction processes involving a large number of molecular species and reactions, the collection of slow reactions may still include a large number of molecular species and reactions. Consequently, the SSA is still computationally expensive. Because the chemical Langevin equations (CLEs) can effectively work for a large number of molecular species and reactions, this paper develops a reduction method based on the CLE by the stochastic averaging principle developed in the work of Khasminskii and Yin [SIAM J. Appl. Math. 56, 1766–1793 (1996); ibid. 56, 1794–1819 (1996)] to average out the fast-reacting variables. This reduction method leads to a limit averaging system, which is an approximation of the slow reactions. Because in the stochastic chemical kinetics, the CLE is seen as the approximation of the SSA, the limit averaging system can be treated as the approximation of the slow reactions. As an application, we examine the reduction of computation complexity for the gene regulatory networks with two-time scales driven by intrinsic noise. For linear and nonlinear protein production functions, the simulations show that the sample average (expectation) of the limit averaging system is close to that of the slow-reaction process based on the SSA. It demonstrates that the limit averaging system is an efficient approximation of the slow-reaction process in the sense of the weak convergence.
- OSTI ID:
- 22657957
- Journal Information:
- Journal of Chemical Physics, Journal Name: Journal of Chemical Physics Journal Issue: 17 Vol. 144; ISSN JCPSA6; ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Explicit time integration of the stiff chemical Langevin equations using computational singular perturbation
Hybrid framework for the simulation of stochastic chemical kinetics