Gillespie's Stochastic Simulation Algorithm for Chemical Reactions
Gillespie's Stochastic Simulation Algorithm for Chemical Reactions
The Gillespie stochastic simulation algorithm is a Monte Carlo algorithm that simulates the time evolution of a chemical system according to its chemical master equation. In this Demonstration, the reaction is simulated with variable forward and backward rate constants and initial particle numbers. The Demonstration shows the key steps of the Gillespie method. At the top left, the time until the next reaction is drawn from an exponential distribution. At the bottom left, the next reaction step (forward or backward) is randomly chosen based on the probabilities at each step of the reaction. The plot on the right shows the time course of each species , , . Use the "reaction event" slider to scroll in time through the reactions that occur.
A+B⇄C
k
f
k
b
A
B
C