Fully ergodic simulations using radial updates

Abstract

A sensible application of the Hybrid Monte Carlo (HMC) method is often hindered by the presence of large - or even infinite - potential barriers. These potential barriers separate the configuration space into distinct sectors and can lead to ergodicity violations that bias measurements. In this work, we address this problem by augmenting HMC with a multiplicative Metropolis-Hastings update in a so-called ''radial direction'' of the fields which enables crossing the potential barriers and ensures ergodicity of the sampling algorithm at comparably low computational cost. We demonstrate the algorithm on a simple toy model and show how it can be applied to the fermionic Hubbard model describing physics ranging from an exactly-solvable two-site system to the $C_{20}{H}_{12}$ perylene molecule. Our numerical results show that the radial updates successfully remove ergodicity violations, while simultaneously reducing autocorrelation times.

19 pages, 13 figures

References in corpus (12)

• Monte Carlo errors with less errors
• Revisiting the Hybrid Quantum Monte Carlo Method for Hubbard and Electron-Phonon Models
• Monte-Carlo simulation of the tight-binding model of graphene with partially screened Coulomb interactions
• Hybrid-Monte-Carlo study of competing order in the extended fermionic Hubbard model on the hexagonal lattice
• Quantum Monte Carlo Calculations for Carbon Nanotubes
• Numerical evidence of conformal phase transition in graphene with long-range interactions
• Markov Chain Monte Carlo Based on Deterministic Transformations
• Avoiding Ergodicity Problems in Lattice Discretizations of the Hubbard Model
• Multiplicative random walk Metropolis-Hastings on the real line
• Single Particle Spectrum of Doped ${\mathrm{C}}_{20}{\mathrm{H}}_{12}$-Perylene
• Exponential speed up in Monte Carlo sampling through Radial Updates
• Minimal Autocorrelation in Hybrid Monte Carlo simulations using Exact Fourier Acceleration