Numerical study of the iterative solution of the one-electron Dirac equation based on 'direct perturbation theory'

The one-electron Dirac equation is solved in an iterative manner starting with the solution of the Schrodinger equation. The method is applied in a basis of atom-centred Gaussian-type functions to the ground state of selected hydrogen-like ions up to Eka Pt109+ and the heavy quasi-molecules Th-2(179+), NiPb(109+)and Th-3(269+) (in D-infinity h and D-3h symmetry). An overall 8-figure accuracy in the absolute relativistic energies is achieved. The iterative procedure converges better than linearly for light systems and linearly for systems containing nuclear charges greater than Z approximate to 40.