Calculation of $K \to \pi\pi$ decay amplitudes with improved Wilson fermion action in lattice QCD

TL;DR

This paper calculates $K o\pi\pi$ decay amplitudes using improved Wilson fermions in lattice QCD, demonstrating that operator mixing issues are manageable and providing initial results at unphysical quark masses.

Contribution

It shows that Wilson fermions can be used for $K o\pi\pi$ decay calculations without operator mixing complications, and provides the first such results with this approach.

Findings

Operator mixing with wrong chirality is absent due to CPS symmetry.

First calculation of $K o\pi\pi$ decay amplitudes with Wilson fermions.

Obtained decay amplitudes at unphysical quark masses with controlled uncertainties.

Abstract

We present our result for the $K\to\pi\pi$ decay amplitudes for both the $\Delta I=1/2$ and $3/2$ processes with the improved Wilson fermion action. Expanding on the earlier works by Bernard {\it et al.} and by Donini {\it et al.}, we show that mixings with four-fermion operators with wrong chirality are absent even for the Wilson fermion action for the parity odd process in both channels due to CPS symmetry. Therefore, after subtraction of an effect from the lower dimensional operator, a calculation of the decay amplitudes is possible without complications from operators with wrong chirality, as for the case with chirally symmetric lattice actions. As a first step to verify the possibility of calculations with the Wilson fermion action, we consider the decay amplitudes at an unphysical quark mass $m_K \sim 2 m_\pi$. Our calculations are carried out with $N_f=2+1$ gauge configurations generated with the Iwasaki gauge action and nonperturbatively $O(a)$-improved Wilson fermion action at $a=0.091\,{\rm fm}$, $m_\pi=280\,{\rm MeV}$, and $m_K=580\,{\rm MeV}$ on a $32^3\times 64$ ($La=2.9\,{\rm fm}$) lattice. For the quark loops in the penguin and disconnected contributions in the $I=0$ channel, the combined hopping parameter expansion and truncated solver method work very well for variance reduction. We obtain, for the first time with a Wilson-type fermion action, that ${\rm Re}A_0 = 60(36) \times10^{ -8}\,{\rm GeV}$ and ${\rm Im}A_0 =-67(56) \times10^{-12}\,{\rm GeV}$ for a matching scale $q^* =1/a$. The dependence on the matching scale $q^*$ for these values is weak.