Analytical eighth-order light-by-light QED contributions from leptons with heavier masses to the anomalous magnetic moment of electron

The important consequences of the recent results of the numerical evaluations of eighth and tenth order QED contributions to the anomalous magnetic moment of electron are commented. The correctness of the results of the numerical evaluation of new eighth order QED corrections to the electron anomaly are supported by the demonstration of their consistency with the new analytical expressions for the QED contributions to $a_e$ from the diagrams with fourth-order light-by-light scattering muon and tau-lepton loops. The consistency of the similar results are demonstrated in the case of eighth order massive dependent contribution to the muon anomalous magnetic moment.

One of the most precise at present experimental results in the modern particle physics is the measurement of the electron anomaly a e = (g e − 2)/2 [1], [2], which gives a e = 1159652180.73(0.28) × 10 −12 The previous stage of the Standard Model theoretical prediction a th e = a QED e + a e (hadrons) + a e (weak) (2) was summarized in the most detailed review on the subject [3]. The perturbative QED contribution to Eq.(2) is defined as a QED e = A 1 + A 2 (m e /m µ ) + A 2 (m e /m τ ) + A 3 (m e /m µ , m e /m τ ) .
Up to recently theoretical and phenomenological applications were based on the following expression for its dominant term A 1 where l is the number of loops of the Feynman diagrams, which are contributing to the corresponding perturbative QED expression. Three first coefficients, presented in Eq.(4) in the numerical form, were evaluated analytically. The first term was calculated by Schwinger [4], the second correction was evaluated by Petermann [5] and Sommerfield [6]. This result was confirmed later on by Terentiev [7], who used different technique. The project of analytical evaluation of all three-loop QED contributions to a e was completed by Laporta and Remiddi [8]. The cited value for the four-loop QED correction to a e was numerically obtained in Refs. [9], [10] by Kinoshita and collaborators. The rough CODATA estimate of the coefficient of the 5-loop term in Eq.(4), namely A 1 | [11], gave the idea what might be theoretical uncertainties in the value of a e due to unknown up to recently total value of the tenth-order QED effects. The project of their direct numerical evaluation, started in 2005 by Kinoshita and Nio in Ref. [12], continued in the series of works of Refs. [13]- [21], was successfully completed in May of 2012 by Ayoama, Hayakawa, Kinoshita and Nio [22].
In general, the coefficient A is defined by the contribution of 12672 diagrams, which were classified into 32 gauge-invariant subsets. The result of evaluation of the final, most complicated Set V, together with more detailed numerical calculation of the eighth-order contribution to a e , was reported recently in Ref. [22]. The long-expected expression for the tenth-order massless contribution [22] is : while at the eight-order level the following new results were obtained [22]: The a (8) e -term differs from Eq.(6) due to the inclusion of complete mass-dependent contributions, numerically evaluated in the process of the works of Refs. [13]- [21]. These massivedependent effects, summarized in in the numerical form in Ref. [22], read  and A 2 (m e /m µ ) are presented in Table II of Ref. [22]. Note, that the sum of 10-th order massive dependent terms result in the small value of the overall massive correction A (10) 2 (m e /m µ ) = −0.00382(39) [22]. It should be stressed, that the determination of the concrete value of A (10) 1 -and a (8) e -terms is of real importance. Indeed, prior the work of Ref. [22] the scientific community faced with the unique case, when the comparison of the low-energy experimental result of Eq.(1) with the massless perturbative QED predictions of Eq.(4), supplemented with analytically known from the works of Refs. [25]- [31] massive-dependent forth and sixth-order contributions into A 2 and A 3 -terms and with the well-known values for a e (hadrons) and a e (weak) (see the review [3]), namely led to the value of the inverse fine coupling constant with theoretical uncertainties ±39 × 10 −8 , comparable with experimental ones ±33 × 10 −8 .
The dominant contribution to the theoretical error, namely ±30 × 10 −18 , was defined by the CODATA estimate of the coefficient of the O(α 5 )-term in Eq.(4). The additional sizable uncertainty in Eq.(13) came from the theoretical error of the numerical evaluation of the eighth-order contributions to Eq.(4) performed in Refs. [9], [10]. There were no additional theoretical errors in the sum of cited above analytically evaluated forth and sixth-order massive-dependent contributions A 2 and A 3 -terms. Their total contribution can be extracted from the summary part of Ref. [31] and reads: The errors in Eq.(14)-Eq.(16) are related to the indicated above uncertainties of the CODATA-2010 values of the ratios of leptons masses.
As to the uncertainties of the important contributions from a e (hadrons) and a e (weak), they were not sensitive at the previous stage of comparing theoretical and experimental predictions for a e .
The results of Ref. [22] and Eq.(5) allowed to solve this intriguing problem and to make theoretical uncertainties in the analog of Eq.(13) less important , than experimental ones.
Indeed, substituting new eighth-and tenth-order QED effects of Ref. [22] into the the procedure of the comparison with the precise experimental result of Eq.(1), the authors of Ref. [22] obtained more precise value of the inverse fine coupling constant: Here the first and second errors are related to the uncertainties of the numerical evaluation of the eighth and tenth-order QED corrections. The third error is determined by the combined uncertainties of the hadronic and electroweak contributions to a e , which start to manifest themselves at this more precise level of perturbative QED calculations, while the fourth and hugest uncertainty is determined by the experimental error in Eq.(1). In view of the importance of the results of the complicated eight and tenth-order numerical calculations it is highly desirable to perform the independent cross-checks if not all of them, but at least of some their parts.
In this work this problem is studied by analyzing the consistency of the numerical results for the massive -dependent contributions into A 2 (m e /m τ )-corrections to the electron anomaly, which are described by the subset of diagrams, formed by external light-by-light scattering muon and τ -lepton subgraphs with extra virtual photon, propagating inside this subgraph. The corresponding analytical expressions for the leading term of heavy-mass expansions of these contributions follow from the obtained in Ref. [32] result of analytical QCD calculations of the hadronic light-by-light scattering contribution to the lepton anomalous magnetic moments a e and a µ with hadronic effects, modeled by the internal light-by-light-scattering quark loop, crossed by virtual gluon. Taking into account that in QED one has C F = 1 and α M S = α(1 + O(α 2 )) (see Ref. [33] for details) we get the corresponding analytical contribution of the leading term in heavy lepton mass expansion of the the eighth-order QED correction to a e from this subset of massive-dependent light-bylight scattering graphs. The results contain the contributions of the Riemann ζ-functions ζ k = ∞ n=1 (1/n) k and polylogarithmic functions a k = Li k (1/2) = ∞ k=1 (1/2 k n k ) and read: where X 1+i = (m 1 /m 1+i ) with i = 1, 2 are the ratios of the electron mass m 1 = m e and the muon mass m 2 = m µ or τ -lepton mass m 3 = m τ . The abbreviations in the parenthesis labels the expressions for the subsets of graphs, which contain the next-to-leading order (NLO) approximation of the light-by-light (lbl) scattering subgraph. Substituting the known values for the transcendental and polylogarithmic functions into Eq.(18) we get: The CODATA-2010 values and errors of X 1 = m e /m µ and X 2 = m e /m τ are fixing the numerical expressions for the terms we are interested in : The results of the numerical calculations for the similar contributions are presented in Table 1 of Ref. [22] and read Taking into account the related uncertainties one can observe good agreement with the numbers of Eq. (19) and Eq. (20), which follow from analytical expression of Eq. (18). This can be considered as the strong check of the results of computer numerical calculations of eightorder contributions to a e , summarized in Ref. [22]. Moreover, minor discrepancy between Eq. (20) and Eq.(22) may indicate, that in the case of larger lepton mass (m τ ) the exact numerical results of Ref. [22] may be sensitive to taking into account higher terms of large mass expansion of the corresponding eight-order light-by-light-type diagrams, contributing to a e . Let us have a look whether the similar feature is manifesting itself in the case of the comparison of the results of analytical and numerical calculation of the similar eight-order light-by-light scattering corrections to the anomalous magnetic moment of muon a µ . In this case the expression , analogous to Eq.(18), reads Using the CODATA-10 value for X 4 = m µ /m τ we get This number is in good agreement with the result of the numerical eighth-order computer calculations, obtained in Ref. [34] in the process of completing numerical evaluation of tenth-order QED contributions to the muon anomaly a µ , namely with the result However, as in the case of τ -lepton light-by-light-scattering eighth order contributions to a e , the numerical comparison of the analytically-based result of Eq.(24) with the result of numerical calculations of Eq.(25) seem to indicate the sensitivity to still unknown higher terms of large mass expansions of the analytically evaluated massive-dependent Feynman graphs.
In any case our considerations demonstrates the reliability of the definite results of the complicated important numerical QED calculatiobns from Ref. [22] and Ref. [34]. More detailed considerations of new 10-th QED order results for a µ using the renormalization-group inspired studies of Ref. [35] may be done in future.
This work grew up from the detailed studies of the results for the eighth-and tenth-order QED contributions to the anomalous magnetic moment of electron and muon, started in January of 2012 during the visit to Theory Division of CERN. I am grateful to the members of Ph-Th Division of CERN for hospitality. It is the pleasure to thank A.E. Dorokhov, M. Nio, M.Passera and T. Kinoshita for discussions of various scientific subjects, related to this work. This work is done within the scientific program of the RFBR grants N 11-02-00112 and N 11-01-00182 and was supported by Grant NS-5590.2012.2.