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.

Exploiting stochastic locality in lattice QCD: hadronic observables and their uncertainties

Optimising stochastic algorithms for hadron correlation function computations in lattice QCD using a localised distillation basis

Leveraging neural control variates for enhanced precision in lattice field theory

A stochastic approach to the reconstruction of spectral functions in lattice QCD

Reconstructing lattice QCD spectral functions with stochastic pole expansion and Nevanlinna analytic continuation

Stochastic and Long-Distance Level Spacing Statistics in Many-Body Localization

Life Outside the Golden Window: Statistical Angles on the Signal-to-Noise Problem

Hadron Spectroscopy with lattice QCD: challenges and opportunities

Locality with staggered fermions

Two-nucleon S-wave interactions at the $SU(3)$ flavor-symmetric point with $m_{ud}\simeq m_s^{\rm phys}$: a first lattice QCD calculation with the stochastic Laplacian Heaviside method

High precision applications of lattice gauge theories in the quest for new physics

Bootstrap-determined p-values in Lattice QCD

Toward inclusive observables with staggered quarks: the smeared $R$~ratio

Multilevel Monte Carlo for quantum mechanics on a lattice

Quantum transport and localization in 1d and 2d tight-binding lattices

On the Locality and Scaling of Overlap Fermions at Coarse Lattice Spacings

Locality in collider tests of Quantum Mechanics with top quark pairs

Dynamical Localization in $\mathbb{Z}_2$ Lattice Gauge Theories