K to ππ decay amplitudes from lattice QCD

We report a direct lattice calculation of the K to ?? decay matrix elements for both the ?I=1/2 and 3/2 amplitudes A0 and A2 on 2+1 flavor, domain wall fermion, 163×32×16 lattices. This is a complete calculation in which all contractions for the required ten, four-quark operators are evaluated, including the disconnected graphs in which no quark line connects the initial kaon and final two-pion states. These lattice operators are nonperturbatively renormalized using the Rome-Southampton method and the quadratic divergences are studied and removed. This is an important but notoriously difficult calculation, requiring high statistics on a large volume. In this paper, we take a major step toward the computation of the physical K??? amplitudes by performing a complete calculation at unphysical kinematics with pions of mass 422 MeV at rest in the kaon rest frame. With this simplification, we are able to resolve Re(A0) from zero for the first time, with a 25% statistical error and can develop and evaluate methods for computing the complete, complex amplitude A0, a calculation central to understanding the ?=1/2 rule and testing the standard model of CP violation in the kaon system.