Lattice QCD calculation of $Kto ell u_ell ell^+ ell^-$ decay width

We develop a methodology for the computation of the $Kto ell u_ell ell^+ ell^-$ decay width using lattice QCD and present an exploratory study here. We use a scalar function method to account for the momentum dependence of the decay amplitude and adopt the infinite volume reconstruction (IVR) method to reduce the systematic errors such as the temporal truncation effects and the finite-volume effects. We then perform a four-body phase-space integral to obtain the decay width. The only remaining technical problem is the possible power-law finite-volume effects associated with the process of $Ktopipi ell u_ellto ell u_ell ell^+ ell^-$, where the intermediate state involves multiple hadrons. In this work, we use a gauge ensemble of twisted mass fermion with a pion mass $m_pi=352$ MeV and a nearly-physical kaon mass. At this kinematics, the $pipi$ in the intermediate state cannot be on shell simultaneously as $2m_pi>m_K$ and the finite-volume effects associate with $pipi$ state are exponentially suppressed. Using the developed methods mentioned above, we calculate the branching ratios for four channels of $Kto ell u_ellell^+ ell^-$, and obtain the results close to the experimental measurements and ChPT predictions. Our work demonstrates the capability of lattice QCD to improve Standard Model prediction in $Kto ell u_ell ell^+ ell^-$ decay width.