Premium
Stochastic models for chemically reacting systems using polynomial stochastic hybrid systems
Author(s) -
Hespanha João Pedro,
Singh Abhyudai
Publication year - 2005
Publication title -
international journal of robust and nonlinear control
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.361
H-Index - 106
eISSN - 1099-1239
pISSN - 1049-8923
DOI - 10.1002/rnc.1017
Subject(s) - bistability , ode , polynomial , nonlinear system , master equation , stochastic differential equation , mathematics , ordinary differential equation , population , statistical physics , differential equation , mathematical analysis , physics , demography , quantum mechanics , sociology , quantum
A stochastic model for chemical reactions is presented, which represents the population of various species involved in a chemical reaction as the continuous state of a polynomial stochastic hybrid system (pSHS). pSHSs correspond to stochastic hybrid systems with polynomial continuous vector fields, reset maps, and transition intensities. We show that for pSHSs, the dynamics of the statistical moments of its continuous states, evolves according to infinite‐dimensional linear ordinary differential equations (ODEs), which can be approximated by finite‐dimensional nonlinear ODEs with arbitrary precision. Based on this result, a procedure to build this types of approximation is provided. This procedure is used to construct approximate stochastic models for a variety of chemical reactions that have appeared in literature. These reactions include a simple bimolecular reaction, for which one can solve the Master equation; a decaying–dimerizing reaction set which exhibits two distinct time scales; a reaction for which the chemical rate equations have a continuum of equilibrium points; and the bistable Schögl reaction. The accuracy of the approximate models is investigated by comparing with Monte Carlo simulations or the solution to the Master equation, when available. Copyright © 2005 John Wiley & Sons, Ltd.