SO is not sufficient to describe the SO coupling within the spin-only ground multiplet 8 H 7/2 of Gd(III). Here, the crystal field levels calculated by H ) are largely overestimated with respect to the H BP SO results (∆E = 0.3, 0.5, 1.0 cm -1 ), leading to µ ≈ +350%.