CASSCF-SO calculations were performed with OpenMolcas 77 on XRD structures of paramagnetic 1-M. Basis sets from the ANO-RCC library were used with VTZP quality on the metal atom, VDZP quality on the coordinating atoms and VDZ quality on all other atoms, employing the second-order DKH transformation. Cholesky decomposition of the two-electron integrals with a threshold of 10 -8 was performed to save disk space and reduce computational demand. The molecular orbitals (MOs) were averaged in state-averaged CASSCF calculations with active spaces of (1,7) (1-Ce), (2,7) (1-Pr) and (3,7) (1-Nd and 1-U) averaging over seven doublets (1-Ce), 21 triplets and 28 singlets (1-Pr), and 35 quartets and 112 doublets (1-Nd and 1-U). The spin-free wavefunctions obtained from these CASSCF calculations were mixed by spin-orbit coupling in the RASSI module. The g-values, magnetization and magnetic susceptibility (isotropic value and tensor) were calculated using SINGLE_ANISO, 122 and the spin-orbit wavefunctions were decomposed into their crystal field wavefunctions, with the quantization axis defined by the g1 direction in the ground doublet. From the magnetic susceptibility tensor we calculate the pseudo-contact shift (𝛿 PCS para ) for the Si and where 𝜒̅ = 1 3