Application of Quadrature Methods for Re-Weighting in Lattice QCD

Abstract

Re-weighting is a useful tool that has been employed in Lattice QCD in different contexts including, tuning the strange quark mass, approaching the light quark mass regime, and simulating electromagnetic fields on top of QCD gauge configurations. In case of re-weighting the sea quark mass, the re-weighting factor is given by the ratio of the determinants of two Dirac operators $D_a$ and $D_b$. A popular approach for computing this ratio is to use a pseudofermion representation of the determinant of the composite operator $\Omega=D_a(D_b^\dagger D_b)^{-1} D_a^\dagger$. Here, we study using quadrature methods together with noise vectors to compute the ratio of determinants. We show that, with quadrature methods each determinant can be computed separately using the operators $\Omega_a=D_a^\dagger D_a$ and $\Omega_b=D_b^\dagger D_b$. We also discuss using bootstrap re-sampling to remove the bias from the determinant estimator.