Search

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

Rutile ($R$) phase VO$_2$ is a quintessential example of a strongly correlated bad-metal, which undergoes a metal-insulator transition (MIT) concomitant with a structural transition to a V-V dimerized monoclinic phase below T$_{MIT} \sim 340K$. It has been experimentally shown that one can control this transition by doping VO$_2$. In particular, doping with oxygen vacancies ($V_O$) has been shown to completely suppress this MIT {\em without} any structural transition. We explain this suppression by elucidating the influence of oxygen-vacancies on the electronic-structure of the metallic $R$ phase VO$_2$, explicitly treating strong electron-electron correlations using dynamical mean-field theory (DMFT) as well as diffusion Monte Carlo (DMC) flavor of quantum Monte Carlo (QMC) techniques. We show that $V_O$'s tend to change the V-3$d$ filling away from its nominal half-filled value, with the $e_{g}^{\pi}$ orbitals competing with the otherwise dominant $a_{1g}$ orbital. Loss of this near orbital polarization of the $a_{1g}$ orbital is associated with a weakening of electron correlations, especially along the V-V dimerization direction. This removes a charge-density wave (CDW) instability along this direction above a critical doping concentration, which further suppresses the metal-insulator transition. Our study also suggests that the MIT is predominantly driven by a correlation-induced CDW instability along the V-V dimerization direction.

QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org .

With recent developments in simulating nonadiabatic systems to high accuracy, it has become possible to determine how much energy is attributed to nuclear quantum effects beyond zero-point energy. In this work we calculate the non-relativistic ground-state energies of atomic and molecular systems without the Born-Oppenheimer approximation. For this purpose we utilize the fixed-node diffusion Monte Carlo method, in which the nodes depend on both the electronic and ionic positions. We report ground-state energies for all systems studied, ionization energies for the first-row atoms and atomization energies for the first-row hydrides. We find the ionization energies of the atoms to be nearly independent of the Born-Oppenheimer approximation, within the accuracy of our results. The atomization energies of molecular systems, however, show small effects of the nonadiabatic coupling between electrons and nuclei., Comment: 9 pages, 9 figures; corrected boron mass and error bars for nonadiabatic dissociation energies