Hybrid Monte Carlo algorithm for lattice QCD with two flavors of dynamical Ginsparg-Wilson quarks

We study aspects concerning numerical simulations of lattice QCD with two flavors of dynamical Ginsparg-Wilson quarks with degenerate masses. A hybrid Monte Carlo algorithm is described and a formula for the fermionic force is derived for two specific implementations. The implementation with the optimal rational approximation method is favored in both CPU time and memory consumption. (C) 1999 Published by Elsevier Science B.V. All rights reserved.Physics, Particles & FieldsSCI(E)10ARTICLE1-2313-32255