Dispersive analysis of the pion transition form factor.

We analyze the pion transition form factor using dispersion theory. We calculate the singly-virtual form factor in the time-like region based on data for the [Formula: see text] cross section, generalizing previous studies on [Formula: see text] decays and [Formula: see text] scattering, and verify our result by comparing to [Formula: see text] data. We perform the analytic continuation to the space-like region, predicting the poorly-constrained space-like transition form factor below [Formula: see text], and extract the slope of the form factor at vanishing momentum transfer [Formula: see text]. We derive the dispersive formalism necessary for the extension of these results to the doubly-virtual case, as required for the pion-pole contribution to hadronic light-by-light scattering in the anomalous magnetic moment of the muon.

Introduction

One of the biggest challenges of contemporary particle physics is the unambiguous identification of signs of beyond-the-standard-model physics. While high-energy experiments are mainly devoted to the search for new particles, high-statistics low-energy experiments can provide such a high precision that standard-model predictions can be seriously scrutinized. A particularly promising candidate for such an enterprise is the gyro-magnetic ratio of the muon; for a review see [1]. Since the muon is an elementary spin-1/2 fermion, the decisive quantity is the deviation of its gyro-magnetic ratio from its classical value. This difference, caused by quantum effects, is denoted by .

From the theory side the potential to isolate effects of physics beyond the standard model is limited by the accuracy of the standard-model prediction. Typically the limiting factor is our incomplete understanding of the non-perturbative sector of the standard model, i.e. the low-energy sector of the strong interaction, which is governed by hadrons as the relevant degrees of freedom instead of the elementary quarks and gluons. In fact, for the hadronic contributions by far dominate the uncertainties for the standard-model prediction. The largest hadronic contribution, hadronic vacuum polarization (HVP), enters at order in the fine-structure constant and can be directly related to one observable quantity, the cross section of the reaction hadrons, by means of dispersion theory. In that way a reliable error estimate of HVP emerges from the knowledge of the experimental uncertainties in the measured cross section. At order there are next-to-leading-order iterations of HVP as well as a new topology, hadronic light-by-light scattering (HLbL) [2]. It was recently shown in [3] that even next-to-next-to-leading-order iterations of HVP are not negligible at the level of accuracy required for the next round of experiments planned at FNAL [4] and J-PARC [5], while an estimate of next-to-leading-order HLbL scattering indicated a larger suppression [6].

With the increasing accuracy of the cross-section measurement for hadrons that can be expected in the near future [7], the largest uncertainty for will then reside in the HLbL contribution. The key quantity here is the coupling of two (real or virtual) photons to any hadronic single- or many-body state. This quantity is not directly related to a single observable. However, it is conceivable to build up the hadronic states starting with the ones most dominant at low energies, in particular the light one- and two-body intermediate states. Based on a dispersive description of the HLbL tensor an initiative has recently been started to relate the one- and two-pion contributions for HLbL scattering to observable quantities [8-10].1

The present work should be understood as an input for this initiative. We focus on the lowest hadronic state, the neutral pion, and its coupling to two (real or virtual) photons (a similar program is currently also being pursued for and ; see [13, 14]). Thus the central object of interest is the pion transition form factor. Its importance for the HLbL contribution to has been stressed early on, see e.g. [1, 15, 16], and triggered many studies of the transition from factor in this context [17-25]. It is defined bywheredenotes the electromagnetic current carried by the quarks and the electric charge of the quark of flavor (in units of the proton charge ).

The normalization of the form factor is given by a low-energy theorem [26-28]. In the chiral limit one findswhich agrees with experiment to a remarkable accuracy; see [29] for a recent review. In (3) denotes the pion decay constant [30].

For the dispersive treatment of the HLbL contribution to as envisaged in [8-10] one needs the pion transition form factor for arbitrary space-like virtualities and of the two photons. We will approach this aim in a multi-step process. In the present work we will formulate the dispersive framework for the general doubly-virtual transition form factor, but restrict the numerical analysis to the singly-virtual case, both in the space- and time-like regions. We will use data on to fix the parameters and predict the cross section for as well as the space-like transition form factor to demonstrate the viability of the approach. While presently low-energy space-like data are scarce [31, 32], new high-statistics data can be expected in the near future from BESIII (see [33, 34]), which makes a calculation of the space-like singly-virtual form factor particularly timely. In a second step, the experimental information from both in space- and time-like kinematics will then serve as additional input for a full analysis of the doubly-virtual form factor.

The basic idea of the dispersive approach for the calculation of the pion transition form factor is its reconstruction from the most important intermediate states in the unitarity relation (see also [34, 35]). At low energies these are the two-pion and three-pion states with isospin and , respectively. Assuming perfect isospin symmetry one of the two photons of the amplitude must be in an isovector and one in an isoscalar state. We shall denote this assignment by the indices and , respectively. Then at low energies the unitarity relation for is dominated by ; see the left diagram in Fig. 1. Additional inelasticities start contributing only at an invariant mass of the isovector photon above , predominantly in the form of four pions, cf. [36]. We will not consider such contributions explicitly in the present work, but estimate their potential impact by variations of the phase shifts in the inelastic region. The crucial building blocks of the dispersive treatment are the charged pion vector form factor , defined byand the amplitude for the reaction. The pion vector form factor with its normalization has been studied in great detail both from the theoretical and experimental side; see e.g. [36-40]. It is closely related to the Omnès function to which we will come back in Sect. 2; see also [41, 42] for more details.

Fig. 1

Two-body unitarity relations for (left) and the amplitude (right). Solid (wiggly) lines denote pions (photons) and the indicates -wave final-state interactions

In contrast, the structure of the amplitude for is much more involved. It will be discussed in detail in Sect. 2. Its two-body unitarity relation, illustrated by the right diagram in Fig. 1, involves the rescattering of pion pairs, which can be resummed in terms of the -wave phase shift within the dispersive approach. While two-body unitarity is exact, we do not consider full three-body unitarity as required by the intermediate states in ; see left diagram in Fig. 2. However, with two-body unitarity fully implemented, the rescattering in generates topologies such as the one shown in the right diagram in Fig. 2, which manifestly contains three-pion cuts. The part of this diagram indicated by the dashed box can be interpreted as a special case of the full amplitude. Therefore, in our framework the structure of the left-hand cut in is approximated by pion-pole terms.

Fig. 2

Three-body unitarity relation for (left) and the approximation inherent in our formalism (right)

The rest of the paper is organized as follows: in Sect. 2 we describe our framework for the determination of the amplitude. In Sect. 3 we formulate the general dispersion relation for the pion transition form factor with arbitrary virtualities for the two photons. In Sect. 4 we specialize the general framework to the case of one on-shell and one time-like photon. As a first application we will determine the cross section of the reaction and compare to the corresponding experimental results. Section 5 is devoted to the analytic continuation into the space-like region as well as the calculation of the slope of the form factor at zero momentum transfer. The Dalitz decay region is discussed in Sect. 6. We close with a summary and outlook in Sect. 7. An appendix is added to discuss the comparison of our results to the simple vector-meson-dominance picture.

The amplitude

Formalism

A key ingredient for the dispersive calculation of the pion transition form factor is the amplitude for the reaction . We definewith , , , , and .

The low-energy limit of is dictated by the chiral anomaly. In the chiral limit this leads to the identification [42-47]A comment is in order to which extent the chiral predictions (3) and (6) have been confronted with experiment so far. has been tested up to in Primakoff measurements of  [48] including chiral [49, 50] and radiative [51] corrections, the former up to two-loop order [52]. Both the world average [30] and the PrimEx result [48] are fully consistent with the chiral tree-level prediction (3), the former even at accuracy, while chiral corrections predict an increase of up to mainly due to mixing [52], in slight tension with the world average. Here, we use (3) directly, given that apart from the very low-energy region the associated uncertainties are sub-dominant.

In contrast to this high accuracy the extractions of both from Primakoff measurements [53] (with chiral and radiative corrections from [54-56]) and  [57] presently allow a test at the level only. In [42] a dispersive framework (see also [55, 58, 59] for earlier work in this direction) was presented that provides a two-parameter description of the cross section valid up to . This opens the possibility to profit from the high-statistics Primakoff data currently analyzed at COMPASS [60] concerning the extraction of to higher accuracy.

We decompose asThis decomposition neglects discontinuities in - and higher partial waves; see [55]. Using the (-channel) partial-wave decompositionwith the Källén function and , we find that the function in (7) is related to the -wave amplitude according to [61]withNote that for positive the evaluation of (8) is straightforward, while some care is needed for the proper analytic continuation of the square roots for negative . Therefore the framework presented here can be immediately applied for instance to the singly-virtual time-like transition form factor, as will be shown in Sect. 4. For the corresponding space-like form factor, to be tackled in Sect. 5, we will refrain from an analytic continuation of the formulas presented here but instead use a dispersion relation to determine the space-like transition form factor from the imaginary part of the time-like one.

For fixed , the quantity , given in (9), only has a right-hand cut starting at . The left-hand cut of the partial wave entirely resides in . Furthermore, the amplitude develops a three-pion cut for , i.e. in kinematics allowing for the physical decay . In this situation, the right- and left-hand cuts in begin to overlap, which leads to a significant complication of the analytic structure, see the corresponding discussion in [61].

The discontinuity of the partial wave along the right-hand cut is given bywhere is the -wave phase shift. Noting that along the right-hand cut, we can recast this relation into the formA once-subtracted dispersive representation solving (12) is given by [61]whereis the Omnès function [62].

An important property of (13) concerns its linearity in the subtraction function , which follows from the fact that is defined in terms of the angular average of itself (9). In this way, takes the role of a normalization, so that in practice (9) and (13) are solved by iteration for , while the full solution is recovered by multiplying with in the end. However, since as a function of implicitly depends on , the subtraction function is not the only source of dependence in the full solution.

For fixed virtualities the solutions of (9) and (13) have been studied in [61] to describe the vector-meson decays .2 In this case the respective subtraction constant is fixed by the overall normalization of the Dalitz plot distribution and hence the corresponding partial decay width. The main complication when extending (13) to arbitrary virtualities of the incoming photon arises from the fact that depends on , a dependence that cannot be predicted within the dispersive framework itself, but has to be determined by different methods. Physically, contains the information how the isoscalar photon couples to hadrons. At low energies, this coupling is dominated by the three-pion state and can be accessed in . For the extraction of we need a representation that preserves analyticity and accounts for the phenomenological finding that the three-pion state is strongly correlated to the very narrow and resonances. We takewith modeled using two relativistic Breit–Wigner functions,In the following we refer to as the spectral function. In (16) is the energy-dependent width of the meson, respectively. We take into account the main decay channels of and viawhere denotes the measured partial decay width for the decay , while the energy-dependent coefficients are given byand the calculation of is performed along the lines described in [61]. For completeness we also include the decay channel of the , which strictly speaking corresponds to a radiative correction. As a consequence the threshold in (15) is actually instead of . However, we checked that as expected the impact of the channel is very small numerically.

The representation (15) can be understood as a dispersively improved Breit–Wigner parametrization [54, 64]: the reconstruction of the real part via a dispersive integral ensures a reasonable behavior of the phase of despite the energy dependence of the widths. We decide to subtract (15) twice: the first subtraction constant is fixed by the chiral anomaly for at the real-photon point (corrected for quark-mass renormalization) [42, 54],The second subtraction serves as an additional background term and is fitted to cross-section data, together with the residues and . Note that the precise form of the spectral function in (16) is irrelevant: the only requirement is to have an analytically rigorous representation of the cross section.

Finally, we give the explicit relation between the amplitude (7) and the cross section (neglecting the electron mass)withandas well as integration boundariesandWe note in passing that for fixed, but arbitrary we can predict the shape of the two-fold differential distribution (21). The knowledge of is only needed for the overall normalization, not for the and dependence.

It has been noted in [61] that the amplitude representation (13) is not accurate enough to give a statistically valid description of the very precise Dalitz plot determination by the KLOE collaboration [66]. For this purpose, a second subtraction was introduced, leading to the representation(only used for in [61]). Similarly, for a twice-subtracted amplitude representation was envisaged theoretically in [42]. For general , the second subtraction will again be -dependent. Provided future measurements allow us to determine such a second subtraction both from cross-section data () and from an Dalitz plot (), the three data points—together with —should permit a smooth interpolation of in a representation similar to (15) (with only a single subtraction). In the absence of such additional high-precision data, we will utilize the singly-subtracted representation (13) of the partial wave for the purpose of this study.

Fits to

Before turning to the fit results, we first summarize the various uncertainty estimates that we have performed in the context of our fits to . First of all, in the calculation of we used three different phase shifts, the phases from [67, 68] and a version of [67] that includes the and the resonances in an elastic approximation to try to mimic the possible impact of inelasticities [61]. In addition, we varied the cutoff in the dispersive integral (13) above which asymptotic behavior is assumed between and , see [41].

Next, our representation for is only adequate below , given that above this energy excited states of and may contribute. The isoscalar vector resonances listed in [30] below with a sizable branching fraction are the and the , with masses and widthsTo estimate the effect of these states, we also consider a version of the fits where additional terms for and are included in (16), identical to the expression for the apart from the channel (we assume branching fraction to for and ). In total, we thus have a three- (five-) parameter representation to be fit to data, with free parameters , , (and , ).

The prime source of data below/above are the SND [69, 70] and CMD2 [71, 72]/the BaBar data sets [73], respectively. Restricting the fit (without and ) to the energy region below , we observed that the SND data set can be described with a reduced close to , while the CMD2 scans can only be accommodated with a significantly worse (around ). We also checked if the respective fit reproduced the correct chiral anomaly by including in (15) as another fit parameter. For SND we indeed obtain , while the fit to CMD2 even produces a negative value of .

One explanation for this apparent tension could be provided by the fact that radiative corrections were not treated in exactly the same way in both experiments. Moreover, the CMD2 scans were restricted to a relatively narrow region around the and masses, limiting the sensitivity to the low-energy region (and thus particularly to the chiral anomaly). Such inconsistencies in the data base were already observed in [74] in the context of the HVP contribution to , where the channel entered with a global reduced of . For the present study we will therefore consider two data sets: first, SND+BaBar and, second, the compilation from [74], in the following denoted by HLMNT. It includes all data sets mentioned so far as well as some older experiments [75-80]. The rationale for doing so is that for the reasons explained above SND/BaBar appear to be the most comprehensive single data sets for low/high energies. Confronting the outcome of fits to the combination of both and to the comprehensive data compilation of [74] should allow for a reasonable estimate of the impact of the uncertainties in the cross section on the prediction for the pion transition form factor.

The result of the three-parameter fit to SND+BaBar below is shown in the left panel of Fig. 3, with fit parameters summarized in Table 1. Since the fits to are hardly distinguishable visually, we only show the curves for the phase shift from [67] and , but give the ranges for the fit parameters found in the full calculation. For these data sets and energy region the reduced is very close to . As alluded to above, the deteriorates substantially when fitting to the full data base of [74], but the central values of the fit parameters remain largely unaffected.

Fig. 3

Fit to the cross-section data of [69, 70] and [73] below (left) and (right), with phase shift from [67] and . The small inserts amplify the regions around the and resonance peaks. Only the fit in the right panel includes in the spectral function. The dashed line indicates the outcome of the fit to the data base of [74]

Table 1

Fit parameters and reduced for the fits to SND+BaBar [69, 70, 73] and HLMNT [74] as described in the main text. The ranges indicate the variation found for the different phase shifts and values of

	β[GeV^-5]	c_ω[GeV^-1]	c_ϕ[GeV^-1]	c_ω^′[GeV^-1]	c_ω^′′[GeV^-1]	χ^2/dof
SND+BaBar, 1.1GeV	5.94…6.21	2.88…2.90	-(0.392…0.406)	–	–	1.01…1.04
HLMNT, 1.1GeV	5.92…6.18	2.81…2.83	-(0.374…0.387)	–	–	6.33…6.36
SND+BaBar, 1.8GeV	7.73…7.78	2.92…2.95	-(0.386…0.400)	-(0.27…0.43)	-(0.70…1.22)	3.18…3.48
HLMNT, 1.8GeV	7.78…7.82	2.88…2.90	-(0.366…0.378)	-(0.19…0.32)	-(0.53…1.02)	7.28…7.62

Extending the fit to higher energies by including and in the spectral function yields a reasonable fit up to , at the expense of a slight deterioration of the data description between the and ; see the right panel of Fig. 3 and Table 1. Again, we observe that the fit result is relatively insensitive to the data set chosen, with larger differences evolving in the region. We will use the outcome of this extended fit to estimate the impact of the high-energy region on the analytic continuation of the transition form factor into the space-like region in Sect. 5.

Dispersion relations for the doubly-virtual transition form factor

We decompose the pion transition form factor into definite isospin components according towhere the first/second index refers to isovector () and isoscalar () quantum numbers of the photon with momentum /. For fixed isoscalar virtuality we can write a once-subtracted dispersion relation in the isovector virtuality [42]where , and is the pion vector form factor (4). Assuming both and to asymptotically fall off like  [41, 61, 81–83] (for fixed ), there is a sum rule for the subtraction function in (28),This sum rule formally converges only with a partial wave based on the singly-subtracted representation (13), with a second subtraction (25) it can at best be evaluated below a certain cutoff. The representation (28) as well as the sum rule (29) have been employed before: for , they yield the vector-meson transition form factors for , including (from the sum rule) the normalization for the real-photon decays [41]. For , one obtains the isovector part of the singly-virtual transition form factor, with the sum rule yielding  [42]. Numerically, these sum rules were found to be saturated at the level [41, 42].

Taken together, (28) and (29) are equivalent to an unsubtracted dispersion relationWe can perform a (necessarily less explicit) subtraction of (28) in as well, defining a subtracted partial waveThe alternative formulation of the dispersive representation, making use of the sum rule (29), then reads

Time-like form factor and

We now specialize the general expressions (27) and (30) to the singly-virtual case for further phenomenological investigation. The transition form factor can be written out explicitly according toHere we have again made use of the sum rule (29) to fix the full transition form factor at to the chiral anomaly . Neglecting the mass of the electron for simplicity, the relation between the cross section and the pion transition form factor is given by

To ensure consistency with the calculation of the amplitude we assume asymptotic behavior of and in (33) above and use a twice-subtracted Omnès representation for (cf. [84])with a radius and the same phase shift as in the respective version of . The isoscalar part, corresponding to the difference in (33), is then calculated by the same methods as in [61] with the normalization fixed from as described in Sect. 2. The isovector part, corresponding to the last term in (33), is completely determined by and can thus be measured in . Here, we use a finite matching point of and fix the normalization to the chiral anomaly [42], but this representation can be improved once the COMPASS data for become available.

Our result for the cross section is shown in Fig. 4. We repeat the calculation for each set of phase shifts and , fitting the isoscalar part in each case both to SND+BaBar and HLMNT. The error band in Fig. 4 represents the uncertainty deduced from scanning over the input quantities in this way. Within uncertainties, the outcome agrees perfectly with the cross section measured by [85-87]. We would like to stress that this result is a prediction solely based on the input quantities described above, most prominently, cross-section data, the -wave phase shift, the pion vector form factor, and the low-energy theorems for and .

Fig. 4

cross section predicted from (left fit with and only, right fit including , ), compared to the data of [85–87]. The inserts again zoom in on the and resonance peaks. The error band represents the variation found by changing the phase shifts, , and the data base as described in the main text

To provide a quantitative measure of the agreement between our result and experiment, we first give the reduced of the mean of our band when comparing to the various data sets; see Table 2. However, the usual does not account for the theory uncertainty, so that it is not surprising that values significantly larger than are obtained. If one assumed the theory band to be statistically distributed with mean values and uncertainties , uncorrelated for each data point , one could consider the difference between theory and experiment with combined error and test the distribution for consistency with zero, leading to a modified ,The corresponding values for this quantity are also summarized in Table 2. Given that in practice correlations between different points of the theory band are not negligible, the statistical interpretation of (36) is not obvious. However, taken together with the observation that curves within the theory band can be constructed with even smaller , it provides quantitative evidence for the consistency of our result with the data. In addition, the comparison of the and for the two fits reveals that, while the fit is deteriorated mostly in the energy region above the , including , improves the agreement with below .

Table 2

Reduced and for the comparison of our result to the data of SND [85, 86] and CMD2 [87] as well as the combined data set. In each case, the upper line refers to the fit with and only, the lower line to the fit including , . and are calculated for all data points below (upper line) and (lower line), respectively

	SND	CMD2	SND+CMD2
χ^2/dof	1.74	4.50	3.12
	1.05	2.37	1.71
χ~^2/dof	0.71	1.42	1.06
	0.56	1.02	0.79

Slope parameter and space-like form factor

We reconstruct the transition form factor in the space-like region again dispersively, making use of the imaginary part determined from the study of the time-like region in the previous sections,If we assume the transition form factor to fulfill even an unsubtracted dispersion relation, this relation implies a sum rule for the chiral anomaly:The slope of the form factor obeysFor the evaluation of these relations we need to specify how to treat the high-energy region of the integrals. Perturbative QCD in the factorization framework of [81] predicts an asymptotic behaviorSince the imaginary part has to vanish at least as fast as the real part, we will assume above a cutoff and estimate the sensitivity to the asymptotic region by varying . We also considered a constant imaginary part above , finding only moderate shifts, but given that such a behavior contradicts [81] we will not include the corresponding variation in the uncertainty bands shown below. Finally, we checked that (37) indeed reproduces the real part in the time-like region, which is non-trivial in view of the imaginary parts generated by three-body cuts in the calculation of the amplitude.

We first turn to the sum rules for and , with results summarized in Table 3. For determined from the fit below , including only and in the spectral function, we find the results given in the first two lines for the slope and the chiral anomaly, respectively. For this fit it does not make sense to increase beyond , given that the fit range in was restricted to this energy region. To estimate the sensitivity to the high-energy region of the dispersive integral, the rest of the table shows the results for the extended fit including in addition and , with three different values for . For each set of parameters we give the ranges corresponding to the variation of the phase shift and as described in Sect. 2. We find very stable results even for the chiral anomaly, whose sum rule is fulfilled at accuracy, although being more sensitive to high energies (it would not converge if we assumed a constant behavior for the imaginary part above ). Averaging over the various fits and data sets we obtain for the slope parameterwhere the error includes the uncertainties from the phase shift, the cutoffs and , the data sets, and the high-energy contribution to the sum rule (estimated via the fits). Our result is appreciably more precise than the value quoted in [30], which is dominated by a monopole fit to the CELLO data [31], or an extraction from an even wider range of space-like data using Padé approximants,  [23].

Table 3

Slope parameter and chiral anomaly from the sum rules (38) and (39). For each fit and data set the upper line refers to the slope in units of , while the lower line gives the sum-rule value for normalized to (3). The ranges correspond to the uncertainty due to the phase shift and

	SND+BaBar	HLMNT
Fit below 1.1GeV	30.4…31.2	30.1…30.9
Λ_π^0=1.1GeV	0.989…1.021	0.976…1.008
Fit below 1.8GeV	30.6…31.4	30.4…31.2
Λ_π^0=1.1GeV	0.992…1.026	0.985…1.019
Fit below 1.8GeV	30.4…31.2	30.3…31.1
Λ_π^0=1.4GeV	0.959…0.987	0.962…0.990
Fit below 1.8GeV	30.3…31.1	30.2…31.0
Λ_π^0=1.8GeV	0.944…0.966	0.947…0.970

Along the same lines, we can also determine the next term in the expansion around ,again with a smaller uncertainty than e.g. from [23]. For a comparison of these numbers to the prediction of vector-meson dominance [88], see 1.

Finally, we use (37) to perform the analytic continuation into the space-like region; see Fig. 5. We follow the convention of the experimental publications to plot . In the case of the CELLO data [31], provided in the original paper for the form factor without the additional factor of , we use the averages given for each bin in the conversion. We also follow the convention to depict the error of the form factor only, and not to propagate an additional uncertainty from the bin size.

Fig. 5

Singly-virtual pion transition form factor in the space-like region, compared to CELLO [31] and CLEO [32] data

As expected, our prediction for the space-like form factor is very accurate at low energies (better than for ), while the uncertainties become more sizable above , reflecting the limited energy range used as input for the time-like calculation. The corresponding error band shown in Fig. 5 comprises the same uncertainty estimates already discussed in the context of the slope parameter (the energy region , which is not reliably described any more in the time-like region, is indicated by the dashed lines in Fig. 5). At low energies the error band is dominated by the variation in the phase shift and ,3 whereas above the treatment of the high-energy region in the dispersive integral becomes increasingly important. The resulting curve is consistent with the existing data base, and will soon be tested by the forthcoming high-statistics low-energy data from BESIII.

Dalitz decay region

So far we have not discussed the third kinematically accessible region of the singly-virtual transition form factor besides and , i.e. the region of the Dalitz decay with , where denotes the electron mass. It is common practice to normalize the corresponding partial decay width to the two-photon decay. The normalized differential decay width is given by [89]Absent high-quality data for this differential decay width we just present our result for the integrated one. In this region of very low momenta it is sufficient to use a polynomial approximation for the transition form factor,Using (41) and (42) the result isin excellent agreement with the experimental value [30]Value and uncertainty in (45) only reflect our form factor calculation and disregard the issue of radiative corrections [90]. The impact of the quadratic term is in the last digit in (45). Note that a pure QED calculation without any form factor yieldsso that the impact of the transition form factor on the integrated decay width is on the level of . High-precision data for the differential decay width (43) will soon become available in the context of dark-photon searches in at NA48/2 [34], but due to the limited sensitivity to the form factor will not improve the PDG value for the slope.

Summary and outlook

We presented the dispersive formalism to analyze the general doubly-virtual pion transition form factor. This includes all effects from elastic rescattering exactly through the respective phase shifts. To determine the isoscalar part that is dominated by intermediate states, we used data on . Furthermore, chiral low-energy theorems on the anomalies and were implemented. As a first step, we carried out the phenomenological analysis of the singly-virtual case. We performed a detailed error analysis and verified our calculation in the time-like region by comparing to data for , yielding very good agreement between theory and experiment. As further applications of the framework, we provided a precise value for the slope parameter, , as well as for the curvature term, . Finally, analytic continuation allowed for a prediction for the transition form factor in the low-energy space-like region that should be compared to the upcoming precise BESIII data.

To extend the calculation to higher energies requires additional input. One could for instance match to the predictions of quark counting rules [81], Regge theory [35], or light-cone sum rules [91, 92]. In the time-like region, with consistency between and demonstrated, one could also fit simultaneously to both reactions to potentially decrease the uncertainties. The most important future extension will concern the generalization to the doubly-virtual case. This can be applied to predict the leptonic neutral pion decay , but most importantly, will help pin down the pion-pole contribution to hadronic light-by-light scattering in . Work in this direction is in progress.