Conservation-laws–preserving algorithms for spin dynamics simulations
Institute for Condensed Matter Physics - 1 Svientsitsky Str., UA-79011 Lviv,
2 Institute for Theoretical Physics, University of Linz - A-4040 Linz, Austria
Accepted: 20 October 2000
We propose new algorithms for numerical integration of the equations of motion for classical spin systems with fixed spatial site positions. The algorithms are derived on the basis of a mid-point scheme in conjunction with the multiple time staging propagation. Contrary to existing predictor-corrector and decomposition approaches, the algorithms introduced preserve all the integrals of motion inherent in the basic equations. As is demonstrated for a lattice ferromagnet model, the present approach appears to be more efficient even over the recently developed decomposition method.
PACS: 02.60.Cb – Numerical simulation; solution of equations / 75.10.Hk – Classical spin models / 75.40.Gb – Dynamic properties (dynamic susceptibility, spin waves, spin diffusion, dynamic scaling, etc.)
© EDP Sciences, 2000