Electron-positron annihilation into hadrons at the higher-loop levels

The strong corrections to the R-ratio of electron-positron annihilation into hadrons are studied at the higher-loop levels. Specifically, the essentials of continuation of the spacelike perturbative results into the timelike domain are elucidated. The derivation of a general form of the commonly employed approximate expression for the R-ratio (which constitutes its truncated re-expansion at high energies) is delineated, the appearance of the pertinent $\pi^2$-terms is expounded, and their basic features are examined. It is demonstrated that the validity range of such approximation is strictly limited to $\sqrt{s}/\Lambda>\exp(\pi/2) \simeq 4.81$ and that it converges rather slowly when the energy scale approaches this value. The spectral function required for the proper calculation of the R-ratio is explicitly derived and its properties at the higher-loop levels are studied. The developed method of calculation of the spectral function enables one to obtain the explicit expression for the latter at an arbitrary loop level. By making use of the derived spectral function the proper expression for the R-ratio is calculated up to the five-loop level and its properties are examined. It is shown that the loop convergence of the proper expression for the R-ratio is better than that of its commonly employed approximation. The impact of the omitted higher-order $\pi^2$-terms on the latter is also discussed.


Introduction
In the studies of a variety of the strong interaction processes a key role is played by the hadronic vacuum polarization function Π (q 2 ), the related function R(s), and the Adler function D(Q 2 ). In particular, these functions govern such processes as the electron-positron annihilation into hadrons, inclusive hadronic decays of τ lepton and Z boson, as well as the hadronic contributions to precise electroweak observables, such as the muon anomalous magnetic moment (g − 2) µ and the running of the electromagnetic fine structure constant. The theoretical analysis of these processes constitutes a decisive self-consistency test of Quantum Chromodynamics (QCD) and entire Standard Model, that, in turn, puts robust restrictions on a possible New Physics beyond the latter.
Additionally, the energy scales relevant to the foregoing strong interaction processes span from the infrared to ultraviolet domain, so that their theoretical investigation provides a native framework for a profound study of both perturbative and intrinsically nonperturbative aspects of hadron dynamics. It is worth noting also that a majority of the aforementioned processes are of a direct relevance to the physics at the currently designed Future Collider Projects, such as the Future Circular Collider FCC-ee [1], Circular Electron-Positron Collider CEPC (its first phase) [2], the International Linear Collider ILC [3], the Compact Linear Collider CLIC [4], as well as the E989 experiment at Fermilab [5], the E34 experiment at J-PARC [6], and others.
In fact, over the past decades the perturbative approach to QCD remains a basic tool for the theoretical exploration of the hadronic physics. However, the QCD perturbation theory can be directly applied to the study of the strong interaction processes only in the spacelike (Euclidean) domain, whereas the proper description of a e-mail: nesterav@theor.jinr.ru arXiv:1707.00668v2 [hep-ph] 18 Dec 2017 hadron dynamics in the timelike (Minkowskian) domain additionally requires the pertinent dispersion relations.
Specifically, the dispersion relation for the R-ratio of electron-positron annihilation into hadrons converts the physical kinematic restrictions on the process on hand into the mathematical form and determines the way how the "timelike" observable R(s) is related to the "spacelike" quantity D(Q 2 ), the corresponding perturbative input being embodied by the so-called spectral function. Since the calculation of the latter at the higher-loop levels constitutes a rather challenging task, one commonly resorts to an approximate form of the R-ratio, namely, its truncated re-expansion at high energies. At the same time, one has to be aware that at any given loop level such re-expansion generates an infinite number of the so-called π 2 -terms (which may not necessarily be small enough to be safely discarded at the higher orders), that also worsen the loop convergence of the resulting approximate R-ratio.
The primary objective of the paper is to explicitly derive a general form of the spectral function required for the proper evaluation of the R-ratio and to study its properties up to the five-loop level. It is also of an apparent interest to calculate the R-ratio itself, to examine its higher-loop convergence, and to elucidate the impact of the omitted higher-order π 2 -terms on its truncated re-expanded approximation.
The layout of the paper is as follows. In Sect. 2 the essentials of continuation of the spacelike perturbative results into the timelike domain are expounded. In Sect. 3 the one-loop expression for the R-ratio is explicated and its approximations are discussed. In Sect. 4 the derivation of a general form of the commonly employed approximate expression for the R-ratio (which constitutes its truncated re-expansion at high energies) is delineated, the appearance of the pertinent π 2 -terms is elucidated, and their basic features are studied. In Sect. 5 the explicit form of the spectral function required for the proper calculation of the R-ratio is obtained and its properties at the higher-loop levels are examined. By making use of the derived spectral function in Sect. 6 the proper expression for the R-ratio is calculated up to the five-loop level and its properties are studied. Additionally, the obtained R-ratio is juxtaposed with its commonly employed approximation and the impact of the omitted higher-order π 2 -terms on the latter is discussed. In the Conclusions (Sect. 7) the basic results are summarized.
2 R-ratio of electron-positron annihilation into hadrons As noted above, the theoretical analysis of certain strong interaction processes relies on the hadronic vacuum polarization function Π (q 2 ), which is defined as the scalar part of the hadronic vacuum polarization tensor For the processes involving final state hadrons the function Π (q 2 ) (1) has the only cut along the positive semiaxis of real q 2 starting at the hadronic production threshold q 2 ≥ 4m 2 π (the discussion of this issue can be found in, e.g., Ref. [7]). In particular, the Feynman amplitude of the respective process vanishes for the energies below the threshold, that expresses the physical fact that the production of the final state hadrons is kinematically forbidden for q 2 < 4m 2 π . In turn, the known location of the cut of function Π (q 2 ) in the complex q 2 plane enables one to write down the pertinent dispersion relation with the once-subtracted Cauchy integral formula being employed. In Eq. (2) ∆ Π (q 2 , q 2 0 ) = Π (q 2 ) − Π (q 2 0 ), whereas R(s) stands for the discontinuity of the hadronic vacuum polarization function across the physical cut  Fig. 1 The integration contour in Eq. (6). The physical cut ζ ≥ 4m 2 π of the Adler function D(−ζ ) (4) is shown along the positive semiaxis of real ζ . This function is commonly identified with the so-called R-ratio of electron-positron annihilation into hadrons R(s) = σ (e + e − → hadrons; s)/σ (e + e − → µ + µ − ; s), with s = q 2 > 0 being the timelike kinematic variable, namely, the center-of-mass energy squared.
In practice one deals with the Adler function D(Q 2 ) [8], which is defined as the logarithmic derivative of the hadronic vacuum polarization function (1) with Q 2 = −q 2 > 0 being the spacelike kinematic variable. Note that the subtraction point q 2 0 entering Eq. (2) does not appear in Eqs. (3) and (4). The widely employed dispersion relation for the Adler function follows immediately from Eqs. (2) and (4), specifically [8] In particular, this dispersion relation enables one to extract the experimental prediction for the Adler function by making use of the corresponding experimental data on electron-positron annihilation into hadrons. However, to obtain the theoretical expression for the R-ratio itself, the relation inverse to Eq. (5) is required. The latter can be obtained by integrating Eq. (4) in finite limits, that yields [9,10] Specifically, this equation relates the R-ratio to the theoretically calculable Adler function and provides a native way to properly account for the effects due to continuation of the spacelike perturbative results into the timelike domain. The integration contour in Eq. (6) lies in the region of analyticity of the integrand, see Fig. 1. Note also that the relation, which expresses the hadronic vacuum polarization function in terms of the Adler function can be obtained in a similar way, namely [11] ∆ where Q 2 and Q 2 0 stand for the spacelike kinematic variable and the subtraction point, respectively. Basically, Eqs. (2)-(7) constitute the complete set of relations, which express the functions on hand in terms of each Table 1 Numerical values of the Adler function perturbative expansion coefficients d j (8). In the last column the numerical estimation of the five-loop coefficient d 5 [21] is listed. other. It is worth mentioning here that, in general, the pattern of applications of the dispersion relations in theoretical particle physics is quite diverse. For example, among the latter are such issues as the refinement of chiral perturbation theory [12,13], the accurate determination of parameters of resonances [14], the assessment of the hadronic light-by-light scattering [15], as well as many others.
It is necessary to outline that the derivation of dispersion relations (2)- (7) is based only on the kinematics of the process on hand and involves neither model-dependent phenomenological assumptions nor additional approximations. In turn, the relations (2)- (7) impose a number of strict physical inherently nonperturbative constraints on the functions Π (q 2 ), R(s), and D(Q 2 ), that should definitely be taken into account when one intends to go beyond the limits of applicability of the QCD perturbation theory. It is worthwhile to note that these nonperturbative restrictions have been merged with the corresponding perturbative input in the framework of dispersively improved perturbation theory (DPT) [16][17][18] (its preliminary formulation was discussed in Ref. [19]). In particular, the DPT enables one to overcome some intrinsic difficulties of the QCD perturbation theory and to extend its applicability range towards the infrared domain, see Ref. [16] and references therein for the details.
In the framework of perturbation theory the Adler function (4) takes the form of the power series in the so-called QCD couplant a ( ) where z = Q 2 /Λ 2 and b m n (the integer superscript m is not to be confused with respective power) stands for the combination of the β function perturbative expansion coefficients, specifically, b 0 The Adler function perturbative expansion coefficients d j were calculated up to the four-loop level (1 ≤ j ≤ 4), see Ref. [20] and references therein, whereas for the five-loop coefficient d 5 only numerical estimation [21] is available so far. The numerical values of the perturbative coefficients d j (8) are listed in Tab. 1. In turn, the β function perturbative expansion coefficients β j have been calculated up to the five-loop level (0 ≤ j ≤ 4), see Ref. [22] and references therein for the details.
In what follows the nonperturbative aspects of the strong interactions will be disregarded and a primary attention will be given to the theoretical description of the R-ratio of electron-positron annihilation into hadrons at moderate and high energies. For this purpose the effects due to the masses of the involved particles can be safely neglected (a discussion of the impact of such effects 1 on the low-energy behavior of the functions Π (q 2 ), R(s), and D(Q 2 ) can be found in, e.g., Refs. [16][17][18][23][24][25]). Additionally, for the scheme-dependent perturbative coefficients β j and d j the MS-scheme will be assumed and for the uncalculated yet five-loop coefficient d 5 its numerical estimation [21] will be employed. Thus, in the massless limit the relation (6) can be represented as (see also Ref. [26]) where stands for the spectral function and d ( ) (Q 2 ) denotes the -loop strong correction to the Adler function. As mentioned above, only perturbative contributions will be retained in Eq. (11) hereinafter, that makes Eq. (10) identical to that of both the foregoing DPT [16][17][18] and the so-called analytic perturbation theory 2 (APT) [26,27].
It has to be noted that, in general, the perturbative spectral function at small values of its argument may be altered by the terms of an intrinsically nonperturbative nature. For instance, the nonperturbative models discussed in Refs. [40][41][42][43] constitute a superposition of the perturbative spectral function with the so-called "flat" terms (which by definition do not affect the corresponding perturbative results at high energies), whereas the models [44] modify the low-energy behavior of the perturbative spectral function proceeding from certain phenomenological assumptions.
It is worthwhile to mention also that a "naive" approach to continue the spacelike perturbative result (8) into the timelike domain consists in merely identifying the timelike kinematic variable (s = q 2 ) with the spacelike However, as thoroughly discussed in Ref. [45], this prescription yields a misleading result, which differs from the proper one (10) even in the deep ultraviolet asymptotic, see also Refs. [21,46,47] as well as [16] and references therein.

R-ratio at the one-loop level
Let us address now the R-ratio of electron-positron annihilation into hadrons at the one-loop level. As discussed in the previous Section, for this purpose the spectral function ρ(σ ) (11), which enters the pertinent integral representation (10), is required. Since the involved strong correction to the Adler function (8) takes a simple form at the one-loop level and the calculation of ρ (1) (σ ) (11) appears to be quite straightforward. Specifically, the one-loop spectral function (11) for the positive values of its argument reads that, in turn, can be represented as  1 (σ ) ≤ 1/π 2 and decreases as 1/y 2 in both ultraviolet (y → ∞) and infrared (y → −∞) asymptotics.
Then, the corresponding one-loop strong correction to the R-ratio can also be easily obtained in an explicit form. Specifically, the integration (10) of the one-loop spectral function (16) where w = s/Λ 2 . The function a TL (s) (17) constitutes the one-loop couplant, which properly accounts for the effects due to continuation of the spacelike perturbative expression (13) into the timelike domain. It is worthwhile to note here that Eq. (17) has first appeared in Ref. [48] and only afterwards was derived in Refs. [9,11,26].  s (Q 2 ) (13) into the timelike domain (12). As one can infer from this Figure, at high energies the two couplants approach each other. At the same time, at moderate energies the deviation between the functions on hand becomes significant, whereas in the infrared domain their behavior turns out to be qualitatively different. Specifically, the function a (1) s (|s|) (13) diverges at low energies due to the infrared unphysical singularities, whereas a (1) (17) is a smooth monotone decreasing function of its argument, which contains no singularities for s > 0.
A somewhat simpler but approximate form of the strong correction to the R-ratio can be obtained by its further re-expansion. Specifically, for this purpose one splits the entire energy range 0 < s < ∞ into three intervals (namely, ln w < −π, −π < ln w < π, and ln w > π) and applies the Taylor expansion to r(s) in each of those intervals. At the one-loop level the implementation of these steps for the expression (17) yields where w = s/Λ 2 . In particular, as one can infer from Fig. 3 B, the re-expansions (18)- (20) may provide an accurate approximation of the function (17) in the aforementioned energy intervals, if the number of retained terms is large enough. However, as one can note, the convergence of the re-expansions (18)  In fact, there is another equivalent way to obtain the re-expansion of the strong correction to the R-ratio at high energies. Specifically, instead of following the lines described above, one can expand the corresponding spectral function ρ ( ) (σ ) (11) and then perform the integration in Eq. (10). In particular, at the one-loop level the Taylor expansion ofρ that, after its integration in Eq. (10), yields the result identical to Eq. (20). It is the latter prescription that will be employed in the next Section for the derivation of an approximate expression for the R-ratio at high energies at an arbitrary loop level.
4 Re-expansion of the R-ratio at high energies: π 2 -terms As outlined in Sect. 2, the strong correction to the R-ratio of electron-positron annihilation into hadrons (10) can be represented as where y = ln(σ /Λ 2 ) and ρ y (y) = ρ Λ 2 exp(y) denotes the corresponding spectral function (11) ρ with Eq. (14) being employed. In Eq. (23) d y (y) = d Λ 2 exp(y) stands for the strong correction to the Adler function being expressed in terms of y = ln(σ /Λ 2 ). Applying to the latter the Taylor expansion one can approximate 4 the spectral function (23) by Therefore, for √ s/Λ > exp(π/2) 4.81 the strong correction to the R-ratio (22) acquires the following form In particular, this equation implies that the strong correction to the R-ratio, being re-expanded at high energies, reproduces the naive continuation of the Adler function into the timelike domain [the first term on the righthand side of Eq. (26)] and additionally produces an infinite number of the so-called π 2 -terms.
Then, at the -loop level the perturbative expression for the strong correction to the Adler function reads (8) d ( ) where a ( ) (y) is the -loop perturbative couplant being expressed in terms of y = ln(σ /Λ 2 ). Since the latter satisfies the renormalization group equation the n-th derivative of the j-th power of the -loop couplant takes the form d n dy n a ( ) Thus, at high energies the -loop strong correction to the R-ratio (22) can be approximated by (see also Ref. [49]) The obtained re-expansion of the strong correction to the R-ratio (30) constitutes the sum of naive continuation of the strong correction to the Adler function into the timelike domain (12) and an infinite number of the π 2 -terms. Equation (30) explicitly proves the fact that at any given loop level the re-expansion of the strong correction to the R-ratio at high energies can be reduced to the form of power series in the naive continuation of the perturbative couplant into the timelike domain a ( ) s (|s|). As one can also note, in the re-expansion (30) the coefficients d j corresponding to various orders of perturbation theory turn out to be all mixed up, i.e., the -loop contribution to Eq. (10) appears to be re-distributed over the higher-order terms.
As discussed earlier, if the number of terms retained in Eq. (30) is large enough, then it can provide a rather accurate approximation of the strong correction to the R-ratio (10) for √ s/Λ > exp(π/2) 4.81. However, one usually truncates the re-expansion (30) at the order , thereby neglecting all the higher-order π 2 -terms 5 , that results in the following expression commonly employed in practical applications: where In this equation d j stand for the Adler function perturbative expansion coefficients (8), whereas δ j embody the contributions of the relevant π 2 -terms.
Equation (30) implies that the π 2 -terms do not appear in the first and second orders of perturbation theory, namely, that makes R ( ) appr (s) (31) identical to the naive expression R ( ) naive (s) (12) for = 1 and = 2. Starting from the third order of perturbation theory (i.e., for ≥ 3) the coefficients δ j (32) are no longer vanishing and constitute a combination of the pertinent perturbative expansion coefficients d j and β j of the first ( −2) orders. Specifically, the third-order and fourth-order coefficients read [21,46,47] and respectively. In these equations d j denote the Adler function perturbative expansion coefficients (8), whereas B j = β j /β j+1 0 stands for the combination of perturbative coefficients of the renormalization group β function.
In turn, at the fifth and sixth orders the coefficients δ j (32) can be represented as [16,21,49] The seventh-and eighth-order coefficients δ j (32) have recently been calculated as well, specifically [16,49]  Numerical values of the coefficients r j of the re-expanded approximate R-ratio (31). The last column employs the numerical estimation of the Adler function perturbative expansion coefficient d 5 [21].  Table 4 The relative weight (1 + |d j /δ j |) −1 ×100% of the π 2 -terms in the coefficients r j of the re-expanded approximate R-ratio (31).
The last column employs the numerical estimation of the Adler function perturbative expansion coefficient d 5 [21].
The explicit expressions for the coefficients δ j (32) at the higher orders can be found in App. C of Ref. [16].  (32), see also Tabs. 3 and 4. As will be discussed below, eventually this results in an essential distortion of the re-expanded approximation R ( ) appr (s) (31) with respect to the naive expression R ( ) naive (s) (12). In particular, the higher-order terms of R ( ) appr (s) (31) turn out to be substantially amplified and even sign-reversed with respect to those of R ( ) naive (s) (12), see Tabs. 1 and 3. It is worthwhile to note also that, as one can infer from Tabs. 2 and 4, the values of coefficients δ j rapidly increase as the order j increases, that makes the loop convergence of R ( ) appr (s) (31) worse than that of both R ( ) (s) (10) and R ( ) naive (s) (12), see Sect. 6 for the details.

Spectral function at the higher-loop levels
It is certainly desirable to leave the truncated re-expanded approximation R ( ) appr (s) (31) aside and have the function R ( ) (s) (10) calculated in a straightforward way beyond the one-loop level. To achieve this objective, the corresponding explicit expression for the involved spectral function ρ ( ) (σ ) (11) is required. Despite the latter becomes rather cumbrous for ≥ 2, the following method enables one to calculate ρ ( ) (σ ) explicitly at an arbitrary 6 loop level.
Specifically, it proves to be convenient to express the spectral function ρ ( ) (σ ) (11) in terms of the so-called "partial" spectral functionsρ ( ) j (σ ) corresponding to the j-th power of the -loop perturbative couplant (9), namely whereρ ( ) To obtain the partial spectral functionρ ( ) j (σ ) (41) for any j ≥ 1 it appears to be enough to calculate only the real and imaginary parts of the -loop couplant a ( ) s (Q 2 ) at the edges of its cut: Here a with n m = n! m! (n − m)! (44) being the binomial coefficient. Then, it is convenient to isolate in Eq. (43) its real and imaginary parts, specifi- where and ( j mod n) denotes the remainder on division of j by n. Therefore, the partial spectral function (41) reads (see also Refs. [41,50])ρ ( ) with σ ≥ 0 and j ≥ 1 being assumed. In particular, the first five relations (47) acquire a compact form In turn, the -loop functions a ( ) Re (σ ) and a ( ) Im (σ ) entering Eq. (47) can also be explicitly calculated in a similar way. Specifically, as mentioned in Sect. 2, the perturbative QCD couplant a ( ) s (Q 2 ) can be represented as the double sum (9) comprised of the functions It is worthwhile to decompose the functionā m n (Q 2 ) (53) at the edges of its cut into the real and imaginary parts, namely where u m n (σ ) and v m n (σ ) are the real functions of their arguments and σ ≥ 0 is assumed. Therefore, the functions a On the left-hand side of Eq. (53) and in Eqs. (54)-(56) the integer superscripts m are not to be confused with respective powers, whereas the coefficients b m n have been specified in Eq. (9). Then, to calculate the functions u m n (σ ) and v m n (σ ) entering Eqs. (55) and (56), it is convenient to split the left-hand side of Eq. (54) into two factors: The first factor on the right-hand side of Eq. (57) reads Proceeding along the same lines as earlier, one can cast the numerator on the right-hand side of this equation      with K(n) being specified in Eq. (46). Hence, the functions u m n (σ ) and v m n (σ ) (54) for m = 0 read with n ≥ 1 being assumed.
In turn, for m ≥ 1 the second factor on the right-hand side of Eq. (57) takes the following form where y = ln(σ /Λ 2 ). Since for real a and b ln(a ± ib) = ln a 2 + b 2 ± iπ where Following the same steps as above, one can cast Eq. (64) to where K(m) is defined in Eq. (46), and m ≥ 1 is assumed.
Therefore, the functions u m n (σ ) and v m n (σ ) (54) read and v m n (σ ) = with n ≥ 1 being assumed. Thus, the explicit expression for the -loop spectral function ρ ( ) (σ ) (40) takes the following form: where the functions u m n (σ ) and v m n (σ ) are specified in Eqs. (69) and (70), respectively. The higher-loop partial spectral functionsρ ( ) j (σ ) (47), which correspond to the j-th power (1 ≤ j ≤ ) of the -loop (2 ≤ ≤ 5) perturbative QCD couplant a implies that at the five-loop level the maximum value of the fifth-order functionρ (5) 5 (σ ) is about three orders of magnitude less than the maximum value of the first-order functionρ (5) 1 (σ ). In turn, as it will be discussed in the next Section, the fact that the functionρ ( ) j+1 (σ ) is subdominant toρ ( ) j (σ ) eventually results in an enhanced higher-loop stability of the proper expression for the R-ratio (10) at moderate and low energies with respect to both its naive form (12) and the commonly employed truncated re-expanded approximation (31).
The plots of the spectral function ρ ( ) (σ ) (71) at the first five loop levels (1 ≤ ≤ 5) are displayed in and Additionally, Fig. 5 A implies that the range of y, where the difference between ρ ( ) (σ ) and ρ ( +1) (σ ) is sizable, is located in the vicinity of y = 0 and becomes smaller at larger . This issue is also elucidated by Fig. 5 B, which shows the relative difference between the -loop and ( + 1)-loop spectral functions (71) at various loop levels. In this equation stands for the j-th order -loop "timelike" effective couplant, that constitutes the proper continuation of the j-th power of -loop QCD couplant a ( ) s (Q 2 ) j into the timelike domain.
As discussed in Sect. 3, at the one-loop level the first-order function (76) acquires a quite simple form, namely (17) A TL,1 (s) = a (1) Despite the fact that the partial spectral functionsρ ( ) j (σ ) (47) become rather cumbrous at the higher loop levels, it appears that the integration in Eq. (76) can be performed explicitly for > 1, too. For example, the two-loop first-order function (76) reads [9] A (2) TL,1 (s) = a (1) whereas the second-order one can be represented as In these equations a TL (s) is given by Eq. (77), B j = β j /β j+1 0 , and W (s) = ln ln 2 w + π 2 , w = s In turn, at the three-loop level the first-order function (76) takes the form where A TL,1 (s) is specified in Eq. (78) and TL (s) +W (s) ln w.
As for the four-loop first-order function (76), it can be represented as TL,1 (s) + 1 ln 2 w + π 2 3 where the functions A  (83), respectively. It is worthwhile to note also that the functions A ( ) TL, j (s) (76) entering Eq. (75) can be computed numerically by making use of the routines included in the freely available program packages [50] [which, being based on a less universal method of calculation of the relevant spectral function (11) than that of Eq. (71), is applicable at first four loop levels only] and [51].
As discussed earlier, the π 2 -terms play a valuable role in the studies of the strong interaction processes in the timelike domain, and their ignorance (complete or partial) may yield misleading results. In particular, this issue is illustrated by Fig. 6, which displays the two-loop timelike effective expansion function of the second order A (2) TL,2 (s) [Eq. (79), solid curve] and the approximations (dashed and dot-dashed curves) corresponding to various orders of its re-expansion for √ s/Λ > exp(π/2) 4.81: 10 ln 2 ln w − 9 ln ln w + 1 − π 2 − − π 4 ln 7 w 3B 1 10 20 Specifically, the dashed curve shows the result of naive continuation of the respective term of the Adler function perturbative expansion into the timelike domain (12) [the first term on the right-hand side of Eq. (85)], whereas the dot-dashed curves additionally include 7 the lowest-order π 2 -terms [numerical labels indicate the highest absolute value of power of ln w retained in Eq. (85)]. Figure 6 implies that the re-expansion (85) converges rather slowly at low and moderate energies. Furthermore, even at relatively high energies A TL,2 (s) (79) considerably differs from a (2) s (|s|) 2 , the latter being the only part of the function (79), which is retained in the two-loop approximation of the R-ratio (31). For example, as one can infer from Fig. 6, for √ s/Λ = 20 the function a (2) s (|s|) 2 exceeds A TL,2 (s) by about 21 %, and to securely achieve 10 % accuracy in the re-expansion (85) the inclusion of the π 2 -terms up to the order of ln −7 w is required.
The issue of the higher-loop stability of the R-ratio of electron-positron annihilation into hadrons is also elucidated by Fig. 7. In particular, this Figure (31) is worse than that of both the proper expression R ( ) (s) (10) and the naive one R ( ) naive (s) (12). As discussed earlier, this is primarily caused by the fact that the convergence range of the re-expanded approximation (31) is strictly limited to √ s/Λ > exp(π/2) 4.81, that, in turn, results in rather large values of the higher-order coefficients δ j (32) embodying the contributions of the corresponding π 2 -terms (30). In particular, as one can infer from Fig. 7, beyond the two-loop level (i.e., for ≥ 3) the curves corresponding to R ( ) (s) (10) are nearly indistinguishable from each other, whereas the curves corresponding to R ( ) appr (s) (31) start to swerve quite above the boundary of its convergence range. For example, at √ s/Λ = 2 exp(π/2) 9.62 the relative difference between the -loop and ( + 1)-loop strong corrections to the proper expression R ( ) (s) (10)  To expound the accuracy of approximation of the R-ratio (10) by its truncated re-expansion R ( ) appr (s) (31), it is worthwhile to mention also the following. At a moderate energy scale of the mass of τ lepton √ s = M τ the relative difference between the -loop strong corrections r ( ) (s) (10) and r ( ) appr (s) (31) turns out to be as high as 26 %, 28 %, 14 %, 2 %, and 7 % at the one-, two-, three-, four-, and five-loop levels, respectively. Moreover, it appears that even at high energies the ignorance of the higher-order π 2 -terms in the truncated re-expanded approximation (31) may produce a considerable effect. In particular, this issue is illustrated by Fig. 8, which displays the quantity at various loop levels. Specifically, as one can infer from Fig. 8, in the energy range planned for the future CLIC experiment [4] the effect of inclusion of the π 2 -terms discarded in the approximate expression R ( ) appr (s) (31) is either comparable to or prevailing over the effect of inclusion of the next-order perturbative correction.

Conclusions
The strong corrections to the R-ratio of electron-positron annihilation into hadrons are studied at the higherloop levels. In particular, the derivation of a general form of the commonly employed approximate expression for the R-ratio (which constitutes its truncated re-expansion at high energies) is delineated, the appearance of the pertinent π 2 -terms is expounded, and their basic features are examined. It is demonstrated that the validity range of such approximation is strictly limited to √ s/Λ > exp(π/2) 4.81 and that it converges rather slowly when the energy scale approaches this value. The spectral function required for the proper calculation of the R-ratio is explicitly derived and its properties at the higher-loop levels are studied. The developed method of calculation of the spectral function enables one to obtain the explicit expression for the latter at an arbitrary loop level. By making use of the derived spectral function the proper expression for the R-ratio is calculated up to the five-loop level and its properties are examined. In particular, it is shown that the loop convergence of the proper expression for the R-ratio is better than that of its commonly employed approximation and that the omitted higher-order π 2 -terms in the latter may produce a considerable effect.