doi:10.1016/j.ccr.2006.01.003
Copyright © 2006 Elsevier B.V. All rights reserved.
Review
Relativistic energy-consistent ab initio pseudopotentials as tools for quantum chemical investigations of actinide systems
Xiaoyan Cao and Michael Dolg
, 
Institute for Theoretical Chemistry, University of Cologne, Greinstr. 4, D-50939 Cologne, Germany
Received 21 September 2005;
accepted 6 January 2006.
Available online 3 March 2006.
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Abstract
The method of relativistic energy-consistent ab initio pseudopotentials is briefly reviewed, paying special attention to the parametrization for actinide elements and the optimization of corresponding valence basis sets. Based on atomic frozen-core data it is argued that only a small-core approach is sufficiently accurate. Calibration calculations for atoms and diatomic molecules are briefly described and an overview is given over recent applications of the method to larger actinide systems. Shortcommings of the currently applied approach are discussed and future developments, towards higher accuracy as well as towards a simplified treatment of actinides in quantum chemical calculations, are outlined.
Keywords: Pseudopotentials; Actinides; Relativistic effects; Electronic structure; Calculation
Fig. 1. (Left) Radial densities of the 5f, 6d and 7s valence spinors of U in the 5f36d17s2 ground state configuration in comparison to the radial densities of the U32+, U14+ and U6+ cores from average-level multi-configuration Dirac–Hartree–Fock calculations [18]. (Right) Non-relativistic orbital energies (Hartree–Fock, HF) and relativistic spinor energies (Dirac–Hartree–Fock, DHF) for valence shells of U from average level calculations [18]. The horizontal dotted lines correspond to a “chemical” core-valence separation based on one-particle energies.
Fig. 2. Third (left bars) and fourth (right bars) ionization potentials of the actinides estimated by pseudopotential (PP) multi-reference averaged coupled-pair functional (MR-ACPF) calculations including spin–orbit corrections and extrapolation to the basis set limit (dotted bars) [19]. Relativistic contributions estimated from the difference of state-averaged multi-configuration Dirac–Hartree–Fock (Dirac–Coulomb–Breit Hamiltonian) and Hartree–Fock calculations (striped bars) [18] and electron correlation contribution estimated from PP MR-ACPF correlation energies extrapolated to the basis set limit (filled bars) [19]. Crosses with error bars denote the experimentally measured values for Th and the semiempirical estimates for U [17].
Fig. 3. Total relativistic contributions estimated from the difference of state-averaged Dirac–Hartree–Fock (Dirac–Coulomb–Breit Hamiltonian) and Hartree–Fock calculations (striped bars) [18] and spin–orbit contributions estimated from pseudopotential calculations with and without spin–orbit operator (filled bars) [19] to the third (left bars) and fourth (right bars) ionization potentials of the actinides.
Fig. 4. Schematic representation of the effective core potential (ECP) approximation in a coordinate system of three axes determining the quality of the many-electron Hamiltonian, the one- and the many-electron basis sets.
Fig. 5. Errors in spin–orbit corrected coupled-cluster singles, double and perturbative triples results for IP4 of Ce and Th. The experimental values are 36.76 and 28.65 eV for Ce and Th, respectively. The errors of uncorrelated calculations are denoted by a circle and square for Ce and Th, respectively. Relativistic energy-consistent small-core ab initio pseudopotentials (32 valence electrons) and large uncontracted basis sets (16s15p12d10f8g8h8i for Ce and 14s13p10d8f6g6h6i for Th) have been applied. The errors for basis sets including functions up to angular quantum number l are linear in 1/l3 for l>3 and can be extrapolated to the basis set limit.
Table 1.
Electronic ground configurations and states of the actinides Ann+ (n=0–4) [17]

f, d and s denote the 5f, 6d and 7s shells. Calculations find a Pa2+
ground state [19].
Table 2.
Relative Dirac–Hartree–Fock/Dirac–Coulomb (DHF/DC) energies (eV) of the 2J+ 1-weighted average of all J levels belonging to a non-relativistic configuration [18] with respect to the value for the U [Rn] 5f3 6d1 7s2 ground state configuration. Only subconfigurations outside the Rn core are listed (conf)

For each f occupation number only the energetically lowest configuration for each ionization level is listed. Contributions of the Breit interaction evaluated in first-order perturbation theory are also given (+B). The frozen-core errors (eV) in the relative energies are given for 32, 14 and 6 valence electron systems. The frozen core was taken from the neutral atom in the ground state configuration. Frozen cores: U Q=6: 1s–5d, 6s, 6p; Q=14: 1s–5d; Q=32: 1s–4d.
Table 3.
Bond lengths Re (Å), vibrational constants ωe (cm−1) and binding energies De (eV) of ThO in the
ground state from energy-consistent small-core pseudopotential calculations [3] in comparison to experimental data

The calculated values are given without/with Boys–Bernardi counter-poise correction of the basis set superposition error.
a Theoretical
De values have been corrected for molecular (0.03 eV) and atomic (Th 0.38, O 0.01 eV) spin–orbit energy lowerings; the zero-point energy (0.06 eV) was subtracted to obtain
D0 [2].
b PP with 60 core-electrons
[2], basis sets Th (14s13p10d8f6g)/[6s6p5d4f3g] ANO
[3], O aug-cc-pVQZ (spdfg); Th 5s,5p and O 1s frozen in CCSD(T).
Table 4.
Overview over recent calibration studies and applications of energy-consistent relativistic small-core actinide pseudopotentials in chronological order
a Z=2,2':6' 2''-terpyridine, 2,6-bis(5,6-dimethyl-1,2,4-triazin-3-yl)pyridine;
Ng=Ar, Kr, Xe.
b For an explanation of the (standard) acronyms characterizing the computational method the reader is referred to the original articles; short hand notations: B3LYP/
PW91=B3LYP, B3PW91; B(3)LYP/
PW91=BLYP, B3LYP, BPW91, B3PW91;
PBE(0)=PBE, PBE0; etc.
c S: structure optimisation; B: binding (and/or reaction) energy; V: vibrational frequencies; I: infrared intensities; E: electronic excitation energies; N: nuclear magnetic resonance data; T: electron transfer rate constants.
Table 5.
Comparison of the U–F bond length Re (Å), the infrared frequencies νi (cm−1) and the bond dissociation energy De (kcal/mol) of UF6 from various pseudopotential density functional theory calculations to experimental data

The intensities of the main peaks are listed in parentheses. For explanation of the (standard) acronyms characterizing the computational method as well as references for experimental values see Refs. [49] and [63].
Table 6.
Bond length (Å) and binding energies (eV) of actinide(III) monohydrate complexes An(H2 O)3+from 5f-in-core (LC) and 5f-in-valence (SC) pseudopotential (PP) Hartree–Fock [70] and all-electron (AE) Dirac–Hartree–Fock [71] calculations

The PP results are given for small/large basis sets. The 5f Mulliken population is also listed. Basis sets: AE An (30s25p19d13f2g), O (9s4p1d), H (4s1p) for large components; PP small basis sets: LC An (8s7p6d2f); PP SC An (14s13p10d8f)/[6s6p5d4f], O and H cc-pVDZ; PP large basis sets: LC An (8s7p6d2f1g); PP SC An (14s13p10d8f6g)/[6s6p5d4f3g], O and H aug-cc-pVQZ.