Lattice QCD with open boundary conditions and twisted-mass reweighting
Lattice QCD simulations at small lattice spacings and quark masses close to their physical values are technically challenging. In particular, the simulations can get trapped in the topological charge sectors of field space or may run into instabilities triggered by accidental near-zero modes of the lattice Dirac operator. As already noted in ref.[1], the first problem is bypassed if open boundary conditions are imposed in the time direction, while the second can potentially be overcome through twisted-mass determinant reweighting[2]. In this paper, we show that twisted-mass reweighting works out as expected in QCD with open boundary conditions and 2+1 2 + 1 flavours of O(a) ( a ) improved Wilson quarks. Further algorithmic improvements are tested as well and a few physical quantities are computed for illustration.
lattice qcd open boundary conditions hybrid monte carlo algorithm determinant reweighting