Multigrid low-mode averaging

TL;DR

This paper introduces a multigrid low-mode averaging technique that maintains a volume-independent number of low modes for variance reduction in lattice QCD, improving efficiency for large-volume simulations.

Contribution

It generalizes low-mode averaging by exploiting local coherence and hierarchical subspaces, enabling effective variance reduction independent of volume size.

Findings

Variance decreases with volume when using multigrid decomposition.

Effective suppression of fluctuations in quark-line connected diagrams at large separations.

Applicable to various quark propagator calculations, especially at large distances.

Abstract

We develop a generalization of low-mode averaging in which the number of low quark modes of the Dirac operator required for a constant variance reduction can be kept independent of the volume by exploiting their local coherence. Typically in lattice QCD simulations, the benefit of translation averaging quark propagators over the space-time volume is spoiled by large fluctuations introduced by the approximations needed to estimate the average. For quark-line connected diagrams at large separations, most of this additional variance can be efficiently suppressed by the introduction of hierarchical subspaces, thanks to the reduced size of the coarse grid operators that act within the subspaces. In this work, we investigate the contributions to the variance of the isovector vector current correlator with $N_{\mathrm f}=2$ non-perturbatively $\mathrm O(a)$-improved Wilson fermions on lattices approximately of size $L=2,3$ and $4$ $\mathrm {fm}$. The numerical results obtained confirm that the variance decreases as the volume is increased when a multigrid decomposition is used with a fixed number of low modes. While the proposed decomposition can be applied to any quark propagator, it is expected to be especially effective for quark-line connected diagrams at large separations, for example, the isovector contribution to the hadronic vacuum polarization or baryonic correlators.