where α, β, and γ are the Stevens constants for each Ln, 𝑂 ̂𝑘 𝑞 are the operator equivalents expressed as polynomials of the total angular momentum operators, <r k ˃ are expectation values of the radial factor r k , and 𝐴 𝑘 𝑞 are numerical parameters that depend on the nature of the ligand shell. Unlike pseudo-axial systems, where only the 𝐴 𝑘 𝑞 <r k ˃ terms with q = 0 (𝐴 2 0 <r 2 ˃, 𝐴 4 0 <r 4 ˃, 𝐴 6 0 <r 6 ˃) are different from zero, the C5 symmetry allows a term with q = 5(𝐴 6 5 <r 6 ˃) that mixes magnetic states. To experimentally determine 𝐴 𝑘 𝑞 r k products, the temperature-dependent direct current magnetic susceptibility measurement can be used. Although dipole interactions between molecules in a crystal are weak, they are sufficient to cause significant inhomogeneous broadening of EPR lines and give rise to decoherence of the quantum spin state; diminishing the efficacy of these systems as qubits.