We present an exact dynamical QCD simulation algorithm for the $O(a)$-improved Wilson fermion with odd number of flavors. Our algorithm is an extension of the non-Hermitian polynomials HMC algorithm proposed by Takaishi and de Forcrand previously. In our algorithm, the systematic errors caused by the polynomial approximation of the inverse of Dirac operator is removed by a noisy-Metropolis test. For one flavor quark it is achieved by taking the square root of the correction matrix explicitly. We test our algorithm for the case of $N_f=1+1$ on a moderately large lattice size ($16^3times48$). The $N_f=2+1$ case is also investigated.
We report on a calculation of the light hadron spectrum and quark masses in three-flavor dynamical QCD using the non-perturbatively O(a)-improved Wilson quark action and a renormalization-group improved gauge action. Simulations are carried out on a 16^3 times 32 lattice at beta=1.9, where a^{-1} simeq 2GeV, with 6 ud quark masses corresponding to m_{pi}/m_{rho} simeq 0.64-0.77 and 2 s quark masses close to the physical value. We observe that the inclusion of dynamical strange quark brings the lattice QCD meson spectrum to good agreement with experiment. Dynamical strange quarks also lead to a reduction of the uds quark masses by about 15%.
We present simulation details and results for the light hadron spectrum in N f = 2 + 1 lattice QCD with the nonperturbatively O(a)-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at a lattice spacing of 0.09 fm on a (2.9fm)^3 box using the PACS-CS computer. We employ the Luschers domain-decomposed HMC algorithm with several improvements to reduce the degenerate up-down quark mass toward the physical value. So far the resulting pseudoscalar meson mass is ranging from 702MeV down to 156MeV. We discuss on the stability and the efficiency of the algorithm. The light harden spectrum extrapolated at the physical point is compared with the experimental values. We also present the values of the quark masses and the pseudoscalar meson decay constants.
We construct positive-definite pseudofermion actions for one fermion flavor in lattice field theory, for Wilson and domain-wall fermions respectively. The positive definiteness of these actions ensures that they can be simulated with the Hybrid Monte Carlo (HMC) method. For lattice QCD with optimal domain-wall quarks, we compare the efficiency of HMC simulations of 2-flavor and (1+1)-flavor, and find that the efficiency ratio is about 3:2.
We study the finite temperature phase structure for three-flavor QCD with a focus on locating the critical point which separates crossover and first order phase transition region in the chiral regime of the Columbia plot. In this study, we employ the Iwasaki gauge action and the non-perturvatively O($a$) improved Wilson-Clover fermion action. We discuss the finite size scaling analysis including the mixing of magnetization-like and energy-like observables. We carry out the continuum extrapolation of the critical point using newly generated data at $N_{rm t}=8$, $10$ and estimate the upper bound of the critical pseudo-scalar meson mass $m_{rm PS,E} lesssim 170 {rm MeV}$ and the critical temperature $T_{rm E}=134(3){rm MeV}$. Our estimate of the upper bound is derived from the existence of the critical point as an edge of the 1st order phase transition while that of the staggered-type fermions is based on its absence.
We compute the Landau gauge quark propagator from lattice QCD with two flavors of dynamical O(a)-improved Wilson fermions. The calculation is carried out with lattice spacings ranging from 0.06 fm to 0.08 fm, with quark masses corresponding to pion masses of 420, 290 and 150 MeV, and for volumes of up to (4.5fm)^4. Our ensembles allow us to evaluate lattice spacing, volume and quark mass effects. We find that the quark wave function which is suppressed in the infrared, is further suppressed as the quark mass is reduced, but the suppression is weakened as the volume is increased. The quark mass function M(p^2) shows only a weak volume dependence. Hypercubic artefacts beyond O(a) are reduced by applying both cylinder cuts and H4 extrapolations. The H4 extrapolation shifts the quark wave function systematically upwards but does not perform well for the mass function.