P. Lopez Rios

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.

Inhomogeneous backflow transformations in quantum Monte Carlo calculations

An inhomogeneous backflow transformation for many-particle wave functions is presented and applied to electrons in atoms, molecules, and solids. We report variational and diffusion quantum Monte Carlo (VMC and DMC) energies for various systems and study the computational cost of using backflow wave functions. We find that inhomogeneous backflow transformations can provide a substantial increase in the amount of correlation energy retrieved within VMC and DMC calculations. The backflow transformations significantly improve the wave functions and their nodal surfaces.

Energies of the first row atoms from quantum Monte Carlo

All-electron variational and diffusion quantum Monte Carlo calculations of the ground state energies of the first row atoms (from Li to Ne) are reported. The authors use trial wave functions of four types: single-determinant Slater-Jastrow wave functions, multideterminant Slater-Jastrow wave functions, single-determinant Slater-Jastrow wave functions with backflow transformations, and multideterminant Slater-Jastrow wave functions with backflow transformations. At the diffusion quantum Monte Carlo level and using their multideterminant Slater-Jastrow wave functions with backflow transformations, they recover 99% or more of the correlation energies for Li, Be, B, C, N, and Ne, 97% for O, and 98% for F.

Quantum Monte Carlo study of the Ne atom and the Ne+ ion

We report all-electron and pseudopotential calculations of the ground-state energies of the neutral Ne atom and the Ne(+) ion using the variational and diffusion quantum Monte Carlo (DMC) methods. We investigate different levels of Slater-Jastrow trial wave function: (i) using Hartree-Fock orbitals, (ii) using orbitals optimized within a Monte Carlo procedure in the presence of a Jastrow factor, and (iii) including backflow correlations in the wave function. Small reductions in the total energy are obtained by optimizing the orbitals, while more significant reductions are obtained by incorporating backflow correlations. We study the finite-time-step and fixed-node biases in the DMC energy and show that there is a strong tendency for these errors to cancel when the first ionization potential (IP) is calculated. DMC gives highly accurate values for the IP of Ne at all the levels of trial wave function that we have considered.

Quantum Monte Carlo study of the first-row atoms and ions

Quantum Monte Carlo calculations of the first-row atoms Li-Ne and their singly positively charged ions are reported. Multideterminant-Jastrow-backflow trial wave functions are used which recover more than 98% of the correlation energy at the variational Monte Carlo level and more than 99% of the correlation energy at the diffusion Monte Carlo level for both the atoms and ions. We obtain the first ionization potentials to chemical accuracy. We also report scalar relativistic corrections to the energies, mass-polarization terms, and one- and two-electron expectation values.

Quantum Monte Carlo study of the phase diagram of solid molecular hydrogen at extreme pressures

Establishing the phase diagram of hydrogen is a major challenge for experimental and theoretical physics. Experiment alone cannot establish the atomic structure of solid hydrogen at high pressure, because hydrogen scatters X-rays only weakly. Instead, our understanding of the atomic structure is largely based on density functional theory (DFT). By comparing Raman spectra for low-energy structures found in DFT searches with experimental spectra, candidate atomic structures have been identified for each experimentally observed phase. Unfortunately, DFT predicts a metallic structure to be energetically favoured at a broad range of pressures up to 400 GPa, where it is known experimentally that hydrogen is non-metallic. Here we show that more advanced theoretical methods (diffusion quantum Monte Carlo calculations) find the metallic structure to be uncompetitive, and predict a phase diagram in reasonable agreement with experiment. This greatly strengthens the claim that the candidate atomic structures accurately model the experimentally observed phases.

Quantum Monte Carlo study of a positron in an electron gas

Quantum Monte Carlo calculations of the relaxation energy, pair-correlation function, and annihilating-pair momentum density are presented for a positron immersed in a homogeneous electron gas. We find smaller relaxation energies and contact pair-correlation functions in the important low-density regime than predicted by earlier studies. Our annihilating-pair momentum densities have almost zero weight above the Fermi momentum due to the cancellation of electron-electron and electron-positron correlation effects.

Framework for constructing generic Jastrow correlation factors

We have developed a flexible framework for constructing Jastrow factors which allows for the introduction of terms involving arbitrary numbers of particles. The use of various three- and four-body Jastrow terms in quantum Monte Carlo calculations is investigated, including a four-body van der Waals-like term, and anisotropic terms. We have tested these Jastrow factors on one- and two-dimensional homogeneous electron gases, the Be, B, and O atoms, and the BeH, H2O, N2, and H2 molecules. Our optimized Jastrow factors retrieve more than 90% of the fixed-node diffusion Monte Carlo correlation energy in variational Monte Carlo for each system studied.

Quantum Monte Carlo study of the energetics of the rutile, anatase, brookite, and columbite TiO

J.R.T., P.L.R., and R.J.N. acknowledge financial support from the Engineering and Physical Sciences Research Council (EPSRC) of the U.K. under Grant No. EP/J017639/1. B.M. acknowledges support from Robinson College, Cambridge, and the Cambridge Philosophical Society for a Henslow Research Fellowship. R.M. is grateful for financial support from MEXT-KAKENHI Grants No. 26287063, No. 25600156, and No. 22104011, and a grant from the Asahi Glass Foundation. Computational resources were provided by the Archer facility of the U.K.s national high-performance computing service (for which access was obtained via the UKCP consortium, EPSRC Grant No. EP/K014560/1), by the Center for Information Science of the JAIST, and by the K-computer (supported by the Computational Materials Science Initiative, CMSI/Japan, under Projects No. hp120086, No. hp140150, and No. hp150014).

_{2}

The interaction between a single hole and a two-dimensional, paramagnetic, homogeneous electron gas is studied using diffusion quantum Monte Carlo simulations. Electron-hole relaxation energies, pair-correlation functions, and electron-hole center-of-mass momentum densities are reported for a range of electron-hole mass ratios and electron densities. We find numerical evidence of a crossover from a collective excitonic state to a trion-dominated state in a density range in agreement with that found in recent experiments on quantum-well heterostructures.