An exact accelerated stochastic simulation algorithm.
An exact method for stochastic simulation of chemical reaction networks, which accelerates the stochastic simulation algorithm (SSA), is proposed. The present "ER-leap" algorithm is derived from analytic upper and lower bounds on the multireaction probabilities sampled by SSA, together with rejection sampling and an adaptive multiplicity for reactions. The algorithm is tested on a number of well-quantified reaction networks and is found experimentally to be very accurate on test problems including a chaotic reaction network. At the same time ER-leap offers a substantial speedup over SSA with a simulation time proportional to the 23 power of the number of reaction events in a Galton-Watson process.
Department of Computer Science, University of California, Irvine, U.C. Irvine, California 92697, USA. emj@uci.edu
Eric
E
Mjolsness
David
D
Orendorff
Philippe
P
Chatelain
Petros
P
Koumoutsakos
D000465:Algorithms; D003198:Computer Simulation; D008956:Models, Chemical; D013269:Stochastic Processes