Light hadron spectrum and quark masses from quenched lattice QCD

We present the details of simulations for the light hadron spectrum in quenched QCD carried out on the CP-PACS parallel computer. Simulations are made with the Wilson quark action and the plaquette gauge action on lattices of size 32(3)x56-64(3)x112 at four values of lattice spacings in the range aapproximate to0.1-0.05 fm and spatial extent L(s)aapproximate to3 fm. Hadronic observables are calculated at five quark masses corresponding to m(PS)/m(V)approximate to0.75-0.4, assuming the u and d quarks are degenerate, but treating the s quark separately. We find that the presence of quenched chiral singularities is supported from an analysis of the pseudoscalar meson data. The physical values of hadron masses are determined using m(pi), m(rho), and m(K) (or m(phi)) as input to fix the physical scale of lattice spacing and the u, d, and s quark masses. After chiral and continuum extrapolations, the agreement of the calculated mass spectrum with experiment is at a 10% level. In comparison with the statistical accuracy of 1%-3% and systematic errors of at most 1.7% we have achieved, this demonstrates a failure of the quenched approximation for the hadron spectrum: the hyperfine splitting in the meson sector is too small, and in the baryon sector the octet masses and mass splitting of the decuplet are both smaller than experiment. Light quark masses are calculated using two definitions: the conventional one and the one based on the axial-vector Ward identity. The two results converge toward the continuum limit, yielding m(ud)=4.29(14)(-0.79)(+0.51) MeV where the first error is statistical and the second one is systematic due to chiral extrapolation. The s quark mass depends on the strange hadron mass chosen for input: m(s)=113.8(2.3)(-2.9)(+5.8) MeV from m(K) and m(s)=142.3(5.8)(-0)(+22.0) MeV from m(phi), indicating again a failure of the quenched approximation. We obtain the scale of QCD, Lambda((M) over bar)((0))=219.5(5.4) MeV with m(rho) used as input. An O(10%) deviation from experiment is observed in the pseudoscalar meson decay constants