“…The appropriate phenomenological spin Hamiltonian, taking into account the Zeeman perturbation in axially distorted surroundings, is given by Equation (2): (2) in which Ŝ x , Ŝ y and Ŝ z are the spin operators, H x , H y and H z are the components of the magnetic field along the coordinate axes x, y and z, g s and g ⊥ are the parallel and perpendicular components of the g tensor, μ B is the Bohr magneton and 2D is the energy gap between ±3/2 and ±1/2 Kramer doublets in the absence of a magnetic field. The experimental data, at temperatures above that giving the maximum susceptibility value, were analysed by means of Equations (3a, 3b, 3c) derived from the Hamiltonian shown in Equation (2): (3a) (3b) (3c) in which the term θ has been added to account for the magnetic coupling between the mononuclear units and g is supposed to be isotropic to avoid overparametrization, while N is the Avogadro number and k the Boltzmann constant.…”