English

Fast, hierarchical, and adaptive algorithm for Metropolis Monte Carlo simulations of long-range interacting systems

Computational Physics 2023-07-27 v1 Statistical Mechanics

Abstract

We present a fast, hierarchical, and adaptive algorithm for Metropolis Monte Carlo simulations of systems with long-range interactions that reproduces the dynamics of a standard implementation exactly, i.e., the generated configurations and consequently all measured observables are identical, allowing in particular for nonequilibrium studies. The method is demonstrated for the power-law interacting long-range Ising model with nonconserved order parameter and a Lennard-Jones system both in two dimensions. The measured runtimes support an average complexity O(NlogN)O(N\log N), where NN is the number of spins or particles. Importantly, prefactors of this scaling behavior are small, which in practice manifests in speedup factors larger than 10410^4. The method is general and will allow the treatment of large systems that were out of reach before, likely enabling a more detailed understanding of physical phenomena rooted in long-range interactions.

Keywords

Cite

@article{arxiv.2207.14670,
  title  = {Fast, hierarchical, and adaptive algorithm for Metropolis Monte Carlo simulations of long-range interacting systems},
  author = {Fabio Müller and Henrik Christiansen and Stefan Schnabel and Wolfhard Janke},
  journal= {arXiv preprint arXiv:2207.14670},
  year   = {2023}
}