Assessing the accuracy of quantum Monte Carlo and density functional theory for energetics of small water clusters

We present a detailed study of the energetics of water clusters (H$_2$O)$_n$ with $n le 6$, comparing diffusion Monte Carlo (DMC) and approximate density functional theory (DFT) with well converged coupled-cluster benchmarks. We use the many-body decomposition of the total energy to classify the errors of DMC and DFT into 1-body, 2-body and beyond-2-body components. Using both equilibrium cluster configurations and thermal ensembles of configurations, we find DMC to be uniformly much more accurate than DFT, partly because some of the approximate functionals give poor 1-body distortion energies. Even when these are corrected, DFT remains considerably less accurate than DMC. When both 1- and 2-body errors of DFT are corrected, some functionals compete in accuracy with DMC; however, other functionals remain worse, showing that they suffer from significant beyond-2-body errors. Combining the evidence presented here with the recently demonstrated high accuracy of DMC for ice structures, we suggest how DMC can now be used to provide benchmarks for larger clusters and for bulk liquid water.

Scheme for adding electron-nucleus cusps to Gaussian orbitals

A simple scheme is described for introducing the correct cusps at nuclei into orbitals obtained from Gaussian basis set electronic structure calculations. The scheme is tested with all-electron variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods for the Ne atom, the H2 molecule, and 55 molecules from a standard benchmark set. It greatly reduces the variance of the local energy in all cases and slightly improves the variational energy. This scheme yields a general improvement in the efficiency of all-electron VMC and DMC calculations using Gaussian basis sets.

All-electron quantum Monte Carlo calculations for the noble gas atoms He to Xe

We report all-electron variational and diffusion quantum Monte Carlo (VMC and DMC) calculations for the noble gas atoms He, Ne, Ar, Kr, and Xe. The calculations were performed using Slater-Jastrow wave functions with Hartree-Fock single-particle orbitals. The quality of both the optimized Jastrow factors and the nodal surfaces of the wave functions declines with increasing atomic number Z, but the DMC calculations are tractable and well behaved in all cases. We discuss the scaling of the computational cost of DMC calculations with Z.

Jastrow correlation factor for atoms, molecules, and solids

A form of Jastrow factor is introduced for use in quantum Monte Carlo simulations of finite and periodic systems. Test data are presented for atoms, molecules, and solids, including both all-electron and pseudopotential atoms. We demonstrate that our Jastrow factor is able to retrieve a large fraction of the correlation energy.