Calculate the energy change (Δe) using the unscaled distances rij, spins Si, anisotropy constants Ki and Jij/n as the exchange integral parameter, taking exp(−Δe/kBT) as the acceptance probability.