Nested stochastic simulation algorithm for chemical kinetic systems with disparate rates

摘要: An efficient simulation algorithm for chemical kinetic systems with disparate rates is proposed. This new quite general, and it amounts to a simple seamless modification of the classical stochastic (SSA), also known as Gillespie [J. Comput. Phys. 22, 403 (1976); J. Chem. 81, 2340 (1977)] algorithm. The basic idea use an outer SSA simulate slow processes computed from inner which simulates fast reactions. Averaging theorems Markov can be used identify variables in system well effective dynamics over time scale, even though itself does not rely on such information. nested easily generalized more than two separated scales. Convergence efficiency are discussed using established error estimates illustrated through examples.