where 𝐽 &' is the exchange interaction between the atomic sites 𝑖 and 𝑗 with corresponding spin direction 𝐒; 𝑘 ) is the uniaxial anisotropy constant along an easy axis 𝐞. The terms within Equation ( ) indicate two main contributions in a magnetic system: the exchange interaction between pairs of local spins and the magnetic uniaxial anisotropy, respectively. The value of 𝐽 &' and 𝑘 ) for fcc Ni and hcp Co are shown in Table . A Landau-Lifshitz-Gilbert-Heun integrator is used to compute the normalized mean magnetization 𝑀 as a function of temperature.