A hierarchical exact accelerated stochastic simulation algorithm
Abstract
A new algorithm, "HiER-leap" (hierarchical exact reaction-leaping), is derived which improves on the computational properties of the ER-leap algorithm for exact accelerated simulation of stochastic chemical kinetics. Unlike ER-leap, HiER-leap utilizes a hierarchical or divide-and-conquer organization of reaction channels into tightly coupled "blocks" and is thereby able to speed up systems with many reaction channels. Like ER-leap, HiER-leap is based on the use of upper and lower bounds on the reaction propensities to define a rejection sampling algorithm with inexpensive early rejection and acceptance steps. But in HiER-leap, large portions of intra-block sampling may be done in parallel. An accept/reject step is used to synchronize across blocks. This method scales well when many reaction channels are present and has desirable asymptotic properties. The algorithm is exact, parallelizable and achieves a significant speedup over the stochastic simulation algorithm and ER-leap on certain problems. This algorithm offers a potentially important step towards efficient in silico modeling of entire organisms.
- Publication:
-
Journal of Chemical Physics
- Pub Date:
- December 2012
- DOI:
- 10.1063/1.4766353
- arXiv:
- arXiv:1212.4080
- Bibcode:
- 2012JChPh.137u4104O
- Keywords:
-
- reaction kinetics theory;
- stochastic processes;
- 82.20.Wt;
- 02.50.Fz;
- 82.20.Uv;
- Computational modeling;
- simulation;
- Stochastic analysis;
- Stochastic theories of rate constants;
- Quantitative Biology - Molecular Networks;
- Computer Science - Computational Engineering;
- Finance;
- and Science;
- Computer Science - Data Structures and Algorithms
- E-Print:
- 22 pages, 3 figures