Constraints on the [Formula: see text] form factor from analyticity and unitarity.

Motivated by the discrepancies noted recently between the theoretical calculations of the electromagnetic [Formula: see text] form factor and certain experimental data, we investigate this form factor using analyticity and unitarity in a framework known as the method of unitarity bounds. We use a QCD correlator computed on the spacelike axis by operator product expansion and perturbative QCD as input, and exploit unitarity and the positivity of its spectral function, including the two-pion contribution that can be reliably calculated using high-precision data on the pion form factor. From this information, we derive upper and lower bounds on the modulus of the [Formula: see text] form factor in the elastic region. The results provide a significant check on those obtained with standard dispersion relations, confirming the existence of a disagreement with experimental data in the region around [Formula: see text].

Introduction

Recent years have seen strongly increased interest in the transition form factors of light mesons [1]. One of the main reasons is the fact that pseudoscalar () pole terms constitute some of the most important contributions to hadronic light-by-light scattering, which is soon to become the biggest stumbling block in a more accurate theoretical determination of the standard model prediction for the muon’s anomalous magnetic moment; see [2] and references therein. The strength of these pole terms is determined by the singly- and doubly-virtual form factors for . However, for the virtuality of one of the photons fixed to the mass of one of the light isoscalar vector resonances or , these form factors are intimately linked to vector-meson transition form factors that can be measured in decays such as or . In all likelihood, these vector-meson conversion decays present one of the few opportunities to measure the doubly-virtual or form factors with very good precision. Furthermore, they present an essential ingredient to a more advanced theoretical understanding of hadronic light-by-light scattering [3, 4].

Recent dispersive treatments [5, 6] of the electromagnetic form factor are in disagreement with experimental data in the region around [7-9], which show strong deviations from even approximate vector-meson-dominance behavior [10]. The main ingredient in the dispersion relation is unitarity, which allows one to express the discontinuity of the form factor in terms of the P partial wave of the process [6, 11] and the pion electromagnetic form factor, quantities determined with precision. Strictly speaking, this relation is valid only in the elastic region, . Due to the strong phase-space (and chiral) suppression of multiparticle intermediate states, the elastic regime approximately extends up to threshold in the P wave, , above which the inelasticity in pion–pion scattering is assumed to be dominated by the intermediate state [11]. However, due to the lack of information on the discontinuity in the inelastic region, elastic unitarity is assumed to be valid also at higher energies in the evaluation of the dispersion integral. This assumption may affect the precision of the theoretical treatment. Having in view the disagreement with some experimental data on below , it is of interest to investigate the form factor in a more model-independent framework, which avoids this assumption.

In the present paper we exploit alternative information on the form factor above . We use a method proposed originally by Okubo [12, 13] (before the advent of QCD), which leads to bounds on form factors by exploiting the positivity of the spectral function of a suitable current–current correlator. This technique, known as method of unitarity bounds, has been resuscitated in the QCD era and was applied to a variety of form factors of heavy and light mesons [14-21] (for a review and more references see [22, 23]). In the present study, we use a dispersion relation for the polarization function of two isovector vector currents, calculated by operator product expansion (OPE) in the Euclidean region, and exploit unitarity for the spectral function. Including and states in the unitarity sum, we derive an upper bound for an integral on the modulus squared of the form factor along the cut from to infinity. From this condition and the known discontinuity in the unphysical region , we derive bounds on the form factor in the latter region and compare them with the experimental data and the results of the standard dispersion relations.

In this work, we have adapted the original method of unitarity bounds, making it suitable for the information available on the form factor, namely its discontinuity across the cut below the inelastic threshold. To solve the corresponding optimization problem we use analytic techniques already applied to scattering amplitudes [24, 25]. We also had to take into account the fact that, unlike most form factors studied up to now, the form factor is not a real analytic function.1 We have therefore made the proper generalization of the formalism to the case of analytic functions that are not of real type. To our knowledge such a formulation appears here for the first time and represents an important generalization of methods existing in the literature. The above modifications may have applications in the study of other form factors as well.

In Sect. 2 we briefly summarize the standard dispersive treatment of the form factor. In Sect. 3 we review the formalism of unitarity bounds, which leads to the integral constraint on the modulus squared of the form factor along the cut. In Sect. 4 we solve the optimization problem and derive upper and lower bounds on for below . Section 5 contains the numerical results of our work and Sect. 6 our conclusions.

Standard dispersive treatment of

We use the definition from [5], where the form factor is defined from the matrix elementwhere is the isovector part of the electromagnetic current, denotes the polarization, and we defined and . The form factor has dimension . The definition adopted in [26] contains a factor in the right-hand side of (1) and the corresponding form factor is dimensionless. In the present paper we will work with the dimensionful version of the form factor.

Unitarity implies that has a cut along the real axis for . Using the conventions of [5], the discontinuity of across the cut in the elastic approximation is given bywhere , is the pion electromagnetic form factor, and the P partial-wave amplitude of the scattering processIn [26] the partial wave was calculated by the method, with the left-hand cut approximated by -exchange. In this model, the phase of exactly compensates the phase of in the discontinuity (2), and as a consequence the form factor is a real analytic function. However, as discussed in [5], once rescattering effects are taken into account, the phase of the partial wave no longer coincides with the P-wave phase shift. More precisely, in the projection onto the P partial wave, the kinematical variables reach regions where the decayis allowed and rescattering between the final pions including three-pion cuts is possible. As a consequence, the discontinuity (2) is not purely imaginary and the form factor is not a real analytic function.

In [5] the pion vector form factor has been reconstructed from an Omnès representation [27] using pion–pion phase shifts [28, 29] as input. The partial wave was calculated in [11] using Khuri–Treiman techniques [30].2 Besides not fulfilling a straightforward variant of Watson’s final-state phase relation [31], it also has the peculiarity of showing a singular behavior at the pseudothreshold , which can be understood perturbatively from the analytic structure of certain non-trivial two-loop Feynman diagrams (see e.g. [32]). Note that this does not imply any singular behavior of near : the form factor remains regular on the upper rim of the cut.

The expression (2) is valid in the region , since above the threshold other intermediate states contribute in the unitarity sum. By neglecting these contributions, the form factor was obtained from a once-subtracted dispersion relation [5, 26]where is known experimentally from the decay rate. The updated value is [33]The recent analysis performed in [5], based on the dispersion relation (5), leads to results that are inconsistent with some experimental data around [7-9]. As discussed above, one questionable point of the theoretical analysis is the extension of the elastic unitarity relation (2) above the (effective) threshold of multiparticle production. In the next section we shall show how some information on the modulus of the form factor for can be derived from independent sources.

Consequences of perturbative QCD, analyticity, and unitarity

Using standard techniques [12-23], we start with the QCD vacuum polarization tensorwhere is the isovector part of the electromagnetic current. The first derivative of the QCD vacuum polarization amplitude satisfies the dispersion relationwith the spectral function given by the unitarity relation Keeping the and intermediate states explicitly (which should dominate the isovector spectral function at low energies), carrying out the two-body phase-space integrals and using the positivity of the spectral function, we obtain the inequalitywhereIt is convenient to write (10) as an integral constraint on the modulus of the form factor aswhereThis quantity can be evaluated for spacelike values using OPE and perturbative QCD for the correlator , and the rich information available on the modulus of the pion form factor. The value of should be taken large enough such as to ensure the validity of the OPE, and in the same time lead to sufficiently strong bounds. As discussed in [20], a reasonable choice satisfying these requirements is . We used perturbative QCD to four loops (see [34] and references therein):where is the strong coupling at . Using as input the value , which covers most of the recent determinations from hadronic decays (see for instance [35-39]), and the coupling’s running we obtain . This yields for the central value with an error of about 1.3 %. We also checked that a higher-order term added in (14), taken as according to [35], changes by about 1.2 %. The power corrections in the OPE, obtained from [40], bring a small contribution. For instance, using for the gluon condensate the standard value given in [41, 42], we obtain for its contribution the value . This leads to , where the uncertainty includes quadratically the effects of the uncertainty and the truncation error.

The integral involving the pion electromagnetic form factor can be calculated using in the low-energy region BaBar experimental data [43] and the bounds on derived in [44], along with data obtained by BaBar up to [43] and a smooth transition to the decrease predicted by QCD (for details see [44, 45]). This gives for the integral appearing in (13) the value , which leads toFrom the inequality (12) and the discontinuity (2) adopted in the elastic region of validity , we shall obtain bounds on the form factor at points below . The mathematical technique will be explained in the next section.

Bounds on below

In this section, we sequentially formulate the extremal problem, present its solution, and also discuss the special case of the optimal solution for the case of a real analytic form factor.

Formulation of an extremal problem

In order to cast the problem into a canonical form, the first step is to map the plane cut along onto the unit disk in the plane, by using a suitable conformal mapping. In particular we shall use the mappingsuch that the . In the -plane the elastic region becomes the segment of the real axis, where , and the upper (lower) edges of the cut become the upper (lower) semicircles.

Further, we shall construct an outer function , i.e. a function analytic and without zeros in , its modulus on being equal to , where is the inverse of (16). The construction of the outer functions is explained in [46] (see also the review [22]). Using the expression (11) of , we obtain for the exact analytic expressionThen the inequality (12), written in terms of the new function defined asbecomes3 where .

Since is real analytic in , is real for , and from the definition (18) it follows that we can writewhere the discontinuity of the form factor is obtained from (2). The function can be expressed in terms of its discontinuity aswhere the function is analytic in , as its discontinuity vanishes:Since we consider in general form factors that are not real analytic, the function is analytic, but its values on the real axis may be complex.

We now express the available information on the form factor as a number of constraints on the function . By inserting (21) in (19) we obtain the conditionwhile using (6) we see that has the valueThe problem is to find the maximal allowed range of at an arbitrary given point in the interval , for functions analytic in and subject both to the boundary condition (23) and the additional constraint (24).

Let us denotewhere is an unknown parameter. A simple reasoning (see for instance [25]) shows that the allowed range of is described by the inequalitywhere is the solution of the minimization problemon the class of functions analytic in , which satisfy the constraint (24) and the additional condition (25) for a given .

Solution of the extremal problem

We solve the constrained minimum norm problem (27) by the technique of Lagrange multipliers. We use the fact that the norm squared of a complex function given on the boundary of the unit disk is expressed asin terms of its Fourier coefficientsTherefore, we write the Lagrangian of the constrained minimization problem (27) aswhere are the negative-frequency coefficients of the function, which can be written by applying Cauchy’s theorem asand are the positive-frequency coefficients, defined through the Taylor expansion of the analytic function :We must find the minimum of the Lagrangian (30) with respect to the complex parameters , which are free, except which is known from (24). The minimum conditionshave the solutionswhich, introduced in the constraint (25), give the Lagrange multiplierand the optimal coefficientsBy inserting these coefficients in (30) and performing the summation of the first term, we find the minimum normThen the inequality (26) can be written in terms of the unknown quantity aswhereBy using the definition (25) of and the relations (18) and (21), we write the inequality (38) aswhich leads to upper and lower bounds on :We have taken into account the fact that calculated from (17) is positive for real . Using the value of from (24) and from (39), with defined in (20), we have all the ingredients to evaluate (41) numerically at an arbitrary point .

Optimal solution for real analytic functions

It is of interest to consider in particular the case when the form factor is a real analytic function, i.e. it satisfies . A similar, more general, optimization problem for real analytic functions was solved previously in [24, 25] with a slightly different method. Below we shall use the method of Lagrange multipliers applied above, adapting it to real analytic functions. In this case the discontinuity of the form factor across the cut is equal to , and from (20) we obtain the discontinuity of for aswhereMoreover, from (21) it follows that the parameter defined in (25) is real and given bywhere denotes the Cauchy principal value. Therefore, from (38) we obtain a quadratic inequality with real coefficients for the parameter where is now real and can be written asFrom (45) we obtain exact upper and lower bounds on the real part of the form factor:We have checked that these bounds follow as particular cases from the more general expressions given in [24, 25]. By combining the bounds (47) with the known value of the imaginary part, we derive bounds on the modulus of the form factor. These bounds are optimal, unlike the bounds given in (41), where the treatment of the modulus in the last step of the derivation amounts to a loss of optimality.

We end this section with a remark that might be useful for improving the bounds. From general arguments [22] and the expressions given above, it follows that the bounds depend monotonically on the value of in the -norm constraint (12): smaller values of lead to narrower allowed intervals for at . Therefore, the bounds can be made tighter in principle by taking into account more intermediate states, besides the pairs, in the unitarity relation (9) for the QCD correlator. Many of these are (re)measured exclusively with excellent precision in the ongoing quest to improve on the hadronic vacuum polarization contribution to the muon’s anomalous magnetic moment (see e.g. [47]). The positive contribution of these states can be subtracted from the QCD value of the correlator as in (13), reducing the value of . To establish the practical effect on the bounds requires further investigations.

Results

In the calculation of the bounds we have employed the discontinuity of the form factor in the range from two different analyses: the recent dispersive treatment reported in [5] and the older work [26]. For the quantity we used the estimate given in (15). It turns out that the value of the form factor at plays a significant role in producing stringent constraints. Although in principle can be complex, we assumed that it only has a small phase which can be neglected [5]. The upper and lower bounds given below were obtained using as input the central value .

We checked that the bounds are quite stable with respect to the variation of the input: by varying inside the error quoted in (6), the upper bounds in the region of interest are changed by at most 2.5 %. Also, the uncertainty of the quantity quoted in (15) affects the bounds by at most 2 %. As the experimental errors are currently rather in the 10–20 % range, we refrain from displaying these small variations in the bounds graphically and only discuss the central results.

Our results are presented in Fig. 1, which shows upper and lower bounds on the modulus squared (normalized to its value at ) in the part of the elastic region accessible experimentally in . For the input from [5], when the form factor is not a real analytic function, the bounds on were calculated using (41). For the input from [26], where rescattering effects are neglected and the form factor is real analytic, we used the optimal bounds (47) on the real part, and combined them with the knowledge of the imaginary part to obtain bounds on the modulus. To assess the loss of optimality inherent in (41), we also show the upper bound calculated with this expression for the input from [26]. One can see that the bounds calculated using (41) are very close to the optimal bounds calculated with (47) for the elastic energy range of interest. For comparison, we also show the result of the dispersive calculation performed in [5], as well as several experimental data from [7-9].

Fig. 1

Upper and lower bounds compared with experimental data on . Solid red line Bounds calculated using the expressions (41) with input from [5]. Dashed blue Optimal bounds calculated using (47) with input from [26]. The dotted blue line is the upper bound calculated with the same input [26], but using the nonoptimal expression (41). The data are from Lepton-G [7], NA60 (2009) [8], and NA60 (2011) [9]. The yellow band is the result of the dispersive calculation performed in [5]

Figure 1 shows that, although the allowed ranges for the ratio calculated with (41) and (47) are rather large, the upper bounds exclude some of the data points from [7-9] in the region above . Specifically, the exclusion is at the level for the last point (at ) from [9], and at about the level for the other points lying above the upper bounds. On the other hand, the dispersive calculation performed in [5] is situated in the allowed range for the modulus derived here. Note also that in [48, 49], has been calculated based on a low-energy effective theory including explicit vector-meson degrees of freedom. We have checked that, although the representation [48, 49] rises more quickly than the dispersive ones [5, 26] and follows the data more closely than most others, it still lies comfortably inside the bounds.

Discussion and conclusions

The present study was motivated by the discrepancies noticed recently between the theoretical calculation of the form factor by a dispersion relation and some of the data reported in [7-9]. Our aim was to avoid the assumptions made on the discontinuity of the form factor above the threshold , where the elastic unitarity (2) is no longer valid. To this end we have resorted to the formalism of unitarity bounds. The central point of the formalism is the derivation of an integral condition on the modulus squared of the form factor from to infinity, which can be calculated using OPE and perturbative QCD in the Euclidean region for a suitable correlator, together with unitarity and positivity of the spectral function. In the present case it was convenient to consider the tensor (7) of two isovector currents.4 From the integral condition, by using techniques of analytic interpolation theory [46], one can derive bounds on the form factor or its derivatives at points inside the holomorphy domain.

In this work we have considered a modified version of the standard formalism, suitable for including the information available on the form factor, namely the discontinuity (2) known in the elastic region. For real analytic functions, the bounds (47) are consistent with the more general results given in [24, 25]. We have also derived bounds on the modulus of the form factor in the case when, due to rescattering effects, it is not an analytic function of real type. This is a generalization of the formalism of unitarity bounds, considered for the first time in this paper. The framework is not specific to the form factor and could easily be adapted to the analysis of other form factors. The derivation uses the maximization of the modulus of a difference of complex numbers in the last step, which implies that strictly speaking the bounds (41) are not optimal. In practice, for the real analytic form factor with the discontinuity from [26], the predictions of (41) are almost indistinguishable from the optimal bounds (47) in the energy region of interest.

The numerical results show that several experimental data around are situated above the upper bounds derived in the present paper, using two different evaluations of the discontinuity (2) in the elastic region. Having in view the model-independent treatment of the region above the inelastic threshold adopted in our analysis, the disagreement signals possible problems with the experimental data. Adding to the seeming inconsistency between data on from and data on from [50-52] (see e.g. the compilation in [1] or the continuation of the calculation [5] above shown in [53]), our results therefore strongly support the desirability of renewed experimental efforts to measure the conversion decay in a more exclusive setting [54, 55].