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.

An adaptive aggregation based domain decomposition multilevel method for the lattice wilson dirac operator: multilevel results

Coarsest-level improvements in multigrid for lattice QCD on large-scale computers

Multilevel Monte Carlo for quantum mechanics on a lattice

Accelerating Lattice QCD Simulations using GPUs

Extending DD-$α$AMG on heterogeneous machines

Improved Dirichlet boundary conditions for lattice gauge-fermion theories

Hierarchically deflated conjugate gradient

Scalable Multilevel Monte Carlo Methods Exploiting Parallel Redistribution on Coarse Levels

Deflation and Flexible SAP-Preconditioning of GMRES in Lattice QCD Simulation

A scalable multilevel domain decomposition preconditioner with a subspace-based coarsening algorithm for the neutron transport calculations

Deflation of eigenvalues for iterative methods in lattice QCD

Deflated GMRES with multigrid for lattice QCD

On the computation of lattice sums without translational invariance

An Additive Overlapping Domain Decomposition Method for the Helmholtz Equation

Parallel Scheme for Multi-Layer Refinement Non-Uniform Grid Lattice Boltzmann Method Based on Load Balancing

Parallelized Domain Decomposition for Multi-Dimensional Lagrangian Random Walk, Mass-Transfer Particle Tracking Schemes

Auxiliary Field Quantum Monte Carlo for Dilute Neutrons on the Lattice

Multigrid low-mode averaging

A subcycling/non-subcycling time advancement scheme-based DLM immersed boundary method framework for solving single and multiphase fluid–structure interaction problems on dynamically adaptive grids

Solving the left-hand cut problem in lattice QCD: from finite volume energy levels

Enhancing Scalability of a Matrix-Free Eigensolver for Studying Many-Body Localization