Stochastic simulation of coupled reaction-diffusion processes (Q1924870)

Reaction-diffusion processes in spatially inhomogeneous systems are studied by an extension of the method proposed by \textit{D. T. Gillespie} [J. Comput. Phys. 22, 403-434 (1976; MR 58.20140)]. The treatment of diffusion enables to relate the diffusion coefficient to the corresponding transition rate probability in the multivariate master equation. Using the Green function the bulk diffusion coefficient in Fick's differential law is mapped. The method provides an approach to simulating stochastic reaction trajectories consistent with this master equation and allows the reaction-diffusion kinetics to be treated for low densities without imposing questionable continuum approximations. The biologically important nonlinear reaction-diffusion process of calcium wave propagation through the cell interior is studied.