Analysis of Brownian Dynamics Simulations of Reversible Bimolecular Reactions

A class of Brownian dynamics algorithms for stochastic reaction-diffusion models which include reversible bimolecular reactions is presented and analyzed. The method is a generalization of the $\lambda$–$\text{newrho}$ model for irreversible bimolecular reactions which was introduced in [arXiv:0903.1298]. The formulae relating the experimentally measurable quantities (reaction rate constants and diffusion constants) with the algorithm parameters are derived. The probability of geminate recombination is also investigated.