Bulk first-order phase transition in three-flavor lattice QCD with O(a)-improved Wilson fermion action at zero temperature -: art. no. 054510

Three-flavor QCD simulation with the O(a)-improved Wilson fermion action is made employing an exact fermion algorithm developed for an odd number of quark flavors. For the plaquette gauge action, an unexpected first-order phase transition is found in the strong coupling regime (beta less than or similar to 5.0) at relatively heavy quark masses (m(PS)/m(V)similar to 0.74-0.87). Strong metastability persists on a large lattice of size 12(3)x32, which indicates that the transition has a bulk nature. The phase gap becomes smaller toward weaker couplings and vanishes at beta similar or equal to 5.0, which corresponds to a lattice spacing a similar or equal to 0.1 fm. These results imply that realistic simulations of QCD with three flavors of dynamical Wilson-type fermions at lattice spacings in the range a=0.1-0.2 fm are not possible with the plaquette gauge action. Extending the study to improved gauge actions, we do not observe evidence for first-order phase transition, at least within the (beta,kappa) range we explored. This suggests the possibility that the phase transition either moves away or weakens with improved gauge actions. Possible origins of the phase transition are discussed.