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.

Small-scale Hamiltonian optimization of interpolating operators for Lagrangian lattice quantum field theory

Strategies for quantum-optimized construction of interpolating operators in classical simulations of lattice quantum field theories

General quantum algorithms for Hamiltonian simulation with applications to a non-Abelian lattice gauge theory

Riemannian quantum circuit optimization for Hamiltonian simulation

Quantum Simulation of Lattice QCD with Improved Hamiltonians

An efficient finite-resource formulation of non-Abelian lattice gauge theories beyond one dimension

Towards determining the (2+1)-dimensional Quantum Electrodynamics running coupling with Monte Carlo and quantum computing methods

Hamiltonian truncation tensor networks for quantum field theories

Bridging lattice-scale physics and continuum field theory with quantum Monte Carlo simulations

Multilevel Monte Carlo for quantum mechanics on a lattice

Strategies for simulating time evolution of Hamiltonian lattice field theories

A new basis for Hamiltonian SU(2) simulations

Efficient Quantum Simulation Algorithms in the Path Integral Formulation

Simulating scalar field theories on quantum computers with limited resources

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

Multiparticle interpolating operators in quantum field theories with cubic symmetry

New basis for Hamiltonian SU(2) simulations

Helical Luttinger liquid on a space-time lattice

Digital quantum simulation of lattice fermion theories with local encoding

Worldline Monte Carlo method for few body nuclear physics