Kramers restricted Hartree-Fock method

Even when a REP is explicitly considered, the total wave function of a ground state closed-shell molecule is approximated by one Slater determinant in the HF method as 

The two-component Hamiltonian of Eq. (1) is invariant under the time reversal operation. In the one-electron case and for a special choice of phases, the time reversal operator T is given by 

As in the nonrelativistic case, all one-electron molecular spinors are expressed in terms of orbital basis and spin functions 

Using the total electronic energy for the KRHF theory in Eq. (14), the gradient of the total energy with respect to nuclear displacements is expressed as 

---

In Table 6.3, the values of De for RfCU are compared with those obtained within various approximations using relativistic effective core potentials (RECP) Kramers-restricted Hartree-Fock (KRHF) (Han et al 1999), averaged RECP including second-order M0ller-Plesset perturbation theory (AREP-MP2) for the correlation part (Han et al. 1999), RECP coupled-cluster single double (triple) [CCSD(T)] excitations (Han et al. 1999), and a Dirac-Fock-Breit (DFB) method (Malli and Styszynski 1998). The AREP-MP2 calculation of De gives 20.4 eV, while the RECP-CCSD(T) method with correlation leads to 18.8 eV. Our value of De of 19.5 eV is just between these calculated values. 

The incorporation of electron correlation effects in a relativistic framework is considered. Three post Hartree-Fock methods are outlined after an introduction that defines the second quantized Dirac-Coulomb-Breit Hamiltonian in the no-pair approximation. Aspects that are considered are the approximations possible within the 4-component framework and the relation of these to other relativistic methods. The possibility of employing Kramers restricted algorithms in the Configuration Interaction and the Coupled Cluster methods are discussed to provide a link to non-relativistic methods and implementations thereof. It is shown how molecular symmetry can be used to make computations more efficient. 

The alternative to the development of new algorithms to handle relativistic Hamiltonians is to search for a way to extend non-relativistic algorithms such that they can handle the additional couplings. Since most implementations are based on a restricted Hartree-Fock scheme the first step is to mimic the spin-restricted excitation operators used in the non-relativistic methods by Kramers restricted excitation operators. This can be done by employing the so-called X-operator formalism [37]. 

---