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.     

Relativistic small-core energy-consistent pseudopotentials for the alkaline-earth elements from Ca to Ra

Small-core ten-valence electron energy-consistent scalar- and two-component relativistic pseudopotentials for the alkaline-earth elements from Ca to Ra are presented. The accuracy and reliability of these pseudopotentials are discussed in terms of their ability to reproduce all-electron calculated and experimental dipole polarizabilities and ionization potentials.     

All-electron and relativistic pseudopotential studies for the group 1 element polarizabilities from K to element 119

Two-component and scalar relativistic energy-consistent pseudopotentials for the group 1 elements from K to element 119 are presented using nine electrons for the valence space definition. The accuracy of such an approximation is discussed for dipole polarizabilities and ionization potentials obtained at the coupled-cluster level as compared to experimental and all-electron Douglas-Kroll results.     

Quasirelativistic f-in-core pseudopotentials and core-polarization potentials for trivalent actinides and lanthanides: molecular test for trifluorides

Calibration studies of actinide and lanthanide trifluorides are reported for actinide and lanthanide scalar-relativistic energy-consistent f-in-core pseudopotentials, respectively, accompanying valence basis sets as well as core-polarization potentials. Results from Hartree–Fock and coupled-cluster singles, doubles, and perturbative triples f-in-core pseudopotential calculations are compared to corresponding data from f-in-valence pseudopotential and all-electron calculations as well as to experimental data. In general, good agreement is observed between the f-in-core and f-in-valence pseudopotential results, whereas due to the lack of experimental data for the actinides only a good agreement of the calculated and experimentally determined bond lengths of the lanthanide systems can be established. Nevertheless, the results indicate that the core-polarization potentials devised here for actinides improve the f-in-core results.     

Improved valence basis sets for divalent lanthanide 4f-in-core pseudopotentials

Improved energy-optimized (6s5p4d) and (7s6p5d) primitive valence basis sets have been derived for energy-consistent scalar-relativistic 4f-in-core pseudopotentials of the Stuttgart-Cologne variety modeling divalent lanthanides with a $4\hbox{f}^{n+1}$ occupation (n = 0?C13 for La?CYb). Segmented contracted basis sets covering the range of polarized double-, triple-, and quadruple-zeta quality, augmented by 2f1g correlation sets, were created for use in molecular calculations. The basis sets contain smaller (4s4p3d) and (5s5p4d) primitive subsets, which are designed in particular for solid state calculations of crystals containing divalent lanthanide ions. Hartree?CFock, density functional theory and coupled cluster results obtained with the new basis sets for lanthanide atomic ionization potentials as well as of geometry optimizations of various test molecules, i.e. selected lanthanide mono- and dihydrides, mono- and difluorides, and monooxides, show a satisfactory agreement with experimental data as well as with corresponding scalar-relativistic all-electron results. Core-polarization potentials are found to improve the results, especially for the atomic first and second ionization potentials.     

Relativistic energy-consistent pseudopotentials--recent developments

The direct adjustment of two-component pseudopotentials (scalar-relativistic + spin-orbit potentials), to atomic total energy valence spectra derived from four-component multiconfiguration Dirac-Hartree-Fock all-electron calculations based on the Dirac-Coulomb-Breit Hamiltonian, has been made a routine tool for an efficient treatment of heavy main-group elements. Both large-core (nsp valence shell) and small-core ((n - 1)spd nsp valence shell) potentials have been generated for all the post-d elements of groups 13-17. At the example of lead and bismuth compounds (PbHal, BiH, BiO, BiHal (Hal = F, Cl, Br, I)), we show how small-core and large-core potentials can be combined in accurate, yet computationally economic, spin-free-state-shifted relativistic electronic structure calculations of molecular ground and excited states.     

Development of new pseudopotential methods: improved model core potentials for the first-row transition metals

We have recently developed new nonrelativistic and scalar-relativistic pseudopotentials for the first-row transition metal and several main-group elements. These improved Model Core Potentials were tested on a variety of transition metal complexes to determine their accuracy in reproducing electronic structures, bond lengths, and harmonic vibrational frequencies with respect to both all-electron reference data as well as experimental data. The new potentials are also compared with the previous model core potentials available for the first-row transition metals. The new potentials do a superior job at reproducing atomic data, reproduce molecular data as well as the previous version, and in conjunction with new main-group pseudopotentials that have L-shell structure of the valence basis set, they are slightly faster.     

Energy-consistent relativistic pseudopotentials for the 4d elements: atomic and molecular applications

Recently reported energy-consistent relativistic pseudopotentials have been used with series of matching correlation consistent basis sets in benchmark calculations of various atomic and molecular properties. The basis set convergence of the 4d metal electron affinities and 5s2-->5s0 excitation energies are reported at the CCSD(T) level of theory, and the effects of valence and 4s4p correlation are investigated. In addition the impact of correlating the low-lying 3d electrons was also studied in all-electron Douglas-Kroll-Hess (DKH) calculations, which also included the ionization potentials and 5s2-->5s1 excitation energies. For all four atomic properties, higher order coupled cluster calculations through CCSDTQ are reported. The final calculated values are generally all within 1 kcal/mol of experiment. A notable exception is the ionization potential of Tc, the currently accepted experimental value of which is suggested to be too high by about 3 kcal/mol. Molecular calculations are also reported for the low-lying electronic states of ZrO and RuF, as well as the ground electronic state of Pd2. The effects of spin-orbit coupling are investigated for these cases in pseudopotential calculations. Wherever possible, the pseudopotential results have been calibrated against DKH calculations with correlation consistent basis sets of triple-zeta quality. In all cases the calculated data for these species are in very good agreement with experiment. In particular, the correct electronic ground state for the RuF molecule (4Phi92) was obtained, which was made possible by utilizing systematic sequences of correlation consistent basis sets.     

Fully relativistic pseudopotentials for alkaline atoms: Dirac-Hartree-Fock and configuration interaction calculations of alkaline monohydrides

Fully relativistic four-component energy-adjusted pseudopotentials and corresponding valence basis sets have been derived for the alkaline atoms Li through Cs, treating them as one-valence electron systems. Core-valence correlation effects are accounted for by a core-polarization potential, deviations of the core-nucleus repulsion from a point charge model by a suitable correction. The results of Dirac-Hartree-Fock and configuration interaction calculations are presented for atomic properties not used in the pseudopotential adjustment, i.e. electron affinities and dipole polarizabilities, as well as for the spectroscopic constants of the ground states of the alkaline monohydrides. The analytic form of the cut-off function for the electric field in the core-polarization term and its effects on atomic and molecular properties is discussed.     

A small-core multiconfiguration Dirac–Hartree–Fock-adjusted pseudopotential for Tl – application to TlX (X = F, Cl, Br, I)

 A relativistic pseudopotential of the energy-consistent variety simulating the Tl²¹⁺ (1s– 4f) core has been generated by adjustment to multiconfiguration Dirac–Hartree–Fock data based on the Dirac–Coulomb–Breit Hamiltonian. Valence ab initio calculations using this pseudopotential have been performed for atomic excitation energies and for spectroscopic constants of the X0⁺ and A0⁺ states of TlX (X = F, Cl, Br, I). Comparison is made to experiment and to four-component density functional results. Received: 22 June 1999 / Accepted: 5 August 1999 / Published online: 15 December 1999     

Pseudopotential study of the ground and excited states of Yb2

The ground state of the van der Waals-type lanthanide dimer Yb₂ has been studied by means of relativistic energy-consistent ab initio pseudopotentials using three different core definitions. Electron correlation was treated by coupled-cluster theory, whereby core-valence correlation effects have been accounted for either explicitly by correlating the energetically highest coreorbitals or implicitly by means of an effective core-polarization potential. Results for the first and second atomic ionization potentials, the atomic dipole polarizability, and the spectroscopic constants of the molecular ground state are reported. Low-lying excited states have been investigated with spin-orbit configuration interaction calculations. It is also demonstrated for the whole lanthanide series that correlation effects due to the atomic-like, possibly open 4f-shell in lanthanides can be modeled effectively by adding a core-polarization potential to pseudopotentials attributing the 4f-shell to the core. Received: 3 April 1998 / Accepted: 27 July 1998 / Published online: 9 October 1998     

Quasirelativistic energy-consistent 5f-in-core pseudopotentials for pentavalent and hexavalent actinide elements

Quasirelativistic energy-consistent 5f-in-core pseudopotentials modeling pentavalent (5f ^n?2 occupation with n = 2–6 for Pa–Am) and hexavalent (5f ^n?3 occupation with n = 3–6 for U–Am) actinides have been adjusted. Energy-optimized (6s5p4d) and (7s6p5d) valence basis sets contracted to polarized double- to quadruple-zeta quality as well as 2f1g correlation functions have been derived. Corresponding smaller basis sets (4s4p3d) and (5s5p4d) suitable for calculations on actinide(V) and actinide(VI) ions in crystalline solids form subsets of these basis sets designed for calculations on neutral molecules. Calculations using the Hartree–Fock and the coupled-cluster method with single and double excitation operators and a perturbative estimate of triple excitations for actinide pentafluorides show satisfactory agreement with calculations using 5f-in-valence pseudopotentials and experimental data, respectively. However, in the hexavalent case the 5f-in-core approximation seems to reach its limitations except for hexavalent uranium (5f⁰), where results for both uranium hexafluoride and the uranyl ion deviate only slightly from the 5f-in-valence reference data.     

Quasirelativistic energy-consistent 4f-in-core pseudopotentials for tetravalent lanthanide elements

Quasirelativistic energy-consistent 4f-in-core pseudopotentials modeling tetravalent lanthanides (4f ^n?1 occupation with n = 1, 2, 3, 8, 9 for Ce, Pr, Nd, Tb, Dy) have been adjusted. Energy-optimized (6s5p4d) and (7s6p5d) valence basis sets contracted to polarized double- to quadruple-zeta quality as well as 2f1g correlation functions have been derived. Corresponding smaller (4s4p3d) and (5s5p4d) basis sets suitable for calculations on lanthanide(IV) ions in crystalline solids form subsets of these basis sets designed for calculations on neutral molecules. Calculations for lanthanide tetrafluorides using the 4f-in-core pseudopotentials at the Hartree–Fock level show satisfactory agreement with calculations using 4f-in-valence pseudopotentials. For cerium tetrafluoride the experimental bond length is well reproduced using the 4f-in-core pseudopotential at the coupled-cluster level with single and double excitation operators and a perturbative estimate of triple excitations. For cerium dioxide 4f-in-core and 4f-in-valence pseudopotential calculations agree quite well, if a proper f basis set instead of f polarization functions is applied.     

Energy-adjusted pseudopotentials for the rare earth elements

Nonrelativistic and quasirelativistic energy-adjusted pseudopotentials and optimized (7s6p5d)/[5s4p3d]-GTO valence basis sets for use in molecular calculations for fixed f-subconfigurations of the rare earth elements, La through Lu, have been generated. Atomic excitation and ionization energies from numerical HF, as well as SCF pseudopotential calculations using the derived basis sets, differ by less than 0.1 eV from numerical HF all-electron results. Corresponding values obtained from CI(SD), CEPA-1, as well as density functional calculations using the quasirelativistic pseudopotentials, are in reasonable agreement with experimental data.     

Quasirelativistic valence ab initio calculation of the potential-energy curves for Cd–rare gas atom pairs

The results of large-scale valence ab initio calculations of the potential-energy curves for the ground states and several excited states of Cd–rare gas (RG) van der Waals molecules are reported. In the calculations, Cd²⁰⁺ and RG⁸⁺ cores are simulated by energy-consistent pseudopotentials, which also account for scalar-relativistic effects and spin-orbit interaction within the valence shell. The potential energies of the Cd–RG species in the ΛS coupling scheme have been evaluated by means of ab initio complete-active-space multiconfiguration self-consistent-field (CASSCF)/CAS multireference second-order perturbation theory (CASPT2) calculations with a total 28 valence electrons, but the spin-orbit matrix has been computed in a reduced configuration interaction space restricted to the CASSCF level. Finally, the Ω potential curves are obtained by diagonalization of the modified spin-orbit matrix (its diagonal elements before diagonalization substituted by the corresponding CASPT2 eigenenergies). The calculated potential curves, especially the spectroscopic parameters derived for the ground states and several excited states of the Cd–RG species are presented and discussed in the context of available experimental data. The theoretical results exhibit very good agreement with experiment. Received: 20 April 2000 / Accepted: 1 September 2000 / Published online: 21 December 2000     

Combined pseudopotential and density functional study of Bis-eta6-benzene d and f element complexes

The electronic and molecular structures of selected zerovalent d and f metal bis-eta6-benzene sandwich complexes M(C6H6)2 (M = Sc, Ti, V, Cr, Y, Zr, Nb, Mo, La, Ce, Gd, Lu, Hf, Ta, Nb, Th) were investigated at the scalar-relativistic level using a combination of energy-consistent ab initio pseudopotentials and gradient-corrected density functionals. The results are compared to previous pseudopotential coupled-cluster and all-electron density functional results based on the zero-order regular approximate Hamiltonian of Chang, Pelissier, and Durand as well as to experimental data.     

An analysis of core effects on shape-consistent pseudopotentials

Large core (seven-valence electrons) shape-consistent averaged relativistic pseudopotentials (AREP) including core effects have been derived for the halogen series (Cl,Br,I,At). The influence of core effects on the spin-orbit splitting of the halogen and alkali atoms is clearly demonstrated within an all-electron four-component atomic reference calculation by means of a perturbation analysis. In particular, it is shown that AREPs extracted at the Dirac-Coulomb-Fock level, which already include spin-orbit polarization effects, give excellent results for atomic spectroscopy and equilibrium distances of halogen dimers. We also show that in our approach the core effects, included by configuration interaction using the numerical GRASP code, are transferred to the averaged orbital one-electron energy, defined in a perturbational way. This leads to a modification of the extracted AREPs by core effects, which is illustrated by calculations of the first atomic excited states using these AREPs. These results support the validity of including core effects directly in the AREPs extracted in a shape-consistent scheme. The transferability to the atomic excited states as well as to the molecular case is also verified.     

Ground-state properties and static dipole polarizabilities of the alkali dimers from K2 n to Fr2 n(n=0,+1) from scalar relativistic pseudopotential coupled cluster and density functional studies

The newly adjusted energy-consistent nine-valence-electron pseudopotentials for K to Fr are used to calculate spectroscopic properties for the neutral and positively charged alkali dimers using coupled cluster and density functional theory. For the neutral dimers the static dipole polarizability was calculated. The coupled cluster results are all in excellent agreement with experimental values. The density functionals used can give quite different spectroscopic properties especially for the dipole polarizability, with the Perdew-Wang PW91 functional performing best.     

Quasirelativistic energy-consistent 5f-in-core pseudopotentials for trivalent actinide elements

Quasirelativistic energy-consistent 5f-in-core pseudopotentials modelling trivalent actinides, corresponding to a near-integral 5f ⁿ occupation (n = 0–14 for Ac–Lr), have been generated. Energy-optimized (6s5p4d), (7s6p5d), and (8s7p6d) primitive valence basis sets contracted to polarized double to quadruple zeta quality as well as 2f1g correlation functions have been derived. Corresponding smaller basis sets (4s4p3d), (5s5p4d), and (6s6p5d) suitable for calculations on actinide(III) ions in crystalline solids form subsets of these basis sets designed for calculations on neutral molecules. Results of Hartree–Fock test calculations for actinide(III) monohydrates and actinide trifluorides show a satisfactory agreement with corresponding calculations using 5f-in-valence pseudopotentials. Even in the beginning of the actinide series, where the 5f shell is relatively diffuse, only quite acceptable small deviations occur as long as the 5f-shell does not participate significantly in covalent bonding. Electronic Supplementary Material Supplementary material is available in the online version of this article at and is accessible for authorized users.