Matrix elements relevant for Delta I=1/2 rule and epsilon-prime from Lattice QCD with staggered fermions

We perform a study of matrix elements relevant for the Delta I=1/2 rule and the direct CP-violation parameter epsilon-prime from first principles by computer simulation in Lattice QCD. We use staggered (Kogut-Susskind) fermions, and employ the chiral perturbation theory method for studying K to 2 Pi decays. Having obtained a reasonable statistical accuracy, we observe an enhancement of the Delta I=1/2 amplitude, consistent with experiment within our large systematic errors. Finite volume and quenching effects have been studied and were found small compared to noise. The estimates of epsilon-prime are hindered by large uncertainties associated with operator matching. In this paper we explain the simulation method, present the results and address the systematic uncertainties.Comment: 40 pages, 17 figures, LATEX with epsf, to be submitted to Phys. Rev. D. Minor errors are corrected, some wording and notation change

Lattice Calculation of Heavy-Light Decay Constants with Two Flavors of Dynamical Quarks

We present results for $f_B$, $f_{B_s}$, $f_D$, $f_{D_s}$ and their ratios in the presence of two flavors of light sea quarks ($N_f=2$). We use Wilson light valence quarks and Wilson and static heavy valence quarks; the sea quarks are simulated with staggered fermions. Additional quenched simulations with nonperturbatively improved clover fermions allow us to improve our control of the continuum extrapolation. For our central values the masses of the sea quarks are not extrapolated to the physical $u$, $d$ masses; that is, the central values are "partially quenched." A calculation using "fat-link clover" valence fermions is also discussed but is not included in our final results. We find, for example, $f_B = 190 (7) (^{+24}_{-17}) (^{+11}_{-2}) (^{+8}_{-0})$ MeV, $f_{B_s}/f_B = 1.16 (1) (2) (2) (^{+4}_{-0})$, $f_{D_s} = 241 (5) (^{+27}_{-26}) (^{+9}_{-4}) (^{+5}_{-0})$ MeV, and $f_{B}/f_{D_s} = 0.79 (2) (^{+5}_{-4}) (3) (^{+5}_{-0})$, where in each case the first error is statistical and the remaining three are systematic: the error within the partially quenched $N_f=2$ approximation, the error due to the missing strange sea quark and to partial quenching, and an estimate of the effects of chiral logarithms at small quark mass. The last error, though quite significant in decay constant ratios, appears to be smaller than has been recently suggested by Kronfeld and Ryan, and Yamada. We emphasize, however, that as in other lattice computations to date, the lattice $u,d$ quark masses are not very light and chiral log effects may not be fully under control.Comment: Revised version includes an attempt to estimate the effects of chiral logarithms at small quark mass; central values are unchanged but one more systematic error has been added. Sections III E and V D are completely new; some changes for clarity have also been made elsewhere. 82 pages; 32 figure

Quenched Lattice QCD with Domain Wall Fermions and the Chiral Limit

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.Comment: 91 pages, 34 figure