This paper analyzes how special relativity changes the equation for the many-body-induced current density starting from the Foldy–Wouthuysen diagonalization of the Dirac–Coulomb Hamiltonian. This current density differs from that obtained with the Gordon decomposition due to the presence of a spin-orbit coupling contribution not considered before for many-body molecular systems. This contribution diverges on atomic nuclei due to the nature of the point charges considered in the nonrelativistic

more »

... pproach, demonstrating that conventionally used nonrelativistic methods are not suitable for dealing with spin effects such as spin-orbit coupling or effects smaller than α2, with α the fine structure constant, and that a fully relativistic approach with a finite charge should be used. Despite the singularity, the spin-orbit coupling current becomes an important contribution to the total current in open-shell systems with high-spin multiplicity and a high atomic number in the nuclear proximity. On long ranges, this contribution is overcome by the Coulomb potential and the derived electric field which decays very quickly for small distances from nuclear charges. An evaluation of this spin-orbit current has been performed in the linear response approach at the HF/DFT level of theory.