Radek Erban and S Jonathan Chapman 2009 Phys. Biol. 6 046001 doi:10.1088/1478-3975/6/4/046001
Radek Erban and S Jonathan Chapman
Show affiliationsSeveral stochastic simulation algorithms (SSAs) have recently been proposed for modelling reaction–diffusion processes in cellular and molecular biology. In this paper, two commonly used SSAs are studied. The first SSA is an on-lattice model described by the reaction–diffusion master equation. The second SSA is an off-lattice model based on the simulation of Brownian motion of individual molecules and their reactive collisions. In both cases, it is shown that the commonly used implementation of bimolecular reactions (i.e. the reactions of the form A + B → C or A + A → C) might lead to incorrect results. Improvements of both SSAs are suggested which overcome the difficulties highlighted. In particular, a formula is presented for the smallest possible compartment size (lattice spacing) which can be correctly implemented in the first model. This implementation uses a new formula for the rate of bimolecular reactions per compartment (lattice site).
87.15.R- Reactions and kinetics
Issue 4 (December 2009)
Received 6 March 2009, accepted for publication 3 July 2009
Published 21 August 2009
Radek Erban and S Jonathan Chapman 2009 Phys. Biol. 6 046001
Yulan Yang and Qingyao Liu 2003 The Astronomical Journal 126 1960
N. R. Sheeley, Jr. and Y.-M. Wang 2007 ApJ 655 1142
Jacques P. Vallée 2008 The Astronomical Journal 135 1301
Dirk G Kurth 2008 Sci. Technol. Adv. Mater. 9 014103
Ravi Sankrit et al. 2007 The Astronomical Journal 133 1383
S B Whitfield et al 2007 J. Phys. B: At. Mol. Opt. Phys. 40 3647
Yaakov Garb et al 2008 Environ. Res. Lett. 3 045015
2008 J. Radiol. Prot. 28 613
B Bastiaensen et al 2007 Class. Quantum Grav. 24 3211