The Standard Model (SM) predictions of the anomalous magnetic moment of the muon,
, with the muon gyromagnetic factor,
and of the running electromagnetic coupling constant, , a crucial ingredient
of electroweak theory, are
limited in precision by hadronic vacuum polarisation (HVP) contributions.
The dominant hadronic terms can be calculated with a combination of experimental
cross-section data, involving annihilation to hadrons, and perturbative QCD. They are
used to evaluate energy-squared dispersion integrals ranging from the
threshold to infinity. The kernels occurring in these integrals
emphasise low photon virtualities, owing to the descent of the cross section,
and, in case of , to an additional suppression. In the latter case, about 73%
of the lowest order hadronic contribution and 59% of the total uncertainty-squared
are given by the final state,
In this work, we reevaluate the lowest-order hadronic contribution, , to the muon magnetic anomaly, and the hadronic contribution, , to the running at the -boson mass using newest cross-section data. The BABAR Collaboration has essentially completed a programme of precise measurements of exclusive hadronic cross sections for all the dominant channels from threshold to an energy of 3–5GeV using the initial-state radiation (ISR) method. Also new results are being produced at the VEPP-2000 facility in Novosibirsk, Russia in the 1–2GeV energy range. The new data complement the available information on exclusive channels allowing to alleviate the need for estimating missing channels with the use of isospin symmetry.
We reevaluate all the experimental contributions using the software package HVPTools [], and add to these narrow resonance contributions evaluated analytically and continuum contributions obtained from perturbative QCD.
2 Input data
Exclusive bare hadronic cross-section measurements are integrated up to 1.8GeV over the relevant dispersion kernels. In the present work 39 channels are included, as compared to only 22 in our latest work from 2011 []. Thanks to the new measurements only very few final states remain to be estimated using isospin symmetry. In the energy range 1.8–3.7GeV and above 5GeV four-loop perturbative QCD is used []. The contributions from the open charm pair production region between 3.7 and 5GeV are again computed using experimental data. For the narrow resonances and Breit-Wigner line shapes are integrated using their currently best known parameters.
The integration of data points belonging to different experiments with their own data densities requires a careful treatment especially with respect to correlated systematic uncertainties within the same experiment and between different experiments. Quadratic interpolation of adjacent data points is performed for each experiment and a local combination between the interpolations is computed in bins of 1MeV. Full covariance matrices are constructed between experiments and channels. Uncertainties are propagated using pseudo-experiment generation and closure tests with known distributions are performed to validate both the combination and integration. Where results from different experiments are locally inconsistent the combined uncertainty is rescaled according to the local value following the well-known PDG approach. At present, for the dominant channel such inconsistencies are limiting the precision of the combination. In most exclusive channels the largest weight in the combination is provided by BABAR measurements.
The following channel-wise discussion focuses on the HVP contribution to as it stronger relies on the low-energy experimental data. We mainly explore the impact of the data released since our last publication [], which provides references to all the older datasets used in the combination.
2.1 The dominant \sansmath channel
The channel dominates both the HVP contribution to and its uncertainty. Recent experiments are generally limited by systematic uncertainties. The main contributors are BABAR [, ] (relative systematic uncertainty of 0.5% per measurement), KLOE-2008 [] (0.8%), KLOE-2010 [] (1.4%), CMD2 [, ] (0.8%), and SND [] (1.5%), For this update we newly included KLOE-2012 [] (0.8%) and the recent BESIII-2015 [] (0.9%). Only BABAR covers the full mass range with high precision.
The three KLOE measurements exhibit statistical correlations due to the common two-pion events used in the 2008 and 2012 results (the 2012 analysis uses the pion over muon pair cross-section ratio), and systematic correlations from common uncertainty sources.
Figure 1 shows the available cross-section measurements in various panels zooming into different centre-of-mass energies (). The green band indicates the HVPTools combination within its uncertainty. The deviations between the combination and the most precise individual measurements are plotted in Fig. 2. Figure 3 (left) shows the local combination weight versus per experiment. The BABAR and KLOE measurements dominate over the entire energy range. Owing to the sharp radiator function, the event yield increases for KLOE towards the mass, hence outperforming BABAR above . The group of experiments labelled “Other exp” in Fig. 3 corresponds to older data with incomplete radiative corrections. Their weights are small throughout the entire energy domain. The computation of the dispersion integral over the full spectrum requires to extend the available data to the region between threshold and , for which we use a fit as described in Ref. [].
A tension between the BABAR and KLOE measurements is observed at and above the peak region (cf. Fig. 2), while the other measurements are consistent with both. We stress the importance of locally assessing the compatibility of the (correlated) cross-section measurements, rather than comparing integrated values where discrepancies could cancel or be diluted. The local uncertainty rescaling we apply (cf. right-hand plot of Fig. 3) increases the combined uncertainty by 15% in the channel.
In spite of this problem, progress in the evaluation of the contribution to has been steady during the last decade. While the central value stayed within
quoted uncertainties, the uncertainty dropped from
Our () estimate using data from ALEPH, OPAL, CLEO, and Belle, [], where the first uncertainty is experimental and the second due to isospin-breaking corrections, is larger than the current -based value. The difference can be reduced by applying off-resonance – mixing corrections [] that come with additional uncertainties. Because of the progress in the data, the input is now less precise and less reliable due to additional theoretical uncertainties. While the versus comparison is interesting in its own right, we do no longer consider the data for the HVP evaluation.
2.2 The channel
Following the treatment described in Ref. [] the contributions from and decaying to three pions are directly evaluated from the cross-section measurements (cf. Fig. 4). Other resonant decays are included in the corresponding , , and spectra, while small remaining non-resonant decay modes are considered separately.
2.3 The four-pion channels
Recent results using the full BABAR data on are now available []. As with other BABAR measurements using the ISR method with the ISR photon measured at large angle, the acceptance for the recoiling hadronic system is large so that the resonance substructure, dominated by the , , and final states, can be fully identified and accurately modelled with a Monte Carlo generator. The systematic uncertainty is 3.1% below 2.7GeV, a considerable improvement over the value of about 10% of preliminary results available so far. Data from some older experiments, both imprecise and inconsistent, are now discarded. As seen in the left hand plot of Fig. 6 the BABAR results lead to a substantial precision improvement in this channel.
The HVP contribution to from threshold to 1.8GeV is , where the total uncertainty of 0.55 is reduced by a factor of 2.3 compared to our 2011 result []. We note that the -based result (the second uncertainty accounts for isospin-symmetry breaking corrections), obtained from a combination of and spectral functions measured by ALEPH [], is larger than the value and twice less precise.
New cross-section data (cf. right hand plot in Fig. 6) were published by BABAR in 2012 [] using the full available data sample and with a reduced systematic uncertainty (2.4%) compared to previous partial results. New measurements from CMD3 between 0.920 and 1.060GeV are also available []. The resulting combined HVP contribution is , with a total uncertainty of 0.31 reduced by a factor of 1.7 compared to our 2011 result [].
For comparison, the ALEPH -based prediction of , [], is consistent, but more than twice less precise than the -based one. The -based evaluation of the sum of the two four-pion channels, , benefits from an anticorrelation due to the contribution in both channels. It is consistent with the -based value of within . The -based cross-section predictions are compared to the data in Fig. 6.
2.4 The \sansmath channels
New cross-section measurements are available for the channel. The BABAR experiment detects both and from threshold up to 2.2GeV [], while CMD3 counts in the resonance region []. Consistency is observed between the two experiments as well as with older measurements from CMD2 and SND. The measured cross sections are shown in the right panel of Fig. 6.
The contribution to up to 1.8GeV amounts to with a total uncertainty of 0.24, which is reduced by a factor of 1.6 over that of our 2011 estimate [].
Recent measurements from SND [] at VEPP-2000 for the channel agree with BABAR [], while both show a discrepancy with former SND data, obtained at VEPP-2M below 1.4GeV, that exceeds the quoted systematic uncertainty. The BABAR and new SND data are displayed in the left hand panel of Fig. 6.
Some concern arises with regard to the cross-section measurements. The BABAR result has a systematic uncertainty of 0.7%, but it is larger by 5.1% (9.6%) than CMD2 (SND) which has a systematic uncertainty of 2.2% (7.1%). Including the BABAR data the contribution to increases from 21.63 to 22.67 with an uncertainty of 0.43. A recent preliminary result from CMD3 [] shows a very large (11%) excess of the cross section over CMD2 and 5% over BABAR. The origin of this large discrepancy is not understood at present []. It raises doubts on the ability to perform a precision measurement of this channel with the standard energy-scan method because the detection efficiency of the low-momentum from decay is difficult to model. Owing to the boost of the final state, the ISR method is expected to be more reliable for the charged kaon detection.
2.5 The \sansmathpions channels
In previous hadronic vacuum polarisation analyses the available exclusive pions data were incomplete. Missing channels were constrained based on assumptions about the process dynamics and isospin symmetry [, ] leading to considerable uncertainty. This procedure became unnecessary since the BABAR experiment produced cross-section results for the three channels contributing to the final state and six channels contributing to . A key ingredient of the BABAR analyses is the detection of neutral kaons, both , through the decay, and interacting in the calorimeters.
Together with previous measurements of and , data for [] complete the picture for the channel (cf. top row panels in Fig. 7). Because that final state is dominated by production below 1.8GeV (with a small contribution from ), it is expected that isospin symmetry provides a good approximation. Indeed, the contribution from the sum of the measured channels, , is in agreement with and has similar precision as the result obtained from the data only together with isospin symmetry.
Of the six channels contributing to only two, and , had been measured by BABAR in 2011. Constraints from isospin symmetry were used to estimate the missing channels [], but because of the complex dynamics involving , and intermediate states, these estimates were plagued by substantial uncertainties. Among the remaining channels, [], [], [], and [] have been measured by BABAR. In addition, the previously measured channels and have been updated with the full data sample []. New measurements of became available from CMD3 [] and are in agreement with the BABAR data. Except for , which can be safely estimated using CP symmetry, all cross sections have now been measured.
2.6 Other channels
Data on many processes with smaller cross sections, mainly from VEPP-2000, have become available and are included in the HVPTools database. This is the case for from CMD3 [], [], [], [], and [] from SND, from SND [] and CMD3 [], nonresonant from CMD3 [], from BABAR [], which extend and improve older measurements in these channels. Except for the cross section above 1.6GeV, results using the ISR technique at BABAR and the scan method at VEPP-2000 are in agreement notwithstanding their different systematic uncertainties. Above 1.8GeV the production of measured by BABAR [] and CMD3 [], by ADONE [] and SND [], and by SND [] are included.
Figure 8 shows the available measurements and their combination of the charm resonance region above the opening of the channel. The individual datasets agree within uncertainties. While Crystal Ball [] and BES [] published bare inclusive cross-section results, PLUTO applied only radiative corrections [] following the formalism of Ref. [], which does not include HVP. Following similar previous cases [], we have applied this correction and assigned a 50% systematic uncertainty to it.
2.7 Estimated missing channels
Even with the large number of exclusive cross-section measurements available below 2GeV, covering particle multiplicities up to six hadrons including and mesons, a few channels with more than two neutral pions are still unmeasured and their contributions must be estimated using isospin symmetry. The treatment of the channels , , and follows our previous approach detailed in Ref. [].
Whereas the cross-section data were already included, the previously neglected smaller contribution from where the does not originate from resonant decay is now taken into account following a BABAR measurement []. Its unmeasured counterpart can be crudely estimated to equal the corresponding contribution with a 100% systematic uncertainty. This estimate is consistent with the upper limit that can be deduced using BABAR’s measurement [].
Altogether the contribution from missing channels to up to 1.8GeV is estimated to be , corresponding to a fraction of only % of the full HVP contribution in this range. The corresponding fraction in our 2011 analysis was much larger, %, illustrating the experimental progress made.
3 Compilation and results
|( excl., isospin)|
|( excl., isospin)|
A compilation of the various contributions to and to , as well as the total results are given in Table 1. 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 and resonances in
Table 1 are obtained by numerically integrating the corresponding
Sufficiently far from the quark thresholds we use four-loop [] perturbative QCD, including quark mass corrections [], to compute the inclusive hadronic cross section. Nonperturbative contributions at were determined from data [] and found to be small. The uncertainties of the contributions given in Table 1 are obtained from the quadratic sum of the uncertainty in (we use from the fit to precision data []), the truncation of the perturbative series (we use the full four-loop contribution as systematic uncertainty), the difference between fixed-order perturbation theory and, so-called, contour-improved perturbation theory [], as well as quark mass uncertainties (we use the values and uncertainties from Ref. []). The former three uncertainties are taken to be fully correlated between the various energy regions (see Table 1), whereas the (smaller) quark-mass uncertainties are taken to be uncorrelated.
The KEDR Collaboration has recently published results from an inclusive scan from to 3.05GeV [], complementing their previous measurements obtained between 3.12 and 3.72GeV []. These data are the most precise and complete in this energy range with a typical systematic uncertainty of 3% for a total of 20 measurements. They constitute a very valuable input to test the validity of the perturbative QCD estimate (cf. Fig. 9). Integrating the dispersion integrals between 2.0 and 3.7GeV gives for the combined data ( in the usual units) and (), compared to the QCD predictions and , respectively. Agreement within is found.
To examine the transition region between the sum of exclusive measurements and QCD we have computed and in the narrow energy interval –. For the former quantity we find and for data and QCD, respectively. The full difference of ( in the case of ) is assigned as additional systematic uncertainty, labelled by “dual” subscripts in Table 1. It accounts for possible low-mass quark-hadron duality violation affecting the perturbative QCD approximation that we use for this interval to avoid systematic effects due to possible unmeasured high-multiplicity channels.
Figure 10 shows the total hadronic annihilation rate versus centre-of-mass
energy as obtained from the sum of exclusive data below 2GeV and from inclusive data between
1.8 and 5GeV.
Muon magnetic anomaly
Adding all lowest-order hadronic contributions together gives
which is dominated by experimental systematic uncertainties (cf. Table 1 for a separation of the total uncertainty into its components). The new result is units larger than our previous evaluation [] and 21% more precise owing to the new and improved experimental data.
Adding to (1) the contributions from higher order hadronic loops, (NLO) and (NNLO) [], hadronic light-by-light scattering, [], as well as QED, [] (see also [] and references therein), and electroweak effects, [], we obtain the complete SM prediction
A compilation of recent SM predictions for compared with the experimental result is given in Fig. 11.
Running electromagnetic coupling at
The sum of all quark-flavour terms from Table 1 gives for the hadronic contribution to the running of
the uncertainty of which is dominated by data systematic effects () and the uncertainty in the QCD prediction ().
The current uncertainty on is sub-dominant in the SM prediction of the -boson mass (the dominant uncertainties are due to the top mass and of theoretical origin), but dominates the prediction of , which, however, is about twice more accurate than the combination of all present measurements [].
4 Conclusions and perspectives
Using newest available 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 -boson mass. For the former quantity we find . The uncertainty of 0.5% on this contribution is now reduced to about half the current uncertainty of the measurement, and has improved by more than a factor of two during the last thirteen years. The discrepancy between measurement and complete Standard Model prediction remains at a non-conclusive level. The forthcoming experiments at Fermilab [] and JPARC [], aiming at up to four times better ultimate precision, have the potential to clarify the situation.
To match the precision of these experiments further progress is needed to reduce the uncertainty on from dispersion relations. New analyses of the dominant channel are underway at the BABAR and CMD3 experiments for which a systematic uncertainty of 0.3% may be reachable. In the 1–2GeV range it is important to improve the precision of the and channels. Independently of the data-driven approach, Lattice QCD calculations of are also progressing albeit not yet reaching competitive precision [].
The determination of is closing in on the estimated uncertainty of the hadronic light-by-light scattering contribution of , which appears irreducible at present. Here only phenomenological models have been used so far and Lattice QCD calculations could have a strong impact [].
- Throughout this paper, final state photon radiation is implied for all hadronic final states.
- Correlations due to systematic uncertainties among the three KLOE datasets are accounted for by matching one-by-one the uncertainties impacting the measurements. For example, the uncertainties due to the luminosity, radiator function and vacuum polarisation are taken to be correlated among the first two measurements, while they are found negligible or not present in the 2012 measurement. The final state radiation uncertainty as well as the trigger, tracking, acceptance and background uncertainties for the data are taken to be correlated among the three measurements, while the background uncertainty for the data only impacts the 2012 measurement. We will replace this approximate treatment by a more accurate one once it is provided by the KLOE Collaboration.
- Unless specified, these and all other related values throughout this paper are given in units of .
- The undressing uses the BABAR programme Afkvac correcting for both leptonic and hadronic VP effects. The correction factors amount to and for the and , respectively.
- Using instead the narrow-width approximation, , gives consistent results.
- We have verified that the integration of the finely binned distribution shown in Fig. 10, together with its covariance matrix, accurately reproduces the and results obtained by summing the exclusive modes below 1.8GeV in Table 1.
- M. Davier, A. Hoecker, B. Malaescu, C.Z. Yuan, and Z. Zhang, Eur. Phys. J. C66, 1 (2010) [arXiv:0908.4300].
- M. Davier, A. Hoecker, B. Malaescu, and Z. Zhang, Eur. Phys. J. C 71, 1515 (2011) [arXiv:1010.4180].
- P.A. Baikov, K.G. Chetyrkin, and J.H. Kühn, Phys. Rev. Lett. 101, 012002 (2008) [arXiv:0801.1821].
- BABAR Collaboration, Phys. Rev. Lett. 103, 231801 (2009) [arXiv:0908.3589].
- BABAR Collaboration, Phys. Rev. D 86, 032013 (2012) [arXiv:1205.2228].
- KLOE Collaboration, Phys. Lett. B 670, 285 (2009) [arXiv:0809.3950].
- KLOE Collaboration, Phys. Lett. B 700, 102 (2011) [arXiv:1006.5313].
- KLOE Collaboration, Phys. Lett. B 720, 336 (2013) [arXiv:1212.4524].
- BESIII Collaboration, Phys. Lett. B 753, 629 (2015) [arXiv:1507.08188].
- CMD2 Collaboration, Phys. Lett. B 578, 285 (2004) [hep-ex/0308008].
- CMD2 Collaboration, JETP Lett. 82, 743 (2005) [hep-ex/0603021]; CMD2 Collaboration, JETP Lett. 84, 413 (2006) [hep-ex/0610016]; CMD2 Collaboration, Phys. Lett. B 648, 28 (2007) [hep-ex/0610021].
- SND Collaboration, JETP Lett. 103, 380 (2006) [hep-ex/0605013].
- M. Davier, A. Hoecker, B. Malaescu, C.Z. Yuan, and Z. Zhang, Eur. Phys. J. C 74, 2803 (2014) [arXiv:1312.1501].
- F. Jegerlehner and R. Szafron, Eur. Phys. J. C 71, 1632 (2011) [arXiv:1101.2872].
- K. Griessinger, (for the BABAR Collaboration), proceedings of the Tau Workshop, IHEP-Beijing (Sep. 2016).
- BABAR Collaboration, Phys. Rev. D 85, 112009 (2012) [arXiv:1201.5677].
- CMD3 Collaboration, Phys. Lett. B 768, 345 (2017) [arXiv:1612.04483].
- BABAR Collaboration, Phys. Rev. D 89, 092002 (2014) [arXiv:1403.7593].
- CMD3 Collaboration, Phys. Lett. B 760, 314 (2016) [arXiv:1604.02981].
- SND Collaboration, Phys. Rev. D 94, 112006 (2016) [arXiv:1608.08757].
- BABAR Collaboration, Phys. Rev. D 88, 032013 (2013) [arXiv:1306.3600].
- CMD3 Collaboration, J. Univ. Sci. Tech. China 46, 507 (2016) [arXiv:1603.03230].
- E.P. Solodov, private communication (2016).
- M. Davier, S. Eidelman, A. Hoecker, and Z. Zhang, Eur. Phys.J. C27, 497 (2003) [hep-ph/0208177]; Eur. Phys. J. C 31, 503 (2003) [hep-ph/0308213].
- BABAR Collaboration, Phys. Rev. D 95, 052001 (2017) [arXiv:1701.08297].
- BABAR Collaboration, arXiv:1704.05009 (2017).
- BABAR Collaboration, Phys. Rev. D 86, 012008 (2012) [arXiv:1103.3001].
- CMD3 Collaboration, Phys. Lett. B 756, 153 (2016) [arXiv:1510.00654].
- CMD3 Collaboration, Phys. Lett. B 723, 82 (2013) [arXiv:1302.0053].
- SND Collaboration, Phys. Rev. D 93, 092001 (2016) [arXiv:1601.08061].
- SND Collaboration, Phys. Rev. D 90, 032002 (2014) [arXiv:1312.7078].
- SND Collaboration, Phys. Rev. D 91, 052013 (2015) [arXiv:1412.1971].
- SND Collaboration, Phys. Rev. D 88, 054013 (2013) [arXiv:1303.5198]; Phys. Rev. D 94, 112001 (2016) [arXiv:1610.00235].
- SND Collaboration, Phys. Rev. D 94, 092002 (2016) [arXiv:1607.00371].
- CMD3 Collaboration, arXiv:1706.06267 (2017).
- BABAR Collaboration, Phys. Rev. D 91, 052013 (2015) [arXiv:1412.1971].
- BABAR Collaboration, Phys. Rev. D 87, 092005 (2013) [arXiv:1302.0055].
- CMD3 Collaboration, Phys. Lett. B 759, 634 (2016) [arXiv:1507.08013]
- ADONE Collaboration, Nucl. Phys. B 517, 3 (1998).
- SND Collaboration, Phys. Rev. D 90, 112007 (2014) [arXiv:1410.3188].
- SND Collaboration, Phys. Rev. D 94, 032010 (2016) [arXiv:1606.06481].
- 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].
- KEDR Collaboration, Phys. Lett. B, 753, 533 (2016) [arXiv:1510.02667].
- KEDR Collaboration, Phys. Lett. B 770, 174 (2017) [arXiv:1610.02827].
- Crystal Ball Collaboration, Z. Phys. C 40, 49 (1988); Crystal Ball Collaboration, SLAC-PUB-5160 (1990).
- A. Bäcker, Preprint DESY F33-77/03 (1977).
- G. Bonneau and F. Martin, Nucl. Phys. B 27, 381 (1971).
- M. Davier, S. Eidelman, A. Hoecker, and Z. Zhang, Eur. Phys. J. C 27, 497 (2003) [hep-ph/0208177].
- BABAR Collaboration, Phys. Rev. D 95, 052001 (2017) [arXiv:1701.08297].
- BABAR Collaboration, Phys. Rev. D 73, 052003 (2006) [hep-ex/0602006].
- BESIII Collaboration, Phys. Lett. B 761, 98 (2016) [arXiv:1604.01924] (and references therein).
- Particle Data Group, Chin. Phys. C 40, 100001 (2016), http://pdglive.lbl.gov.
- P. Baikov, K.G. Chetyrkin, and J.H. Kühn, Phys. Rev. Lett. 101, 012002 (2008) [arXiv:0801.1821].
- K.G. Chetyrkin, J.H. Kühn and M. Steinhauser, Nucl. Phys. B 482, 213 (1996) [hep-ph/9606230].
- M. Davier and A. Hoecker, Phys. Lett. B 419, 419 (1998) [hep-ph/9801361].
- Gfitter Group, Eur. Phys. J. C 74, 3046 (2014) [arXiv:1407.3792].
- F. Le Diberder and A. Pich, Phys. Lett. B 286, 147 (1992).
- Particle Data Group, J. Phys. G 37, 075021 (2010).
- Muon – Collaboration, Phys. Rev. D 73, 072003 (2006) [hep-ex/0602035].
- A. Hoecker and W. Marciano, “The Muon Anomalous Magnetic Moment”, in: Review of Particle Physics, Physarticle Data Group, Chin. Phys. C, 40, 100001 (2016).
- K. Hagiwara et al., J. Phys. G 38, 085003 (2011) [arXiv:1105.3149].
- F. Jegerlehner and A. Nyffeler, Phys. Rept. 477, 1 (2009) [arXiv:0902.3360].
- A. Kurz et al., EPJ Web Conf. 118, 01033 (2016) [arXiv:1511.08222].
- J. Prades, E. de Rafael, and A. Vainshtein, Adv. Ser. Direct. High Energy Phys. 20, 303 (2009) [arXiv:0901.0306].
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. Lett. 109, 111808 (2012) [arXiv:1205.5370].
- C. Gnendiger, D. Stökinger, and H. Stökinger-Kim, Phys. Rev. D 88, 053005 (2013) [arXiv:1306.5546]; M. Knecht, S. Peris, M. Perrottet, and E. de Rafael, JHEP 0211, 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).
- M. Steinhauser, Phys. Lett. B 429, 158 (1998) [hep-ph/9803313].
- Muon – Collaboration, FERMILAB-FN-0992-E, FERMILAB-DESIGN-2014-02, arXiv:1501.06858 (2015).
- E34 Muon –/EDM experiment at JPARC, http://g-2.kek.jp/portal/index.html.
- HPQCD Collaboration, arXiv:1601.03071 (2016); RBC and UKQCD Collaborations, Phys. Rev. Lett. 116, 232002 (2016) [arXiv:1512.09054]; T. Blum, Phys. Rev. Lett. 91, 052001 (2003).
- T. Blum et al., Phys. Rev. D 93, 014503 (2016) [arXiv:1510.07100].