Using improved Operator Product Expansion in Borel-Laplace Sum Rules with ALEPH $\tau$ decay data, and determination of pQCD coupling

We use improved truncated Operator Product Expansion (OPE) for the Adler function, involving two types of terms with dimension $D=6$, in the double-pinched Borel-Laplace Sum Rules and Finite Energy Sum Rules for the V+A channel strangeless semihadronic $\tau$ decays. The generation of the higher order perturbative QCD terms of the $D=0$ part of the Adler function is carried out using a renormalon-motivated ansatz incorporating the leading UV renormalon and the first two leading IR renormalons. The trunacted $D=0$ part of the Sum Rules is evaluated by two variants of the fixed-order perturbation theory (FO), by Principal Value of the Borel resummation (PV), and by contour-improved perturbation theory (CI). For the experimental V+A channel spectral function we use the ALEPH $\tau$-decay data. We point out that the truncated FO and PV evaluation methods account correctly for the renormalon structure of the Sum Rules, while this is not the case for the truncated CI evaluation. We extract the value of the ${\overline {\rm MS}}$ coupling $\alpha_s(m_{\tau}^2) = 0.3235^{+0.0138}_{-0.0126}$ [$\alpha_s(M_Z^2)=0.1191 \pm 0.0016$] for the average of the two FO methods and the PV method, which we consider as our main result. If we included in the average also CI extraction, the value would be $\alpha_s(m_{\tau}^2) = 0.3299^{+0.0232}_{-0.0225}$ [$\alpha_s(M_Z^2)=0.1199^{+0.0026}_{-0.0028}$]. This work is an extension and improvement of our previous work [Eur.Phys.J.C81 (2021) 10, 930] where we used for the truncated OPE a more naive (and widely used) form and where the extracted values for $\alpha_s(M_Z^2)$ were somewhat lower.


I. INTRODUCTION
One of the most important challenges of physics is the determination of fundamental parameters. In the case of strong interactions, the main parameter is the strong running coupling a(Q 2 ) ≡ α s (Q 2 )/π which depends on the squared (spacelike) momentum scale Q 2 (≡ −q 2 ) characteristic of the process. This coupling is well determined at high energies 1 GeV 2 < Q 2 ( M 2 Z ) because of the high precision of the experiments and the theory. On the other hand, determinations at lower energies are a good test of the consistency of the theory. We study here the quantities related with the semihadronic decay width of τ lepton, which is at low scales Q m τ (∼ 1 GeV) and for which high precision experimental results are available from the ALPEH collaboration [2][3][4][5]. Both the experimental and the theoretical results are related with the two-point correlation function of (ud) quark currents Π(Q 2 ); the experimental results are given in terms of the spectral function ω(σ) ∝ ImΠ(−σ − i ), while the theoretical results are usually given in terms of contour integrals of the derivative of the correlation function, D(Q 2 ) ∝ dΠ(Q 2 )/d ln Q 2 , known as the Adler function.
Specific contour integrals involving the Adler function give us various τ -decay sum rules, including the inclusive strangeless decay ratio R τ = Γ(τ → ν τ hadrons)/Γ(τ → ν τ e −ν e ) 1 . The Adler function was calculated by perturbation techniques of QCD up to O(α 4 s ), cf. [12][13][14]. Determinations of the coupling from such sum rules, i.e., at low momenta ∼ m τ , allow us to see the reliability of the theory because we can compare the extracted value of α s with the known determinations at high energies through renormalisation group equation (RGE) evolution (see for instance [15][16][17]).
(dimension D = 0), and the nonperturbative corrections (D > 0), and the latter are often small such as in the case of R τ (r τ ) [9,18]. One of the novel aspects covered in the present work, in comparison with the previous one [1], is that we will deal with the nonperturbative contributions more carefully. In particular here the D (≡ 2k) ≥ 6 condensate contributions to the Adler function d(Q 2 ) have two parts, which eliminate the renormalon ambiguity originating from the two infrared (IR) renormalon contributions of the Adler function at u = k (k = 3, 4, . . .), one of the type ∼ 1/(k − u) γ k +1 and the other ∼ 1/(k − u) γ k (where γ k = 1 + kβ 1 /β 2 0 is unity in the large-β 0 limit).
The D > 6 condensates will not be included in our analysis, as we use an extended Adler function that is based on a renormalon-motivated ansatz that includes only the first two IR renormalons u = 2, 3 (and the first UV renormalon u = −1) [19]. Evaluation of the perturbative part (D = 0) of sum rules is based on (re)summations and truncations. These (re)summations require integration of the perturbation series on a circular contour in the complex Q 2 -plane with radius |Q 2 | = σ max (∼ m 2 τ ), which involve the QCD running coupling parameter a(Q 2 ) along this circle. Due to truncations, different ways of evaluation give different results. We will apply four of them. The first two are the fixed-order (FO) approach and its variant FO, where the coupling in the Adler function is Taylor-expanded around the spacelike point Q 2 = σ max > 0, and in the FO the series in powers of a(σ max ) is performed and truncated, while in the FO the series in logarithmic derivatives of a(σ max ) is performed and truncated. The other is the contourimproved (CI) evaluation, which evolves a(Q 2 ) through RGE along the contour Q 2 = σ max e iφ . The third is the so called principal-value (PV) evaluation, where the inverse Borel transformation is appplied to the singular part of the Borel transform of the Adler function and the principal value prescription is applied for the integration across the IR renormalon singularities [20][21][22].
The extraction of α s (m 2 τ ) from the τ -decay data was performed in the past in the literature with the FO and CI approaches, giving different results. The tendency shows that the CI gives higher values of α s (m 2 τ ) than the FO, even if the duality violation effects are explicitly taken into account with a model [4,14,23,24]. The discrepancy between the CI and FO was discussed in [25,26] in the large-β 0 (LO) approximation and then beyond-LO (bLO) in a renormalon-motivated model. The authors of [25,26] argued that the truncated FO approach in the r τ sum rule takes correctly into account some cancellations of the leading renormalon contributions of the Adler function. They presented theoretical arguments for this cancellations in the LO approximation. Further, they argued that such cancellation do not take place in the truncated CI approach. Beyond LO (bLO), such cancellations were demonstrated in an approach with a modified Borel transform in a specific renormalisation scheme [27]. In our previous work [1], we presented arguments that such cancellations take place in sum rules at the bLO level in any renormalisation scheme when the D = 0 contribution to the sum rule is written in terms of the series in logarithmic derivatives of a(σ max ) (cf. Appendix of [1]).
The main goal of this work is to determine the value of the running coupling α s (m 2 τ ), by using a renormalonmotivated extension [19] of the Adler function for the D = 0 contribution and the OPE with the terms D = 4 and D = 6 Eq. (1). The sum rules we use for this determination are the Borel-Laplace sum rules, which are in the double-pinched form in order to suppress the duality-violating effects. 2 The D = 0 contribution is evaluated by the (truncated) methods FO, FO, PV and CI. The truncation index of these contributions is then determined by considering the double-pinched finite energy sum rules (FESRs) a (2,0) and a (2,1) and finding the value of the truncation index for which these quantities become locally stable under the variation of the index. This work can be regarded as a continuation and improvement of our previous analysis [1], where now the improved form (1) of the OPE terms is used; further, the variation of the higher order perturbation coefficients of the renormalon-motivated extension of the Adler function is performed now in a more realistic way. This paper is organized as follows. In Sec. II we present the improved D ≥ 6 terms of the OPE for the Adler function, and summarise the general form of the sum rules for the (strangeless) semihadronic τ -decay data. In Sec. III we present briefly the renormalon-motivated extension of the D = 0 Adler function, including the variation of the parameters of the model which reflect the variation of the higher order perturbation coefficients. In Sec. IV the specific sum rules used later in the analysis are presented. In Sec. V we summarise the various (truncated) methods used for the evaluation of the D = 0 contribution to the sum rules. Section VI is the main part of this work, it contains the numerical analysis of the extraction of the value of α s (m 2 τ ) and of the D = 4, 6 condensates, for the full V+A channel. In Sec. VII we summarise the results and present our conclusions. A summary of our results is given in [35].

II. SUM RULES AND ADLER FUNCTION
The Adler function D(Q 2 ) is defined as a logarithmic derivative of the quark current polarisation function Π(Q 2 ) We will consider the total (V+A)-channel, i.e., Π(Q 2 ) will be the total (V+A)-channel polarisation function In this V+A sum (3), the term Π gives negligible contribution (to the sum rules) because ImΠ where J µ are the up-down quark vector and axial vector currents, J µ =ūγ µ d andūγ µ γ 5 d for J = V, A, respectively. We recall that q 2 ≡ −Q 2 is the square of the momentum transfer, q 2 = (q 0 ) 2 − q 2 . The usually used theoretical expression of the polarisation function has the following OPE [36] form: The D = 2k term in this expansion is related to a renormalon singularity of the Borel transform B[d](u) of the Adler function d(Q 2 ) D=0 . Namely, if the perturbation expansion of d(Q 2 ) D=0 in powers of a(µ 2 ) ≡ α s (µ 2 )/π is 3 the expansion of the Borel transform B[d](u; κ) is In the leading-β 0 (LB) approximation (i.e., where β 1 = 0), we have 3 Here, κ ≡ µ 2 /Q 2 (0 < κ ∼ 1) is the dimensionless parameter for the renormalisation scale µ 2 . 4 The β j coefficients in our convention appear in the following form of the RGE for a(Q 2 ): .. The one-loop and two-loop coefficients β 0 and β 1 are universal (i.e., scheme independent) in mass independent schemes, β 0 = (11 − 2N f /3)/4 (= 9/4 for N f = 3) and β 1 = (1/16)(102 − 38N f /3). The other coefficients β j (j ≥ 2) depend on (and characterise) the renormalisation scheme. 5 In the general case, the power in ∼ 1/(k − u) γ k is γ k = 1 + kβ 1 /β 2 0 − γ is the effective leading-order anomalous dimension of the D-dimensional OPE operator O D (D = 2k), and the corresponding term in the OPE (6) would have in such a case, [19] [cf. also the later discussion in the second paragraph after Eq. (21)]. γ k = 1, i.e., the singularity ∼ 1/(k − u) γ k corresponds in the LB approximation to a simple single pole. However, the Borel transform of the D = 0 Adler function is known to all orders in the LB approximation [43][44][45], and it turns out that only the u = 2 singularity is a single pole, but all other IR singularities (at u = 3, 4, . . .) are combinations of double and single poles. Consequently, when going beyond the LB approximation, the renormalon poles of the Adler function at u = k ≥ 3 become ∼ 1/(k − u) γ k +1 and ∼ 1/(k − u) γ k (and lower singularities), if we assume that the beyond-LB effects come from the RGE-evolution (beyond one-loop) of the QCD coupling and that the leading-order anomalous dimensions of the corresponding operators remain of the LB-type. 6 The corresponding OPE terms which cancel the ambiguities from these singularities at u = k (k ≥ 3) are then ∼ 1/(Q 2 ) k /a(Q 2 ) and 1/(Q 2 ) k [19]. 7 This implies that the OPE expansion of the Adler function, improved with respect to the expansion (6), has the following form: which has two different condensates for each D = 2k ≥ 6. This form of the OPE terms was noted already in our previous work [1] [Eq. (59) there], but the implementation was left pending. In the expansion (9), we neglected terms O(a), i.e., terms ∼ O 2k a(Q 2 ), because they are relatively small for the considered momenta (|Q 2 | ≈ 3 GeV 2 ). It can be checked that this expansion then corresponds to the following expansion of the polarisation function: We will use this OPE expansion (9)-(10) in the sum rules. According to the general principles of Quantum Field Theory, the considered polarisation function Π(Q 2 ; µ 2 ) and its logarithmic derivative D(Q 2 ) are holomorphic (i.e., analytic) functions of Q 2 in the complex Q 2 -plane with the exception of the real negative axis (−∞, −m 2 π ). Then if g(Q 2 ) is a (arbitrary) holomorphic function of Q 2 , then the Cauchy theorem can be applied to the integral dQ 2 g(Q 2 )Π(Q 2 ; µ 2 ) along a closed path in the complex Q 2 -plane that consists, for example, of the circle of finite radius |Q 2 | = σ max (≡ σ m ) and lines above and below the negative axis avoiding thus the enclosure of the values Q 2 < 0 where the integrand is not holomorphic (cf. Fig. 1). This then implies In the sectors of the path C 1 we have ε → +0. 6 Cf. discussion in the next Section on this point. 7 This can also be formulated in the following way [19]: the corresponding dimension D ≡ 2k operators have the leading-order anomalous where the integration on the right-hand side of Eq. (11b) is counterclockwise (Q 2 = σ m e iφ , −π < φ < π), and ω(σ) is proportional to the discontinuity (spectral) function of the (V+A)-channel polarisation function The quantity ω(σ) was measured by OPAL Collaboration [40,46], and to an even higher precision by ALEPH Collaboration [2][3][4][5], in semihadronic strangeless decays of the τ lepton. In our analysis we will use the ALEPH data; they are presented in Fig. 2. Integration by parts allows us to replace the theoretical polarisation function in the sum where D th (Q 2 ) is given by the OPE expansion (9), and the (holomorphic) function G is the following integral of g: which is independent of the form of path from −σ m to Q 2 in the Q 2 -complex plane because g(Q 2 ) is holomorphic. The quark vector and axial vector current polarisation functions Π (k) J (Q 2 ) and the Adler function D(Q 2 ) are quantities which are holomorphic (analytic) functions of Q 2 in the complex Q 2 -plane with the exception of the negative semiaxis, i.e., they are spacelike quantities. On the other hand, the spectral function (12) and the sum rules (13) are timelike observables, being functions of the squared energy σ > 0 (= −Q 2 ) or σ m ≡ σ max > 0. Several other timelike quantities exist that have the form of integrals of D(Q 2 ) and are of phenomenological interest [47], such as: (a) e + e − → hadrons production ratio R(s) [48,49] (with any s > 0); (b) the leading order hadronic vacuum polarisation contribution to the muon anomalous magnetic moment (g µ − 2) [50,51], where the dominant squared timelike momenta are s ∼ m 2 µ (∼ 0.01 GeV 2 ) [52,53].

III. RENORMALON-MOTIVATED EXTENSION OF THE ADLER FUNCTION
Here we summarise the main results of the renormalon-motivated extension of the Adler function. We refer for details to [1,19]. The expansion of the D = 0 part of the Adler function in powers of a ≡ α s /π has the form (7), where the first four expansion coefficients (d 0 = 1; d 1 ; d 2 ; d 3 ) are exactly known [12][13][14]. If we reorganise this expansion into the expansion in terms of the (related) logarithmic derivatives we obtain By the use of the MS scheme RGE (which is known up to five-loops [54]) the logarithmic derivatives a n and the powers a k can be related and the expansion coefficients d n and d k are related analogously The coefficients k s (n + 1 − s) and k s (n + 1 − s) are specific κ-independent combinations of the β-function coefficients c j = β j /β 0 , cf. [19].
which contains all the information about the Adler function coefficients d n (and thus d n ) but, in contrast to B[d](u), has the simple one-loop type renormalisation scale dependence As a consequence, the renormalon structure of B[ d](u) has no explicit effects coming from the beyond one-loop RGE running of the coupling, i.e., no terms ∝ β 1 in the power indices of the singularities [cf. Eq. (17)]. We adopt the approximation that the renormalon structure of B[ d](u) has the form as obtained in the LB approximation [43,44]: . . near the UV renormalon locations u = −1, −2, . . . As shown in [19], this then We point out that the general structure of the IR singularity at O D /β 0 (using the notation of [19]) where γ O6 /β 0 = −1 or 0 only in the LB approximation, but becomes a set of noninteger numbers in the exact case (beyond-LB) [55], namely γ (1) O6 /β 0 ≈ 0.197; 0.247; . . ., corresponding to a set of nine operators with dimension D = 6 (cf. also [56][57][58]). The corresponding D = 6 OPE terms of the Adler function would then be nine terms ∼ a(Q 2 ) γ (1) O 6 /β0 /(Q 2 ) 3 , instead of the two D = 6 terms appearing in Eq. (9). We will not pursue this extension in this work, but wish to point out that this indicates the challenging way how one may extend the present renormalon-motivated model to structures which contain improved beyond-LB effects in the IR renormalons u = p ≥ 3 and the corresponding OPE terms of dimension D = 2p ≥ 6.
Thus the ansatz we make for B[ d](u) (with κ = 1) in the MS scheme includes the singularities at the locations u = 2, 3 and u = −1 [19] where the value of the parameter α in the MS scheme is fixed, α = −0.255 [19], by the knowledge of the subleading coefficientĉ (D=4) 1 coming from the D = 4 subleading-order anomalous dimension and Wilson coefficient [45,59]. The other five parameters ( K and the residues d IR 2,1 , d IR 3,2 , d IR 3,1 , d UV 1,2 ) are determined by the knowledge of the first five coefficients d n (and thus d n ), n = 0, 1, 2, 3, 4. However, for n = 4, the coefficient d 4 of the Adler function is not yet exactly known. There are, however, estimates for the value of this coefficient. In Ref. [19], a similar ansatz as (22) was used in the lattice-related MiniMOM renormalisation scheme, not containing d IR 3,1 /(3 − u) term, and the resulting estimate d 4 = 338.19 (in MS scheme, with κ = 1) was extracted upon reexpansion. On the other hand, the effective charge (ECH) method [60] gives an estimate d 4 = 275 [14,61]; a similar estimate d 4 = 277 ± 51 was obtained in [62] based on Padé approximants; d 4 = 283 was obtained in [25] by extrapolating the approximately geometric series behaviour of FOPT of r In Table I we present the resulting parameters K and d IR i,j (i = 2, 3; j = 1, 2) and d UV 1,2 for the central and the border values of d 4 Eq. (23). In Table II, where γ p = 1 + pβ 1 /β 2 0 ; γ p = 1 − pβ 1 /β 2 0 , and the residues in the expression (24) and a linear polynomial of u. That model was later referred to and compared with Padé-related approaches in Ref. [62].

IV. SPECIFIC SUM RULES USED
Once having the expansion coefficients of d(Q 2 ) D=0 , we can apply various sum rules, i.e., various weight functions g(Q 2 ) [cf. Eqs. (13)-(14 We will primarily consider the double-pinched Borel-Laplace transforms B(M 2 ) where M 2 is a complex squared energy parameter 9 The general sum rule Eq. (13) has on the left-hand side the experimental value, and on the right-hand side the theoretical value. Using the (double-pinched) Borel-Laplace weight function g(Q 2 ) Eq. (26a), the left-hand side of the sum rule is written as and the right-hand side as The last terms here are the contributions of the dimension D = 2k condensates of the OPE of the Adler function (9) where k = 3, 4, . . .. The coupling a(Q 2 ) at the OPE term O 2k /a(Q 2 ) will be taken to run according to one-loop RGE along the contour and 10 J s (A) are the integrals The expressions for J s (A) in terms of sums, for any complex A and any positive integer s, are given in the Appendix. Further, explicit closed expressions are given there for J s (A) with s = 1, 2, 3 (which are relevant for the D = 6 contribution). On the other hand, one can use FESRs with (double-pinched) momenta a (2,n) which are associated with the following weight functions g (2,n) (n = 0, 1, . . .): The experimental and the theoretical parts of these FESR moments are then (we subtract unity for convenience) We will consider in particular the first two moments a (2,0) and a (2,1) , up to D = 6 terms k≥2 a (2,0) The expressions for a (2,n) th (σ m ) D=2k for general integer n ≥ 0 (and k ≥ 2) are given in the Appendix. We note that a (2,1) (σ = m 2 τ ) D=0 is the D = 0 part of the canonical QCD (and strangeless and massless) τ -decay ratio r which appears in the contour integrals in the sum rules, Eqs. (28) and (33b), is written as truncated Taylor expansion in powers of a(κσ m ) up to (and including) a(κσ m ) Nt . We point out that the Adler function d(Q 2 ) D=0 is a spacelike quantity (defined for general complex Q 2 ); the sum rules are timelike quantities (defined only for positive quantities σ = σ m > 0), but are written in the FO approach in terms of powers of a(κσ m ) where Q 2 = κσ m > 0 is a spacelike point in the complex Q 2 -plane. 11 The expression (35) has some dependence on the renormalisation scale parameter κ due to truncation, (d/d ln κ)d(Q 2 ; κ) D=0,pt ∼ a N t +1 .  (36) in the countour integrals is written as truncated Taylor expansion in logarithmic derivatives a k+1 (κσ m ) up to (and including) a Nt (κσ m ). D=0 are thus practically free of renormalon growth when n increases. We refer for additional explanation to [1] (Sec. IV.B there).

VI. RESULTS OF FITTING THE BOREL-LAPLACE SUM RULES
In this Section we first fit the theoretical double-pinched Borel-Laplace sum rules, cf. Eqs. (13) and (27)- (29), to the ALEPH experimental data (V+A channel) as explained in Sec. II. The theoretical Borel-Laplace sum rules B th (M 2 ; σ m ) are evaluated with various evaluation methods and various truncation indices N t in the D = 0 contribution, cf. Sec. V for explanations. Subsequently, the resulting predictions for the (double-pinched) FESR momenta a (2,0) (σ m ) and a (2,1) (σ m ) [cf. Eqs. (33)- (34)] are compared with the experimental data, for various truncation indices N t , and the optimal N t is fixed where the relative stability of these FESRs under variation of N t is reached. We point out that in the analysis, the higher order contributions of the D = 0 Adler function contributions are generated (estimated) by the renormalon-motivated ansatz mentioned in Sec. III. Further, throughout the analysis, the OPE expansion (9) of the Adler function is performed up to D(≡ 2k) = 6 terms. Going beyond D = 6 terms is not well motivated in the present analysis, as the D = 0 Adler function is generated by a renormalon-motivated ansatz which includes IR renormalons up to u = +3 and not beyond. This means that an assumption is made that the higher IR renormalons (u = +4, etc.) give suppressed contributions to the D = 0 Adler function; such an assumption would then also suggest that the D ≥ 8 OPE contributions to the Adler function are in general suppressed.
[1] for more explanation. is given in units of 10 −4 GeV 6 because it is in general smaller than O . The minimisation of χ 2 Eq. (39) is performed by varying these four parameters simultaneously. In most cases the fits are very good, namely χ 2 10 −3 , , cf. Fig. 3.  In Table III we present the results of this analysis. 13 The uncertainties in the Table were obtained by combining various theoretical uncertainties and the experimental uncertainty, as will be explained below in more detail for the case of the parameter α s (m 2 τ ). 13 From the above values of O 4 V +A , the corresponding values for the gluon condensate are obtained by using the relation aGG = 6 O 4 + 6f 2 π m 2 π , where 6f 2 π m 2 π ≈ 0.00199 GeV 4 .
The extracted values for α s , with uncertainties from various sources given separately, are The central values were extracted for the truncation index N t = 8, 5, 8, 4 for the methods FO, FO, PV and CI, respectively, cf. Table III, 14 and these values of N t were obtained by looking for the local stability of the resulting FESRs a (2,j) (σ m ) (j = 0, 1) under variation of N t as explained earlier (see also Figs. 4 later below). In the results (40), at the symbol '(κ)' is the variation when the renormalisation scale parameter κ (≡ µ 2 /Q 2 ) is vared from κ = 1 up to κ max = 2 and down to κ min = 2/3, respectively. 15 At the symbol '(N t )' is the maximal variation when the truncation number is varied around its central value N t : in the cases of FO, FO, and PV the variation was in the interval N t = 8 +2 −3 , 5 ± 2, 8 +2 −3 , respectively (thus the case N t = 5 that is independent of the renormalon model is always included); in the case of CI, N t = 4 +2 −1 , i.e., N t = 2 was not considered, for being an extreme truncation). At the symbol '(d 4 )' is the variation when the coefficient d 4 varies according to Eq. (23) [cf. also Table I]. At the symbol '(exp)', the variations are (rough) estimates of the experimental uncertainties, and were obtained by the method explained in [1]. 16 In the PV method, there is an additional (fourth) source of theoretical uncertainty '(amb)', which is an estimate of uncertainty due to the Borel integration ambiguity for the Adler function [1].
In the results Eqs. (40) and Table III, the total uncertainties were obtained by adding the mentioned various uncertainties in quadrature. We see from Eqs.   (40). We can see that in general, and for a reasonably wide range of N t , the resulting predicted values of a (2,0) (σ m ) and a (2,1) (σ m ) are well compatible with the experimental results. One may ask how the momenta a (2,0) (σ m ) and a (2,1) (σ m ) behave under variation of N t when the parameters α s , and O (2) 6 are not varied but kept fixed. In that case, only the D = 0 contribution varies with N t , while D = 4, 6 contributions become N t -independent. In Figs. 5(a)-(b) we present these results, using for α s the central value of each method, i.e., the corresponding central values in Table III. We can see that, in contrast to the results of  At the end of this Section, we wish to address the following question. In the considered MS scheme we know the β-function (17) only up to the five-loop level (i.e., up to the coefficient β 4 [54]), but we nonetheless use in the FO and PV approaches the relatively high truncation index N t = 8. In principle, due to the connections (19b), the Adler function perturbation coefficients d n generated from the perturbation coefficients d n−s of B[ d](u) involve the coefficients k s (n + 1 − s) (s = 1, . . . , n − 1) which in turn are functions of β j (j = 1, . . . , n − 1) [19]. This implies that d 6 involves β 5 and d 7 (at a Nt = a 8 ) involves β 5 and β 6 coefficients which are taken to be zero here. Such an effect of truncation of the β-function, i.e, setting β 5 = β 6 = 0, in principle changes the coefficients d n (n ≥ 6) deduced from the Borel transform B[ d](u) Eq. (22) of the renormalon-motivated model. In [19], these effects were investigated in the case of the so-called Lambert MiniMOM (LMM) scheme, and it was shown there that the truncation of the LMM scheme at four-loop level (β j = 0 for j ≥ 4) gives numerically very similar results for the relevant parameters C (D) j,k and ratios of residues d X p,k / d X p,k , leading to similar expressions for the Borel transform B[d](u) of the Adler function in the two cases and thus to similar values of d n (n ≥ 5). 17 For these reasons, we believe that in the renormalon-motivated model considered here (in the MS scheme), the values of the coefficients d n in FO and δd n in PV approach (n = 6, 7) are not affected significantly by the setting β 5 = β 6 = 0.

VII. FINAL RESULTS AND CONCLUSIONS
We applied double-pinched Borel-Laplace sum rules to the V+A channel semihadronic strangeless τ decay data of ALEPH. We used for D = 0 contribution to the Adler function the renormalon-motivated extension [19] for the higher order terms. We applied four different methods of evaluation (FO, FO, CI, PV). The optimal truncation index N t (of the D = 0 Adler contribution) in the Borel-Laplace sum rules was fixed in such a way that the predicted FESR momenta a (2,0) and a (2,1) showed local stability under the variation of N t . As argued in detail in our previous work, the FO, FO and PV methods of evaluation of the sum rules lead to truncated series 18 that have the leading renormalon contribution of the Adler function (double UV renormalon pole at u = +1) suppressed in them; and that the CI method in the sum rules does not have this property. 19 We refer for details on these points to [1], and especially Appendix A there.
The four different methods (FO, FO, CI, PV) give the main results presented, for each method, in Eqs. (40) and Table III. On the grounds argued above, we exclude the CI methods from our average, and the arithmetic average of the three methods (FO, FO and PV) then gives  Table IV we present, for comparison, the values of α s (m 2 τ ) which were extracted from the ALEPH τ -decay data by various groups, who used various sum rules and various methods of evaluation. We refer to [1] for discussion and comparison of various of the approaches and results presented in the Table. None of those works used OPE in the form (9) but rather in the form (6).
If, on the other hand, we included the CI method results in the average, the value of the coupling would significantly increase  18 In the PV case, the truncated series refers to the series δd(σme iφ ) in Eq. (37) which is free from the leading renormalon singularities and gives the corresponding contour integrals (sum rules) also free from the leading renormalon singularities. 19 These arguments, presented in a very general way in [1] (especially Appendix A there), are also compatible with somewhat related arguments presented in [25][26][27]. A construction of Borel transforms of CI sum rules was proposed and investigated in [64] (cf. also [65]).
Since our results have relatively high truncation index N t = 8 (where a Nt is the power where truncation is made) for the FOPT and PV approaches, one may wonder how much the obtained results depend on the specific renormalonmotivated Adler function. The FOPT is independent of the renormalon-motivated extension for N t = 5 [the coefficient d 4 is taken according to Eq. (23)], and this is true to a large degree also for the PV approach. As can be seen from the results (40a) and (40e), the extracted values of the coupling change only little if N t = 5 is taken instead of N t = 8: the central value of α s (m 2 τ ) changes from 0.3228 to 0.3171 in the FOPT case, and from 0.3269 to 0.3277 in the PV case; the average central value α s (m 2 τ ) = 0.3235 Eq. (41) changes to α s (m 2 τ ) = 0.3219 when N t = 5 in all three methods is taken [and α s (M 2 Z ) changes from 0.1191 to 0.1189]. Therefore, we can conclude that the influence of the extension of the Adler function beyond the order a 5 (by our renormalon-motivated model) does not have a numerically important role in the determination of α s . The case N t = 5 is properly included in the uncertainties of the results Eqs. (40).
The effects of the quark hadron duality violations (DVs) may be important for specific moments in the sum rules, as argued in [42,67] where a DV model [23] was used. On the other hand, we extracted our results from the Borel-Laplace sum rules whose weight function Eq. (26a) is suppressed at the Minkowskian point Q 2 = −σ m through the doublepinching factor (1+Q 2 /σ m ) 2 and, in addition, the exponential factor exp(Q 2 /M 2 ). The exponential factor suppression is especially strong and extended around the mentioned Minkowskian point when |M 2 | is small (|M 2 | ≈ 1 GeV 2 ), i.e., at those values of |M 2 | where the nonperturbative effects (D = 4, 6 terms) in the Borel-Laplace sum rule are large. This indicates that the applied Borel-Laplace sum rules are designed to suppress the DV effects more strongly when these DV effects are stronger.
The results were obtained by an analysis based on programs written by us in Mathematica that are freely available from [68]. A (A > 0); C ± ε are the linear paths parallel to the negative axis (z = −s ± iε, −A ≤ s ≤ −δ; −C δ is the (clockwise) circular path of radius δ. The limit ε δ → +0 is taken, i.e., first ε → +0 and then δ → +0.
When this is combined with the integral where Γ(a, A) is the incomplete Gamma function (we have here a = 0), and γ E is the Euler-Mascheroni constant (γ E ≈ 0.5772). When A is not positive real but complex, these same formulas remain valid by complex continuation in A. It can be checked by numerical evaluation (integration) of the integrals J s (A) Eq. (31) for complex values of A that the explicit expressions (A7) are valid. The general expression for the contribution of the dimension D = 2k (k ≥ 2) operators to the (double-pinched) FESR momenta a (2,n) (σ m ) (n ≥ 0) is where we took one-loop running of 1/a(Q 2 ) on the contour Q 2 = σ m e iφ around σ m , Eq. (30). Evaluation of this integral then gives a(σ m ) n + 3 n + 1 δ k,n+2 + δ k,n+3 (−1) n + O 2k β 0 (−1) k n + 3 n + 1