Optimal finite-differences discretization for the diffusion equation from the perspective of large-deviation theory

Abstract

When applying the finite-differences method to numerically solve the one-dimensional diffusion equation, one must choose discretization steps Δx, Δt in space and time, respectively. By applying large-deviation theory on the discretized dynamics, we analyze the numerical errors due to the discretization, and find that the (relative) errors are especially large in regions of space where the concentration of particles is very small. We find that the choice Δt=Δx2/(6D) , where D is the diffusion coefficient, gives optimal accuracy compared to any other choice (including, in particular, the limit Δt→0 ), thus reproducing the known result that may be obtained using truncation error analysis. In addition, we give quantitative estimates for the dynamical lengthscale that describes the size of the spatial region in which the numerical solution is accurate, and study its dependence on the discretization parameters. We then turn to study the advection–diffusion equation, and obtain explicit expressions for the optimal Δt and other parameters of the finite-differences scheme, in terms of Δx, D and the advection velocity. We apply these results to study large deviations of the area swept by a diffusing particle in one dimension, trapped by an external potential ∼|x| . We extend our analysis to higher dimensions by combining our results from the one dimensional case with the locally one-dimension method.