“…From this Hamiltonian, the energy eigenvalues, ε n , and eigenvectors, |n , are numerically determined to calculate the magnetization components. Details of the calculations of magnetization and entropy were described elsewhere [9]. In short, for a given magnetic field direction, one can use relations (9) and (10) to obtain the magnetization components˛= (x, y, z) as well as the directional magnetic entropy of the magnetic system under investigation.…”