A typical MC Metropolis algorithm refers to single spins (magnetic moments) placed in the so-called nodes with defined distances r. This well-known algorithm lies in a random walk over the system nodes and change of the spin direction with the acceptance probability equal to unity if it decreases the energy of the system. If the energy increases, the spin-flip can be accepted with the Metropolis probability exp(−ΔE/kBT) where ΔE = Enew − Eold is the difference in energy between the new and old spin configurations and kBT is the thermal energy.