Recovering four-component solutions by the inverse transformation of the infinite-order two-component wave functions

Abstract

The two-component Hamiltonian of the infinite-order two-component (IOTC) theory is obtained by a unitary block-diagonalizing transformation of the Dirac-Hamiltonian. Once the IOTC spin orbitals are calculated, they can be back transformed into four-component solutions. The transformed four component solutions are then used to evaluate different moments of the electron density distribution. This formally exact method may, however, suffer from certain approximations involved in its numerical implementation. As shown by the present study, with sufficiently large basis set of Gaussian functions, the Dirac values of these moments are fully recovered in spite of using the approximate identity resolution into eigenvectors of the p(2) operator.