Riemannian manifold hybrid Monte Carlo in lattice QCD

2021/12/08 by Tuan Nguyen, Peter Boyle, Nguyen, Tuan +7
Physics and Astronomy · #FOS: Physical sciences #High Energy Physics - Lattice (hep-lat) #hep-lat

paper · pdf · doi:10.48550/arxiv.2112.04556

7 pages, 4 figures, 2021 Lattice Proceedings

arxiv created 2021/12/08 · arxiv updated 2021/12/10

Abstract

Critical slowing down presents a critical obstacle to lattice QCD calculation at the smaller lattice spacings made possible by Exascale computers. Inspired by the concept of Fourier acceleration, we study a version of the Riemannian Manifold HMC (RMHMC) algorithm in which the canonical mass term of the HMC algorithm is replaced by a rational function of the SU(3) gauge covariant Laplacian. We have developed a suite of tools using Chebyshev filters based on the SU(3) gauge covariant Laplacian that provides the power spectra of both the gauge and fermion forces and determines the spectral dependence of the resulting RMHMC evolution of long- and short-distance QCD observables. These tools can be used to optimize the RMHMC mass term and to monitor the resulting acceleration mode-wise.

Related

• Riemannian Manifold HMC with fermions
• Ensemble Quasi-Newton HMC
• Fourier acceleration, the HMC algorithm and renormalizability
• Playing with the kinetic term in the HMC
• Modified Cholesky Riemann Manifold Hamiltonian Monte Carlo: Exploiting Sparsity for Fast Sampling of High-dimensional Targets