Idea Behind The Algorithm
Traditional continuous and deterministic biochemical rate equations do not accurately predict cellular reactions since they rely on bulk reactions that require the interactions of millions of molecules. They are typically modeled as a set of coupled ordinary differential equations. In contrast, the Gillespie algorithm allows a discrete and stochastic simulation of a system with few reactants because every reaction is explicitly simulated. When simulated, a Gillespie realization represents a random walk that exactly represents the distribution of the master equation.
The physical basis of the algorithm is the collision of molecules within a reaction vessel. It is assumed that collisions are frequent, but collisions with the proper orientation and energy are infrequent. Therefore, all reactions within the Gillespie framework must involve at most two molecules. Reactions involving three molecules are assumed to be extremely rare and are modeled as a sequence of binary reactions. It is also assumed that the reaction environment is well mixed.
Read more about this topic: Gillespie Algorithm
Famous quotes containing the word idea:
“A cow does not know how much milk it has until the milkman starts working on it. Then it looks round in surprise and sees the pail full to the brim. In the same way a writer has no idea how much he has to say till his pen draws it out of him. Thoughts will then appear on the paper that he is amazed to find that he possessed. “How brilliant!” he says to himself. “I had no idea I was so intelligent.” But the reader may not be so im pressed.”
—Gerald Branan (1894–1987)