Abstract
Two-dimensional, fully numerical, finite-difference approach to the second-order Dirac equation for one-electron diatomics is revisited. Instead of using the transformation of variables to get rid of the singularities of σ-type spinors the Gaussian nuclear charge distribution model is employed. The approach is tested on the several lowest σ, π, δ and φ states of H, Kr and Th
systems with the reference data of 12-digit accuracy provided by the GRASP
and QRHF atomic codes.
The method is also tested on H, Kr
and Th
systems. The energies of
and
states of H
are converged to 9 significant figures and the relativistic correction to 6. In the case of the Th
system the calculations for the
state were performed for a range of internuclear separations
a.u. and compared with the eXact 2-Component (X2C) and 4-component results. The agreement to within 1 milihartree or better is observed for all the separations greater than 0.2 a.u. Alas, the
orbital energy of the system with
a.u., could only be calculated with 0.1% accuracy. This problem also seems to haunt Sundholm's implementation of the method and we indicate the reasons behind the low accuracy at extremely short internuclear distances.
GRAPHICAL ABSTRACT
![](/cms/asset/9716e2e3-b721-470c-99c7-b6ad00deb8c6/tmph_a_2092563_uf0001_ob.jpg)
Acknowledgments
We would like to thank Jacek Karwowski for his comments on the draft version of the manuscript. We are also very much indebted to the referee for the host of valuable remarks that helped improve the manuscript.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Correction Statement
This article has been corrected with minor changes. These changes do not impact the academic content of the article.
Notes
1 The following values of atomic masses (in amu) were used in the DIRAC2D and Molcas programs: H – 1.0078250350, Kr – 83.911507, Rn – 222.017571, Th – 232.0380508 (in the DIRAC program the mass of the hydrogen nucleus is set to 1.007825).
2 Note that the term is missing in Equation 11a,b of the paper by Sundholm et al. [Citation13].
3 The calculations were realised via the complete active space self-consistent field method (CASSCF) [Citation49,Citation50]
4 According to Bağcı and Hoggan [Citation51] the relativistic energy of 1s state of the hydrogen atom with the point nucleus can be calculated via the Rayleigh-Ritz method as
hartree and – according to this paper – is equal to the binding energy that can be obtained from the analytical formula (as given by Eg. 62 of that paper), when the velocity of light equal to
a.u. is used. However, when the formula is used with this very velocity of light it leads to a slightly different value of the energy, namely,
hartree. This value can be obtained by means of the bc utility implementing an arbitrary precision calculator.