Hartree-Fock theory relativistic methods

In this paper we present the first application of the ZORA (Zeroth Order Regular Approximation of the Dirac Fock equation) formalism in Ab Initio electronic structure calculations. The ZORA method, which has been tested previously in the context of Density Functional Theory, has been implemented in the GAMESS-UK package. As was shown earlier we can split off a scalar part from the two component ZORA Hamiltonian. In the present work only the one component part is considered. We introduce a separate internal basis to represent the extra matrix elements, needed for the ZORA corrections. This leads to different options for the computation of the Coulomb matrix in this internal basis. The performance of this Hamiltonian and the effect of the different Coulomb matrix alternatives is tested in calculations on the radon en xenon atoms and the AuH molecule. In the atomic cases we compare with numerical Dirac Fock and numerical ZORA methods and with non relativistic and full Dirac basis set calculations. It is shown that ZORA recovers the bulk of the relativistic effect and that ZORA and Dirac Fock perform equally well in medium size basis set calculations. For AuH we have calculated the equilibrium bond length with the non relativistic Hartree Fock and ZORA methods and compare with the Dirac Fock result and the experimental value. Again the ZORA and Dirac Fock errors are of the same order of magnitude. 

Accounting for relativistic effects in computational organotin studies becomes complicated, because Hartree-Fock (HF), density functional theory (DFT), and post-HF methods such as n-th order Mpller-Plesset perturbation (MPn), coupled cluster (CC), and quadratic configuration interaction (QCI) methods are non-relativistic. Relativistic effects can be incorporated in quantum chemical methods with Dirac-Hartree-Fock theory, which is based on the four-component Dirac equation. " Unformnately the four-component Flamiltonian in the all-electron relativistic Dirac-Fock method makes calculations time consuming, with calculations becoming 100 times more expensive. The four-component Dirac equation can be approximated by a two-component form, as seen in the Douglas-Kroll (DK) Hamiltonian or by the zero-order regular approximation To address the electron cor-... 

The UF6 molecule has also been studied extensively using a more elaborate method, namely configuration interaction, to assign the experimental photoelectron spectrum (de Jong and Nieuwpoort 1998). The qualitative analysis of chemical bonding exhibits that the U-F bond is more ionic in the relativistic framework (de Jong and Nieuwpoort 1998). The 6s orbital of uranium remains atom-like in the molecule due to relativistic contraction and does not contribute to chemical bonding, while it contributed in nonrelativistic Hartree-Fock theory. 

In molecular quantum chemistry Gaussian-fxmction-based computations, effective core potentials were originally derived from a reference calculation of a single atom within the nonrelativistic Hartree-Fock or relativistic Dirac-Fock (see Sect. 8.3) approximations, or from some method including electron correlations (Cl, for instance). A review of these methods, as well as a general theory of ECPs is provided in [480,481]. 

This expression excludes self-interaction. There have been a number of attempts to include into the Hartree-Fock equations the main terms of relativistic and correlation effects, however without great success, because the appropriate equations become much more complex. For a large variety of atoms and ions both these effects are fairly small. Therefore, they can be easily accounted for as corrections in the framework of first-order perturbation theory. Having in mind the constantly growing possibilities of computers, the Hartree-Fock self-consistent field method in various... 

---