Stochastic Simulation of Continuous Time Random Walks: Minimizing Error in Time-Dependent Rate Coefficients for Diffusion-Limited Reactions

A reaction limited by standard diffusion is simulated stochastically to illustrate how the continuous time random walk (CTRW) formalism can be implemented with minimum statistical error. A step-by-step simulation of the diffusive random walk in one dimension reveals the fraction of surviving reactants P(t) as a function of time, and the time-dependent unimolecular reaction rate coefficient K(t). Accuracy is confirmed by comparing the time-dependent simulation to results from the analytical master equation, and the asymptotic solution to that of Fickian diffusion. An early transient feature is shown to arise from higher spatial harmonics in the Fourier distribution of walkers between reaction sites. Statistical ‘shot’ noise in the simulation is quantified along with the offset error due to the discrete time derivative, and an optimal simulation time interval Δt0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta t_0$$\end{document} is derived to achieve minimal error in the finite time-difference estimation of the reaction rate. The number of walkers necessary to achieve a given error tolerance is derived, and W=107\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W = 10^7$$\end{document} walkers is shown to achieve an accuracy of ±0.2%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pm 0.2\%$$\end{document} when the survival probability reaches P(t)∼13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P(t) \sim \frac{1}{3}$$\end{document}. The stochastic method presented here serves as an intuitive basis for understanding the CTRW formalism, and can be generalized to model anomalous diffusion-limited reactions to prespecified precision in regimes where the governing wait-time distributions have no analytical solution.