We analyze the problem of eliminating finite-size errors from quantum Monte Carlo (QMC) energy data. We demonstrate that both (i) adding a recently proposed [S. Chiesa et al., Phys. Rev. Lett. 97, 076404 (2006)] finite-size correction to the Ewald en
ergy and (ii) using the model periodic Coulomb (MPC) interaction [L. M. Fraser et al., Phys. Rev. B 53, 1814 (1996); P. R. C. Kent et al., Phys. Rev. B 59, 1917 (1999); A. J. Williamson et al., Phys. Rev. B 55, 4851 (1997)] are good solutions to the problem of removing finite-size effects from the interaction energy in cubic systems, provided the exchange-correlation (XC) hole has converged with respect to system size. However, we find that the MPC interaction distorts the XC hole in finite systems, implying that the Ewald interaction should be used to generate the configuration distribution. The finite-size correction of Chiesa et al. is shown to be incomplete in systems of low symmetry. Beyond-leading-order corrections to the kinetic energy are found to be necessary at intermediate and high densities, and we investigate the effect of adding such corrections to QMC data for the homogeneous electron gas. We analyze finite-size errors in two-dimensional systems and show that the leading-order behavior differs from that which has hitherto been supposed. We compare the efficiency of different twist-averaging methods for reducing single-particle finite-size errors and we examine the performance of various finite-size extrapolation formulas. Finally, we investigate the system-size scaling of biases in diffusion QMC.
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 quan
tum 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.
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 orbi
tals. 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.
We present density-functional theory (DFT) and quantum Monte Carlo (QMC) calculations designed to resolve experimental and theoretical controversies over the optical properties of H-terminated C nanoparticles (diamondoids). The QMC results follow the
trends of well-converged plane-wave DFT calculations for the size dependence of the optical gap, but they predict gaps that are 1-2 eV higher. They confirm that quantum confinement effects disappear in diamondoids larger than 1 nm, which have gaps below that of bulk diamond. Our QMC calculations predict a small exciton binding energy and a negative electron affinity (NEA) for diamondoids up to 1 nm, resulting from the delocalized nature of the lowest unoccupied molecular orbital. The NEA suggests a range of possible applications of diamondoids as low-voltage electron emitters.
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.
We report diffusion quantum Monte Carlo calculations of three-dimensional Wigner crystals in the density range r_s=100-150. We have tested different types of orbital for use in the approximate wave functions but none improve upon the simple Gaussian
form. The Gaussian exponents are optimized by directly minimizing the diffusion quantum Monte Carlo energy. We have carefully investigated and sought to minimize the potential biases in our Monte Carlo results. We conclude that the uniform electron gas undergoes a transition from a ferromagnetic fluid to a body-centered-cubic Wigner crystal at r_s=106+/-1. The diffusion quantum Monte Carlo results are compared with those from Hartree-Fock and Hartree theory in order to understand the role played by exchange and correlation in Wigner crystals. We also study floating Wigner crystals and give results for their pair-correlation functions.
We report calculations of the energies of excitons and biexcitons in ideal two-dimensional bilayer systems within the effective-mass approximation with isotropic electron and hole masses. The exciton energies are obtained by a simple numerical integr
ation technique, while the biexciton energies are obtained from diffusion quantum Monte Carlo calculations. The exciton binding energy decays as the inverse of the separation of the layers, while the binding energy of the biexciton with respect to dissociation into two separate excitons decays exponentially.
Quantum Monte Carlo (QMC) methods have been used to obtain accurate binding-energy data for pairs of parallel thin metallic wires and layers modeled by 1D and 2D homogeneous electron gases. We compare our QMC binding energies with results obtained wi
thin the random phase approximation, finding significant quantitative differences and disagreement over the asymptotic behavior for bilayers at low densities. We have calculated pair-correlation functions for metallic biwire and bilayer systems. Our QMC data could be used to investigate van der Waals energy functionals.
