The vector form factor at zero momentum transfer on the lattice
Introduction
The most precise determination of the Cabibbo angle, or equivalently of the CKM matrix element [1], is obtained from () decays. The key observation which allows to reach a good theoretical control on these transitions is the Ademollo–Gatto theorem [2], which states that the form factors, and , at zero four-momentum transfer, are renormalized only by terms of at least second order in the breaking of the SU(3) flavor symmetry. The estimate of these smallish corrections, i.e., of the difference of from unity, is presently the dominant source of theoretical uncertainty in the extraction of .
Chiral perturbation theory (CHPT) provides a natural and powerful tool to analyse the amount of SU(3) (and isospin) breaking due to light quark masses. As shown by Leutwyler and Roos [3], within CHPT one can perform a systematic expansion of the type , where . Because of the Ademollo–Gatto theorem, the first non-trivial term in the chiral expansion, , does not receive contributions of local operators appearing in the effective theory and can be computed unambiguously in terms of , and ( in the case [3]). The problem of estimating can thus be re-expressed as the problem of finding a prediction for This quantity is difficult to be evaluated since it depends on unknown coefficients of chiral operators, with . Using a general parameterization of the SU(3) breaking structure of the pseudoscalar meson wave functions, Leutwyler and Roos estimated . Very recently, Bijnens and Talavera [4] showed that, in principle, the leading contribution to Δf could be constrained by experimental data on the slope and curvature of ; however, the required level of experimental precision is far from the presently available one. For the time being we are therefore left with the Leutwyler–Roos result, and the large scale dependence of the loop calculations [4], [5] seems to indicate that its 0.008 error might well be underestimated [6].
The theoretical error on due to the Leutwyler–Roos estimate of is already comparable with the present experimental uncertainty (see, e.g., Ref. [7]). When the high-statistics results from KLOE and NA48 will be available, this theoretical error will become the dominant source of uncertainty on . Given this situation, it is then highly desirable to obtain independent estimates of at the level (or below). The purpose of the present work is to show that this precision can be achieved using lattice QCD.
The strategy adopted in order to reach the challenging goal of a error, is based on the following three main steps:
- 1.
Evaluation of the scalar form factor at . Applying a method originally proposed in Ref. [8] to investigate heavy-light form factors, we extract from the relation where all mesons are at rest. The double ratio and the kinematical configuration allow to reduce most of the systematic uncertainties and to reach a statistical accuracy on well below 1%.
- 2.
Extrapolation of to . By evaluating the slope of the scalar form factor, we extrapolate from to . We note that in order to obtain at the percent level the precision required for the slope can be much lower, because it is possible to choose values of very close to .
For each set of quark masses we calculate two- and three-point correlation functions of mesons with various momenta in order to study the dependence of both and . The latter turns out to be well determined on the lattice, whereas the former does not. We improve the precision in the extraction of by constructing a new suitable double ratio which provides an accurate determination of the ratio . We will define this ratio in Section 4. Fitting the -dependence of with different functional forms, we finally extrapolate to . The systematic error induced by this extrapolation, which is strongly reduced by the use of small values for , is estimated by the spread of the results obtained with different extrapolation functions.
- 3.
Subtraction of the leading chiral logs and chiral extrapolation. The Ademollo–Gatto theorem holds also within the quenched approximation [9], which has been adopted in this work. The leading chiral corrections to , denoted by where the superscript q refers to the quenched approximation, are finite and can be computed unambiguously in terms of the couplings of the quenched CHPT (qCHPT) Lagrangian [10], [11]. For these reasons, in order to get rid of some of the quenched artifacts we define the quantity where and extrapolate it to the physical kaon and pion masses. The ratio : (i) is finite in the SU(3)-symmetric limit; (ii) does not depend on any subtraction scale; (iii) is free from the dominant quenched chiral logs. We emphasize that the subtraction of in Eq. (3) does not imply necessarily a good convergence of (q)CHPT at order for the meson masses used in our lattice simulations. The aim of this subtraction is to define the quantity Δf in such a way that its chiral expansion starts at order independently of the values of the meson masses. In the presence of sizable local contributions, we expect to have a smooth chiral behavior and to be closer to its unquenched analog than the SU(3)-breaking quantity . Extrapolating the values of to the physical meson masses, we finally obtain where the systematic error does not include an estimate of quenched effects beyond . Our result (4) is in good agreement with the estimate obtained by Leutwyler and Roos in Ref. [3].
The systematic error quoted in Eq. (4) is mainly due to the uncertainties resulting from the functional dependence of the scalar form factor on both and the meson masses. This error can be further reduced by using larger lattice volumes (leading to smaller lattice momenta) as well as smaller meson masses. In our estimate of Δf discretization effects start at and are also proportional to , as the physical SU(3)-breaking effects. In other words, our result is not affected by the whole discretization error on the three-point correlation function, but only by its smaller SU(3)-breaking part. Discretization errors on Δf are estimated to be few percent of the physical term, i.e., well within the systematic uncertainty quoted in Eq. (4). For a more refined estimate of these effects, calculations at different values of the lattice spacing are required. Finally, we stress again that the effects of quenching onto the terms of are not estimated and thus not included in our final systematic error.
Section snippets
Notations and lattice details
The form factors of the weak vector current are defined by where is a Clebsch–Gordan coefficient, equal to 1 () for neutral (charged) kaons. As usual, we express in terms of the so-called scalar form factor, By construction and the differences between and channels are only due to isospin-breaking effects. In the following we shall concentrate
Calculation of
Following a procedure originally proposed in Ref. [8] to study the heavy-light form factors, the scalar form factor has been calculated very efficiently at (i.e., ) from the double ratio of three-point correlation functions with both mesons at rest: When the vector current and the two interpolating fields are separated far enough from each other, the contribution of the ground states
Momentum dependence of the form factors and extrapolation to
In this section we perform the extrapolation of the scalar form factor from to . To this end we need to evaluate the slope of , which in turn means to study the -dependence of the scalar form factor. We stress that, in order to obtain at the percent level, the precision required for the slope can be much lower, since the values of used in our lattice calculations are quite close to (see Table 2). We find indeed that a precision is enough for our purposes.
The
Extraction of Δf
In order to determine the physical value of , we need to extrapolate our results of Table 3 to the physical kaon and pion masses. As discussed in the introduction, the problem of the chiral extrapolation is substantially simplified if we remove the effect of the leading chiral logs. To this purpose we consider the quantity where represents the leading non-local contribution determined by pseudoscalar meson loops within CHPT. This procedure is well defined thanks to the
Estimate of and impact on
The value of Δf given in Eq. (37) is in excellent agreement with the estimate of Δf made by Leutwyler and Roos using a general parameterization of the SU(3) breaking structure of the pseudoscalar-meson wave functions [3]. It must be stressed, however, that our result is the first computation of this quantity using a non-perturbative method based on QCD. Combining our estimate of Δf with the physical value of from Eq. (28), we finally obtain to be compared
Conclusions
We presented a quenched lattice study of the vector form factor at zero-momentum transfer. Our calculation is the first one obtained by using a non-perturbative method based only on QCD, except for the quenched approximation. Our main goal is the determination of the second-order SU(3)-breaking quantity , which is necessary to extract from decays. In order to reach the required level of precision we employed the double ratio method originally proposed in Ref. [8] for the
Acknowledgments
The work of G.I. and F.M. is partially supported by IHP-RTN, EC contract No. HPRN-CT-2002-00311 (EURIDICE).
References (31)
- et al.
Eur. Phys. J. C
(2002) - et al.
Phys. Rev. D
(1992) - et al.
Nucl. Phys. B
(1997) - et al.
Phys. Rev. D
(2001)et al.Nucl. Phys. B (Proc. Suppl.)
(2002) - et al.
Phys. Rev. D
(2002)- et al.
Phys. Lett. B
(1987) - et al.
Phys. Rev. D
(2000) Phys. Rev. Lett.
(1963)et al.Prog. Theor. Phys.
(1973)- et al.
Phys. Rev. Lett.
(1964)
Z. Phys. C
Nucl. Phys. B
Phys. Rev. D
Cited by (139)
Beta decays and non-standard interactions in the LHC era
2013, Progress in Particle and Nuclear PhysicsWhen hadrons become unstable: A novel type of non-analyticity in chiral extrapolations
2011, Physics Letters, Section B: Nuclear, Elementary Particle and High-Energy PhysicsFLAG Review 2021
2022, European Physical Journal CFLAG Review 2021
2021, arXivFLAG Review 2019: Flavour Lattice Averaging Group (FLAG)
2020, European Physical Journal C