# Full text of "$K\toππ$ Decays with Domain Wall Fermions: Lattice Matrix Elements"

## See other formats

K ^ TTTT Decays with Domain Wall Fermions: Lattice Matrix Elements T. Blum*^ RBC Collaboration * ^RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973-5000, USA We present a lattice calculation of the K ^ n and K matrix elements of the AS' — 1 effective weak Hamiltonian which can be used to determine e' /e and the AI = 1/2 rule for K decays in the Standard Model. The matrix elements for K —> tttv decays are related to A" tt and K using lowest order chiral perturbation theory. We also present results for the kaon B parameter, Bk- Our quenched domain wall fermion simulation was done at /3 = 6.0 (a~^ ^ 2 GeV), lattice size 16^ x 32 x 16, and domain wall height Ms = 1.8. 1. Introduction Recent measurements of direct CP violation (e'/e 7^ 0) in if ^ TTTT decays at FNAL and CERN allow an important test of the Standard Model (in particular, the CKM mixing paradigm). The ef- fective weak Hamiltonian governing strangeness changing K decays has been computed to next- to-leading order in QCD and QED by the Mu- nich and Rome groups; the remaining piece of the puzzle is the hadronic matrix elements of the operators of this effective weak Hamiltonian. The recent advance of domain wall (and over- lap) fermions which maintain chiral symmetry to a high degree of accuracy allows for a new attempt at this old problem. Chiral sym- metry of domain wall fermions provides a sig- nificant advantage when computing light quark QCD observables since the lattice artifacts that arise when this symmetry is explicitly broken are greatly reduced. Mixing and renormaliza- tion of operators, which is already complicated in the continuum, is readily handled with domain wall fermions]^. However, calculations using im- proved Wilson fermions were also reported at this meeting and it is still unclear which method will prove most advantageous. *This work was done in collaboration with N. Christ, C. Cristian, C. Dawson, G. Fleming, X. Liao, G. Liu, S. Ohta, A. Soni, P. Vranas, R. Mawhinney, M. Wingate, L. Wu, and Y. Zhestkov. We thank RIKEN, Brookhaven National Laboratory and the U.S. Department of Energy for providing the facilities essential for this work. In this study, we present preliminary results for the K ^ n and X — > matrix elements of these operators, which when combined with lowest or- der chiral perturbation theory, yield the desired K TTTT matrix elements Q. The contribution of R. Mawhinney in these proceedings takes up this point Here we are concerned with the simpler K ^ TT and K ^ lattice matrix elements. Also, the CP-PACS collaboration has presented a very similar calculation at this meeting]^]. 2. Theoretical Framework The AS* = 1 effective weak Hamiltonian is gen- erated from the fundamental Standard Model La- grangian by integrating out the top quark and W boson. The resulting effective Hamiltonian is then evolved to a low scale (/i <ti Mw) appropri- ate for lattice calculations using the renormaliza- tion group equations. The effective Hamiltonian above the charm threshold is: Gp H AS=1 V2 1 {Cli^l)iQl{^i) - QM) C2(/i)(Q2(M)-Q2c(M))) VudV- (1) where fi is the renormalization scale, Q are a basis of local four-quark operators which are closed under renormalization, C'i{ii) the cor- responding Wilson coefficients, and Vqq' the 2 Cabibbo-Kobayashi-Maskawa mixing matrix ele- ments which are fundamental parameters of the Standard Model. The effective operators renormalized at the scale /i are Ql = SaluPLdaUp-fyPLUp = SalvPhda ^ qpltiliyP(L,R)q(3 u,d,s,c... = SalvPhdp ^ qi3ltJ.lvP(L,B.)qa u.d,s^c... 3 ^ 13,5 Q 4,6 u.d.s.c. , Qlc = SalvPhdaCp^^PLCfj Q2c = SaluPLdpCp-fuPhCa, (2) with color indices a and /3, P(l,r) = 1 T 75j the sums are taken over active quark flavors at the scale /x, and summation over v is implied. Ql, 2,1c, 2c are often referred to as current-current operators, Q3-6 QCD penguin operators, and Qt-iq electroweak penguin operators. It is use- ful to split the above operators according to their isospin, Qr = Qf'''^ +Qf'^\ As mentioned earlier, lowest order chiral per- turbation theory relates K tttt matrix elements to a linear combination of K it and X — s- 0. For all Qi except the electroweak penguins Qj^s and all particles at rest, we haveQ (^+7r-|Q|i^°) = (n+\Q\K+) = 4i(m|. — ml)ai P - 4i(m|- — / where niM is the meson mass for unphysical pseu- doscalar states with vris = md- For Qi which transform in a (27,1) chiral multiplet a2 — 0. Note that each matrix element vanishes linearly with the meson mass squared. This is an im- portant prediction of chiral perturbation theory. and therefore QCD, and provides a solid test of the chiral symmetry properties of domain wall fermions. The strength of the above approach is that it allows less computationally demanding K IT and A' matrix elements to be cal- culated on the lattice. A significant drawback to this approach is that it manifestly does not con- tain information on the final state interactions of the pions (for calculation of e' the final state s- wave scattering phases from experiment can be put in by hand, however). In the case of the electroweak penguins Qj^s^ the contribution in lowest order chiral perturba- tion theory to the K ^ tt matrix elements is con- stant in the chiral limit. The K ^ matrix ele- ment, however, still vanishes. Since on the lattice 02 is quadratically diver- gent, the process of combining the second and third lines in Eq. 4 to obtain an expression for K ^ TTTT requires a delicate cancellation of this quadratic divergence. See R. Mawhinney's con- tribution for details. Finally, lattice counterparts of the operators in Eq. 1^ must be matched to the continuum and renormalized since they are logarithmically di- vergent after power divergences have been sub- tracted. In addition, operators in the same symmetry multiplets mix through renormaliza- tion group running from M\y down to the low scale /i. This poses a serious challenge for lat- tice calculations. In our calculation this prob- lem is handled remarkably well with the nonper- turbative renormalization method of the Rome- Southampton group which was explained in the talk by C. Dawsong. 3. Simulation details We have calculated matrix elements on 200 quenched gauge configurations at /3 = 6.0, with lattice four volume 16^ x 32, domain wall fermion extra dimension size Lg = 16, and domain wall height Afs = 1.8. We have calculated with light quark masses m/ = 0.01 — 0.05, and charm quark masses rric = 0.1 — 0.4. The physical kaon state made from degenerate quarks corresponds to nif Ri 0.02, and nic ~ 0.5 for the physical charm quark. 3 We extract matrix elements from three-point correlation functions. The external pseudoscalar states are interpolated from wall sources near the time direction boundaries, t = 5 and 27, and the operator is inserted between them. When the operator is far from either boundary, the de- sired lowest mass states dominate the correlation function. The forward and backward (in time) quark propagators used to interpolate the K and TT states are linear combinations of propagators computed with periodic and anti-periodic bound- ary conditions which amounts to doubling the gauge field configuration in the time direction. The closed fermion loops necessary for operators that have self contractions are computed from a complex Gaussian random source spread over time slices 14-17. All results are given as aver- ages over these four time slices. We have performed several important checks of our computer code. Most importantly, a com- pletely independent check code was written to compare with our two production versions (the check code and one of the production codes run on the QCDSP supercomputer and the other pro- duction code, based on the MILC code, runs on the NERSC T3E). Output from each code gen- erated on the same configuration agreed up to machine precision. As a final useful check, the left-left operators in Eq. |^ go into themselves under a Fierz trans- formation. Thus color-mixed contractions can be compared to corresponding color-diagonal ones. We find perfect agreement in all cases. The following results were obtained on the RIKEN BNL and Columbia University QCDSP supercomputers. 4. Results Fig. I shows (7r|Q^^/^^|X) as a function of quark mass, w/ = 77is — m^. An uncorrelated linear extrapolation yields a zero intercept, within statistical errors, which is in agreement with chi- ral perturbation theory. For strictly low energy QCD observables, we expect quantities to vanish ai irif = —rrires^- Since the AI — 1/2 opera- tors have contributions from physics scales near the (high energy) lattice cut-off, this is no longer true. Thus, for this matrix element, the statisti- cal errors are not small enough to resolve these systematic effects. However, in R. Mawhinney's contribution, we see that such effects are visible in the subtracted operator^. Presumably this is due to the strong statistical correlations between the A' — > TT matrix element and the subtraction term. In Fig. || we show a similar plot for Qq. 0.000 -0.001 - -0.002 - -0.003 - -0.004 0.00 0.01 0.02 0.03 0.04 0.05 Figure 1. The K it matrix element of the bare operator Q^2^'^^ . Here, an uncorrelated linear extrapolation has a non-zero intercept of roughly three standard de- viations. Note that it vanishes for m/ > 0. Thus explicit chiral symmetry breaking effects are visi- ble, though small. In Fig. ||we show (ttIQ^/^^IA-) which exhibits noticeable nonlinearity and does not vanish in the chiral limit, as expected. Finally, we show an example of a X — > ma- trix element in Fig. |. Note that for the AT ^ matrix elements we use ^ and fit to the form i^\Qi\K) I (^\s^^d\K) = const + {rus — md)r]i where chiral perturbation theory predicts const — 0. This ratio is useful since it is exactly the coef- ficient of the subtraction operator used to remove the quadratic divergence in {Tr\Qi\K). Since the quark masses enter as a difference, we expect ex- 4 0.4 0.3 0.2 0.1 0.0 <7T|Qe|K> 0.00 0.01 0.02 0.03 0.04 0.05 -0.010 -0.015 - -0.020 - -0.025 0.00 0.01 0.02 0.03 0.04 0.05 Figure 2. The K ^ tt matrix element of the bare Figure 3. The K ^ tt matrix element of the bare operator Qq (charm contribution not included) . operator Q^g^^^ . plicit chiral symmetry breaking effects which do not depend on the quark mass to cancel. From the fit depicted in Fig. ^ wc find that the constant term is zero within errors. Note that this ratio is extremely well resolved, and quite linear. We take this opportunity to quote our value for the kaon B parameter, B'j^'^\2 GeV) = 0.538(8), with Zll/Z\ = 0.928(6) computed nonperturba- tively in the RI scheme and matched to the MS- NDR scheme 0. The errors are statistical only, and the error on Bk is obtained by adding the er- rors on the matrix element and the renormaliza- tion factor in quadrature. Our value is lower than the one quoted by Taniguchi at this meeting||], probably due to the fact that our renormalization constant is lower than the pcrturbative one used in that study. Note our result is for f3 — 6.0 with a single lattice size. 0.00 -0.10 0.00 0.02 0.04 0.06 Figure 4. The K ^ matrix element of the bare operator Qq (charm contribution not included). REFERENCES 1. A. Ah Kh an, lat/00070ll T. et al. Blum, (CP-PACS), et al. (RBC), hcp- hep- lat/0007038 , and references therein. 2. C. Dawson, in these proceedings. 3. G. MartincUi, in these proceedings. 4. C. Bernard, et at, PRD32 (1985) 2343. 5. R. D. Mawhinney, in these proceedings. 6. J. Noaki, in these proceedings. 7. M. Crisafulh, et al, PLB369 (1996) 325. 8. Y. Taniguchi, in these proceedings.