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.

Linear-scaling quantum Monte Carlo with non-orthogonal localized orbitals

Toward linear scaling auxiliary field quantum Monte Carlo with local natural orbitals

Linear Scaling Calculation of Maximally Localized Wannier Functions with Atomic Basis Set.

Linear Scaling Local Correlation Approach for Solving the Coupled Cluster Equations of Large Systems.

Robust linear-scaling optimization of compact localized orbitals in density functional theory

Multilevel Monte Carlo for quantum mechanics on a lattice

A Blocked Linear Method for Optimizing Large Parameter Sets in Variational Monte Carlo

Self-consistent field theory based molecular dynamics with linear system-size scaling

Linear scaling computation of the Fock matrix VII. Periodic Density Functional Theory at the $Γ$-point

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

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

Linear-Scaling Subspace-Iteration Algorithm with Optimally Localized Nonorthogonal Wave Functions for Kohn-Sham Density Functional Theory

Efficient Linear-Scaling Quantum Transport Calculations on Graphics Processing Units and Applications on Electron Transport in Graphene

An Efficient Linear Scaling Procedure for Constructing Localized Orbitals of Large Molecules Based on the One-Particle Density Matrix

Electronic structure calculations and molecular dynamics simulations with linear system-size scaling

Automated quantum conductance calculations using maximally-localised Wannier functions

A quantum Monte Carlo algorithm for arbitrary spin-1/2 Hamiltonians

The sign problem in full configuration interaction quantum Monte Carlo: Linear and sub-linear representation regimes for the exact wave function

Local Second Order Mo̷ller-Plesset Theory with a Single Threshold Using Orthogonal Virtual Orbitals: A Distributed Memory Implementation

A Sub-Linear Scaling Algorithm for Computing the Electronic Structure of Materials

Monte Carlo based techniques for quantum magnets with long-range interactions