Reweighting towards the chiral limit

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 164, (1.85 fm)4 volume configurations from mq≅20 MeV to mq≅5 MeV quark masses, reaching the ε-regime. We illustrate the strength of the method by calculating the low energy constant F from the ε-regime pseudoscalar correlator.