Theoretical modeling of the nuclear magnetic resonance (NMR) chemical shifts in paramagnetic transition metal complexes has undergone considerable progress in the last fifty years. The foundation was laid early in the 1970s by McGarvey to model the paramagnetic chemical shifts in terms of the EPR spin Hamiltonian parameters which to date are desired by experimentalists to interpret the paramagnetic chemical shifts . The most recent and widely used model equations based on the EPR spin Hamiltonian parameters of the paramagnetic center, as summarized in the review , are upgrades of Kurland and McGarvey's work, with the primary emphasis being on evaluating the tensorial parameters using QC methods . The model equations require system-specific modeling and evaluation of the parameters from the QC methods. However there is another way to probe the paramagnetic center that is used routinely by other groups within the magnetic resonance community, particularly those con-cerned with paramagnetic biomolecules and solidstate NMR of materials . In this later case, the QC evaluation of the chemical shielding tensor (σ Ki j ) of nucleus K in a paramagnetic complex follows eq. 1,