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.

Mitigating spikes in fermion Monte Carlo methods by reshuffling measurements

A solution for infinite variance problem of fermionic observables

Post-measurement Quantum Monte Carlo

Addressing general measurements in quantum Monte Carlo

An integral algorithm of exponential observables for interacting fermions in quantum Monte Carlo simulation

Quantum Monte Carlo for gauge fields and matter without the fermion determinant

Multilevel Monte Carlo for quantum mechanics on a lattice

Systematic lowering of the scaling of Monte Carlo calculations by partitioning andsubsampling

Computation of Dynamical Correlation Functions for Many Fermion Systems with Auxiliary-Field Quantum Monte Carlo

On the Infinite Variance Problem in Fermion Models

Non-equilibrium quantum Monte Carlo algorithm for stabilizer Renyi entropy in spin systems

A Simple and Efficient Joint Measurement Strategy for Estimating Fermionic Observables and Hamiltonians

Optimal Fermionic Joint Measurements for Estimating Non-Commuting Majorana Observables

Elective-momentum Ultrasize Quantum Monte Carlo Method

High-efficiency quantum Monte Carlo algorithm for extracting entanglement entropy in interacting fermion systems

Quantum-assisted Monte Carlo algorithms for fermions

Measuring Rényi entanglement entropy with high efficiency and precision in quantum Monte Carlo simulations

Self-learning quantum Monte Carlo method in interacting fermion systems

Overcoming Ergodicity Problems of the Hybrid Monte Carlo Method using Radial Updates

Fermionic physics from ab initio path integral Monte Carlo simulations of fictitious identical particles