There have been several successful applications of the Dirac–Hartree–Fock (DHF) equations to the calculation of numerical electronic wave functions for diatomic molecules (Laaksonen and Grant 1984a, 1984b, Sundholm 1988, 1994, Kullie et al. 1999). However, the use of numerical techniques in relativistic molecular calculations encounters the same difficulties as in the nonrelativistic case, and to proceed to general applications beyond simple diatomic and linear molecules it is necessary to resort to an analytic approximation using a basis set expansion of the wave function. The techniques for such calculations may to a large extent be based on the methods developed for nonrelativistic calculations, but it turns out that the transfer of these methods to the relativistic case requires special considerations. These considerations, as well as the development of the finite basis versions of both the Dirac and DHF equations, form the subject of the present chapter. In particular, in the early days of relativistic quantum chemistry, attempts to solve the DHF equations in a basis set expansion sometimes led to unexpected results. One of the problems was that some calculations did not tend to the correct nonrelativistic limit. Subsequent investigations revealed that this was caused by inconsistencies in the choice of basis set for the small-component space, and some basic principles of basisset selection for relativistic calculations were established. The variational stability of the DHF equations in a finite basis has also been a subject of debate. As we show in this chapter, it is possible to establish lower variational bounds, thus ensuring that the iterative solution of the DHF equations does not collapse. There are two basically different strategies that may be followed when developing a finite basis formulation for relativistic molecular calculations. One possibility is to expand the large and small components of the 4-spinor in a basis of 2-spinors. The alternative is to expand each of the scalar components of the 4-spinor in a scalar basis. Both approaches have their advantages and disadvantages, though the latter approach is obviously the easier one for adapting nonrelativistic methods, which work in real scalar arithmetic.
