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 ${\cal O} (p^4)$ counterterms of the weak chiral lagrangian can be determined in closed form. This involves in a decisive way the decay $K^\pm \to \pi ^\pm \pi ^0 l^+ l^-$, 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 ${\cal{O}}(p^4)$ 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 $K_S\to \pi^+\pi^-\gamma^*$, which will be measured by LHCb and will help reduce (together with the related $K_L$ decay) the experimental uncertainty on the radiative weak chiral couplings. A precise experimental determination of the ${\cal{O}}(p^4)$ 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 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 even with the NA48/2 data. In particular, 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 fact that its measurement is rather challenging, 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 measurements on other 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.
2 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] (2) K + → π + γγ :ĉ = 1.56 ± 0.23 ± 0.11 [8] .
Whether one can also extract information on the (much suppressed) electric counterterms depends on the size of the interference term. Besides N 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 rather robust. We will then show how 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 E 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 As in [2], we have chosen to work with the set of variables (E γ ,T c , q), originally defined in [16], with the redefinitionT c = T c + m π .
The kinematic functions f i are given by 1 q c is defined as the lower integration limit of the integral in q. Accordingly, q c = √ 2m e corresponds to performing the whole integral in q, while increasing q c progressively eliminates the low-q region, where the Bremsstrahlung dominates.
Values for the counterterms N are taken from experiment (see the previous section). 2 Using the experimental 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 E . It is expected that N 17 be small (a 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  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 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 at 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 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 E is important to revert this sign at already 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 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 a 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 [20,1].
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 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 1 collect the strong and weak O(p 4 ) corrections. In (19) we have used N (0) 17 , which was already present in K + → π + π 0 γ * , and defined 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 [20,1]: • 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 nonnegligible contribution.
• The weak electric counterterm combinations are the same in both decay modes. However, N E and N E exchange their roles: for K S the impact of N E that can be safely dismissed.
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. One therefore finds that the long-distance contributions to K S → π + π − γ * can be determined with remarkable accuracy, namely BR(K S → π + π − e + e − ) = 4.74 · 10 −5 Brems.  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 BR(K S → π + π − µ + µ − ) = 4.17 · 10 −14 Brems.
Using the same formalism employed in the previous section for K + → π + π 0 γ * , in fig. 3 we show the Dalitz plot for 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 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).
Adding the measurements on K ± → π ± γγ, one can also conclude that The determination of N 16 and N 17 requires information on K ± → π ± π 0 γ and K ± → π ± π 0 γ * . The latter 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 .
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, with vector meson dominance (VDM) as a solid guiding principle. The weak counterterms, as opposed to the strong ones, counterterm combinations 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 nonnegligible 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 [23] for a more comprehensive list), together with the predictions of a number of models: the weak deformation model (WDM) [24], factorization model (FM) [25], holographic electroweak model (HEW) [26] and the resonance model studied in ref. [27].
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 contrast to the WDM model, the resonance model of ref. [27] 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. (28). An exception is the counterterm entering K L → π + π − γ * , which is predicted to have a marginal vectorial component and accordingly should be expected to be small. A determination of N 16 and N 17 through K + → π + π 0 γ * should resolve this issue.