Non-hermitian exact local bosonic algorithm for dynamical quarks Original Research Article

We present an exact version of the local bosonic algorithm for the simulation of dynamical quarks in lattice QCD. This version is based on a non-hermitian polynomial approximation of the inverse of the quark matrix. A Metropolis test corrects the systematic errors. Two variants of this test are presented. For both of them, a formal proof is given that this Monte Carlo algorithm converges to the right distribution. Simulation data are presented for different lattice parameters. The dynamics of the algorithm and its scaling in lattice volume and quark mass are investigated.