Complete Tenth-Order QED Contribution to the Muon
We report the result of our calculation of the complete tenth-order QED terms of the muon . Our result is in units of , which is about 4.5 s.d. larger than the leading-logarithmic estimate . We also improved the precision of the eighth-order QED term of , obtaining in units of . The new QED contribution is , which does not resolve the existing discrepancy between the standard-model prediction and measurement of .
The anomalous magnetic moment of the muon has been studied extensively both experimentally and theoretically since it provides one of the promising paths in exploring possible new physics beyond the standard model. For this purpose it is crucial to know the prediction of the standard model as precisely as possible.
In the standard model, can be divided into electromagnetic, hadronic, and electroweak contributions
At present (hadronic) is the largest source of theoretical uncertainty. The uncertainty comes mostly from the hadronic vacuum-polarization (v.p.) term, being the fine-structure constant. The lattice QCD simulations have attempted to evaluate this contribution Blum (2003); Gockeler et al. (2004); Aubin and Blum (2007); Feng et al. (2011); Boyle et al. (2012); Della Morte et al. (2012). At present, most accurate evaluations must rely on the experimental information. Three types of measurements are available for this purpose: (1) , (2) , (3) . These processes have been investigated intensely by many groups Davier et al. (2011); Jegerlehner and Szafron (2011); Hagiwara et al. (2011).We list here one of them Hagiwara et al. (2011):
which overlaps other values based on the data Davier et al. (2011); Jegerlehner and Szafron (2011) and makes the standard-model prediction closest to the experiment (1). The next-to-leading-order (NLO) hadronic vacuum-polarization contribution is also known Hagiwara et al. (2011):
The hadronic light-by-light scattering contribution (l-l) is of similar size as , but has a much larger theoretical uncertainty Melnikov and Vainshtein (2004); Bijnens and Prades (2007); Prades et al. (2009); Nyffeler (2009)
where the uncertainty covers almost all values obtained in different publications.
Since this uncertainty is 30 times smaller than the experimental precision of (1), it can be regarded as known precisely.
The primary purpose of this letter is to report the complete numerical evaluation of all tenth-order QED contribution to . It leads to a sizable reduction of the uncertainty of the previous estimate by the leading-log approximations Kinoshita and Nio (2006a); Kataev (2006). We have also improved the numerical precision of the eighth-order QED contribution including the newly evaluated tau-lepton contribution. Together they represent a significant reduction in the theoretical uncertainty of the QED part of .
The QED contribution to can be evaluated by the perturbative expansion in :
where is finite thanks to the renormalizability of QED and can be written as
is independent of mass and universal for all leptons. , and are known exactly Schwinger (1948); Petermann (1957); Sommerfield (1958); Laporta and Remiddi (1996). Mass dependence is known analytically for and for Samuel and Li (1991); Li et al. (1993); Laporta (1993a); Laporta and Remiddi (1993); Czarnecki and Skrzypek (1999). We reevaluated them using the latest values of the muon-electron mass ratio and/or the muon-tau mass ratio Mohr et al. (2012). In the same order of terms as shown on the right-hand-side of (8), the results are summarized as follows:
The value of has been obtained mostly by numerical integration Laporta (1993b); Kinoshita and Nio (2006b); Aoyama et al. (2007); ?. They arise from 13 gauge-invariant sets whose representative diagrams are shown in Fig. 1. We have reevaluated some of them for further check and improvement of numerical precision. The results for the mass-dependent terms are summarized in Table 1.
From the data listed in Table 1 and the value of from Refs. Kinoshita and Nio (2006b); Aoyama et al. (2007); ?; Aoyama et al. (2012a), we obtain the following value for the eighth-order QED contribution :
Over the period of more than nine years we have numerically evaluated all 32 gauge-invariant sets of diagrams that contribute to Kinoshita and Nio (2006a); Aoyama et al. (2008b, 2010a, c); ?; ?; ?; ?; ?; ?; Aoyama et al. (2012a), whose representative diagrams are shown in Fig. 2. The results for mass-dependent terms are summarized in Table 2. Some simple diagrams were evaluated analytically or in the asymptotic expansion in Kataev (1992); ?; Broadhurst et al. (1993); Laporta (1994); Aguilar et al. (2008); Baikov et al. (2008). The results are consistent with our numerical ones.
From the data listed in this Table and the value of from Ref. Aoyama et al. (2012a), we obtain the complete tenth-order result:
The uncertainty is attributed entirely to the statistical fluctuation in the Monte-Carlo integration of Feynman amplitudes by VEGAS Lepage (1978). This is 20 times more precise than the previous estimate, , obtained in the leading-logarithmic approximation Kinoshita and Nio (2006a). This is mainly because we had underestimated the magnitude of the contribution of the Set III(a). Note also that (11) is about s.d. larger than the leading-log estimate. The numerical values of for are summarized in Table 3.
In order to evaluate (QED) using (7), a precise value of is needed. At present, the best non-QED is the one obtained from the measurement of Bouchendira et al. (2011), combined with the very precisely known Rydberg constant and Mohr et al. (2012):
where the first three uncertainties are due to the eighth-order term, tenth-order term, and the hadronic and electroweak terms, involved in the evaluation of . The fourth uncertainty comes from the measurement of . At present the difference between (12) and (13) is much smaller than the current uncertainty in the measurement of so that one may use either one of these two. However, some caution must be exercised to employ to calculate , when more accurate experiment of becomes available, because theoretical calculation of is strongly correlated with that of .
Note that the uncertainties of the lepton mass ratios, the eighth-order term, the tenth-order terms, and are improved by factors 1.7, 1.3, 20, and 1.5, respectively, compared with given in Eq. (99) of Ref. Jegerlehner and Nyffeler (2009).
In view of the rather large value of one might wonder how large might be. As a matter of fact it is not difficult to estimate its size. For this purpose note that the dominant contribution to comes from the Group IV(a) and the dominant contribution to comes from the Set VI(a). Both are integrals obtained by inserting several second-order vacuum-polarization loops into the virtual photon lines of the sixth-order diagram which contains a light-by-light scattering electron loop. Analogously the leading contribution to the twelfth-order term will come from insertion of three ’s in , namely,
noting that and the factor 10 accounts for the possible ways of insertion of . Including the contribution of other diagrams, the size of the 12th-order term might be as large as . This is larger than the uncertainty of the 10th-order term in (14) so that it would be desirable to obtain at least a crude evaluation of this term.
We have therefore
The size of discrepancy between theory and experiment has not changed much, since the tenth-order QED contribution is not a significant source of theoretical uncertainties. Let us emphasize, however, that the complete calculation of enables us to concentrate on improving the precision of the hadronic contributions.
Acknowledgements.We thank J. Rosner for a helpful comment. This work is supported in part by the JSPS Grant-in-Aid for Scientific Research (C)20540261 and (C)23540331. T. K.’s work is supported in part by the U. S. National Science Foundation under Grant No. NSF-PHY-0757868. T. K. thanks RIKEN for the hospitality extended to him while a part of this work was carried out. Numerical calculations were conducted on RSCC and RICC supercomputer systems at RIKEN.
- preprint: RIKEN-QHP-26
- G. W. Bennett et al. (Muon g-2), Phys. Rev. Lett. 92, 161802 (2004).
- B. L. Roberts, Chin. Phys. C34, 741 (2010).
- B. Lee Roberts (Fermilab P989 Collaboration), Nucl. Phys. Proc. Suppl. 218, 237 (2011).
- H. Iinuma (J-PARC New g-2/EDM experiment Collaboration), J. Phys. Conf. Ser. 295, 012032 (2011).
- T. Blum, Phys. Rev. Lett. 91, 052001 (2003).
- M. Gockeler et al. (QCDSF Collaboration), Nucl. Phys. B688, 135 (2004).
- C. Aubin and T. Blum, Phys. Rev. D75, 114502 (2007).
- X. Feng, K. Jansen, M. Petschlies, and D. B. Renner, Phys. Rev. Lett. 107, 081802 (2011).
- P. Boyle, L. Del Debbio, E. Kerrane, and J. Zanotti, Phys. Rev. D85, 074504 (2012).
- M. Della Morte, B. Jager, A. Juttner, and H. Wittig, JHEP 1203, 055 (2012).
- M. Davier, A. Hoecker, B. Malaescu, and Z. Zhang, Eur. Phys. J. C71, 1515 (2011).
- F. Jegerlehner and R. Szafron, Eur. Phys. J. C71, 1632 (2011).
- K. Hagiwara, R. Liao, A. D. Martin, D. Nomura, and T. Teubner, J. Phys. G G38, 085003 (2011).
- K. Melnikov and A. Vainshtein, Phys. Rev. D 70, 113006 (2004).
- J. Bijnens and J. Prades, Mod. Phys. Lett. A22, 767 (2007).
- J. Prades, E. de Rafael, and A. Vainshtein, in Lepton Dipole Moments, edited by B. L. Roberts and W. J. Marciano (World Scientific, Singapore, 2009) pp. 303–319.
- A. Nyffeler, Phys. Rev. D 79, 073012 (2009).
- K. Fujikawa, B. Lee, and A. Sanda, Phys. Rev. D 6, 2923 (1972).
- A. Czarnecki, B. Krause, and W. J. Marciano, Phys. Rev. Lett. 76, 3267 (1996).
- M. Knecht, S. Peris, M. Perrottet, and E. De Rafael, J. High Energy Phys. 11, 003 (2002).
- A. Czarnecki, W. J. Marciano, and A. Vainshtein, Phys. Rev. D 67, 073006 (2003).
- T. Kinoshita and M. Nio, Phys. Rev. D 73, 053007 (2006a).
- A. Kataev, Phys. Rev. D74, 073011 (2006).
- J. S. Schwinger, Phys. Rev. 73, 416 (1948).
- A. Petermann, Helv. Phys. Acta 30, 407 (1957).
- C. M. Sommerfield, Ann. Phys. (N.Y.) 5, 26 (1958).
- S. Laporta and E. Remiddi, Phys. Lett. B379, 283 (1996).
- M. A. Samuel and G.-w. Li, Phys. Rev. D 44, 3935 (1991).
- G. Li, R. Mendel, and M. A. Samuel, Phys. Rev. D 47, 1723 (1993).
- S. Laporta, Nuovo Cim. A106, 675 (1993a).
- S. Laporta and E. Remiddi, Phys. Lett. B301, 440 (1993).
- A. Czarnecki and M. Skrzypek, Phys. Lett. B449, 354 (1999).
- P. J. Mohr, B. N. Taylor, and D. B. Newell, (2012), arXiv:1203.5425 .
- S. Laporta, Phys. Lett. B312, 495 (1993b).
- T. Kinoshita and M. Nio, Phys. Rev. D 73, 013003 (2006b).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. Lett. 99, 110406 (2007).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. D 77, 053012 (2008a).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, (2012a), arXiv:1205.5368 .
- T. Aoyama, M. Hayakawa, T. Kinoshita, M. Nio, and N. Watanabe, Phys. Rev. D 78, 053005 (2008b).
- T. Aoyama, K. Asano, M. Hayakawa, T. Kinoshita, M. Nio, and N. Watanabe, Phys. Rev. D 81, 053009 (2010a).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. D 78, 113006 (2008c).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. D 82, 113004 (2010b).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. D 83, 053002 (2011a).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. D 83, 053003 (2011b).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. D 84, 053003 (2011c).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. D 85, 033007 (2012b).
- T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, Phys. Rev. D 85, 093013 (2012c).
- A. L. Kataev, Phys. Lett. B284, 401 (1992).
- A. L. Kataev, Phys. Lett. B710, 710 (2012).
- D. J. Broadhurst, A. L. Kataev, and O. V. Tarasov, Phys. Lett. B298, 445 (1993).
- S. Laporta, Phys. Lett. B328, 522 (1994).
- J.-P. Aguilar, D. Greynat, and E. De Rafael, Phys. Rev. D 77, 093010 (2008).
- P. A. Baikov, K. G. Chetyrkin, and C. Sturm, Nucl. Phys. Proc. Suppl. 183, 8 (2008).
- G. P. Lepage, J. Comput. Phys. 27, 192 (1978).
- R. Bouchendira, P. Clade, S. Guellati-Khelifa, F. Nez, and F. Biraben, Phys. Rev. Lett. 106, 080801 (2011).
- D. Hanneke, S. Fogwell, and G. Gabrielse, Phys. Rev. Lett. 100, 120801 (2008).
- D. Hanneke, S. Fogwell Hoogerheide, and G. Gabrielse, Phys. Rev. A 83, 052122 (2011).
- F. Jegerlehner and A. Nyffeler, Phys. Rept. 477, 1 (2009).