On the calculation of second-order magnetic properties using subsystem   approaches in the relativistic framework
  We report an implementation of the nuclear magnetic resonance (NMR) shielding ($\sigma$), isotope-independent indirect spin-spin coupling ($K$) and the magnetizability ($\xi$) tensors in the frozen density embedding (FDE) scheme using the four-component (4c) relativistic Dirac--Coulomb (DC) Hamiltonian and the non-collinear spin density functional theory (SDFT). The formalism takes into account the magnetic balance between the large and the small components of molecular spinors and assures the gauge-origin independence of NMR shielding and magnetizability results. This implementation has been applied to hydrogen-bonded HXH$\cdots$OH$_2$ complexes (X = Se, Te, Po) and compared with the supermolecular calculations and with the approach based on the integration of the magnetically induced current density vector. A comparison with the approximate Zeroth-Order Regular Approximation (ZORA) Hamiltonian indicates non-negligible differences in $\sigma$ and $K$ in the HPoH$\cdots$OH$_2$ complex, and calls for a thourough comparison of ZORA and DC in the description of environment effects on NMR parameters for molecular systems with heavy elements. 
Malgorzata Olejniczak and Radovan Bast and Andre Severo Pereira Gomes
physics.chem-ph