Four-component relativistic theory for nuclear magnetic shielding: Magnetically balanced gauge-including atomic orbitals

It is recognized only recently that the incorporation of the magnetic balance condition is absolutely essential for four-component relativistic theories of magnetic properties. Another important issue to be handled is the so-called gauge problem in calculations of, e.g., molecular magnetic shielding tensors with finite bases. It is shown here that the magnetic balance can be adapted to distributed gauge origins, leading to, e.g., magnetically balanced gauge-including atomic orbitals (MB-GIAOs) in which each magnetically balanced atomic orbital has its own local gauge origin placed on its center. Such a MB-GIAO scheme can be combined with any level of theory for electron correlation. The first implementation is done here at the coupled-perturbed Dirac-Kohn-Sham level. The calculated molecular magnetic shielding tensors are not only independent of the choice of gauge origin but also converge rapidly to the basis set limit. Close inspections reveal that (zeroth order) negative energy states are only important for the expansion of first order electronic core orbitals. Their contributions to the paramagnetism are therefore transferable from atoms to molecule and are essentially canceled out for chemical shifts. This allows for simplifications of the coupled-perturbed equations.