text:
# 1 Introduction

## 1 Introduction

Event-shape distributions, as well as other observables that are not completely inclusive, tend to have significant contributions from soft and collinear gluon emission. Perturbatively, this appears in the form of logarithmically enhanced terms, the well-known Sudakov logs. Another, more general source of enhancement of perturbative coefficients, which is associated with the running coupling, is infrared renormalons. These appear due to integration over particularly small momenta and therefore reflect the sensitive of the observable to large distance physics. The perturbative expansion of an event-shape distribution has quite a rich structure incorporating both sources of enhancement. Resummation techniques exist for both Sudakov logs and renormalons, separately. Yet, in many cases, both are important. The Dressed Gluon Exponentiation (DGE) method advocated here takes into account both these aspects.

It is understood how infrared renormalons are related to power-corrections through sensitivity to the phase-space boundaries . It is also known that Sudakov logs appear due to the same type of sensitivity . However, the precise relation between the two phenomena was not addressed before, except in the specific case of Drell-Yan production , where the emphasis was on the identification of the leading power-correction. DGE allows one to address this relation directly, and thus analyse the distribution in the peak region where the hierarchy between the renormalons starts to break down.

Consider a generic event-shape variable $`y`$ which is infrared and collinear safe , where the limit $`y0`$ corresponds to the extreme two-jet configuration. The logarithmically enhanced part of the differential cross-section has the form :

$$\frac{1}{\sigma }\frac{d\sigma }{dy}|_{\mathrm{log}}\frac{1}{y}\underset{n=1}{\overset{\mathrm{}}{}}\underset{m=0}{\overset{2n1}{}}c_{n,m}\left(\mathrm{ln}\frac{1}{y}\right)^ma^n,$$
(1)
where $`a\alpha _s(Q^2)/\pi `$. When $`y1`$ a naive fixed-order calculation is not adequate and the large logs must be resummed. However, this does not imply that the full perturbative calculation has to be carried out to all orders: the factorisation of soft and collinear radiation allows one to calculate the log-enhanced part of the cross-section, roughly speaking, “exponentiating the one-loop result”.

It is useful to define the integrated cross-section

$$R(y)_0^y\frac{1}{\sigma }\frac{d\sigma }{dy}𝑑y$$

where the singularity in the two-jet limit, $`y0`$, is cancelled by virtual corrections, so that $`R(0)=0`$. The logarithmically enhanced part of $`R(y)`$ has the form $`R(y)|_{\mathrm{log}}1+_{n=1}R_n(y)a^n`$ with

$$R_n_{}(y)=\underset{m=1}{\overset{2n}{}}r_{n,m}L^m$$

where $`L\mathrm{ln}\frac{1}{y}`$ and $`r_{n,m}=c_{n,m1}/m`$. The full $`R(y)`$ canThere are various different ways to make this separation. The one we use here is motivated by the so-called log-R matching scheme. See . be written  as

$$\mathrm{ln}R(y)=\mathrm{ln}R(y)|_{\mathrm{log}}+\mathrm{ln}R(y)|_{\mathrm{non}\mathrm{log}},$$
(2)
where the first term contains only logarithmic terms that diverge in the two-jet limit whereas the second is finite in this limit.

Exponentiation means that $`\mathrm{ln}R(y)|_{\mathrm{log}}`$ has the following perturbative expansion,

$`\mathrm{ln}R(y)|_{\mathrm{log}}`$ $`=`$ $`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}{\displaystyle \underset{m=1}{\overset{n+1}{}}}G_{n,m}L^ma^n`$
$`=`$ $`Lg_1(aL)+g_2(aL)+ag_3(aL)+a^2g_4(aL)+\mathrm{}`$
where the function $`g_1(aL)`$ resums the leading-logs (LL) $`G_{n,n+1}L^{n+1}a^n`$ to all orders, $`g_2(aL)`$ resums the next-to-leading logs (NLL) $`G_{n,n}L^na^n`$, etc. Exponentiation is a universal feature which follows from the factorisation property of QCD matrix elements (see e.g. ), provided that the phase-space integration also factorizes. The latter condition makes the calculation of the exact coefficients $`G_{n,m}`$ rather difficult. However, thanks to an intensive research effort in the last two decades, there is today a detailed understanding of Sudakov logs in a large class of QCD observables. In particular, in the case of event-shape observables in $`e^+e^{}`$ annihilation, like the thrust, the jet mass parameters, the C-parameter, and jet broadening parameters, the state of the art is resummation of $`\mathrm{ln}R(y)|_{\mathrm{log}}`$ up to NLL accuracy .

Current phenomenology of event shape distributions (see  for a recent review) is based on NLL resummation, neglecting $`g_k(aL)`$ for $`k3`$ in (1). The resummed cross-section is matched with next-to-leading order calculations that are available numerically . This combination allows a good qualitative description of the observables, with the main caveat being the fact that “hadronization corrections”, which are not under control theoretically, are required to bridge the gap between the calculation and the measurements.

Before alluding to non-perturbative physics it is worthwhile to re-examine the perturbative calculation. It was already shown in the example of the average thrust, that a major part of the discrepancy between the next-to-leading order calculation and the data can be explained by resummation of higher-order perturbative terms associated with infrared renormalons . Physically the Single Dressed Gluon (SDG) resummation of  allows a better description of the branching of a single gluon emitted from the primary quark–anti-quark pair. Following the spirit of BLM , this resummation takes into account the fact that the physical scale of the emitted gluon is its virtuality, rather than the centre of mass energy $`Q`$, which is usually used as the scale of the coupling. As emphasized in , the essence of this approach is the idea that the perturbative expansion should be reorganized in a “skeleton expansion” in analogy with the Abelian theory. Then, the separation between running coupling effects and conformal effects is unique and physical. We shall see here, that in the specific context of the Sudakov region, one has access to sub-leading terms in the “skeleton expansion”, corresponding to multiple gluon emission, through exponentiation.

In order to examine the approximations in the standard NLL perturbative treatment of the thrust distribution, let us first recall two main features of the resummation (1). The first is the fact that the calculation of the functions $`g_k(aL)`$ is based on an integral over the running coupling. The second is that each $`g_k(aL)`$ resums a convergent series in $`aL`$ (for $`aL1`$). At first sight these facts seem contradictory: usually integrals over the running coupling translate into renormalon factorial divergence. The resolution is that the expansion (1) itself is divergent. The functions $`g_k(aL)`$ increase factorially with $`k`$, thus endangering the validity of the truncation at a certain logarithmic accuracy, e.g. NLL. Consequently, when examining the logs $`G_{n,m}L^m`$ at a fixed order $`n`$ in eq. (1), one finds that the sub-leading logs are enhanced by a relative factor of $`m!`$ compared to the leading ones. In the region where the perturbative treatment holds ($`aL1`$) this numerical factor can over-compensate the power of the log and thus invalidate the expansion.

Another (related) argument for keeping track of sub-leading logs is the fact that integrals over the running coupling yield at once all powers of $`L`$. Therefore, truncation at a certain logarithmic accuracy implies renormalization scheme and scale dependence. A renormalization-group invariant calculation of $`\mathrm{ln}R(y)|_{\mathrm{log}}`$ requires keeping all the logs. Of course, the exact calculation of all the logs is far beyond reach. However, a specific class of logs is fully accessible: these are the logs that are leading in $`\beta _0`$, $`G_{n,m}\beta _0^{n1}`$, namely the ones associated with a single dressed gluon (SDG) emission. Within this class all the sub-leading logs, i.e. the functions $`g_k(aL)`$, are calculable. Moreover, since the exponentiation kernel is primarily associated with a single gluon, the entire logarithmically enhanced cross-section, which includes any number of gluons, can be generated in some approximation by exponentiating the SDG cross-section. The result is fully consistent with NLL accuracy, provided an adequate choice of the coupling (the “gluon bremsstrahlung” scheme ) is made. Our exponentiation procedure is based on the standard assumption that gluons emitted from the primary quark and anti-quark are completely independent and contribute additively to the thrust (this assumption holds for soft gluons). Given this assumption the SDG cross-section exponentiates under the Laplace transform at any logarithmic accuracy. Thus, through integrals over the running coupling we are able to calculate $`\mathrm{ln}R(y)|_{\mathrm{log}}`$ with the same formal logarithmic accuracy (NLL), yet avoiding any truncation of sub-leading logs, which would inevitably violate renormalization group invariance. In this way we resum at once Sudakov logs and renormalons.

Comparison of the perturbatively calculated event-shape distribution with data requires some phenomenological model to deal with hadronization. Traditionally, hadronization effects are included using detailed Monte-Carlo hadronization models such as HERWIG  or PYTHIA . An attractive alternative that became popular recently is to parametrise the effects in terms of a small number of non-perturbative parameters which control power-suppressed contributions. The basic assumption, which we adopt, is that the form of the most important non-perturbative corrections can be deduced from the ambiguities of the perturbative result . It was shown  that the primary non-perturbative effect is a shift of the perturbatively calculated distribution to larger values of $`y`$:

$$\frac{1}{\sigma }\frac{d\sigma }{dy}|_{\text{PT}}(y)\frac{1}{\sigma }\frac{d\sigma }{dy}|_{\text{PT}}(y\lambda _1\mathrm{\Lambda }/Q)$$
(4)
where $`\lambda _1`$ is a non-perturbative parameter. According to the model of , this parameter is expressed as an observable dependent factor which is perturbatively calculable times the “effective infrared coupling” (see ). The shift (4) allows a fair description of the data at all measured energies and in a wide range of $`y`$ values for several event-shapes in terms of a single parameter – the “effective infrared coupling”. The success of this perturbatively motivated approach suggests that the perturbative calculation itself, if pushed further, may be much closer to the data. It should be emphasised though that non-perturbative effects are essential in the two-jet limit when $`aL𝒪(1)`$, and in particular in determining the precise location of the distribution peak and its shape. In this region the above-mentioned approach, which uses a single non-perturbative parameter, is insufficient. In  a more general approach was suggested, namely to replace the infinite sum of power-corrections of the form $`\lambda _n\mathrm{\Lambda }^n/(yQ)^n`$ which exponentiate by a single shape-function whose central moments are $`\lambda _n`$.

As a consequence of infrared renormalons the resummation of $`\mathrm{ln}R(y)|_{\mathrm{log}}`$ contains power-suppressed ambiguities which signal the necessity of non-perturbative corrections. The fact that renormalons occur in the exponent strongly suggest that non-perturbative power corrections also exponentiate, in accordance with the shift and the shape-function approaches. We show that the renormalon ambiguity contains valuable information on the non-perturbative corrections, which is not restricted to the level of the leading power correction. In particular, it is possible to deduce the $`y`$ dependenceThe idea that the dependence of power-corrections on external variables can be deduced from the functional form of the renormalon ambiguity is not new. It was used e.g. in  for the case of higher-twist in deep-inelastic structure functions and in  for the case of power corrections to fragmentation functions in $`e^+e^{}`$ annihilation. In the latter reference (see also ) it was justified by the ‘ultraviolet dominance’ assumption. Here we use this idea at the level of the exponent. of each power correction $`1/Q^n`$ from the corresponding residue of the Borel transform of the exponent. This allows us to make definite predictions concerning the form of the non-perturbative corrections. In the specific case of the thrust distribution the main conclusion from this analysis is that the even central moments of the shape-function are suppressed. Moreover, in the case of the second central moment, an additional cancellation occurs between large-angle contributions, which are included in the shape-function of refs. , and a collinear contribution which was not addressed before.

Another important aspect which our approach incorporates is the complementarity between perturbative and non-perturbative contributions, namely renormalon resummation and explicit parametrisation of power corrections . It was clearly shown in  that a consistent treatment of non-perturbative power-corrections as well as a reliable determination of $`\alpha _s`$ requires renormalon resummation. Such resummation is performed here for the first time at the level of the differential cross-section.

This paper deals with a specific observable, the thrust distribution. Nevertheless, DGE is a general method that can be applied in any case where Sudakov logs appear. It makes a direct link between Sudakov resummation and renormalons. This paper is organised as follows: we begin in section 2 by recalling the calculation of the thrust distribution in the single dressed gluon approximation . In principle, this calculation does not require any specific kinematic approximations, however, it is technically simpler if one neglects  non-inclusive decay  of the gluon into opposite hemispheres. Here we adopt this approximation and provide further evidence for its validity at small values of $`1`$ thrust. We then observe that log-enhanced terms dominate the single dressed gluon result in a large range of thrust values. This justifies concentrating on the log-enhanced cross-section of the form (1). In section 3 we calculate the resummed cross-section by exponentiating the logarithmically enhanced part of the SDG result. We also analyse there the renormalon structure of the exponent and the relation with power-correction models mentioned above. In section 4 we compare our approach with known results. First we refer to the NLL calculation which is based on solving an evolution equation for the jet mass distribution . We show that the two calculations coincide up to NLL accuracy, provided an appropriate choice of the coupling is made . Next we briefly recall the matching procedure of the resummed result with the full NLO calculation. In section 5 we study the phenomenological implication of the suggested resummation by fitting to data in a large range of energies and in section 6 we summarise our conclusions.

## 2 Thrust distribution from Single Dressed Gluon

The purpose of this section is to calculate the thrust distribution in case of a single gluon emission. The result will be used in section 3 to calculate the physical multi-gluon cross-section, through exponentiation. We begin by recalling the calculation of the Single Dressed Gluon characteristic function (sec. 2.1). Then, in sec. 2.2, we isolate the logarithmically enhanced terms and study their structure and finally, in sec 2.3, we construct the Borel representation of the result.

### 2.1 Single Dressed Gluon (SDG) characteristic function

In renormalon calculations  the all-order perturbative result, which resums all the terms that are leading in $`\beta _0`$, can be expressed as a time-like momentum integral of an observable dependent Single Dressed Gluon (SDG) characteristic function, times an observable independent effective-charge:

$$_0^{\mathrm{}}\frac{dϵ}{ϵ}\dot{}(ϵ)\overline{A}_{\text{eff}}(ϵQ^2)=_0^{\mathrm{}}\frac{dϵ}{ϵ}\left[(ϵ)(0)\right]\overline{\rho }(ϵQ^2)$$
(5)
where $`\overline{A}(k^2)\beta _0\overline{a}(k^2)=\beta _0\overline{\alpha }_s(k^2)/\pi `$, and the bar stands for a specific renormalization-scheme. In the large $`\beta _0`$ limit (large $`N_f`$ limit) $`\overline{A}(k^2)`$ is related to the $`\overline{\mathrm{MS}}`$ coupling $`A(k^2)\beta _0\alpha _s^{\overline{\text{MS}}}(k^2)/\pi `$ by

$$\overline{A}(k^2)=\frac{A(k^2)}{1\frac{5}{3}A(k^2)}.$$
(6)
Going beyond this limit is discussed in sec. 4, where we identify $`\overline{A}`$ with the “gluon bremsstrahlung” coupling , thus making the DGE result exact to NLL accuracy.

The two integrals in (5) are related by integration by parts:

$$\dot{}(ϵ)ϵ\frac{d}{dϵ}(ϵ),$$
(7)
and the function $`\overline{\rho }(\mu ^2)`$ is identified as the time-like discontinuity of the coupling,

$$\overline{\rho }(\mu ^2)=\frac{1}{2\pi i}\mathrm{Disc}\left\{\overline{A}(\mu ^2)\right\}\frac{1}{2\pi i}\left[\overline{A}\left(\mu ^2+iϵ\right)\overline{A}\left(\mu ^2iϵ\right)\right].$$
(8)
The “time-like coupling” $`\overline{A}_{\text{eff}}(\mu ^2)`$ in (5) obeys

$$\mu ^2\frac{d\overline{A}_{\text{eff}}(\mu ^2)}{d\mu ^2}=\overline{\rho }(\mu ^2).$$
(9)
For example, in the one-loop case $`\overline{A}(k^2)=1/(\mathrm{ln}k^2/\overline{\mathrm{\Lambda }}^2)`$ and

$$\overline{A}_{\text{eff}}(\mu ^2)=\frac{1}{2}\frac{1}{\pi }\mathrm{arctan}\left(\frac{1}{\pi }\mathrm{ln}\frac{\mu ^2}{\overline{\mathrm{\Lambda }}^2}\right).$$
(10)
Eq. (5) which represents the all-order perturbative sum is not yet defined to power accuracy. The precise definition will be given in sec. 3.3 using Borel summation. A delicate issue, which was addressed in  is the fact that the dispersive integral (5) over the time-like perturbative coupling (e.g. (10)) is formally well-defined but it differs from the Borel sum by power corrections, which are not entirely related to large distances. In  the dispersive integral is given non-perturbative meaning by replacing the perturbative coupling by a non-perturbative coupling which is regular in the space-like region (i.e. it does not have a Landau singularity) and obeys a dispersion relation. Here we do not make this replacement. We refrain from attaching any non-perturbative meaning to (5). Large distance power-corrections will be deduced directly from the ambiguity of the Borel sum (sec. 3.5). Eventually, the normalisation of these power corrections can be expressed as an integral over an infrared-finite coupling, which is conjectured  to be universal.

The characteristic function $`(ϵ)`$ itself can be calculated either using the Abelian large $`N_f`$ limit followed by naive non-Abelianization , or in the dispersive approach , referring to $`\mu ^2=ϵQ^2`$ as the squared gluon mass. In case of inclusive observables, such as the total cross-section in $`e^+e^{}`$ annihilation, the two calculations coincide. The two differ , however, in the case of observables such as event-shape variables, that are not completely inclusive with respect to the decay products of the gluon. Here we shall use the massive-gluon characteristic function for the thrust distribution, which was calculated in . It corresponds to an inclusive integration over the decay products of the gluon, so that the possibility that a gluon decays into partons that end up in opposite hemispheres is not taken into account. It is important to emphasise, though, that the steps that follow can be repeated in very much the same way starting with the large $`N_f`$ result, once available. In the case of the thrust, the inclusive approximation is expected to be good (some evidence is given below). Note that for some event-shapes, e.g. the heavy jet mass, this is probably not so.

Let us now briefly recall the calculation of the thrust distribution characteristic function of ref. . At the SDG level, the distribution of $`t=1`$ thrust is

$$\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}=\frac{C_F}{\beta _0}_0^1\frac{dϵ}{ϵ}\overline{A}_{\text{eff}}(ϵQ^2)\dot{}(ϵ,t)=\frac{C_F}{\beta _0}_0^1\frac{dϵ}{ϵ}\overline{\rho }(ϵQ^2)\left[(ϵ,t)(0,t)\right].$$
(11)
The characteristic function $`(ϵ,t)`$ is obtained from the following integral over phase-space,

$$(ϵ,t)=_{\mathrm{phase}\mathrm{space}}𝑑x_1𝑑x_2(x_1,x_2,ϵ)\delta \left(1T(x_1,x_2,ϵ)t\right)$$
(12)
where $`C_Fa`$ is the squared tree level matrix element for the production of a quark–anti-quark pair and a gluon of virtuality $`\mu ^2ϵQ^2`$, and

$$(x_1,x_2,ϵ)=\frac{1}{2}\left[\frac{(x_1+ϵ)^2+(x_2+ϵ)^2}{(1x_1)(1x_2)}\frac{ϵ}{(1x_1)^2}\frac{ϵ}{(1x_2)^2}\right].$$
(13)
The integration variables $`x_{1,2}`$ represent the energy fraction of each of the quarks in the centre-of-mass frame. The energy fraction of the gluon is $`x_3=2x_1x_2`$.

We use the following  definition of the thrust

$$T=\frac{_i\left|\stackrel{}{p}_i\stackrel{}{n}_T\right|}{_iE_i}=\frac{_i\left|\stackrel{}{p}_i\stackrel{}{n}_T\right|}{Q}.$$
(14)
In case of three partons (a quark, an anti-quark and a “massive” gluon) it yields ,

$$1T(x_1,x_2,ϵ)=\mathrm{min}\{1x_1,\mathrm{\hspace{0.17em}1}x_2,\mathrm{\hspace{0.17em}1}\sqrt{(2x_1x_2)^24ϵ}\}.$$
(15)
Evaluating (12) one obtains  the characteristic function for the thrust distribution,

$$(ϵ,t)=\{\begin{array}{ccc}_Q^l(ϵ,t)+_G(ϵ,t)\hfill & & ϵ<t<\sqrt{ϵ}\hfill \\ _Q^h(ϵ,t)+_G(ϵ,t)\hfill & & \sqrt{ϵ}<t<\frac{2}{3}\sqrt{1+3ϵ}\frac{1}{3}\hfill \end{array}$$
(16)
where the dominant contribution $`_Q(ϵ,t)`$ corresponds to the phase-space regions where one of the primary quarks carries the largest momentum ($`T=x_{1,2}`$) and $`_G(ϵ,t)`$ corresponds to the region where the gluon momentum is the largest. The superscripts $`l`$ and $`h`$ on $`_Q(ϵ,t)`$ denote low and high $`t`$ values, respectively. These functions are given by

$`_Q^h(ϵ,t)`$ $`=`$ $`{\displaystyle \frac{1}{t}}\left[(1t+ϵ)^2+(1+ϵ)^2\right]\mathrm{ln}{\displaystyle \frac{t}{qt}}+\left(32{\displaystyle \frac{q}{t}}+{\displaystyle \frac{1}{2}}{\displaystyle \frac{1}{t}}+{\displaystyle \frac{1}{2}}tq\right)`$
$`+`$ $`\left(42{\displaystyle \frac{q}{t}}+3{\displaystyle \frac{1}{t}}{\displaystyle \frac{q}{t^2}}+{\displaystyle \frac{1}{qt}}\right)ϵ`$
$`_Q^l(ϵ,t)`$ $`=`$ $`{\displaystyle \frac{1}{t}}\left[(1t+ϵ)^2+(1+ϵ)^2\right]\mathrm{ln}{\displaystyle \frac{ϵ}{t(qt)}}+\left(12{\displaystyle \frac{q}{t}}+{\displaystyle \frac{1}{2}}{\displaystyle \frac{1}{t}}q\right)`$ (17)
$`+`$ $`\left(3{\displaystyle \frac{1}{t}}+{\displaystyle \frac{1}{qt}}{\displaystyle \frac{q}{t^2}}+2{\displaystyle \frac{1}{t^2}}+22{\displaystyle \frac{q}{t}}\right)ϵ+\left(2+{\displaystyle \frac{1}{2t}}\right){\displaystyle \frac{ϵ^2}{t^2}}`$
$`_G(ϵ,t)`$ $`=`$ $`{\displaystyle \frac{1t}{q^2}}\left[\left((1+ϵ)^2+(1+ϵq)^2\right)\mathrm{ln}{\displaystyle \frac{qt}{t}}+(2tq)\left(q+{\displaystyle \frac{ϵ}{t}}+{\displaystyle \frac{ϵ}{qt}}\right)\right]`$
where $`q\sqrt{(1t)^2+4ϵ}`$.

### 2.2 Logarithmically Enhanced Terms

Next we examine the SDG perturbative expression (11) with the characteristic function of eqs. (16) and (2.1). For this purpose it is convenient to expand the time-like coupling $`\overline{A}_{\text{eff}}(ϵQ^2)`$ in terms of a fixed-scale space-like coupling $`\overline{A}(Q^2)`$,

$$\overline{A}_{\text{eff}}(ϵQ^2)=\underset{m=0}{\overset{\mathrm{}}{}}\left(\overline{A}(Q^2)\right)^{m+1}\underset{j=0}{\overset{m}{}}K_{m,j}\left(\mathrm{ln}\frac{1}{ϵ}\right)^j,$$
(18)
where, in general, $`K_{m,j}`$ depend on the coefficients of the $`\beta `$ function which controls the evolution of $`\overline{A}(Q^2)`$. Performing the integral (11) over $`\dot{}(ϵ,t)ϵ\frac{d}{dϵ}(ϵ,t)`$ order by order one gets the sum

$$\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}=\frac{C_F}{\beta _0}\underset{m=0}{\overset{\mathrm{}}{}}\left(\overline{A}(Q^2)\right)^{m+1}\underset{j=0}{\overset{m}{}}K_{m,j}h_j(t)$$
(19)
where the functions $`h_j(t)`$ are the log-moments of the characteristic function,

$$h_j(t)=_0^1\frac{dϵ}{ϵ}\left(\mathrm{ln}\frac{1}{ϵ}\right)^j\dot{}(ϵ,t).$$
(20)
The analytic calculation of $`h_j(t)`$ with the full characteristic function of eqs. (16) and (2.1) is rather involved. Alternatively, concentrating in the small $`t`$ region, we distinguish in $`h_j(t)`$ between logarithmically enhanced terms which are not suppressed by any power of $`t`$ and the rest, which we eventually neglect.

As we shall see, the logarithmically enhanced terms, which have the form of eq. (1), can be easily computed analytically to all orders. Using (16) for $`t<1/3`$ we obtain

$$h_j(t)=_0^t\frac{dϵ}{ϵ}\left(\mathrm{ln}\frac{1}{ϵ}\right)^j\dot{}_G(ϵ,t)+_0^{t^2}\frac{dϵ}{ϵ}\left(\mathrm{ln}\frac{1}{ϵ}\right)^j\dot{}_Q^h(ϵ,t)+_{t^2}^t\frac{dϵ}{ϵ}\left(\mathrm{ln}\frac{1}{ϵ}\right)^j\dot{}_Q^l(ϵ,t).$$
(21)
We notice that a term of the form $`ϵ^n/t^k`$ in $`\dot{}(ϵ,t)`$ contributes to $`h_j(t)`$ a term proportional to $`t^{nk}`$ (times a polynomial in $`\mathrm{ln}\frac{1}{t}`$) from the collinear gluon emission limit of phase-space<sup>§</sup><sup>§</sup>§A detailed analysis of the phase-space can be found in . See figure 1 there. where $`ϵ=t`$, and a term proportional to $`t^{2nk}`$ (times a polynomial in $`\mathrm{ln}\frac{1}{t}`$) from the large-angle soft gluon emission limit of phase-space where $`ϵ=t^2`$. Taking the double expansion of (2.1) at small $`t`$ and $`ϵ`$ we have

$`\dot{}_Q^h(ϵ,t)`$ $``$ $`{\displaystyle \frac{ϵ}{t^2}}+{\displaystyle \frac{𝒪\left(ϵ\right)}{t}}`$ (22)
$`\dot{}_Q^l(ϵ,t)`$ $``$ $`{\displaystyle \frac{ϵ^2}{t^3}}{\displaystyle \frac{ϵ}{t^2}}+{\displaystyle \frac{2}{t}}+𝒪\left(1\right)`$
$`\dot{}_G(ϵ,t)`$ $``$ $`{\displaystyle \frac{ϵ}{t}}+𝒪\left(1\right)`$
where we ignored higher powers of $`ϵ`$ multiplying a given power of $`1/t`$, which are not relevant. From (21) and (22) we conclude that only $`\dot{}_Q^l(ϵ,t)`$ contributes to the logarithmically enhanced cross-section and thus

$$h_j(t)|_{\mathrm{log}}=_{t^2}^t\frac{dϵ}{ϵ}\left(\mathrm{ln}\frac{1}{ϵ}\right)^j\left[\frac{ϵ^2}{t^3}\frac{ϵ}{t^2}+\frac{2}{t}\right].$$
(23)
Performing the integral over $`ϵ`$ yields

$`h_j(t)|_{\mathrm{log}}`$ $`=`$ $`{\displaystyle \frac{1}{t}}\left[2{\displaystyle \frac{2^{j+1}1}{j+1}}L^{j+1}e^L\mathrm{\Gamma }(j+1,L)\mathrm{\hspace{0.17em}2}^{j1}e^{2L}\mathrm{\Gamma }(j+1,2L)\right]`$ (24)
$`=`$ $`{\displaystyle \frac{1}{t}}\left[2{\displaystyle \frac{2^{j+1}1}{j+1}}L^{j+1}{\displaystyle \underset{k=0}{\overset{j}{}}}{\displaystyle \frac{j!}{k!}}\left(1+2^{kj1}\right)L^k\right]`$
where $`L\mathrm{ln}\frac{1}{t}`$.

Having evaluated the log-moments we can substitute them into (19) to obtain the logarithmically enhanced SDG differential cross-section to arbitrarily high order. The result (24) has a few features that deserve attention. First we note that the leading-logs in $`h_j(t)`$ depend on both collinear and large-angle soft gluon limits of phase-space, while the sub-leading logs are associated with the collinear limit alone. Next, we observe that the leading-logs in $`h_j(t)`$ originate in the $`ϵ`$ independent term $`2/t`$ in $`\dot{}_Q^l(ϵ,t)`$, while the sub-leading logs originate in the two $`ϵ`$ dependent terms, $`ϵ^2/t^3ϵ/t^2`$. This $`ϵ`$ dependence is reflected in factorial increase of the coefficients (see below) which is associated with the sensitivity of the thrust to collinear gluon emission.

It is important to notice that the factorial growth is associated with sub-leading logs. Considering the different logs $`L^k`$ in a given log-moment $`h_j(t)`$ we find that with decreasing power $`k`$ the numerical coefficient increases as $`j!/k!`$. Note, in particular, that the next-to-leading logs $`L^j`$ (which are the only sub-leading logs taken into account in the standard analysis) have a small coefficient ($`\frac{3}{2}`$), while further sub-leading logs (which are usually ignored) are multiplied by large numbers. It is useful to analyse the square brackets in (24) by considering two extreme limits for the function $`F(L)e^L\mathrm{\Gamma }(j+1,L)`$. In the limit $`L0`$, one finds factorial growth $`F(L)j!`$. In the limit $`L\mathrm{}`$ one can keep only the leading logarithmic term in $`F(L)`$, i.e. $`F(L)L^j`$, recovering the standard NLL result. The fact that the latter formally holds may be of little relevance in practice given that the entire perturbative treatment is adequate only when $`L1/\overline{A}(Q^2)`$. One can now examine the validity of these two extreme approximations for phenomenologically interesting values of $`L`$. Consider for example the case $`t=0.05`$, $`L=\mathrm{ln}(1/0.05)3`$ at order $`j=3`$: the value is $`F(L)=77.8`$ while the NLL gives $`L^j=26.9`$ and the naive factorial gives $`j!=6`$. Clearly, none of these extreme approximations is valid. One must retain all the logs.

One may wonder to what extent it is justified to consider only the logarithmically enhanced part of the log moments. Factorial increase is certainly not restricted to the logarithmically enhanced part. The immediate answer is, of course, that the neglected part is parametrically smaller, by a factor of $`t`$ compared to the one kept. This answer is not complete: as we know, parametric suppression can be compensated by large numerical coefficients. To be convinced that this approximation is valid at physically interesting values of $`t`$ let us now compare the logarithmic part of the log-moments $`h_j(t)`$ with the full function which we evaluate numerically. Writing $`h_j(t)=h_j(t)|_{\mathrm{log}}+\mathrm{\Delta }h_j(t)`$, where the logarithmically enhanced part $`h_j(t)|_{\mathrm{log}}`$ is given by (24) and $`\mathrm{\Delta }h_j(t)`$ are suppressed by a relative factor of $`t`$, the first few log-moments are

$`h_1(t)`$ $`=`$ $`{\displaystyle \frac{1}{t}}\left[3L^2{\displaystyle \frac{3}{2}}L{\displaystyle \frac{5}{4}}+\mathrm{\Delta }h_1(t)\right]`$
$`h_2(t)`$ $`=`$ $`{\displaystyle \frac{1}{t}}\left[{\displaystyle \frac{14}{3}}L^3{\displaystyle \frac{3}{2}}L^2{\displaystyle \frac{5}{2}}L{\displaystyle \frac{9}{4}}+\mathrm{\Delta }h_2(t)\right]`$
$`h_3(t)`$ $`=`$ $`{\displaystyle \frac{1}{t}}\left[{\displaystyle \frac{15}{2}}L^4{\displaystyle \frac{3}{2}}L^3{\displaystyle \frac{15}{4}}L^2{\displaystyle \frac{27}{4}}L{\displaystyle \frac{51}{8}}+\mathrm{\Delta }h_3(t)\right]`$
$`h_4(t)`$ $`=`$ $`{\displaystyle \frac{1}{t}}\left[{\displaystyle \frac{62}{5}}L^5{\displaystyle \frac{3}{2}}L^45L^3{\displaystyle \frac{27}{2}}L^2{\displaystyle \frac{51}{2}}L{\displaystyle \frac{99}{4}}+\mathrm{\Delta }h_4(t)\right]`$
The $`h_j(t)`$ and thus $`\mathrm{\Delta }h_j(t)`$ are calculated by integrating (20) numerically where $`\dot{}(ϵ,t)`$ is the derivative of the full characteristic function (2.1). The comparison for the first four log-moments is presented in fig. 1.

We see that while the NLL approximation is quite far from the actual values of $`h_j(t)`$ even at rather small $`t`$, there is a fairly wide range in which the logarithmic approximation $`h_j(t)|_{\mathrm{log}}`$ is good. At large $`t`$ values the gap between the two increases. However, since our interest here is in improving the calculation in the small $`t`$ region, we consider this agreement satisfactory. It should be kept in mind, though, that our resummation procedure is not the appropriate way to improve the calculation of the cross-section at $`t\begin{array}{c}>\hfill \\ \hfill \end{array}0.2`$. We comment that for higher log-moments $`j5`$ the approximation deteriorates even at smaller $`t`$ values; however, the resulting discrepancy at the level of the SDG cross-section (19) is not large. We shall revisit the significance of terms that are suppressed by $`t`$ in the phenomenological analysis in sec. 5.

Finally, using (6), the first few orders of the SDG calculation (19) of the log-enhanced cross-section are given in the $`\overline{\mathrm{MS}}`$ scheme by

$`{\displaystyle \frac{1}{\sigma }}{\displaystyle \frac{d\sigma }{dt}}(Q^2,t)|_{\text{SDG}}={\displaystyle \frac{C_F^{}}{t\beta _0}}\{(2L1.5)A(Q^2)+(3L^2+1.8333L3.75)A(Q^2)^2`$ (25)
$`+(4.6667L^3+8.5000L^28.5242L5.6481)A(Q^2)^3`$
$`+(7.5000L^4+21.833L^315.859L^240.585L+2.025)A(Q^2)^4`$
$`+(12.400L^5+48.500L^429.339L^3175.72L^247.164L+40.22)A(Q^2)^5+\mathrm{}\}`$
The leading and next-to-leading logs are consistent with previous calculations . In addition, at NLO we can compare the next-to-next-to-leading log ($`3.75`$ in (25)) with the exact numerical computation . To make an explicit comparison, one should separate the contributions of different colour factors in the NLO coefficient (this was done before using the EVENT program  in ) and then fit the numerical results to a parametric form which includes a log-enhanced part, and terms which are suppressed by powers of $`t`$. The result of such a fit, where the first two logs are fixed to their exact values, gives the following log-enhanced NLO expression in the large $`\beta _0`$ limit,

$`{\displaystyle \frac{1}{\sigma }}{\displaystyle \frac{d\sigma }{dt}}(Q^2,t)|_{\text{NLO}}={\displaystyle \frac{C_F^{}}{t\beta _0}}\left\{(2L1.5)A(Q^2)+(3L^2+1.8333L3.77\pm 0.07)A(Q^2)^2\right\}.`$ (26)
Details concerning the fit and the error estimate can be found in the Appendix. The closeness of the NNLL coefficients is encouraging. It strongly suggests that the impact of the non-inclusive splitting of the gluon into partons that end up in opposite hemispheres , which is not taken into account in our calculation, is small, at least at the logarithmic level. A similar comparison, which goes beyond the logarithmic level, has been performed for the average thrust in  and the higher moments $`t^n`$ in . It was found that the non-inclusive effect is rather small for the average thrust ($`4.4\%`$) while it increases for the higher moments. This implies that the effect is mainly important in the large $`t`$ region. Being interested here in the small $`t`$ region and in particular in the logarithmically enhanced part of the cross section, we shall neglect this effect altogether. The reader should keep in mind, however, that the evidence provided above concerns directly only the Abelian part of the coefficients. Only through the assumption that running coupling effects dominate (making the large $`\beta _0`$ approximation legitimate) this evidence becomes relevant for the non-Abelian non-inclusive effect. Moreover, in the context of power-corrections within the infrared-finite coupling approach, the impact of the non-inclusive effect in the non-Abelian case was shown to be more significant than in the Abelian one .

### 2.3 Borel representation of the SDG result

A convenient way to deal with perturbation theory to all orders is Borel summation. In this section we construct the scheme-invariant Borel representation of the logarithmically enhanced part of the SDG result. This will be used in the next section to calculate the multiple emission cross-section. We start with the scheme-invariant Borel representation of the coupling ,

$$\overline{A}(k^2)=_0^{\mathrm{}}𝑑z\mathrm{exp}\left(z\mathrm{ln}k^2/\mathrm{\Lambda }^2\right)\overline{A}_B(z).$$
(27)
In the case of one-loop running coupling $`\overline{A}_B(z)=\mathrm{exp}\left(z\mathrm{ln}\mathrm{\Lambda }^2/\overline{\mathrm{\Lambda }}^2\right)`$. Taking the time-like discontinuity  we get

$$\overline{A}_{\text{eff}}(Q^2)=_0^{\mathrm{}}𝑑z\mathrm{exp}\left(z\mathrm{ln}Q^2/\mathrm{\Lambda }^2\right)\frac{\mathrm{sin}\pi z}{\pi z}\overline{A}_B(z).$$
(28)
In eq. (11) we use the part of $`\dot{}(ϵ,t)`$ that contributes to log-enhanced terms (as in (23)),

$$\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}=\frac{C_F}{\beta _0}_{t^2}^t\frac{dϵ}{ϵ}\overline{A}_{\text{eff}}(ϵQ^2)\left[\frac{ϵ^2}{t^3}\frac{ϵ}{t^2}+\frac{2}{t}\right].$$
(29)
Substituting $`\overline{A}_{\text{eff}}(ϵQ^2)`$ from (28) and integrating over $`ϵ`$ we obtain the Borel representation of the logarithmically enhanced SDG cross-section,

$$\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}=\frac{C_F}{\beta _0}\frac{1}{t}_0^{\mathrm{}}𝑑zB_{\text{SDG}}(z,t)\mathrm{exp}\left(z\mathrm{ln}Q^2/\mathrm{\Lambda }^2\right)\frac{\mathrm{sin}\pi z}{\pi z}\overline{A}_B(z)$$
(30)
with

$$B_{\text{SDG}}(z,t)=\frac{2}{z}\mathrm{exp}\left(2zL\right)\left(\frac{2}{z}+\frac{1}{1z}+\frac{1}{2z}\right)\mathrm{exp}\left(zL\right)$$
(31)
where $`L\mathrm{ln}\frac{1}{t}`$. Note that the first exponent in the Borel transform is associated with the large-angle soft gluon limit of phase-space, while the second is associated with the collinear limit. The Borel integral is well-defined at the perturbative level thanks to the cancellation of the pole at $`z=0`$ between these two contributions. Moreover, provided that the coupling does not enhance the singularity (e.g. the case of one-loop coupling), the poles at $`z=1`$ and $`z=2`$ are regulated by the $`\frac{\mathrm{sin}\pi z}{\pi z}`$ factor, which originates in the analytic continuation to the time-like region.

Note that the Borel integral converges at infinity only provided that $`2L<\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2`$, or $`Qt>\overline{\mathrm{\Lambda }}`$, reflecting large-angle gluon emission sensitivity. A weaker condition appears from the collinear limit: $`L<\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2`$, or $`Q^2t>\overline{\mathrm{\Lambda }}^2`$. These convergence conditions signal the fact that the small $`t`$ region is not under control in perturbation theory, and that power corrections of the form $`\mathrm{\Lambda }^n/(Qt)^n`$ (and eventually also $`\mathrm{\Lambda }^{2n}/(Q^2t)^n`$) should be included . We will return to this point in more detail in the next section.

## 3 Exponentiation

The calculation of the previous section concerns a single gluon emission from the primary quark–anti-quark pair. Close to the two-jet limit, resummation of multiple soft and collinear gluons is essential. This resummation is achieved by exponentiating the SDG cross-section. As mentioned in the introduction, our basic assumption is that emission of several soft or collinear gluons can be regarded as independent. The total contribution to the thrust is then computed as the sum of the individual contributions of these emissions: $`t=_{k=1}^nt_k`$. For large-angle soft gluons, which give the dominant contribution at small $`t`$, additivity holds since the contribution of each gluon to $`t`$ is proportional to its transverse momentum. This assumption fails, of course, for sufficiently hard gluons. In general, the contribution of several gluons to the thrust, as any other event-shape, are correlated. Such a correlation emerges from the modification of the phase-space limits in the integration over the gluon momenta (the quarks do recoil) as well as from additional interactions between the emitted gluons. These correlations are expected to be important only at large values of $`1`$ thrust, away from the two-jet configuration. Since they may have an effect at the logarithmic level, our calculation of sub-leading logs (beyond NLL accuracy) that are sub-leading in $`\beta _0`$Such correlations do not affect the logs that are leading in $`\beta _0`$. These terms are not modified by the exponentiation and their calculation in the previous section is exact barring the inclusive approximation made in the definition of the thrust. is only partial. It remains for future work to quantify the effect of these correlations. Here we simply assume that it is small.

### 3.1 Exponentiation Formula

The result of the previous section, $`\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}`$, has a probabilistic interpretation as the probability distribution, of the random variable $`t(0,1)`$, given that only one gluon is emitted. Since, in reality any number of gluons $`n`$ can be emitted, we sum over these possibilities

$$\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{res}}=\underset{n=0}{\overset{\mathrm{}}{}}\frac{1}{n!}\underset{k=1}{\overset{n}{}}_0^1𝑑t_k\frac{1}{\sigma }\frac{d\sigma }{dt}(t_k)|_{\text{SDG}}\left[\delta \left(t\underset{k=1}{\overset{n}{}}t_k\right)\delta (t)\right]$$
(32)
where the singularity due to zero gluon emission probability $`n=0`$ is cancelled by virtual corrections, $`\delta (t)`$. The combinatorial $`n!`$ is introduced to avoid multiple counting – the expression is symmetric for any permutation of the $`n`$ gluons. Using the Laplace conjugate variable $`\nu `$, the $`\delta `$ function can be written in a factorized form as (see e.g. ),

$$\delta (t\underset{k=1}{\overset{n}{}}t_k)=_𝒞\frac{d\nu }{2\pi i}\mathrm{exp}\left\{\nu \left(\underset{k=1}{\overset{n}{}}t_kt\right)\right\}.$$
(33)
where the contour $`𝒞`$ runs parallel to the imaginary axis (to the right of all the singularities of the integrand, if any). Changing the order of the $`\nu `$ integration and the sum over $`n`$ in (32) we get the exponentiated form,

$`{\displaystyle \frac{1}{\sigma }}{\displaystyle \frac{d\sigma }{dt}}(Q^2,t)|_{\text{res}}`$ $`=`$ $`{\displaystyle _𝒞}{\displaystyle \frac{d\nu }{2\pi i}}e^{\nu t}{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{n!}}\left[{\displaystyle _0^1}{\displaystyle \frac{1}{\sigma }}{\displaystyle \frac{d\sigma }{dt}}(t_k)|_{\text{SDG}}\left(e^{\nu t_k}1\right)dt_k\right]^n`$ (34)
$`=`$ $`{\displaystyle _𝒞}{\displaystyle \frac{d\nu }{2\pi i}}e^{\nu t}\mathrm{exp}\left\{{\displaystyle _0^1}{\displaystyle \frac{1}{\sigma }}{\displaystyle \frac{d\sigma }{dt}}(Q^2,\stackrel{~}{t})|_{\text{SDG}}\left(e^{\nu \stackrel{~}{t}}1\right)d\stackrel{~}{t}\right\}`$
Defining

$$𝒮_{\text{PT}}(Q^2,\mathrm{ln}\nu )e^{\nu t}_{\text{SDG}}_0^1\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,\stackrel{~}{t})|_{\text{SDG}}\left(e^{\nu \stackrel{~}{t}}1\right)d\stackrel{~}{t},$$
(35)
we have

$$\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{res}}=_𝒞\frac{d\nu }{2\pi i}e^{\nu t}\mathrm{exp}\left\{e^{\nu t}_{\text{SDG}}\right\},$$
(36)
and finally, integrating over $`t`$,

$$R(Q^2,t)|_{\text{res}}=_𝒞\frac{d\nu }{2\pi i\nu }e^{\nu t}\mathrm{exp}\left\{e^{\nu t}_{\text{SDG}}\right\}.$$
(37)
Note that the terms which are leading in $`\beta _0`$ in the coefficients of $`\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{res}}`$, or in the log-enhanced part of the exact result, are the same as the ones in $`\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}`$. These coefficients correspond to a single emission.

### 3.2 The Exponent in the Conjugate Variable

Starting with the Borel representation (30) with the Borel function (31) we now calculate the exponent $`e^{\nu t}_{\text{SDG}}`$ as defined in (35),

$`e^{\nu t}_{\text{SDG}}`$ $`=`$ $`{\displaystyle \frac{C_F}{\beta _0}}{\displaystyle _0^{\mathrm{}}}𝑑zB_\nu (z)\mathrm{exp}\left(z\mathrm{ln}Q^2/\mathrm{\Lambda }^2\right){\displaystyle \frac{\mathrm{sin}\pi z}{\pi z}}\overline{A}_B(z)`$ (38)
with

$$B_\nu (z)_0^1\frac{dt}{t}\left[\frac{2}{z}e^{2z\mathrm{ln}\frac{1}{t}}\left(\frac{2}{z}+\frac{1}{1z}+\frac{1}{2z}\right)e^{z\mathrm{ln}\frac{1}{t}}\right]\left(e^{\nu t}1\right).$$
(39)
To evaluate $`B_\nu (z)`$ we use the integral

$$_0^1\frac{dt}{t}e^{z\mathrm{ln}\frac{1}{t}}\left(e^{\nu t}1\right)=\nu ^z\gamma (z,\nu )+\frac{1}{z}$$
(40)
where

$$\gamma (z,\nu )\mathrm{\Gamma }(z)\mathrm{\Gamma }(z,\nu ).$$
(41)
Note that (40) is regular at $`z=0`$ due to cancellation between the $`\gamma `$ function and the simple pole – the first originates in the $`e^{\nu t}`$ part, i.e. in real contributions, while the latter in the unity, namely in virtual corrections. We thus find

$$B_\nu (z)=\frac{2}{z}\left[e^{2z\mathrm{ln}\nu }\gamma (2z,\nu )+\frac{1}{2z}\right]\left(\frac{2}{z}+\frac{1}{1z}+\frac{1}{2z}\right)\left[e^{z\mathrm{ln}\nu }\gamma (z,\nu )+\frac{1}{z}\right].$$
(42)
Similarly to (31), this Borel function is well defined at the perturbative level thanks to cancellation of the pole at $`z=0`$ between the two terms in (42); the first corresponds to the large-angle soft gluon emission and the second to collinear gluon emission. In contrastAt first sight it seems surprising that the Laplace-transformed $`B_\nu (z)`$ has a completely different renormalon structure from $`B_{\text{SDG}}(z,t)`$. The reason is, of course, that $`B_\nu (z)`$, thanks to the integration over $`t`$, is sensitive to the small $`t`$ region, where the exponents $`\mathrm{exp}\left(2zL\right)`$ and $`\mathrm{exp}\left(zL\right)`$ in (31) are large and impose convergence constraints on (30). with $`B_{\text{SDG}}(z,t)`$ in (31), $`B_\nu (z)`$ has a rich renormalon structure, making the integral (38) ambiguous at the level of power accuracy. Note that $`\gamma (z,\nu )`$ has the same singularity structure as $`\mathrm{\Gamma }(z)`$: a simple pole at any non-negative integer. It follows that $`B_\nu (z)`$ in (42) has two sets of infrared renormalons: the first at all positive integers and half-integers due to large-angle emission sensitivity and the second at all positive integers – due to collinear sensitivity. In the latter, the renormalons at $`z=1`$ and $`z=2`$ become double poles. It is important to note that the $`\mathrm{sin}(\pi z)`$ factor in (38) regulates simple poles and transforms double poles into simple poles at all integer $`z`$ values, while it does not affect the singularities at half integer $`z`$ values. We will return to this point in the context of power-corrections in sec. 3.5.

An alternative representation of $`e^{\nu t}_{\text{SDG}}`$ in terms of a characteristic function can be obtained using eq. (29), namely

$`e^{\nu t}_{\text{SDG}}`$ $`=`$ $`{\displaystyle \frac{C_F}{\beta _0}}{\displaystyle _0^1}𝑑t(e^{\nu t}1){\displaystyle _{t^2}^t}{\displaystyle \frac{dϵ}{ϵ}}\overline{A}_{\text{eff}}(ϵQ^2)\left[{\displaystyle \frac{ϵ^2}{t^3}}{\displaystyle \frac{ϵ}{t^2}}+{\displaystyle \frac{2}{t}}\right]`$ (43)
$`=`$ $`{\displaystyle \frac{C_F}{\beta _0}}{\displaystyle _0^1}{\displaystyle \frac{dϵ}{ϵ}}\overline{A}_{\text{eff}}(ϵQ^2)\dot{}_\nu (ϵ),`$
where the order of integrations over $`t`$ and $`ϵ`$ has been changed. The characteristic function is given by,

$`\dot{}_\nu (ϵ)`$ $`=`$ $`\left(2+\nu ϵ{\displaystyle \frac{1}{2}}\nu ^2ϵ^2\right)\mathrm{Ei}(1,\sqrt{ϵ}\nu )+\left(\nu ϵ2+{\displaystyle \frac{1}{2}}\nu ^2ϵ^2\right)\mathrm{Ei}(1,\nu ϵ)`$
$`+`$ $`\left(\sqrt{ϵ}{\displaystyle \frac{1}{2}}ϵ+{\displaystyle \frac{1}{2}}\nu ϵ^{3/2}\right)e^{\sqrt{ϵ}\nu }+\left({\displaystyle \frac{3}{2}}{\displaystyle \frac{1}{2}}\nu ϵ\right)e^{\nu ϵ}`$
$`+`$ $`2\mathrm{ln}(\sqrt{ϵ}\nu )2\mathrm{ln}(\nu ϵ)+{\displaystyle \frac{1}{2}}ϵ+\sqrt{ϵ}{\displaystyle \frac{3}{2}}.`$
Given a model<sup>*</sup><sup>*</sup>*An example is provided by . for the coupling which is regular in the infrared , the integral (43) could be readily evaluated. We prefer, however, to complete the perturbative treatment without any additional assumptions on non-perturbative physics. For this purpose we find the Borel representation (42) most suitable. In the next section we calculate this integral explicitly, but before doing so it is worthwhile to examine more closely the structure of the sum, term by term.

Since we are interested in large $`\nu `$ we can replace $`\gamma (z,\nu )`$ in (42) by $`\mathrm{\Gamma }(z)`$. The logarithmically enhanced terms are not modified by such a replacement. Expanding the integrand in (38) in powers of $`z`$, we obtain

$`B_\nu (z)={\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(\mathrm{ln}\nu )^n}{n!}}`$ $`[\mathrm{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}2}{\displaystyle \underset{m\mathrm{max}\{1,n1\}}{}}(12^m)c_{mn}z^{m1}`$
$`+{\displaystyle \underset{l0}{}}(1+2^{l1}){\displaystyle \underset{m\mathrm{max}\{n+l,n+1\}}{}}c_{mnl1}z^{m1}],`$
where the first term in the square brackets originates in the $`2/t`$ part in $`\dot{}(ϵ,t)`$ from both the soft and collinear limits of phase-space, while the second originates in the $`ϵ`$ dependent terms in $`\dot{}(ϵ,t)`$ from the collinear limit alone. The numbers $`c_k`$ are defined by $`\mathrm{\Gamma }(z)=_{k=1}^{\mathrm{}}c_kz^k`$. This means that $`c_1=1`$, $`c_0=\gamma _E`$, $`c_1=\pi ^2/12+\gamma _E^2/2`$, etc. Higher $`c_k`$ contain higher $`\zeta _i`$ numbers ($`ik+1`$), yet numerically $`c_k`$ ($`k1`$) are all close to $`1`$.

Next, performing the Borel integral term by term we get a series in the coupling $`\overline{A}(Q^2)`$. It is then straightforward to rewrite the sum as

$$e^{\nu t}_{\text{SDG}}=\frac{C_F}{\beta _0}\underset{k=1}{\overset{\mathrm{}}{}}\overline{A}(Q^2)^{k2}f_k\left(\overline{A}(Q^2)\mathrm{ln}\nu \right)$$
(46)
and perform the summation over all powers of $`\xi \overline{A}(Q^2)\mathrm{ln}\nu `$. The functions $`f_k`$ are given by

$$f_k(\xi )=\underset{u=0}{\overset{(k1)/2}{}}\frac{(\pi ^2)^u}{(2u+1)!}\left[g_k^{(u)}(\xi )M_{\mathrm{col}}(k)+g_k^{(u)}(2\xi )M_{\mathrm{la}}(k)\right]$$
(47)
where

$`M_{\mathrm{col}}(k)`$ $`=`$ $`2c_{k2}+{\displaystyle \underset{m=0}{\overset{k2}{}}}(1+2^{m1})c_{k3m}`$ (48)
$`M_{\mathrm{la}}(k)`$ $`=`$ $`2^{k1}c_{k2}`$
and

$`g_k^{(u)}(\xi )=\xi ^{n_0}{\displaystyle \frac{\mathrm{\Gamma }(k2+n_0)}{\mathrm{\Gamma }(1+n_0)}}\mathrm{Hypergeom}([1,k2+n_0],[1+n_0],\xi )`$
with $`n_0=\mathrm{max}(0,2u+3k)`$. In (48), $`M_{\mathrm{col}}(k)`$ and $`M_{\mathrm{la}}(k)`$ correspond to the collinear and large-angle contributions to the Borel function, respectively. Explicitly, this yields

$`f_1(\xi )`$ $`=`$ $`2(1\xi )\mathrm{ln}(1\xi )(12\xi )\mathrm{ln}(12\xi )`$
$`f_2(\xi )`$ $`=`$ $`2\gamma _E(\mathrm{ln}(1\xi )\mathrm{ln}(12\xi )){\displaystyle \frac{3}{2}}\mathrm{ln}(1\xi )`$ (49)
and

$`\begin{array}{cc}f_3(\xi )=0.804/(1\xi )2.32/(12\xi )\hfill & f_4(\xi )=0.779/(1\xi )^22.68/(12\xi )^2\hfill \\ f_5(\xi )=2.32/(1\xi )^35.00/(12\xi )^3\hfill & f_6(\xi )=6.12/(1\xi )^415.32/(12\xi )^4\hfill \\ f_7(\xi )=23.8/(1\xi )^561.12/(12\xi )^5\hfill & f_8(\xi )=120.9/(1\xi )^6305.52/(12\xi )^6,\hfill \end{array}`$ (53)
etc. Here the function $`f_1`$ resums the leading-logs (LL), $`f_2`$ resums the next-to-leading logs (NLL) and so on. The effect of the analytic continuation to the time-like region is taken into account through the summation over $`u`$ in (47). Note that we omitted a constant term in eq. (53). Such a constant does not necessarily exponentiate and it can be shifted into the “remainder function” $`\mathrm{ln}R(y)|_{\mathrm{non}\mathrm{log}}`$ in (2).

The first two functions $`f_1(\xi )`$ and $`f_2(\xi )`$ are already known from . When comparing $`f_2(\xi )`$ in (3.2) to  one should be careful concerning the definition of the coupling. Strictly speaking, our calculation is done in the large $`\beta _0`$ limit, so $`\beta _1`$ terms as well as other termsSuch terms originate  in the singular part of the NLO splitting function. which are sub-leading in $`\beta _0`$ are missing. However, these terms can be generated upon replacing the large $`\beta _0`$ coupling, which is one-loop, by a two-loop coupling in the “gluon bremsstrahlung” scheme  (see sec. 4) inside the integral (38). Such a replacement will be made in sec. 3.3 and then in the phenomenological analysis in sec. 5.

The functions $`f_k(\xi )`$, for $`k2`$ are singular at $`\xi =1/2`$ and at $`\xi =1`$; the first is associated with large-angle emission $`\nu \overline{\mathrm{\Lambda }}/Q1`$ and the second with collinear emission $`\nu \overline{\mathrm{\Lambda }}^2/Q^21`$. The singularity (at both limits) becomes stronger as the order $`k`$ increases, in particular $`f_k1/(12\xi )^{k2}`$ close to $`\xi =1/2`$ and $`f_k1/(1\xi )^{k2}`$ close to $`\xi =1`$, for any $`k3`$. Since soft gluon resummation has a universal structure, we expect that this singularity structure in $`\xi `$ will be common to different physical processes. An example is provided by a recent calculation of NNLL for deep-inelastic scattering and the Drell-Yan process . There, the singularity structure of $`f_k(\xi )`$ for $`k=1`$ through $`3`$ is indeed similar to (3.2) and (53). In particular $`f_3(\xi )`$ has simple poles.

As usual in perturbation theory the expansion (46) is divergent. In principle it can be given sense e.g. by truncation at the minimal term. Note that this truncation is cumbersome since the minimal term will correspond to different orders depending on the value of $`\nu `$. The standard procedure  amounts to truncating this expansion at NLL, so $`f_k(\xi )`$ for $`k3`$ are just neglected. To see whether this is legitimate let us now examine the numerical values of $`f_k(\xi )`$.

The functions $`f_k(\xi )`$ are plotted in fig. 2. It is clear from the plot that the expansion (46) is divergent. In particular, the ratios $`f_k(\xi )/f_{k+1}(\xi )`$ increase with $`k`$, reflecting factorial growth.

The actual effect of these increasing coefficients depends on the value of the expansion parameter, here $`\overline{A}(Q^2)=\beta _0a(Q^2)`$. $`\overline{A}(Q^2)`$ is defined in the “gluon bremsstrahlung” scheme  (the reason for this choice is explained in section 4). In fig. 3 we show the individual terms in (46) at two representative centre-of-mass energy scales: $`Q=91`$ and $`22\mathrm{GeV}`$.

First note that the NNLL ($`f_3`$) gives a significant contribution in all cases. It should be emphasised that the leading part of its contribution is taken into account in the standard NLL calculation in the process of matching with the NLO result. Next note that $`\mathrm{N}^3\mathrm{LL}`$ ($`f_4`$) also gives quite an important correction: $`8\%`$ at $`\mathrm{M}_\mathrm{Z}`$ and $`18\%`$ at $`22\mathrm{GeV}`$. A closer examination shows that the hierarchy between the contributions of sub-leading logs depends strongly on both the energy (the value of the coupling) and $`1/\nu `$. At sufficiently large $`1/\nu `$ values, i.e. far enough from the two-jet limit, the expansion (46) appears to be convergent at the first few orders. Consider for example the case $`1/\nu =0.075`$ at $`22\mathrm{GeV}`$: the contribution of $`f_4`$ (full line) is smaller than that of $`f_3`$ and $`f_5`$ (dotted) is yet smaller. But then the trend changes: $`f_6`$ (dashed) contributes slightly more than $`f_5`$ and the series starts diverging. In this particular example the contribution at the minimal term is about $`7.8\%`$. Thus a resummed perturbative calculation should not be expected to be better than that. For smaller $`1/\nu `$ values, and likewise at lower energies, the expansion diverges faster: the order at which the minimal term is reached becomes smaller while its contribution becomes larger. At some point, close to the distribution peak, the expansion becomes meaningless. Technically, this behaviour is a direct consequence of the enhanced singularity of the functions $`f_k`$ and their factorially increasing numerical coefficients. This is the way infrared renormalons show up in this framework. Physically it is well expected that the perturbative expansion would gradually become less predictive and eventually irrelevant as a quantitative evaluation of the cross-section, since it misses non-perturbative soft gluon contributions or, in other words, hadronization effects. We shall return to this issue in sec. 3.5. Until then, our goal is to complete the perturbative evaluation of the cross section. To achieve this we leave the expansion (46) and revert to Borel summation (38) which allows us, in rough terms, to resum the renormalons first, and then complete the summation over the logs.

### 3.3 Evaluation of the Borel sum

Our next task is to evaluate explicitly the Borel sum (38). An exact analytic integration is difficult. However, it is possible to make a systematic approximation of the Borel function and then perform the integration analytically.

We begin by writing the Borel function $`B_\nu (z)`$ as

$$B_\nu (z)=\frac{2}{z}V(2z)\left(\frac{2}{z}+\frac{1}{1z}+\frac{1}{2z}\right)V(z)$$
(54)
with

$$V(z)e^{z\mathrm{ln}\nu }\mathrm{\Gamma }(z)+\frac{1}{z}.$$
(55)
Since $`V(z)`$ has only simple poles at positive integer $`z`$ values, it is natural to approximate it by a rational function. Clearly, we would like the approximation to coincide with the function at least at the first few orders in the expansion around $`z=0`$. One possibility is then to use Padé approximants based on some truncated sum. However, since we know the singularity structure, it is better to use an approximation of the formA similar approximation for the Borel function was used in  in a different application. Comparison with Padé approximants showed consistent results.

$$V(z)\underset{p=1}{\overset{p_{\mathrm{max}}}{}}\frac{r_p}{pz},$$
(56)
where $`r_p`$ are set such that the expansion of this approximant coincides with that of (55) up to order $`z^{p_{\mathrm{max}}1}`$. This means that $`r_p`$, for every $`1pp_{\mathrm{max}}`$, is a polynomial of order $`p_{\mathrm{max}}`$ in $`\mathrm{ln}\nu `$,

$$r_p=\underset{q=0}{\overset{p_{\mathrm{max}}}{}}\rho _{p,q}\left(\mathrm{ln}\nu \right)^q$$
(57)
where $`\rho _{p,q}`$ are (irrational) numerical coefficients.

Next, $`B_\nu (z)`$ can be written as a sum of simple and double poles,

$$B_\nu (z)\underset{n=1}{\overset{p_{\mathrm{max}}}{}}\frac{\overline{r}_{n/2}^{\mathrm{la}}}{n/2z}+\underset{p=1}{\overset{p_{\mathrm{max}}}{}}\frac{\overline{r}_p^{\mathrm{col}}}{pz}\frac{r_1}{\left(1z\right)^2}\frac{r_2}{\left(2z\right)^2}$$
(58)
where $`\overline{r}_p^{\mathrm{la}}`$ and $`\overline{r}_p^{\mathrm{col}}`$ are linear combinations of $`\left\{r_p\right\}_{p=1}^{p_{\mathrm{max}}}`$, corresponding to the large-angle (first term) and collinear (second term) parts of (54), respectively. Simple poles appear at all positive integers and half integers $`1/2pp_{\mathrm{max}}/2`$ due to the large-angle part, and at all integers $`3pp_{\mathrm{max}}`$ due to the collinear part. The double poles at $`z=1`$ and $`z=2`$ appear due to the collinear part. As an example, consider the case of $`p_{\mathrm{max}}=3`$ where one has

$$B_\nu (z)\frac{2r_1}{\frac{1}{2}z}+\frac{3r_1\frac{1}{2}r_3}{1z}+\frac{\frac{2}{3}r_3}{\frac{3}{2}z}+\frac{r_1r_3}{2z}+\frac{\frac{5}{6}r_3}{3z}\frac{r_1}{(1z)^2}\frac{r_2}{(2z)^2}.$$

Finally, given a specific coupling and the function (58) the Borel integral (38) can be evaluated analytically. The simple pole integrals are of the form

$`I_p(s)`$ $`=`$ $`{\displaystyle _0^{\mathrm{}}}𝑑z\mathrm{exp}\left(z\mathrm{ln}Q^2/\mathrm{\Lambda }^2\right){\displaystyle \frac{1}{pz}}{\displaystyle \frac{\mathrm{sin}\pi z}{\pi z}}\overline{A}_B(z)`$
$`=`$ $`{\displaystyle _0^{\mathrm{}}}𝑑z\mathrm{exp}\left(zs\right){\displaystyle \frac{1}{pz}}{\displaystyle \frac{\mathrm{sin}\pi z}{\pi z}}A_B(z)`$
where $`A_B(z)`$ is the Borel transform of the $`\overline{\mathrm{MS}}`$ coupling and $`s\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2`$. Regarding $`p`$ as a complex number, $`I_p(s)`$ can be written in terms of $`\stackrel{~}{I}_p(s)`$,

$$\stackrel{~}{I}_p(s)=_0^{\mathrm{}}𝑑z\mathrm{exp}\left(zs\right)\frac{1}{pz}A_B(z),$$
(60)
by analytically continuing in $`s`$, $`\overline{I}_p(s)=\frac{1}{2\pi i}\left(\stackrel{~}{I}_p(si\pi )\stackrel{~}{I}_p(s+i\pi )\right)`$, and

$$I_p(s)=\frac{1}{p}\left[\overline{I}_p(s)\underset{p=0^{}}{lim}\overline{I}_p(s)\right].$$
(61)
The double pole integrals

$$I_p^{(1)}(s)=_0^{\mathrm{}}𝑑z\mathrm{exp}\left(zs\right)\frac{1}{\left(pz\right)^2}\frac{\mathrm{sin}\pi z}{\pi z}A_B(z)$$
(62)
can be simply obtained using $`I_p^{(1)}(s)=\frac{d}{dp}I_p(s)`$.

In case of infrared renormalons in space-like observables $`\stackrel{~}{I}_p(s)`$ should be evaluated for real $`s`$ and $`p`$ where it is ill-defined. Then Cauchy principal value can be defined taking the real part of (60). Our definition for the time-like renormalon integral is motivated by this regularisation. Here there is no need to take the real part: thanks to the analytic continuation procedure, (60) is evaluated above and below the cut in a symmetric way (the relevant complex parameter is $`ps`$), yielding a real value for $`I_p(s)`$ and $`I_p^{(1)}(s)`$. The ambiguity of these quantities is directly related to the ambiguity in $`\stackrel{~}{I}_p(s)`$. The definition of the Borel-sum through the space-like renormalon integral also allows to introduce a cut-off regularisation .

$`\stackrel{~}{I}_p(s)`$ can be written in terms of known analytic functions. The simplest case is the one-loop running coupling, where $`A_B(z)1`$. Here we have

$$\stackrel{~}{I}_p(s)|_{\mathrm{one}\mathrm{loop}}=\mathrm{Ei}(1,ps)e^{ps}.$$
(63)
The more difficult two-loop case was dealt with in  (see sec. 3.4 there) using results from . It is useful to express the coupling using the Lambert W function ,

$$\overline{A}(Q^2)|_{\mathrm{two}\mathrm{loop}}=\frac{1}{\delta }\frac{1}{1+w(s)}$$
(64)
where $`\delta \beta _1/\beta _0^2`$ and <sup>§</sup><sup>§</sup>§The analytic continuation in the calculation of $`\overline{I}_p(s)`$ should be done carefully: one should use the correct branch $`n`$ of the Lambert W function, $`W_n`$. In practice this is simple since for $`N_c=3`$, $`N_f6`$ the time-like axis is always contained in the $`n=\pm 1`$ branches (see  for more details). Also note that with the current definition of the two-loop coupling the Landau branch point is at $`Q^2=\overline{\mathrm{\Lambda }}^2`$.  $`w(s)W_1\left(e^{s/\delta 1}\right)`$. In these variables the two-loop renormalon integral (60) is simply

$$\stackrel{~}{I}_p(s)|_{\mathrm{two}\mathrm{loop}}=\left[w(s)p\delta \right]^{p\delta }e^{w(s)p\delta }\mathrm{\Gamma }(p\delta ,w(s)p\delta ).$$
(65)
Altogether the Borel integral (38) is given by

$$e^{\nu t}_{\text{SDG}}\frac{C_F}{\beta _0}\left[\underset{n=1}{\overset{p_{\mathrm{max}}}{}}\overline{r}_{n/2}^{\mathrm{la}}I_{n/2}(s)+\underset{p=1}{\overset{p_{\mathrm{max}}}{}}\overline{r}_p^{\mathrm{col}}I_p(s)r_1I_1^{(1)}(s)r_2I_2^{(1)}(s)\right],$$
(66)
where the (polynomial) dependence on $`\mathrm{ln}\nu `$ is fully contained in the coefficients $`\overline{r}_p`$ and $`r_p`$, given by eqs. (55) through (57), while the dependence on the external energy scale $`s\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2`$ is fully contained in the renormalon integrals $`I_p(s)`$. Eventually, the log-resummation order $`p_{\mathrm{max}}`$ should be set large enough, so that higher orders can be safely neglected. In sec. 5 we show that in practice $`p_{\mathrm{max}}6`$ is enough to make the result stable.

### 3.4 Evaluation of the inverse Laplace transform

The last stage of the calculation is to evaluate the inverse Laplace transform (37),

$$R(Q^2,t)|_{\text{res}}=_𝒞\frac{d\nu }{2\pi i\nu }e^{\nu t}\mathrm{exp}\left\{𝒮_{\text{PT}}(Q^2,\mathrm{ln}\nu )\right\}$$
(67)
with $`𝒮_{\text{PT}}(Q^2,\mathrm{ln}\nu )e^{\nu t}_{\text{SDG}}`$. We apply the same technique used in  (see section 5 there). However, at a difference with the latter, we must keep track of sub-leading logs so we refrain from making any approximation. Defining $`u=\nu t`$ and expanding $`𝒮(Q^2,\mathrm{ln}\nu )=𝒮(Q^2,L+\mathrm{ln}u)`$ (where $`L=\mathrm{ln}\frac{1}{t}`$) in powers of $`\mathrm{ln}u`$, one has

$`R(Q^2,t)|_{\text{res}}`$ $`=`$ $`{\displaystyle _𝒞}{\displaystyle \frac{du}{2\pi iu}}\mathrm{exp}\left\{u+𝒮(Q^2,L+\mathrm{ln}u)\right\}`$
$`=`$ $`e^{𝒮(Q^2,L)}{\displaystyle _𝒞}{\displaystyle \frac{du}{2\pi iu}}\mathrm{exp}\{u+d_1(Q^2,L)\mathrm{ln}u+{\displaystyle \underset{k=2}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{k!}}d_k(Q^2,L)(\mathrm{ln}u)^k)\}`$
where $`d_k(Q^2,L)\frac{d^k}{dL^k}𝒮(Q^2,L)`$. Finally, we calculate $`R(Q^2,t)|_{\text{res}}`$ by expanding

$$\mathrm{exp}\{\underset{k=2}{\overset{\mathrm{}}{}}\frac{1}{k!}d_k(Q^2,L)(\mathrm{ln}u)^k)\}$$

in powers of $`\mathrm{ln}u`$ and using the integral

$$_𝒞\frac{du}{2\pi iu}(\mathrm{ln}u)^n\mathrm{exp}\left\{u+g\mathrm{ln}u\right\}=\frac{d^n}{dg^n}\frac{1}{\mathrm{\Gamma }(1g)}.$$
(69)
Our final expression for $`R(Q^2,t)|_{\text{res}}`$ includes several special functions: derivatives of the $`\mathrm{\Gamma }`$ function from (69) and the expressions for $`I_p(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2)`$, written in terms of the exponential integral function in the one-loop case (63), or the incomplete $`\mathrm{\Gamma }`$ function and the Lambert W function in the two-loop case (65). We therefore chose to evaluate it using the computer algebra program Maple.

### 3.5 Power-Corrections

Since hadronization is not described by perturbation theory it is a priori clear that there are non-perturbative corrections to the resummed cross-section. In this section we study the structure of these corrections. We already saw that the perturbative sum is, in general, ambiguous due to infrared renormalons. Having used a specific regularization of the renormalon singularities we are bound to address the question of the related power-suppressed ambiguity. The appearance of infrared renormalons does not only signal the existence of non-perturbative power-corrections, but also suggests their functional dependence on the energy and the event-shape variable. One should not expect, of course, that the entire non-perturbative contribution would be deducible from a perturbative result. Appropriate non-perturbative parameters will be introduced and eventually fixed by the data. Upon including these non-perturbative corrections the result becomes well-defined to power accuracy.

If the Borel sum (38) is supplemented by a non-perturbative contribution, eq. (67) becomes

$$R(Q^2,t)=_𝒞\frac{d\nu }{2\pi i\nu }e^{\nu t}\mathrm{exp}\left\{𝒮_{\text{PT}}(Q^2,\mathrm{ln}\nu )\right\}\mathrm{exp}\left\{𝒮_{\text{NP}}(Q^2,\mathrm{ln}\nu )\right\}.$$
(70)
In sec. 3.4 we evaluated the inverse Laplace integral ignoring the non-perturbative part. Thanks to the factorized form of (70) non-perturbative corrections can be included by convoluting  the perturbative result (67) with a non-perturbative function:

$$R(Q^2,t)=_0^t𝑑\stackrel{~}{t}R(Q^2,t\stackrel{~}{t})|_{\text{res}}f_{\text{NP}}(Q^2,\stackrel{~}{t})R(Q^2,t)|_{\text{res}}f_{\text{NP}}(Q^2,t)$$
(71)
where

$$f_{\text{NP}}(Q^2,t)=_𝒞\frac{d\nu }{2\pi i}e^{\nu t}\mathrm{exp}\left\{𝒮_{\text{NP}}(Q^2,\mathrm{ln}\nu )\right\}.$$
(72)
In order to learn about the functional form of $`𝒮_{\text{NP}}(Q^2,\mathrm{ln}\nu )`$, we shall now examine the renormalon structure of $`𝒮_{\text{PT}}(Q^2,\mathrm{ln}\nu )`$. Considering the integrand of (38) with $`B_\nu (z)`$ of eq. (42) (see description following the equation), we identify potential singularities at all positive integers and half integers. On the other hand, the $`\mathrm{sin}\pi z`$ factor which originates in the analytic continuation to the time-like region introduces an explicit zero at all integer $`z`$ values. It is important to realize that the complete cancellation of some of the singularities may be an artifact of the large $`\beta _0`$ approximation. A similar situation occurs in any time-like observable, e.g. in the total cross-section . There, having an operator product expansion, one can relate the singularity strength of the BorelNote that the $`p\beta _1/\beta _0^2`$ term appears in the standard Borel transform but not in the scheme-invariant one . pole $`1/\left(pz\right)^{1+\kappa }`$ to the anomalous dimension of the corresponding operator(s) $`𝒪_p`$ of dimension $`2p`$: $`\kappa =p\beta _1/\beta _0^2+\gamma _0/\beta _0`$, where $`\left[d/d\mathrm{ln}\mu ^2+\gamma \right]𝒪_p=0`$ with $`\gamma =\gamma _0(\alpha _s/\pi )+\mathrm{}`$. Thus, in general, a large-$`\beta _0`$ calculation does not yield the correct singularity strength: not only the residue, but also $`\kappa `$ is expected to be modified when sub-leading terms in $`\beta _0`$ are included. Being aware of this fact (and having no operator product expansion to supply additional information) it is reasonable to allow for modifications of order $`1`$ in the singularities of $`B_\nu (z)`$,

$$b(z)=\frac{r_p}{\left(pz\right)^m}\stackrel{~}{b}(z)=\frac{\stackrel{~}{r}_p}{\left(pz\right)^{m+\kappa _p}}$$
(73)
where we assume $`\left|\kappa _p\right|<1`$. Thus, instead of simple and double poles ($`m=1,2`$), the Borel integrand has cuts which survive the suppression by the analytic continuation and introduce some ambiguity into $`𝒮_{\text{PT}}(Q^2,\mathrm{ln}\nu )`$. Note that, we will not consider the possibility of a modified $`\nu `$ dependence of the exact Borel function compared to the calculated one. We allow only a modification of the overall coefficient and the singularity strength.

Consider first a scheme-invariant Borel integral in the space-like region,

$$\stackrel{~}{I}(s)=_0^{\mathrm{}}e^{zs}\stackrel{~}{b}(z)𝑑z,$$
(74)
where the Borel function behaves as (73) with $`m=1`$ near $`z=p`$. The Borel integral can be written as a sum of two terms,

$`\stackrel{~}{I}(s)`$ $`=`$ $`{\displaystyle _0^p}𝑑ze^{zs}{\displaystyle \frac{1}{\left(pz\right)^{1+\kappa }}}+{\displaystyle _p^{\mathrm{}}}𝑑ze^{zs}{\displaystyle \frac{1}{\left(pz\right)^{1+\kappa }}}`$
$`=`$ $`(s)^\kappa \gamma (\kappa ,sp)e^{sp}(s)^\kappa \mathrm{\Gamma }(\kappa )e^{sp}=(s)^\kappa \mathrm{\Gamma }(\kappa ,sp)e^{sp}.`$
Next we use the fact that the ambiguity $`\mathrm{\Delta }\stackrel{~}{I}(s)`$ is purely imaginary and originates uniquely in the second integral where $`zp`$, i.e. in the complete $`\mathrm{\Gamma }`$ function partA similar mathematical analysis appears in  in the context of a renormalon integral with a two-loop running coupling.. Taking $`1/\pi `$ times the imaginary part, we obtain

$$\mathrm{\Delta }\stackrel{~}{I}(s)=\frac{\mathrm{sin}\left(\pi \kappa \right)}{\pi }s^\kappa \mathrm{\Gamma }(\kappa )e^{sp}=\frac{1}{\mathrm{\Gamma }(1+\kappa )}s^\kappa e^{sp}.$$
(75)
The ambiguity for the $`m=2`$ case can be calculated using $`\stackrel{~}{I}^{(1)}(s)=\frac{d}{dp}\stackrel{~}{I}(s)`$: one finds that $`\mathrm{\Delta }\stackrel{~}{I}^{(1)}(s)=\mathrm{\Delta }\stackrel{~}{I}(s)s`$. Substituting $`s=\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2`$ we see that the ambiguity $`\mathrm{\Delta }\stackrel{~}{I}(s)`$ scales as a power $`(Q^2/\overline{\mathrm{\Lambda }}^2)^p`$, up to a logarithmic factor, $`\left(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2\right)^\kappa `$, which remains undetermined as long as $`\kappa `$ is unknown. Next, the ambiguity of the corresponding Borel integral in the time-like region

$$I(s)=_0^{\mathrm{}}e^{zs}\stackrel{~}{b}(z)\frac{\mathrm{sin}\pi z}{\pi z}𝑑z,$$
(76)
can be deduced by analytically continuing (75) in $`Q^2Q^2`$, i.e. $`ss\pm i\pi `$, and taking ($`\frac{1}{\pi p}`$ times) the imaginary part,

$`\mathrm{\Delta }I(s)`$ $`=`$ $`{\displaystyle \frac{1}{\pi p}}{\displaystyle \frac{\left(s^2+\pi ^2\right)^{\kappa /2}}{\mathrm{\Gamma }(1+\kappa )}}e^{sp}\left[\mathrm{cos}\left(\kappa \mathrm{arctan}{\displaystyle \frac{\pi }{s}}\right)\mathrm{sin}(\pi p)\mathrm{sin}\left(\kappa \mathrm{arctan}{\displaystyle \frac{\pi }{s}}\right)\mathrm{cos}(\pi p)\right]`$ (77)
$``$ $`{\displaystyle \frac{1}{\pi p}}{\displaystyle \frac{1}{\mathrm{\Gamma }(1+\kappa )}}s^\kappa e^{sp}\left[\mathrm{sin}(\pi p){\displaystyle \frac{\kappa \pi }{s}}\mathrm{cos}(\pi p)\right],`$
where the expression in the second line is obtained expanding in $`1/s`$ and neglecting $`1/s^2`$ terms and higher. Similarly, we obtain

$`\mathrm{\Delta }I^{(1)}(s){\displaystyle \frac{1}{\pi p}}{\displaystyle \frac{1}{\mathrm{\Gamma }(1+\kappa )}}s^\kappa e^{sp}\left[s\mathrm{sin}(\pi p)(1+\kappa )\pi \mathrm{cos}(\pi p)\right].`$ (78)
Given the general singularity structure of the SDG Borel function $`B_\nu (z)`$ around $`z=p`$,

$$b(z)=\frac{r_p^{(1)}}{\left(pz\right)^2}+\frac{r_p}{pz},$$
(79)
our result for the ambiguity ((77) and (78)) suggests the following assignments of non-perturbative contributions,

$$I_{\text{NP}}(s)=\frac{\rho _{2p}}{\pi p}\frac{e^{sp}s^{\kappa _p}}{\mathrm{\Gamma }(1+\kappa _p)}\left[\left(r_p^{(1)}s+r_p\right)\mathrm{sin}\left(\pi p\right)\left(\left(1+\kappa _p\right)r_p^{(1)}+\kappa _pr_p/s\right)\pi \mathrm{cos}\left(\pi p\right)\right],$$
(80)
where $`s=\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2`$ and $`\rho _{2p}`$ are dimensionless phenomenological parameters to be fixed, eventually, by the data. On general grounds we expect that $`\rho _{2p}𝒪(1)`$. It is important to keep in mind that the actual values of $`\rho _{2p}`$ depend on the regularization prescription used for the Borel sum. For example, when the perturbative sum is defined using a cut-off (on the Euclidean momentum), $`\rho _{2p}`$ can be interpreted as moments of an “infrared finite effective coupling” $`\overline{A}_{\text{NP}}(k^2)`$ , which is assumed to coincide with the perturbative coupling $`\overline{A}_{\text{PT}}(k^2)`$ above some scale $`\mu _I`$ ($`\overline{\mathrm{\Lambda }}\mu _IQ`$),

$$\rho _{2p}^{\mu _I}=\left(\frac{\mu _I^2}{\overline{\mathrm{\Lambda }}^2}\right)^pm_p^{\text{NP}}(\mu _I^2)$$
(81)
where

$$m_p^{\text{NP}}(\mu _I^2)_0^{\mu _I^2}p\frac{dk^2}{k^2}\left(\frac{k^2}{\mu _I^2}\right)^p\overline{A}_{\text{NP}}(k^2).$$
(82)
Here we define the perturbative sum $`\stackrel{`}{\mathrm{a}}`$ la Borel. In this case $`\rho _{2p}`$ is

$$\rho _{2p}=\left(\frac{\mu _I^2}{\overline{\mathrm{\Lambda }}^2}\right)^p\left[m_p^{\text{NP}}(\mu _I^2)m_p^{\text{PT}}(\mu _I^2)\right]$$
(83)
where

$$m_p^{\text{PT}}(\mu _I^2)\mathrm{P}.\mathrm{V}._0^{\mu _I^2}p\frac{dk^2}{k^2}\left(\frac{k^2}{\mu _I^2}\right)^p\overline{A}_{\text{PT}}(k^2),$$
(84)
and P.V. indicates the principal value prescription. Note that in this case $`\rho _{2p}`$ are independent of $`\mu _I`$. A detailed comparison of the two regularization prescriptions can be found in . The explicit expressions for $`m_p^{\text{PT}}(\mu _I^2)`$ in the one-loop and two-loop cases can be obtained by substituting $`s\mathrm{ln}\mu _I^2/\overline{\mathrm{\Lambda }}^2`$ in eqs. (63) and (65) (taking the real part and multiplying by $`p`$), respectively.

Let us now examine in more detail the singularities in (42). First note that double poles at half integer values do not occur. Thus the term $`r_p^{(1)}s\mathrm{sin}\left(\pi p\right)`$ in (80) is actually not relevant. Next, to determine the residues $`r_p`$ of $`B_\nu (z)`$ we use the fact that<sup>\**</sup><sup>\**</sup>\**Since our interest is in large $`\nu `$ values, corresponding to small $`t`$, we can replace $`\gamma (z,\nu )`$ by $`\mathrm{\Gamma }(z)`$, which is a good approximation at any $`\nu >1`$. $`\mathrm{Residue}\left\{\mathrm{\Gamma }(z)\right\}=(1)^{p+1}/p!`$, and therefore $`\mathrm{Residue}\left\{\mathrm{\Gamma }(2z)\right\}=(1)^{2p+1}/\left(2(2p)!\right)`$. Consider the first term in (42) which is related to large-angle emission. Here $`r_p^{la}=\frac{1}{p(2p)!}\nu ^{2p}(1)^{2p}`$ and there are no double poles, so

$`I_p^{la}={\displaystyle \frac{\rho _{2p}}{\pi p^2(2p)!\mathrm{\Gamma }(1+\kappa _p)}}\left({\displaystyle \frac{\nu \overline{\mathrm{\Lambda }}}{Q}}\right)^{2p}\left[\left(\mathrm{ln}{\displaystyle \frac{Q^2}{\overline{\mathrm{\Lambda }}^2}}\right)^{\kappa _p}\mathrm{sin}\left(\pi p\right)+\kappa _p\pi \left(\mathrm{ln}{\displaystyle \frac{Q^2}{\overline{\mathrm{\Lambda }}^2}}\right)^{\kappa _p1}\mathrm{cos}\left(\pi p\right)\right].`$ (85)
Summing over $`p`$ this gives

$$e^{\nu t}_{\text{NP}}^{\mathrm{large}\mathrm{angle}}=\underset{p=\frac{1}{2},1,\frac{3}{2},2,\frac{5}{2},\mathrm{}}{}\lambda _{2p}\frac{1}{(2p)!}\left(\frac{\nu \overline{\mathrm{\Lambda }}}{Q}\right)^{2p}=\underset{n=1}{\overset{\mathrm{}}{}}\lambda _n\frac{1}{n!}\left(\frac{\nu \overline{\mathrm{\Lambda }}}{Q}\right)^n$$
(86)
where we defined

$`2p\mathrm{odd}`$ $`\lambda _{2p}={\displaystyle \frac{C_F}{\beta _0}}{\displaystyle \frac{\rho _{2p}}{p^2}}{\displaystyle \frac{1}{\mathrm{\Gamma }(1+\kappa _p)}}\left(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2\right)^{\kappa _p}{\displaystyle \frac{\mathrm{sin}\left(\pi p\right)}{\pi }}`$
$`2p\mathrm{even}`$ $`\lambda _{2p}={\displaystyle \frac{C_F}{\beta _0}}{\displaystyle \frac{\rho _{2p}}{p^2}}\kappa _p{\displaystyle \frac{1}{\mathrm{\Gamma }(1+\kappa _p)}}\left(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2\right)^{\kappa _p1}\mathrm{cos}\left(\pi p\right).`$ (87)
The second term in (42), associated with collinear emission, is treated in a similar manner. Note that here there are also double poles at $`p=1`$ and $`2`$ and consequently the residues have a more complicated dependence on $`\nu `$. The result is

$`e^{\nu t}_{\text{NP}}^{\mathrm{col}}`$ $`=`$ $`\lambda _2^{(1)}{\displaystyle \frac{\nu \overline{\mathrm{\Lambda }}^2}{Q^2}}\lambda _4^{(1)}{\displaystyle \frac{\nu ^2\overline{\mathrm{\Lambda }}^4}{Q^4}}`$ (88)
$``$ $`\left[\left(4\gamma _E\mathrm{ln}\nu \right)\nu 1\right]\lambda _2{\displaystyle \frac{\overline{\mathrm{\Lambda }}^2}{Q^2}}\left[\left({\displaystyle \frac{3}{2}}+\gamma _E+\mathrm{ln}\nu \right)\nu ^21\right]\lambda _4{\displaystyle \frac{\overline{\mathrm{\Lambda }}^4}{Q^4}}`$
$`{\displaystyle \underset{p3}{\overset{\mathrm{}}{}}}\lambda _{2p}\left({\displaystyle \frac{2}{p}}+{\displaystyle \frac{1}{1p}}+{\displaystyle \frac{1}{2p}}\right){\displaystyle \frac{(1)^{p+1}}{(p1)!}}\left({\displaystyle \frac{\nu \overline{\mathrm{\Lambda }}^2}{Q^2}}\right)^p,`$
where the first line corresponds to the double poles at $`p=1`$ and $`2`$, with

$$\lambda _{2p}^{(1)}=\frac{C_F}{\beta _0}\frac{\rho _{2p}}{p}\frac{(1+\kappa _p)}{\mathrm{\Gamma }(1+\kappa _p)}\left(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2\right)^{\kappa _p};$$

the second line corresponds to the simple poles at $`p=1`$ and $`2`$, and the last line to the remaining simple poles at integers $`p3`$. $`\lambda _{2p}`$ are given in (3.5).

Finally, the total non-perturbative contribution in (70) is given by

$$𝒮_{\text{NP}}(Q^2,\mathrm{ln}\nu )e^{\nu t}_{\text{NP}}^{\mathrm{large}\mathrm{angle}}+e^{\nu t}_{\text{NP}}^{\mathrm{col}}.$$
(89)
Note that in both (86) and (88) we use the same non-perturbative parameters $`\rho _{2p}`$. One could introduce a more general parametrisation of power-corrections by defining separate parameters for the large-angle and collinear regions of phase-space. This, however, seems redundant from the current point of view: compensating the ambiguity in the perturbative sum does not require separate parameters. Moreover, the delicate relation between the large-angle and collinear limits of phase-space at the perturbative level suggests that the two are correlated in the same way also at the non-perturbative level.

It is now straightforward to see the relation with previously suggested models for power-corrections . As mentioned in the introduction, it was shown  that the primary non-perturbative effect in the region $`t\overline{\mathrm{\Lambda }}/Q`$ is a shift of the resummed distribution to larger values of $`t`$:

$$\frac{1}{\sigma }\frac{d\sigma }{dt}|_{\text{res}}(Q^2,t)\frac{1}{\sigma }\frac{d\sigma }{dt}|_{\text{res}}(t\lambda _1\overline{\mathrm{\Lambda }}/Q).$$
(90)
This is indeed the effect if one keeps only the first term in the sum in (86): then the only non-perturbative contribution to $`𝒮_{\text{NP}}(Q^2,\mathrm{ln}\nu )`$ in (70) is $`\lambda _1\nu \overline{\mathrm{\Lambda }}/Q`$, and

$$f_{\mathrm{la}}=\delta (t\lambda _1\overline{\mathrm{\Lambda }}/Q).$$
(91)
Assuming that the average thrust is dominated by the two-jet region, the corresponding correction is $`tt+\lambda _1\overline{\mathrm{\Lambda }}/Q`$, namely it is the same parameter that controls the shift of the distribution and the leading power-correction to $`t`$. This $`1/Q`$ correction to the average thrust was discussed extensively in the literature .

As explained in , the shift (90) is no longer adequate when $`t\overline{\mathrm{\Lambda }}/Q`$. The shape-function  allows parametrisation of all the terms in (86) making the following identification:

$$_0^{\mathrm{}}𝑑\stackrel{~}{t}\mathrm{exp}(\nu \stackrel{~}{t})(Q/\overline{\mathrm{\Lambda }})f_{\mathrm{la}}(\stackrel{~}{t}Q/\overline{\mathrm{\Lambda }})=\mathrm{exp}\left\{e^{\nu t}_{\text{NP}}^{\mathrm{large}\mathrm{angle}}\right\}.$$
(92)
The single argument function $`(Q/\overline{\mathrm{\Lambda }})f_{\mathrm{la}}(tQ/\overline{\mathrm{\Lambda }})`$ is therefore a particular example of the two-argument non-perturbative function $`f_{\text{NP}}(Q^2,t)`$ defined in eqs. (71) and (72). Similarly to , the moments of the large-angle shape-function $`\sigma _m_0^{\mathrm{}}f_{\mathrm{la}}(ϵ)ϵ^m𝑑ϵ`$ are related to the parameters $`\lambda _{2p}`$. Writing (92) as

$$_0^{\mathrm{}}\mathrm{exp}\left\{\frac{\nu ϵ\overline{\mathrm{\Lambda }}}{Q}\right\}f_{\mathrm{la}}(ϵ)𝑑ϵ=\mathrm{exp}\left\{\underset{n=1}{\overset{\mathrm{}}{}}\lambda _n\frac{1}{n!}\left(\frac{\nu \overline{\mathrm{\Lambda }}}{Q}\right)^n\right\}$$
(93)
and expanding the exponents on both sides one has:

$$\sigma _0=1,\sigma _1=\lambda _1,\sigma _2=\lambda _1^2\lambda _2,\sigma _3=\lambda _1^33\lambda _1\lambda _2+\lambda _3.$$
(94)
The function $`f_{\mathrm{la}}(ϵ)`$, like the parameters $`\lambda _{2p}`$, depends on the regularization prescription of the perturbative calculation. In the absence of a non-perturbative calculation (ref.  gives a field theoretic definition) the shape-function and thus $`\lambda _n`$ can only be determined by fitting the data.

Finally, to deal with the collinear power-corrections, one should include (88). Here, the most important corrections arise from the terms that are enhanced by $`\mathrm{ln}\nu `$. Note that the remaining terms can be simply recast in the form

$$_0^{\mathrm{}}𝑑\stackrel{~}{t}\mathrm{exp}(\nu \stackrel{~}{t})(Q^2/\overline{\mathrm{\Lambda }}^2)f_{\mathrm{col}}^{\mathrm{sub}\mathrm{leading}}(\stackrel{~}{t}Q^2/\overline{\mathrm{\Lambda }}^2)$$
(95)
similarly to the shape-function of (92). These corrections become important only at $`t\overline{\mathrm{\Lambda }}^2/Q^2`$ and can be safely neglected in the region $`t\overline{\mathrm{\Lambda }}/Q`$. To treat the terms that are enhanced by $`\mathrm{ln}\nu `$ we first expand the exponent,

$`\mathrm{exp}\left\{e^{\nu t}_{\text{NP}}^{\mathrm{col}}\right\}`$ $`=`$ $`1+\lambda _2{\displaystyle \frac{\nu \overline{\mathrm{\Lambda }}^2}{Q^2}}\left(\mathrm{ln}\nu +\gamma _E\right)+\left[{\displaystyle \frac{\lambda _2^2}{2}}\left(\mathrm{ln}\nu +\gamma _E\right)^2\lambda _4\left(\mathrm{ln}\nu +\gamma _E\right)\right]{\displaystyle \frac{\nu ^2\overline{\mathrm{\Lambda }}^4}{Q^4}}`$ (96)
$`+`$ $`\left[{\displaystyle \frac{\lambda _2^3}{6}}\left(\mathrm{ln}\nu +\gamma _E\right)^3\lambda _2^2\lambda _4\left(\mathrm{ln}\nu +\gamma _E\right)^2\right]{\displaystyle \frac{\nu ^3\overline{\mathrm{\Lambda }}^6}{Q^6}}+\mathrm{},`$
where we kept the sub-leading $`\gamma _E`$ constants for convenience (see below). Next we use the fact that

$$\mathrm{\Gamma }(ϵ)\nu ^ϵ=_0^{\mathrm{}}𝑑t\frac{1}{t}t^ϵe^{\nu t}$$

to write the following convolution integral in the “+” prescription, with a generic test function $`h(t)`$ whose inverse Laplace transform is $`H(\nu )`$,

$`\left[\mathrm{\Gamma }(ϵ)\nu ^ϵ{\displaystyle \frac{1}{ϵ}}\right]H(\nu )`$ $`=`$ $`{\displaystyle _0^{\mathrm{}}}𝑑te^{\nu t}{\displaystyle _0^t}𝑑\stackrel{~}{t}\left[{\displaystyle \frac{1}{\stackrel{~}{t}}}e^{ϵ\mathrm{ln}\stackrel{~}{t}}\right]_+h(t\stackrel{~}{t})`$
$`=`$ $`{\displaystyle _0^{\mathrm{}}}𝑑te^{\nu t}{\displaystyle _0^t}𝑑\stackrel{~}{t}{\displaystyle \frac{1}{\stackrel{~}{t}}}e^{ϵ\mathrm{ln}\stackrel{~}{t}}\left[h(t\stackrel{~}{t})h(t)\right]`$
Now one can extract the inverse Laplace transform of the terms in (96), in the distribution sense, by expanding both sides of eq. (3.5) in powers of $`ϵ`$. In particular, one has

$$_0^{\mathrm{}}𝑑t\left[\frac{1}{t}\right]_+e^{\nu t}=(\mathrm{ln}\nu +\gamma _E)_0^{\mathrm{}}𝑑t\left[\frac{\mathrm{ln}t}{t}\right]_+e^{\nu t}=\frac{1}{2}(\mathrm{ln}\nu +\gamma _E)^2+\frac{\pi ^2}{12}.$$

Thus,

$`f_{\mathrm{col}}(Q^2,t)`$ $`=`$ $`{\displaystyle _𝒞}{\displaystyle \frac{d\nu }{2\pi i}}e^{\nu t}\mathrm{exp}\left\{e^{\nu t}_{\text{NP}}^{\mathrm{col}}\right\}`$
$`=`$ $`\delta (t)\lambda _2\left[{\displaystyle \frac{\overline{\mathrm{\Lambda }}^2}{Q^2t}}\right]_+\delta ^{(1)}(t)+\left[\lambda _2^2{\displaystyle \frac{\overline{\mathrm{\Lambda }}^4}{Q^4}}\left({\displaystyle \frac{\mathrm{ln}t}{t}}{\displaystyle \frac{\pi ^2}{12}}\delta (t)\right)+\lambda _4{\displaystyle \frac{\overline{\mathrm{\Lambda }}^4}{Q^4}}{\displaystyle \frac{1}{t}}\right]_+\delta ^{(2)}(t)+\mathrm{}`$
$`=`$ $`\left[1+\lambda _2{\displaystyle \frac{\overline{\mathrm{\Lambda }}^2}{Q^2t^2}}\left(2\lambda _2^2\mathrm{ln}{\displaystyle \frac{1}{t}}+𝒪\left(1\right)\right){\displaystyle \frac{\overline{\mathrm{\Lambda }}^4}{Q^4t^3}}+𝒪\left({\displaystyle \frac{\overline{\mathrm{\Lambda }}^6}{Q^6t^4}}\right)\right]\delta (t),`$
where $`\delta ^{(n)}(t)`$ is the $`n`$-th derivative of the $`\delta `$-function<sup>††</sup><sup>††</sup>††In the last line, the appropriate integration prescription at $`t0`$ is understood.. Clearly, the leading correction $`\frac{\lambda _2}{Q^2t^2}`$ is as important around $`t\overline{\mathrm{\Lambda }}/Q`$ as the large-angle corrections contained in the shape-function (92). The sub-leading terms in (3.5) are less important and can be ignored in a first approximation.

In conclusion, the most important power correction at $`t\overline{\mathrm{\Lambda }}/Q`$ is the one associated with the simple renormalon pole at $`z=\frac{1}{2}`$. Its exponentiation amounts to a shift of the resummed distribution by $`\lambda _1\overline{\mathrm{\Lambda }}/Q`$. Closer to the two-jet region sub-leading power corrections of the form $`\lambda _n(\overline{\mathrm{\Lambda }}/tQ)^n`$ become important. Strictly based on the large $`\beta _0`$ renormalon calculation, only odd values associated with renormalon poles at half integers, are relevant. However, since the large $`\beta _0`$ limit is not expected to predict the singularity strength of the exact Borel transform, one should expect $`\lambda _n(\overline{\mathrm{\Lambda }}/tQ)^n`$ with even $`n`$ to survive. Assuming that the exact Borel function differs from the calculated one only by a modification of the singularity strength, with the $`\nu `$ dependence left unchanged, we predict that the even terms will be suppressed by $`1/\mathrm{ln}(Q^2/\overline{\mathrm{\Lambda }}^2)`$ compared to the odd ones (see 3.5). Allowing both odd and even terms in the sum in (93) one recovers the shape-function description as far as the large-angle gluon emission is concerned. Nevertheless, at the same time a collinear correction proportional to $`\lambda _2(\overline{\mathrm{\Lambda }}/tQ)^2`$ appears in (3.5). The latter has a sign apposite to the relevant term in the large-angle shape-function, with a slightly different functional form, thus reducing significantly the overall effect of the non-perturbative parameter $`\lambda _2`$.

Finally, we comment that an alternative way to derive the form of the power corrections from the large $`\beta _0`$ calculation is to use the characteristic function (3.2). Expanding the latter at small $`ϵ`$ one obtains

$`\dot{}_\nu \left(ϵ\right)`$ $`=`$ $`2\sqrt{ϵ}\nu +\left({\displaystyle \frac{1}{2}}\nu \mathrm{ln}\left({\displaystyle \frac{1}{ϵ}}\right)3\nu {\displaystyle \frac{1}{2}}\nu ^2\right)ϵ+\left(\nu +{\displaystyle \frac{1}{2}}\nu ^2+{\displaystyle \frac{1}{9}}\nu ^3\right)ϵ^{3/2}`$
$`+`$ $`\left({\displaystyle \frac{1}{4}}\nu ^2\mathrm{ln}\left({\displaystyle \frac{1}{ϵ}}\right){\displaystyle \frac{1}{12}}\nu ^3{\displaystyle \frac{1}{48}}\nu ^4\right)ϵ^2+\left({\displaystyle \frac{1}{6}}\nu ^3+{\displaystyle \frac{1}{72}}\nu ^4+{\displaystyle \frac{1}{300}}\nu ^5\right)ϵ^{5/2}`$
$`+`$ $`\left({\displaystyle \frac{5}{36}}\nu ^3+{\displaystyle \frac{1}{48}}\nu ^4{\displaystyle \frac{1}{480}}\nu ^5{\displaystyle \frac{1}{2160}}\nu ^6\right)ϵ^3`$
$`+`$ $`\left({\displaystyle \frac{1}{360}}\nu ^5+{\displaystyle \frac{1}{3600}}\nu ^6+{\displaystyle \frac{1}{17640}}\nu ^7\right)ϵ^{7/2}+\mathrm{}.`$
Here only non-analytic terms are associated with ambiguity of the perturbative sum (43), and therefore with power corrections . One can immediately identify the terms $`(\nu \sqrt{ϵ})^n`$ where $`n=1,3,5,\mathrm{}`$ as the half-integer renormalon poles associated with large-angle soft gluons. These are the most important power corrections. As for integer powers $`ϵ^p`$, the only non-analytic terms are at $`p=1`$ and $`p=2`$, due to the logarithm. These ambiguities, however, are negligible for any $`t\overline{\mathrm{\Lambda }}^2/Q^2`$ since each power of $`ϵ`$ is accompanied by $`\nu `$ rather than by $`\nu ^2`$. Thus the picture is consistent with what we already learned based on the Borel analysis in the case that the singularities at integer $`z`$ are just simple poles and are therefore cancelled by $`\mathrm{sin}(\pi z)`$ factor from the analytic continuation.

## 4 Relation to fixed logarithmic accuracy and fixed-order calculations

In order to prove consistency of DGE with the exact exponent to NLL accuracy we will now compare our calculation to that of . This will lead to identification of coupling $`\overline{A}`$ as the “gluon bremsstrahlung” coupling . Next we briefly recall the way to match the resummed result to the fixed-order calculation.

### 4.1 NLL resummation from an evolution equation

The calculation of the differential cross-section by DGE requires to identify the coupling $`\overline{A}`$. The analogy with the Abelian theory  suggests that there is a unique all-order definition of $`\overline{A}`$ (analogous to the Gell-Mann–Low effective charge) which does not depend on the observable considered, such that (5) is the leading term in the skeleton expansion. Since a systematic skeleton expansion has not yet been constructed in QCD other considerations must be used. The analogy with the Abelian theory allows one to fix the large $`\beta _0`$ coupling, as done in (6). To go beyond this limit we require consistency of the exponent to NLL accuracy. We will see that this amounts to identifying the coupling $`\overline{A}`$ as the “gluon bremsstrahlung” coupling .

The standard way of calculating the logarithmically enhanced cross-section  is based on an evolution equation for the jet mass distribution $`J(Q^2,k^2)`$, where the kernel is the splitting function of a gluon off a quark, $`P(a,x)`$. At leading order the splitting function is

$$P(a,x)=C_Fa\left[\frac{x}{1x}+\frac{1}{2}(1x)\right]_+=C_Fa\left[\frac{1}{1x}\frac{1}{2}(1+x)\right]_+$$
(99)
where $`a(k^2)`$ is the coupling ($`\alpha _s(k^2)/\pi `$) and $`x`$ is the longitudinal momentum fraction in the branching. Fixing the coupling to be the “gluon bremsstrahlung” coupling ,

$$a\stackrel{~}{a}a_{\overline{\text{MS}}}+\left[\frac{5}{3}\beta _0+\left(\frac{1}{3}\frac{\pi ^2}{12}\right)C_A\right]a_{\overline{\text{MS}}}^2+\mathrm{}$$
(100)
this splitting function is correct to next-to-leading order, as far as the singular $`1/(1x)`$ terms are concerned. At NLL accuracy, the evolution equation is ,

$$\frac{dJ_\nu (Q^2)}{d\mathrm{ln}Q^2}=_0^1𝑑xP(\stackrel{~}{a}((1x)Q^2),x)\left[e^{\nu (1x)}1\right]J_\nu (Q^2)$$
(101)
where $`J_\nu (Q^2)`$ is the Laplace transform of $`J(Q^2,k^2)`$ with respect to $`k^2/Q^2`$. The solution of this differential equation is

$$\mathrm{ln}J_\nu (Q^2)=_0^1𝑑u\left[e^{\nu u}1\right]\frac{1}{u}_0^{1u}𝑑xP(\stackrel{~}{a}((1x)uQ^2),x).$$
(102)
To this accuracy $`1\text{thrust}`$, $`t`$, is the sum of invariant masses of the two hemispheres, so

$$R(Q^2,t)=_0^{\mathrm{}}𝑑k^2𝑑\overline{k}^2J(Q^2,k^2)J(Q^2,\overline{k}^2)\theta \left(tQ^2k^2\overline{k}^2\right),$$
(103)
and therefore,

$$R(Q^2,t)|_{\text{res}}=_𝒞\frac{d\nu }{2\pi i\nu }e^{\nu t}\mathrm{exp}\left\{2\mathrm{ln}J_\nu (Q^2)\right\}.$$
(104)
From here, a straightforward calculation gives the leading and next-to-leading logs in the thrust distribution.

Let us now compare this calculation with ours. Specifically, $`2\mathrm{ln}J_\nu (Q^2)`$ in (104) should be compared with $`e^{\nu t}_{\text{SDG}}`$ in (36). The two quantities are defined in (102) and (35), respectively. The first observation is that the integration variable $`u`$ in (102) can be identified with $`t`$, making the two similar provided one identifies

$$\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}\frac{2}{t}_0^{1t}𝑑xP(\stackrel{~}{a}((1x)tQ^2),x),$$
(105)
or, explicitly (see 29),

$$_{t^2}^t\frac{dϵ}{ϵ}\overline{A}_{\text{eff}}(ϵQ^2)\left[\frac{2}{t}\frac{ϵ}{t^2}\frac{ϵ^2}{t^3}\right]\frac{2}{t}_0^{1t}𝑑x\beta _0\stackrel{~}{a}((1x)tQ^2)\left[\frac{1}{1x}\frac{1}{2}(1+x)\right].$$
(106)
If one identifies $`\overline{A}_{\text{eff}}(ϵQ^2)`$ with $`\beta _0\stackrel{~}{a}((1x)tQ^2)`$, and changes variables $`ϵ=(1x)t`$, the r.h.s. becomes

$$_{t^2}^t\frac{dϵ}{ϵ}\overline{A}_{\text{eff}}(ϵQ^2)\left[\frac{2}{t}\frac{2ϵ}{t^2}+\frac{ϵ^2}{t^3}\right].$$

One thus finds that the first term in the square brackets on both sides is the same. Recall that this term is responsible for the leading logs. The other terms, contributing to sub-leading logs, are not the same. However, the NLL are: to NLL accuracy it is enough to consider the conformal part  by fixing the scale of the coupling. The latter gives, in both cases, $`\frac{1}{t}\frac{3}{2}\overline{A}(Q^2)`$ from the collinear limit $`ϵ=t`$.

We thus learned that DGE is consistent with the exact exponent to NLL accuracy provided one identifies the running coupling with the “gluon bremsstrahlung” coupling  of eq. (100). In evaluating the Borel sum (38) we therefore use

$$\overline{A}(k^2)=\frac{A(k^2)}{1\left[\frac{5}{3}+\left(\frac{1}{3}\frac{\pi ^2}{12}\right)C_A/\beta _0\right]A(k^2)}.$$
(107)
instead of (6). This simply amounts to fixing the definition of $`\overline{\mathrm{\Lambda }}`$ (see (27)).

The qualitative features of the SDG characteristic-function based calculation appear also in the splitting-function based calculation: there appear $`ϵ/t^2`$ and $`ϵ^2/t^3`$ terms which contribute to the log-enhanced cross-section due to the collinear limit $`ϵt`$. As we saw in the previous sections these terms lead to factorial enhancement of sub-leading logs (already prior to exponentiation). However, the actual values of sub-leading logs predicted by the splitting-function beyond the NLL are different. Performing the integral on the r.h.s. of (106) explicitly we arrive at the following expansion for $`\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)`$ prior to exponentiation:

$`{\displaystyle \frac{1}{\sigma }}{\displaystyle \frac{d\sigma }{dt}}(Q^2,t)={\displaystyle \frac{C_F^{}}{t\beta _0}}\{(2.L1.5)A(Q^2)+(3L^2+1.8333L4.2500)A(Q^2)^2`$ (108)
$`+(4.6667L^3+8.5000L^29.5242L8.8151)A(Q^2)^3`$
$`+(7.5000L^4+21.833L^317.359L^250.085L9.957)A(Q^2)^4`$
$`+(12.400L^5+48.500L^431.339L^3194.72L^295.096L+10.97)A(Q^2)^5+\mathrm{}`$
This expansion can be compared directly with (25). The LL and NLL are of course the same. Sub-leading logs are similar in many cases; note in particular the NNLL. Nevertheless, at lower logs there are big variations, as one can deduced directly from (106).

The fact that the two calculations differ beyond NLL is of no surprise. Most importantly, the kinematic approximations made in the splitting-function based calculation (e.g. the matrix element, the definition of the observable), modify the logs beyond NLL accuracy.

It should be emphasised that some approximations have been made also in the characteristic-function based calculation. Although this calculation is exact at the level of three partons with an off-shell gluon, there are approximations beyond this level: the massive-gluon characteristic-function does not take into account the possible branching of the gluon into opposite hemispheres, the so-called non-inclusive contribution. When the gluon is roughly collinear with one of the quarks (the only limit contributing to sub-leading logs in $`\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}`$) this type of branching is probably not significant. Indeed, the comparison between $`\frac{1}{\sigma }\frac{d\sigma }{dt}(Q^2,t)|_{\text{SDG}}`$ and the leading term in $`\beta _0`$ in the (numerically evaluated) exact NLO coefficient in sec. 2.2 confirms this assertion. Another important aspect is the validity of the exponentiation formula (37) beyond NLL accuracy. As discussed in sec. 3, correlations in the emission of multiple soft and collinear gluons may influence the sub-leading logs. This effect was neglected.

### 4.2 Matching to the fixed-order calculation

At large values of $`t=1\text{thrust}`$, the resummed cross-section (be it DGE or NLL) must be complemented by the exact fixed order result. To do this we will use the $`\mathrm{log}`$-R matching scheme. The basic idea of matching is to replace the low orders in the resummed expression by the exact low order terms. Schematically one writes  $`R_{\mathrm{matched}}=R^{\text{NLO}}+R_{\text{res}}R_{\text{res}}^{\text{NLO}}`$ where $`R_{\text{res}}`$ is the resummed cross-section and $`R_{\text{res}}^{\text{NLO}}`$ is its expansion to NLO. There is no general argument for matching the integrated cross-section itself (R-matching) as compared to matching some function of $`R`$. Exponentiation may hold beyond the level of the logarithmically enhanced cross-section, so it is natural to match the logarithm of $`R`$ (log-R matching). In practice the difference between different matching procedures is small (see e.g. ) and for simplicity we only consider $`\mathrm{log}`$-R matching.

Using the known NLO calculation ,

$$R(Q^2,t)1+R_1(t)a+R_2(t)a^2+\mathrm{}$$

where $`R_1(t)`$ is known analytically and $`R_2(t)`$ is parametrised based on a numerical computation by EVENT2 , we write

$$R(Q^2,t)=\mathrm{exp}\left\{R_1a(Q^2)+\left(R_2\frac{1}{2}R_1^2\right)a(Q^2)^2+\mathrm{ln}R(Q^2,t)|_{\text{res}}\mathrm{ln}R(Q^2,t)|_{\text{res}}^{\text{NLO}}\right\}.$$
(109)
Here $`\mathrm{ln}R(Q^2,t)|_{\text{res}}`$ is our resummed result for the logarithmically enhanced integrated cross-section, given by eqs. (67) and (66) and $`\mathrm{ln}R(Q^2,t)|_{\text{res}}^{\text{NLO}}`$ is its expansion, up to NLO in $`\overline{\mathrm{MS}}`$,

$`\mathrm{ln}R(Q^2,t)|_{\text{res}}^{\text{NLO}}`$ $`=`$ $`[C_FL^2+{\displaystyle \frac{3}{2}}C_FL]a(Q^2)+[C_F\beta _0L^3^{}`$
$`+`$ $`\left({\displaystyle \frac{11}{12}}C_F\beta _0{\displaystyle \frac{1}{3}}\pi ^2C_F^2+\left({\displaystyle \frac{1}{12}}\pi ^2{\displaystyle \frac{1}{3}}\right)C_AC_F\right)L^2`$
$`+`$ $`\left({\displaystyle \frac{15}{4}}C_F\beta _0+\left(4\zeta _3+{\displaystyle \frac{1}{2}}\pi ^2\right)C_F^2+\left({\displaystyle \frac{1}{2}}{\displaystyle \frac{1}{8}}\pi ^2\right)C_AC_F\right)L`$
$`+`$ $`(3\zeta _3{\displaystyle \frac{1}{180}}\pi ^4{\displaystyle \frac{3}{16}}\pi ^2)C_F^2]a(Q^2)^2.`$
The matching (109) is done in the $`\overline{\mathrm{MS}}`$ scheme.

## 5 The phenomenological implication

We saw that a careful treatment of running coupling effects (or renormalons) in the logarithmically enhanced cross-section requires resummation of sub-leading logs. Our investigation in section 3.2 leaves no doubt: sub-leading logs contribute significantly at all physically interesting values of the thrust. On the other hand, it is well known that successful fits to the distribution were obtained in  with a much simpler treatment of the perturbative calculation, based on the standard NLL resummation. Non-perturbative corrections whose magnitude is not under control theoretically are anyway important, and must be included in any fit. One may argue that a more complete perturbative treatment is not necessary, and that it amounts simply to a redefinition of the non-perturbative correction. The answer is, of course, that as soon as quantitative predictions are expected one should start with a reliable perturbative calculation. Physically, the transition between the perturbative and non-perturbative regimes is smooth, and one must make sure that an arbitrary separation of the two does not have too strong an impact on the results. Legitimate separations are cut-off regularization of momentum integrals and principal-value Borel summation . It was shown in  that these two procedures can be fully consistent with each other provided that renormalon resummation is performed. This means in particular that the extracted value of $`\alpha _s`$ is the same using the two prescriptions. Such consistency cannot be achieved if a truncated perturbative expansion, e.g. NLO, is used instead. Our results in sec. 3.2 imply that NLL truncation of the exponent is just as dangerous. By performing renormalon resummation, such a truncation is avoided and a more reliable comparison with data can be achieved.

Here we shall perform global fits to the thrust distribution at all energies, based on our perturbative calculation of section 3. As discussed in section 3.3 we use Borel summation. An equivalent calculation can be done using cut-off regularization. The translation of our results to this language is straightforward and may be useful once comparison can be made to a non-perturbative calculation of the shape-function with a hard cut-off as an ultraviolet regulator.

The perturbative calculation must be supplemented by non-perturbative corrections in order to be compared with the data. Our guiding principle is that the functional form of the non-perturbative corrections is dominated by the ambiguities identified in the calculation of the perturbative sum. We assume that the SDG based calculation is sufficient to identify these ambiguities. Of course, there are inherent limitations to this calculation: splitting of a gluon into partons that end up in opposite hemispheres is not taken into account, nor are correlations between the emitted gluons. These effects are not expected to be important in the two-jet region, so it is reasonable to neglect them.

### 5.1 Perturbative calculation

Before attempting any fits to the data involving non-perturbative corrections, it is useful to examine the impact of the resummation simply by looking at the perturbatively calculated cross-section in fig. 4. Needless to say, the fixed-order calculations, LO and NLO, do not describe the physical cross section in the two-jet region. The large difference between the NLO and the LO signals large higher-order corrections. So does the renormalization scale dependence of the result (not shown in the plot). The qualitatively different behaviour of the NLL result is clearly seen: thanks to exponentiation it vanishes at $`t=0`$. Nevertheless, as fig. 4 shows, the renormalization scale dependence is still appreciable<sup>*</sup><sup>*</sup>*It has been shown  that renormalon resummation in single scale observables can be imitated by an appropriate choice of the renormalization scale, similarly to the original formulation of the BLM approach . One might wonder to what extent DGE can be imitated by simply tuning the scale of the coupling in the standard NLL formula. Physically, the naive scale-setting approach is bound to fail for multi-scale observables such as differential cross-sections, since the typical scale of emitted gluons is not anymore a fixed fraction of the center-of-mass energy, but rather some function of all the external parameters. In our case, it is a function of both $`Q`$ and $`t`$. Indeed, choosing lower scales for the coupling in the NLL formula one cannot reproduce the DGE result.. We stress that the factor $`2`$ between the renormalization points is completely arbitrary. This difference should not be considered as a reliable estimate of the error due to running coupling effects that have been neglected, since the physical scale, (e.g. the transverse momentum) in the two-jet region is significantly smaller than $`Q`$. As anticipated in sec. 3.2, the impact of the additional resummation of sub-leading logs by DGE is quite significant. Having performed this resummation, the renormalization scale dependence (at the level of the logarithmically enhanced cross-section) is avoided.

In fig. 4, as in all forthcoming calculations we use the so-called log-R matching scheme  , defined in Eq. (109), to complement the calculation of the logarithmically enhanced cross-section. Renormalization scale dependence appears only beyond the level of the logarithmically enhanced cross-section and it is therefore suppressed by $`t`$. Numerically, the residual scale dependence (see table 1) is very small as long as the log-enhanced cross-section dominates. The difference is much too small to be observed in the resolution of fig. 4. At large $`t`$ ($`t\begin{array}{c}>\hfill \\ \hfill \end{array}0.2`$) the terms which are suppressed by $`t`$ start making an impact and eventually the scale dependence becomes comparable to its values in the NLL and the NLO results.

Next, we confront the question of reliability of the approximation technique we used in section 3.3 in evaluating the Borel sum. In the calculation presented in fig. 4 we evaluated the Borel sum in (66) using $`p_{\mathrm{max}}=8`$. This means that the function $`V(z)`$ in eq. (55) was approximated by $`8`$ poles. Note that $`p_{\mathrm{max}}`$ is also the largest power of $`\mathrm{ln}\nu `$ respected by this approximation (see (57)). In addition, the renormalon integrals $`\stackrel{~}{I}_p(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2)`$ were evaluated using the two-loop running coupling in the “gluon bremsstrahlung” scheme (see section 4). How sensitive are the results to these approximations? To answer this question we repeat the calculation varying the choice of $`p_{\mathrm{max}}`$ and the running coupling. Considering first the choice of $`p_{\mathrm{max}}`$, we show in fig. 5 approximants of increasing order.

The convergence is quite remarkable. In the upper plot the $`p_{\mathrm{max}}=5`$ curve (dashed) can hardly be distinguished from the $`p_{\mathrm{max}}=8`$ curve. Indeed, as the lower plot shows, the difference between them is very small and it exceeds $`1\%`$ only around the distribution peak. For our fits, we choose $`p_{\mathrm{max}}=8`$ which guarantees a truncation errorThis truncation error must not be confused with the ambiguity of the perturbative sum. of less than $`1\%`$ down to $`t=0.01`$ at $`\mathrm{M}_\mathrm{Z}`$ (or $`0.01\mathrm{M}_\mathrm{Z}/Q`$ at $`Q`$).

The stability of the calculation with respect to the running coupling can be checked by repeating the calculation with a different renormalization-group equation. In principle, one may also consider changing the definition of $`\mathrm{\Lambda }`$, thus replacing the “gluon bremsstrahlung” coupling by some other coupling that matches the same Abelian limit. This was done in a similar context in the case of the average thrust . It was found there that such a replacement introduces a non-negligible variation in the calculated observable, but the final impact on the extracted value of $`\alpha _s`$ was rather small (of order $`1\%`$). In the context of the distribution, the choice of the “gluon bremsstrahlung” coupling is almost imposed on us by consistency with the NLL exponentiation kernel. Therefore we shall not make such a modification. The simplest variation in the renormalization-group equation is to replace it by the one-loop. In this case the renormalon integrals $`\stackrel{~}{I}_p(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2)`$ are evaluated by eq. (63) instead of eq. (65). Such a comparison is presented in fig. 6.

The difference between the two curves is quite significant, in particular in the peak region. Note, however, that contrary to the two-loop case, the one-loop resummed result does not include all the next-to-leading logs (the terms proportional to $`\beta _1`$ are omitted). The difference between two- and three-loopIn the spirit of  and our discussion in sec. 4, the appropriate value for $`\beta _2`$ should probably be set based on the NNLO coefficient of the singular term $`1/(1x)`$ in the splitting function. running coupling is expected to be much smaller. However, it may well be non-negligible close to the peak, so it is certainly worth consideration in the future. In the following we will use the two-loop running coupling.

### 5.2 Non-perturbative corrections

It was realized before  that non-perturbative power corrections of the form $`\lambda _n(\mathrm{\Lambda }/tQ)^n`$ to any power $`n`$ become relevant around the distribution peak. As shown in section 3.5, such corrections are indeed expected based on the form of the ambiguity of the perturbative result. If all the powers are as important, one would need an infinite number of parameters – or a completely arbitrary shape-function – to bridge the gap between the perturbative result and the data. In practice one assumes that there is a region in $`t`$ where a power corrections expansion such as (93) is dominated by its first few terms, or in other words, that the shape-function can be described by its first few moments. It is therefore essential that $`\lambda _n`$ do not grow too fast with $`n`$. Having made this assumption, one can put a cut $`t_L`$ on the fitted range in $`t`$ for any given energy, and fit the data for $`t>t_L`$ using a reasonably flexible parametrisation of the shape-function. This cut should scale like the distribution peak, i.e. as $`1/Q`$,

$$t_L(Q)=t_L(M_Z)\frac{M_Z}{Q}.$$
(111)
In any case, the dependence of the results on the precise location of this cut and on the specific functional form of the shape-function must not be too large.

The simplest possibility is to place the cut at large enough values of $`t`$ such that the only relevant non-perturbative parameter will be $`\lambda _1`$, which controls the shift of the distribution (90). This approach was suggested in . The advantage is that the perturbative calculation is better under control and the “model dependence” in the parametrisation of non-perturbative corrections is minimized. The price, however, is that one does not use a large part of the data, and in particular, one excludes the region where the distribution changes sharply and is therefore expected to be most constraining and most informative.

Here we use various different minimal $`t`$ cuts and perform fits where a flexible shape-function is convoluted with the perturbative distribution as well as fits with just a shift. For the shape-function-based fits the minimal cut can be quite low, thus using almost all available data. At the low end we choose $`t_L(\mathrm{M}_\mathrm{Z})=0.01`$. This is equivalent to $`Qt\begin{array}{c}>\hfill \\ \hfill \end{array}4\overline{\mathrm{\Lambda }}`$ and it guarantees convergence of the Borel function approximation within $`1\%`$, as shown in fig. 5. Less ambitious cuts $`t_L(\mathrm{M}_\mathrm{Z})=0.02`$ and $`0.03`$ are useful to check the stability, which is indicative of whether the power-correction expansion is convergent enough. For the shift-based fit higher cuts must be used. Also here stability of the results is important to guarantee that sub-leading power corrections are not too relevant in the fitted range.

For simplicity we shall assume that the parameters of the shift or the shape-function are $`Q`$ independent, even though we have shown that logarithmic dependence is expected in general. In principle, a fit that allows logarithmic dependence of $`\lambda _n`$ according to eq. (3.5) is more appropriate. But, since $`\kappa _p`$ are not known, the price will be having more parameters which can only be constrained by very precise data at several different energies. The available data is not constraining enough for such a fit.

#### 5.2.1 Fits with a shifted distribution

Let us consider first a fit based on shifting the perturbative distribution. We use experimental data from all energies in the range $`14\mathrm{GeV}Q189\mathrm{GeV}`$ (summarised in table 5 below). When available, the systematic errors have been added in quadrature to the statistical errors. Table 2 summarizes the fit results for $`\alpha _s`$ and $`\lambda _1`$ as a function of the low cut.

First of all we note that the fits are very good: the $`\chi ^2/\mathrm{dof}`$ values are unusually low. Next, we see that the results depend on the low cut. Particularly stable results, as far as the extracted value of $`\alpha _s`$ is concerned, are obtained for $`t_L(\mathrm{M}_\mathrm{Z})\begin{array}{c}>\hfill \\ \hfill \end{array}0.06`$. Above this value consecutive cuts yield $`\alpha _s`$ values that are within errors of each other. This stability means that sub-leading power corrections in this region are small enough to be neglected. Thus we regard the result $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})0.110`$ as reliable. Too high cuts $`t_L(\mathrm{M}_\mathrm{Z})0.10`$ are disfavoured because the number of fitted points decreases significantly (note that the propagated experimental error increases accordingly). Cuts like $`t_L(\mathrm{M}_\mathrm{Z})=0.05`$ have the advantage of increased sensitivity to $`\lambda _1`$. Lower cuts are probably too sensitive to sub-leading power corrections. This is reflected in the stronger dependence on the cut and in the deterioration of $`\chi ^2/\mathrm{dof}`$. In the following we shall use $`t_L(\mathrm{M}_\mathrm{Z})=0.05`$ for illustrations. It has the additional advantage of consistency with the shape-function-based fit (see below).

The table and fig. 7 show that there is a strong correlation between $`\alpha _s`$ and $`\lambda _1`$. A similar correlation was found in . This means that quantitative discussion of the values of non-perturbative parameters, like $`\lambda _n`$ is relevant only for a fixed $`\alpha _s`$.

In all the fits discussed so far we used an upper cut $`t_H=0.32`$ which is $`Q`$ independent. We checked that the location of this cut does not influence the fit. For example, keeping the lower cut fixed $`t_L(\mathrm{M}_\mathrm{Z})=0.05`$ and changing the upper cut between $`t_H=0.29`$ and $`t_H=0.35`$ (corresponding to a total of $`154`$ to $`184`$ fitted points, respectively) the central value of $`\alpha _s`$ changes by just $`0.2\%`$ (and that of $`\lambda _1`$ by $`4\%`$).

Figures 8 through 10 show the fitted distribution and the original perturbative distribution together with the experimental data. For large $`Q`$ we show only the small $`t`$ region, which is more interesting. The agreement between the shifted distribution and the data is good. Nevertheless, it is clear that around the peak additional corrections are required (in the case of $`\mathrm{M}_\mathrm{Z}`$, for example). Indeed we saw that such corrections are expected based on the renormalon ambiguity.

#### 5.2.2 Fits with a shape-function

Consider now a fit where the perturbative distribution is convoluted with a shape-function. According to sec. 3.5, this means, in fact, two convolutions: one with the large-angle shape-function $`f_{\mathrm{la}}(ϵ)`$, defined in eqs. (92) and (93), and another which corrects for the collinear contribution where we only include the leading term in (3.5),

$$f_{\mathrm{col}}(Q^2,t)=\delta (t)\lambda _2\left[\frac{\overline{\mathrm{\Lambda }}^2}{Q^2t}\right]_+\delta ^{(1)}(t).$$

The main difference with previous works where a shape-function was used  is in the starting point: the perturbative distribution there is the standard NLL result (with some cut-off regularization) while we start with the DGE result avoiding the truncation of sub-leading logs. The latter provides an analytic continuation of the perturbative treatment beyond its original range of applicability.

In spite of the significant difference in the perturbative distribution as compared to , we consider their ansatz for the shape-function appropriate also in our case, barring the fact that it is positive definite. To allow in principle the suppression of even moments as implied by our analysis in sec. 3.5, we consider a function that is flexible enough:

$$f_{\mathrm{la}}(ϵ)=n_0ϵ^q(1+k_1ϵ+k_2ϵ^2)e^{b_1ϵb_2ϵ^2}$$
(112)
where $`n_0`$ is fixed by the normalization requirement $`f(ϵ)𝑑ϵ=1`$ and $`q,k_1,k_2,b_1`$ and $`b_2`$ are free parameters to be fixed by the fit to the data. We assume $`q>0`$ so that $`f_{\mathrm{la}}(0)=0`$ and $`b_2>0`$ (or $`b_2=0`$ and $`b_1>0`$) so that the shape-function decreases at large $`ϵ`$.

Note that with a shape-function of this form, the central moments $`\lambda _n`$ do not have a definite sign. In particular both positive and negative values for $`\lambda _1`$ are possible. When a cut-off regularization is used it is intuitively clear that the shift is to the right: soft gluons widen the jets. Indeed $`\lambda _1`$ was found to be positive in both the shift-based fit in  and in the shape-function  based fit. The shift parameter $`\lambda _1`$ was shown to be roughly consistent with the leading power correction to average event-shapes , which is always positive. The intuitive argument is not anymore valid when principal value regularization is used instead. According to , in this regularization $`\lambda _1`$ is still positive. We stress, however, that there is no general reason to exclude a priori the possibility of a negative $`\lambda _1`$.

Performing a fit with $`f_{\mathrm{la}}(ϵ)`$ and the collinear correction (112) with a lower cut $`t_L(\mathrm{M}_\mathrm{Z})=0.01`$ and an upper cut $`t_H=0.32`$, we obtained a very good fit with $`\chi ^2/\mathrm{dof}=0.84`$ (for a total of 246 fitted data points). The fit results are shown in figures 8 through 10. The entire range of $`t`$ is shown for all energies in fig. 11.

The parameters in this case are

(117)
corresponding to

$`\lambda _1`$ $`=`$ $`2.34\pm 0.21`$ (118)
$`\lambda _2`$ $`=`$ $`0.37\pm 1.38.`$
It turns out that the fit has a strong preference for a linear term in the exponent, so that if both $`b_1`$ and $`b_2`$ are free, then $`b_2`$ is set by the fit to zero. In order to try another functional form, thus testing the sensitivity of the fit results to the properties of the shape-function, we also fitted the data in the same range with the function where the linear term in the exponent is missing: $`b_10`$. In this case we got a best fit with $`\chi ^2/\mathrm{dof}=0.88`$, i.e. slightly worse than the previous one, with the following parameters,

(123)
corresponding to

$`\lambda _1`$ $`=`$ $`2.92\pm 0.14`$ (124)
$`\lambda _2`$ $`=`$ $`3.52\pm 0.57.`$
We see that most of the parameters of the shape-function are fairly well constrained by the fit. Unfortunately, this does not imply that the central moments of the shape-function are well determined. As explained above, due to the correlation between $`\alpha _s`$ and the non-perturbative parameters, it makes sense to fix $`\alpha _s`$ in the fit when comparing the results for the two possible exponents. Doing so with $`\alpha _s=0.1093`$ we still obtain good fits in both cases. For a linear term in the exponent we get $`\lambda _1=2.35\pm 0.05`$ ($`\chi ^2/\mathrm{dof}=0.84`$) and for a square in the exponent we get $`\lambda _1=2.45\pm 0.04`$ ($`\chi ^2/\mathrm{dof}=0.92`$). Now the values of $`\lambda _1`$ are much closer. Note that also the propagated errors become smaller, since $`\alpha _s`$ and $`\lambda _1`$ are so strongly correlated.

Figure 7 shows how $`\lambda _1`$ varies when the fits with the shift and with the shape-function (with a linear term in the exponent) are performed with a fixed value of $`\alpha _s`$. As intuitively expected, forcing a larger value of $`\alpha _s`$, $`\lambda _1`$ decreases. For $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})0.109`$ similar values are obtained for $`\lambda _1`$ from the shift and from the shape-function (note that this depends on the chosen lower cut), however, for higher values of the coupling the values of $`\lambda _1`$ start differing. In this case the shape-function is not so well approximated by a $`\delta `$-function (a shift) and its higher central moments $`\lambda _n`$ play a role.

To allow a quantitative discussion concerning the higher central moments let us fix also $`\lambda _1`$ and compare the fits with the two exponents for the same $`\alpha _s`$ and $`\lambda _1`$. We choose them to be the values obtained in our best bit $`(\text{117})`$: $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})=0.1093`$ and $`\lambda _1=2.34`$.

Quantitative information on higher moments is theoretically important for several reasons. First, it would be interesting to see the effect of the collinear contribution. At the same time it would be interesting to check our prediction that $`\lambda _2`$ (like other even central moments of the shape-function) is suppressed. Even ignoring the strong dependence on $`\alpha _s`$, these two ingredients are hard to address by fitting the data: if $`\lambda _2`$ is small then the collinear correction is not important, and conversely – if the collinear correction is included a major part of the effect of $`\lambda _2`$ is cancelled, so $`\lambda _2`$ would be very hard to access. To deal with this we make three different fits with the same fixed values of $`\alpha _s`$ and $`\lambda _1`$ mentioned above: with the collinear correction, without it, and finally a fit where $`\lambda _2`$ is fixed to zero. The results are summarised in table 3.

The first conclusion is that $`\lambda _2`$ is not well determined by the fits: even for fixed values of $`\alpha _s`$ and $`\lambda _1`$ the two functional forms yield somewhat different results. Nevertheless, a vanishing $`\lambda _2`$ is certainly not inconsistent with the data: the fit does not deteriorate much by fixing it to zero. One also sees that the collinear correction is not important: including it or not leaves $`\lambda _2`$ almost identically the same. Had the value of $`\lambda _2`$ been significantly different from zero, such a situation would have been impossible, but for a small (and not very well constrained) $`\lambda _2`$ it is quite reasonable.

Table 3 also gives the $`\lambda _3`$ values. It is clear that $`\lambda _3`$ is not constrained enough unless $`\alpha _s`$, $`\lambda _1`$ and $`\lambda _2`$ are all fixed. Still, one thing that can be said in general (varying the other parameters, including $`\alpha _s`$) is that negative values for $`\lambda _3`$ are strongly preferred. Thus, a positive-definite ansatz for the shape-function (in the principal-value regularisation!) is not adequate. Note that the preferred signs of $`\lambda _1`$ (positive) and $`\lambda _3`$ (negative) are compatible with $`\rho _{2p}`$ being all positive (see eq. (83) and (3.5)).

The shape-function-based fit can be considered an alternative way to extract the value of $`\alpha _s`$. We saw already, that in spite of the fact that the determination of the non-perturbative parameters is far from being complete, $`\alpha _s`$ does not change much. In particular, for the two functions used, $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})`$ varies only from $`0.1082`$ to $`0.1093`$. A further stability check is to what extent it is sensitive to the location of the low cut $`t_L(Q)`$. The answer is provided by table 4.

The sensitivity to $`t_L`$ in the determination of $`\alpha _s`$ is rather small. In particular it is smaller than in the case of a shift-based fit. There is still some tendency that $`\alpha _s`$ becomes larger for higher (thus less ambitious) cuts. While higher cuts are safer from the point of view suppressing sub-leading power corrections, they are less constraining as far as the multi-parameter shape-function is concerned. Note also that the propagated experimental error on $`\alpha _s`$ and $`\lambda _1`$ increases for higher cuts. “Ambitious cuts” where the peak region itself is included in the fit , e.g. $`t_L(\mathrm{M}_\mathrm{Z})=0.010`$, are therefore favoured. We conclude that according to the shape-function-based fit $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})`$ ranges between $`0.108`$ and $`0.111`$. These variations should be considered part of the theoretical uncertainty.

Finally, we summarise in table 5 the contribution to $`\chi ^2`$ in the fit by each experiment separately. The middle columns correspond to a fit of the shifted distribution for $`\alpha _s`$ and $`\lambda _1`$ with a lower cut $`t_L(\mathrm{M}_\mathrm{Z})=0.05`$ and an upper cut $`t_H=0.32`$. The columns on the right correspond to a fit for $`\alpha _s`$ and the parameters of the shape-function with a lower cut $`t_L(\mathrm{M}_\mathrm{Z})=0.01`$ and upper limit $`t_H=0.32`$. The OPAL and SLD data sets at $`91.2`$ GeV are not included in the fit since they are not consistent with the other data sets at this energy for small $`t`$. We note that there are no particular trends in the $`\chi ^2/\mathrm{point}`$ as a function of energy.

### 5.3 Moments of the thrust distribution

Given the assumption that the first few moments of the thrust $`t^n`$ are dominated by the two-jet configuration, and therefore by the logarithmically enhanced cross-section, one can relate  the parameters $`\lambda _n`$ to the power corrections of the moments:

$`t_{\mathrm{two}\mathrm{jet}}`$ $``$ $`t_{\text{PT}}+\lambda _1\overline{\mathrm{\Lambda }}/Q`$ (125)
$`t^2_{\mathrm{two}\mathrm{jet}}`$ $``$ $`t^2_{\text{PT}}+2\lambda _1t_{\text{PT}}\overline{\mathrm{\Lambda }}/Q+\lambda _2\overline{\mathrm{\Lambda }}^2/Q^2.`$
These predictions have been compared in  with the power correction from a single gluon emission, based on the characteristic function (2.1),

$`t_{\text{SDG}}`$ $``$ $`t_{\text{PT}}+\lambda \overline{\mathrm{\Lambda }}/Q+𝒪(\overline{\mathrm{\Lambda }}^3/Q^3)`$ (126)
$`t^2_{\text{SDG}}`$ $``$ $`t^2_{\text{PT}}+𝒪(\overline{\mathrm{\Lambda }}^3/Q^3).`$
The assumption that the two-jet region is dominant probably holds for the average thrust. The immediate conclusion is that the shift of the distribution $`\lambda _1`$ must coincide with the leading power-correction for the average, $`\lambda `$. The latter was determined from a fit in  using the principal value Borel sum regularization. For $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})=0.110`$, which is the best fit value in the case of the average thrust, $`\lambda \overline{\mathrm{\Lambda }}=0.62\pm 0.12`$ GeV. This number can be compared with our current fits to the distribution with a fixed $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})=0.110`$. We get $`\lambda _1\overline{\mathrm{\Lambda }}=0.44\pm 0.02`$ GeV in a shift-based fit with a cut $`t_L(\mathrm{M}_\mathrm{Z})=0.05`$ (corresponding to $`\chi ^2/\mathrm{dof}=0.67`$ for 177 points) and $`\lambda _1\overline{\mathrm{\Lambda }}=0.48\pm 0.01`$ GeV with a shape-function-based fit with a cut $`t_L(\mathrm{M}_\mathrm{Z})=0.01`$ (corresponding to $`\chi ^2/\mathrm{dof}=0.85`$ for 246 points). The power correction extracted from the distribution is somewhat lower than the corresponding power correction extracted from $`t`$, however the discrepancy is not large. Note that there is a difference between the resummation procedures used in the two cases. In the case of the distribution DGE is used, where the exponentiation plays a major role but terms that are suppressed by $`t`$ are neglected. On the other hand in case of the average a single dressed gluon renormalon sum is used, where terms that are suppressed by $`t`$ (and higher powers) are included but multiple gluon emission is neglected (apart from the NLO term).

An even more difficult case is that of the second moment $`t^2`$. As explained in , the absence of the term $`2\lambda _1t_{\text{PT}}\overline{\mathrm{\Lambda }}/Q`$ in (126) is an artifact of using the single gluon approximation. On the other hand the absence of the term $`\lambda _2\overline{\mathrm{\Lambda }}^2/Q^2`$ is physically meaningful (a SDG analysis is the appropriate tool to identify such a power-correction, if it exists). Our current finding that the effect of $`\lambda _2`$ is cancelled at the level of the distribution is consistent with that of eq. (126). It is interesting to note, however, that eq. (126) was derived in  based on the entire characteristic function, while our current conclusion concerning $`\lambda _2`$ is based on the logarithmically enhanced cross-section alone. If two-jet dominance applies to $`t^2`$, fixing $`\lambda _1`$ and $`\lambda _2`$ by the distribution would allow the determination of the power-corrections to $`t^2`$ based on (125). Unfortunately, this is too optimistic: as one can verify, e.g. by looking at the data for the distribution scaled by $`t^2`$ (see ), at $`\mathrm{M}_\mathrm{Z}`$ this observable gets a significant contribution from the region $`t\begin{array}{c}>\hfill \\ \hfill \end{array}0.2`$. As we saw at various stages of the current investigation (see e.g. fig. 1 and table 1) the logarithmically enhanced cross-section ceases to dominate at $`t\begin{array}{c}>\hfill \\ \hfill \end{array}0.2`$. It follows that the second moment of the thrust depends on contributions we neglected here. Perturbative and non-perturbative physics of three jets<sup>§</sup><sup>§</sup>§Recently, non-perturbative effects in three-jet observables were addressed ., including terms that are sub-leading in $`t`$ at NNLO, must be taken into account for a reliable analysis of this quantity.

### 5.4 Theoretical uncertainty

As always, there are unproven assumptions and various limitations to our approach. Here we address several issues shortly. A quantitative analysis is restricted to those approximations we have made where a particularly significant impact on the extracted value of $`\alpha _s`$ is expected.

First in line is the uncertainty in the perturbative calculation. Here there are two major issues: the first is the uncertainty due to missing NNLO calculations. The significance of NNLO corrections was emphasized recently in . In our analysis it is reflected for instance in the large renormalization scale dependence in the region $`t1/3`$. Once such a calculation is available, matching with the DGE would allow a significant reduction in the uncertainty, in particular, concerning the extraction of $`\alpha _s`$. The second issue concerns exponentiation: the DGE performed here is only approximate beyond NLL. An exact NNLL calculation would allow one to test certain aspects of this approach and complement the resummation.

Concentrating on the two-jet region we neglected effects that are explicitly suppressed by $`t`$, but renormalon factorial growth is not restricted to the logarithmically enhanced terms. Large numerical coefficients appear also in the “remainder function”. Fig. 1 shows that at $`t0.2`$, such sub-leading terms become important. Eventually, DGE may be generalized to included them. Currently it is not known whether these terms exponentiate or not (and therefore they are treated differently in different “matching schemes”). One can give an estimate of the possible impact of such terms by making a modification of the log, as suggested in , e.g. $`\mathrm{ln}1/t\mathrm{ln}(1/t1)`$. Such a modification ensures that the logarithmically enhanced cross-section will vanish at $`t=1/2`$, like the physical cross-section, rather than at $`t=1`$. Some results for a shift-based fit are summarized in table 6.

The general quality of these fits is worse than in the case of the non-modified log, table 2. In addition there is a stronger dependence of $`\alpha _s`$ on both cuts. These two facts suggest that such a modification is not quite the appropriate way to take into account the fact that the physical cross-section vanishes at $`t=1/2`$. Still, as a rough measure of the impact of the terms we neglected these results are useful. For the favoured cuts the difference in $`\alpha _s`$ is quite appreciable. It changes from $`0.110`$ to $`0.114`$. A similar study for the shape-function-based fit (with $`t_L(\mathrm{M}_\mathrm{Z})`$) shows a somewhat smaller change in $`\alpha _s`$: the central value changes from $`0.109`$ to $`0.112`$. Altogether we see that NNLO (and higher order) terms which are not logarithmically enhanced have a relatively significant impact on $`\alpha _s`$. We consider this as the largest source of uncertainty.

One aspect in which the DGE may be improved is the running coupling formula. Our entire analysis was performed with a two-loop running coupling. The comparison with the one-loop case in fig. 6 shows a relatively large variation. We expect the difference between three-loop and two-loop to be much smaller, but it may be non-negligible.

Next, there are uncertainties in the way hadronization corrections are treated. Here we assumed that the dominant non-perturbative corrections appear in a way that matches the ambiguity in the perturbative calculation. The success this general approach has had in event-shapes and elsewhere is encouraging. However, here we take it one step further: we assume that the $`\nu `$ dependence of each renormalon ambiguity (not only the first!) is reflected in the dependence of the corresponding non-perturbative corrections. Of course, any assumption we took in the perturbative calculation may have an impact on the prediction. In particular,

* Having neglected effects that are explicitly suppressed by $`t`$ in the perturbative calculation, we cannot expect the power correction analysis to be reliable beyond $`𝒪(1/t)`$ accuracy.
* Since exponentiation is based on two-jet kinematics, non-perturbative effects that are associated with three jets (where the recoil of the quark must be taken into account) are not included.
* The normalization of the thrust variable by the sum of energies ($`Q`$) in (14) (rather than by the sum of the absolute values of the three-momenta) was shown to affect the coefficient of the leading power-correction in the case of the average thrust . There is a similar impact on the $`1/Qt`$ correction in the case of the distribution, and probably also on higher power corrections.
* Having neglected the non-inclusive decay of the gluon in the perturbative calculation in sec. 2, the effect is missed also on the non-perturbative level. See discussion below eq. (25).
* In our framework it is very natural to treat non-perturbative corrections as if they all exponentiate, since the perturbative calculation indicating their necessity is a calculation of the exponent. This is most clearly put in eq. (70): the ambiguity appears in the exponent and so should be the power corrections. It is probably true that in the peak region, the leading power corrections indeed exponentiate, but in general there are power corrections that do not.

There are also hadronization effects that cannot be associated with ambiguities in the perturbative calculation. For example, the latter does not contain information on the spectrum of the hadrons produced. In a recent work Salam and Wicke  showed that hadron mass effects introduce important corrections, for average values of event-shape variables.

Another fact we ignored here is that in reality also heavy quarks are produced. Since the data contains all events whereas in the calculation the quarks are treated as massless, some systematic error is expected.

## 6 Conclusions

Many infrared safe observables in QCD are sensitive to soft and collinear gluon radiation. This sensitivity appears at the perturbative level in the form of large Sudakov logs. In order to deal with this situation one must perform resummation of the logarithmically enhanced cross-section. The resummation has two aspects: exponentiation and integrals over the running coupling. The Dressed Gluon Exponentiation (DGE) differs from previous resummation methods in the way the second aspect is addressed. The approach makes a direct link between the resummation of Sudakov logs and that of infrared renormalons. We showed that a careful treatment of running coupling integrals makes it possible to take into account a large class of logs, which is fully consistent with both the exact exponent up to NLL accuracy and renormalization group invariance. The latter is realized only when all the powers of the log are resummed. Moreover, sub-leading logs are factorially enhanced compared to the leading ones. As a consequence the standard approximation based on keeping the leading and the next-to-leading logs alone is not always justified.

DGE is a method to calculate a well defined class of log-enhanced terms based on exponentiation of the Single Dressed Gluon (SDG) renormalon-sum. The basic idea is that the exponent can be represented in terms of an integral over an observable-dependent characteristic function times the running coupling, similarly to the standard dispersive approach in renormalon calculations . This representation is natural since the exponentiation kernel is primarily associated with a single gluon emission. The method is quite general and can readily be applied to other physical quantities. The first stage is to calculate analytically the large $`N_f`$ or massive gluon characteristic function, and then to identify the phase-space limits in the integration over the gluon virtuality as well as the specific terms in the characteristic function that contribute to logarithmically enhanced terms in the perturbative coefficients. This allows to obtain the logarithmically enhanced part of the SDG cross-section to all orders. At the next stage the physical cross section is obtained by exponentiating the SDG result. In the case considered here this exponentiation is rather straightforward, given the assumption that gluons are emitted independently from the quarks and contribute additively to $`1`$ thrust. In making this assumption we ignored correlations in multi-gluon emission that occur at higher orders. As a result, the exponentiation is only approximate as far as sub-leading logs (NNLL and beyond) are concerned. DGE cannot replace the construction of the exact exponent as a solution of the evolution equation up to a given logarithmic accuracy. It is rather complementary to it, similarly to the way resummation in general is complementary to fixed-order calculations.

One of the main results of this work is the observation that sub-leading Sudakov logs are factorially enhanced compared to the leading logs. This observation was made first at the level of the SDG cross-section, before exponentiation. At this level factorial growth appears only due to the collinear limit of phase space. The Borel function (in the large $`\beta _0`$ limit) has simple poles at $`z=1`$ and $`2`$. However, the Borel integral is still well defined thanks to the analytic continuation. A much more dramatic factorial growth appears at the level of the exponentiated result in the Laplace conjugate variable $`\nu `$. Here both the collinear and large-angle (soft) limits of phase space contribute to infrared renormalons. As a result of the large-angle contributions, the first infrared renormalon appears at $`z=1/2`$, thus enhancing the divergence of the perturbative expansion as well as the relative significance of sub-leading logs compared to the leading logs. In the standard exponentiation formula  the most prominent enhancement of sub-leading logs, associated with the $`z=1/2`$ renormalon, appear due to the leading term in the splitting function, the same term that generates the leading logs.

The traditional formulation of resummation is based on writing the exponent as an expansion in the coupling, with coefficients $`g_n(\xi )`$ that are functions of the combination $`\xi L\beta _0\alpha _s/\pi `$, as in (1) (or in (46)). This way of organising the expansion is convenient when working with a fixed logarithmic accuracy. In this formulation, the enhancement of sub-leading logs translates into the properties of $`g_n(\xi )`$: these functions become more singular at $`\xi =1/2`$ and $`\xi =1`$ as $`n`$ increases, and in addition they increase factorially. As was shown in fig. 3, the convergence of the expansion (46) depends strongly on the coupling (and thus on $`Q^2`$). In the case of the thrust sub-leading logs are quite important at all relevant energies. When using DGE and keeping all the logs the standard way of organizing the expansion loses its attractiveness: as a consequence of the factorial divergence of the functions $`g_n(\xi )`$ at large $`n`$ the series must be truncated. Truncation at the minimal term would minimise the error, but this is not practical since the enhanced singularity of increasing order functions implies that the truncation order must itself be a function of $`L`$. An elegant alternative was presented in section. 3.3: the exponent can be written in terms of a finite set of analytic functions $`I_p(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2)`$ with coefficients $`\overline{r}_p`$ that are polynomials in $`L`$. Each function $`I_p(\mathrm{ln}Q^2/\overline{\mathrm{\Lambda }}^2)`$ corresponds to a single renormalon integral over a simple or a double pole. In this way running coupling effects are resummed first and the logs only later. As was shown in fig. 5, the convergence of this resummation method is quite remarkable.

The ambiguities of the perturbative result obtained by DGE can be used as a template for the non-perturbative power corrections. In particular, we saw that it is possible to identify the $`t`$ dependence associated with each renormalon singularity in the exponent. The leading power-correction for $`t\overline{\mathrm{\Lambda }}/Q`$ scales as $`1/Q`$. It corresponds to a simple renormalon pole at $`z=1/2`$, with an ambiguity which is proportional to $`1/t`$. This results, upon exponentiation, in a shift of the distribution, as predicted by . The physical origin of this correction is the sensitivity of the thrust to large-angle soft gluon emission. The same sensitivity gives rise to other renormalon singularities that are further away from the origin. The corresponding ambiguities scale as $`\overline{\mathrm{\Lambda }}^n/Q^n`$ and are therefore less important at large $`t`$. However, at small $`t`$ they behave as $`1/t^n`$, implying non-perturbative corrections of the form $`\lambda _n\overline{\mathrm{\Lambda }}^n/(Qt)^n`$. As emphasised by Korchemsky and Sterman, these corrections are all important in the distribution peak region. Since they exponentiate, it is natural to describe them by a non-perturbative shape-function. The parameters $`\lambda _n`$ then correspond to the central moments of this function.

Our calculation shows that the ambiguities $`\overline{\mathrm{\Lambda }}^n/Q^n`$ with $`n=2,4,\mathrm{}`$ vanish in the large $`\beta _0`$ limit. This strongly suggests that the parameters $`\lambda _n`$ for even $`n`$ are suppressed according to eq. (3.5). In addition, the appearance of a double renormalon pole at $`z=1`$ leads to a collinear power correction proportional to $`\lambda _2\overline{\mathrm{\Lambda }}^2/(Qt)^2`$, which has the opposite sign to the $`\lambda _2`$ term in the large-angle contribution. This leads to an almost complete cancellation of the overall effect of $`\lambda _2`$.

The immediate phenomenological implication of these findings is that by shifting the distribution one should be able to fit the data in a wider range (closer to the peak) than a priori expected: there is no new effect beyond that of the shift so long as $`\frac{\lambda _3}{\lambda _1}\frac{\overline{\mathrm{\Lambda }}^2}{Q^2t^2}1`$, instead of $`\frac{\lambda _2}{\lambda _1}\frac{\overline{\mathrm{\Lambda }}}{Qt}1`$. The phenomenological analysis indeed shows very successful shift-based fits, even for rather low $`t`$ values (this was first observed in ). Clearly, in the peak region itself, a shift is not enough, and a convolution with a shape-function must replace it.

Finally, our main conclusions from the phenomenological analysis can be summarised as follows: we obtained good fits based on the DGE result both with a shifted distribution and with a shape function. The two methods are fairly consistent with each other. This demonstrates that the data in the peak region can be used to extract $`\alpha _s`$ provided the appropriate resummation is performed and the power corrections are included. In the fits, the central value of $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})`$ ranges between $`0.108`$ and $`0.115`$. The most important factor which limits the determination of $`\alpha _s`$ is the absence of a NNLO calculation. In particular, terms that are explicitly suppressed by $`t`$ at NNLO are important, as can be deduced from the fact that the modification of the log: $`\mathrm{ln}1/t\mathrm{ln}(1/t1)`$ can change the extracted $`\alpha _s^{\overline{\text{MS}}}(\mathrm{M}_\mathrm{Z})`$ by $`4\%`$. There are also non-negligible differences ($`3\%`$) between different fitting procedures, such as dependence on the low $`t`$ cut, the use of a shift or a shape-function, and the parametric form of the latter. Other factors contributing to the uncertainty are discussed in sec. 5.4.

The correlation between the non-perturbative parameters $`\lambda _n`$ and $`\alpha _s`$ is strong. Consequently a quantitative discussion is restricted to the case where $`\alpha _s`$ is fixed. Then the shift parameter, or the first moment of the shape function, $`\lambda _1`$ is reasonably well determined (see fig. 7). Even when a modified log is used, the value of $`\lambda _1`$ as a function of $`\alpha _s`$ remains roughly the same. The sub-leading power-corrections are harder to constrain by the data. To determine the higher central moments of the shape-function it is useful to fix also $`\lambda _1`$. As an example, fixing $`\alpha _s`$ and $`\lambda _1`$ to their best fit values we find that $`\lambda _2`$ is small and that the fit is not affected much by fixing it to zero (see table 3). Thus our theoretical prediction that the effect of $`\lambda _2`$ is small can be consistent with the data. A clear preference of the fits for a shape-function which is not positive-definite was observed.

Acknowledgments

We would like to thank Stefano Catani, Gregory Korchemsky and Gavin Salam for very interesting and useful discussions and Otmar Biebel and Stefan Kluth for useful correspondence.

## Appendix A Numerical calculation of the NNLL part of the NLO coefficient

To calculate the NNLL part of the NLO coefficient we have used the EVENT2  program. Writing the thrust-distribution on the form

$`{\displaystyle \frac{1}{\sigma _0}}{\displaystyle \frac{d\sigma }{dL}}(Q^2,L)|_{\text{NLO}}`$ $`=`$ $`𝒜(L)a(Q^2)+(L)a(Q^2)^2,`$ (127)
where $`L=\mathrm{ln}(1/t)`$, $`a=\alpha _s^{\overline{\text{MS}}}/\pi `$, and $`𝒜(L)`$ is the analytically known leading order coefficient, the program calculates the next-to-leading order coefficient $`(L)`$ in the $`\overline{\text{MS}}`$ scheme using Monte Carlo integration. Note that the normalisation in Eq. (127) is to the Born cross-section $`\sigma _0`$ and not to the total cross-section $`\sigma =\left(1+\frac{3}{4}C_Fa(Q^2)+\mathrm{}\right)\sigma _0`$. The program is written such that the different colour factor parts can be calculated separately (this was done before using the EVENT program  in ) in the following way,

$`(L)`$ $`=`$ $`C_F^2_{C_F}(L)+C_FC_A_{C_A}(L)+C_FT_FN_F_{T_F}(L).`$
Since we are interested in the large $`L`$ (small $`t`$) region we decreased the invariant mass cut-off on any pair of partons, $`m_{ij}^2>\text{CUTOFF}Q^2`$, in the program to $`10^{14}`$ from the default $`10^8`$ in order not to restrict the phase-space for small $`t`$. To increase the sampling of the small $`t`$ region we also changed the parameters NPOW1 and NPOW2 to 4 instead of the default values $`2`$. The calculations presented below are based on $`10^9`$ events (samplings of the integral). We used a linear binning in $`L`$ and the bin size was chosen to be $`\mathrm{\Delta }L=0.5`$.

For each bin in $`L`$ the Monte Carlo integration gives the integrated cross-section. From these integrals we then subtract the known NLL parts (except the constant contributions which will be dealt with separately) and normalise to the bin size giving the following integrals for each bin,

$`I_i`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Delta }L}}{\displaystyle _{L_{\mathrm{min}}}^{L_{\mathrm{max}}}}\left(_i(L)_i^{\text{NLL}}(L)\right)𝑑L,`$
where $`i=C_F,C_A,T_F`$, and

$`_{C_F}^{\text{NLL}}(L)`$ $`=`$ $`2L^3{\displaystyle \frac{9}{2}}L^2+\left({\displaystyle \frac{13}{4}}\pi ^2\right)L,`$
$`_{C_A}^{\text{NLL}}(L)`$ $`=`$ $`{\displaystyle \frac{11}{4}}L^2+\left({\displaystyle \frac{\pi ^2}{6}}{\displaystyle \frac{169}{72}}\right)L,`$
$`_{T_F}^{\text{NLL}}(L)`$ $`=`$ $`L^2+{\displaystyle \frac{11}{18}}L.`$
The results for the different integrals are shown in fig. 12.

If it were not for sub-leading corrections, the integrals $`I_i`$ would coincide with the NNLL coefficients $`_{i,1}`$. In principle one could take into account sub-leading terms by fitting a more general form to the integrals $`I_i`$. We choose instead to estimate a systematic error from sub-leading corrections for each bin in $`L`$. The estimate is based on assuming that the most important form of sub-leading corrections has the same power of logs as the leading log multiplied by one factor of $`t`$. The size of the coefficient multiplying this term is then assumed to be the same as $`_{i,1}`$ estimated from the integral $`I_i`$. The systematic error for each bin is thus estimated by, e.g.

$`\mathrm{\Delta }I_{C_F}`$ $`=`$ $`{\displaystyle \frac{I_{C_F}}{\mathrm{\Delta }L}}{\displaystyle _{L_{\mathrm{min}}}^{L_{\mathrm{max}}}}L^3\mathrm{exp}(L)𝑑L.`$
These errors are illustrated as bands in fig. 12. When making the fit these systematic errors are added in quadrature to the statistical errors from the Monte Carlo integration.

From fig. 12 it is clear that for $`L\begin{array}{c}<\hfill \\ \hfill \end{array}7`$ the sub-leading corrections start to become important. It is also clear that for $`L\begin{array}{c}>\hfill \\ \hfill \end{array}12`$ the statistical errors start to increase significantlyFor even larger $`L`$ the phase-space limitation imposed by the invariant mass cut-off becomes clearly noticeable.. Thus we choose the fit range $`7.1<L<12.1`$. Fitting a constant in this range gives the following results,

$`_{C_F,1}`$ $`=`$ $`2.95\pm 0.27\pm 0.11,`$
$`_{C_A,1}`$ $`=`$ $`1.855\pm 0.081\pm 0.093,`$
$`_{T_F,1}`$ $`=`$ $`1.255\pm 0.022\pm 0.010.`$
where the first error is from the fit and the second from the maximal variation when varying the upper and lower limits for the fit by $`\pm 1`$.

From these results we can now get the NNLL term $`𝒢_{21}`$ in the order $`a^2`$ contribution to $`\mathrm{ln}R`$,

$$\mathrm{ln}R=\left[𝒢_{12}L^2+𝒢_{11}L\right]a+\left[𝒢_{23}L^3+𝒢_{22}L^2+𝒢_{21}L\right]a^2+\mathrm{}.$$

For the $`C_A`$ and $`T_F`$ parts these are just the same as the $`_{i,1}`$ but for the $`C_F`$ part one has to subtract the contribution coming from the expansion of the exponent and the different normalisation used in Eq. (127), $`(\pi ^23)/4`$. Adding the errors in quadrature and writing the result as

$`𝒢_{21}`$ $`=`$ $`𝒢_{\beta _0,1}C_F\beta _0+𝒢_{C_F,1}C_F^2+𝒢_{C_A,1}C_AC_F`$
we get

$`𝒢_{C_F,1}`$ $`=`$ $`1.23\pm 0.29,`$
$`𝒢_{C_A,1}`$ $`=`$ $`1.60\pm 0.14,`$
$`𝒢_{\beta _0,1}`$ $`=`$ $`3.766\pm 0.070.`$
These results can be compared with Eq. (4.2) where

$`𝒢_{C_F,1}`$ $``$ $`0.13,`$
$`𝒢_{C_A,1}`$ $``$ $`0.73,`$
$`𝒢_{\beta _0,1}`$ $`=`$ $`3.75,`$
in order to examine to what accuracy the DGE generates the sub-leading logs. The coefficient which is directly calculated in our approach is $`𝒢_{\beta _0,1}`$. In this case the agreement is very good. As expected, the other coefficients which are generated by the exponentiation do not agree. It is encouraging, however, that their values are smaller and that they have the correct signs. Thus, including the corresponding terms at NNLO and beyond is certainly an improvement.

For easy comparison with earlier results  we also write the result in the expansion

$$\mathrm{ln}R=\left[G_{12}L^2+G_{11}L\right]a/2+\left[G_{23}L^3+G_{22}L^2+G_{21}L\right]a^2/4+\mathrm{},$$

with $`G_{21}=G_{T_F,1}+G_{C_F,1}+G_{C_A,1}`$. Our results (for $`N_F=5`$) are then given by

$`G_{C_F,1}`$ $`=`$ $`9\pm 2,`$
$`G_{C_A,1}`$ $`=`$ $`30\pm 2,`$
$`G_{T_F,1}`$ $`=`$ $`16.7\pm 0.3.`$
The result for the sum $`G_{21}=22\pm 3`$ agrees within errors with earlier results, but the results for the colour factors $`C_A`$ and $`T_F`$ are significantly different . The reason for these discrepancies can probably be found in the different ranges in $`L`$ used for the fit, namely $`7.1<L<12.1`$ and $`3.1<L<6.7`$ respectively.