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 323×56–643×112 at four values of lattice spacings in the range a≈0.1–0.05 fm and spatial extent Lsa≈3 fm. Hadronic observables are calculated at five quark masses corresponding to mPS/mV≈0.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π, mρ, and mK (or mφ) 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 mud=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: ms=113.8(2.3)-2.9+5.8 MeV from mK and ms=142.3(5.8)-0+22.0 MeV from mφ, indicating again a failure of the quenched approximation. We obtain the scale of QCD, ΛMS̅ (0)=219.5(5.4) MeV with mρ used as input. An O(10%) deviation from experiment is observed in the pseudoscalar meson decay constants