J. R. Trail

Smooth relativistic Hartree-Fock pseudopotentials for H to Ba and Lu to Hg

We report smooth relativistic Hartree-Fock pseudopotentials (also known as averaged relativistic effective potentials) and spin-orbit operators for the atoms H to Ba and Lu to Hg. We remove the unphysical extremely nonlocal behavior resulting from the exchange interaction in a controlled manner, and represent the resulting pseudopotentials in an analytic form suitable for use within standard quantum chemistry codes. These pseudopotentials are suitable for use within Hartree-Fock and correlated wave function methods, including diffusion quantum Monte Carlo calculations.

Norm-conserving Hartree-Fock pseudopotentials and their asymptotic behavior.

We investigate the properties of norm-conserving pseudopotentials (effective core potentials) generated by inversion of the Hartree-Fock equations. In particular, we investigate the asymptotic behavior as r-->infinity and find that such pseudopotentials are nonlocal over all space, apart from a few special cases such as H and He. Such extreme nonlocality leads to a lack of transferability and, within periodic boundary conditions, an undefined total energy. The extreme nonlocality must therefore be removed, and we argue that the best way to accomplish this is a minor relaxation of the norm-conservation condition. This is implemented, and pseudopotentials for the atoms H-Ar are constructed and tested.

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.

Exact and approximate Kohn-Sham potentials in ensemble density-functional theory

We construct exact Kohn-Sham potentials for the ensemble density-functional theory (EDFT) from the ground and excited states of helium. The exchange-correlation (XC) potential is compared with the quasi-localdensity approximation and both single-determinant and symmetry-eigenstate ghost-corrected exact exchange approximations. Symmetry-eigenstate Hartree exchange recovers distinctive features of the exact XC potential and is used to calculate the correlation potential. Unlike the exact case, excitation energies calculated from these approximations depend on ensemble weight, and it is shown that only the symmetry-eigenstate method produces an ensemble derivative discontinuity. Differences in asymptotic and near-ground-state behavior of exact and approximate XC potentials are discussed in the context of producing accurate optical gaps.

Diffusion quantum Monte Carlo study of three-dimensional Wigner crystals

We report diffusion quantum Monte Carlo calculations of three-dimensional Wigner crystals in the density range rs=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 rs=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.

Excitations and benchmark ensemble density functional theory for two electrons

A new method for extracting ensemble Kohn-Sham potentials from accurate excited state densities is applied to a variety of two-electron systems, exploring the behavior of exact ensemble density functional theory. The issue of separating the Hartree energy and the choice of degenerate eigenstates is explored. A new approximation, spin eigenstate Hartree-exchange, is derived. Exact conditions that are proven include the signs of the correlation energy components and the asymptotic behavior of the potential for small weights of the excited states. Many energy components are given as a function of the weights for two electrons in a one-dimensional flat box, in a box with a large barrier to create charge transfer excitations, in a three-dimensional harmonic well (Hookes atom), and for the He atom singlet-triplet ensemble, singlet-triplet-singlet ensemble, and triplet bi-ensemble.

Spectroscopic data for the LiH molecule from pseudopotential quantum Monte Carlo calculations

Quantum Monte Carlo and quantum chemistry techniques are used to investigate pseudopotential models of the lithium hydride (LiH) molecule. Interatomic potentials are calculated and tested by comparing with the experimental spectroscopic constants and well depth. Two recently developed pseudopotentials are tested, and the effects of introducing a Li core polarization potential are investigated. The calculations are sufficiently accurate to isolate the errors from the pseudopotentials and core polarization potential. Core-valence correlation and core relaxation are found to be important in determining the interatomic potential.

Pseudopotentials for correlated electron systems

A scheme is developed for creating pseudopotentials for use in correlated-electron calculations. Pseudopotentials for the light elements H, Li, Be, B, C, N, O, and F, are reported, based on data from high-level quantum chemical calculations. Results obtained with these correlated electron pseudopotentials (CEPPs) are compared with data for atomic energy levels and the dissociation energies, molecular geometries, and zero-point vibrational energies of small molecules obtained from coupled cluster single double triple calculations with large basis sets. The CEPPs give better results in correlated-electron calculations than Hartree-Fock-based pseudopotentials available in the literature.

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).