An Exact Algorithm for Any-flavor Lattice QCD with Kogut-Susskind Fermion
arXiv:hep-lat/0208058 · doi:10.1016/S0010-4655(03)00353-9
Abstract
We propose an exact simulation algorithm for lattice QCD with dynamical Kogut-Susskind fermion in which the N_f-flavor fermion operator is defined as the N_f/4-th root of the Kogut-Susskind (KS) fermion operator. The algorithm is an extension of the Polynomial Hybrid Monte Carlo (PHMC) algorithm to KS fermions. The fractional power of the KS fermion operator is approximated with a Hermitian Chebyshev polynomial, with which we can construct an algorithm for any number of flavors. The error which arises from the approximation is corrected by the Kennedy-Kuti noisy Metropolis test. Numerical simulations are performed for the two-flavor case for several lattice parameters in order to confirm the validity and the practical feasibility of the algorithm. In particular tests on a 16^4 lattice with a quark mass corresponding to m_{PS}/m_V ~ 0.68 are successfully accomplished. We conclude that our algorithm provides an attractive exact method for dynamical QCD simulations with KS fermions.
35 pages, 12 figures, uses elsart.cls, v2: typos corrected, comparison to the R-algorithm on large size lattice is added
References in corpus (8)
- The QCD spectrum with three quark flavors
- Light Hadron Spectroscopy with Two Flavors of Dynamical Quarks on the Lattice
- Numerical Methods for the QCD Overlap Operator: I. Sign-Function and Error Bounds
- Effects of non-perturbatively improved dynamical fermions in QCD at fixed lattice spacing
- Polynomial Hybrid Monte Carlo algorithm for lattice QCD with an odd number of flavors
- Dynamical Fermions with Fat Links
- Odd-flavor Hybrid Monte Carlo Algorithm for Lattice QCD
- Light Quark Masses with Wilson Fermions