$K\toππ$ Decays with Domain Wall Fermions: Lattice Matrix Elements

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.