Systematic reduction of sign errors in many-body calculations of atoms and molecules

The self-healing diffusion Monte Carlo algorithm (SHDMC) [Phys. Rev. B {\bf 79}, 195117 (2009), {\it ibid.} {\bf 80}, 125110 (2009)] is shown to be an accurate and robust method for calculating the ground state of atoms and molecules. By direct comparison with accurate configuration interaction results for the oxygen atom we show that SHDMC converges systematically towards the ground-state wave function. We present results for the challenging N$_2$ molecule, where the binding energies obtained via both energy minimization and SHDMC are near chemical accuracy (1 kcal/mol). Moreover, we demonstrate that SHDMC is robust enough to find the nodal surface for systems at least as large as C$_{20}$ starting from random coefficients. SHDMC is a linear-scaling method, in the degrees of freedom of the nodes, that systematically reduces the fermion sign problem.Comment: Final version accepted in Physical Review Letters. The review history (referees' comments and our replies) is included in the source

Comparative study of density functional theories of the exchange-correlation hole and energy in silicon

We present a detailed study of the exchange-correlation hole and exchange-correlation energy per particle in the Si crystal as calculated by the Variational Monte Carlo method and predicted by various density functional models. Nonlocal density averaging methods prove to be successful in correcting severe errors in the local density approximation (LDA) at low densities where the density changes dramatically over the correlation length of the LDA hole, but fail to provide systematic improvements at higher densities where the effects of density inhomogeneity are more subtle. Exchange and correlation considered separately show a sensitivity to the nonlocal semiconductor crystal environment, particularly within the Si bond, which is not predicted by the nonlocal approaches based on density averaging. The exchange hole is well described by a bonding orbital picture, while the correlation hole has a significant component due to the polarization of the nearby bonds, which partially screens out the anisotropy in the exchange hole.Comment: 16 pages, 5 figures, RevTeX, added conten

Diffusion Quantum Monte Carlo Calculations of Excited States of Silicon

The band structure of silicon is calculated at the Gamma, X, and L wave vectors using diffusion quantum Monte Carlo methods. Excited states are formed by promoting an electron from the valence band into the conduction band. We obtain good agreement with experiment for states around the gap region and demonstrate that the method works equally well for direct and indirect excitations, and that one can calculate many excited states at each wave vector. This work establishes the fixed-node DMC approach as an accurate method for calculating the energies of low lying excitations in solids.Comment: 5 pages, 1 figur

Two-dimensional limit of exchange-correlation energy functional approximations in density functional theory

We investigate the behavior of three-dimensional (3D) exchange-correlation energy functional approximations of density functional theory in anisotropic systems with two-dimensional (2D) character. Using two simple models, quasi-2D electron gas and two-electron quantum dot, we show a {\it fundamental limitation} of the local density approximation (LDA), and its semi-local extensions, generalized gradient approximation (GGA) and meta-GGA (MGGA), the most widely used forms of which are worse than the LDA in the strong 2D limit. The origin of these shortcomings is in the inability of the local (LDA) and semi-local (GGA/MGGA) approximations to describe systems with 2D character in which the nature of the exchange-correlation hole is very nonlocal. Nonlocal functionals provide an alternative approach, and explicitly the average density approximation (ADA) is shown to be remarkably accurate for the quasi-2D electron gas system. Our study is not only relevant for understanding of the functionals but also practical applications to semiconductor quantum structures and materials such as graphite and metal surfaces. We also comment on the implication of our findings to the practical device simulations based on the (semi-)local density functional method.Comment: 21 pages including 9 figures, to be published in Phys. Rev.

Quantum Monte Carlo calculations of the one-body density matrix and excitation energies of silicon

Quantum Monte Carlo (QMC) techniques are used to calculate the one-body density matrix and excitation energies for the valence electrons of bulk silicon. The one-body density matrix and energies are obtained from a Slater-Jastrow wave function with a determinant of local density approximation (LDA) orbitals. The QMC density matrix evaluated in a basis of LDA orbitals is strongly diagonally dominant. The natural orbitals obtained by diagonalizing the QMC density matrix resemble the LDA orbitals very closely. Replacing the determinant of LDA orbitals in the wave function by a determinant of natural orbitals makes no significant difference to the quality of the wave function's nodal surface, leaving the diffusion Monte Carlo energy unchanged. The Extended Koopmans' Theorem for correlated wave functions is used to calculate excitation energies for silicon, which are in reasonable agreement with the available experimental data. A diagonal approximation to the theorem, evaluated in the basis of LDA orbitals, works quite well for both the quasihole and quasielectron states. We have found that this approximation has an advantageous scaling with system size, allowing more efficient studies of larger systems.Comment: 13 pages, 4 figures. To appear in Phys. Rev.

Finite size errors in quantum many-body simulations of extended systems

Further developments are introduced in the theory of finite size errors in quantum many-body simulations of extended systems using periodic boundary conditions. We show that our recently introduced Model Periodic Coulomb interaction [A. J. Williamson et al., Phys. Rev. B 55, R4851 (1997)] can be applied consistently to all Coulomb interactions in the system. The Model Periodic Coulomb interaction greatly reduces the finite size errors in quantum many-body simulations. We illustrate the practical application of our techniques with Hartree-Fock and variational and diffusion quantum Monte Carlo calculations for ground and excited state calculations. We demonstrate that the finite size effects in electron promotion and electron addition/subtraction excitation energy calculations are very similar.Comment: 15 pages, 6 figures. To appear in Phys. Rev.