Diffusion Monte Carlo method: numerical analysis in a simple case
Abstract
The Diffusion Monte Carlo method is devoted to the computation of electronic ground-state energies of molecules. In this paper, we focus on implementations of this method which consist in exploring the configuration space with a {\bf fixed} number of random walkers evolving according to a Stochastic Differential Equation discretized in time. We allow stochastic reconfigurations of the walkers to reduce the discrepancy between the weights that they carry. On a simple one-dimensional example, we prove the convergence of the method for a fixed number of reconfigurations when the number of walkers tends to while the timestep tends to 0. We confirm our theoretical rates of convergence by numerical experiments. Various resampling algorithms are investigated, both theoretically and numerically
Cite
@article{arxiv.math/0703407,
title = {Diffusion Monte Carlo method: numerical analysis in a simple case},
author = {Tony Lelievre and Mohamed El Makrini and Benjamin Jourdain},
journal= {arXiv preprint arXiv:math/0703407},
year = {2007}
}