An equation-free probabilistic steady state approximation: Dynamic application to the stochastic simulation of chemically reacting systems

Howard M. Salis, Yiannis Kaznessis

Research output: Chapter in Book/Report/Conference proceedingConference contribution

1 Scopus citations

Abstract

'Small' systems of chemical or biochemical reactions typically require a mesoscopic mathematical representation, such as a jump Markov process. Kinetic Monte Carlo, also known as stochastic simulation, is a numerical method for generating trajectories of a jump Markov process. However, the computational cost of these simulations scale with the number of occurrences of reaction or diffusion events in the system. We present an equation-free probabilistic steady state approximation that speeds up the simulation of jump Markov processes which contain a separation in timescales. The method partitions the system into subsets of fast/discrete and slow/discrete reactions; detects the convergence of the unknown distribution to a stable quasi steady state; samples from the underlying distribution; and uses those samples to compute the waiting time of the next slow/discrete reaction and the state of the system at that time. The method is equation-free because the steady state distribution is not computed via the Master equation or any differential equation. Instead, by sampling directly from the unknown distribution, we only assume the quasi steady state distribution is ergodic. We present multiple examples and demonstrate the accuracy and efficiency of the proposed method. The presented ideas apply to both homogeneous and heterogeneous systems, but, for simplicity, we will focus on homogeneous ones.

Original languageEnglish (US)
Title of host publication05AIChE: 2005 AIChE Annual Meeting and Fall Showcase, Conference Proceedings
Pages11469
Number of pages1
StatePublished - 2005
Event05AIChE: 2005 AIChE Annual Meeting and Fall Showcase - Cincinnati, OH, United States
Duration: Oct 30 2005Nov 4 2005

Other

Other05AIChE: 2005 AIChE Annual Meeting and Fall Showcase
CountryUnited States
CityCincinnati, OH
Period10/30/0511/4/05

All Science Journal Classification (ASJC) codes

  • Engineering(all)

Fingerprint Dive into the research topics of 'An equation-free probabilistic steady state approximation: Dynamic application to the stochastic simulation of chemically reacting systems'. Together they form a unique fingerprint.

Cite this