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.

Temperature Dependence of the Upper Critical Field in High-Temperature Superconductors: Localization Effects

Multigrid Monte Carlo Revisited: Theory and Bayesian Inference

Minimal Autocorrelation in Hybrid Monte Carlo simulations using Exact Fourier Acceleration

Toward Large-Scale AFQMC Calculations: Large Time Step Auxiliary-Field Quantum Monte Carlo

Yet Another Paper on Sn1987a: Large Angle Oscillations, and the Electron Neutrino Mass

A new algorithm for Monte Carlo simulation of Ising spin systems

The effect of branchless collisions and population control on correlations in Monte Carlo power iteration

Analysis of autocorrelation times in Neural Markov Chain Monte Carlo simulations

Towards Large-Scale AFQMC Calculations: Large Time Step Auxiliary-Field Quantum Monte Carlo

Multilevel Monte Carlo for quantum mechanics on a lattice

Hybrid Monte Carlo Simulation with Fourier Acceleration of the $N=2$ Principal Chiral Model in two Dimensions

Scheduling massively parallel multigrid for multilevel Monte Carlo methods

Multilevel Asymptotic-Preserving Monte Carlo for Particle Simulations

Monte Carlo coupled Multi-Physics with spatially continuous material properties

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Quadratic Speed-up in Infinite Variance Quantum Monte Carlo

Some recent developments in auxiliary-field quantum Monte Carlo for real materials

Parallelization, processor communication and error analysis in lattice kinetic Monte Carlo

Improving computational efficiency of Monte-Carlo simulations with variance reduction

Mitigating Spatial Error in the iterative-Quasi-Monte Carlo (iQMC) Method for Neutron Transport Simulations with Linear Discontinuous Source Tilting and Effective Scattering and Fission Rate Tallies

Parallel Tempered Metadynamics: Overcoming potential barriers without surfing or tunneling