@article{oai:tsukuba.repo.nii.ac.jp:00007787,
 author = {Aoki, S. and Burkhalter, R. and Fukugita, M. and Hashimoto, S. and Ishikawa, K-I. and Ishizuka, N. and Iwasaki, Y. and 岩崎, 洋一 and Kanaya, K. and Kaneko, T. and Kuramashi, Y. and Okawa, M. and Onogi, T. and Tominaga, S. and Tsutsui, N. and Ukawa, A. and 宇川, 彰 and Yamada, N. and Yoshie, T.},
 issue = {9},
 journal = {Physical review D},
 month = {Apr},
 note = {We present a polynomial hybrid Monte Carlo (PHMC) algorithm for lattice QCD with odd numbers of flavors of O(a)-improved Wilson quark action. The algorithm makes use of the non-Hermitian Chebyshev polynomial to approximate the inverse square root of the fermion matrix required for an odd number of flavors. The systematic error from the polynomial approximation is removed by a noisy Metropolis test for which a new method is developed. Investigating the property of our PHMC algorithm in the Nf=2 QCD case, we find that it is as efficient as the conventional HMC algorithm for a moderately large lattice size (16^3×48) with intermediate quark masses (mPS/mV~0.7-0.8). We test our odd-flavor algorithm through extensive simulations of two-flavor QCD treated as an Nf=1+1 system, and comparing the results with those of the established algorithms for Nf=2 QCD. These tests establish that our PHMC algorithm works on a moderately large lattice size with intermediate quark masses (16^3×48,mPS/mV~0.7-0.8). Finally we experiment with the (2+1)-flavor QCD simulation on small lattices (4^3×8 and 8^3×16), and confirm the agreement of our results with those obtained with the R algorithm and extrapolated to a zero molecular dynamics step size.},
 title = {Polynomial hybrid Monte Carlo algorithm for lattice QCD with an odd number of flavors},
 volume = {65},
 year = {2002}
}