There has been much recent progress in the understanding and reduction of the computational cost of the hybrid Monte Carlo algorithm for lattice QCD as the quark mass parameter is reduced. In this letter we present a new solution to this problem, where we represent the fermionic determinant using $n$ pseudofermion fields, each with an $n$th root kernel. We implement this within the framework of the rational hybrid Monte Carlo algorithm. We compare this algorithm with other recent methods in this area and find it is competitive with them.

• Received 22 August 2006

DOI:https://doi.org/10.1103/PhysRevLett.98.051601