A numerical study of quenched QCD for light quarks is presented using O(a) improved fermions. Particular attention is paid to the possible existence and determination of quenched chiral logarithms. A `safe region to use for chiral extrapolations appears to be at and above the strange quark mass.
Quenched QCD simulations on three volumes, $8^3 times$, $12^3 times$ and $16^3 times 32$ and three couplings, $beta=5.7$, 5.85 and 6.0 using domain wall fermions provide a consistent picture of quenched QCD. We demonstrate that the small induced effects of chiral symmetry breaking inherent in this formulation can be described by a residual mass ($mres$) whose size decreases as the separation between the domain walls ($L_s$) is increased. However, at stronger couplings much larger values of $L_s$ are required to achieve a given physical value of $mres$. For $beta=6.0$ and $L_s=16$, we find $mres/m_s=0.033(3)$, while for $beta=5.7$, and $L_s=48$, $mres/m_s=0.074(5)$, where $m_s$ is the strange quark mass. These values are significantly smaller than those obtained from a more naive determination in our earlier studies. Important effects of topological near zero modes which should afflict an accurate quenched calculation are easily visible in both the chiral condensate and the pion propagator. These effects can be controlled by working at an appropriately large volume. A non-linear behavior of $m_pi^2$ in the limit of small quark mass suggests the presence of additional infrared subtlety in the quenched approximation. Good scaling is seen both in masses and in $f_pi$ over our entire range, with inverse lattice spacing varying between 1 and 2 GeV.
We compute hadron masses and the lowest moments of unpolarized and polarized nucleon structure functions down to pion masses of 300 MeV, in an effort to make unambiguous predictions at the physical light quark mass.
We investigate two-point correlation functions of left-handed currents computed in quenched lattice QCD with the Neuberger-Dirac operator. We consider two lattice spacings a~0.09,0.12 fm and two different lattice extents L~ 1.5, 2.0 fm; quark masses span both the p- and the epsilon-regimes. We compare the results with the predictions of quenched chiral perturbation theory, with the purpose of testing to what extent the effective theory reproduces quenched QCD at low energy. In the p-regime we test volume and quark mass dependence of the pseudoscalar decay constant and mass; in the epsilon-regime, we investigate volume and topology dependence of the correlators. While the leading order behaviour predicted by the effective theory is very well reproduced by the lattice data in the range of parameters that we explored, our numerical data are not precise enough to test next-to-leading order effects.
We investigate the chiral properties of quenched domain-wall QCD (DWQCD) at the lattice spacings $a^{-1} simeq 1$ and 2 GeV for both plaquette and renormalization-group (RG) improved gauge actions. In the case of the plaquette action we find that the quark mass defined through the axial Ward-Takahashi identity remains non-vanishing in the DWQCD chiral limit that the bare quark mass $m_fto 0$ and the length of the fifth dimension $N_stoinfty$, indicating that chiral symmetry is not realized with quenched DWQCD up to $a^{-1} simeq 2$ GeV. The behavior is much improved for the RG-improved gauge action: while a non-vanishing quark mass remains in the chiral limit at $a^{-1}simeq 1$ GeV, the result at $a^{-1}simeq 2$ GeV is consistent with an exponentially vanishing quark mass in the DWQCD chiral limit, indicating the realization of exact chiral symmetry. An interpretation and implications are briefly discussed.
We compute charm and bottom quark masses in the quenched approximation and in the continuum limit of lattice QCD. We make use of a step scaling method, previously introduced to deal with two scale problems, that allows to take the continuum limit of the lattice data. We determine the RGI quark masses and make the connection to the MSbar scheme. The continuum extrapolation gives us a value m_b^{RGI} = 6.73(16) GeV for the b-quark and m_c^{RGI} = 1.681(36) GeV for the c-quark, corresponding respectively to m_b^{MSbar}(m_b^{MSbar}) = 4.33(10) GeV and m_c^{MSbar}(m_c^{MSbar}) = 1.319(28) GeV. The latter result, in agreement with current estimates, is for us a check of the method. Using our results on the heavy quark masses we compute the mass of the Bc meson, M_{Bc} = 6.46(15) GeV.