Introduction to the Diffusion Monte Carlo Method

Quantum Monte Carlo (QMC) methods represent a powerful family of computational techniques for tackling complex quantum many-body problems and performing calculations of stationary state properties. QMC is among the most accurate and powerful approaches to the study of electronic structure, but its application is often hindered by a steep learning curve, hence it is rarely addressed in undergraduate and postgraduate classes. This tutorial is a step towards filling this gap. We...

Direct sampling of multi-dimensional systems with quantum Monte Carlo methods allows exact account of many-body effects or particle correlations. The most straightforward approach to solve the Schr\"odinger equation, Diffusion Monte Carlo, has been used in several benchmark cases for other methods to pursue. Its robustness is based on direct sampling of a positive probability density for diffusion in imaginary time. It has been argued that the corresponding real time diffusio...

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-dimen...

Quantum Monte Carlo methods are first-principle approaches that approximately solve the Schr\"odinger equation stochastically. As compared to traditional quantum chemistry methods, they offer important advantages such as the ability to handle a large variety of many-body wave functions, the favorable scaling with the number of particles, and the intrinsic parallelism of the algorithms which are particularly suitable to modern massively parallel computers. In this chapter, we ...

The computational cost of a Monte Carlo algorithm can only be meaningfully discussed when taking into account the magnitude of the resulting statistical error. Aiming for a fixed error per particle, we study the scaling behavior of the diffusion Monte Carlo method for large quantum systems. We identify the correlation within the population of walkers as the dominant scaling factor for large systems. While this factor is negligible for small and medium sized systems that are t...

We present an overview of the variational and diffusion quantum Monte Carlo methods as implemented in the CASINO program. We particularly focus on developments made in the last decade, describing state-of-the-art quantum Monte Carlo algorithms and software and discussing their strengths and their weaknesses. We review a range of recent applications of CASINO.

In this work we develop tools that enable the study of non-adiabatic effects with variational and diffusion Monte Carlo methods. We introduce a highly accurate wave function ansatz for electron-ion systems that can involve a combination of both fixed and quantum ions. We explicitly calculate the ground state energies of H$_{2}$, LiH, H$_{2}$O and FHF$^{-}$ using fixed-node quantum Monte Carlo with wave function nodes that explicitly depend on the ion positions. The obtained e...

This topical review describes the methodology of continuum variational and diffusion quantum Monte Carlo calculations. These stochastic methods are based on many-body wave functions and are capable of achieving very high accuracy. The algorithms are intrinsically parallel and well-suited to petascale computers, and the computational cost scales as a polynomial of the number of particles. A guide to the systems and topics which have been investigated using these methods is giv...

We study one-dimensional (1D) and two-dimensional (2D) Helium atoms using a new time-dependent quantum Monte Carlo (TDQMC) method. The TDQMC method employs random walkers, with a separate guiding wave attached to each walker. The ground state is calculated by a self-consistent solution of complex-time Schroedinger equations for the guiding waves and of equations for the velocity fields of the walkers. Our results show that the many-body wavefunction and the ground state energ...

We provide a pedagogical introduction to the two main variants of real-space quantum Monte Carlo methods for electronic-structure calculations: variational Monte Carlo (VMC) and diffusion Monte Carlo (DMC). Assuming no prior knowledge on the subject, we review in depth the Metropolis-Hastings algorithm used in VMC for sampling the square of an approximate wave function, discussing details important for applications to electronic systems. We also review in detail the more soph...