Optimal phase space sampling for Monte Carlo simulations of Heisenberg spin systems

J. D. Alzate-Cardona*, D. Sabogal-Suárez, R. F.L. Evans, E. Restrepo-Parra

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review


We present an adaptive algorithm for the optimal phase space sampling in Monte Carlo simulations of 3D Heisenberg spin systems. Based on a golden rule of the Metropolis algorithm which states that an acceptance rate of 50% is ideal to efficiently sample the phase space, the algorithm adaptively modifies a cone-based spin update method keeping the acceptance rate close to 50%. We have assessed the efficiency of the adaptive algorithm through four different tests and contrasted its performance with that of other common spin update methods. In systems at low and high temperatures and anisotropies, the adaptive algorithm proved to be the most efficient for magnetization reversal and for the convergence to equilibrium of the thermal averages and the coercivity in hysteresis calculations. Thus, the adaptive algorithm can be used to significantly reduce the computational cost in Monte Carlo simulations of 3D Heisenberg spin systems.

Original languageEnglish
Article number095802
Number of pages10
JournalJournal of Physics Condensed Matter
Issue number9
Early online date17 Jan 2019
Publication statusPublished - 6 Mar 2019

Bibliographical note

© 2019 IOP Publishing Ltd


  • Heisenberg model
  • Monte Carlo
  • phase space sampling

Cite this