1998/05/18 by M. Krech, Alex Bunker, D. P. Landau · 4 citations
Computer Science · Mathematics · Physics and Astronomy · #Matrix Theory and Algorithms #Numerical methods for differential equations #Physics of Superconductivity and Magnetism #cond-mat.mtrl-sci #cond-mat.stat-mech
paper · pdf · doi:10.1016/s0010-4655(98)00009-5
published as Computer Physics Communications 111 (1998) 1-13 · 9 pages RevTeX with 8 figures
arxiv created 1998/05/18 · openalex publication_date 1998/06/01 · arxiv updated 2009/11/30 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We have proposed new algorithms for the numerical integration of the equations of motion for classical spin systems. In close analogy to symplectic integrators for Hamiltonian equations of motion used in Molecular Dynamics these algorithms are based on the Suzuki-Trotter decomposition of exponential operators and unlike more commonly used algorithms exactly conserve spin length and, in special cases, energy. Using higher order decompositions we investigate integration schemes of up to fourth order and compare them to a well established fourth order predictor-corrector method. We demonstrate that these methods can be used with much larger time steps than the predictor-corrector method and thus may lead to a substantial speedup of computer simulations of the dynamical behavior of magnetic materials.