Introduction to the Diffusion Monte Carlo Method

In this review we discuss, from a unified point of view, a variety of Monte Carlo methods used to solve eigenvalue problems in statistical mechanics and quantum mechanics. Although the applications of these methods differ widely, the underlying mathematics is quite similar in that they are stochastic implementations of the power method. In all cases, optimized trial states can be used to reduce the errors of Monte Carlo estimates.

The diffusion quantum Monte Carlo method is extended to solve the old theoretical physics problem of many-electron atoms and ions in intense magnetic fields. The feature of our approach is the use of adiabatic approximation wave functions augmented by a Jastrow factor as guiding functions to initialize the quantum Monte Carlo prodecure. We calcula te the ground state energies of atoms and ions with nuclear charges from Z= 2, 3, 4, ..., 26 for magnetic field strengths relevant...

We describe QWalk, a new computational package capable of performing Quantum Monte Carlo electronic structure calculations for molecules and solids with many electrons. We describe the structure of the program and its implementation of Quantum Monte Carlo methods. It is open-source, licensed under the GPL, and available at the web site http://www.qwalk.org

Diffusion Monte Carlo (DMC) calculations typically yield highly accurate results in solid-state and quantum-chemical calculations. However, operators that do not commute with the Hamiltonian are at best sampled correctly up to second order in the error of the underlying trial wavefunction, once simple corrections have been applied. This error is of the same order as that for the energy in variational calculations. Operators that suffer from these problems include potential en...

A recently developed self-healing diffusion Monte Carlo algorithm [PRB 79, 195117] is extended to the calculation of excited states. The formalism is based on an excited-state fixed-node approximation and the mixed estimator of the excited-state probability density. The fixed-node ground state wave-functions of inequivalent nodal pockets are found simultaneously using a recursive approach. The decay of the wave-function into lower energy states is prevented using two method...

In order to find the equilibrium geometries of molecules and solids and to perform ab initio molecular dynamics, it is necessary to calculate the forces on the nuclei. We present a correlated sampling method to efficiently calculate numerical forces and potential energy surfaces in diffusion Monte Carlo. It employs a novel coordinate transformation, earlier used in variational Monte Carlo, to greatly reduce the statistical error. Results are presented for first-row diatomic m...

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a di...

The implementation and reliability of a quadratic diffusion Monte Carlo method for the study of ground-state properties of atoms are discussed. We show in the simple yet non-trivial calculation of the binding energy of the Li atom that the method presented is effectively second-order in the time step. The fulfilment of the expected quadratic behavior relies on some basic requirements of the trial wave function used for importance sampling, in the context of the fixed-node app...

Treating the fermionic ground state problem as a constrained stochastic optimization problem, a formalism for fermionic quantum Monte Carlo is developed that makes no reference to a trial wavefunction. Exchange symmetry is enforced by nonlocal terms appearing in the Green's function corresponding to a new kind of walker propagation. Complemented by a treatment of diffusion that encourages the formation of a stochastic nodal surface, we find that an approximate long-range exte...

In this paper we propose an ab initio method to solve quantum many-body problems of molecular dynamics where both the electronic and the nuclear degrees are represented by ensembles of trajectories and guiding waves in physical space. Both electrons and nuclei can be treated quantum mechanically where the guiding waves obey a set of coupled Schrodinger equations (quantum-quantum description) or, alternatively, coupled Schroedinger-Newtonian equations are solved for the quantu...