Total: 1
We present an implementation of spin-orbit coupling (SOC) for the computation of nuclear magnetic resonance (NMR) chemical shielding tensors within linear response theory. Our implementation in the Vienna {\it Ab initio} Simulation Package (VASP) is tailored to solid-state systems by employing periodic boundary conditions and the gauge-including projector augmented waves (GIPAW) 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.