Energy balance and energy correction in dynamics of classical spin systems

Energy-correction method is proposed as an addition to mainstream integrators for equations of motion of systems of classical spins. This solves the problem of non-conservation of energy in long computations and makes mainstream integrators competitive with symplectic integrators for spin systems that for different-site interactions conserve the energy explicitly. The proposed method is promising for spin systems with single-site interactions for which symplectic integrators do not conserve energy and thus have no edge against mainstream integrators. From the energy balance in the spin system with a phenomenological damping and Langevin fields, a formula for the dynamical spin temperature in the presence of single-site anisotropy is obtained.