## A new evaluation of the hadronic vacuum polarisation contributions to the muon anomalous magnetic moment and to α ( mZ 2 )

M. Davier , A. Hoecker , B. Malaescu , and Z. Zhang 1 2 3 1

- 1 IJCLab, Universit´ Paris-Saclay et CNRS/IN2P3, Orsay, France e
- 2 CERN, CH-1211, Geneva 23, Switzerland
- 3 LPNHE, Sorbonne Universit´ e, Paris Diderot Sorbonne Paris Cit´, CNRS/IN2P3, Paris, France e

Abstract. We reevaluate the hadronic vacuum polarisation contributions to the muon magnetic anomaly and to the running of the electromagnetic coupling constant at the Z -boson mass. We include newest e + e -→ hadrons cross-section data together with a phenomenological fit of the threshold region in the evaluation of the dispersion integrals. The precision in the individual datasets cannot be fully exploited due to discrepancies that lead to additional systematic uncertainty in particular between BABAR and KLOE data in the dominant π + π -channel. For the muon ( g -2) / 2, we find for the lowest-order hadronic contribution (694.0 ± 4.0) · 10 -10 . The full Standard Model prediction differs by 3.3 σ from the experimental value. The five-quark hadronic contribution to α m ( 2 Z ) is evaluated to be (276.0 ± 1.0) · 10 -4 .

## 1 Introduction

The Standard Model (SM) predictions of the anomalous magnetic moment of the muon, a µ = ( g µ -2) / 2, with g µ the muon gyromagnetic factor, and of the running electromagnetic coupling constant, α s ( ), an important ingredient of electroweak theory, are limited in precision by hadronic vacuum polarisation (HVP) contributions. The dominant hadronic terms can be calculated with the use of experimental cross-section data, involving e + e -annihilation to hadrons, and perturbative QCD to evaluate energy-squared dispersion integrals ranging from the π γ 0 threshold to infinity. The kernels occurring in these integrals emphasise low photon virtualities, owing to the 1 /s descent of the cross section, and, in case of a µ , to an additional 1 /s suppression. About 73% of the lowest order hadronic contribution to a µ and 58% of the total uncertainty-squared are given by the π + π -( γ ) final state, 1 while this channel amounts to only 12% of the hadronic contribution to α s ( ) at s = m 2 Z [1].

In this work, we reevaluate the lowest-order hadronic contribution, a had,LO µ , to the muon magnetic anomaly, and the hadronic contribution, ∆α had ( m 2 Z ), to the running α m ( 2 Z ) at the Z -boson mass using newest e + e -→ hadrons cross-section data and updated techniques. In particular, we perform a phenomenological fit to supplement less precise data in the low-energy domain up to 0.6 GeV. We also reconsider the systematic uncertainty in the π + π -channel in view of discrepancies among the most precise datasets.

1 Throughout this paper, final state photon radiation is implied for all hadronic final states.

All the experimental contributions are evaluated using the software package HVPTools [2]. To these are added narrow resonance contributions evaluated analytically, and continuum contributions computed using perturbative QCD.

## 2 Combination of experimental inputs

The integration of data points belonging to different experiments with their own data densities requires a careful treatment to avoid biases and to properly account for correlated systematic uncertainties within the same experiment and between different experiments, as well as within and between different channels. Quadratic interpolation (splines) of adjacent data points is performed for each experiment, and a local combination in form of a weighted average of the interpolations is computed in bins of 1 MeV, or in narrower bins for the ω and φ resonances.

The uncertainties on the combined dataset, the data integration and the phenomenological fit are computed using large numbers of pseudo-experiments. These are generated taking into account all measurement uncertainties and their correlations. While this treatment guarantees a proper propagation of uncertainties, the resulting precision of the combination still depends on the chosen test statistic: a poor choice (e.g., an arithmetic instead of a weighted average) would lead to poor precision, while an aggressive choice (e.g., exploiting the available correlation information globally over the full spectrum, thereby benefiting from constraints among

different energy regimes ) could lead to an optimistic 2 precision claim with the risk of undercoverage with respect to the (unknown) truth. To avoid either case, we employ a test statistic that only relies on local measurement uncertainties and correlations to combine datasets in a given bin. 3 As stated above, the uncertainty in each combined bin and the correlation among bins are evaluated using pseudo-experiments generated with the full correlation information. Correlations between channels are accounted for by propagating individually the common systematic uncertainties. 4

Where results from different datasets are locally inconsistent, the combined uncertainty is rescaled according to the local χ 2 value and number of degrees of freedom following the PDG prescription [6]. Such inconsistencies are currently limiting the precision of the combination in the dominant π + π -channel as well as in the K K + -channel (see discussions below). In most exclusive channels the largest weight in the combination is provided by BABAR data.

Closure tests with known distributions have been performed in the dominant π + π -channel to validate both the combination and integration procedures.

## 3 Input data

Exclusive bare hadronic cross-section measurements for 32 channels are integrated up to 1.8 GeV over the rele-

2 Systematic uncertainties are based on estimates which are impacted by imponderables regarding size and correlation among measurements, in particular uncertainties due to theoretical modelling. Systematic uncertainties are often evaluated in relatively wide mass ranges, the event topology may evolve between measurements performed at different centre-of-mass energies (affecting for example the acceptance and tracking efficiency) as does the background composition, systematic uncertainties due to trigger and tracking may be correlated, etc. It is therefore important to treat systematic uncertainties and their correlations with care and avoid the use of long-range correlations to constrain measurements among different centre-of-mass energies. Ambiguities in systematic uncertainties and their correlations have been studied in other experimental areas and different 'configurations'/'scenarios' of uncertainties were proposed [3-5].

3 This information on the uncertainties and correlations is used on slightly wider ranges, of typically up to a couple of 100 MeV, when averaging regions are defined in order to account for the difference between the point-spacing and bin-sizes for the various experiments [2]. In this procedure the systematic uncertainties are not constrained, but rather directly propagated from each input measurement to the averaging regions and then to the fine bins.

- 4 A number of 15 such uncertainties are accounted for in the current study. Typical examples are the luminosity uncertainties, if the data stem from the same experimental facility but measure different channels, and uncertainties related to radiative corrections.

vant dispersion kernels. This analysis uses all the available public data with recent additions [7-16]. References for data already included in our 2017 analysis are provided in the corresponding paper [1] as well as earlier publications [17, 18].

In the energy range 1.8-3.7 GeV and above 5 GeV four-loop perturbative QCD is used [19]. The contributions from the open charm pair production region between 3.7 and 5 GeV are again computed using experimental data. For the narrow resonances J/ψ and ψ (2 S ) Breit-Wigner line shapes are integrated using their currently best known parameters [6].

The following discussion of individual channels focuses on the HVP contribution to a µ as it relies more strongly on the low-energy experimental data. We mainly explore the impact of the data released since our last publications [1, 17]. If not stated otherwise, all numerical results for a µ are quoted in units of 10 -10 .

## 3.1 The π + π -channel

Data from the BABAR [20, 21] and KLOE [22-24] experiments dominate the measurement of the π + π -channel. Their sub-percent precision is not matched by the other experiments (CMD-2, SND, and BESIII ). 5 New data in this channel stem from CLEO [10] using large angle initial state radiation (ISR) and taking into account up to one additional photon, following the BABAR method [21]. Relatively large statistical uncertainties and a systematic uncertainty of 1.5% are, however, insufficient to improve the precision of the combined π + π -contribution. Recently, a combination of the three KLOE measurements was proposed [26]. We do not use this combination as the KLOE measurements correspond to different ISR topologies and normalisation procedures. 6 Eigenvector decomposition of the statistical and systematic covariance matrices of the three most recent series of KLOE measurements [26] is used. Each eigenvector multiplied by the square-root of the corresponding eigenvalue is treated as an uncertainty source that is fully correlated between the KLOE data points, while the individual sources are treated as independent among each other. Pseudo-experiments are generated in the usual way to propagate correlated uncertainties among the KLOE measurements.

5 There is a small inconsistency between the bin-by-bin statistical uncertainties and the diagonal values of the statistical covariance matrix of the π + π -data published by BESIII [25].

6 Using the KLOE combination [26] we find for a had,LO µ [ ππ ] between the π + π -threshold and 1.8 GeV a value of 506.6 ± 2.4, which is to be compared with 506.7 ± 2.3 as obtained from the HVPTools combination. For both calculations the usual local √ χ / 2 ndof uncertainty rescaling method was applied. Without the rescaling the corresponding results are 506.6 ± 2.0 and 506.7 ± 2.0, respectively. The similarity of the results is maintained when using a phenomenological fit up to 0.6 GeV (see later in text).

Fig. 1. Bare cross section of e + e -→ π + π -versus centre-of-mass energy for different energy ranges. The error bars of the data points include statistical and systematic uncertainties added in quadrature. The green band shows the HVPTools combination within its 1 σ uncertainty.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-3486ec46-8b55-4335-aad2-044d168f274a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-83caf27c-4060-4bb6-b2e8-b3ba10e5193d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-3f01808e-59bf-49a1-91cf-90d087e94edd.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-e4cf6a7a-335d-44e1-81f8-93ddd4ae7b07.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-49ce5c24-ffb8-417d-8b04-ac91600a8269.jpg)

Fig. 2. Comparison between individual e + e -→ π + π -cross-section measurements from BABAR [20, 21], KLOE 08 [22], KLOE10 [23], KLOE12 [24], BESIII [25], CLEO [10], CMD-2 03 [27], CMD-2 06 [28], SND [29], and the HVPTools combination. The error bars include statistical and systematic uncertainties added in quadrature.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-5846f3bd-2b72-4b32-abaa-3ca5b2144977.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-5c7a00b2-192e-4bbd-88c4-33d85bfd485b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-328b0c9f-dc2d-4ad3-a058-8373028f3777.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-89985c65-1a7b-4207-8b85-1328816f71aa.jpg)

Fig. 3. Left: relative local weight per experiment contributing to the e + e -→ π + π -cross-section combination versus centre-of-mass energy. Right: local scale factor versus centre-of-mass energy applied to the combined π + π -cross-section uncertainty to account for inconsistency in the individual measurements.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-0b66185f-3ecd-4588-b2e3-007b6056f5f4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-a5305711-c370-46ad-a681-a86147bc90d3.jpg)

Figure 1 shows the available e + e -→ π + π -crosssection measurements in various panels zooming into different energy ranges. The green band indicates the HVPTools combination within its 1 σ uncertainty. Comparisons between the combination and the most precise individual measurements are plotted in Fig. 2. Figure 3 (left) shows the local combination weight versus √ s per experiment. The BABAR and KLOE measurements dominate over the entire energy range. Owing to the sharp radiator function, the KLOE event yield increases towards the φ (1020) mass leading to a better precision than BABAR in the 0.8 -1.0 GeV region. The group of experiments labelled 'Other exp' in the left panel of Fig. 3 corresponds to older data with incomplete radiative corrections. Their weights are small throughout the entire energy domain. The right hand panel of Fig. 3 shows the scale factor versus centreof-mass energy that is locally applied to the combined π + π -cross-section uncertainty to account for inconsistencies among the individual measurements. Significant inconsistencies are found between the most precise BABAR and KLOE datasets.

The computation of the dispersion integral over the full π + π -spectrum requires to extend the available data to the region between threshold and 0.3 GeV, for which we use a fit as described below.

## Phenomenological fit

The bare e + e -→ π + π -annihilation cross section is related to the pion form factor F 0 π ( s ) (excluding vacuum polarisation) by

$$\sigma ^ { ( 0 ) } ( e ^ { + } e ^ { - } & \to \pi ^ { + } \pi ^ { - } ) = & \quad \text{in tl} \\ \frac { \pi \alpha ^ { 2 } } { 3 s } \beta _ { 0 } ^ { 3 } ( s ) \cdot | F _ { \pi } ^ { 0 } ( s ) | ^ { 2 } \cdot F S R ( s ) \,, \quad ( 1 ) \quad \text{note} \\ \text{dep} \text{c}$$

where α is the electromagnetic coupling constant, β 0 ( s ) = √ 1 -4 m /s 2 π is a threshold kinematic factor and FSR( ) is the final state radiation contribution. s

The pion form factor is an analytic function of s in the complex plane, except on the real axis above 4 m 2 π . It can be parameterised as a product of two functions [30]

$$F _ { \pi } ^ { 0 } = G ( s ) \cdot J ( s ) \,, \quad \quad \ \ ( 2 ) \quad \ \ e \alpha$$

where

$$G ( s ) = 1 + \alpha _ { V } s + \frac { \kappa s } { m _ { \omega } ^ { 2 } - s - i m _ { \omega } \Gamma _ { \omega } } \,, \quad ( 3 ) \quad \text{distr} _ { 1 }$$

and, exploiting the unitarity constraint which identifies arg( F 0 π ) with the P-wave π + π -phase shift δ 1 ( s ),

$$\cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdclots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ J ( s ) & = e ^ { 1 - \delta _ { 1 } ( s _ { 0 } ) / \pi } \cdot \left ( 1 - \frac { s _ { 1 } ( s _ { 0 } ) } { \pi } \right ) ^ { \frac { s _ { 0 } } { \pi } } \left ( 1 - \frac { s _ { 0 } } { \pi } \right ) ^ { \frac { s _ { 0 } } { \pi } } \left ( 1 - \frac { s _ { 0 } } { \pi } \right ) ^ { - 1 } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \frac { s _ { 0 } } { \pi } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \frac { s _ { 0 } } { \pi } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \frac { s _ { 0 } } { \pi } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \frac { s _ { 0 } } { \pi } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \frac { s _ { 0 } } { \pi } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \frac { s _ { 0 } } { \pi } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 1 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac { s _ { 0 } \right ) } { \pi } \right ) ^ { \binom { 1 } { 2 } } \left ( \frac < \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

The last term in Eq. (3) accounts for ρ -ω mixing. The function J s ( ) is taken from Refs. [31, 32]. Owing to ρ dominance, the phase shift δ 1 ( s ) can be parameterised by [33]

$$\stackrel { \sim } { \Pr } \quad \begin{array} { c } \lambda \vdots \cdots \\ \text{pre-} \\ \text{ersus} \end{array} \quad \cot \delta _ { 1 } ( s ) = \frac { \sqrt { s } } { 2 k ^ { 3 } ( s ) } \left ( m _ { \rho } ^ { 2 } - s \right ) \left ( \frac { 2 m _ { \pi } ^ { 3 } } { m _ { \rho } ^ { 2 } \sqrt { s } } + B _ { 0 } + B _ { 1 } \omega ( s ) \right ) \\ \text{users} \quad \text{with} \quad \end{array}$$

with

$$k ( s ) = \frac { \sqrt { s - 4 m _ { \pi } ^ { 2 } } } { 2 } \,, \quad \omega ( s ) = \frac { \sqrt { s } - \sqrt { s _ { 0 } - s } } { \sqrt { s } + \sqrt { s _ { 0 } - s } } \,.$$

The six free parameters α V , κ , m ω , m ρ , B 0 and B 1 are determined by the fit to the π + π -data restricted to the region up to 1 GeV to stay below the threshold of significant inelastic channels. The width of the ω resonance is fixed to its PDG value of 8.49 MeV [6], and √ s 0 = 1.05 GeV. The results of the fit are given in Table 1. To derive an estimate for the model uncertainty, we independently vary √ s 0 to 1.3 GeV and remove the linear term B ω s 1 ( ) from Eq. (5) since the resulting value of B 1 from the nominal fit is consistent with zero.

The fit is performed using as test statistic a diagonal χ 2 function that accounts for the statistical and systematic uncertainties of the experimental measurements. 7 The same uncertainty rescaling in case of local discrepancies among datasets as for the HVPTools based combination is applied. Correlations are ignored in the test statistic, but accounted for in the uncertainty propagation through a series of pseudo-experiments for each of which the full fit procedure is repeated. This is a conservative procedure, as exploiting correlations in the test statistic would improve the precision of the fit. Currently, the most precise measurements are dominated by systematic uncertainties, whose size and mass dependence as well as correlation among each other and among data points rely on estimates with somewhat limited precision, as discussed in section 2. Since there are also clear indications of a significant underestimate of the size of uncertainties in the discrepant dataset(s), we prefer not to exploit this information in the constrained fit. Pseudo-experiments are also used to assess the goodness-of-fit on the data, which yields a p-value of 0.27. 8 We have checked the reliability of this procedure by generating a set of pseudo-experiments and evaluating the p-value for each of them. The expected distribution of p-values reconstructed this way is indeed uniform between 0 and 1.

A graphical comparison of the fit result with the data is shown in Fig. 4. In the energy range between 0.3

7 For the fit we use the original data provided by each experiment instead of the HVPTools combination.

8 The p-values for each individual dataset read 0.042 (BABAR), 0.097 (KLOE), 0.449 (CMD), 0.675 (TOF), 0.718 (DM1), 0.756 (CMD-2), 0.796 (SND), and 0.984 (CLEO). The p-values for both OLYA and BESIII are close to 1.

Fig. 4. Same as Fig. 2 except that the comparison is made with respect to the fit instead of the combination. The black and red curves show the results of two alternative fits where the data from KLOE and BABAR, respectively, were excluded.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-9cd9edd6-7263-4785-9bd6-85f878caa142.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-c3f16acc-ec5d-40f6-b899-d66d9f9871a1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-e13eefca-f8e8-4e95-a789-770e012c1c43.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-915e2440-52c1-450b-94d0-2d52ae49ed4c.jpg)

Table 1. Results of the fit to all π + π -data. The diagonal elements give the fitted parameter values and their uncertainties, while the off-diagonal elements give the correlation coefficients.

|             | α V           | κ [10 - 4 ]   | B 0           | B 1           | m ρ [MeV]   | m ω [MeV]   |
|-------------|---------------|---------------|---------------|---------------|-------------|-------------|
| α V         | 0.133 ± 0.020 | 0.52          | - 0.45        | - 0.97        | 0.90        | - 0.25      |
| κ [10 - 4 ] |               | 21.6 ± 0.5    | - 0.33        | - 0.57        | 0.64        | - 0.08      |
| B 0         |               |               | 1.040 ± 0.003 | 0.40          | - 0.40      | 0.29        |
| B 1         |               |               |               | - 0.13 ± 0.11 | - 0.96      | 0.20        |
| m ρ [MeV]   |               |               |               |               | 774.5 ± 0.8 | - 0.17      |
| m ω [MeV]   |               |               |               |               |             | 782.0 ± 0.1 |

and 0.6 GeV, the result of the fit yields for a had,LO µ [ ππ ] a contribution of 109.8 ± 0.4 ± 0.4, where the first error is experimental and the second the model uncertainty. The latter is obtained by adding linearly the absolute values of following two variations: the √ s 0 variation of -0.13 ± 0.10 and the difference of without and with the B ω s 1 ( ) term of 0.24 ± 0.14, where the uncertainty of each variation accounts for the correlation between the integral results. The corresponding result based on data integration is 109.6 ± 1.0. Taking into account the correlation of 72% between the experimental uncertainties, the difference between the two evaluations amounts to

0.2 ± 0.9. Similarly, for ∆α had ( m 2 Z ) the difference is 0.020 ± 0.028. The fit therefore gives compatible but more precise results than the direct data integration.

Other studies using constraints from unitarity and analyticity with the aim to improve the precision of the π + π -HVP contribution to the muon g -2 exist in the literature. 9 The treatment followed in Ref. [35] is similar to ours with, however, a more elaborate theoretical analysis. Differences are also present in the treatment of

9 In Ref. [34] an analyticity-based phenomenological fit has been used for the π + π -π 0 channel.

experimental data, Ref. [35] using a χ 2 computed globally, including correlations across all the experimental data points and bins in the full mass range of interest. However, in order to avoid too low p-values, some bins of the KLOE measurements were removed in that study and energy rescaling parameters were introduced to fit each measured energy spectrum. A different approach is followed in Ref. [36, 37] where the low-mass contribution was obtained from input data at a fixed mass, followed by a simple average used to combine inputs in the √ s range between 0.65 and 0.71 GeV (and then to combine values from different experiments). Instead of a direct evaluation of the correlations from the published information, they were assumed to be the same between all the combination inputs and an attempt was made to evaluate them based on the resulting χ 2 value. It is possible to directly compare the results for the mass region between threshold and 0.63 GeV. In the present analysis a value of 133.2 ± 0.5 ± 0.4 is found, which agrees with the other results, 132.8 ± 0.4 ± 1.0 [35] and 132.9 ± 0.8 [36, 37].

It is also interesting to compare the results given in Table 1 with other analyses. The value obtained for κ corresponds to a branching fraction of ω into π + π -of (2.09 ± 0.09) · 10 -2 , in agreement with the result extracted from the fit of Ref [35], (1.95 ± 0.08) · 10 -2 . Both values disagree with the PDG average [6], (1.51 ± 0.12) · 10 -2 , dominated by the result of Ref. [30] which uses fits to essentially the same data. The fitted ω mass is found to be lower than the PDG average [6] obtained from 3 π decays by (0.65 ± 0.12 ± 0.12 PDG ) MeV, in agreement with previous fits of the ρ -ω interference in the 2 π spectrum (see for instance Refs. [20, 35]).

## The π + π -contribution

The evaluation of the complete a had,LO µ [ ππ ] integral for the π + π -contribution from threshold to 1.8 GeV, using the fit up to 0.6 GeV and the HVPTools data combination above, gives 507.0 ± 1.9. The choice of the ranges is justified by the good agreement between fit and combined data integration in the 0.6-1.0 GeV region with, however, no advantage in precision for the fit. The correlation among the two contributions is found to be 62% using pseudo-experiments.

Removing BABAR or KLOE from the dataset gives 505.1 ± 2.1 and 510.6 ± 2.2, respectively, with an absolute difference of 5.5 that is significantly larger than the individual uncertainties. Figure 5 shows a comparison among the most precise a had,LO µ [ ππ ] evaluations in the interval 0.6-0.9 GeV. The results of all other experiments fall in-between the BABAR and KLOE results, with insufficient precision to resolve the discrepancy.

Figure 6 compares the HVPTools combination and the fits without using the BABAR and KLOE data, respectively, with the fit result for the full dataset. In light of this discrepancy, which is not fully captured by the

]

Fig. 5. Comparison of results for a had,LO µ [ ππ ], evaluated between 0.6 GeV and 0.9 GeV for the various experiments. In case of CMD-2 all available measurements have been combined using HVPTools. For KLOE the result from the public combination [26] is displayed.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-45392570-b850-4525-bc18-e2a12af038e7.jpg)

1

Fig. 6. The HVPTools combination (green band) relative to the result of the fit to all individual π + π -data (blue band) versus centre-of-mass energy. The black and red curves show the results of two alternative fits where the data from KLOE and BABAR, respectively, were excluded.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-db935888-660c-4e6d-835c-21aab200c63c.jpg)

local uncertainty rescaling procedure, we add as additional systematic uncertainty half of the full difference between the complete integrals without BABAR and KLOE, respectively, and we place the central value of the a had,LO µ [ ππ ] contribution half-way between the two results. To avoid double counting, the local uncertainty rescaling between BABAR and KLOE is not applied, but that between these and the other π + π -datasets is kept. This procedure results in a total π + π -contribution of a had,LO µ [ ππ ] = 507.9 ± 0.8 ± 3.2, where the first uncertainty is statistical and the second systematic (dominated by the new uncertainty of 2.8).

## 3.2 The K + K -channel

Tensions among datasets are also present in the K K + -channel (see top panel of Fig. 7 for a display of the available measurements). A discrepancy between BABAR and SND was observed for masses between 1.05 and 1.4 GeV, which has been resolved with the most recent SND result [38] so that older SND data are discarded.

Concerns also arise regarding data on the φ (1020) resonance. Previously, a 5.1% difference between CMD2 at VEPP-2M and BABAR was observed, with the CMD-2 data being lower. New results from CMD-3 at VEPP-2000 [39] exhibit the opposite effect: they are 5.5% higher than BABAR (cf. middle panel in Fig. 7). The discrepancy of almost 11% between the CMD-2 and CMD-3 datasets, which largely exceeds the quoted systematic uncertainty of 2.2%, of which only 1.2% accounts for uncertainties in the detection efficiency, is claimed to originate from a better understanding of the detection efficiency of low-energy kaons in the CMD3 data. 10 Given the yet unresolved situation, we keep both CMD-2 and CMD-3 datasets, which due to the uncertainty rescaling procedure in presence of discrepancies leads to a deterioration of the precision by about a factor of two of the combined data (cf. bottom panel of Fig. 7). 11

## 3.3 Other channels

Recent measurements have been included in the data combinations: π γ 0 from SND [12], π + π -2 π 0 from BABAR [8], π + π -3 π 0 from BABAR [13], ηπ + π -from BABAR [11], ηπ + π -π 0 from CMD-3 [7] and SND [15], φη from CMD-3 [16], and K K π S L 0 from SND [9]. The π γ 0 and π + π -π 0 contributions include small additions of 0.12 ± 0.01 and 0.01 ± 0.00, respectively, to cover the threshold region up to the lowest-energy data measurements [40, 41].

Only very few final states remain to be estimated using isospin symmetry. Already in 2017, a significant step was achieved with the BABAR measurements of all the final states contributing to the KKπ and KKππ channels, so that previous isospin-based estimates became obsolete. Now the only significant (albeit small) contribution obtained with the use of isospin constraints is that for the π + π -4 π 0 channel. The part excluding

10 In comparison with the CMD-2/3 and SND measurements, the ISR method of BABAR benefits from highermomentum kaons with better detection efficiency owing to the boost of the final state.

11 We have verified that the local χ 2 rescaling procedure covers the global discrepancy among the CMD-2 and CMD3 data by removing alternatively one or the other dataset from the K K + -combination. The difference of 0.45 resulting between the two a µ values is covered by the uncertainty rescaling (a similar conclusion is reached for ∆α had ( m 2 Z )). There is therefore no need to introduce an additional global systematic uncertainty as for the π + π -case.

Fig. 7. Top panel: bare cross sections for e + e -→ K K + -. See text for a description of the data used. Middle: comparison between individual e + e -→ K K + -cross-section measurements from BABAR [42], CMD-2 [43], CMD-3 [39] and SND [38], and the HVPTools combination. Bottom: local scale factor versus centre-of-mass energy applied to the combined K K + -cross-section uncertainty to account for inconsistency in the individual measurements.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-a3e7aa27-0c25-44f5-8a06-42da87c60a90.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-b20e046a-5673-45c4-8ee3-ac113e6ebef8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-bac23e84-20c5-42c5-930a-c959d1db9b29.jpg)

η 3 π , which is obtained from measured processes and reevaluated in this analysis, amounts to a fraction of the total dispersion integral of only 0.016% with an assigned systematic uncertainty of 100%. One could also question the completeness of the set of exclusive processes considered below 1.8 GeV, including up to 6-pion and KK +3 pions. A recent measurement of the 3 π + 3 π -π 0 cross section by CMD-3 [14] allows one to estimate a very small total 7-pion contribution, included in this analysis, of only 0.002%. Although such highmultiplicity channels appear to be contributing negligibly below 1.8 GeV their importance is likely to increase above.

All other contributions are identical to the ones described in our previous analysis [1], except for (i) a reevaluation of the contribution from ω decay modes not reconstructed in other exclusive channels, and (ii) a better estimate of the KKπ π π + -0 contribution, excluding φη which is dominated by the KKω final state [44].

## 4 Compilation and results

A compilation of the various contributions to a had,LO µ and to ∆α had ( m 2 Z ), as well as the total results are given in Table 2. The experimental uncertainties are separated into statistical, channel-specific systematic, and common systematic contributions that are correlated with at least one other channel. The contributions from the J/ψ and ψ (2 S ) resonances in Table 2 are obtained by numerically integrating the corresponding undressed 12 Breit-Wigner lineshapes. The uncertainties in the integrals are dominated by the knowledge of the corresponding electronic width Γ R → ee for which we use the values 5.53 ± 0.10 keV for R = J/ψ and 2.34 ± 0.04 keV for R = ψ (2 S ) [6].

Sufficiently far from the quark thresholds we use four-loop [19] perturbative QCD, including O ( α 2 S ) quark mass corrections [45], to compute the inclusive hadronic cross section. Nonperturbative contributions at 1.8 GeV were determined from data [46] and found to be small. The uncertainties of the R QCD contributions given in Table 2 are obtained from the quadratic sum of the uncertainty in α S (we use α S ( m 2 Z ) = 0.1193 ± 0.0028 from the fit to Z precision data [47]), the truncation of the perturbative series (we use the full fourloop contribution as systematic uncertainty), the difference between fixed-order perturbation theory (FOPT)

12 The undressing uses the BABAR programme AFKVAC , correcting for both leptonic and hadronic VP effects. The hadronic part is obtained from a numerical integration over cross section data for the continuum, supplemented by analytical expressions for the contributions of narrow resonances including both their real and imaginary components. The resulting correction factors reduce the J/ψ and ψ (2 S ) contributions to a had,LO µ by about 4% and are known to a precision of better than 10 -3 .

and, so-called, contour-improved perturbation theory (CIPT) [48], as well as quark mass uncertainties (we use the values and uncertainties from Ref. [6]). The former three uncertainties are taken to be fully correlated between the various energy regions (see Table 2), whereas the (smaller) quark-mass uncertainties are taken to be uncorrelated.

To examine the transition region between the sum of exclusive measurements and QCD we have computed a had,LO µ and ∆α had ( m 2 Z ) in the narrow energy interval 1.8-2.0 GeV. For the former quantity we find 7.65 ± 0.31 and 8.30 ± 0.09 for data and QCD, respectively. The full difference of 0.65 (0.28 · 10 -4 in the case of ∆α had ( m 2 Z )) is assigned as additional systematic uncertainty, labelled by 'dual' subscripts in Table 2. It accounts for possible low-mass quark-hadron duality violation effects in the perturbative QCD approximation that we use for this interval to avoid systematic effects due to unmeasured high-multiplicity channels.

Figure 8 shows the total hadronic e + e -annihilation rate R versus centre-of-mass energy as obtained from the sum of exclusive data below 2 GeV and from inclusive data between 1.8 and 5 GeV. 13 Also indicated are the perturbative QCD prediction above 1.5 GeV and the analytical narrow J/ψ and ψ (2 S ) resonances.

## Muon magnetic anomaly

Adding all lowest-order hadronic contributions together gives

$$a _ { \mu } ^ { \text{had,LO} } = 6 9 4. 0 \pm 4. 0 \,,$$

which is dominated by experimental systematic uncertainties (cf. Table 2 for a separation of the total uncertainty into its components), with an uncertainty of 2.8 originating from the BABAR versus KLOE discrepancy in the π + π -channel. The new result is 0.9 units larger than our previous evaluation [1], 693.1 ± 3.4, mostly because we symmetrised the new BABAR/KLOE systematic uncertainty. The total uncertainty is increased by 18%. The result without the additional BABAR/KLOE systematic uncertainty is 693.1 ± 3.2.

Adding to (6) the contributions from higher order hadronic loops, -9.87 ± 0.09 (NLO) and 1.24 ± 0.01 (NNLO) [52], hadronic light-by-light scattering, 10.5 ± 2.6 [53], as well as QED, 11658471.895 ± 0.008 [54] (see also [55] and references therein), and electroweak effects, 15.36 ± 0.10 [56], 14 we obtain the complete SM prediction

$$\psi ^ { \prime } _ { \epsilon } \quad a _ { \mu } ^ { \text{SM} } = 1 1 6 5 9 \, 1 8 3. 1 \pm 4. 0 \pm 2. 6 \pm 0. 1 \, \left ( 4. 8 _ { \text{tot} } \right ), \quad ( 7 )$$

13 We have verified that the integration of the finely binned R distribution shown in Fig. 8, together with its covariance matrix, accurately reproduces the a had,LO µ and ∆α had ( m 2 Z ) results obtained by summing the exclusive modes below 1.8 GeV in Table 2.

14 When adjusting [57] the new full 2-loop calculation in Ref. [58] to physical quark masses, it reproduces the value obtained in [56].

Table 2. Compilation of the contributions to a had,LO µ and ∆α had ( m 2 Z ) as obtained from HVPTools, and the phenomenological fit for the π + π -contribution below 0.6 GeV. Where three (or more) uncertainties are given, the first is statistical, the second channel-specific systematic, and the third common systematic, which is correlated with at least one other channel. For the contributions computed from QCD, only total uncertainties are given, which include effects from the α S uncertainty, the truncation of the perturbative series at four loops, the FOPT vs. CIPT ambiguity, and quark mass uncertainties. Except for the latter uncertainty, all other uncertainties are taken to be fully correlated among the various energy regions where QCD is used. The additional uncertainty dubbed 'dual' estimates possible quark-hadron duality violating effects in the QCD estimate between 1.8 and 2.0 GeV. The uncertainties in the Breit-Wigner integrals of the narrow resonances J/ψ and ψ (2 S ) are dominated by the the respective electronic width measurements [6]. The uncertainties in the sums (last line) are obtained by quadratically adding all statistical and channel-specific systematic uncertainties, and by linearly adding correlated inter-channel systematic uncertainties.

| Channel                           | a had,LO µ [10 - 10 ]             | ∆α had ( m 2 Z ) [10 - 4 ]   |
|-----------------------------------|-----------------------------------|------------------------------|
| π 0 γ                             | 4.41 ± 0.06 ± 0.04 ± 0.07         | 0.35 ± 0.00 ± 0.00 ± 0.01    |
| ηγ                                | 0.65 ± 0.02 ± 0.01 ± 0.01         | 0.08 ± 0.00 ± 0.00 ± 0.00    |
| π + π -                           | 507.85 ± 0.83 ± 3.23 ± 0.55       | 34.50 ± 0.06 ± 0.20 ± 0.04   |
| π + π - π 0                       | 46.21 ± 0.40 ± 1.10 ± 0.86        | 4.60 ± 0.04 ± 0.11 ± 0.08    |
| 2 π + 2 π -                       | 13.68 ± 0.03 ± 0.27 ± 0.14        | 3.58 ± 0.01 ± 0.07 ± 0.03    |
| π + π - 2 π 0                     | 18.03 ± 0.06 ± 0.48 ± 0.26        | 4.45 ± 0.02 ± 0.12 ± 0.07    |
| 2 π + 2 π - π 0 ( η excl.)        | 0.69 ± 0.04 ± 0.06 ± 0.03         | 0.21 ± 0.01 ± 0.02 ± 0.01    |
| π + π - 3 π 0 ( η excl.)          | 0.49 ± 0.03 ± 0.09 ± 0.00         | 0.15 ± 0.01 ± 0.03 ± 0.00    |
| 3 π + 3 π -                       | 0.11 ± 0.00 ± 0.01 ± 0.00         | 0.04 ± 0.00 ± 0.00 ± 0.00    |
| 2 π + 2 π - 2 π 0 ( η excl.)      | 0.71 ± 0.06 ± 0.07 ± 0.14         | 0.25 ± 0.02 ± 0.02 ± 0.05    |
| π + π - 4 π 0 ( η excl., isospin) | 0.08 ± 0.01 ± 0.08 ± 0.00         | 0.03 ± 0.00 ± 0.03 ± 0.00    |
| ηπ + π -                          | 1.19 ± 0.02 ± 0.04 ± 0.02         | 0.35 ± 0.01 ± 0.01 ± 0.01    |
| ηω                                | 0.35 ± 0.01 ± 0.02 ± 0.01         | 0.11 ± 0.00 ± 0.01 ± 0.00    |
| ηπ + π - π 0 (non- ω , φ )        | 0.34 ± 0.03 ± 0.03 ± 0.04         | 0.12 ± 0.01 ± 0.01 ± 0.01    |
| η 2 π + 2 π -                     | 0.02 ± 0.01 ± 0.00 ± 0.00         | 0.01 ± 0.00 ± 0.00 ± 0.00    |
| ωηπ 0                             | 0.06 ± 0.01 ± 0.01 ± 0.00         | 0.02 ± 0.00 ± 0.00 ± 0.00    |
| ωπ 0 ( ω → π 0 γ )                | 0.94 ± 0.01 ± 0.03 ± 0.00         | 0.20 ± 0.00 ± 0.01 ± 0.00    |
| ω 2 π ( ω → π 0 γ )               | 0.07 ± 0.00 ± 0.00 ± 0.00         | 0.02 ± 0.00 ± 0.00 ± 0.00    |
| ω (non-3 π , πγ , ηγ )            | 0.04 ± 0.00 ± 0.00 ± 0.00         | 0.00 ± 0.00 ± 0.00 ± 0.00    |
| K + K -                           | 23.08 ± 0.20 ± 0.33 ± 0.21        | 3.35 ± 0.03 ± 0.05 ± 0.03    |
| K S K L                           | 12.82 ± 0.06 ± 0.18 ± 0.15        | 1.74 ± 0.01 ± 0.03 ± 0.02    |
| φ (non- KK , 3 π , πγ , ηγ )      | 0.05 ± 0.00 ± 0.00 ± 0.00         | 0.01 ± 0.00 ± 0.00 ± 0.00    |
| KKπ                               | 2.45 ± 0.05 ± 0.10 ± 0.06         | 0.78 ± 0.02 ± 0.03 ± 0.02    |
| KK 2 π                            | 0.85 ± 0.02 ± 0.05 ± 0.01         | 0.30 ± 0.01 ± 0.02 ± 0.00    |
| KKω                               | 0.00 ± 0.00 ± 0.00 ± 0.00         | 0.00 ± 0.00 ± 0.00 ± 0.00    |
| ηφ                                | 0.33 ± 0.01 ± 0.01 ± 0.00         | 0.11 ± 0.00 ± 0.00 ± 0.00    |
| ηKK (non- φ )                     | 0.01 ± 0.01 ± 0.01 ± 0.00         | 0.00 ± 0.00 ± 0.01 ± 0.00    |
| ω 3 π ( ω → π 0 γ )               | 0.06 ± 0.01 ± 0.01 ± 0.01         | 0.02 ± 0.00 ± 0.00 ± 0.00    |
| 7 π (3 π + 3 π - π 0 +estimate)   | 0.02 ± 0.00 ± 0.01 ± 0.00         | 0.01 ± 0.00 ± 0.00 ± 0.00    |
| J/ψ (BW integral)                 | 6.20 ± 0.11                       | 7.00 ± 0.13                  |
| ψ (2 S ) (BW integral)            | 1.56 ± 0.05                       | 2.48 ± 0.08                  |
| R data [3.7 - 5.0] GeV            | 7.29 ± 0.05 ± 0.30 ± 0.00         | 15.79 ± 0.12 ± 0.66 ± 0.00   |
| R QCD [1.8 - 3.7 GeV] uds         | 33.45 ± 0.28 ± 0.65 dual          | 24.27 ± 0.18 ± 0.28 dual     |
| R QCD [5.0 - 9.3 GeV] udsc        | 6.86 ± 0.04                       | 34.89 ± 0.18                 |
| R QCD [9.3 - 12.0 GeV] udscb      | 1.20 ± 0.01                       | 15.53 ± 0.04                 |
| R QCD [12.0 - 40.0 GeV] udscb     | 1.64 ± 0.00                       | 77.94 ± 0.13                 |
| R QCD [ > 40.0 GeV] udscb         | 0.16 ± 0.00                       | 42.70 ± 0.05                 |
| R QCD [ > 40.0 GeV] t             | 0.00 ± 0.00                       | - 0.72 ± 0.01                |
| Sum                               | 1.0 ± 3.5 ± 1.6 ± 0.1 ψ ± 0.7 QCD | ± 0.23 ± 0.15 ψ ± 0.55 QCD   |

Fig. 8. The total hadronic e + e -annihilation rate R as a function of centre-of-mass energy. Inclusive measurements from BES [49] and KEDR [50, 51] are shown as data points, while the sum of exclusive channels from this analysis is given by the narrow blue bands. Also shown for the purpose of illustration is the prediction from massless perturbative QCD (solid red line).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-e90eb424-7a3c-4c1a-b941-3d02f25030cb.jpg)

where the uncertainties account for lowest and higher order hadronic, and other contributions, respectively. The result (7) deviates from the experimental value, a exp µ = 11659209.1 ± 5.4 ± 3.3 [55, 59], by 26.0 ± 7.9 (3.3 σ ).

uncertainties are due to the top mass and of theoretical origin), but dominates the prediction of sin 2 θ /lscript eff , which, however, is about twice more accurate than the combination of all present measurements [47].

A compilation of recent SM predictions for a µ compared with the experimental result is given in Fig. 9.

## Running electromagnetic coupling at m 2 Z

The sum of all quark-flavour terms from Table 2 gives for the hadronic contribution to the running of α m ( 2 Z )

$$\Delta \alpha _ { \text{had} } ( m _ { Z } ^ { 2 } ) = ( 2 7 5. 3 \pm 1. 0 ) \cdot 1 0 ^ { - 4 } \,, \quad ( 8 ) \quad \text{zero}.$$

the uncertainty of which is dominated by data systematic effects (0.7 · 10 -4 ) and the uncertainty in the QCD prediction (0.6 · 10 -4 ). The use of the same inputs with different integration kernels in the calculations induces a correlation of +44% between the a had,LO µ and ∆α had ( m 2 Z ) uncertainties. The result without the new BABAR/KLOE systematic uncertainty is 275.2 ± 0.9.

Adding to (8) the four-loop leptonic contribution, ∆α lep ( m 2 Z ) = (314.979 ± 0.002) · 10 -4 [61], one finds

$$\alpha ^ { - 1 } ( m _ { Z } ^ { 2 } ) = 1 2 8. 9 4 7 \pm 0. 0 1 3 \,. \quad \ ( 9 ) \quad \underset { \text{me} } { \text{com} }$$

The current uncertainty on α m ( 2 Z ) is sub-dominant in the SM prediction of the W -boson mass (the dominant

## 5 Conclusions and perspectives

Using newest available e + e -→ hadrons cross-section data we have reevaluated the lowest-order hadronic vacuum polarisation contribution to the Standard Model prediction of the anomalous magnetic moment of the muon, and the hadronic contribution to the running electromagnetic coupling strength at the Z -boson mass. For the former quantity we find a had,LO µ = (694.0 ± 4.0) · 10 -10 . In spite of new data and the use of a more precise fit to evaluate the threshold region up to 0.6 GeV, the uncertainty on this contribution has increased to 0.6% since our last evaluation [1], due to the addition of a new systematic uncertainty to account for a global discrepancy between π + π -data from BABAR and KLOE. Resolving this discrepancy would allow to reduce the a had,LO µ uncertainty by 20%. 15

The discrepancy between measurement and complete Standard Model prediction remains at a nonconclusive 3.3 σ level. The new Fermilab g -2 experiment currently in operation [62] aims at up to four

15 The contribution of the π + π -channel to the total a had,LO µ uncertainty-squared is 71%.

]

Fig. 9. Compilation of recent data-driven results for a SM µ (in units of 10 -10 ), subtracted by the central value of the experimental average [55, 59]. The blue vertical band indicates the experimental uncertainty, with the darker inlet representing the experimental systematic uncertainty. The representative SM predictions are taken from KNT 2019 [41], J 2018 [60], and this work (DHMZ 2019).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00921v3-figure-1f6179b5-cdb5-4d39-9b29-af33d1aae461.jpg)

times better ultimate precision and has the potential to clarify the situation.

To match the precision of the new experiment further progress is needed to reduce the uncertainty on a had,LO µ from dispersion relations. New analyses of the dominant π + π -channel are underway at the BABAR, CMD-3 and SND experiments for which a systematic uncertainty below 0.5% may be reachable. It is also important to improve the precision of the π + π -π 0 and K K + -channels. The new Belle-2 experiment at the KEK Super-B factory will also contribute to measuring hadronic cross sections via the ISR method once the detector performance is fully understood and sufficient statistics has been accumulated.

Independently of the data-driven approach, lattice QCD calculations of a had,LO µ are also progressing albeit not yet reaching competitive precision [63].

The determination of a had,LO µ is closing in on the estimated uncertainty of the hadronic light-by-light scattering contribution a had,LBL µ of 2.6 · 10 -10 , which appears irreducible at present. Here only phenomenological models have been used so far and lattice QCD calculations could have a strong impact [64], as well as a new promising dispersive approach [65].

## References

- 1. M. Davier, A. Hoecker, B. Malaescu and Z. Zhang, Eur. Phys. J. C 77, 827 (2017) [arXiv:1706.09436].
- 2. M. Davier, A. Hoecker, B. Malaescu, C.Z. Yuan, and Z. Zhang, Eur. Phys. J. C 66, 1 (2010) [arXiv:0908.4300].
- 3. ATLAS Collaboration, Eur. Phys. J. C 75, 17 (2015) [arXiv:1406.0076].
- 4. ATLAS Collaboration, JHEP 09, 020 (2017) [arXiv:1706.03192].
- 5. ATLAS Collaboration, JHEP 05, 195 (2018) [arXiv:1711.02692].
- 6. Particle Data Group, Phys. Rev. D 98, 030001 (2018) and 2019 update. http://pdglive.lbl. gov .
- 7. CMD-3 Collaboration, Phys. Lett. B 773, 150 (2017) [arXiv:1706.06267].
- 8. BABAR Collaboration, Phys. Rev. D 96, 092009 (2017) [arXiv:1709.01171].
- 9. M.N. Achasov et al. , Phys. Rev. D 97, 032011 (2018) [arXiv:1711.07143].
- 10. T. Xiao, S. Dobbs, A. Tomaradze, K.K. Seth and G. Bonvicini, Phys. Rev. D 97, 032012 (2018) [arXiv:1712.04530].
- 11. BABAR Collaboration, Phys. Rev. D 97, 052007 (2018) [arXiv:1801.02960].
- 12. M.N. Achasov et al. , Phys. Rev. D 98, 112001 (2018) [arXiv:1809.07631].
- 13. BABAR Collaboration, Phys. Rev. D 98, 112015 (2018) [arXiv:1810.11962].
- 14. CMD-3 Collaboration, Phys. Lett. B 792, 419 (2019) [arXiv:1902.06449].
- 15. M.N. Achasov et al. , Phys. Rev. D 99, 112004 (2019) [arXiv:1903.09307].
- 16. V.L. Ivanov et al. , Phys. Lett. B 798 (2019) 134946 [arXiv:1906.08006].
- 17. M. Davier, A. Hoecker, B. Malaescu, and Z. Zhang, Eur. Phys. J. C 71, 1515 (2011) [arXiv:1010.4180].
- 18. M. Davier, S. Eidelman, A. Hoecker, and Z. Zhang, Eur. Phys.J. C 27, 497 (2003) [hep-ph/0208177]; Eur. Phys. J. C 31, 503 (2003) [hep-ph/0308213].
- 19. P.A. Baikov, K.G. Chetyrkin, and J.H. K¨hn, Phys. u Rev. Lett. 101, 012002 (2008) [arXiv:0801.1821].
- 20. BABAR Collaboration, Phys. Rev. D 8, 032013 (2012) [arXiv:1205.2228].
- 21. BABAR Collaboration, Phys. Rev. Lett. 103, 231801 (2009) [arXiv:0908.3589].
- 22. KLOE Collaboration, Phys. Lett. B 670, 285 (2009) [arXiv:0809.3950].
- 23. KLOE Collaboration, Phys. Lett. B 700, 102 (2011) [arXiv:1006.5313].
- 24. KLOE Collaboration, Phys. Lett. B 720, 336 (2013) [arXiv:1212.4524].
- 25. BESIII Collaboration, Phys. Lett. B 753, 629 (2015) [arXiv:1507.08188].
- 26. KLOE-2 Collaboration, JHEP 03, 173 (2018) [arXiv:1711.03085].
- 27. CMD-2 Collaboration, Phys. Lett. B 578, 285 (2004) [hep-ex/0308008].
- 28. CMD-2 Collaboration, JETP Lett. 82, 743 (2005) [hep-ex/0603021]; CMD-2 Collaboration, JETP Lett. 84, 413 (2006) [hep-ex/0610016]; CMD-2 Collaboration, Phys. Lett. B 648, 28 (2007) [hepex/0610021].

- 29. SND Collaboration, JETP Lett. 103, 380 (2006) [hep-ex/0605013].
- 30. C. Hanhart, S. Holz, B. Kubis, A. Kup´´, A. Wirzba sc and C.W. Xiao, Eur. Phys. J. C 77, 98 (2017) Erratum: [Eur. Phys. J. C 78, 450 (2018)] [arXiv:1611.09359].
- 31. J.F. De Troconiz and F.J. Yndurain, Phys. Rev. D 65, 093001 (2002) [hep-ph/0106025].
- 32. J.F. de Troconiz and F.J. Yndurain, Phys. Rev. D 71, 073008 (2005) [hep-ph/0402285].
- 33. R. Garcia-Martin, R. Kaminski, J.R. Pelaez, J. Ruiz de Elvira and F.J. Yndurain, Phys. Rev. D 83, 074004 (2011) [arXiv:1102.2183].
- 34. M. Hoferichter, B. L. Hoid and B. Kubis, JHEP 08, 137 (2019) [arXiv:1907.01556].
- 35. G. Colangelo, M. Hoferichter and P. Stoffer, JHEP 02, 006 (2019) [arXiv:1810.00007].
- 36. B. Ananthanarayan, I. Caprini and D. Das, Phys. Rev. D 98, 114015 (2018) [arXiv:1810.09265].
- 37. B. Ananthanarayan, I. Caprini and D. Das, Rom. J. Phys. 64 (2019) 401, [arXiv:1907.01767].
- 38. SNDCollaboration, Phys. Rev. D 94, 112006 (2016) [arXiv:1608.08757].
- 39. E.A. Kozyrev et al. , Phys. Lett. B 779, 64 (2018) [arXiv:1710.02989].
- 40. K. Hagiwara, A.D. Martin, D. Nomura, and T. Teubner, Phys. Rev. D 69, 093003 (2004), [arXiv:hep-ph/0312250].
- 41. A. Keshavarzi, D. Nomura,and T. Teubner, Phys. Rev. D 101, 014029 (2020), [arXiv:1911.00367].
- 42. BABAR Collaboration, Phys. Rev. D 88, 032013 (2013) [arXiv:1306.3600].
- 43. CMD-2 Collaboration, Phys. Lett. B 669, 217 (2008) [arXiv:0804.0178].
- 44. BABAR Collaboration, Phys. Rev. D 76, 092005 (2007); Erratum: Phys. Rev. D 77, 119902 (2008) [arXiv:0708.2461].
- 45. K.G. Chetyrkin, J.H. K¨ uhn and M. Steinhauser, Nucl. Phys. B 482, 213 (1996) [hep-ph/9606230].
- 46. M. Davier and A. Hoecker, Phys. Lett. B 419, 419 (1998) [hep-ph/9801361].
- 47. Gfitter Group, Eur. Phys. J. C 74, 3046 (2014) [arXiv:1407.3792].
- 48. F. Le Diberder and A. Pich, Phys. Lett. B 286, 147 (1992).
- 49. BES Collaboration, Phys. Rev. Lett. 84, 594 (2000) [hep-ex/9908046]; Phys. Rev. Lett. 88, 101802 (2002) [hep-ex/0102003]; BES Collaboration, Phys. Lett. B 641, 145 (2006) [hep-ex/0605105]; Phys. Lett. B 677, 239 (2009) [arXiv:0903.0900].
- 50. KEDRCollaboration, Phys. Lett. B 753, 533 (2016) [arXiv:1510.02667].
- 51. KEDRCollaboration, Phys. Lett. B 770, 174 (2017) [arXiv:1610.02827].
- 52. A. Kurz T. Liu, P. Marquard, and M. Steinhauser, Phys. Lett. B 734, 144 (2014) [arXiv:1403.6400].
- 53. J. Prades, E. de Rafael, and A. Vainshtein, Adv. Ser. Direct. High Energy Phys. 20, 303 (2009) [arXiv:0901.0306].
- 54. T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. Lett. 109, 111808 (2012) [arXiv:1205.5370].
- 55. A. Hoecker and W. Marciano, 'The Muon Anomalous Magnetic Moment', in: Review of Particle Physics, Particle Data Group, Phys. Rev. D 98, 030001 (2018).
- 56. C. Gnendiger, D. St¨kinger, and H. St¨kinger-Kim, o o Phys. Rev. D 88, 053005 (2013) [arXiv:1306.5546]; M. Knecht, S. Peris, M. Perrottet, and E. de Rafael, JHEP 11, 003 (2002) [hep-ph/0205102]; A. Czarnecki, W.J. Marciano, and A. Vainshtein, Phys. Rev. D 67, 073006 (2003), Erratum-ibid. D 73, 119901 (2006) [hep-ph/0212229]; R. Jackiw and S. Weinberg, Phys. Rev. D 5, 2396 (1972).
- 57. W.J. Marciano, private communication, Aug. 2019.
- 58. T. Ishikawa, N. Nakazawa, Y. Yasui, Phys. Rev. D 99, 073004 (2019) [arXiv:1810.13445].
- 59. Muon g -2 Collaboration, Phys. Rev. D 73, 072003 (2006) [hep-ex/0602035].
- 60. F. Jegerlehner, EPJ Web of Conferences 166, 00022 (2018) [arXiv:1705.00263].
- 61. C. Sturm, Nucl. Phys. B 874, 698 (2013) [arXiv:1305.0581].
- 62. Muon g -2 Collaboration, FERMILAB-DESIGN2014-02, arXiv:1501.06858 (2015).
- 63. BMW Collaboration, Phys. Rev. Lett. 121, 022002 (2018) [arXiv:1711.04980], arXiv:2002.12347 (2020); HPQCD Collaboration, arXiv:1902.04223 (2019); RBC and UKQCD Collaborations, Phys. Rev. Lett. 121, 022003 (2018) [arXiv:1801.07224]; ETM Collaboration, JHEP 10, 157 (2017) [arXiv:1707.03019], Phys. Rev. D 98, 114504 (2018) [arXiv:1808.00887], Phys. Rev. D 99, 114502 (2019) [arXiv:1901.10462]; A. G´rardin e et al. , Phys. Rev. D 100, 014510 (2019) [arXiv:1904.03120].
- 64. T. Blum et al. , arXiv:1907.00864 (2019).
- 65. M. Hoferichter et al. , JHEP 10, 141 (2018) [arXiv:1808.04823].
M. Davier A. Hoecker B. Malaescu Z. Zhang
A new evaluation of the hadronic vacuum polarisation contributions to the muon anomalous magnetic moment and to $\mathbf{\boldsymbolα(m_Z^2)}$
We reevaluate the hadronic vacuum polarisation contributions to the muon
magnetic anomaly and to the running of the electromagnetic coupling constant at
the $Z$-boson mass. We include newest $e^+e^- \to$ hadrons cross-section data
together with a phenomenological fit of the threshold region in the evaluation
of the dispersion integrals. The precision in the individual datasets cannot be
fully exploited due to discrepancies that lead to additional systematic
uncertainty in particular between BABAR and KLOE data in the dominant
$\pi^+\pi^-$ channel. For the muon $(g-2)/2$, we find for the lowest-order
hadronic contribution $(694.0 \pm 4.0)\cdot10^{-10}$. The full Standard Model
prediction differs by $3.3\sigma$ from the experimental value. The five-quark
hadronic contribution to $\alpha(m_Z^2)$ is evaluated to be
$(276.0\pm1.0)\cdot10^{-4}$.