Comment on “Quantum Monte Carlo study of the dipole moment of CO” [J. Chem. Phys. 110, 11700 (1999)]

Zero-variance zero-bias quantum Monte Carlo estimators for the electron density at a nucleus

We derive new quantum Monte Carlo (QMC) estimators for the electronic density at the position of a point nucleus using the zero-variance and zero-bias principles. The resulting estimators are highly efficient, and are significantly simpler to implement and use than alternative methods, as they contain no adjustable parameters. In addition, they can be used in both variational and diffusion QMC calculations. Our best estimator is used to calculate the most accurate available estimates of the total electron density at the nucleus for the first-row atoms Li-Ne, the Ar atom, and the diatomic molecules B(2), N(2), and F(2).

Methods for calculating forces within quantum Monte Carlo simulations

Atomic force calculations within the variational and diffusion quantum Monte Carlo methods are described. The advantages of calculating diffusion quantum Monte Carlo forces with the 'pure' rather than the 'mixed' probability distribution are discussed. An accurate and practical method for calculating forces using the pure distribution is presented and tested for the SiH molecule. The statistics of force estimators are explored and violations of the central limit theorem are found in some cases.

Continuum variational and diffusion quantum Monte Carlo calculations

This topical review describes the methodology of continuum variational and diffusion quantum Monte Carlo calculations. These stochastic methods are based on many-body wavefunctions and are capable of achieving very high accuracy. The algorithms are intrinsically parallel and well suited to implementation on petascale computers, and the computational cost scales as a polynomial in the number of particles. A guide to the systems and topics which have been investigated using these methods is given. The bulk of the article is devoted to an overview of the basic quantum Monte Carlo methods, the forms and optimization of wavefunctions, performing calculations under periodic boundary conditions, using pseudopotentials, excited-state calculations, sources of calculational inaccuracy, and calculating energy differences and forces.

Accurate forces in quantum Monte Carlo calculations with nonlocal pseudopotentials

Calculating accurate forces within variational and diffusion Monte Carlo (VMC and DMC) methods is a very challenging problem. We derive expressions for the contribution to the Hellmann-Feynman force from nonlocal pseudopotentials for use within the VMC and DMC methods. Equilibrium bond lengths and harmonic vibrational frequencies are calculated from the Hellmann-Feynman forces and compared with those obtained from the energies at the Hartree-Fock, VMC, and pure DMC levels. Results for five small molecules show that the equilibrium bond lengths obtained from the force and energy calculations differ by less than 0.007 A at the DMC level.

Energetics and dipole moment of transition metal monoxides by quantum Monte Carlo

The transition metal (TM) oxygen bond appears very prominently throughout chemistry and solid-state physics. Many materials, from biomolecules to ferroelectrics to the components of supernova remnants, contain this bond in some form. Many of these materials' properties depend strongly on fine details of the TM-O bond, which makes accurate calculations of their properties very challenging. Here the authors report on highly accurate first principles calculations of the properties of TM monoxide molecules within fixed-node diffusion Monte Carlo and reptation Monte Carlo.

Dielectric response of periodic systems from quantum Monte Carlo calculations

We present a novel approach that allows us to calculate the dielectric response of periodic systems in the quantum Monte Carlo formalism. We employ a many-body generalization for the electric-enthalpy functional, where the coupling with the field is expressed via the Berry-phase formulation for the macroscopic polarization. A self-consistent local Hamiltonian then determines the ground-state wave function, allowing for accurate diffusion quantum Monte Carlo calculations where the polarization's fixed point is estimated from the average on an iterative sequence, sampled via forward walking. This approach has been validated for the case of an isolated hydrogen atom and then applied to a periodic system, to calculate the dielectric susceptibility of molecular-hydrogen chains. The results found are in excellent agreement with the best estimates obtained from the extrapolation of quantum-chemistry calculations.

Accurate, efficient, and simple forces computed with quantum Monte Carlo methods

Computation of ionic forces using quantum Monte Carlo methods has long been a challenge. We introduce a simple procedure, based on known properties of physical electronic densities, to make the variance of the Hellmann-Feynman estimator finite. We obtain very accurate geometries for the molecules H(2), LiH, CH(4), NH(3), H(2)O, and HF, with a Slater-Jastrow trial wave function. Harmonic frequencies for diatomics are also in good agreement with experiment. An antithetical sampling method is also discussed for additional reduction of the variance.