Reweighting towards the chiral limit

Abstract

We propose to perform fully dynamical simulations at small quark masses by reweighting in the quark mass. This approach avoids some of the technical difficulties associated with direct simulations at very small quark masses. We calculate the weight factors stochastically, using determinant breakup and low mode projection to reduce the statistical fluctuations. We find that the weight factors fluctuate only moderately on nHYP smeared dynamical Wilson-clover ensembles, and we could successfully reweight 16^4, (1.85fm)^4 volume configurations from m_q = 20MeV to m_q = 5MeV quark masses, reaching the epsilon-regime. We illustrate the strength of the method by calculating the low energy constant F from the epsilon-regime pseudo-scalar correlator.