Gaussian functions, spin-orbit operators

The quasirelativistic (QR) PP of Hay and Wadt [61] use two-component wave functions, but the Hamiltonian includes the Darwin and mass-velocity terms and omits the spin-orbit effects. The latter are then included via the perturbation operator after the wave functions have been obtained. The advantage of die method is the possibility to calculate quite economically rather large systems. The method is implemented in the commercial system Gaussian 98 It has extensively been applied to calculations of transition-element and actinide systems [62],... 

---