Coulomb operator Hartree-Fock calculations

The expression for the lowest order contribution to the parity violating potential within the Dirac Hartree-Fock framework is identical to that within the relativistically parameterised extended Hiickel approach in eq. (146). The difference is, however, that in DHF typically atomic basis sets with fixed radial functions are employed (see [161]) and that the molecular orbital coefficients are obtained in a self-consistent Dirac Hartree-Fock procedure. Computations of parity violating potentials along these lines have occasionally been called fully relativistic, although this term is rather unfortunate. In the four-component Dirac Hartree-Fock calculations by Quiney, Skaane and Grant [155] as well as in those by Schwerdtfeger, Laerdahl and coworkers [65,156,162,163] the Dirac-Coulomb operator has been employed, which is for systems with n electrons given by... 

By combining the separation into Coulomb and exchange terms with the multipole technique of Section 9.14, it is possible to arrive at an operation count that scales linearly with the size of the system, allowing us to cany out direct Hartree-Fock calculations for very large systems. However,... 

In many calculations beyond the Hartree-Fock level a first step is the transformation of at least some integrals. For the simplest such calculation, second-order perturbation theory, integrals with two indices transformed into the occupied MO basis axe required. Such integrals appear in many situations, including the MO basis formulation of coupled-perturbed Hartree-Fock theory. We can represent the first phase of this transformation as obtaining Coulomb and exchange operators ... 

One guesses at an initial set of wave functions, , and constructs the Hartree-Fock Hamilton S which depends on the through the definitions of the Coulomb and exchange operators, (/ and One then calculates the new set of , and compares it (or the energy or the density matrix) to the input set (or to the energy or density matrix computed from the input set). This procedure is continued until the appropriate self-consistency is obtained. 

Shape-consistent pseudopotentials including spin-orbit operators based on Dirac-Hartree-Fock AE calculations using the Dirac-Coulomb Hamiltonian have been generated by Christiansen, Ermler and coworkers [161-170]. The potentials and corresponding valence basis sets are also available on the internet under http //www.clarkson.edu/ pac/reps.html. A similar, quite popular set for main group and transition elements based on scalar-relativistic Cowan-Griffin AE calculations was published by Hay and Wadt [171-175]. 

These are important points for any quantitative work, electron-electron interactions must be taken into account, and the theories underpinning computation of MOs do this at various levels of accuracy. Approaches such as Hartree-Fock or density functional theory adapt the Hamiltonian operator to include electron-electron terms in an averaged way so electrons see the Coulomb field of each other averaged over the calculated density associated with each MO (see the Further Reading section in this chapter). 

According to the approach of Beattie and Landsberg [31], the intrinsic Auger lifetime is calculated by perturbation method. The perturbation operator of Auger mechanism, i.e., of Coulomb interaction of two electrons, is calculated by subtracting Hartree-Fock single-electron Hamiltonian from the complete Hamiltonian of the system (and actually the simplest Hamiltonian that still sees the Coulomb interaction of the Auger process). It has the form of a screened Coulomb potential... 

---