Hadronic vacuum polarization and vector-meson resonance parameters from e + e - → π 0 γ.

We study the reaction e + e - → π 0 γ based on a dispersive representation of the underlying π 0 → γ γ ∗ transition form factor. As a first application, we evaluate the contribution of the π 0 γ channel to the hadronic-vacuum-polarization correction to the anomalous magnetic moment of the muon. We find a μ π 0 γ | ≤ 1.35 GeV = 43.8 ( 6 ) × 10 - 11 , in line with evaluations from the direct integration of the data. Second, our fit determines the resonance parameters of ω and ϕ . We observe good agreement with the e + e - → 3 π channel, explaining a previous tension in the ω mass between π 0 γ and 3 π by an unphysical phase in the fit function. Combining both channels we find M ¯ ω = 782.736 ( 24 ) MeV and M ¯ ϕ = 1019.457 ( 20 ) MeV for the masses including vacuum-polarization corrections. The ϕ mass agrees perfectly with the PDG average, which is dominated by determinations from the K ¯ K channel, demonstrating consistency with 3 π and π 0 γ . For the ω mass, our result is consistent but more precise, exacerbating tensions with the ω mass extracted via isospin-breaking effects from the 2 π channel.

Introduction

The vector mesons and are narrow states compared to other hadronic resonances in the low-energy QCD spectrum. In the case of the , this is because two-body decays are either forbidden by G parity () or require electromagnetic interactions (, ), so that the dominant decay proceeds into . In contrast, for the a G-parity conserving two-body decay into is possible, but suppressed by very small phase space, while the decay into is small due to the Okubo–Zweig–Iizuka rule [1-3]. Accordingly, the most precise information on the mass of the comes from [4-8], which indeed dominates the PDG average [9]. For the determination of the mass, the reaction is the primary source of information [5, 10], but here the three-particle nature of the decay complicates a reliable extraction of the resonance parameters. In particular, there is a significant tension with the mass determination from [11], which together with [12] leads to a scale factor in the PDG average. In this work, we consider the reaction using a dispersive representation of the transition form factor (TFF), which together with our previous work on the channel [13] allows us to present a combined determination of the and resonance parameters within the same framework consistent with the constraints from analyticity, unitarity, and crossing symmetry as well as low-energy theorems.

These constraints, as incorporated in the dispersive representation of the TFF [14-16], are not only valuable for a reliable extraction of resonance parameters, but also define a global fit function for the cross section that allows one to check the consistency of the data sets with these general principles. Applications to the [17-20] and [13] channels have provided such analyses for the two dominant channels in the hadronic-vacuum-polarization (HVP) contribution to the anomalous magnetic moment of the muon . Here, we will study the channel in the same spirit. Since the total contribution is about an order of magnitude smaller than the one of the channel, very large relative changes would be required to notably influence the Standard Model prediction [13, 15–17, 19–39] and thus the tension with the BNL measurement [40]. However, in view of recent results from lattice QCD [41] that suggest large modifications of the hadronic cross section at low energies [42, 43], any further corroboration of the phenomenological HVP evaluation, especially for the channels relevant below such as , is certainly worthwhile – in anticipation of improved measurements at Fermilab [44] and J-PARC [45].

The paper is organized as follows: in Sect. 2 we review the dispersive formalism for the pion TFF and the cross section, which is then applied in Sect. 3 to fit the available data sets. In Sect. 4 we discuss the consequences for the HVP contribution to , in Sect. 5 the combined analysis of the and resonance parameters from and . We close with a summary in Sect. 6.

Time-like pion transition form factor and cross section

Based on the unitarity relation and its crucial building blocks, a once-subtracted dispersive representation for the time-like singly-virtual TFF was constructed in [14],where , is the pion vector form factor, and is the partial-wave amplitude for [14, 46–48], as a generalization of previous studies of the TFFs [49, 50]. In particular, was studied in [14] as a first step towards the doubly-virtual space-like TFF [15, 16], which determines the strength of the pion-pole contribution in a dispersive approach to hadronic light-by-light scattering [51-54], to demonstrate the consistency between and data. Similarly, the and TFFs become relevant for the description of the left-hand cuts in the two-pion contributions [55-60].

denotes the normalization at , as determined at leading order by the Wess–Zumino–Witten anomaly [61, 62]This value, obtained from the pion decay constant [9], agrees with the recent PrimEx-II measurement of the neutral-pion life time [63], which implies . The relation between the cross section and the pion TFF, calculated from the dispersion relation (1), readswhere and we neglected the mass of the electron. Strictly speaking, the dispersion relation (1) applies to the pure QCD process without further radiative correction, so that (3) describes the bare cross section excluding vacuum-polarization (VP) corrections. Accordingly, the mass parameters for and extracted from the fit do not include these VP corrections, in contrast to the PDG convention, see Sect. 5. We use the VP routine from [27] to remove VP from the experimental cross sections.

The isoscalar contribution, corresponding to in the integrand of (1), was calculated in [14] using the previously determined partial wave , where the normalization function was fixed from a fit to data; the isovector part, the last term in (1), was determined using a finite matching point of and a normalization at fixed to the chiral anomaly for the amplitude [64-66]. We will implement the same constraint here, i.e., including quark-mass corrections [47, 67]We stress that in contrast to , whose anomaly-constraint (2) has been confirmed by PrimEx-II at the level of , the chiral prediction for has only been tested experimentally with precision, from Primakoff measurements [68] and [69]. In the remainder of this paper, we assume that follows the precedent, so that the remaining uncertainty in (4), from the quark-mass renormalization, becomes subleading compared to other sources of systematic uncertainty in the dispersive representation of the TFF. In view of open questions regarding the role of subleading terms in the chiral expansion of the amplitude [32, 70–73], a more stringent test of would be highly desirable, which could be achieved with data on taken in the COMPASS Primakoff program [74], using the dispersive framework proposed in [47, 48].

As already remarked in [14], the normalization function could also be determined by a fit to instead of the channel. We follow this approach in the present work and consider an update of this once-subtracted analysis based on the improved parameterization for developed in [15, 16], including a conformal polynomial to be able to describe the inelastic effects that were found to be relevant in above the resonance [13]. For the details of the calculation of we refer to [13, 15, 16], but reiterate the free parameters that enter the dispersive representation for the normalization function : apart from the and resonance parameters, these are their residues and , as well as, potentially, further free parameters in the conformal polynomial. For the evaluation of the final dispersion relation (1), we choose an integration cutoff above which an asymptotic behavior is assumed for both and [75-79]. The isovector part is updated as well in line with the isoscalar contribution.

The systematic uncertainties of the dispersive representation are taken into account as follows: the pion vector form factor is calculated with different variations of the Omnès function [80] using different phase shifts [81, 82] as in [16]; in the meantime, the integration cutoffs in the solution of the Khuri–Treiman equations [83] and in the solution of the pion TFF (1) are varied in the range ; lastly, the asymptotic behavior of the imaginary part of the conformal polynomial is varied as in [13]. The central values of the cross sections are obtained by the best fits to the data sets scanning over the variations of these quantities. The systematic uncertainties are defined as the maximum deviations of all the variations from the central cross sections.

Fits to data

Data sets and normalization uncertainties

In addition to the cross section measurements [11, 84, 85] already included in [14], there are two new data sets, the most accurate new data determined from the whole data sample of the SND experiment [86] and another one that explored a new region between 1.4 and [87]. The full data sets that we consider in our analysis are listed in Table 1. These measurements were performed at the VEPP-2M collider with the SND [84-87] and CMD-2 [11] detectors.

Table 1

Summary of the data sets. For [87] only data points for are included, as the cross section in the region was found to be consistent with zero. In the last column we indicate the size of the systematic errors that we interpret as a normalization-type uncertainty and therefore assume to be correlated

Experiment	Region of s (GeV)	# Data points	Normalization uncertainty
SND 2000 [84]	[0.99, 1.03]	12	3.3%
SND 2003 [85]	[0.60, 0.97]	30	All systematics
SND 2016 [86]	[0.63, 1.35]	60	All systematics
SND 2018 [87]	[1.08, 1.35]	5	All systematics
CMD-2 2005 [11]	[0.60, 1.31]	46	6.0%

As first observed in [88], a naive treatment of normalization-type systematic uncertainties would lead to a bias in the fit. For the data sets in Table 1, the systematic uncertainties of [11, 84] are explicitly given in percentages and therefore interpreted as normalization uncertainties. Likewise, we assume that the systematic uncertainties of [85-87] can be attributed primarily to effects in the same category and thus treat all the systematics uncertainties as correlated. Accordingly, we employ the iterative solution strategy introduced in [89] to treat the normalization uncertainties in a consistent manner and consider both fits with diagonal and full covariance matrices to better monitor the role of the correlations, in analogy to the strategy in [13].

Fits to SND

First, we perform fits to the SND data sets [84-87], with the results shown in Table 2. We display the best results for both the diagonal fit and also the fully correlated one. Only the fit uncertainties are displayed in Table 2 at this step, as we will add the systematic uncertainties of our approach later. Fit errors are already inflated by the scale factorto account for potential inconsistencies between the data sets following the PDG prescription [9].

Table 2

Fits to the combined SND data sets [84–87], for diagonal uncertainties and full covariance matrices. All errors refer to fit uncertainties only

	Diagonal	Full
χ^2/dof	116.9/100	151.3/100
	=1.17	=1.51
p value	0.12	7×10^-4
M_ω(MeV)	782.55(3)	782.58(3)
Γ_ω(MeV)	8.73(7)	8.68(6)
M_ϕ(MeV)	1019.18(5)	1019.18(6)
Γ_ϕ(MeV)	4.24(16)	4.27(17)
c_ω(GeV^-1)	2.95(2)	2.95(3)
c_ϕ(GeV^-1)	-0.378(11)	-0.382(13)
10^4×ξ	3.5(1.3)	4.0(1.0)
10^11×aμπ^0γ|_≤1.35GeV	44.05(24)	44.14(57)

In contrast to [13], we do not include the or other excited vector mesons in the fits since their residues come out consistent with zero, in such a way that their inclusion does not improve the quality of the fit. This strategy is consistent with the observation of a negligible cross section above in [87]. Similarly, the data points above the region are scarce, so that additional free parameters in the conformal polynomial in the parameterization of also do not improve the fits. Therefore, we will use the conformal polynomial to implement the chiral low-energy theorem (with S-wave singularities removed), but do not add additional free parameters.

The accuracy of the center-of-mass energy determination of the data set [85] is worse than the accuracy of the mass value. Therefore, an energy-scale bias was introduced in [85]. A separate fit to [85] indeed produces a smaller mass that is not compatible with the most precise measurement [86]. Therefore, we allow for an energy rescaling for [85],The introduced scaling indeed leads to a considerable improvement of the fits, and its value around comes out in agreement with the energy-bias uncertainties. Similar rescalings within the quoted energy uncertainties were also found to improve the fit quality for the [17] and [13] channels. In the case of , the data set from [85] is the only one for which we see a need for such a rescaling.

We observe that the correlated fit produces larger uncertainties for the parameters and the HVP contribution compared to the diagonal one. Otherwise, the central values of the parameters of both fits are in good agreement within uncertainties. Besides, we find that the correlated fit has a worse description than the diagonal fit, which is a general observation of the iterative fit strategy [89] concerning normalization uncertainties. In fact, this effect may be overestimated here because all systematic uncertainties of [85-87] were assumed to contribute in that category, so that the description could likely be improved if more details on the systematic uncertainties were available. At present, the relatively large of the correlated fit is mainly driven by [87]: a fit to this data set alone gives a and a p value of . The fact that the p value drops by another factor of 3 in the combined SND fit thus points to some minor tensions among [84-87].

Fits to CMD-2

Next, we turn to the fits to the CMD-2 data [11]. Although there is only a single data set, it covers almost the entire relevant energy region. The results are given in Table 3, in the same form as the SND fits, the only exception being the exclusion of the rescaling parameter. For comparison, the fit uncertainties are also inflated by the scale factor (5).

Table 3

Fits to the CMD-2 data set [11]

	Diagonal	Full
χ^2/dof	42.50/40	57.39/40
	=1.06	=1.43
p value	0.36	0.04
M_ω(MeV)	782.53(14)	782.68(9)
Γ_ω(MeV)	8.25(28)	8.41(19)
M_ϕ(MeV)	1019.18(7)	1019.18(6)
Γ_ϕ(MeV)	3.90(21)	3.90(17)
c_ω(GeV^-1)	2.91(7)	2.92(13)
c_ϕ(GeV^-1)	-0.342(13)	-0.341(17)
10^11×aμπ^0γ|_≤1.35GeV	44.88(99)	44.48(3.05)

As for the SND fits, we again find internal consistency for the parameters of the diagonal and the correlated fits. A minor difference concerns the mass and width of the , which display relativity large upward shifts once the correlations are included.

Even once accounting for VP corrections, see Sect. 5, our result for the mass is substantially smaller than in [11], which quotes . A key difference to our formalism is that the vector-meson-dominance ansatz from [11] (see also [90]) permits a complex phase between the and contributions, which cannot be physical because it violates analyticity and unitarity, e.g., by introducing an imaginary part below the respective thresholds. In our fits, we do not see a conflict with the mass extracted from cross sections, and thus conclude that the result from [11] is likely affected by the unphysical phase.

Compared to the SND fits, we observe that the width of the comes out appreciably smaller, albeit with rather large fit uncertainties. This observation will also be reflected in the determination of the width of the in the combined fit presented in the next section.

Combined fits

Finally, our combined SND and CMD-2 fit results are presented in Table 4, including all the data sets listed in Table 1. We take the correlated full fit as our central value, and define our systematic uncertainties as the maximum deviations from the different fit variations discussed in Sect. 2. In all cases, the uncertainties are statistics dominated, in part because a main source of systematic uncertainty from the channel [13], the degree of the conformal polynomial, does not become relevant here given that the observed cross section becomes negligibly small around , with few data points above the resonance.

Table 4

Fits to the combined data sets as shown in Table 1

	Diagonal	Full
χ^2/dof	173.3/146	238.6/146
	=1.19	=1.63
p value	0.06	2×10^-6
M_ω(MeV)	782.55(3)	782.58(3)
Γ_ω(MeV)	8.71(7)	8.65(6)
M_ϕ(MeV)	1019.20(4)	1019.21(4)
Γ_ϕ(MeV)	4.08(13)	4.07(13)
c_ω(GeV^-1)	2.95(2)	2.93(3)
c_ϕ(GeV^-1)	-0.363(9)	-0.358(10)
10^4×ξ	3.5(1.3)	4.1(1.0)
10^11×aμπ^0γ|_≤1.35GeV	44.04(23)	43.82(58)

The combined fit, although dominated by the SND data, reflects some inconsistencies between SND and CMD-2. Most prominently, the downward shift of the width of the in comparison to Table 2 is due to the CMD-2 data [11]. The coupling is also affected and shifted to a smaller value compared to the SND fits. Comparing the residues and to the fit [13], and , we observe reasonable agreement, which indeed is better for than for . Taken together with the fact that also the width from the CMD-2 data drives the combined fit away from the value, we conclude that indeed the interchannel consistency is better for the SND data sets. Figure 1 illustrates our final preferred fit, with close-up views of the and regions in Fig. 2.

Fig. 1

The final fit to the data sets as listed in Table 1 (with VP removed everywhere), where the gray band indicates the full uncertainty and the black band indicates the fit uncertainty

Fig. 2

Fit around the and resonance regions. The black band represents the fit uncertainties, and the gray band indicates the total uncertainty

The final result for the and parameters readswith systematic errors in the second brackets derived as described above. We stress that these resonance parameters do not include VP corrections, see Sect. 5 for a more detailed discussion.

Consequences for the anomalous magnetic moment of the muon

The HVP contribution to the anomalous magnetic moment of the muon reads [91, 92]where the R-ratiois a substitute for the (bare) hadronic cross section and the kernel function is known analytically in terms of the center-of-mass energy s and the muon mass . By convention, the hadrons in the final state of the cross section include photons, so that the channel is actually the first to contribute and sets the integration threshold in (8) to . Based on the fits presented in the previous section our central result for the HVP contribution from the channel becomeswhere the second uncertainty is systematic.1 In comparison to the most recent direct-data-integration analyses, our result is in good agreement with [19], with a slight improvement in the uncertainty thanks to the incorporation of the general QCD constraints. The small difference to [28] partly originates from the application of the trapezoidal rule to scarce data in the tails of the resonance, similarly to the case of . Higher-order interpolations to the data combination of [28] indeed move the HVP contribution towards (10). Our analysis does not support values as low as [93], which is based on a Breit–Wigner description of and . The analysis [19] has updated [26] to account for the threshold contribution , which was already included in [27, 28]. It was determined in [94] based upon a combination of the chiral-anomaly term and -meson dominance [95]. This result is in line with our finding for the threshold region, . Indeed, the agreement between the prediction and the cross section of the first few data points was already observed in [94]. Although these small differences are negligible at the current level of accuracy required for HVP, it is reassuring that the dispersive analysis also corroborates current estimates for the channel, making significant changes in HVP in the energy region up to increasingly unlikely. Other radiative effects beyond , , and the infrared-enhanced contributions in are negligibly small compared to the current uncertainty of the full , see, e.g., [57].

and resonance parameters

Table 6

Comparison of and resonance parameters from to the PDG values, including VP corrections

	e^+e^-→3π,π^0γ	PDG
M¯_ω(MeV)	782.736(24)	782.65(12)
Γ¯_ω(MeV)	8.63(5)	8.49(8)
M¯_ϕ(MeV)	1019.457(20)	1019.461(16)
Γ¯_ϕ(MeV)	4.22(5)	4.249(13)

Our final results for the and resonance parameters as determined from are contrasted to the results from [13] in Table 5. There is good agreement throughout, leading to the combination in the last column. Since the channel is statistics-dominated for all quantities, see (7), the combination is straightforward despite the fact that the systematic errors related to the dispersive representation are correlated. Likewise, the statistical correlations among the resonance parameters (and with the residues) from the respective fits have a negligible impact on the combination. and require a small scale factor (defined in accordance with the PDG conventions [9]). The slight tension for can be traced back to the CMD-2 data set [11], see Sect. 3.3. However, we conclude that within uncertainties the and channels yield a consistent picture for the and resonance parameters.

Table 5

and resonance parameters from [13], (this work), and their combination. The final uncertainties for and include a scale factor . All parameters do not include VP corrections, see Table 6 for the comparison to the PDG parameters

	e^+e^-→3π	e^+e^-→π^0γ	Combination
M_ω(MeV)	782.631(28)	782.584(28)	782.607(23)
Γ_ω(MeV)	8.71(6)	8.65(6)	8.69(4)
M_ϕ(MeV)	1019.196(21)	1019.205(55)	1019.197(20)
Γ_ϕ(MeV)	4.23(4)	4.07(13)	4.22(5)

To be able to compare our results to the PDG conventions, we need to restore the VP corrections that have been removed in the definition of the bare cross sections, which we will denote by a bar over the corresponding quantities. As argued in [13], this leads to the shiftswhere the couplings are related to the respective widths, e.g., , and the uncertainties have been propagated from the PDG values [9] (with potential differences to our determinations being higher-order effects). While otherwise shifts in the widths are negligible, there is an effect enhanced by – mixingwhere we have assigned a generous uncertainty because the estimate relies on a narrow-resonance assumption for the .

The resulting parameters, in comparison to the PDG values, are shown in Table 6. First, one sees that the mass agrees perfectly, with competitive uncertainties. This is an important observation because it demonstrates consistency between and . The latter includes the BaBar measurements [6, 7], which, in contrast to all data sets for considered in this work as well as all the data sets relevant for the and parameters, have not been taken in energy-scan mode (at the VEPP-2M collider), but using initial-state radiation. The width also agrees within uncertainties, but not at the level of accuracy that can be achieved in the channel.

For the mass, its PDG value is dominated by the weighted average of determinations from ( [5], [10]), ( [11]), and ( [12]), where the spread among these determinations drives the scale factor and thus an uncertainty much larger than we obtain from .

As described in Sect. 3.3, we believe that the large value for the mass determined from in [11] originates from an unphysical phase in the vector-meson-dominance model used for the extraction. For the reaction, the uncertainties are more difficult to assess than in the processes because the shape of the background processes is unknown and because the width of the signal, , is dominated by the experimental resolution and much larger than the intrinsic width. Energy scans in , for which the entire amplitude can be reconstructed from general principles and whose energy resolution lies well below the width, should thus yield a much more reliable probe of the resonance parameters.

The mass can also be extracted via – mixing in , and it has been known for a while [96] that without further constraints such fits prefer significantly smaller values for than both the PDG average and our determination from . This conclusion was recently confirmed in [17] within a dispersive approach, leading to , in significant tension with Table 5. However, given the high accuracy required in the channel, additional imaginary parts from the radiative channels , , etc. may actually become relevant [97]. Before their impact is better understood, we would thus consider the mass determination from to be more reliable.

As for the width, our value is consistent with earlier determinations from the channel ( [5], [10]), but lies above the PDG average by . This tension is partly driven by an extraction from the reaction ( [98]), but mostly due to an earlier measurement of by the ND collaboration ( [99]). However, it should be noted that the error quoted in [99] is only statistical, while the modern data sets [5, 10] provide a complete error estimate. Moreover, without access to the original data for from [99] it is impossible to assess its weight in global fits to the data base [13]. In such a situation we do not believe it is adequate to keep the ND measurement in the average for and would therefore consider our determination from modern data sets to be more reliable than the current PDG average.

Summary

We have studied the cross section for in a dispersive framework, which implements constraints from analyticity, unitarity, and crossing symmetry as well as low-energy theorems for the amplitude and the transition form factor for . The relation between this form factor and the cross section forms the basis for the subsequent data analysis.

As the next step, we considered the full data sets for from SND and CMD-2. An iterative fit algorithm was applied to eliminate the D’Agostini bias. Some tensions among different data sets exist and the resulting scale factor of the global fit turns out to be larger compared to those of similar analyses of the and reactions, which in part can be traced back to assumptions necessary for the details of the systematic uncertainties. However, we did not find any data set that needed to be excluded because of severe tensions nor did we identify problematic outliers in the data sets.

As a first application, we evaluated the contribution to HVP, with our central result given in (10). In general, the outcome is in good agreement with analyses using a direct integration of the data, with a slightly reduced uncertainty thanks to the global fit function defined by the dispersive representation. In combination with previous work on and , the three largest channels below have now been subject to scrutiny using constraints from analyticity, unitarity, and low-energy theorems.

Finally, we studied the resulting and resonance parameters first from and then in combination with . Contrary to previous analyses, we find good agreement between the two channels, suggesting that a previous tension could be due to unphysical complex phases in a vector-meson-dominance model employed for the channel. Comparing the combined determinations to the current PDG averages, see Table 6, we observe that for the mass, the value obtained from agrees perfectly at a similar level of precision, demonstrating consistency between extractions from and , the latter dominating the PDG average. The width also comes out consistent, but with larger uncertainty than from the channel. For the , we find that the combination of and determines its mass at a level not far from the mass, and argue that the resulting values both for the mass and the width are more reliable than the current PDG averages. However, the tension with the mass determination from the channel persists, suggesting that an improved understanding of isospin-breaking effects therein will become necessary.