We introduce a factorization of the fermion determinant in lattice QCD with Wilson-type fermions that leads to a bosonic action which is local in the block fields. The interaction among gauge fields on distant blocks is mediated by multiboson fields located on the boundaries of the blocks. The resultant multiboson domain-decomposed hybrid Monte Carlo passes extensive numerical tests carried out by measuring standard gluonic observables. The combination of the determinant factorization and of the one of the propagator, that we put forward recently, paves the way for multilevel Monte Carlo integration in the presence of fermions. We test this possibility by computing the disconnected correlator of two flavor-diagonal pseudoscalar densities, and we observe a significant increase of the signal-to-noise ratio due to a two-level integration.

• Received 15 September 2016

DOI:https://doi.org/10.1103/PhysRevD.95.034503