NMR chemical shielding for solid-state systems using spin–orbit coupled ZORA GIPAW

We present an implementation of spin–orbit coupling (SOC) for the computation of nuclear magnetic resonance chemical shielding tensors within linear response theory. Our implementation in the Vienna Ab initio Simulation Package is tailored to solid-state systems by employing periodic boundary conditions and the gauge-including projector augmented waves approach. Relativistic effects are included on the level of the zeroth-order regular approximation (ZORA). We discuss the challenges posed by the PAW partial wave basis in describing SOC regarding chemical shielding tensors. Our method is in good agreement with existing local-basis ZORA implementations for a series of Sn, Hg, and Pb molecules and cluster approximations for crystalline systems.