Symplectic integrators for classical spin systems

We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully utilized for other Hamiltonian systems, e. g. for molecular dynamics or non-linear wave equations. Our procedure rests on a decomposition of the spin Hamiltonian into a sum of two completely integrable Hamiltonians and on the corresponding Lie-Trotter decomposition of the time evolution operator. In order to make this method widely applicable we provide a large class of integrable spin systems whose time evolution consists of a sequence of rotations about fixed axes. We test the proposed symplectic integrator for small spin systems, including the model of a recently synthesized magnetic molecule, and compare the results for variants of different order.

Abstract licence: CC0

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.