Consequently, (6) ΔEi′=n3Δei where (7) Δei=−n−1∑jJijΔSi⋅Sj−Ki(ΔSi^⋅ni^)2−gμBμ0Hi⋅ΔSi+D∑jΔSi⋅Sj−3(ΔSi⋅eij)(Sj⋅eij)rij3 is the energy change of, let us say, a reduced “e-system”, similar to the original one but with the exchange energy divided by the n parameter. It is interesting that the energy change of the super-cell ΔE’ is equivalent to the energy change of the single spin in n3 e-systems simultaneously. Through this conclusion, the re-scaled system can be simulated with the thermodynamic balance determined for the reduced one, which is more effective, taking into account the fact that exp(−ΔE’/kBT) << exp(−Δe/kBT). It should be emphasized that the most important assumption is the coherent rotation of all spins in the super-cells. This causes a restriction in the temperature range, which has to be significantly lower than the Curie point. Therefore, the value kBT can be considered a kind of “annealing temperature”, enabling effective relaxation of the system and leading to minimization of its free energy.