Abstract:Monte Carlo simulations of quantum field theories on a lattice become increasingly expensive as the continuum limit is approached since the cost per independent sample grows with a high power of the inverse lattice spacing. Simulations on fine lattices suffer from critical slowdown, the rapid growth of autocorrelations in the Markov chain. This causes a strong increase in the number of lattice configurations that have to be generated to obtain statistically significant results. This paper discusses hierarchical sampling methods to tame the growth in autocorrelations. Combined with multilevel variance reduction, this significantly reduces the computational cost of simulations for given tolerances $\epsilon_{\text{disc}}$ on the discretisation error and $\epsilon_{\text{stat}}$ on the statistical error. For observables with lattice errors of order $\alpha$ and integrated autocorrelation times that grow like $\tau_{\mathrm{int}}\propto a^{-z}$, multilevel Monte Carlo (MLMC) reduces the cost from $\mathcal{O}(\epsilon_{\text{stat}}^{-2}\epsilon_{\text{disc}}^{-(1+z)/\alpha})$ to $\mathcal{O}(\epsilon_{\text{stat}}^{-2}\vert\log \epsilon_{\text{disc}} \vert^2+\epsilon_{\text{disc}}^{-1/\alpha})$ or $\mathcal{O}(\epsilon_{\text{stat}}^{-2}+\epsilon_{\text{disc}}^{-1/\alpha})$. Higher gains are expected for simulations of quantum field theories in $D$ dimensions. The efficiency of the approach is demonstrated on two model systems, including a topological oscillator that is badly affected by critical slowdown from topological charge freezing. On fine lattices, the new methods are orders of magnitude faster than standard Hybrid Monte Carlo sampling. For high resolutions, MLMC can be used to accelerate even the cluster algorithm for the topological oscillator. Performance is further improved through perturbative matching which guarantees efficient coupling of theories on the multilevel hierarchy.

Self-learning Monte Carlo Method and Cumulative Update in Fermion Systems

Self-learning quantum Monte Carlo method in interacting fermion systems

Symmetry-enforced self-learning Monte Carlo method applied to the Holstein model

Self-learning Monte Carlo method: Continuous-time algorithm

Effective Ruderman–Kittel–Kasuya–Yosida-like Interaction in Diluted Double-exchange Model: Self-learning Monte Carlo Approach

Self-learning Monte Carlo with equivariant Transformer

Wide Sampling and Efficient Updating Monte Carlo Algorithms for Dimer Models

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

Accessing Excitation of Many-body Systems via Single-Mode Approximation within Quantum Monte Carlo Simulations

Elective-momentum Ultrasize Quantum Monte Carlo Method

An Exact Monte Carlo Method for Continuum Fermion Systems

Some aspects of simulation algorithms for dynamical fermions

Monte Carlo simulation of electrostatic interactions in inhomogeneous dielectric media: Correct sampling for the local lattice simulation algorithm

Efficient swap algorithms for molecular dynamics simulations of equilibrium supercooled liquids

Boosting Determinant Quantum Monte Carlo with Submatrix Updates: Unveiling the Phase Diagram of the 3D Hubbard Model

A multiscale hybrid Maxwellian-Monte-Carlo Coulomb collision algorithm for particle simulations

A quantitative study on the approximation error and speed-up of the multi-scale MCMC (Monte Carlo Markov chain) method for molecular dynamics

Machine learning force-field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation

Sign-Problem-Free Fermionic Quantum Monte Carlo: Developments and Applications

Multilevel Monte Carlo for quantum mechanics on a lattice