A Monte Carlo method for large reaction–diffusion systems

Abstract

We present a Monte Carlo algorithm which is specially designed for studying the stochastic behavior of large reaction–diffusion systems. We compare its performance and results to those obtained by a standard method for a solvable model system for which we calculate the mean concentration, the variance, and the spatial autocorrelation function of concentration.