Basis functions large-component

The molecular orbitals in the nonrelativistic and one-component calculations and the large component in the Dirac-Fock functions were spanned in the Cd s Ap9d)l[9slp6d basis of [63] and the H (5s 2p)/[35 l/>] set [61]. Contraction coefficients were taken from corresponding atomic SCF calculations. The basis for the small components in the Dirac-Fock calculations is derived by the MOLFDIR program from the large-component basis. The basis set superposition error is corrected by the counterpoise method [64]. The Breit interaction was found to have a very small effect and is therefore not included in the results. 

The inclusion of relativistic effects is essential in quantum chemical studies of molecules containing heavy elements. A full relativistic calculation, i.e. based upon Quantum Electro Dynamics, is only feasible for the smallest systems. In the SCF approximation it involves the solution of the Dirac Fock equation. Due to the four component complex wave functions and the large number of basis functions needed to describe the small component Dirac spinors, these computations are much more demanding than the corresponding non-relativistic ones. This limits Dirac Fock calculations, which can be performed using e.g. the MOLFDIR package [1], to small molecular systems, UFe being a typical example, see e.g. [2]. 

The IAM model further assumes the atoms in a crystal to be neutral. This assumption is contradicted by the fact that molecules have dipole and higher electrostatic moments, which can indeed be derived from the X-ray diffraction intensities, as further discussed in chapter 7. The molecular dipole moment results, in part, from the nonspherical distribution of the atomic densities, but a large component is due to charge transfer between atoms of different electronegativity. A population analysis of an extended basis-set SCF wave function of HF, for example, gives a net charge q of +0.4 electron units (e) on the H atom in HF for CH4 the value is +0.12 e (Szabo and Ostlund 1989). 

The l/REP(r), U ARKP(r), and terms At/f EP(r) in 11s0 of Eqs. (23), (31), and (34) or Eq. (6), respectively, are derived in the form of numerical functions consistent with the large components of Dirac spinors as calculated using the Dirac-Fock program of Desclaux (27). These operators have been used in their numerical form in applications to diatomic systems where basis sets of Slater-type functions are employed (39,42,43). It is often more convenient to represent the operators as expansions in exponential or Gaussian functions (32). The general form of an expansion involving M terms is... 

---