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.

Field-Level Inference with Microcanonical Langevin Monte Carlo

Tracking the Multifield Dynamics with Cosmological Data: A Monte Carlo approach

The future of cosmological likelihood-based inference: accelerated high-dimensional parameter estimation and model comparison

Hamiltonian and Langevin Monte Carlo

Multilevel Monte Carlo for quantum mechanics on a lattice

A comparison of Bayesian sampling algorithms for high-dimensional particle physics and cosmology applications

Constrained cosmological simulations of the Local Group using Bayesian hierarchical field-level inference

Multigrid Monte Carlo Revisited: Theory and Bayesian Inference

Multilevel Sequential${}^2$ Monte Carlo for Bayesian Inverse Problems

Field-level simulation-based inference with galaxy catalogs: the impact of systematic effects

Microcanonical Hamiltonian Monte Carlo

Multilevel Importance Sampling for Rare Events Associated With the McKean--Vlasov Equation

Diffusion-HMC: Parameter Inference with Diffusion Model driven Hamiltonian Monte Carlo

ELUCID - Exploring the Local Universe with reConstructed Initial Density field I: Hamiltonian Markov Chain Monte Carlo Method with Particle Mesh Dynamics

Mean-Field Simulation-Based Inference for Cosmological Initial Conditions

Euclid-era cosmology for everyone: Neural net assisted MCMC sampling for the joint 3x2 likelihood

Fluctuation without dissipation: Microcanonical Langevin Monte Carlo

Fast likelihood-free inference in the LSS Stage IV era

Field-level multiprobe analysis of the CMB, integrated Sachs-Wolfe effect, and the galaxy density maps

Variational Hybrid Monte Carlo for Efficient Multi-Modal Data Sampling

Markov Chain Monte Carlo for Bayesian Parametric Galaxy Modeling in LSST