Density-functional calculations of relativistic spin-orbit effects on nuclear magnetic shielding in paramagnetic molecules

Terms arising from the relativistic spin-orbit effect on both hyperfine and Zeeman interactions are introduced to density-functional theory calculation of nuclear magnetic shielding in paramagnetic molecules. The theory is a generalization of the former nonrelativistic formulation for doublet systems and is consistent to O(alpha4), the fourth power of the fine structure constant, for the spin-orbit terms. The new temperature-dependent terms arise from the deviation of the electronic g tensor from the free-electron g value as well as spin-orbit corrections to hyperfine coupling tensor A, the latter introduced in the present work. In particular, the new contributions include a redefined isotropic pseudocontact contribution that consists of effects due to both the g tensor and spin-orbit corrections to hyperfine coupling. The implementation of the spin-orbit terms makes use of all-electron atomic mean-field operators and/or spin-orbit pseudopotentials. Sample results are given for group-9 metallocenes and a nitroxide radical. The new O(alpha4) corrections are found significant for the metallocene systems while they obtain small values for the nitroxide radical. For the isotropic shifts, none of the three beyond-leading-order hyperfine contributions are negligible.