Abstract
For systems of classical spins interacting with the bath via damping and thermal noise, an approach is suggested to replace the white noise by a pulse noise acting at regular time intervals , within which the system evolves conservatively. The method is working well in the typical case of a small dimensionless damping constant and allows a considerable speedup of computations by using high-order numerical integrators with a large time step (up to a fraction of the precession period), while keeping to reduce the relative contribution of noise-related operations. In cases when precession can be discarded, can be increased up to a fraction of the relaxation time that leads to a further speedup. This makes equilibration speed comparable with that of the Metropolis Monte Carlo method. The pulse-noise approach is tested on single-spin and multispin models.
- Received 14 October 2016
DOI:https://doi.org/10.1103/PhysRevE.95.013306
©2017 American Physical Society