Minimal Autocorrelation in Hybrid Monte Carlo simulations using Exact Fourier Acceleration
arXiv:2404.09723 · doi:10.1016/j.cpc.2025.109624
Abstract
The hybrid Monte Carlo (HMC) algorithm is a ubiquitous method in computational physics with applications ranging from condensed matter to lattice QCD and beyond. However, HMC simulations often suffer from long autocorrelation times, severely reducing their efficiency. In this work two of the main sources of autocorrelations are identified and eliminated. The first source is the sampling of the canonical momenta from a sub-optimal normal distribution, the second is a badly chosen trajectory length. Analytic solutions to both problems are presented and implemented in the exact Fourier acceleration (EFA) method. It completely removes autocorrelations for near-harmonic potentials and consistently yields (close-to-) optimal results for numerical simulations of the Su-Schrieffer-Heeger and the Ising models as well as in lattice gauge theory, in some cases reducing the autocorrelation by multiple orders of magnitude. EFA is advantageous for and easily applicable to any HMC simulation of an action that includes a quadratic part.
arXiv comment15 + 8 pages, 4 + 1 figures, 3 algorithms; v2: removed gauge fixing; v3: minor revision, published
References in corpus (25)
- MCMC using Hamiltonian dynamics
- Normalizing Flows: An Introduction and Review of Current Methods
- Ab initio calculation of the neutron-proton mass difference
- Monte Carlo errors with less errors
- Critical slowing down and error analysis in lattice QCD simulations
- HMC algorithm with multiple time scale integration and mass preconditioning
- Quantum Monte Carlo Simulations of the 2D Su-Schrieffer-Heeger Model
- Valence-bond solid to antiferromagnet transition in the two-dimensional Su-Schrieffer-Heeger model by Langevin dynamics
- Langevin Simulations of a Long Range Electron Phonon Model
- The Semimetal-Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice
- Optimised Trotter Decompositions for Classical and Quantum Computing
- Exploring the HMC trajectory-length dependence of autocorrelation times in lattice QCD
- Quantum Monte Carlo Calculations for Carbon Nanotubes
- Testing a Fourier Accelerated Hybrid Monte Carlo Algorithm
- SmoQyDQMC.jl: A flexible implementation of determinant quantum Monte Carlo for Hubbard and electron-phonon interactions (version 2.0 release)
- Accelerating Hybrid Monte Carlo simulations of the Hubbard model on the hexagonal lattice
- Schur complement solver for Quantum Monte-Carlo simulations of strongly interacting fermions
- First-principle quantum Monte-Carlo study of charge carrier mobility in organic molecular semiconductors
- The phase diagram of the Sigma Model and its Implications for Chiral Hierarchies
- The Ising Model with Hybrid Monte Carlo
- Gauge-Fixed Fourier Acceleration
- Lattice sine-Gordon model
- Overcoming Ergodicity Problems of the Hybrid Monte Carlo Method using Radial Updates
- Hybrid Monte Carlo Simulation with Fourier Acceleration of the Principal Chiral Model in two Dimensions
- Fourier Acceleration in a Linear Sigma Model with Spontaneous Symmetry Breaking