# Non-hermitian Exact Local Bosonic Algorithm for Dynamical Quarks

## TL;DR

The paper introduces an exact local bosonic algorithm for lattice QCD simulations of dynamical quarks, utilizing a non-hermitian polynomial approximation and a Metropolis test to ensure accuracy and convergence.

## Contribution

It presents a novel exact algorithm based on non-hermitian polynomial approximation, with proven convergence and practical simulation results.

## Key findings

- Algorithm converges to the correct distribution.
- Performance scales with lattice volume and quark mass.
- Simulation data demonstrate effectiveness of the method.

## Abstract

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.