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.     

Relativistic effects determined using the Douglas-Kroll contracted basis sets and correlation consistent basis sets with small-core relativistic pseudopotentials

The coupled cluster approximation with single, double, and quasiperturbative triple excitations [CCSD(T)] was used in combination with the Douglas-Kroll contracted correlation consistent basis sets [cc-pVnZ-DK, where n = D(2), T(3), Q(4), and 5] and small-core relativistic pseudopotentials (PP) with correlation consistent polarized valence basis sets (cc-pVnZ-PP and aug-cc-pVnZ-PP) to investigate the impact of scalar relativistic corrections on energetic and structural properties of small molecules containing third-row (Ga-Kr) atoms. These molecules were taken from the Gaussian-2 extended test set for third-row atoms. Atomization energies, ionization energies, electron affinities, and proton affinities for molecules in the test set were determined and compared with nonrelativistic results which were obtained in a recent study in which the standard and augmented correlation consistent basis sets were used in combination with CCSD(T). Several schemes were used to extrapolate the energies to the complete basis set limit.     

Auxiliary basis sets for density fitting-MP2 calculations: Nonrelativistic triple-zeta all-electron correlation consistent basis sets for the 3d elements Sc-Zn

Auxiliary basis sets for density fitting second-order Moller-Plesset perturbation theory (DF-MP2) have been optimized for use with the triple-zeta nonrelativistic all-electron correlation consistent orbital basis sets, cc-pVTZ-NR and aug-cc-pVTZ-NR, for the 3d elements Sc-Zn. The relative error in using these auxiliary basis sets is found to be around four orders of magnitude smaller than that from utilizing triple-zeta orbital basis sets rather than corresponding quadruple-zeta basis sets, in calculation of the correlation energy for a test set of 54 small to medium sized transition metal complexes.     

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.     

All-electron scalar relativistic basis sets for the 6p elements

New segmented all-electron relativistically contracted (SARC) basis sets have been developed for the elements ₈₁Tl–₈₆Rn, thus extending the SARC family of all-electron basis sets to include the 6p block. The SARC basis sets are separately contracted for the second-order Douglas–Kroll–Hess and the zeroth-order regular approximation scalar relativistic Hamiltonians. Their compact size and segmented construction are best suited to the requirements of routine density functional theory (DFT) applications. Evaluation of the basis sets is performed in terms of incompleteness and contraction errors, orbital properties, ionization energies, electron affinities, and atomic polarizabilities. From these atomic metrics and from computed basis set superposition errors for a series of homonuclear dimers, it is shown that the SARC basis sets achieve a good balance between accuracy and size for efficient all-electron scalar relativistic DFT applications.     

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.     

Molecular core-valence correlation effects involving the post-d elements Ga-Rn: benchmarks and new pseudopotential-based correlation consistent basis sets

Correlation consistent basis sets that are suitable for the correlation of the outer-core (n-1)spd electrons of the post-d elements Ga-Rn have been developed. These new sets, denoted by cc-pwCVXZ-PP (X=D,T,Q,5), are based on the previously reported cc-pVXZ-PP sets that were built in conjunction with accurate small-core relativistic pseudopotentials (PPs) and designed only for valence nsp correlation. These new basis sets have been utilized in benchmark coupled cluster calculations of the core-valence correlation effects on the dissociation energies and spectroscopic properties of several small molecules. As expected, the most important contribution is the correlation of the (n-1)d electrons. For example, in the case of the group 13 homonuclear diatomics (Ga(2),In(2),Tl(2)), this leads to a dissociation energy increase compared to a valence-only treatment from 1.5 to 3.2 kcal/mol, bond length shortenings from -0.076 to -0.125 A?, and harmonic frequency increases of 7-8?cm(-1). Even in the group 15 cases (As(2),Sb(2),Bi(2)), the analogous effects of (n-1)d electron correlation are certainly not insignificant, the largest values being +4.4?kcal/mol, -0.049 A?, and +9.6?cm(-1) for the effects on D(e), r(e), and ω(e), respectively. In general, the effects increase in magnitude down a group from 4p to 6p. Correlation of the outer-core (n-1)p electrons is about an order of magnitude less important than (n-1)d but larger than that of the (n-1)s. The effect of additional tight functions for Hartree-Fock and valence sp correlation was found to be surprisingly large, especially for the post-4d and post-5d elements. The pseudopotential results for the molecules containing post-3d elements are also compared to the analogous all-electron calculations employing the Douglas-Kroll-Hess Hamiltonian. The errors attributed to the PP approximation are found to be very small.     

Density functional theory and the correlation consistent basis sets: the tight d effect on HSO and HOS

The HSO and HOS isomers have been revisited using the DFT functionals, B3LYP, B3PW91, and PBE, in combination with tight d-augmented correlation consistent basis sets, cc-pV(x+d)Z and aug-cc-pV(x+d)Z for second-row atoms. Structures, vibrationally averaged structures, relative energies, harmonic and anharmonic frequencies, enthalpies of formation of HSO and HOS, and the barrier for the HSO/HOS isomerization have been determined. These results were compared with results from previous DFT and ab initio studies in which the standard correlation consistent basis sets were used. The relative energies of the two isomers converge more rapidly and smoothly with respect to increasing basis set size for the tight d-augmented sets than for the standard basis sets. Our best calculations, B3PW91/aug-cc-pV(5+d)Z, for the relative energy of the isomers are in excellent agreement with previous CCSD(T) results given by Wilson and Dunning.     

Polarization consistent basis sets. V. The elements Si-Cl

Polarization consistent basis sets, optimized for density functional calculations, are proposed for the elements Si-Cl. Their performance for atomization energies, equilibrium geometries, harmonic vibrational frequencies, and associated infrared intensities is compared with other commonly used basis sets. Atomization energies can be predicted to within 0.01 kJ/mol per atom of the basis set limit by extrapolation of the pc-2, -3, and -4 results. Equilibrium bond distances and harmonic vibrational frequencies can be calculated to within 10(-5) A and 0.5 cm(-1), respectively, of the basis set limit. The pc-n basis sets are shown to give comparable or better accuracy than other alternatives, while containing fewer or equal number of primitive basis functions.     

Valence basis sets for lanthanide 4f-in-core pseudopotentials adapted for crystal orbital ab initio calculations

Crystal orbital adapted Gaussian (4s4p3d), (5s5p4d) and (6s6p5d) valence primitive basis sets have been derived for calculating periodic bulk materials containing trivalent lanthanide ions modeled with relativistic energy-consistent 4f-in-core lanthanide pseudopotentials of the Stuttgart-Koeln variety. The calibration calculations of crystalline A-type Pm₂O₃ using different segmented contraction schemes (4s4p3d)/[2s2p2d], (4s4p3d)/[3s3p2d], (5s5p4d)/[2s2p2d], (5s5p4d)/[3s3p3d], (5s5p4d)/[4s4p3d], (6s6p5d)/[2s2p2d], (6s6p5d)/[3s3p3d] and (6s6p5d)/[4s4p4d] are discussed at both Hartree–Fock (HF) and density functional theory (DFT) levels for the investigation of basis set size effects. Applications to the geometry optimization of A-type Ln₂O₃ (Ln = La-Pm) show a satisfactory agreement with experimental data using the lanthanide valence basis sets (6s6p5d)/[4s4p4d] and the standard set 6-311G* for oxygen. The corresponding augmented sets (8s7p6d)/[6s5p5d] with additional diffuse functions for describing neutral lanthanide atoms were applied to calculate atomic energies of free lanthanide atoms for the evaluation of cohesive energies for A-Ln₂O₃ within both conventional Kohn-Sham DFT and the a posteriori-HF correlation DFT schemes.     

Comments on relativistic basis sets

 Adding the tight and diffuse Gaussian-type functions (GTFs), Faegri's variationally determined double-zeta-quality basis sets for molecular relativistic calculations are examined. An example atom is Cm. When the tight s-type GTF is added the total energy increases, whereas when diffuse GTFs are added the total energy decreases. The reasons for these findings are clarified. It is also pointed out that not only the Faegri's sets but also other variationally determined basis sets would show similar behavior so far as the expansion terms are not sufficient. Received: 22 July 2002 / Accepted: 21 October 2002 / Published online: 31 January 2003 Correspondence to: H. Tatewaki e-mail: htatewaki@nsc.nagoya.cu.ac.jp     

Systematically convergent correlation consistent basis sets for molecular core-valence correlation effects: the third-row atoms gallium through krypton

The family of correlation consistent polarized valence basis sets has been extended in order to account for core-core and core-valence correlation effects within the third-row, main group atoms gallium through krypton. Construction of the basis sets is similar to that of the atoms boron through argon, where either the difference between core-correlated and valence-only correlation energies were calculated via configuration interaction (CISD) computations on the ground electronic states of the atoms (named cc-pCVnZ) or the sets were optimized with respect to the core-valence correlation energy and a small weight of core-core correlation energy (cc-pwCVnZ). Due to the correlation of 3d orbitals, added shells of higher angular momentum exponents compared to the valence sets are necessary to describe the core region. The pattern of added core-correlating functions is (1s1p1d1f) for double-zeta, (2s2p2d2f1g) for triple-zeta, (3s3p3d3f2g1h) for quadruple-zeta, and (4s4p4d4f3g2h1i) for quintuple-zeta. Atomic and molecular results show good convergence to the CBS limit, with the cc-pwCVnZ sets showing improved convergence compared to the cc-pCVnZ ones for molecular core-valence correlation effects. After testing the basis sets on the homonuclear diatomics Ga2-Kr2 with coupled cluster wave functions, it is concluded that a treatment of core-valence correlation effects is essential for high-accuracy ab initio investigations of third-row-containing molecules. Though the basis sets are optimal for 3s3p3d correlation, preliminary atomic and molecular results show the basis sets to be efficient with respect to 3d-only correlation, and these potentially could be used with 3d-only correlation for more qualitative studies on larger species.     

A pseudopotential-based composite method: the relativistic pseudopotential correlation consistent composite approach for molecules containing 4d transition metals (Y-Cd)

The correlation consistent composite approach (ccCA) has proven to be an effective first-principles-based composite approach for main group and first-row transition metal species. By combining relativistic pseudopotentials and ccCA, accurate energetic and thermodynamic data for heavier elements, including transition metals, is obtainable. Relativistic pseudopotential ccCA (rp-ccCA) was formulated and tested on 25 molecules from the G3∕05 set that contain 4p elements (Ga-Kr). A 32.5% time savings was obtained using rp-ccCA, relative to ccCA employing all-electron basis sets. When implementing rp-ccCA to compute dissociation energies and enthalpies of formation for molecules from the 4p block, rp-ccCA results in a mean absolute deviation of 0.89 kcal?mol(-1) from experimental data. rp-ccCA was also applied to a set of 30 4d transition metal-containing molecules, ranging from diatomics to Mo(CO)(6), and enthalpies of formation for these species were obtained with a mean absolute deviation of 2.89 kcal mol(-1) in comparison to experimental data. Based on quality of the experimentally available enthalpies of formation, where the average value of reported experimental error bars is 3.43 kcal mol(-1), rp-ccCA is within transition metal chemical accuracy for the 4d molecule set. rp-ccCA is a pseudopotential-based composite method for transition metals and is shown to yield accurate thermodynamic results for molecules containing heavy elements Ga-Kr and Y-Cd.     

Truncation of the correlation consistent basis sets: extension to third-row (Ga-Kr) molecules

The systematic reduction of the commonly used correlation consistent basis sets [cc-pVnZ where n=D(2), T(3), Q(4), and 5] as a means to reduce computational cost has been extended to hydrogen-containing third-row (Ga-Kr) molecules of the G2 test suite. Coupled cluster with singles, doubles, and quasiperturbative triple excitations [CCSD(T)] calculations were performed using both the full correlation consistent basis sets and a series of truncated basis sets in order to assess the impact of basis set reduction upon the structures and energies of the species. The impact that truncation of the basis sets for hydrogen has upon extrapolation of energies to the complete basis set limit also has been examined, and the cost savings that can be achieved are discussed. Overall, basis set reduction can be accomplished which preserves the systematic convergence behavior of the full correlation consistent basis sets.     

New relativistic ANO basis sets for transition metal atoms

New basis sets of the atomic natural orbital (ANO) type have been developed for the first, second, and third row transition metal atoms. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive and negative ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies, electron affinities, and excitation energies for all atoms and polarizabilities for spherically symmetric atoms. These calculations include spin-orbit coupling using a variation-perturbation approach. Computed ionization energies have an accuracy better than 0.2 eV in most cases. The accuracy of computed electron affinities is the same except in cases where the experimental values are smaller than 0.5 eV. Accurate results are obtained for the polarizabilities of atoms with spherical symmetry. Multiplet levels are presented for some of the third row transition metals.     

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.     

Accurate relativistic energy-consistent pseudopotentials for the superheavy elements 111 to 118 including quantum electrodynamic effects

Energy-consistent two-component semi-local pseudopotentials for the superheavy elements with atomic numbers 111-118 have been adjusted to fully relativistic multi-configuration Dirac-Hartree-Fock calculations based on the Dirac-Coulomb Hamiltonian, including perturbative corrections for the frequency-dependent Breit interaction in the Coulomb gauge and lowest-order quantum electrodynamic effects. The pseudopotential core includes 92 electrons corresponding to the configuration [Xe]4f(14)5d(10)5f(14). The parameters for the elements 111-118 were fitted by two-component multi-configuration Hartree-Fock calculations in the intermediate coupling scheme to the total energies of 267 up to 797 J levels arising from 31 up to 62 nonrelativistic configurations, including also anionic and highly ionized states, with mean absolute errors clearly below 0.02 eV for averages corresponding to nonrelativistic configurations. Primitive basis sets for one- and two-component pseudopotential calculations have been optimized for the ground and excited states and exhibit finite basis set errors with respect to the finite-difference Hartree-Fock limit below 0.01 and 0.02 eV, respectively. General contraction schemes have been applied to obtain valence basis sets of polarized valence double- to quadruple-zeta quality. Results of atomic test calculations in the intermediate coupling scheme at the Fock-space coupled-cluster level are in good agreement with those of corresponding fully relativistic all-electron calculations based on the Dirac-Coulomb-Breit Hamiltonian. The results demonstrate besides the well-known need of a relativistic treatment at the Dirac-Coulomb level also the necessity to include higher-order corrections for the superheavy elements.