Covariant approximation averaging

INSPIRE: 1279923
arXiv: 1402.0244
DOI: 10.1103/PhysRevD.91.114511

Authors: Shintani, Eigo, Arthur, Rudy, Blum, Thomas, Izubuchi, Taku, Jung, Chulwoo, Lehner, Christoph

Submitted: 2 February 2014

Subjects:

Journal reference: Phys.Rev.D 91 (11) 114511 (2015)

Citations: 190

Abstract

We present a new class of statistical error reduction techniques for Monte Carlo simulations. Using covariant symmetries, we show that correlation functions can be constructed from inexpensive approximations without introducing any systematic bias in the final result. We introduce a new class of covariant approximation averaging techniques, known as all-mode averaging (AMA), in which the approximation takes account of contributions of all eigenmodes through the inverse of the Dirac operator computed from the conjugate gradient method with a relaxed stopping condition. In this paper we compare the performance and computational cost of our new method with traditional methods using correlation functions and masses of the pion, nucleon, and vector meson in Nf=2+1 lattice QCD using domain-wall fermions. This comparison indicates that AMA significantly reduces statistical errors in Monte Carlo calculations over conventional methods for the same cost.