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

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: Cattaneo, Luca, Ettori, Federico, Cerri, Giovanni, Biscari, Paolo, Puppin, Ezio
Natura: Preprint
Pubblicazione: 2026
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866910057518596096
author Cattaneo, Luca
Ettori, Federico
Cerri, Giovanni
Biscari, Paolo
Puppin, Ezio
author_facet Cattaneo, Luca
Ettori, Federico
Cerri, Giovanni
Biscari, Paolo
Puppin, Ezio
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2603_17032
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Rejection-free Glauber Monte Carlo for the 2D Random Field Ising Model via Hierarchical Probabilistic Counters
Cattaneo, Luca
Ettori, Federico
Cerri, Giovanni
Biscari, Paolo
Puppin, Ezio
Statistical Mechanics
Disordered Systems and Neural Networks
Computational Physics
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.
title Rejection-free Glauber Monte Carlo for the 2D Random Field Ising Model via Hierarchical Probabilistic Counters
topic Statistical Mechanics
Disordered Systems and Neural Networks
Computational Physics
url https://arxiv.org/abs/2603.17032