Deterministic versus Stochastic Chemical Kinetics
Deterministic versus Stochastic Chemical Kinetics
Traditional phenomenological kinetics uses first-order ordinary differential equations (ODEs) to describe the rate of change of the species in a chemical reaction. The solution to these equations provides expressions for the concentration of each species as a function of time. The use of ODEs to describe rate laws implies that concentrations are continuous quantities. As a consequence, this model is deterministic and provides exact solutions with no uncertainty involved. Because Avogadro's number is so enormous , even a nanomole of material involves over 100 trillion molecules, so in a typical lab-scale reaction it is reasonable to assume that concentrations are continuous variables. However, when concentrations or volumes become tiny (as in the inside of a human cell), this assumption can break down.
(~)
23
10
In truth, molecules are discrete entities: one-tenth of a water molecule does not exist in the physical world. As a consequence, chemical reactions involving small numbers of molecules are stochastic processes better described by probability methods. Stochastic simulation algorithms provide an alternative approach to chemical kinetics, in which trajectories consistent with these probability functions can be simulated via Monte Carlo methods.
The model allows you to compare the results predicted by a traditional phenomenological method to a stochastic simulation algorithm known as the Gillespie (or direct) method [1].