On some aspects of isospin breaking in the decay K±→π0π0e±νe(-)\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 ^0 \pi ^0 e^\pm \mathop {\nu _e}\limits ^{_{(-)}}$$\end{document}

Two aspects of isospin breaking in the decay K±→π0π0e±νe(-)\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 ^0 \pi ^0 e^\pm \mathop {\nu _e}\limits ^{_{(-)}}$$\end{document} are studied and discussed. The first addresses the possible influence of the phenomenological description of the unitarity cusp on the extraction of the normalization of the form factor from data. Using the scalar form factor of the pion as a theoretical laboratory, we find that this determination is robust under variations of the phenomenological parameterizations of the form factor. The second aspect concerns the issue of radiative corrections. We compute the radiative corrections to the total decay rate for K±→π0π0e±νe(-)\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 ^0 \pi ^0 e^\pm \mathop {\nu _e}\limits ^{_{(-)}}$$\end{document} in a setting that allows comparison with the way radiative corrections were handled in the channel K±→π+π-e±νe(-)\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 ^+ \pi ^- e^\pm \mathop {\nu _e}\limits ^{_{(-)}}$$\end{document}. We find that, once radiative corrections are included, the normalizations of the form factor as determined experimentally from data in the two decay channels come to a better agreement. The remaining discrepancy can easily be accounted for by other isospin-breaking corrections, mainly those due to the difference between the masses of the up and down quarks.


Introduction
The program of analyzing K 4 decays of the charged kaon conducted by the NA48/2 Collaboration at the CERN SPS has so far been very successful. In the π + π − channel of the electron mode, = e (the decay K ± → π + π − e ± (−) ν e will henceforth be referred to as K +− e4 ), it has led, besides a more precise determination of the corresponding brancha e-mail: bernard@ipno.in2p3.fr b e-mail: descotes@th.u-psud.fr c e-mail: knecht@cpt.univ-mrs.fr ing ratio and hadronic form factors [1], to a very accurate determination of the ππ S-wave scattering lengths a 0 0 and a 2 0 [2,3], which constitutes a stringent test of the QCD prediction obtained within the framework of chiral perturbation theory [4][5][6][7]. These data have also contributed to improvement of the accuracy of the determination of some of the low-energy constants of chiral perturbation theory [8][9][10].
More recently, the results concerning an analysis of the data obtained in the π 0 π 0 channel of the electron mode (the decay K ± → π 0 π 0 e ± (−) ν e will henceforth be referred to as K 00 e4 ) have also become available [11]. Although the number of events is lower (∼ 6.5 · 10 4 events in the K 00 e4 mode vs. ∼ 10 6 events for K +− e4 ), this allows for some cross checks at the level of the structure of one of the form factors, that is identical for the two channels in the isospin limit. The normalization of this common form factor, as measured in the two channels, reads [1,11] |V us | f s [ Ignoring for the time being the correction factor δ E M (we will discuss radiative corrections below), the difference of the two values, as compared to the value measured in the K +− e4 channel, amounts to 6.5 % in relative terms. This might be considered as a small difference, but given the uncertainties, it is, in statistical terms, quite significant. Adding all errors in quadrature 1 this gives = 1.065 ± 0.010. (1.2) It seems difficult to ascribe a variation of 6.5 % to the radiative correction factor δ E M alone. While in some regions of phase space radiative corrections can reach the ±10 % level, they usually sum up to ±1 % in the decay rate. The radiative corrections to the K +− e4 decay mode have been discussed in several places [12][13][14] at one-loop precision in the lowenergy expansion. But no comparable study has been done for the K 00 e4 decay mode. There exists an older, less systematic, analysis [15] that covers the corrections due to virtual photon exchanges and real photon emission, which could provide the relevant contributions at a first stage, but its practical use is somewhat limited, since the expressions given there are not always very explicit, and, moreover, need to be checked. Furthermore, not all radiative corrections occurring in the charged K +− e4 channel [12][13][14] have been taken into account in the analysis of the experimental data. These additional radiative corrections could affect f s [K 00 e4 ] and f s [K +− e4 ] in different ways, and make up for another part of the discrepancy.
If δ E M alone does not explain the discrepancy (1.2), one has to look for other sources of isospin-breaking effects. These can be due to the difference between the up and down quark masses m u and m d , conveniently described by the parameter R, with 1/R = (m d − m u )/(m s − m), where m s is the mass of the strange quark, whereas m denotes the average mass of the up and down quarks, m = (m u + m d )/2. For instance, at lowest order in the chiral expansion, one has [12,16]  Finally, there are also isospin-breaking effects induced by the mass difference between charged and neutral pions. Most notable from this point of view is the presence of a unitarity cusp [11] in the form factor describing the amplitude of the K 00 e4 mode. The interpretation of this cusp is by now well understood, and as in the case of the K ± → π 0 π 0 π ± decay [18][19][20][21][22], it arises from the contribution of a π + π − intermediate state in the unitarity sum [for a general discussion of the properties of the K e4 form factors from the point of Footnote 1 continued The ratio in Eq. (1.2), however, does not depend on τ K ± anymore, which lowers the contribution of the "external" uncertainties to Eq. (1.2). At the level of precision shown, this does not impinge on the uncertainty in Eq. (1.2). We are indebted to B. Bloch-Devaux for drawing our attention to this point. view of analyticity and unitarity, see Ref. [23] and references therein].
This cusp contains information on the combination a 0 0 −a 2 0 that describes the amplitude for the process π + π − → π 0 π 0 at threshold. Although this information probably cannot be extracted from the K 00 e4 data in a way as statistically significant as the determination from the cusp in K ± → π 0 π 0 π ± [24], it is nevertheless important to include a correct description of this cusp in the parameterization of the form factor used to analyze the data. This necessity has been demonstrated in full detail in the case of the K ± → π 0 π 0 π ± decay, and it is to be expected that the same attention to these matters should be paid also in the analysis of the K 00 e4 data. Failure to do so may introduce a systematic bias which would make the comparison with the information on the form factor extracted from the K +− e4 data spurious to some extent.
It is the purpose of the present note to address some of these issues. In a first step, we investigate the possible influence that various parameterizations of the form factors could have on the outcome of the analysis. In order to control inputs and outputs fully, we choose to work with a simplified model, where the exact form factors are known from a theoretical point of view, and where one can assess the effects of various choices of parameterizations for the form factors used in order to analyze the numerically generated data (which we will henceforth refer to as pseudo-data). This framework is provided by the scalar form factors of the pions, defined as π 0 ( p 1 )π 0 ( p 2 )| m(uu + dd)(0)|Ω = +F π 0 S (s) [s ≡ ( p 1 + p 2 ) 2 ], π + ( p + )π − ( p − )| m(uu + dd)(0)|Ω Expressions of these form factors, with isospin-breaking contributions due to the difference of masses between charged and neutral pions included, have been recently obtained in [25] up to and including two loops in the low-energy expansion. We will use these expressions in order to generate pseudo-data, which we can then submit to analysis, using various parameterizations for the form factors, inspired by those in use for the analyses of the K +− e4 and K 00 e4 experimental data. The reason for working with the scalar form factors is at least twofold. First, the form factors, with isospin-breaking effects included, are known at two loops in both channels, whereas in the K e4 case, only the form factors in the channel with two charged pions have been studied at the same level of accuracy as far as isospinbreaking corrections are concerned [23] (see Ref. [26] for a systematic study at one loop). Second, the K e4 form factors depend on two more kinematical variables, besides the di-pion invariant mass. The scalar form factors depend only on the latter, and offer therefore a simple kinematical environment, so that the issues we wish to focus on can be addressed without unnecessary additional complications.
In a second step, we address the issue of radiative corrections to the total decay rate of the decay K ∓ → π 0 π 0 e ∓ (−) ν e . Our intent here is not to develop a full one-loop calculation, at the same level of precision as those that exist for the decay channel into two charged pions [12][13][14]. We rather aim at providing a simple estimate for the radiative corrections to the total decay rate, much in the spirit of Refs. [15] and [27] or, on a more general level, of Ref. [28]. This will allow us to assess how much of the discrepancy (1.2) has to be ascribed to other isospin-breaking effects in the form factors, such as discussed above.
The remainder of this study is then organized in the following way. First, we give (Sect. 2) a theoretical discussion of the structure of the scalar form factors of the pions using the explicit expressions obtained in Ref. [25]. We will thus adapt the discussion of Ref. [11] to the case at hand. Working on this analogy will allow us to give an assessment of some additional assumptions regarding the structure of the form factors implicitly made in Ref. [11]. Next, we generate pseudo-data (Sect. 3) using the known twoloop expressions of the form factors, which we then analyze using various phenomenological parameterizations, which do not necessarily comply with the outcome of Sect. 2. The purpose here is to discuss in a quantitative way the possible systematic biases that can be induced by these different choices. The last part of Sect. 3 addresses the determination of the combination a 0 0 − a 2 0 of S-wave scattering lengths. Radiative corrections, aiming at an estimate of the correction factor δ E M in Eq. (1.1), are discussed in Sect. 4. We first compute radiative corrections to the K 00 e4 decay rate in a similar set-up to the one used for the treatment of the data in the K +− e4 channel, in order to obtain a meaningful comparison between the two channels. Then we compute the effects of additional photonic corrections, not included in this treatment. Finally, we end our study with a summary and conclusions. Two appendices contain technical details relevant for the discussions in Sects. 2 and 4, respectively.

Describing the cusp: theory
According to the general analysis of Ref. [20], the occurrence of both π 0 π 0 and π + π − intermediate states at different thresholds leads to the following structure for the scalar form factor of the neutral pion F π 0 S (s) (M π stands for the charged-pion mass, whereas M π 0 is the mass of the neutral pion): where F π 0 0 (s) is a function of s that is smooth as long as no other threshold, corresponding to higher intermediate states, is reached. Here δ(s) represents a phase. It can be chosen arbitrarily, as long as it is also a smooth function of s. The cusp at s = 4M 2 π observed in the differential decay rate corresponding to this simplified (as compared to K e4 ) situation then results from this decomposition, since (2.2) Apart from the dependence with respect to the second kinematical variable s e , the empirical parameterization used for the fit to the K 00 e4 data, Eq. (9.1) in Ref. [11], complies with this general representation provided (the variable S π used in this reference corresponds to the variable s used here): A more theoretically based parameterization, adapted from the simple discussion of the cusp in K ± → π 0 π 0 π ± given in Ref. [19], is considered in Sect. 9.4 of Ref. [11], though not used for the data analysis. As compared to Eq. (2.2), its validity also rests on additional assumptions, which, once transposed to the present situation, read: (1 ) The phase δ(s) can be chosen such as to make the two functions F π 0 0 (s) and F π 0 1 (s) simultaneously real, so that Eq. (2.2) takes the simpler form . (2.4) (2 ) F π 0 1 (s) is related to the scalar form factor F π S (s) of the charged pion, multiplied by a combination of the two S-wave ππ scattering lengths a 0 0 and a 2 0 in the isospin limit, In view of the discussion in Ref. [11], F π S (s) should be identified with the phase-removed scalar form factor of the charged pion. The latter is given by . Our purpose in this section is twofold. First, we will rewrite the two-loop representation of the form factor F π 0 S (s) obtained in Ref. [25] in the form (2.2) that makes the cusp structure explicit. Second, we will assess to which extent the additional features mentioned above and assumed in Ref. [11] are actually reproduced by the structure of the form factors at two loops in the low-energy expansion. In particular, we will establish the precise relation between F π S (s) and F π S (s) in Eq. (2.5) at this order. In what follows, and unless otherwise stated, it will always be understood that s ≥ 4M 2 π 0 . Furthermore, in practice s ≥ 4M 2 π will actually mean 4M 2 π ≤ s ≤ M 2 K , where M K is the mass of the charged kaon, so that we need not worry about thresholds other than those produced by two-pion intermediate states.

The cusp in the one-loop form factors
We start with the study of the cusp using the one-loop expression of the form factor F π 0 S (s), In this expression, a π 0 S denotes a subtraction constant, that we need not specify further for the time being. The loop functionsJ 0 (s) andJ (s) are given bȳ with The functions ϕ 00 0 (s) and ϕ x 0 (s) denote the lowest-order real parts of the S-wave projections of the amplitudes of the processes π 0 π 0 → π 0 π 0 and π 0 π 0 → π + π − , respectively. Their expressions read with [23,25,29] a 00 = 2 3 a 0 0 + 2a 2 and Δ π ≡ M 2 π − M 2 π 0 . In the range of s under consideration, the functionJ 0 (s) is complex, but both its real and imaginary parts are smooth, whereasJ (s) may be rewritten as . (2.12) The two functionsJ [0] (s) andJ [1] (s) are smooth, and read , (2.14) In these expressions, the definition of σ (s) has been extended below s = 4M 2 π by 2 .
(2.15) According to Eqs. (2.1) and (2.2),J (s) exhibits a cusp structure at s = 4M 2 π . One thus obtains, at this order, the decomposition of the form (2.1) for F π 0 S (s), with (2.16) provided one factorizes the global phase Therefore, up to so far unspecified higher-order corrections, the one-loop expression of the form factor can be brought into the form (2.1). Both functions F π 0 0 (s) and F π 0 1 (s)/σ (s) are real and smooth for s ≥ 4M 2 π 0 at this stage. However, the expression for F π 0 1 (s) in (2.16) does not quite comply with Eq. (2.5). Whereas at this stage F π S (s) is equal to the constant F π S (0), which, at this order, can be identified with the phaseremoved form factor, the combination of scattering lengths that occurs in Eq. (2.5) corresponds to (s) is the expression of ϕ x 0 (s) in the isospin limit. Thus, at this order, the expression (2.5) misses both the dependence with respect to s in ϕ x 0 (s), and the isospin-breaking corrections in the scattering lengths.
For later reference, we briefly extend the discussion to the scalar form factor of the charged pion. At one loop, it is given by Besides the subtraction constant a π S , which differs from a π 0 S (they become identical in the isospin limit), this expression involves the lowest-order real part of the S-wave projection of the amplitude for the scattering process π + π − → π + π − , where [23,25,29] (2.20) After having factorized the global phasẽ At this stage, both functions F π 0 (s) and F π 1 (s)/σ (s) are real and smooth for s ≥ 4M 2 π 0 .
2.2 The cusp in the two-loop form factor F π 0 S (s) of the neutral pion Let us now go through the same analysis, but with the twoloop expression of the form factor. The expressions of the pion scalar form factors at two loops and in presence of isospin breaking have been worked out in Ref. [25] using a recursive construction based on general properties like relativistic invariance, unitarity, analyticity, and chiral counting. The scalar form factor of the neutral pion can be written as 3 In this formula, the functions ξ x (s), . . . are polynomials of at most second order in the variable s. Their expressions can be found in Ref. [25], except for Δ 1 ξ x (s), which reads (2.24) It is also useful to be aware of the relations ξ (2;0) In order to achieve the decomposition (2.1), we need to extend the decomposition of the functionJ (s) in Eqs. (2.12) and (2.13) to the other functions, denoted generically byK α n (s), which appear in the expression (2.23). This may be done as follows. First, we may observe that, likē J 0 (s) orJ (s), these functions can also be defined by a dispersive representation of the form [forJ 0 (s) ≡K 0 0 (s), Explicit expressions for the functions k α n (s) are given in A. For the set of functions K 0 n (s) andK ∇ n (s), one has s thr = 4M 2 π 0 . These functions will therefore each develop an imaginary part for s ≥ s thr = 4M π 0 , ImK α n (s) = k α n (s)θ (s − 4M 2 π 0 ), while the real part displays a cusp at s = 4M π 0 , but is smooth for s ≥ 4M π 0 . The situation is different for the remaining functions,K n (s),K x n (s), andK x (s), for which s thr = 4M 2 π , so that, in a generic way, they have the following structure: . (2.27) In general, the function k α n (s)/σ (s), although real, is not smooth for the whole range s ≥ 4M 2 π 0 , but only for s ≥ 4M 2 π . Suppose one can find a functionk α n (s) such that it coincides with k α n (s) for s ≥ 4M 2 π , and such thatk α n (s)/σ (s) is real and smooth for all s ≥ 4M 2 π 0 . Then one can perform the decomposition , (2.28) in terms of two real and smooth functionsK [0] (s) and K [1] (s), given byK Such a decomposition can indeed be achieved for the various functions considered here, as discussed in detail in Appendix A. The decomposition (2.1) of the form factor now follows immediately, with and Both functions e iδ(s) F π 0 0 (s) and e iδ(s) F π 0 1 (s)/σ (s) are smooth for s ≥ 4M 2 π 0 , but complex. It remains to discuss the phase δ(s). If we want to make the function F π 0 0 (s) real, while keeping it smooth, then its choice is unique, withψ 00 0 (s) given by Note thatψ 00 0 (s) differs from the quantity ψ 00 0 (s) defined in Eq. (4.6) of Ref. [25] by the presence of the functionL(s) instead of L(s) in the last term between square brackets, see Eq. (A.2). This makes σ 0 (s)ψ 00 0 (s) a smooth function for s ≥ 4M 2 π 0 , whereas σ (s)L(s), and hence ψ 00 0 (s), displays a cusp at s = 4M 2 π . Making use of Eq. (2.27), the removal of the phase δ(s) indeed leads to a real and smooth expression for the function F π 0 0 (s): As far as F π 0 1 (s) is concerned, we may even proceed in a more direct way by noticing that, up to higher-order corrections, Eq. (2.31) rewrites as with (for the notation, see Appendix A) It is possible to give a more precise interpretation of the combination ϕ x 0 (s) +ψ x 0 (s) that occurs in (2.37). To this end, let us recall from Ref. [25] that the = 0 partial-wave projection f x 0 (s) for the scattering amplitude of the process π 0 π 0 → π + π − is given, at the order of one loop and for s ≥ 4M 2 π 0 , by where ψ x 0 (s) is defined in Eq. (4.15) of Ref. [25] [the contribution Δ 2 ψ x 0 (s), of second order in isospin breaking, is numerically quite small, and is omitted for simplicity]. It differs fromψ x 0 (s) by the replacement ofk 2 (s) by k 2 (s) in Eq. (2.36). Then applying the decomposition (2.28) to f x 0 (s), one finds with f x We may summarize this theoretical study of the cusp in the scalar form factor of the neutral pion with a couple of remarks: -It is, in general, not possible to choose the phase δ(s) in Eq. (2.1) such as to make both F π 0 0 (s) and F π 0 1 (s) real simultaneously. A relative phase remains; see Eqs.
(2.37) and (2.38). At lowest order, this phase is given by the S-wave projection of the inelastic rescattering of a pair of neutral pions through a pair of charged pions, cf. Eq. (2.21).
-The structure of F π 0 1 (s) is more complicated than just the product of the scattering length corresponding to this rescattering amplitude times the phase-removed scalar form factor of the charged pion. At the order we have been working, it involves the decomposition (2.40) of the Swave projection of this amplitude times the part F π 0 (s) of the decomposition (2.1) of F π S (s). This is different from the phase-removed form factor, as already seen at order one loop: (2.41) Note, however, that this difference only concerns the region 4M 2 π 0 ≤ s ≤ 4M 2 π , which contributes very little to the total decay rate as defined by Eqs. (3.1) and (3.2) below.

Description of the two-loop form factor F π S (s) of the charged pion
We now briefly address the scalar form factor of the charged pion. The issue here is not to describe the cusp, which occurs below the physical threshold at s = 4M 2 π , but to provide the expressions that will be used in the sequel. Again, we will rely on the results obtained in Ref. [25], and rewrite the form factor F π S (s) at two loops in a way that is adapted to our purposes. In particular, we will consider the phase-removed form factor, which reads, in the relevant domain s ≥ 4M 2 π , The functionsK α n (s) that appear in this expression have again a dispersive representation of the form displayed in Eq. (2.26). The absorptive parts are in part given in Appendix A. For the remaining ones, we have for the functions corresponding to s thr = 4M 2 π 0 , and, for those whose dispersive integrals start at s thr = 4M 2 π , with the definitions of the functions σ (s), σ 0 (s) and L 0 (s) for s ≥ 4M 2 π 0 given in Eqs. (2.8), (2.11), and (2.15).

Generation and analysis of the pseudo-data
In order to study the effect that particular choices of phenomenological parameterizations of the form factors can have on the output, we will first generate numerical data sets for the scalar form factors of the neutral and charged pions. The pseudo-data in question consist of the (unnormalized) decay distribution defined by The total decay rate is obtained by integrating the distributions (3.1), convoluted with the K 4 phase space, over the whole physical range (we now consider the electron mode only, and set m e = 0): Since we consider the scalar form factor instead of K 4 form factors, the integration with respect to s involves the phase space only. The overall normalization factor has been chosen to be the one of the K e4 decay, up to the factor 1/|F π S (0)| 2 , introduced for convenience.

Form factors and input parameters used for the generation of pseudo-data
The form factors involved in the preceding expressions are considered as known exactly and are constructed as follows.
For F π 0 S;data (s), we will basically use the decomposition of Eq. (2.1), with F π 0 0 (s) given by Eq. (2.34), and F π 0 1 (s) given by Eqs. (2.37) and (2.38). Unfortunately, for some of the functions involved in Eq. (2.34), like ReK ∇ n (s) or ReK x n (s), explicit analytical expressions are not known. For a numerical approach, we could use their dispersive representation, as given by Eq. (2.26). We have, however, found it more convenient to start from expressions upon which we have full analytical control. For that purpose, one may replace the functions ReK ∇ n (s) by the corresponding functions ReK 0 n (s), and likewise ReK x n (s) by ReK n (s). We also drop the contribution proportional to Δ 1 ξ x (s). In the range of s we are interested in, 4M 2 π 0 ≤ s ≤ M 2 K , the difference induced in F π 0 0 (s) by these changes is numerically very small. For the scalar form factor of the neutral pion, the resulting expression then reads with and Finally, in the case of F π S;data (s), we use the expression (2.42) of the phase-removed form factor, replacing the functions K x 0 n (s) byK 0 n (s), and the functionsK Δ n (s) byK n (s), respectively. This then gives (2;±) In the sequel, we will generate pseudo-data using the expressions presented in this subsection, considered to provide exact descriptions of the form factors. In particular, it is understood that higher-order contributions are considered as vanishing. As already mentioned, in order to work within a framework where we deal with fully analytical expressions of the form factors, we have made some approximations as compared to the two-loop expressions discussed in the preceding section. Numerically, these differences are small, but most important is that the approximations we have made preserve the general features of the form factors as described after Eq. (2.40).
For the numerical generation of the pseudo-data, we need to fix the values of the various parameters that occur in the expressions of the form factors. As we want to compare different methods of analysis of the K 4 form factors, we only aim at choosing values that are representative of the expected situation in these decays, with some limited arbitrariness in this choice. In the following, we consider the case We fix the subtraction constants of the form factors by requesting that the charged scalar form factor has the typical values r 2 π = 0.60 fm 2 and c π S = 10 GeV −4 [30], leading to a π S = 2.63 GeV −2 b π S = 2.96 GeV −4 (3.9) Using Ref. [25], one can compute the isospin-breaking shift between a π 0 S and a π S . Assuming that c π S = c π 0 S , we obtain a π 0 S = 2.60 GeV −2 b π 0 S = 3.24 GeV −4 (3.10) leading to r 2 π 0 = 0.59 fm 2 . For the remaining parameters, we use the same input values as in Ref. [25]. We normalize the charged scalar form factor to unity at s = 0 and rescale the neutral one accordingly: As an illustration, we quote the values obtained for the total decay rate of Eq. (3.2) from these inputs, for V us = 0.2255: We also show, on the left panel of Fig. 1, the various contributions to the form factor F π 0 S;data (s) obtained with our input values. For comparison, the right panel shows the equivalent results in the case of the parameterization of the K 00 e4 form factor M(s) discussed in Section 9.4 of Ref. [11]. It corresponds, for M 1 (s), to the parameterization of Eq. (3.13) below with f s = 0 and the remaining parameters f s / f s and f s / f s fixed at the central values given in Table 3 of Ref. [3]. For M 0 (s), we take the expression of F π 0 0 (s) given in Eq. (3.14). The parameters, f s , f s0 / f s , f s0 / f s , and f s0 / f s it involves are determined from a fit to the phase-space distribution given by Eq. (9.1) of Ref. [11] at s e = 0, and with the values of the remaining coefficients a, b, and d taken at their central values as shown in Table 1 of that same reference. The overall normalization N is fixed such that the distribution is equal to unity at s = 4M 2 π , s e = 0. One can observe similar features in both plots on Fig. 1 [the absence of an imaginary part in the right-hand plot has already been discussed at the beginning of Sect. 2], suggesting that our subsequent analysis, based on the scalar form factors of the pion, has also some bearing on the K e4 form factor. Note that due to our choice of phase space in Eqs. (3.1) and (3.2), in practice the region of interest in Fig. 1 corresponds to s ≤ 0.15 GeV.
In summary, the form factors used in order to generate the pseudo-data are defined by the expressions (3.5), (3.6) and (3.7) (with vanishing higher-order corrections), together with the values (central values for all parameters, no error bars) of the parameters specified above. The form factors thus defined will be referred to as the "exact" form factors, considered to represent the "truth" to which we will fit different model parameterizations of the form factors, in order to obtain a quantitative determination of the possible biases different parameterizations can have on the output of the analysis.

Phenomenological parameterizations
In order to mimic the situation in K e4 decays, the pseudodata generated with the exact scalar form factors of the pion will now be analyzed using approximate phenomenological parameterizations. For the analysis itself, we will consider a framework close to the experimental set-up for the K +− e4 [1][2][3] and the K 00 e4 [11] decay channels. From here on, we therefore also use q 2 ≡ s/4M 2 π − 1 in addition to s, the square of the center-of-mass energy of the di-pion system. The region below the cusp corresponds to q 2 ≤ 0, while positive values of q 2 describe the region above the cusp.
For the charged-pion form factor, this means that we consider a parameterization of the form In the case of the neutral-pion form factor, we consider two parameterizations:  [11] for the K 00 e4 form factor, with the parameters chosen as described in the text after Eq. (3.12) -Model 1: (3.14) -Model 2: In order to check the influence of possible higher-order terms in the q 2 expansion, as compared to the parameterization considered in Ref. [11], we have introduced a coefficient f s , which can be set to zero or kept as a free variable in the fit.
The first parameterization with f s = 0 reproduces exactly the one that was considered in Sect. 9.4 of Ref. [11]. The second parameterization incorporates more information gathered from the theoretical discussion in Sect. 2, while remaining sufficiently simple. Although we have chosen not to distinguish them, the parameters f s , f s , f s , and f s appearing in Eqs. (3.14) and (3.15) are not, a priori, identical to those occurring in the expression (3.13). One issue of the analysis we will present is precisely to determine to which extent e.g. f s in Eq. (3.14) should be expected to agree with f s in Eq. (3.13). Our first task is therefore to provide reference values for the various parameters. This is done by performing, in Eqs. The comparison with our various fits will illustrate how different the parameters extracted from the fit and those describing the real Taylor expansion can be, and they will thus give information on the possible bias introduced by the fitting procedure. For convenience, in the following f s0 , f s0 . . . will be called "neutral" parameters, whereas f s , f s . . . are referred to as the "charged" parameters.

Fitting procedures
In order to stay close to the NA48/2 experimental set-up, we will thus assume that we have measurements of |F π 0 S (s)| 2 at the 12 points corresponding to the barycenters of the experimental bins, and we assign a statistical uncertainty derived from the number of events collected in each bin [∼2900 events in the first two bins, and ∼5900 events in all the other ones], without any correlations between the bins. As the parameterizations given in Eqs. (3.14) and (3.15) depend on the S-wave ππ scattering lengths, our χ 2 will also include an uncertainty on these quantities in order to mock up the fact that in the real analysis these scattering lengths are determined from the charged form factor. Here we will use the experimental information on these quantities, namely the latest NA48/2 combination of K → 3π and K 4 results  where we combine statistical and systematic uncertainties in quadrature. One should notice that the central values are Table 1 Results of the different models (1-2) and methods (A-B-C), compared to the reference values, obtained from the Taylor expansion of the exact form factors (3.4), (3.5), and (3.6) with the values of the parameters discussed in the text. The value of the ratio r does not result from the fit, but is obtained once the fit has been performed; see the last paragraph of Sect. 3.3. The column label A1 (A2) refers to the fit method A using model 1 (2), and so on A1 A2 Reference χ 2 /N dof (9.1 ± 4.3)/9 (6.9 ± 3.6)/7 close (but not identical) to the "true" values used to generate our pseudo-data. We consider the following methods to determine the coefficients of the above models.
-Method A: fit of |F π 0 S (s)| 2 for all points to determine all the (neutral, charged) parameters (setting f s = 0 to ensure a reasonable convergence of the fit), assuming the equality of the neutral and charged normalization ( f s = f s0 ). -Method B: fit of |F π 0 S (s)| 2 for all points to determine all the (neutral, charged) parameters, setting f s = 0 and keeping the normalizations f s and f s0 distinct. -Method C: fit of |F π 0 S (s)| 2 to determine the neutral parameters, injecting information on charged parameters by adding to the χ 2 a contribution corresponding to a fit of the charged form factor F π S (s) to the polynomial expression (3.13), effectively identifying the charged parameters in the models for F π 0 S (s) with the parameters occurring in the charged scalar form factor.
In method C, we generate pseudo-data points for the chargedpion scalar form factor with energies corresponding to the barycenters given in Ref. [1], and we use the relative uncertainties for F S (combined in quadrature) quoted for each bin in this reference, without correlations. In agreement with Ref. [1], we add an overall 0.62 % relative uncertainty, completely correlated between all the charged bins. The curvature of the charged form factor F π S (s) being more pronounced than that of the scalar K ± e4 form factor, a q 6 term must be included in the polynomial in order to obtain a good description of the form factor over the whole kinematic range.
We give the resulting χ 2 min (obtained from the best-fit values of each method). Even though each model provides through its fit a value of f s0 , one can also determine the latter by considering the branching ratio. In this case, f s0 is determined by integrating the decay distribution obtained by using as inputs the slope parameters determined from the different methods of fitting, and fixing the normalization by comparison with the total decay rate Γ π 0 defined in Eq. (3.2), and evaluated with the exact form factor F π 0 S;data (s). The corresponding numerical value is given in Eq. (3.12). We denote by r the ratio between the value of f s0 determined this way from the branching ratio, and the true value computed from the exact form factor, i.e. the reference value f s0 = 1.381.

Discussion of the results
In order to obtain an estimate of the uncertainty attached to the coefficients of models 1-2 using methods A and B, we will perform fits of the models on a series of 10000 pseudoexperiments, generated by assuming that the data are random variables with a mean given by our theoretical model for the neutral scalar form factor and a standard deviation given by the relative uncertainty of the corresponding form factor measured in K e4 decays by the NA48 experiment [1,11]. We will then determine the mean and the variance of the resulting distribution for each coefficient of the parameterization considered. The results are gathered in Table 1. The column labeled "Reference" provides a comparison with the coefficients obtained from the Taylor expansions of the form factors given in Eqs. (3.5) and (3.7), as described after Eq. (3.15).
As shown by the comparison between models 1 and 2 for methods A and B, the higher powers of q 2 are only weakly constrained. Model 1 is very rough and provides a very poor description of the charged form factor (modeling it as a simple constant), which explains the very bad χ 2 min for method C1. Only f s0 and f s0 / f s0 can be determined with a good accuracy, but there is no significant bias introduced by the fitting procedure with respect to the reference values. Despite its shortcomings, method A gives good results for the neutral parameters. As expected, compared to method B, method C provides a much better accuracy on the neutral parameters since the charged ones are constrained in this method. As shown by the ratio r , both methods yield accurate values of f s0 (at the few percent level) obtained by integrating over the phase space to consider the branching ratio, even methods that do not attempt at describing the q 2 < 0 region correctly. This can easily be understood: both methods are constrained to describe correctly |F π 0 S (s)| 2 for small q 2 > 0 (as can be seen by their agreement concerning f s0 and f s0 ), but they may differ for q 2 < 0 (which exhibit larger uncertainties). However, this region is very narrow (4M 2 π 0 ≤ s ≤ 4M 2 π ) and its contribution is further sup-pressed by phase space. Therefore, the impact of this region on the estimation of the branching ratio is very small, and the latter is completely dominated by the region s ≥ 4M 2 π where all parameterizations agree (the uncertainties reflecting mainly the uncertainties of the inputs and the lack of data at large q 2 ). From this discussion, one thus expects that using the fit function given by Eq. (9.1) of Ref. [11], and described at the beginning of Sect. 2, will lead to similar results for the ratio r , despite the fact that this model complies with the expected structure of the cusp only if one imposes strong assumptions, and should be considered as a mere phenomenological parameterization to reproduce data smoothly. One indeed obtains r = 1 ± 0.01 and χ 2 /N dof = (8.0 ± 4.1)/9 for method A, and r = 1 ± 0.01 and χ 2 /N dof = (8.1 ± 4.1)/8 for method B, illustrating once more that a smooth parameterization of the curve above the cusp in good agreement with the data is enough to obtain an accurate and unbiased value for the normalization f s0 .
The outcome of this discussion is that the measurement of |F π 0 S | 2 allows for an accurate determination of f s0 (at the percent level), in the current experimental setting. As shown by the ratio r , the value of f s0 obtained from the computation of the branching ratio is equal (within uncertainties) to its true value for all methods and parameterizations considered here. Even though one has to keep in mind that this observation is done using the pion scalar form factors rather than the actual K 4 form factors, it nevertheless suggests that the fit procedure adopted in Ref. [11] does not bias the determination of f s0 , and thus cannot explain the surprisingly higher value of f s extracted by the NA48/2 Collaboration from the K 00 e4 channel, as compared to the value for f s determined from the K +− e4 channel.

Constraining the scattering lengths
The presence of a cusp similar to the one observed in the three-body K + → π + π 0 π 0 decay suggests that it should, in principle, be possible to extract information on the scattering lengths from an accurate measurement of the K 00 e4 differential decay rate. At leading order, the cusp is related to the difference of scattering lengths a 0 0 − a 2 0 . Going to higher orders in F π 0 1 (i.e. including ϕ x 0 ) will also add a [weaker] dependence on a 2 0 . The scattering lengths can be determined only once the relative normalization of form factors involved in F π 0 0 and F π 0 1 is fixed, which requires the determination of the charged parameters in some way. We define two methods for this purpose. Method D is exactly as method A, without including any experimental information on a 0 0 and a 2 0 in the χ 2 [i.e. removing them from the χ 2 as described in Eq. (3.16)], and similarly for method E with respect to method C. We proceed as before, but now also fitting the scattering lengths. We consider model 2, as model 1 yielded poor results in the previous section for method C. In order to discuss the potential impact future experimental improvements could have, we consider also a situation where all statistical errors are reduced by 10 (but the number of bins is unchanged) for the neutral channel, keeping the uncertainties unchanged for the charged channel.
The results gathered in Table 2 show that the current statistical uncertainties yield a relative uncertainty on a 0 0 − a 2 0 of around 80 % for D2 and 40 % for E2. For D2, the charged parameters are only very poorly constrained, but this does not prevent the fit to be reasonable. Reducing the statistical uncertainties by 10 (for the neutral part) yields a significant reduction in the uncertainties, leading to a relative uncertainty on a 0 0 − a 2 0 of 27 % for D2 and 10 % for E2. At this level of accuracy, there is no significant bias in the value of a 0 0 − a 2 0 extracted through these various approaches. As expected, no relevant information can be obtained on a 2 0 , due to the very small sensitivity of the neutral-pion channel to this quantity. To illustrate this point, if instead we fix the value of a 2 0 to its central value in Eq. (3.16), our results concerning the uncertainty on a 0 0 − a 2 0 and the quality of the fit remain unchanged.
From this discussion, one can hope to get some information on a 0 0 −a 2 0 using model 2, should a larger data set become available for K 00 e4 in the future. One has, however, to keep in mind that we have assumed the equality between the charged and neutral normalizations in the polynomials for the neutral scalar form factors in the case of method D, as well as the equality between the charged parameters in the polynomials for the charged and neutral scalar form factors in the case of method E. These assumptions are certainly reasonable considering the current uncertainties involved, but one might need to reassess them in the presence of more accurate data. In this context, it is also interesting to notice that the current result from the DIRAC experiment Ref. [31] is |a 0 0 − a 2 0 | = 0.253 ± 0.011, i.e. a 4.3 % uncertainty, so that a substantial increase of the statistical sample of K 00 e4 decays is needed in order to reach a comparable accuracy.

Radiative corrections to the K 00 edecay rate
In this section we now discuss radiative corrections, which were addressed differently in the analyses of the K +− e4 and K 00 e4 channels so far. In the latter case, no radiative corrections were applied to the decay rate measured in Ref. [11].
This accounts for the unspecified factor δ E M in Eq. (1.1). It is thus natural to ask how much of the observed 6.5 % discrepancy [see Eq. (1.2)] in the normalization of the form factor measured in the two channels is due to this correction factor. Our aim here is not to provide a complete discussion of radiative corrections in the K e4 decay channels at a level of sophistication that would match the treatment of isospin breaking due to the difference between masses of the charged and neutral pions. We rather want to work out these corrections in a somewhat simpler framework, trying to reproduce a treatment of radiative corrections in the neutral channel similar to the one that was applied in the charged channel, in order to make the comparison as meaningful as possible.

Treatment of radiative corrections in
Let us recall how radiative corrections are treated in the charged channel [2,3]. First, virtual photon exchange between all possible pairs of charged external lines are considered, and the corresponding Sommerfeld-Gamow-Sakharov factors are applied. The corrections induced by emission of real photons are treated with PHOTOS [32][33][34][35]. The latter also implements the wave-function renormalization on the external charged legs. The couplings of photons to mesons are treated as point-like interactions, given by scalar QED. The result is then free from infrared singularities. Furthermore, one neglects the contributions that vanish when the electron mass goes to zero, which is a sensible limit to consider for the K e4 decay channels. Apart from the Sommerfeld-Gamow-Sakharov factors, some contributions that would arise within a more systematic approach, provided by the effective low-energy theory of QCD and QED for light quarks and leptons [37][38][39], as applied in Refs. [12][13][14] to the channel with two charged pions, are not considered. These include, for instance, all structure-dependent corrections, where the photon is emitted from the tree-level K e4 vertices or from internal charged lines. The outcome of such a truncated calculation is affected by an ultraviolet divergence, which is removed by renormalizing the coupling |V us | 2 G 2 F [note that in Eq. (12) of Ref. [35] the factor (α/π ) has been inadvertently omitted], This same correction factor also appears in Ref. [27], with the ultraviolet cut-off Λ taken equal to M W . From this last reference, we also see that the factor 9/4 decomposes as 9/4 = 3 × (1/2) − 1/4 + 1/2 + 1/2, where the first contribution comes from the wave-function renormalization of the three charged mesons, the second from the (charged) lepton wave-function renormalization, and the last two ones from the virtual photon loops between the charged kaon and the charged lepton on the one hand, and between the two charged pions on the other hand (the remaining divergent contributions of this type, i.e. a photon line connecting the external kaon to the charged-pion lines, or the charged lepton with each of the two pions, cancel pairwise).
In the case of the K 00 e4 channel, we therefore expect that this factor becomes 1 × (1/2) − 1/4 + 1/2 = 3/4. Since it differs from the previous one, it cannot be absorbed by the renormalization of the same prefactor |V us | 2 G 2 F as before. It seems more natural instead to absorb these ultraviolet divergences into the normalizations of the form factors This is also more in line with the effective theory approach mentioned above, where the form factors are also corrected by (different) contributions from the low-energy constants K i [38] or X i [39], which are renormalized by the ultraviolet divergences coming from the photon loops. Using instead Eq. (4.1) in both cases would leave a remaining cut-off dependent contribution to the K 00 e4 amplitude. For a typical value of Λ = 1 GeV, this would modify Eq. (5.1) at the per mille level.

Radiative corrections à la PHOTOS for the K 00
e4 decay rate In the following, we will try to estimate the potential impact of PHOTOS on K 00 e4 rather than pursuing an effective field theory approach. If we want to reproduce the analog of the PHOTOS treatment [35] of radiative corrections for the K 00 e4 decay rate, we need to consider the wave-function renormalization of the charged lepton and of the kaon in (scalar) QED, and the vertex correction corresponding to diagram (a) in Fig. 2. Using a Pauli-Villars regularization, and taking the photon propagator in the Feynman gauge, we reproduce the expressions of Eq. (6) in Ref. [27] for the former. In order to evaluate and discuss the contribution from diagram (a) in Fig. 2, we choose to describe the tree-level K 00 4 vertex as with constant form factors [Bose symmetry forbids a contribution of the form G 00 ( p 1 − p 2 ) μ with G 00 constant], so that the lowest-order amplitude reads (4.4) In the limit m e → 0, we obtain A(K 00 e4 ; 2(a)) = e 2 A 0 (K 00 e4 ) 1 16π 2 ln The various loop functions occurring in this expression are defined in Appendix B. Adding to it the wave-function renormalizations on the charged external lines gives the following result, in the framework adopted here, for the radiatively corrected amplitude (m γ denotes a small photon mass, introduced as an infrared regulator, to be sent to zero once an infrared-safe observable has been constructed): We make a few comments about this result: -Although the result (4.5) holds in the Feynman gauge ξ = 1, we have also computed the wave-function renormalizations and A(K 00 e4 ; 2(a)) in an arbitrary linear and covariant ξ -gauge, and we have checked that the final result (4.6) does actually not depend on the gauge-fixing parameter ξ .
-In order to reproduce the analog of the PHOTOS treatment [27,35] of radiative corrections for the K 00 e4 decay rate, one needs to add the emission of soft photons from the charged external lines, diagrams (a) and (c) of Fig.  3 (in the K +− e4 case, there are two more diagrams where the photon is emitted from the charged-pion lines) so that the result is free of infrared singularities at order α. These corrections will be discussed later on. At this stage, we simply note that the infrared-divergence of Eq. (4.6) is equal to with the function τ ( p 1 , p 2 ) defined in Eq. (B.21). -The factor 3 8 α π ln Λ 2 M 2 π discussed in Eq. (4.2) is indeed to be found in Eq. (4.6), provided that one sets R 00 to zero. The only remaining contribution comes from F 00 , which is proportional to f s0 at this level. This indeed corresponds to the situation considered in Ref. [27]. In the absence of radiative corrections, the form factor R 00 (or R +− in the charged channel) does not contribute to the decay distribution for m e = 0. In this case, one may as well take R 00 = 0 from the beginning. But once radiative corrections are switched on, taking m e = 0 or R 00 = 0 are no longer equivalent options. As shown by the second contribution in Eq. (4.6), there is a correction to F 00 that is induced by R 00 , and this contribution is not considered in Ref. [27], and is hence also missing in Ref. [35].
-At lowest order, and in the isospin limit, one has [36] R 00 = F 00 (4.8) Actually, as shown in the second expression, the vertex in the diagram (a) in Fig. 2 only accounts for the contribution R 00 = (2/3) · F 00 . The second factor comes from the diagram ( f ) in Fig. 2. 4.3 Additional non-factorizable radiative corrections to the K 00 e4 decay rate After these preliminary remarks concerning the PHOTOStype treatment of radiative corrections in the K +− e4 and K 00 e4 channels, let us now address radiative corrections in the channel with two neutral pions in a somewhat more systematic manner. This will allow us to estimate the size of the radiative corrections that are not included in the experimental analysis, as described in the previous section. We keep on considering the limit where m e vanishes, so that in the absence of radiative corrections the amplitude reads simply For our purpose, it is convenient to distinguish between two types of radiative corrections, which we call factorizable and non-factorizable. Factorizable radiative corrections are defined by the contributions where both ends of the virtual photon line connect to a charged mesonic line or to the vertex with the leptonic current, or when both ends connect to the electron line. These factorizable contributions will not modify the structure of the matrix element, but will change the form factors F 00 and R 00 . We find it convenient to express them as 4 where we have factored out the wave-function renormalization factors computed in QED for Z e , and in scalar QED for Z K : In contrast to the preceding subsection, we use now dimensional regularization, with the minimal subtraction of the combination (4.12) 4 In the K +− e4 case, there are additional structures in A fact (K +− e4 ), due to the possibility, already at tree level, for a virtual photon to emit a pair of charged pions; see [12,14]. Notice in this respect that the contributions in Fig. 6 and in Fig. 8 of Ref. [14] vanish in our case.
It is understood that F 00 in Eq. (4.10) now includes all the remaining factorizable photonic corrections, together with the contributions from the low-energy constants L i [37] and K i [38]. These will take care of the UV divergences due to the meson loops and to the photon loops, respectively, so that the product F 00 √ Z K is actually UV finite. It, however, inherits the infrared divergence contained in √ Z K . Let us next consider the non-factorizable contributions. As far as the corrections to F 00 are concerned, they are represented by the diagrams shown in Fig. 2. One finds that the contributions coming from the diagrams (b), (c), and (d) are proportional to the lepton mass, and thus vanish in the limit m e → 0. There are therefore only three diagrams to compute in this approximation. Consistently dropping terms that vanish as m e → 0, one finds that and Apart from the change of regularization, the expression for A(K 00 e4 ; 2(a)) in Eq. (4.13) reproduces the one of Eq. (4.5) obtained previously, provided one takes R 00 = (2/3) · F 00 , as discussed at the end of Sect. 4.2, and makes use of the identities where the ellipses denote terms that vanish in the limit m e → 0.
Adding up the contributions discussed so far, one obtains an expression for the radiative corrections at order O(α) that still contains both infrared and ultraviolet divergences. The latter will be taken care of by the contributions from the counterterms X i introduced in [39]. Their contribution reads The low-energy constant X 1 is not renormalized, whereas X 6 = X r 6 (μ) − 5λ. Collecting the divergent pieces from the various contributions leads to (we recall that at this stage F 00 √ Z K has already been made UV-finite through the contributions of the low-energy constants L i and K i ) (4.18) which vanishes, as it should.
As to the infrared divergences, collecting the IR-divergent pieces contained in the contributions computed so far, one obtains , (4.19) with the function τ ( p 1 , p 2 ) defined in Eq. (B.21). Besides the wave-function renormalizations, such divergences only arise from the contribution of C (k − p ) 2 ; M K , m e in A(K 00 e4 ; 2(a)). Notice that this infrared divergence coincides with the one of the result (4.6), given in Eq. (4.7). The construction of an infrared-safe observable at order O(α) requires also to consider the process with the emission of one soft photon. The corresponding differential decay rate is given by (4.20) The amplitude A(K 00 e4γ ) for the radiative decay K ± → π 0 π 0 e ± (−) ν e γ can be expanded in powers of the photon energy, The Low approximation consists in keeping A −1 (K 00 e4γ ) alone. This is enough in order to study the emission of only soft photons and to discuss the issue of infrared divergences. Explicitly, one has [q μ is the momentum of the emitted (real) photon, ε * (q) the corresponding polarization vector] One may then perform the integration over the undetected soft photon. In the soft-photon approximation, the photon momentum in the delta-function of the phase-space integration is neglected, and one takes Expressions for the corresponding integrals can be found in [12,14]. The integration is limited to photon energies |q| below the experimental detection threshold ΔE in the kaon rest frame. As far as the infrared divergences are concerned, one has (4.25) Therefore, the contributions proportional to ln m γ cancel in the sum dΓ (K 00 e4 ) + dΓ (K 00 e4γ ). For later convenience, we rewrite Eq. (4.24) in a way that explicitly displays the IRsingular part: We then add the contribution to the amplitudes involving virtual photons, such as to make them infrared finite. To this end, we define the function We can now add the virtual and real contributions described up to now which should be involved in a PHOTOS-like treatment of this decay. We include them as a correction of the form (1 + δ EM ) to the determination of the form factor from the measurement of the branching ratio. To this end, we compute the total decay rate including the soft-photon emission where Γ (K 00 e4 ) includes corrections at first order in the finestructure constant α, and write it in terms of the decay rate Γ 0 (K 00 e4 ) without radiative corrections in the form Let us first discuss the corrections computed in subsection B above. In order to obtain a result that is as close as possible to the treatment of radiative corrections in the K +− e4 channel, we absorb the UV-divergent factor of Eq. (4.6) in F 00 and take R 00 equal to zero. Then Γ (K 00 e4 ) is computed by performing the phase-space integration of We take ΔE = 11.7 MeV, the value corresponding to the NA48/2 experiment [11], for the real-photon detection threshold in the kaon rest frame. This gives then This value has the expected size. Moreover, it goes into the right direction, in the sense that it reduces the discrepancy in Eq. (1.2) from 6.5 to 4.6 %. As a test of the stability of the result (4.33) we may also evaluate the non-factorizable radiative corrections corresponding to all the diagrams in Fig. 2. This amounts to taking the expressions in Eqs. (4.13), (4.14), and (4.15) for the evaluation of Γ (K 00 e4 ) (let us stress again that these equations have been obtained in a regularization scheme differing from the one discussed in Sect. 4.2). We absorb the ultraviolet divergences, as well as the contribution (4.17) into F 00 , in order to build a quantity both UV and IR finite. Constant terms have been discarded as they could be included in the contribution of the counterterms X 1 and X 6 . The resulting expression for ΔA(K 00 e4 ) then reads This value is quite close to the one obtained in Eq. (4.33), so that in the present case the treatment of radiative corrections à la PHOTOS seems to yield stable results even after the inclusion of non-factorizable contributions.

Summary and conclusion
The present study is devoted to isospin-breaking effects in the semileptonic decay of the charged kaon into two neutral pions, K ± → π 0 π 0 e ± (−) ν e . Because of the smallness of the electron mass and of the limited experimental precision, this decay can be described in terms of a single form factor. This form factor also occurs in the description of the decay into two charged pions, K ± → π + π − e ± (−) ν e , and up to isospinbreaking contributions, the two determinations should agree. The present study focuses mainly on two aspects related to this issue: (1) to ascertain quantitatively to which extent the phenomenological parameterizations used in order to analyze the data could impinge on the resulting value of the normalization or on the shape of the form factor measured in the decay K ± → π 0 π 0 e ± (−) ν e , and (2) to obtain a quantitative estimate of the radiative corrections to the total decay rate, which again might affect the normalization of the form factor.
Concerning the first issue, we have considered the form factors of the pion as a case study. As a first step, we have discussed the structure of the form factors, and their properties linked to the presence of a cusp, using exact expressions of the form factors valid up to two loops in the low-energy expansion. We have clearly established that the phenomenological parameterizations used in order to analyze the data did not agree with the general properties that can be inferred form these exact expressions. In a second step, we have generated pseudo-data from these form factors, which we have then analyzed using several phenomenological parameterizations. The outcome of this study is that the determination of the normalization of the form factor is actually not sensitive to the parameterizations used. As a side product, we see that the higher orders in the Taylor expansion of form factors are not accurately determined by a direct fit to simplified (polynomial) formulas, as shown in Table 1. Although our study was carried out for the scalar form factor of the neutral pion, we expect that the conclusion also holds for the K 00 e4 form factor. This expectation rests on the following observations: (1) the two form factors have a similar shape, in particular as far as the cusp is concerned, the curvature of the scalar form factor being somewhat more pronounced, see Fig. 1; however, as just stated, the coefficients of the higher powers in q 2 are only weakly constrained by the fits. (2) In both form factors the cusp is described by a constant times the function σ (s) at first approximation. Refining the description of the cusp, as in model 2 for instance, does not impinge on the determination of the normalization of the form factor. This is related to the fact that, in both cases, the region between threshold and cusp is very narrow.
We have also considered the possibility to constrain the ππ S-wave scattering lengths from the measurement of the decay distribution. We have found that, unfortunately, with the sample of events presently available, the statistical uncertainties remain large. A statistical sample comparable to the one available in the K +− e4 channel would be required in order to reach a precision close to that obtained by the Dirac experiment.
The second issue addressed in this paper consists in radiative corrections. We have determined the correction factor δ E M to the total decay rate in Eq. (1.2). In order to make a meaningful comparison with the value of the normalization of the form factor extracted from the K +− e4 channel, we have used a simplified framework, including only those corrections that were also included in the latter case Note that the error bar in this equation is purely from experimental origin, and does not include the systematic uncertainties from the methods used for the evaluation of radiative corrections in both channels. Such additional uncertainties can stem, for instance, from the regularization dependence of the PHOTOS(-like) treatment of radiative corrections, and from neglecting the dependence in the cut-off Λ discussed in Sect. 4.1.
We have also considered additional photonic corrections estimated within a different regularization scheme, and we have found that they do not modify the previous estimate in a significant way. A few comments are in order: -The analysis of radiative corrections we have performed provides an adequate estimate of the global factor δ E M that modifies the total decay rate. It need not be suitable for an analysis of radiative corrections to the phase-space distribution itself. -Other isospin-breaking corrections, among them factorizable exchanges of virtual photons, but also effects due to m u − m d or to the mass differences between charged and neutral pions and/or kaons, are not covered by our analysis. They could affect the normalization of the form factors measured in the two channels in different ways.
A mode elaborate study is needed in order to reach a quantitatively meaningful interpretation of the result in Eq. In view of the value given in Eq. (5.1), the corrections from higher orders to this ratio should therefore be small. -Conversely, using the relation (1.3) in regard to the result (5.1), and discarding yet to be computed corrections to the former, we obtain R = 32 +9 −6 .
The discussion of the radiative corrections presented here is clearly only a first step. In view of the statistical accuracy of the data, a full model-independent calculation of these corrections in the neutral as well as in the charged channels is certainly mandatory before a definite conclusion can be reached concerning the observed difference in the normalization of the form factors between the neutral and the charged channels. This task is clearly beyond the scope of the present note, and is left for future work.