Rejection-free Glauber Monte Carlo for the 2D Random Field Ising Model via Hierarchical Probabilistic Counters

Abstract

We present an efficient Monte Carlo algorithm for the simulation of the two-dimensional Random Field Ising Model (RFIM). The method combines the event-driven, rejection-free character of the Bortz Kalos-Lebowitz (BKL) algorithm with Glauber transition probabilities, introducing hierarchical probabilistic counters to perform spin selection in O(log N) operations. This enables efficient sampling of the system's dynamics, especially in the low-temperature and low-disorder regime, where traditional Metropolis updates suffer from critical slowing down. Furthermore, this approach allows a proper dynamical simulation of the Ising system's behavior even in the presence of a Random Field (RF), unlike the BKL method. RFIM simulations with Gaussian field distributions reproduce the expected reduction of the pseudo-critical temperature with increasing disorder. Benchmarking shows speedups exceeding two orders of magnitude compared to the Metropolis algorithm in the low-temperature regime. The proposed method provides an efficient and dynamically faithful tool for studying both equilibrium and non-equilibrium phenomena in disordered spin systems.

arXiv comment16 pages, 11 figures

References in corpus (9)

• Exchange Monte Carlo Method and Application to Spin Glass Simulations
• Population Annealing with Weighted Averages: A Monte Carlo Method for Rough Free Energy Landscapes
• Criticality of the random field Ising model in and out of equilibrium: a nonperturbative functional renormalization group description
• Review of recent developments in the random-field Ising model
• Efficient numerical methods for the random-field Ising model: Finite-size scaling, reweighting extrapolation, and computation of response functions
• The Hobbyhorse of Magnetic Systems: The Ising Model
• Fast, hierarchical, and adaptive algorithm for Metropolis Monte Carlo simulations of long-range interacting systems
• A power-law upper bound on the correlations in the 2D random field Ising model
• Dynamical Properties of Random Field Ising Model