Hadronic vacuum-polarization contribution to various QED observables

Due to precision tests of quantum electrodynamics (QED), determination of accurate values of fundamental constants, and constraints on new physics, it is important in a consistent way to evaluate a number of QED observables such as the Lamb shift in hydrogen-like atomic systems. Even in a pure leptonic case, those QED variables are in fact not pure QED ones since hadronic effects are involved through intermediate states while accounting for higher-order effects. One of them is hadronic vacuum polarization (hVP). Complex evaluations often involve a number of QED quantities, for which treatment of hVP is not consistent. The highest accuracy for a calculation of the hVP term is required for the anomalous magnetic moment of a muon. However, a standard data-driven treatment of hVP, based on a dispersion integration of experimental data on electron-positron annihilation to hadrons and some other phenomena, leads to a contradiction with the experimental value of aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_\mu $$\end{document}. This experimental value can be considered as an indirect determination of the hVP contribution to aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_\mu $$\end{document} and the scatter of theory and experiment allows one to obtain a conservative estimation of the related hVP contribution. In this paper, we derive exact and approximate relations between the leading-order (LO) hVP contributions to various observables. Using those relations, we obtain for them a consistent set of the results, based on the scatter of aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_\mu $$\end{document} values. While calculating the LO hVP term, we have to remember that next-to-LO (NLO) hVP corrections are often comparable with the uncertainty of the LO term. Special attention is payed to hVP contribution to simple atoms. In particular, we discuss the NLO contribution to the Lamb shift in ordinary and muonic hydrogen and other two-body atoms for Z≤10\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z\le 10$$\end{document}. We also consider the NLO contribution of the muonic vacuum polarization to the Lamb shift in hydrogen-like atoms. With the aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_\mu $$\end{document} puzzle unresolved, one may still require present-days values of the hVP contributions to various observable for comparison to experiment etc. the presence of contradicting values and a lack of consistency means an additional uncertainty for aμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_\mu $$\end{document} and for key contributions to it, including the LO hVP one. We present here an estimation of such a propagated uncertainty in hVP contributions to different QED observables and recommend a consistent set of the related LO hVP contributions.


Introduction
Precision low-energy tests of quantum electrodynamics are not entirely pure QED tests. Even in the case of purely leptonic systems, theory requires a certain hadronic input and, in particular, an input from hadronic vacuum polarization (hVP).
There is a number of QED observables for which we need such hVP contributions. The most important one is the anomalous magnetic moment of a muon a μ . A calculation of the leading hVP contribution to a μ at the level better than one percent is required for a comparison of theory and experiment (see, e.g., [1]). Other observables include the anomalous magnetic moment of an electron a e , muonium hyperfine (Mu HFS) intera e-mail: savely.karshenboim@mpq.mpg.de (corresponding author) val ν Mu (of the 1s state), the Lamb shift in hydrogen E L (H) for various states, the most important of which is 1s, and the Lamb shift in muonic hydrogen E L (μH) for the 2s state or, more accurately for the 2p − 2s interval. Contributions to the Lamb shift in other ordinary and muonic light atoms such as [μ]D, [μ] 3 He + , and [μ] 4 He + are also required for theoretical evaluation and interpretation of experimental data.
While the leading-order (LO) term for Δa μ (hVP) has been recently accurately calculated in many papers, the LO contributions for Δa e (hVP) ( [2,3]) and for Δν Mu (hVP) ( [2]) are found in a few publications only. There are also not very accurate calculations of the related contributions to Lamb shift in H [4] and μH (see, e.g., [5]). A number of results on the next-to-leadingorder (NLO) hVP contributions to a μ and LO and NLO hVP contributions to some other observables, which do not require a high accuracy, exist as well.
Even if the required accuracy of the other applications is not as high as that for a μ , there is still a question of a consistent evaluation of all the hVP contributions once one makes a 'grand QED evaluation' with many observables. Such an evaluation would be helpful for a determination of values of fundamental constants (see, e.g., [6]) or for a 'gross QED test' to constrain new physics.
Here we consider relations between the leading hVP contributions to various QED observables. Some of the relations are exact, others are approximate. All the corrections can be presented in a generic form, since they are proportional to an integral over a dispersion variable s. Here K(s) is a result of an integration over a QED part of the diagrams, usually expressed in terms of elementary functions, and is a quantity experimentally determined from e + e −annihilation to hadrons or extracted from similar phenomena. Accurate calculations may require interpolations or extrapolations of the existing data as well as theoretical partly model-dependent input in the area where the data are not very accurate or not available at all. Until recently the dispersion integration in (1) was considered as the only way to find LO hVP contributions, however recently lattice calculations started delivering results with comparable accuracy. The purpose of our work here is not to find an hVP contribution by itself, say, to Mu HFS interval from scratch, but to check how a change in R(s), due to a possible update of the experimental data or corrections, could affect the existing results. We take advantage of exact and approximate relations between K's for various observables which we derive in this paper. Such approximate relations may be not sufficient to find a contribution by itself, however for a small correction (compared to a previously obtained value) the approximate relations should be sufficient. They can be also used to consider the uncertainty due to a scatter of the data and discrepancies. One more application of the approximate relations is due to higher-order hVP contributions, for which the high accuracy is not required.
The variety of possible results for the R function is not entirely due to its successful improvement. There is a scatter of the experimental results obtained in different experiments and/or evaluated by different methods. Therefore, different existing parameterizations of the R function could be in part inconsistent.
There are two very different groups of problems in hVP calculations with the R function for various QED observables.
• For pure leptonic cases (a μ , a e , and ν Mu ), the hVP contribution should be taken into account. The required fractional accuracy for the leading hVP term differs from an application to application. The highest one is for Δa μ (hVP). One of the problems is an internal theoretical problem-the R function, determined from different groups of data or with help of corrections being approached differently, has inconsistencies. That should increase the uncertainty. The spread is present in the calculation of Δa μ (hVP) (see, e.g., [1]), but not for two other leptonic observables. For them the uncertainty due to the scatter of the results should be estimated.
A bigger problem is a discrepancy between theory and experiment for a μ . There are four options to resolve the discrepancy. One is 'new physics'. Once a model of new physics is introduced to force agreement between theory and experiment for a μ , it should be considered for the other observables. The second option is a possible problem with the experiment on a μ [7]. The third one is due to a possible incorrectness of the hVP calculations (questions about the required [raw] data for R(s) and corrections for measurements, say, radiative corrections). The fourth option is a problem with a contribution of hadronic light-by-light scattering (hLbL). Other possibilities are unlikely. For all the options without involvement of new physics, we can consider either corrections (one can correct various values to reach agreement or choose which of the values is a 'correct' one) or an increase in the related uncertainties. Note, in the case of Mu HFS interval we need not only Δν Mu (hVP), but also a value of a μ to calculate the leading term, and a certain hLBL correction to find one of the higher-order hadronic contributions. Their treatment should be consistent with the interpretation of the a μ discrepancy.
Let us consider one of the options-the a μ experiment is correct, while the hLbL contribution is reasonably estimated. Then, we can extract which would likely increase the other hVP values by a few percent (but less that for the related results 'from experiment') and set their uncertainty to a level of a few percent (essentially higher that for the results 'from experiment').
Use of either of those two values of Δa μ (hVP) for other observables should push the related hVP contributions for other observables beyond the present uncertainty of their direct calculations. A correct interpretation of any hVP contribution requires establishment of a relation between Δa μ (hVP) and the hVP contribution to an observable quantity of interest. • Observables with hadronic nuclei, such as the Lamb shift in ordinary and muonic atoms, should be treated differently. We cannot measure the nuclear radii, but we can extract their value from the electric-charge form factors. Meanwhile, we basically measure not the electric form factor G E (q 2 ), but a value G E (q 2 )/(1 − π h (q 2 )) (see, e.g., a discussion in [4]). Here π h (q 2 ) is a certain factor related to hadronic vacuum polarization. That means that any 'incorrect' treatment of the hadronic VP means rather an effective redefinition of the electric form factor and the related rms charge radius. At the present level of accuracy, such a 'redefinition' of the charge radius affects neither a determination of R ∞ , nor internal consistency of evaluation of data related to H, μH, D, μD, etc., as far as the hVP factor (1 − π h (q 2 )) is introduced in all the observables in a consistent way and as far as all the related nuclei are characterized by the data on the elastic electronnucleus scattering. In such a case, the hVP contributions can be completely ignored and that would not lead to any problem in the interpretation.
However, the situation with compound nuclei is not that simple as for the hydrogen nucleus, a proton. For a proton, one has no alternative to the scattering data and the measured quantity is expressed through G E (q 2 )/(1 − π h (q 2 )). For heavier nuclei, starting with the deuteron, one can construct an efficient nuclear model. In this case, the model 'speaks' rather directly in terms of G(q 2 ) and the hVP correction should be added. Often the radius itself is obtained from scattering, while a model describes various higher-order nuclear-structure effects, but still should produce a reasonable value of the radius as a consistency test.
Once the overall set of available data includes both H and D atoms and their muonic analogues, a certain consistent treatment of hVP is required.
Numerically, the mentioned effects are of minor importance for an evaluation of the current data; however, their importance may change with increase in the accuracy and increase in the nuclear charge Z. It is also important to be in control of those small corrections to be on a safe side.
With a dispersion-integral substitution for the Euclidean momentum q 2 one can present the result for a LO hVP contribution of the observable X i as Here, C i is a normalization constant (see Table 3) and m i is a characteristic mass for the observable X i . For X being the Lamb shift for an ns state in hydrogen (or other H-like atoms) and a e , m i = m e , while for a μ , ν Mu and the Lamb shift in μH (and other muonic atoms) m i = m μ . As mentioned, often a parameterization for the dispersion density ρ(s) in terms of the R function [see (2)] is utilized. The most important property of the expression (5) is that the QED kernels K i (s) are known [8][9][10][11] or can be relatively easily derived for the leading hVP contributions to various QED observables, while the R function can be expressed through a cross section of the e + e − annihilation to hadrons, which is measured (or through similar effects such as the rate of hadronic modes of τ decay). That makes the outcome for different observables X i correlated. Comparing various K i we are to identify those correlations.
The kernels for the leading-order hVP contributions are known. The result for the anomalous magnetic the kernel for the LO hVP contribution to the 1s Mu HFS interval (m = m μ ) (see Fig. 2) reads [11] (see also [12]) while the result for the Lamb shift in H (m = m e ) and μH (m = m μ ) (see Fig. 3) is obvious (Here and throughout the paper we apply units in which = c = 1.) Sometimes calculations (within a 'big evaluation') deal with different R functions for different observables X i (see, e.g., [6]). (Indeed, nobody applies by themselves different R's parameterizations within one paper, but frequently results, which have already been obtained in various original papers with different R's for different X's, are utilized within a single review or compilation without any reservations.) Such an evaluation is inconsistent. A consistent one would suggest the use of the same R for all the observables involved.
We intend to find the corrections to calculation of various hVP contributions. One of the purposes is to be able to correct and to make consistent the results for different X's for the case when the original evaluation has been performed with different R's. To do that we should be able to find what would happen with values of the hVP contributions to two observables if we [slightly] change the R function.
Having R(s), it is possible to find the LO hVP contribution to any observable directly. We intend to consider relations between the LO hVP contributions for two reasons. Indeed, the relations will simplify an adaptation of a result of a calculation on one observable to another. (We remind that there is quite a large number of different evaluations of LO hVP for a μ and a relatively few calculations for other observables. We prepare an 'instrument' to adjust any result for a μ to a calculation for another observable.) As we mentioned, a direct calculation is possible for all observables and such relations would only simplify the adaptation of direct results for a μ (LO hVP) to other quantities. In the meantime, we also have to deal with indirect results on a μ (LO hVP) which follow from a comparison of theory and experiment (see Eqs. 3 and 4)). Those indirect results are not related to any set of data on R (or related only in part). In this case, we need relations to convert a result on a μ to other observables. Similar situation exists with adaptation of lattice results a μ (LO hVP). Those are direct results, but they are not based on a dispersion presentation (6).
With the so-called a μ puzzle still unresolved, there is no consensus choice of a value of a μ . It is uncertain due to the discrepancy and the uncertainty is essentially larger than [claimed] uncertainties of individual theoretical or experimental determinations. Since we have no preferences on the resolution of the discrepancy, such an uncertainty should be considered as an uncertainty of all the involved key contributions including the one due to hVP. Once we consider Δa μ (hVP) more uncertain than its individual determinations [(cf. (4)], we have to acknowledge that until a consensus on the solution of the a μ puzzle is reached the experimental part of the dispersion integrand for Δa μ (hVP) [(cf. (6)] is compromised and an additional uncertainty should be assigned to it. That in its turn means an additional uncertainty for all calculations with R(s). We demonstrate below that certain relations between hVP contributions to different QED observables could serve this purpose. We stress that we consider an error in R(s) just as a possible solution of the problem on the same ground as the other possibilities (an experimental error in a μ , a theoretical error in hLbL, or new physics). Table 1 LO hVP contribution to various QED observables. Here, EF = 8/3 (Zα) 4 m 3 r /(memμ) is the leading nonrelativistic (NR) HFS term, which includes only the Dirac's values of the lepton magnetic moments (see Sect. 4); mr is the reduced mass of the related two-body system (muonium eμ, hydrogen ep, muonic hydrogen μp). Z is the nuclear charge which is equal to unity for all the mentioned two-body atoms; however, the result may be generalized to other hydrogen-like atoms with Z = 1 as we do for the Lamb shift in this paper

Xi
Ci Ki(mi, s) mi

Slope of the hadronic vacuum polarization Π h (0) and LO hVP contributions to various observables
As seen from Table 1, some kernels are of a simple form and in particular a few of them are ∝ m 2 /s. Let us introduce a notation for the slope of hVP polarization at zero The slope is given for the vacuum polarization as a function of the Euclidean momentum. The slope for the polarization as a function of the standard Minkovsky four-momentum has the opposite sign.
In the terms of the slope, we find The LO hVP for several other observables, such as the 1s − 2s frequency in positronium and its ground-state HFS interval, is also expressed in terms of Π h (0) (see Sect. A for detail). Being too small, the positronium LO hVP contributions are out of our interest here. Results of the hVP contribution to the Lamb shift in ordinary and muonic hydrogen, published in various recent compilations, are consistent, but based on old and relatively inaccurate calculations such as [4]. Meanwhile, a few of recent and accurate evaluations of the hVP contribution to a e have been available in the literature, which allows us to improve the Lamb shift results (see below). (We remind that our concern is not an accuracy of the calculations by itself, but their reliability (for the individual ones) and their overall consistency (for comparison of the results for different observables).) To find a value of Π h (0), we use calculations of the hVP contribution to g e − 2 in (11) The latter was recently found in a few publications [2,3]. The results of those direct calculations of the integral over the experimental data for the a e contribution are given in Table 2. In the same table, we give the extracted value for Π h (0) and the related results for the hVP contributions to the Lamb shift in H and μH.
The results from [2] to [3] are in agreement; however, they are based very much on the same data and we cannot treat them as uncorrelated. Because of that we choose their average as a central value and use their individual uncertainty as the uncertainty of the average. Note, at this stage we discuss a result of 'direct calculations' only, i.e., the result performed with an explicit value of the R function and for a quantity directly related to Π h (0).
We have studied above the most simple kernels from Table 1. The other are more complicated. To obtain simple relations, we may rely on their asymptotics. Even K ae (s) has a more complicated shape [see (7)] than it is given in Table 1. However, the dispersion integral starts 1 with s 0 = (2m π ) 2 , which leads to a  Table 3 The leading hVP contributions in the limit mμ, me → 0. Here we also include the cases where K(s) is simple and coincides with its asymptotics at s m 2 e , that is valid for the whole integration area. The corrections beyond the leading term due to the further m 2 e /s expansion are below 10 −5 for the all area of the s integration, so we may practically consider the identity (14) as an exact one, since the uncertainty of R(s) is always much higher than 10 −5 . For this reason, we have given the expression above as an exact one in Table 1.
A similar equation for the case of a μ is applicable, but as an approximate one (see Table 3). Its accuracy with actual s is roughly at 15% level for the most of the integration area (see the plot in Fig. 4). The lowest accuracy is at the lowest s, close to the 2π threshold. The error there is about 35%. However, that area does not contribute too much to the s integral, which is related to both its central value and the uncertainty.
To understand the accuracy of an approximation for the kernel K(s), we have to look at the function R(s), a convolution with which we are to calculate. Roughly related to NLO. One of such channels is with π 0 γ intermediate states, which starts at s = m 2 π . The effect of the π 0 γ channel is marginal and is ignored in this paper for the sake of simplicity. In higher orders, there should also be pure photon states produced via hadronic loops which means s0 = 0. Those effects are also negligible for practical applications.

Fig. 5
The R-function determined from the electronpositron annihilation to hadrons as published in our paper [5]). The plot was kindly provided to us by Simon Eidelman and we are grateful to him for his permission to use it in [5] and in this paper.) the behavior of the R function is presented in Fig. 5. The dominant contributions come from s m 2 ρ . Note, the ω and φ contributions are at s ≈ m 2 ρ from the point of view of the accuracy of the approximation with a smooth-behaving kernel K. Meanwhile (m μ /m ρ ) 2 0.02 which does not make the limit m μ → 0 sufficient for accurate calculations, but still makes it useful for a consideration of the correlations between different contributions.
Comparing the kernel K aµ and its asymptotics (see Fig. 4), we find that except of the area essentially below the ρ meson mass, the difference between the kernel K aµ and its asymptotics K aµ: lead can be described as A use of K aµ: lead allows us to express any a μ (LO hVP) result in terms of Π h (0). However, the accuracy of (15) may be not entirely sufficient for the whole area of s integration. We have to specifically consider an area of low s where the approximation is not that accurate. Another area where 15% in (15) is not sufficient is related to very high s. The ratio at s → ∞ is unity; however, the essential result comes from s in a wide area around m 2 ρ for which the unity is not achievable. As is well known for a while, the experimental result on a μ [7] does not agree with theory (see [1] for detail).
If any corrections to theory should be done, they are to be small. We intend to apply the relation in (15) rather for corrections to results (or differences between the results) than to the complete values. We are interested in the area below √ s = 2 GeV. A reason for that is that above this value one can not only measure ρ(s), but also calculate it within pQCD [14,15]. The results are in good agreement (see, e.g., [16,17]). We do not expect any big corrections there.
Applying that relation to the corrections should only allow us to maintain a high accuracy of the results. Since we consider possible corrections to theory, we have to choose theoretical values which will be 'corrected', i.e., our starting point. There are somewhat different theoretical results, but for analyzing the shifts it is not important which one to use. We take advantage that the accurate calculations were performed in [3] on a μ (LO hVP) and a e (LO hVP) in a consistent way and utilize the related results from [3] Δa μ (LO hVP) = 689. 5 Δa e (LO hVP) = 1.849(11) × 10 −12 (17) as our starting point. We remind that the latter one corresponds to Π h (0) m 2 e = 1.028(6) × 10 −6 (cf. Table 2). Denoting a shift (from the reference values) as δ, we find We intend to apply the relation (15) not to the complete contribution to a μ , which is a convolution with R(s), but to its possible changes due to changes in R(s). Such an approximation on K(s) with focussing on area relatively close to m ρ is valid for relatively large changes in R(s). 'Large change' means either a shift in a μ larger than or comparable with the uncertainty of a μ (hVP) or a serious reconsideration of the uncertainty by increasing or decreasing it by a factor, say, of two. Probability of all such changes should follow the value of the uncertainty of partial hVP contribution, i.e., by a product of K(s) and the [absolute] uncertainty R(s) [cf. (6)].  Table 5 The area, which dominates in the uncertainty budget, is expected to be the most probable area from which the change would come. That validates the approximation. Once we consider a small change in a μ , it may come from any area including those where the approximation in (15) is not applicable.

The a μ discrepancy and its consequences for various LO hVP
The main motivation to obtain the approximate relations above is the discrepancy in determination of a μ . One can interpret them as a discrepancy between a theoretical result, where the uncertainty comes from hadronic effects, and an experimental one [7], which may be considered as an indirect determination of the hadronic contribution.
Theoretical predictions for the LO hVP contribution to the anomalous magnetic moment of a muon have a certain scatter (see, e.g., [1,3,18]). The scatter is essentially smaller than the discrepancy of theory and experiment and we ignore it here. As a theoretical value, we use the result in (16) obtained in [3], which is chosen by us because in that publication a result on Δa e (LO hVP) is also presented. As mentioned above, a scatter of theoretical results on Δa μ (LO hVP) is not our concern in this paper, so we can start with any theoretical result on Δa μ (LO hVP).
The most realistic value, however, is neither theoretical nor experimental, but (4) found from scatter of the theoretical and experimental results (see Fig. 6). In the evaluations below, we give values of various hVP contributions, related to direct theoretical calculations; based on Δa from expt μ (LO hVP), extracted from experiment, and the one due to the scatter of the data. The latter is the estimation recommended in this paper.
On base of the approximate relations between the LO hVP contributions to a μ and a e [see (15)] and exact relations between the LO hVP contributions to a e , Π h (0), and to the Lamb shift in H and μH [see (10), (11), (12), (13), and Table 2], we find the results for the LO hVP contribution to various observables as summarized in Table 4 and Fig. 7.
Note that there is no consensus on the theoretical value. Results from scan over s in e + e − annihilation are currently less accurate (around the ρ-meson mass etc.). The data on τ -decay should be corrected due to violation of the isotopic invariance (introduction of which in [20] explained in the principle then discrepancy between decay data and annihilation by a ρ (0) −γ mixing). However, there exist a few different parameterizations for such a correction, which lead to somewhat controversial results [21][22][23]. Overall analysis sometimes utilizes τ -decay data (see, e.g., [3,23]), while in other cases it does not (see, e.g., [18,24]).
Use of recent and more accurate initial-state-radiation (ISR) data shows a certain discrepancy between different experiments (see, e.g., ISR data for the π + π − channel from KLOE [25][26][27], BaBar [28], and BESIII [29] and their comparison in [18]). Nevertheless, the scatter of theoretical data is essentially smaller than the discrepancy of theory with the experiment. Currently, it is not that important which theoretical value is used as input as far as eventually we plan to calculate the value of Δa from scat μ (LO hVP).  Table 2), indirect (from experiment and from scatter of theory and experiment), and lattice (from [19]) values for the LO hVP contribution to the anomalous magnetic moment of an electron and to Π h (0). The lattice value is discussed below In our analysis, in this paper we are limited to publications which consider the LO hVP contribution to [at least] two observables. One of such papers [3] presents results on the contributions to a μ and a e while applying τ -decay data, while the other paper [2], important for our consideration, gives results for a e and Mu HFS interval ignoring τ data. Consideration of both options made our analysis a representative one.
One has also to mention the result for LO hVP from lattice calculations. Until recently their accuracy was not comparable with that of theory and experiment. Significant progress in the field has produced a situation when the uncertainty of the lattice results is comparable with the uncertainty of Δa from scat μ (LO hVP). The lattice results are summarized in Table 5. The results are obtained by independent groups but using somewhat similar methods. Below we do not apply the lattice results on Δa μ (LO hVP) to obtain the LO hVP contribution to other variables, but consider them as a  [36] confirmation that Δa from scat μ (LO hVP), which we utilize here, is a reasonable estimation. Once the lattice results on a μ will reach a somewhat higher accuracy, they may be applied to other QED observables along with the technique that was developed in this paper.
There is also a lattice calculation for Δa e (LO hVP). The result [19] of 1.89(6)×10 −12 is in perfect agreement with a result from scatter in Table 4 (see also Fig. 7). Tables 1 and 3 on the LO hVP contributions contain results for a number of observables considered above in detail. One of important QED observables has not yet been discussed. This remaining case is considered here. It is for the muonium 1s HFS interval. Prior a discussion of the LO hVP term, we have to mention the leading NR contribution to the HFS interval, which is called the Fermi energy. It has several slightly different definitions, such as

LO hVP contribution to muonium HFS interval
where the upper version of the Fermi energy is a result of NR interaction of two full magnetic moments of an electron and a muon, while the one in the bottom is with their Dirac's values. It is convenient to consider the former as the leading term, while expressing the recoil effect (including hVP which is one of them) in units related to the latter. (That is not important for LO hVP, but a difference in definitions of Fermi energy as a prefactor for hVP affects definition of NLO hVP.) Indeed, the anomalous magnetic moments of leptons have their hVP contributions of order α 2 and the related contributions to the muonium HFS interval would be α 2 (Zα) 4 m e with various suppression of the mass factors. However, once the leading NR term, E F , is calculated with using the actual values of a e , a μ ,  (21) and (22). We remind that mρ 0.78 GeV, 2mπ 0.28 Gev those contributions are included in E F and we are interested in finding the hVP contributions beyond them. The related LO hVP contribution to the muonium 1s HFS interval is given in terms of (6) and (1) in Table 1 [see also (8)] [11]. It is of order α(Zα) 5 m e with a characteristic suppression mass factor of (m e /m ρ ) 2 . The related asymptotics at s m is of the form The results at the limit of m i → 0 are summarized in Table 3 and one can note that the asymptotics in (21) has a logarithmically enhanced term in contrast to all the other asymptotics considered above [cf. (14) and (9)]. That does not allow us to easily relate the results for the hVP contributions to a μ , a e and to the Lamb shift with that for the Mu HFS interval. To find such relations, we need an additional simplifying approximation.
To make such, we note that the dominant contribution comes from an area around s = m 2 ρ (see Fig. 5). Let us introduce The expression is less accurate as an approximation compared to the one in (21); however, it is easily related to the asymptotics of other kernels [cf. (14) and (9)]. As to the error of the approximation, it is plotted in   [2,3] or experimental result in [7]; the italic ones are from combinations of original data including the evaluation in this paper

Ref.
Ia µ I  (23) and the R function is well described by a step function obtained from perturbative QCD [14,15]. The latter circumstance makes reliable the high-s end of the s integration in the dispersion integral, because of consistency between experimental data and perturbative QCD theory. Due to that, we expect that the major shift in R(s) (if any) may happen between the 2π-threshold and, say, √ s = 2 GeV. That allows us to write for the difference between a result of a direct calculation and one with a true R(s) function, based on certain value of a μ (which may be obtained from experiment, estimated from scatter of experiment and theory, or taken from a lattice calculation), as an indirect estimation, the following We remind that we are interested in relatively large changes in a μ (see discussion above).
As 'reference' points (to calculate shift δI Mu ) we consider the results from [2] to [3]. In [2] they found a e (LO hVP) and ν Mu (LO hVP), which allows one to relate a result on I ae to the corrected value of I Mu . Since (23) sets a relation between shifts in I aµ and I Mu , we have to find first a shift in I aµ from a known shift in I ae . The results of [3], which include Δa e (LO hVP) and Δa μ (LO hVP), are suitable for that (see the previous section for details). In terms of the I integrals, the results of our evaluations are summarized in Table 6.
The results due to the a μ -discrepancy are

NLO hVP contribution to the Lamb shift in ordinary and muonic hydrogen
The NLO hVP contributions to a number of QED observables are known, but not to all of them. In particular, the NLO contribution for a μ is known [23,37] (see also, e.g., [1]) as well as the ones for a e (see, e.g., [2,3,37]) and Mu HFS interval [38,39]. (For some variables such as a μ , a e even NNLO hVP is known [40]; for a μ , a e , and Mu hfs interval, the hLbL contributions are of the same orders as NLO hVP and needed for a complete NLO hadronic term (see, e.g., [1,3,23,41,42] for detail).) Here we consider the NLO hVP contribution to the Lamb shift in ordinary and muonic hydrogen. NLO in H and μH includes a term, which is due to the twoloop vertex with one of those two loops being hVP (see Fig. 9). It is of relative order α compared to the related LO hVP term. In the case of the Lamb shift of an ns state in hydrogen that is a complete NLO correction.
That NLO term is of order α 2 (Zα) 4 m(m/m ρ ) 2 . It comes from a slope of the Dirac form factor F 1 (0) at zero momentum transfer and from the value of the Pauli form factor F 2 (0) at zero momentum due to the mentioned two-loop vertex, depicted in Fig. 9. The latter, F 2 (0), is the anomalous magnetic moment of the orbiting lepton. This contribution relates to LO hVP to the related anomalous magnetic moment and it is indeed known. However, since two contributions to the Lamb shift of an s state (due to F 1 (0) and F 2 (0) ), which are technically very similar, we consider both of them as a single corrections within the same formalism (see below).
There is also a contribution to the Lamb shift at l = 0 in ordinary and muonic hydrogen of the same order as the mentioned NLO hVP contribution for the s states. It is completely described by the anomalous magnetic moment in the same matter as the leading QED term for l = 0, once we deal with a complete value of the anomalous magnetic moment of the orbiting lepton a L ΔE (a) (25) To find the LO hVP contribution to the Lamb shift of l = 0 states, which is of order α 2 (Zα) 4 m(m/m ρ ) 2 , i.e., of the same order as the NLO hVP contribution for s states, it is sufficient to utilize the LO hVP contribution of the anomalous magnetic moment, Δa L (LO hVP), into the identity above. Since this contribution is known (see above), we do not consider it in this paper and focus on the Lamb shift in ns states.
The integral presentation (in terms of s) for the vertex contribution is of the same shape in H and μH. While that is a complete result for hydrogen, for its We remind that in H and μH the LO hVP contribution to the Lamb shift (see Fig. 3) is proportional to |Ψ (0)| 2 , the squared value of the wave function at origin. In the muonic case, there is a number of different sources for the α correction to |Ψ (0)| 2 in addition to the two-loop-vertex term discussed above. One of them is from an NR correction to |Ψ (0)| 2 due to the Uehling potential which is known in a semi-analytic form (see, e.g., [43]), while the other is due to an NR correction with electronic vacuum polarization (eVP) on the same photon as hVP. Those NR specific μ-atom NLO hVP corrections are presented in Fig. 3. The results could be generalized for light two-body muonic atoms. The Uehling correction to the wave function at origin is given for the 1s and 2s states of atoms A ≤ 10 in [44].
Let us consider now all the mentioned effects in more detail, starting with the two-loop vertex. A two-loop vertex, where one loop is a VP one, delivers pure lepton contributions as well as hVP ones. The simplest leptonic one is the eVP contribution to the Lamb shift in hydrogen, which is well known [45][46][47][48]. A similar contribution is due to muonic vacuum polarization (μVP) in μH (which can be found with a simple scaling substitution m e → m μ ). There are also asymmetric in masses pure leptonic contributions, such as eVP in μH and μVP in H, which are known only in part.
Let us find C vert2 (m) and K vert2 (m, s) for s states, defined in general in (6), where m is the mass of the orbiting particle. Following [49], we obtain for the Lamb shift of the ns state with an orbiting particle with mass where The identities above are for the sum of two mentioned earlier contributions (due to F 1 (0) and F 2 (0)).
Here s is a dispersion variable [cf. (5)], which can be introduced both for leptonic and hadronic VP loops. In the leptonic case, we use the parameterization where m l is the lepton mass in the VP loop.
Performing the x-integration, we find where It is useful to consider asymptotic behavior at s m 2 . The result reads Using the kernel obtained above, we can consider various leptonic VP contributions to the Lamb shift of s states in H and μH. Such contributions are eVP and μVP for H and eVP and μVP for μH. The contribution, for which the loop particle and the orbiting one are the same, agrees with the known theoretical result [48].
which agrees with the result from [50].
The μVP correction to the Lamb shift in hydrogen (see Fig. 11), which is a new result, is obtained as There are different possibilities to obtain a fast preliminary estimation. In particular, we consider two simple methods below. One is to use a simplistic model for the dispersion function R(s), which we have previously used for other occasions (cf. [38,39]). A related result is given in Table 7. For μH, we have applied an exact kernel, since the asymptotics in (32) is not very accurate for μH. Therefore, our estimation is obtained with an exact kernel K(s) and approximate function ρ(s).
Another estimation can be obtained with complete R(s) but very approximate K(s). To begin with, we For numerical estimations, which are given in the first line of Table 7, we utilize results from [3] for the righthand side of the identities in (33). We have obtained two independent estimations for the Lamb shift in both ordinary and muonic hydrogen, and their scatter characterizes the uncertainty of those estimations.
The estimation can be improved once we use an approximation of the kernel as a combination of kernels for Mu HFS interval and for a e (or a μ ) with the coefficients good for asymptotics s m 2 . In this case, we can use the identity (32) with ln(s/m 2 ) without a substitution s → m 2 ρ , which is not very accurate for a broad range of s.
In particular, one can express the asymptotics of this kernel in terms of asymptotics for the kernels for a μ , a e and for Mu HFS interval as Using the same coefficients, we apply now the relations for the complete kernels, which makes those relations approximate. The result for their integration over s is known (see above) and, following the scatter results, we obtain To check the accuracy of (35) and introduce corrections, we have studied the ratios of interest. In previous sections, we were interested in small corrections to known values. It was sufficient to study corrections to R(s) in a relatively narrow area. Here we are to calculate the whole NLO contributions and therefore we are rather interested in the complete area of s integration. The accuracy for the contribution to the Lamb shift in hydrogen is determined by the ratio Using the relations above, we obtain central values given in Table 7. We still do not give any uncertainty. The most uncertain estimation can be obtained from the scatter value. Even its uncertainty is smaller that the uncertainty of the approximate relations we apply. The uncertainty of the NLO: vert2 term is much smaller than the one of LO hVP and can be ignored.

The hVP contribution to the Lamb shift in hydrogen and muonic hydrogen
The NLO hVP contributions have been known for a μ (see, e.g., [18]), a e (see, e.g, [2]), and Mu HFS interval [38,39], but not for the Lamb shift in both ordinary and muonic hydrogen. In principle, any accurate evaluation of LO should be accompanied with an NLO calculation, since the LO uncertainty is often comparable to the NLO contribution. The NLO contributions to Lamb shift in H and μH have been already in part considered above.
In the case of hydrogen, NLO consists of the vert2 term only. It is of relative order of α to the LO hVP term. A contribution with an extra Zα exists, but it has an additional suppression factor of m e /m ρ . Theory of hVP for the 1s Lamb shift in hydrogen is summarized in Table 8. We consider NLO not as the result by itself, but rather to compare it with the uncertainty of LO. For the 1s state in hydrogen it is 18 Hz.
The NLO contribution is the same (in units of the LO term) for any s state. Since the NLO contribution has an additional factor of α, it is the same in units of the   LO term for all the low-Z H-like atoms. In particular, for an arbitrary s state of a low-Z H-like atom it is The LO hVP term also scales as Z 4 /n 3 . Theory of the Lamb shift in muonic hydrogen is somewhat different from that in hydrogen. In particular, additionally to the vert2 term there are NR Uehling contributions. One of them is due to a modification of the wave function at origin, that has been considered for a number of occasions (see the left diagram in Fig. 10). The other one is with eVP on the same photon as hVP (see the right diagram in Fig. 10). Both are similar to the corrections in muonic atoms for HFS interaction and for finite-nuclear-size contributions (FNS) to the Lamb shift (cf., e.g., [44,51]). We apply here the related results following [44]. The evaluation of soft contributions due to the Uehling potential is given in Appendix in Sect. B for the nuclear charge Z = 1 − 10.
Summary on theory of hVP for the 2p−2s Lamb shift in μH is given in Table 9. Similarly to the case of hydrogen above, here we consider NLO not as the result by itself, but rather to compare it with the LO uncertainty. The results on other low-Z two-body muonic atoms are given in Appendix in Sect. B.
As mentioned in a number of publications (see, e.g., [4,5,52]), there is a certain interference in treatment of the FNS and hVP contributions to the Lamb shift in H and μH (and other atoms). In a number of situations, such as, e.g., in determination of the form factors, the hVP correction is not subtracted while evaluating the scattering data and therefore instead of form factors one measures rather For the expansion at low q 2 , that means that instead of the charge radius squared we effectively deal with If only the LO hVP contribution to the Lamb shift of either ordinary or muonic atom is important, it is more advantageous not to consider hVP explicitly at all, which means an effective redefinition of the rms nuclear charge radius following (40). Mathematically speaking, FNS and hVP do not make a single block. While fitting, they have a somewhat different structure and it may be good to separate them. Corrections, e.g., a soft QED correction for μ atoms, are also not of exactly the same form. The combinatoric coefficient is different (see Sect. B in Appendix). A similar situation is with the radiative corrections to electron-proton scattering, e.g., eVP corrections for one-photon exchange have a different (by a factor of two) coefficient for form factor and for hVP. Because of logarithmic enhancement for q 2 m 2 e , such a correction can be more important for scattering than for the bound problem.
In the case of light compound nuclei (deuteron, triton, helion, α-particle), the situation is more tricky. The radius can be determined either from muonic atoms or from scattering. In the meantime, there are various nuclear models (see, e.g., [53]) which may be efficiently used for calculation of the nuclear-structure contributions. A value of the radius used in theory of compound nuclei is usually not the one obtained from the models, but they still have to produce such a radius from the model for a consistency test. The hVP effects are definitely to be added to nuclear-physics-model calculations of the radius. However, at the present level of accuracy that is not really important. The most accurate applications relate H and μH and their isotopes, where the mentioned effects are rather marginal.
Here, we have calculated hVP corrections in ordinary and μ atoms and it is up to those who consider their complete theory to include them or not. However, the decision should be explicitly explained in order to avoid confusion and to allow ones to maintain consistency with other calculations.

Conclusions
In area of precision low-energy tests of QED, there is a number of discrepancies, such as the proton-radius puzzle and a controversy with a μ . These discrepancies create problems at two different levels. One that attracts most of the attention is roughly speaking a question, which of controversial results are right and which are wrong. In meantime, there is an additional problem in the shadow of the previous one. Presence of a controversy means that certain data and certain standard treatments are compromised. Meanwhile, they are not applied exclusively to the controversial quantities, but are often also used for a number of others.
For instance, the proton-radius puzzle tells us that probably the use of the data on the form factors of a proton in a wide region, as it is usually done, is incorrect. (If it is, it is not clear whether we have a problem with the data or their interpretation, including the fitting 'as usual'.) That means not only that the related value of the proton charge radius is [probably] wrong, but also that the magnetic radius and various other quantities such as the so-called Zemach and Friar radii, cannot be considered as reliable, once they are deduced from the same evaluation that delivers a probably incorrect value of R p . The Zemach radius [54] is required to obtain theoretical predictions for the HFS interval in ordinary and muonic hydrogen, while the Friar radius [55][56][57] is a key element of a presentation of a higherorder FNS contribution, which plays an important role for a determination of R p from the Lamb shift in μH.
(The related higher-order FNS contribution is essentially smaller than a discrepancy for the radius from μH spectroscopy and e − p scattering, but essentially larger than the uncertainty of R p due to the μH experiment [58] and of a pure QED part of the related theory (see, e.g., [5]).) See [59,60] for more details.
A number of recent indications from hydrogen spectroscopy [61,62] and electron-proton low-momentumtransfer scattering [63] have pointed out that the muonic value of the proton charge radius [58] is in general correct. (Not all the recent results support the muonic value of R p , e.g., the spectroscopic result [64] confirms the 'old' hydrogen value, which is consistent with a result of a [former] 'usual' evaluation from e − p scattering [65].) However, they do not help for understanding what to do with the overall scattering data and their fits for the electric form factor, that is required for a calculation of the Friar term in order to eventually obtain a muonic value of R p .
The problem with a μ is somewhat different. If theory is compromised, that means that a calculation of LO hVP is compromised at the first place. (The [dominant] pure QED part of theory [66] is much more accurate and reliable than the LO hVP contribution, while the higher-order hadronic effects are smaller than LO hVP and therefore also more reliable (see [1] for detail).) Standard theoretical evaluations of the LO hVP term are based on the dispersion integral in (1) over experimental data (presented with the R(s) function). The same function R(s) is applied for other observables.
The discrepancy for a μ should be considered as the most accurate test of calculations of the LO hVP contribution through a dispersion integral. The discrepancy indicates that the dispersion calculation is compromised. It is the most important not just to update various hVP results from an occasion to an occasion, but to check whether a more recent value is consistent with the previous one, and to be sure that all the results for different observables, applied within a certain analysis, are consistent.
Disagreement between experiment and theory on a μ opens a question, how the hVP results for a μ should be interpreted. There are four options: (i) experiment [7] is not correct; LO hVP theory (see, e.g., [3]) is not; the evaluation of hLbL is not; or certain new physics is present.
Technically, the option of new physics as a solution of the discrepancy between theory and experiment leads to the same consequences for the application as a possible [relatively] large error in the hLbL contribution. The experimental value of a μ is applicable in such a case as well as the R function and all the direct results on the LO hVP contributions to various observables; however, one has to consider an additional contribution to various QED observables, which has not been accounted for (or has not been correctly accounted for) previously. Such a contribution (e.g., due to a conceptual error in the hLbL term or due to new physics) for an observable of interest may happen to be larger than the uncertainties due to LO hVP (theoretical) and a μ (experimental). In both cases (hLbL and new physics), the contribution is to be model dependent. The required new physics should be introduced with a certain model or one has to consider a completely new approach to hLbL, drastically different from the previous calculations. Eventually, that leads to a bunch of model-dependent corrections to various observables. Until such a model appears we have to focus our attention on two other options: either the a μ measurement [7] has a systematic error, or the data on R(s), applied for a standard calculation of LO hVP, have to be reconsidered. The data may involve the measurements errors, the errors in interpretation (such as radiative corrections), or models to extend the data beyond the experimental reach (to low-luminosity kinematic areas and/or channels).
Without new data or a hint, we are not capable to choose between hVP theory or experiment 2 , and therefore as an estimation of the accuracy of LO hVP we consider the scatter of the data from its direct and indirect determination. The estimation, based on the scatter, pretty well agrees with the lattice results [19,[30][31][32][33][34][35][36] (see Fig. 6). As a consistent set of hVP contributions, we choose the values given in Table 10. They are based on scatter of a μ (see Fig. 6) and details of their evaluation have been given in the paper. The muonium 1s HFS interval deserves a somewhat more accurate consideration. Technically speaking, the largest hVP contribution to Mu HFS interval is not the one considered here as the LO hVP term, but another one due to a μ . They are of the same order in α but have different suppression mass factors, such as m 2 μ /m 2 ρ (for a μ ) and m e m μ /m 2 ρ (for a 'direct' contribution to Mu HFS). A small LO hVP contribution to a e is also of the same order in α but has a suppression factor of m 2 e /m 2 ρ . Once we define the leading NR term of the muonium HFS interval, the Fermi energy E F , as the result of interaction of full moments of electron and muon as given in (19), the a μ -related hVP contribution is a part of the Fermi energy, while the leading hVP correction to it is Δν Mu (LO hVP). However, we have to consider the leading NR term and the LO hVP correction to the Mu HFS interval in a consistent way. Since QED theory of a μ is set, any uncertainty and error is understood as uncertainty or error of direct or indirect determination of the contribution of hVP effects. A simultaneous inclusion in the same theoretical expression of an experimental value of a μ from [7] and Δν Mu (LO hVP) from a dispersion integration (1) (such as the result from [2]) is inconsistent.
Our recommended set of consistent values for the involved quantities, based on the scatter for a μ , is where we also give a correlation coefficient.
In Table 10, we present for the Lamb shift in hydrogen and muonic hydrogen a sum LO+NLO hVP, since NLO contribution is an original result of this paper. One may be interested in LO and NLO separately and not only for ordinary and muonic hydrogen, but also for other atoms. The NLO hVP in ordinary atoms and a hard part of the contribution to muonic atoms is discussed in Sect. 5, while the soft part of the muonicatom correction is present in Appendix in Sect. B. Here we summarize the results for A = 1 − 4 for ordinary and muonic atoms that are given in Tables 11 and 12, respectively.
One may also be interested in their intercorrelations. As concerns correlations for the Lamb shift of the ns states in ordinary and muonic atoms, all the LO hVP contributions are proportional to the same value, namely Π h (0) and therefore they are correlated with a correlation coefficient r = 1. (We remind that the LO hVP contribution vanishes for l = 0.) We do not consider here few-electron atoms, such as He, or simple molecules, such as HD + ; however, we expect that LO hVP for most of low-energy variables is to be described with contact interactions (an electron-electron or electron-nucleus one) and therefore it is to be proportional to Π h (0), which makes our finding directly applicable for such atomic and molecular systems. The exception is the hyperfine interaction (cf. Sect. 4) for which correlations make a more complicated form.
Having in mind various projects on accurate measurements of pure leptonic properties, such as the anomalous magnetic moment of electron [68] and muon [67,69], energy intervals in muonium [70,71] and positronium [72,73], progress in related theoretical QED calculations [74][75][76][77][78], and in determination of involved fundamental constants (such as the fine structure constant α [79]), a better understanding of related hVP contributions, that is one of the limiting factor of theory of pure leptonic systems, is important.
After the work was completed a number of new publications become available [80][81][82][83][84][85][86]. They somewhat improve the situation on certain issues but do not affect the conclusions of the paper. One of the results of this paper is the relations between LO hVP contributions to different observable. The data may play a role of starting point. The most important problem with the data is a discrepancy of theoretical and experimental results on a μ . A minor change of a dispersion result for the theoretical value is marginal comparing with the size of the discrepancy. Once we deal with a half sum of theoretical and experimental values as the central value of our prediction and their half difference as an estimation of the uncertainty, any change in actual theoretical values is negligible.
Still, new data-driven evaluations of the LO hVP contribution to a μ was presented in [80][81][82]. Two of them [80,81] utilize e + e − annihilation data only, while the other [82] uses also the τ -decay data. The consensus on use or non-use of τ data has not yet reached.
There is a new result on LO hVP from lattice [83]. There has been also a progress in calculation of hadronic LbL contribution [84][85][86]. Within our logics the latter contributes to the indirect determination of the LO hVP by striping the experimental value of a μ out of all the other theoretical contributions. The hLbL one is a dominant contribution to the theoretical uncertainty of such a procedure.
Technically, the results from [80] could be useful for realization of our strategy line, because that seems to be the first paper with appropriate description of R(s) where the LO hVP contribution is calculated with the same R(s) for a μ , a e , and Mu HFS. That should sim- plify the evaluations, however, the eventual results (see Table 10) would not change much.
Funding Open Access funding enabled and organized by Projekt DEAL.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All the data we used are published elsewhere and we give the related citations. We present the most important input data with appropriate citations explicitly in our paper, while our output results are also explicitly given in various tables of our paper.].  also a suppression factor of (me/mρ) 2 , where mρ represents the scale of hadronic masses in the dispersion integral of (6) where most of the contribution comes from.) The α 5 me hVP contributions are exclusively for the S states, being proportional to the squared value of the wave function at origin. They are also proportional to Π h (0), which makes the hVO positronium results to be of interest for this paper.

Open Access
One of the LO hVP contributions in Ps is a standard Lamb-shift-type term (see Fig. 3) (cf. Table 1; we remind that in positronium m = me, mr = me/2) The other LO hVP contribution is one to the energy of a triplet S state (see Fig. 14) where is the annihilation-channel contribution to the leading term of the Ps 1S HFS interval. The triplet-state contribution affects both the HFS interval and the energy of center of mass of the nS states. An accurate result on the LO hVP contribution to the Ps 1s HFS interval was published in [87]. The result is related to m 2 e Π h (0) = 1.027(2) , which has a smaller uncertainty than those based on calculations in [2] and [3] (see Table 2), while the latter have been obtained from more sophisticated evaluations. The central Table 13 The numerical results for the wave-function correction for the 1s and 2s states in light muonic atoms (Z ≤ 10) (cf. [51,94]) following (B3). See also [44] Nucleus κ Z C value of the result in [87] is reasonable, while the estimation of the theoretical uncertainty in [87] is considered by us rather as an overoptimistic one. Besides, as we explain in this paper, the very use of a purely theoretical result based on the dispersion integration in (1) is questionable because of a discrepancy between theory and experiment on aμ.

Appendix B: NLO hVP contribution to the Lamb shift in two-body muonic atoms with Z ≤ 10
Let us make a brief 'parallel' calculation of the Uehling correction to the leading FNS term (following [44]) and the LO hVP term in a muonic atom. It is useful to present both contributions in the form (B2) It is generic for any δ-function kernels, which is the case for the leading terms for FNS and hVP (see [44]). The contribution is present only for the ns states. The other correction, C X: kern 1 is a correction to the δfunction-like kernel and, in principle, it depends on the kernel. In general, it may not vanish for l = 0 states in muonic atoms and, in particular, for both FNS and hVP it does not.
The numerical results on the coefficient are summarized in Table 13 for the 1s and 2s states. As we already mentioned, C 1 (nl) = 0 for l = 0 and, in particular, C 1 (2p) = 0. The correction due to the modification of the contact term (i.e., the kernel of the interaction) is specific for any contact term. In the case of the FNS contribution in order α(Zα) 4 m [43] (cf. [51,57,95] and [96][97][98][99][100][101][102]) the result is of the form where A(κ) has been defined above. The result for the hVP: kern term is similar to the known FNS: kern one, but has a different combinatoric coefficient for the 1s Lamb shift and for the n = 2 one (i.e., for the 2p−2s interval) in light muonic atoms are given in Table 14.
A complete result for the hVP contribution (LO+NLO), which includes the Uehling correction discussed in this appendix and the 'hard' NLO: vert2 term considered in Sect. 5 (cf. Table 9) for the 2p − 2s Lamb shift in lightest muonic atoms, is summarized in Table 12 of Sect. 7.