Soft pseudopotentials for efficient quantum Monte Carlo calculations: From Be to Ne and Al to Ar

Correlation Consistent Basis Sets and Core Polarization Potentials for Al-Ar with ccECP Pseudopotentials

New correlation consistent basis sets for the second-row atoms (Al–Ar) to be used with the neon-core correlation consistent effective core potentials (ccECPs) have been developed. The basis sets, denoted cc-pV(n+d)Z-ccECP (n = D, T, Q), include the “tight”-d functions that are known to be important for second-row elements. Sets augmented with additional diffuse functions are also reported. Effective core polarization potentials (CPPs) to account for the effect of core–valence correlation have been adjusted for the same elements, and two different forms of the CPP cutoff function have been analyzed. The accuracy of both the basis sets and the CPPs is assessed through benchmark calculations at the coupled-cluster level of theory for atomic and molecular properties. Agreement with all-electron results is much improved relative to the basis sets that originally accompanied the ccECPs; moreover, the combination of cc-pV(n+d)Z-ccECP and CPPs is found to be a computationally efficient and accurate alternative to including core electrons in the correlation treatment.

High Accuracy Transition Metal Effective Cores for the Many-Body Diffusion Monte Carlo Method

Practical applications of the real-space diffusion Monte Carlo (DMC) method require the removal of core electrons, where currently localization approximations of semilocal potentials are generally used in the projector. Accurate calculations of complex solids and large molecules demand minimizing the impact of approximated atomic cores. Prior works have shown that the errors from such approximations can be sizable in both finite and periodic systems. In this work, we show that a class of differential pseudopotentials, known as pseudo-Hamiltonians, can be constructed for the 3d transition metal atoms, entirely removing the need for any localization scheme in the DMC projector. As a proof of principle, we demonstrate the approach for the case of Co. In order to minimize errors in the pseudo-Hamiltonian at the many-body level, we generalize the recently proposed correlation-consistent pseudopotential generation scheme to successively close semilocal representations of the differential potentials. Our generation scheme successfully produces potentials tailored specifically for real space projector quantum Monte Carlo methods with low error at the many-body level, i.e., with many-body scattering properties very close to relativistic all-electron results. In particular, we show that the agreement with respect to atomic and molecular quantities reach chemical accuracy in many cases─on par with the most accurate semilocal pseudopotentials available. Further, our pseudo-Hamiltonian generation scheme utilizes standard quantum chemistry codes designed only to work with semilocal pseudopotentials, enabling straightforward generation of pseudo-Hamiltonians for additional elements in future works.

Achieving high-quality single-atom nitrogen doping of graphene/SiC(0001) by ion implantation and subsequent thermal stabilization

We report a straightforward method to produce high-quality nitrogen-doped graphene on SiC(0001) using direct nitrogen ion implantation and subsequent stabilization at temperatures above 1300 K. We demonstrate that double defects, which comprise two nitrogen defects in a second-nearest-neighbor (meta) configuration, can be formed in a controlled way by adjusting the duration of bombardment. Two types of atomic contrast of single N defects are identified in scanning tunneling microscopy. We attribute the origin of these two contrasts to different tip structures by means of STM simulations. The characteristic dip observed over N defects is explained in terms of the destructive quantum interference.

Spin multiplicity and symmetry breaking in vanadium-benzene complexes

We present accurate quantum Monte Carlo (QMC) calculations which enabled us to determine the structure, spin multiplicity, ionization energy, dissociation energy, and spin-dependent electronic gaps of the vanadium-benzene system. From total and ionization energy we deduce a high-spin state with vastly different energy gaps for the two spin channels. For this purpose we have used a multistage combination of techniques with consecutive elimination of systematic biases except for the fixed-node approximation in QMC calculations. Our results significantly differ from the established picture based on previous less accurate calculations and point out the importance of high-level many-body methods for predictive calculations of similar transition metal-based organometallic systems.

Electron Structure Quantum Monte Carlo

The diffusion quantum Monte Carlo method (DMC) is capable of calculating accurately the electronic energy for molecules and molecular aggregates. An overview is given on recent developments for the optimization of the guide function that determines the accuracy of the method. Furthermore, the versatility of DMC is shown with applications to Rydberg states, transition metal compounds, and weakly interacting systems.

Quantum Monte Carlo study of porphyrin transition metal complexes

Diffusion quantum Monte Carlo (DMC) calculations for transition metal (M) porphyrin complexes (MPo, M=Ni,Cu,Zn) are reported. We calculate the binding energies of the transition metal atoms to the porphin molecule. Our DMC results are in reasonable agreement with those obtained from density functional theory calculations using the B3LYP hybrid exchange-correlation functional. Our study shows that such calculations are feasible with the DMC method.

Toward the exact solution of the electronic Schrödinger equation for noncovalent molecular interactions: worldwide distributed quantum monte carlo calculations

Quantum Monte Carlo (QMC) calculations on the stacked (st) and Watson/Crick (wc) bound adenine/thymine (A/T) and cytosine/guanine (C/G) DNA base pair complexes were made possible with the first large scale distributed computing project in ab initio quantum chemistry, Quantum Monte Carlo at Home (QMC@HOME). The results for the interaction energies (wc-A/T = 15.7 kcal/mol, wc-C/G = 30.2 kcal/mol, st-A/T = 13.1 kcal/mol, st-C/G = 19.6 kcal/mol) are in very good agreement with the best known coupled-cluster based estimates. The accuracy of these values is further supported by calculations on the S22 benchmark set of noncovalently bound systems, for which we obtain a small mean absolute deviation of 0.68 kcal/mol. Our results support previous claims that the stacking energies are of comparable magnitude to the interactions of the commonly discussed hydrogen-bonded motif. Furthermore, we show that QMC can serve as an advantageous alternative to conventional wave function methods for large noncovalently bound systems. We also investigated in detail all technical parameters of the QMC simulations and recommend a careful optimization procedure of the Jastrow correlation factors in order to obtain numerically stable and reliable results.

Weak intermolecular interactions calculated with diffusion Monte Carlo

The performance of fixed node diffusion Monte Carlo (FNDMC) for weakly interacting molecules is investigated. The effect of Gaussian basis sets on the asymptotic description of the molecular orbitals which is crucial for a successful importance sampling is analyzed for the example of the hydrogen atom. We find that accurate reference binding energies of the water, the ammonia, and the T-shaped as well as the parallel-displaced benzene dimer are correctly reproduced by FNDMC. The binding energies for the benzene dimers are -3.00(0.38) and -3.58(0.38) kcal/mol, respectively. The description of the methane dimer which has the smallest binding energy and a quite large intermolecular distance requires a more flexible basis set of diffuse quadruple-zeta quality in order to prevent sampling errors.

Energy-consistent pseudopotentials for quantum Monte Carlo calculations

The authors present scalar-relativistic energy-consistent Hartree-Fock pseudopotentials for the main-group elements. The pseudopotentials do not exhibit a singularity at the nucleus and are therefore suitable for quantum Monte Carlo (QMC) calculations. They demonstrate their transferability through extensive benchmark calculations of atomic excitation spectra as well as molecular properties. In particular, they compute the vibrational frequencies and binding energies of 26 first- and second-row diatomic molecules using post-Hartree-Fock methods, finding excellent agreement with the corresponding all-electron values. They also show their pseudopotentials give superior accuracy than other existing pseudopotentials constructed specifically for QMC. Finally, valence basis sets of different sizes (VnZ with n=D,T,Q,5 for first and second rows, and n=D,T for third to fifth rows) optimized for our pseudopotentials are also presented.

Energetics and dipole moment of transition metal monoxides by quantum Monte Carlo

The transition metal (TM) oxygen bond appears very prominently throughout chemistry and solid-state physics. Many materials, from biomolecules to ferroelectrics to the components of supernova remnants, contain this bond in some form. Many of these materials' properties depend strongly on fine details of the TM-O bond, which makes accurate calculations of their properties very challenging. Here the authors report on highly accurate first principles calculations of the properties of TM monoxide molecules within fixed-node diffusion Monte Carlo and reptation Monte Carlo.

Pfaffian pairing wave functions in electronic-structure quantum Monte Carlo simulations

We investigate the accuracy of trial wave functions for quantum Monte Carlo based on Pfaffian functional form with singlet and triplet pairing. Using a set of first row atoms and molecules we find that these wave functions provide very consistent and systematic behavior in recovering the correlation energies on the level of 95%. In order to get beyond this limit we explore the possibilities of multi-Pfaffian pairing wave functions. We show that a small number of Pfaffians recovers another large fraction of the missing correlation energy comparable to the larger-scale configuration interaction wave functions. We also find that Pfaffians lead to substantial improvements in fermion nodes when compared to Hartree-Fock wave functions.

Diffusion Monte Carlo method with lattice regularization

We introduce an efficient lattice regularization scheme for quantum Monte Carlo calculations of realistic electronic systems. The kinetic term is discretized by a finite difference Laplacian with two mesh sizes, a and a', chosen so that the electrons can diffuse in a configuration space which is in practice indistinguishable from the continuum, and the different length scales in the system can be efficiently taken in account. The regularized Hamiltonian goes to the continuous limit for a --> 0 and allows the inclusion of nonlocal potentials in a consistent variational scheme, substantially improving the accuracy upon previous nonvariational approaches.

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.

Performance of diffusion Monte Carlo for the first dissociation energies of transition metal carbonyls

Fixed node diffusion Monte Carlo (FNDMC) calculations are carried out for the first ligand dissociation energies of the prototype transition metal carbonyls Cr(CO)6, Fe(CO)5, Ni(CO)4, and Fe(CO)4N2. Since Hartree-Fock theory performs particularly badly for these type of compounds they are difficult to treat with conventional ab initio methods. We find that a Kohn-Sham determinant from a standard density functional provides a balanced description of the fermionic nodal hyper surfaces of all compounds involved in the dissociation reaction. With one exception, the experimental dissociation enthalpies are reproduced by FNDMC within the statistical accuracy of the method.

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.