The estimated \(B_l^m\)  parameters predict the five Kramers' doublets given in Table REF .
All four transitions from the CEF ground state to the excited states occur for energy transfers below 100 meV, which is covered by our experiments.
However, only three flat INS modes are observed experimentally, because the energies of two excited states almost coincide and yield a single broad peak at 38 meV.
Knowing the composition of the CEF ground state allows us to determine the corresponding \(g\)  factor anisotropy [1]},
\(\begin{split}g_{z} = 2g_J |\langle \pm \omega _0|J_z|\pm \omega _0\rangle | = 2.78, \\g_{xy} = g_J |\langle \pm \omega _0|J_\pm |\mp \omega _0\rangle | =1.22,\end{split}\) 
