Sample animation
The animation displays the time-dependent probability density of a
very simple interacting chemical system consisting of two species
that comes to rest at steady-state. Although the density itself is
in equilibrium, any realization of the system would be fluctuating
according to this distribution. Note the negative correlation
which comes from the interaction between the two species; if the
number of molecules of one of the species increases, then the
other species will start to disappear from the system due to the
interaction. Note also the different scales present already in
this simple example: the transport along the line x = y is very
fast while the slow diffusion along x y = constant is much slower.
(Plot info: each axis tick is 10 molecules and the contour levels
have been chosen by Matlab)
References
S. Engblom: Spectral Approximation of Solutions to the Chemical
Master Equation, in J. Comput. Appl. Math.
229(1):208--221,
2009: (doi)
S. Engblom: Galerkin Spectral Method applied to the
Chemical Master Equation, in Commun. Comput. Phys.
5(5):871--896,
2009: (abstract), (pdf).
Stefan Engblom
Last modified: Thu May 10 11:04:40 MEST 2007