Closing in on the radiative weak chiral couplings

We point out that, given the current experimental status of radiative kaon decays, a subclass of the O(p4)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{O} (p^4)$$\end{document} counterterms of the weak chiral lagrangian can be determined in closed form. This involves in a decisive way the decay K±→π±π0l+l-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K^\pm \rightarrow \pi ^\pm \pi ^0 l^+ l^-$$\end{document}, currently being measured at CERN by the NA48/2 and NA62 collaborations. We show that consistency with other radiative kaon decay measurements leads to a rather clean prediction for the O(p4)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{{O}}(p^4)$$\end{document} weak couplings entering this decay mode. This results in a characteristic pattern for the interference Dalitz plot, susceptible to be tested already with the limited statistics available at NA48/2. We also provide the first analysis of KS→π+π-γ∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_S\rightarrow \pi ^+\pi ^-\gamma ^*$$\end{document}, which will be measured by LHCb and will help reduce (together with the related KL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K_L$$\end{document} decay) the experimental uncertainty on the radiative weak chiral couplings. A precise experimental determination of the O(p4)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{{O}}(p^4)$$\end{document} weak couplings is important in order to assess the validity of the existing theoretical models in a conclusive way. We briefly comment on the current theoretical situation and discuss the merits of the different theoretical approaches.


Introduction
Radiative kaon decays have been, and still are, one of the most powerful tools we have to explore the physics behind S = 1 processes, most prominently by performing chiral tests and studies on CP violation. Chiral tests are a direct probe of the Standard Model at low energies. In this regime, physics is described by chiral perturbation theory (ChPT), where nonperturbative effects are parametrized in terms of a set of low-energy couplings. Keeping only the octet contributions, the S = 1 Lagrangian up to NLO order is given by a e-mail: oscar.cata@uni-siegen.de where u 2 = U is the chiral matrix containing pions and kaons, is the electric charge matrix and = uλ 6 u † is the chiral spurion for the S = 1 transitions. The second term collects the NLO operators W j . Out of the 37 operators present, only combinations of W 14 , ..., W 18 and W 28 , ..., W 31 are relevant for radiative kaon decays.
As opposed to their analogs in the bottom and charm sector, radiative kaon decays are generically dominated by Bremsstrahlung, while resonance effects are suppressed. This makes the extraction of the weak counterterms quite challenging. The best strategy to determine the chiral couplings is to study the interference term between the Bremsstrahlung and resonance contributions, which is substantially bigger than the resonance piece. Since the interference term is linear in the counterterms, this has the additional advantage that one is sensitive to both their magnitude and sign.
This strategy has been applied, e.g., to the decay modes K ± → π ± γ * , K S → π 0 γ * , K ± → π ± π 0 γ or K ± → π ± γ γ . This way three independent combinations of counterterms are known. In order to fully determine the set of weak counterterms N 14 , ..., N 18 an extra decay mode has to be measured. From this perspective, K + → π + π 0 γ * stands out as the most promising candidate: besides counterterm combinations already present in K ± → π ± γ * and K ± → π ± π 0 γ , it contains an additional combination. Furthermore, it has already been measured by NA48/2, which collected of the order of 5000 events thereof, statistics that will be improved by the new data currently being collected at NA62.
The main observation of this note is that the undetermined counterterm in K + → π + π 0 γ * turns out to be sizeable. Our observation is not based on model estimates, but arises from combining the different experimental results on radiative kaon decays. If this expectation on the size of the counterterm is combined with the strategy laid out in [1,2], which shows how cuts in the dilepton invariant mass can diminish the dominance of Bremsstrahlung effects, a measurement could be performed already with the NA48/2 data. In this note we provide a detailed analysis of the interference term to aid the counterterm determination.
In the last part of the paper we comment on the decay mode K S → π + π − e + e − . In spite of the challenges that it poses, there are prospects to measure it at LHCb in the near future. For this decay mode the relevant counterterm combinations turn out to be entirely predicted from the already measured radiative kaon decay modes, so that one can have a rather precise estimate of how the interference term should look. LHCb could therefore provide a very relevant consistency check of the counterterm structure, and potentially improve the precision of the weak chiral couplings.
We conclude this note with some remarks on the present theoretical understanding of the radiative weak chiral couplings.

Experimental status of the radiative chiral weak counterterms
The radiative kaon decay modes K ± → π ± γ * , K S → π 0 γ * , K ± → π ± π 0 γ and K + → π + γ γ have been observed and measured by the NA48 collaboration over the last 15 years. Besides the determination of the branching ratios, an extraction of the relevant combinations of chiral weak counterterms for those decays has also been possible. The chiral counterterms are related to the slope of the differential decay rate, which can be most easily accessed through the interference term between Bremsstrahlung and electric emission. Using this technique, the experimental values for the slopes of the different modes were found to be K ± → π ± γ * : a + = −0.578 ± 0.016 [3,4], K + → π + γ γ :ĉ = 1.56 ± 0.23 ± 0.11 [8].
The results are linked to the different weak counterterm combinations through (see [9][10][11][12]) where |M K | = 1.81 × 10 −8 GeV is the amplitude for K + → π + π 0 , f π = 93 MeV is the pion decay constant and G 8 = 9.1 × 10 −6 GeV −2 . The results for the counterterms are summarized in Table 1, where we have used the values taken from [13] and [14], respectively. We remark that while N (0) E and N 0 are scale-independent, N (1) E and N S depend on scale. All quantities have been evaluated at the kaon mass scale.
3 Prospects for K + → π + π 0 e + e − NA48/2 has collected roughly 5 × 10 3 K + → π + π 0 e + e − events [15]. With this rather limited statistics one is experimentally sensitive to the Bremsstrahlung and magnetic terms, which are known from K + → π + π 0 (through Low's theorem) and from K + → π + π 0 γ , respectively. Whether one can also extract information on the (much suppressed) electric counterterms depends on the size of the interference term. Besides N (0) E and N (1) E , defined in the previous section, K + → π + π 0 e + e − also contains the counterterm combination which is so far undetermined. The interplay of the three counterterms will determine whether NA48/2 is sensitive to the chiral counterterms. In this section we will show that N (2) E is expected to be positive and sizeable. This expectation is mostly determined by consistency with the experimental values of the previous section and is therefore rather robust. We will then show how N (0) E and the estimated N E define a rather characteristic pattern for the interference Dalitz plot, which can be tested by using cuts on the invariant dilepton mass.
The decay K + → π + π 0 e + e − was first studied in [1,2]. The amplitude for this process can be factorized into leptonic and hadronic pieces as where The last form factor contains the magnetic emission, already measured in K + → π + π 0 γ [7]. An independent measurement through K + → π + π 0 e + e − would be interesting, and a detailed analysis on how to do that was already presented in [2]. Here we are mostly concerned with the weak chiral couplings associated with the electric emission. They enter the form factors as where δ 2 0 and δ 1 1 are strong phases associated with final state interactions and M K is the amplitude for the K + → π + π 0 decay. The counterterm combinations N depend on scale. However, in both cases this dependence turns out to be very mild and well within the experimental errors. We will accordingly omit it, as was done in [2].
It should be noted that Eq. (14) differs from a NLO computation using chiral perturbation theory (see e.g. [1]): the Bremsstrahlung contribution is here evaluated using Low's theorem, which is an exact current algebra result, whereas the weak and strong countertems are evaluated at the classical level, neglecting loop corrections. Since M K is taken from experiment, we are including the effects of final state interactions (pion rescattering), which are known to be important.
The decay rate is overwhelmingly dominated by the Bremsstrahlung component. However, as initially pointed out in [1], this dominance decreases as the dilepton invariant energy q 2 is increased. Judicious cuts in q 2 can therefore allow experimental analyses to reach information on the weak counterterms. This strategy was developed in detail in [2], suggesting an optimal cut at q c ∼ 50 MeV. 1 The analysis done in [2] was mostly meant for a first experimental analysis of K + → π + π 0 e + e − , with the determination of the Bremsstrahlung contribution and the magnetic piece as main goals (see [15]). At the time, the electric contribution seemed inaccessible without NA62 data. Accordingly, sample values were chosen for the counterterms, mostly to show the shape of the Dalitz plot and get an order of magnitude estimate for the electric contributions. Presently there is an active effort to study the electric counterterms with NA48/2 data. This motivates to refine the analysis presented in [2].
The relevant piece for the extraction of the weak chiral couplings is the interference term. Integrating over the angular variables, the differential decay rate can be parametrized as where δ = δ 1 0 − δ 2 0 π/18 and |M K | = 1.81 × 10 −8 GeV. As in [2], we have chosen to work with the set of variables (E γ ,T c , q), originally defined in [16], with the redefinition T c = T c + m π . The kinematic functions f i are given by Fig. 1 Dalitz plots for the interference differential decay rate of K + → π + π 0 e + e − in the (E γ , T c ) plane for q = 20 MeV (left panel) and q = 50 MeV (right panel). Numbers are given in units of 10 −20 GeV −1 .
The contour plot is 'spikier' the lower the q values are, a pattern mostly dictated by the structure of the Bremsstrahlung term Values for the counterterms N E are taken from experiment (see the previous section). 2 Using the experimen-2 Ref. [2] used N (0) E = −0.0022. The experimental value bears a positive sign, which we take in this note. All tables and plots in [2] were consistently done with the negative choice for N (0) E . tal information of Table 1, one gets the following prediction which is entirely based on experimental input. The value and sign of N 17 is unknown and can only be determined from an actual measurement of N (2) E . It is expected that N 17 be small (a number of theoretical models predict N 17 to be actually zero). In the following we will assume that N 17 is small enough that it can be neglected in Eq. (17). Given the experimental error quoted in (17), this is true if |N 17 | 0.002.
In Fig. 1 we show the Dalitz plot for the interference differential decay rate in the (E γ , T c ) plane for different values of the invariant dilepton mass q. As expected, the plots have a smoother behavior the larger q is, since the Bremsstrahlung peak is located at small q. In Table 2 we provide the values of the interference term (compared to the Bremsstrahlung) for different choices of the counterterms for different cuts q c . We use the shorthand notation where a subscript (1, 1, 1) indicates that all counterterms take on their estimated values. When a counterterm is neglected, it is indicated by a nul value. Given the definition of the variable q c , the first row in Table 2 corresponds to the fully integrated rate in q 2 , while the subsequent rows cut the lower q 2 regions. Notice that N (0) E determines the overall sign of the interference term, while N (2) E is important to revert this sign already at low q c values. N (1) E has a marginal effect. It is interesting to note that a sizeable N (2) E tends to make the interference term smaller, though not dramatically. As the table shows, different patterns for the counterterm values lead to distinct behaviors of the interference term with q c . This could be used to extract a value for N (2) E or, more generally, to fit for In Fig. 2 we plot the results of the last four columns of Table 2, including also the error band, estimated by the 1 sigma shift of the experimental input.
4 Long-distance contributions to K S → π + π − e + e − The branching ratio for K S → π + π − e + e − has already been measured by the NA48 collaboration [17,18]. In order to probe the chiral O( p 4 ) structure of this decay one needs to reach the percent precision, which might be possible at LHCb [19]. Here we provide a first analysis of this decay mode, which bears parallelisms with K L → π + π − e + e − , already studied in [1,20].
A particularly interesting aspect of K S → π + π − e + e − is that the magnetic piece is a CP-violating effect. This means that the interference term is, despite being suppressed, the relevant contribution after the Bremsstrahlung. The charge radius contribution can also be shown to be CP-violating. Thus, to a very good approximation, one finds that the hadronic form factors can be described by where ∼ 10 −3 e iπ/2 parametrizes (indirect) CP violation and M K is the matrix element of the K S → π + π − decay, which can be written as Experimentally, one finds that |M K | = 3.92 × 10 −7 GeV. The first term in both F 1 and F 2 is the Bremsstrahlung piece. In analogy with the analysis of the previous section, we evaluate it using Low's theorem, such that the important effect of final state interactions is taken into account. The remaining pieces in F 1 and F 2 collect the strong and weak  17 , which was already present in K + → π + π 0 γ * , and defined F 1,2 also contain the combination N N 17 ), which is -suppressed and has been omitted in (19). The direct emission is estimated with the O( p 4 ) strong and weak counterterms, neglecting the loop corrections.
It is instructive to compare the results in Eqs. (19) with the corresponding expressions for K L → π + π − e + e − given in [1,20]: • In Eq. (19) the Bremsstrahlung piece is the dominant contribution while the magnetic piece is a CP-violating effect. For K L one finds the reverse situation, which makes this latter decay mode especially suited to determine the magnetic counterterm N 19 + N 31 combination.

• An interesting feature of Eq. (19) is the SU (3)-violating
pole structure that comes with L 9 . In K L such a structure is also present but violates CP and is completely negligible. Such a term is also present in K + → π + π 0 γ * , but there it represents a tiny O(m 2 π + −m 2 π 0 ) isospin violation, which we neglected altogether in the previous section. In contrast, for K S → π + π − γ * it is a non-negligible contribution.
• The weak electric counterterm combinations are the same in both decay modes. However, N E exchange their roles: for K S the impact of N (4) E is negligible, while for K L it is N (3) E that can be safely dismissed.
The dominance of N (3) E over N (4) E is relevant. While the latter combination is presently unknown, the former can be actually predicted from the experimental values in Table 1 to be where the uncertainty is purely experimental, in terms of the measurements of K ± → π ± γ * and K ± → π ± π 0 γ . We remark that the previous result does not depend on any assumptions about N 17 . One therefore finds that the longdistance contributions to K S → π + π − γ * can be determined with remarkable accuracy, namely where in the interference term, by virtue of the I = 1 2 rule, the relevant phase difference is to a very good approximation given by δ 0 0 − δ 1 1 34 0 [21]. This number is in excellent agreement with the PDG average [22]: Similarly, one can predict that Brems.  Fig. 4 Left panel: values of N 14 and N 15 as given by K ± → π ± γ * (blue band) and K S → π 0 γ * (violet band). Right panel: values for N 16 and N 17 extracted from K ± → π ± π 0 γ (blue band) and K ± → π ± π 0 e + e − (yellow band) measurements. The latter is an educated estimate (see main text) Using the same formalism employed in the previous section for K + → π + π 0 γ * , in Fig. 3 we show the Dalitz plot for the interference differential decay rate of the decay K S → π + π − e + e − in the (E γ , T c ) plane for fixed dilepton invariant masses q = 20 MeV and q = 50 MeV. Note that, in contrast to K + → π + π 0 e + e − , the interference term is positive all over the physical region.

Concluding remarks
The experimental situation in radiative kaon decays is at a point where the relevant (electric) weak chiral counterterms N 14 , ..., N 18 can be individually determined. The outcome of the measurements on K ± → π ± γ * and K S → π 0 γ * is plotted in Fig. 4 and gives Adding the measurements on K ± → π ± γ γ , one can also conclude that The determination of N 16 and N 17 requires information from K ± → π ± π 0 γ and an extra radiative decay mode, either K ± → π ± π 0 γ * or K L → π + π − γ * . So far only the combination N 16 + N 17 is known. K ± → π ± π 0 γ * is currently being analyzed by the NA48/2 collaboration, which gathered 5000 events thereof. With this rather limited statistics, an extraction of N 16 and N 17 is only feasible if (a) there is some expectation for the size of the counterterm combinations to be tested; and (b) a strategy is devised to compensate the overwhelming dominance of the Bremsstrahlung contributions.
In this paper we have provided an estimate for the counterterms to be measured in K ± → π ± π 0 e + e − . Our estimates are, to the extent possible, based on the existing experimental knowledge of radiative kaon decays, and only very mild theoretical assumptions are used. We have shown that they generate a distinct shape in the Dalitz plot and, if judicious cuts on the dilepton invariant mass are used, one can largely compensate the relatively low statistics currently available. While the larger statistics of NA62 will definitely allow for a more precise determination of the counterterm combinations, the strategy proposed here is meant to make an extraction with NA48/2 data possible.
On the right panel of Fig. 4 we show estimates for N 16 and N 17 , where the external yellow band corresponds to a value of N (2) E = 0.089 (11), obtained assuming that N 17 is negligible. The internal yellow band corresponds to the same central value but assuming that a precision of 4 × 10 −3 can be reached. The experimental determination would fall out of the yellow band if N 17 turns out to be sizeable, roughly |N 17 | 2 × 10 −3 .
While K ± → π ± π 0 γ * stands out as the most promising mode for an extraction of N 17 (together with the already existing information on K ± → π ± π 0 γ ), one should also bear in mind that the weak counterterms in K L ,S → π + π − γ * satisfy the following relation: Therefore, N 17 can be accessed through both charged and neutral radiative kaon decays. Given the kaon program at LHCb, this latter possibility turns out to be of interest. In this paper we have performed a study of K S → π + π − γ * which, together with the existing studies of K L → π + π − γ * , should make an extraction of N 17 from neutral radiative kaon decays feasible.
We will close this note with some comments on the theoretical understanding of the weak counterterm values, i.e., on the electroweak interactions at low energies. The situation is not as solid as for the strong sector, where the bulk of the nonperturbative effects entering the Gasser-Leutwyler coefficients L i is due to resonance exchange [23], with vector meson dominance (VDM) as a solid guiding principle. The weak counterterms, as opposed to the strong ones, are sensitive to the whole range of energies. Resonance models with VDM are therefore based on the assumption that the low-energy region is dominant. Given the large number of parameters to fit, additional simplifying assumptions are used in order to end up with predictive schemes. Failure to reproduce the experimental numbers can therefore be attributed to (a) a non-negligible contribution from short distances inside the weak counterterms; and (b) assumptions on the long-distance contributions which are not supported phenomenologically.
While the predictions of the different models should be taken with caution, they have merits that sometimes are not fully appreciated. In table 3 we have listed a number of counterterm combinations relevant for different radiative kaon decays (see [24] for a more comprehensive list), together with the predictions of a number of models: the weak deformation model (WDM) [25], factorization model (FM) [26], holographic electroweak model (HEW) [27] and the resonance model studied in Ref. [28].
The WDM/FM/HEW models express their predictions in terms of the strong counterterms L 9 , L 10 and H 1 . Leaving aside H 1 , which has not been determined, the experimental values for L 9 and L 10 indicate that The counterterm combinations listed on the left have been actually measured (see Table 1) and agree with the experimental numbers. Notice that the O(10 −3 ) predictions are interpreted as accidental cancellations, proportional to (L 9 + L 10 ), which is known to be a small number. Guessing the resulting sign when such cancellations occur is stretching the models beyond their range of applicability. 3 The models also predict that the undetermined counterterms in K + → π + π 0 e + e − and K S → π + π − e + e − discussed in this paper should be large and positive, in agreement with the results of Eqs. (17) and (22). In contrast to the WDM model, the resonance model of ref. [28] allows one to gauge the relative weight of vector and axial resonances as long-distance components of the weak chiral couplings and thereby test how robust is the VDM assumption. From Table 3 it is clear that if VMD is assumed, then there is a reasonable qualitative agreement with the conclusions of Eq. (29). In particular, notice that the cancellations of the form (L 9 + L 10 ) in the third column of Table 3 are interpreted as vectorial cancellations in the fourth column. Notice also that the counterterms enter-ing K L ,S → π + π − γ * are predicted to be the same in the third column but have different axial-vector contributions in the fourth column. A determination of N 16 and N 17 through K + → π + π 0 γ * should clarify the importance of axialvectors in those counterterm combinations.