# Kaon semileptonic decays near the physical point

###### Abstract

The CKM matrix element can be extracted from the experimental measurement of semileptonic decays. The determination depends on theory input for the corresponding vector form factor in QCD. We present a preliminary update on our efforts to compute it in lattice QCD using domain wall fermions for several lattice spacings and with a lightest pion mass of about . By using partially twisted boundary conditions we avoid systematic errors associated with an interpolation of the form factor in momentum-transfer, while simulated pion masses near the physical point reduce the systematic error due to the chiral extrapolation.

Kaon semileptonic decays near the physical point

RBC and UKQCD Collaborations

\abstract@cs

The 30th International Symposium on Lattice Field Theory, June 24 - 29, 2012, Cairns, Australia

## 1 Introduction

In the Standard Model (SM) the unitary Cabibbo-Kobayashi-Maskawa (CKM) matrix contains information on the strength of flavour-changing weak decays as well as information on correlations between different such processes. Inconsistencies in the CKM-picture would indicate the presence of new physics beyond the SM. One therefore tries to determine all CKM-matrix elements as precisely as possible by studying flavour changing processes both experimentally (e.g. at NA62 and LHCb at CERN) and theoretically. Here we concentrate on the determination of the matrix element from the study of semileptonic kaon () decays which enables a test of first-row CKM-unitarity, . is neglected because it is smaller than the errors in the other two terms on the l.h.s., and . is known very precisely from neutron -decay and matching the corresponding precision for is crucial in searching for deviations from CKM-unitarity and for possible signs of new physics.

## 2 Kaon semileptonic decays

To date, one of the most precise determinations of comes from semileptonic decays (cf. [4]): The product can be determined with a precision at the per mil level from the experimental decay rate [5] and lattice QCD provides the value for . Currently the level of precision set by experiment sets the precision goal for lattice simulations.

The form factor is defined from the vector part of the strangeness-changing weak current () according to

(2.0) |

where is the momentum transfer. We define the scalar form factor

(2.0) |

with . In the SU(3) flavour limit where , by vector current conservation. It can be expanded in terms of the meson masses as where [6]. is a known function of meson masses and the SU(3) pseudoscalar decay constant. In lattice simulations therefore effectively only the small higher order correction is computed and extrapolated to the physical point.

We compute the matrix element in (2) in terms of the ground-state contribution to suitable ratios of Euclidean two- and three-point functions. In a finite lattice box with periodic boundary conditions for the quark fields, the matrix element can in this way only be computed for meson momenta corresponding to the Fourier modes, i.e. with for a spatial volume . The form factor at is then computed by interpolating between the data for the form factor at these Fourier-points [1] thereby introducing a dependence of the final result on the interpolation model; different ansätze may lead to different results for the form factor, introducing systematic uncertanities.

## 3 Partially twisted boundary condition

In contrast to periodic boundary conditions, twisted boundary conditions allow access to hadron momenta other than [8]. Here we employ partially twisted boundary conditions [9] by applying the twist only to the valence quarks,

(3.0) |

( is the twist angle in the -direction). In our simulation a charged meson of mass with one of the valence quarks twisted with angle then obeys the dispersion relation

(3.0) |

up to exponentially supressed finite volume corrections [9]. For matrix elements like , where the initial and final state mesons carry twists and , respectively, the momentum transfer can be written as [10]

(3.0) |

By adjusting the twists on the initial and final meson (here ), we can evaluate directly at . The commonly used implementation is to contruct ratios of correlation functions, (see e.g. [10]), where either the kaon or pion is kept at rest and where the twist on the other meson is tuned in order to obtain results [10]:

(3.0) |

The form factor can then be evaluated directly at as

(3.0) |

This expression is derived by considering only the time-component of the vector current matrix element (2). Using all the other components of the weak vector current, we can construct an over-constrained system of equations

(3.0) |

which we solve to obtain the form factors directly at . Note that distributing the twist over all spatial directions maximises the number of non-trivial equations. Having eliminated the systematic error due to the interpolation, the remaining dominant source of systematic uncertainty is due to the chiral extrapolation of the lattice data to the physical point [2].

## 4 Results

In [1] we presented results for the form factor from simulations of domain wall fermions (DWF) with the Iwasaki Gauge action, an inverse lattice spacing of 1.7GeV and in the range –. Here we extend the analysis by adding new ensembles 32Fine [11], 32Coarse [12] (cf. table 1) providing data for considerably lighter pion masses down to and two additional lattice spacings.

Label | Size | |||
---|---|---|---|---|

24Coarse | Iwasaki | 1.75(4) | , , , | |

32Fine | Iwasaki | 2.31(4) | , , | |

32Coarse | Iwasaki+DSDR | 1.37(1) | , |

Closer to the physical point the twist-angles prescribed by eqs. (3) are very large, in particular for the case where the kaon is moving and the pion is at rest. This manifests itself in larger statistical fluctuations for the ratio . In order to find a better choice of kinematics we analysed the situation further: From eq. (3), the slope of with respect to is given by

(4.0) |

for the matrix element of the time-component of the vector current, , and

(4.0) |

for the spatial components .

The solutions for all equations for () and () on the 32Coarse ensemble with are shown in the l.h.s. and r.h.s. plots in Fig 1, respectively. While all solutions have a negative slope for the case where only the kaon is twisted, there are solutions with opposite slopes in the case where the pion is twisted. Because the solution is given by the intersection of the individual constraints, the kinematical situation where the pion is moving (twisted) and the kaon is at rest provides the best result. The statistical errors are also smaller in this case. Motivated by these observations, we computed all correlation functions once again for a third choice of kinematics with , where both the kaon and the pion are twisted. This leads to a good constraint for as shown in the left plot in Fig. 2. The result is in agreement with the result obtained by solving all simultaneous equations for the cases where either the pion or the kaon are twisted as shown in the right plot of Fig. 2 (obtained by combining the plots in Fig. 1).

We now turn our attention to the mass- and momentum-dependence of the results. In principle, using partially twisted boundary conditions, one is independent of the momentum interpolation. We wish however to include our earlier data sets, the three heavier ensembles in the 24Coarse set, for which we have not generated data directly at . From [1] we know that in these cases the interpolation introduces hardly any model-dependence. In addition to the momentum-dependence a fit ansatz should also incorporate the -symmetry-breaking nature of and the strange quark mass dependence since the simulated strange quark mass is not exactly at the physical point. Our ansatz with four fit parameters [1] is

(4.0) |

Fig 3 (left) summarizes the results of the analysis for all ensembles. The data points are the ones for the simulated, i.e. unphysical strange-quark mass. After correcting towards the physical strange-quark mass using the ansatz in eq. (4.0), all data points line up on the fit-curve in the r.h.s. plot in Fig. 3. Our preliminary result for at the physical point is indicated by the light blue square. At this early stage of the analysis we find that the statistical error at the physical point has been reduced with respect to our earlier results [2].

In our fit ansatz we use the decay constant in the chiral limit, , in the NLO term and have a form for the NNLO term which together make the ansatz consistent with the Ademollo-Gatto theorem and symmetry under interchange of and . Since we do not know the precise value for , we repeat the global fit for and . This variation in is also assumed to account partially for our NNLO term not being the full chiral perturbation theory expression. The result for using is taken as our best value and the values at and provide an estimate of the systematic error from the chiral extrapolation.

In Fig. 4 we compare our new preliminary result for with other recent Lattice determinations. We emphasize that our result is preliminary, but we expect to reach a new level of precision when the analysis is complete.

## 5 Conclusion

We have presented an update on our simulations aiming at the precise prediction of the kaon semileptonic form factor at vanishing momentum transfer. The novelties are simulations for pion masses down to and for three different lattice spacings. We identified preferred choices for the kinematics when using partially twisted boundary conditions in order to simulate directly at . In particular, we now understand why twisting either only the pion or both the kaon and the pion can constrain the form factor better than twisting only the kaon. At this early stage of the analysis the inclusion of the new ensembles in our global fit leads to a reduction of the statistical error at the physcial point compared to our earlier result [2]. This indicates significant progress in the study of K semileptonic form factors which would directly impact our knowledge of the CKM Matrix and provide for improved constraints for new physics beyond the Standard Model.

## 6 Acknowledgement

We thank our colleagues in RBC and UKQCD within whose programme this calculation was performed. The computations were done using the STFC’s DiRAC facility at Swansea, JUGENE at the Jülich Supercomputing Centre and the DiRAC facility at Edinburgh. KS acknowledges support by the European Union under the Grant Agreement number 238353 (ITN STRONGnet); CTS and JMF by STFC grant ST/J000396/1 and ST/H008888/1; P.A.B, by STFC grants ST/K000411/1, ST/J000329/1 and ST/H008845/1; AJ by European Research Council (FP7/2007-2013)/Grant agreement 27975; and JMZ by the Australian Research Council under grant FT100100005.

## References

- [1] P. A. Boyle et al. [RBC/UKQCD], Phys. Rev. Lett. 100 (2008) 141601, [arXiv:0710.5136].
- [2] P. A. Boyle et al. [RBC/UKQCD], Eur. Phys. J. C 69 (2010) 159 [arXiv:1004.0886 [hep-lat]].
- [3] R. Arthur et al. [RBC/UKQCD], arXiv:1208.4412 [hep-lat].
- [4] G. Colangelo et al. [FLAG], Eur.Phys.J. C71 1695 (2011).
- [5] M. Antonelli et al., Eur.Phys.J. C69 399-424 (2010).
- [6] J. Gasser and H. Leutwyler, Nucl. Phys. B250 (1985) 517-538.
- [7] H. Leutwyler and M. Roos, Z. Phys. C 25, 91 (1984).
- [8] P. F. Bedaque and J.-W. Chen, Phys. Lett. B616 (2005) 208–214, [hep-lat/0412023].
- [9] C. T. Sachrajda and G. Villadoro, Phys. Lett. B609 (2005) 73–85, [hep-lat/0411033].
- [10] P. A. Boyle et al. [RBC/UKQCD], JHEP 05 (2007) 016, [hep-lat/0703005].
- [11] Y. Aoki et al. [RBC/UKQCD], Phys. Rev. D83 (2001) 074508.
- [12] R. Arthur et al. [RBC/UKQCD], arXiv:1208.4412.
- [13] M. Ademollo and R. Gatto, Phys. Rev. Lett. 13, 264 (1964).