Soft-gluon effective coupling and cusp anomalous dimension

TL;DR

This paper extends the soft-gluon effective coupling in QCD beyond next-to-leading logarithms, establishing a relation with the cusp anomalous dimension and providing explicit expressions up to high orders.

Contribution

It introduces two all-order proposals for the soft-gluon effective coupling linked to the cusp anomalous dimension, expanding the theoretical framework in QCD resummation.

Findings

Derived explicit expressions for the couplings up to ${\cal O}(\as^2)$ in $d$ dimensions.

Confirmed previous ${\cal O}(\as^3)$ results for one coupling.

Established a relation between the couplings and the cusp anomalous dimension up to ${\cal O}(\as^4)$.

Abstract

We consider the extension of the CMW soft-gluon effective coupling in the context of soft-gluon resummation for QCD hard-scattering observables beyond the next-to-leading logarithmic accuracy. We present two proposals of a soft-gluon effective coupling that extend the CMW coupling to all perturbative orders in the \ms\ coupling $\as$. Although both effective couplings are well-defined in the physical four-dimensional space time, we examine their behaviour in $d=4 -2\ep$ space time dimensions. We uncover an all-order perturbative relation with the cusp anomalous dimension: the (four dimensional) cusp anomalous dimension is equal to the $d$-dimensional soft-gluon effective coupling at the conformal point $\ep=\beta(\as)$, where the $d$-dimensional QCD $\beta$-function, $\beta(\as) - \ep$, vanishes. We present the explicit expressions of the two soft-gluon couplings up to ${\cal O}(\as^2)$ in $d$ dimensions. In the four-dimensional case we compute the two soft couplings up to ${\cal O}(\as^3)$. For one of the two couplings, we confirm the ${\cal O}(\as^3)$ result previously presented by other authors. For the other coupling, we obtain the explicit relation with the cusp anomalous dimension up to ${\cal O}(\as^4)$. We comment on Casimir scaling at ${\cal O}(\as^4)$.

Equations

$$d{w}_i^{DL}=C_i\frac{{\alpha }_{\mathrm{S}}}{\pi }\frac{d\omega }{\omega }\frac{d{\theta }^2}{{\theta }^2}\simeq C_i\frac{{\alpha }_{\mathrm{S}}}{\pi }\frac{dz}{1-z}\frac{d{q}_T^2}{q_T^2},$$

$$C_i\frac{{\alpha }_{\mathrm{S}}}{\pi }\to {\mathcal{A}}_i^{CMW}({\alpha }_{\mathrm{S}}(q_T^2))=C_i\frac{{\alpha }_{\mathrm{S}}^{CMW}(q_T^2)}{\pi }=C_i\frac{{\alpha }_{\mathrm{S}}(q_T^2)}{\pi }\left(1+\frac{{\alpha }_{\mathrm{S}}(q_T^2)}{2\pi }K\right),$$

$$K=\left(\frac{67}{18}-\frac{{\pi }^2}{6}\right)C_A-\frac{5}{9}{n}_F,$$

$$P_{ii}({\alpha }_{\mathrm{S}};z)=\frac{1}{1-z}{A}_i({\alpha }_{\mathrm{S}})+\dots ,(z<1),$$

$$A_i({\alpha }_{\mathrm{S}})=\sum _{n=1}^{\infty }{\left(\frac{{\alpha }_{\mathrm{S}}}{\pi }\right)}^n{A}_i^{(n)}.$$

$$A_i^{(1)}$$

$$A_i^{(2)}$$

$$A_i^{(3)}$$

$$\displaystyle+\left(-\frac{55}{96}+\frac{1}{2}\zeta_{3}\right)C_{F}\,n_{F}-\frac{1}{108}n_{F}^{2}\;\Bigg{]}\;,$$

$${\alpha }_{\mathrm{S}}^u{\mu }_0^{2ϵ}{S}_{ϵ}={\alpha }_{\mathrm{S}}({\mu }_R^2){\mu }_R^{2ϵ}Z({\alpha }_{\mathrm{S}}({\mu }_R^2);ϵ),S_{ϵ}=(4\pi {)}^{ϵ}{e}^{-ϵ{\gamma }_E},$$

$$Z({\alpha }_{\mathrm{S}},ϵ)=1-{\alpha }_{\mathrm{S}}\frac{{\beta }_0}{ϵ}+{\alpha }_{\mathrm{S}}^2\left(\frac{{\beta }_0^2}{{ϵ}^2}-\frac{{\beta }_1}{2ϵ}\right)+\mathcal{O}({\alpha }_{\mathrm{S}}^3),$$

$$\beta ({\alpha }_{\mathrm{S}})=-{\beta }_0{\alpha }_{\mathrm{S}}-{\beta }_1{\alpha }_{\mathrm{S}}^2+\mathcal{O}({\alpha }_{\mathrm{S}}^3),$$

$$12\pi {\beta }_0=11C_A-2n_F,24{\pi }^2{\beta }_1=17C_A^2-5C_A{n}_F-3C_F{n}_F.$$

$${\tilde{\mathcal{A}}}_{T,i}({\alpha }_{\mathrm{S}}({\mu }^2);ϵ)$$

$${\tilde{\mathcal{A}}}_{0,i}({\alpha }_{\mathrm{S}}({\mu }^2);ϵ)$$

$${\mathcal{A}}_{T,i}({\alpha }_{\mathrm{S}})\equiv {\tilde{\mathcal{A}}}_{T,i}({\alpha }_{\mathrm{S}};ϵ=0),{\mathcal{A}}_{0,i}({\alpha }_{\mathrm{S}})\equiv {\tilde{\mathcal{A}}}_{0,i}({\alpha }_{\mathrm{S}};ϵ=0).$$

$${\tilde{\mathcal{A}}}_i({\alpha }_{\mathrm{S}};ϵ)=\sum _{n=1}^{\infty }{\left(\frac{{\alpha }_{\mathrm{S}}}{\pi }\right)}^n{\tilde{\mathcal{A}}}_i^{(n)}(ϵ),{\mathcal{A}}_i({\alpha }_{\mathrm{S}})=\sum _{n=1}^{\infty }{\left(\frac{{\alpha }_{\mathrm{S}}}{\pi }\right)}^n{\mathcal{A}}_i^{(n)}.$$

$${\tilde{\mathcal{A}}}_i^{(n)}(ϵ)={\mathcal{A}}_i^{(n)}+\sum _{k=1}^{\infty }{ϵ}^k{\tilde{\mathcal{A}}}_i^{(n;k)}.$$

$${\tilde{\mathcal{A}}}_{T,i}^{(1)}(ϵ)={\tilde{\mathcal{A}}}_{0,i}^{(1)}(ϵ)=C_{i}c(ϵ),$$

$$c(ϵ)\equiv \frac{e^{ϵ{\gamma }_E}}{\Gamma (1-ϵ)}=1-\frac{{\pi }^2}{12}{ϵ}^2-\frac{1}{3}{\zeta }_3{ϵ}^3+\mathcal{O}({ϵ}^4),$$

$${\mathcal{A}}_{T,i}^{(2)}={\mathcal{A}}_{0,i}^{(2)}=A_i^{(2)}.$$

$${\tilde{\mathcal{A}}}_{T,i}({\alpha }_{\mathrm{S}};ϵ=\beta ({\alpha }_{\mathrm{S}}))={\tilde{\mathcal{A}}}_{0,i}({\alpha }_{\mathrm{S}};ϵ=\beta ({\alpha }_{\mathrm{S}}))=A_i({\alpha }_{\mathrm{S}}).$$

$$A_i({\alpha }_{\mathrm{S}}({\mu }_F^2))=\frac{d}{d\ln {\mu }_F^2}{\mathcal{P}}_{ϵ}\left\{{\int }_0^{{\mu }_F^2}\frac{d{q}_T^2}{q_T^2}{\tilde{\mathcal{A}}}_i(({\alpha }_{\mathrm{S}}(q_T^2);ϵ)\right\},$$

$${\tilde{\mathcal{A}}}_{T,i}^{(2)}(ϵ)$$

$$h(ϵ)={\gamma }_E+\psi (1-ϵ)+2\psi (1+2ϵ)-2\psi (1+ϵ),$$

$${\tilde{\mathcal{A}}}_{0,i}^{(2)}(ϵ)$$

$$r(ϵ)=\frac{2}{1+ϵ}{}_3{F}_2(1,1,1-ϵ;2-2ϵ,2+ϵ;1)-\frac{1}{1-ϵ}{}_3{F}_2(1,1,1-ϵ;2-2ϵ,2-ϵ;1),$$

$${\tilde{\mathcal{A}}}_{T,i}^{(2)}(ϵ)$$

$${\tilde{\mathcal{A}}}_{0,i}^{(2)}(ϵ)$$

$${\mathcal{A}}_i^{(3)}=A_i^{(3)}-({\beta }_0\pi {)}^2{\tilde{\mathcal{A}}}_i^{(1;2)}+({\beta }_0\pi ){\tilde{\mathcal{A}}}_i^{(2;1)}.$$

$${\mathcal{A}}_{T,i}^{(3)}$$

$${\mathcal{A}}_{0,i}^{(3)}$$

Full text

ZU-TH 17/19

ICAS 42/19

**Soft-gluon effective coupling

and cusp anomalous dimension**

Stefano Catani(a), Daniel de Florian*(b)* and Massimiliano Grazzini*(c)***

*(a)*INFN, Sezione di Firenze and Dipartimento di Fisica e Astronomia,

Università di Firenze, I-50019 Sesto Fiorentino, Florence, Italy

(b) International Center for Advanced Studies (ICAS), ECyT-UNSAM,

Campus Miguelete, 25 de Mayo y Francia, (1650) Buenos Aires, Argentina

*(c)*Physik Institut, Universität Zürich, CH-8057 Zürich, Switzerland

Abstract

We consider the extension of the CMW soft-gluon effective coupling [1] in the context of soft-gluon resummation for QCD hard-scattering observables beyond the next-to-leading logarithmic accuracy. We present two proposals of a soft-gluon effective coupling that extend the CMW coupling to all perturbative orders in the ${\overline{\rm MS}}$ coupling $\alpha_{\mathrm{S}}$. Although both effective couplings are well-defined in the physical four-dimensional space time, we examine their behaviour in $d=4-2\epsilon$ space time dimensions. We uncover an all-order perturbative relation with the cusp anomalous dimension: the (four dimensional) cusp anomalous dimension is equal to the $d$-dimensional soft-gluon effective coupling at the conformal point $\epsilon=\beta(\alpha_{\mathrm{S}})$, where the $d$-dimensional QCD $\beta$-function, $\beta(\alpha_{\mathrm{S}})-\epsilon$, vanishes. We present the explicit expressions of the two soft-gluon couplings up to ${\cal O}(\alpha_{\mathrm{S}}^{2})$ in $d$ dimensions. In the four-dimensional case we compute the two soft couplings up to ${\cal O}(\alpha_{\mathrm{S}}^{3})$. For one of the two couplings, we confirm the ${\cal O}(\alpha_{\mathrm{S}}^{3})$ result previously presented by other authors. For the other coupling, we obtain the explicit relation with the cusp anomalous dimension up to ${\cal O}(\alpha_{\mathrm{S}}^{4})$. We comment on Casimir scaling at ${\cal O}(\alpha_{\mathrm{S}}^{4})$.

April 2019

A well known feature of QCD is that perturbative computations of hard-scattering processes are sensitive to soft-gluon effects. These effects manifest themselves in hard-scattering observables that are evaluated close to the exclusive boundary of the phase space. In such kinematical configurations, real-radiation contributions in the inclusive final state are strongly suppressed and they cannot balance virtual-radiation effects (which are always kinematically allowed). The unbalance leads to large logarithmic radiative corrections (the argument of the logarithms is the distance from the exclusive boundary). These contributions are often refereed to as logarithmically-enhanced terms of Sudakov type.

Real emission has a logarithmic spectrum for radiation that is soft and/or collinear to the direction of the observed hard jets (partons). This spectrum produces perturbative contributions that have at most two powers of logarithms for each power of the QCD coupling $\alpha_{\mathrm{S}}$. In the case of many observables, the double-logarithmic (DL) terms can be resummed to all orders in $\alpha_{\mathrm{S}}$ in exponentiated form. For most of these observables resummation can be extended in exponentiated form to higher (or arbitrary) logarithmic orders. This feature depends on the hard-scattering process and, especially, on the specific kinematical properties of the observable under consideration. In the case of observables that fulfil exponentiation, it is natural to device a resummed perturbative expansion by systematically organizing the exponent in classes of subsequent logarithmic accuracy: leading logarithmic (LL) terms, next-to-leading logarithmic (NLL) terms, next-to-next-to-leading logarithmic (NNLL) terms and so forth.

The explicit computation and resummation of Sudakov-type logarithms can be performed by using both traditional perturbative QCD methods and techniques based on Soft Collinear Effective Theory. Very many observables are nowadays known up to NLL accuracy (see, e.g., the reviews in Refs. [2, 3] and references therein for a large, though still incomplete, list of NLL results), and several observables are known to NNLL or higher logarithmic accuracy (a list of results can be found in Ref. [4]).

Direct inspection of NLL resummed results shows that they have a high degree of universality, with general structures and ingredients that have a ‘minimal’ dependence on the hard-scattering process and on the specific observable to be treated. Roughly speaking, the resummed logarithmic contributions are embodied in a ‘generalized’ Sudakov form factor and they are combined with ‘hard’ (non-logarithmic) factors, which are computable at fixed perturbative orders. The universality structure of NLL resummation is evident in the context of process-independent and observable-independent formulations of resummation that have been explicitly worked out [5, 6, 7] for large classes of (properly specified) observables. In particular, within such formulations, the Sudakov form factor is obtained by integration (over an observable-dependent phase space) of a universal kernel that is explicitly evaluated up to NLL accuracy. In view of these NNL results, progress is being carried out [8] to extend such observable-independent resummation program to NNLL accuracy. Obviously, an improved understanding of NNLL contributions and their possible universality features is also relevant in the context of resummed calculations for specific observables, independently of any observable-independent treatment.

A relevant feature of the NLL results that we have briefly recalled is that the ‘dominant’ (soft and collinear) part of the NLL kernel of the generalized Sudakov form factor is obtained simply and in complete form through the use of the QCD coupling $\alpha_{\mathrm{S}}^{CMW}$ [1] in the Catani–Marchesini–Webber (CMW) scheme (or bremsstrahlung scheme). The CMW coupling $\alpha_{\mathrm{S}}^{CMW}$ has the meaning of an effective (physical) coupling for inclusive radiation of soft and collinear gluons. The purpose of the present paper is to extend the definition of $\alpha_{\mathrm{S}}^{CMW}$ beyond NLL accuracy. A definition of such extension has been proposed in Ref. [4]. Moreover, the authors of Ref. [4] present the relation between the effective coupling and the customary ${\overline{\rm MS}}$ renormalized coupling up to ${\cal O}(\alpha_{\mathrm{S}}^{3})$, and they embody the effective coupling in the context of an explicit formulation of NNLL resummation for generic two-jet observables in $e^{+}e^{-}$ annihilation [8]. We think that there is no unique extension of $\alpha_{\mathrm{S}}^{CMW}$ beyond NLL accuracy. By unique extension, we mean an extension with the same universality features as those of $\alpha_{\mathrm{S}}^{CMW}$ at NLL accuracy (we postpone additional comments on this). Therefore, in the present paper, besides considering the definition of Ref. [4], we propose a different definition of the soft-gluon effective coupling, and we present some results for both definitions.

We recall that the CMW coupling plays a role in two other contexts directly or indirectly related to Sudakov resummation. The coupling $\alpha_{\mathrm{S}}^{CMW}$ can be used in Monte Carlo event generators (see, e.g., Ref. [9]) to improve the logarithmic accuracy of corresponding parton shower algorithms [1]. The dispersive approach to power-behaved terms in QCD hard processes [10, 11] uses $\alpha_{\mathrm{S}}^{CMW}$ to combine contributions from the low-momentum (non-perturbative) region with perturbative contributions at next-to-leading order (NLO) in the high-momentum region. The extension of $\alpha_{\mathrm{S}}^{CMW}$ beyond NLL accuracy can be useful also for higher-order studies within these two contexts.

The outline of the paper is as follows. We first recall (in a sketchy way) the role of the CMW coupling in NLL resummed calculations. Then we introduce two definitions of the soft-gluon effective coupling at arbitrary perturbative orders. We illustrate various perturbative results for both effective couplings, and we present some brief comments on their derivation (details about the derivation of the results will appear in a separate publication). The results include an all-order relation with the cusp anomalous dimension and perturbative expressions up to ${\cal O}(\alpha_{\mathrm{S}}^{4})$. Finally, we conclude with a summary and some general comments.

At the lowest perturbative order, the probability of radiation of a single soft gluon that is soft ($\omega\ll E$) and collinear ($\theta\ll 1$) to the direction of a massless hard parton is given by the well-known DL spectrum

[TABLE]

where $\omega$ is the energy of the soft gluon, $E$ is the energy of the radiating hard parton and $\theta$ is the gluon emission angle. To DL accuracy, $dw_{i}^{DL}$ can be equivalently expressed in terms of the longitudinal-momentum fraction $1-z$ ($1-z\simeq\omega/E$) and transverse momentum $q_{T}$ ($q_{T}\simeq\omega\theta$) of the soft gluon. The subscript $i$ denotes the type of radiating parton ($i=q,{\bar{q}},g$), and $C_{i}$ is the corresponding quadratic Casimir coefficient. We have $C_{i}=C_{F}$ if $i=q$ (or $i={\bar{q}})$) and $C_{i}=C_{A}$ if $i=g$, with $C_{F}=(N_{c}^{2}-1)/(2N_{c})$ and $C_{A}=N_{c}$ in $SU(N_{c})$ QCD with $N_{c}$ colours.

The integration of the spectrum in Eq. (1) over the observable-dependent phase space produces large DL terms (in the vicinity of the exclusive boundary) and infrared divergent contributions that are cancelled by one-loop virtual-radiation effects. In the case of Sudakov sensitive observables that fulfil exponentiation, DL resummation is achieved by simply using $dw_{i}^{DL}$ as integration kernel in the exponent of the observable-dependent Sudakov form factor.

The intensity of soft-gluon radiation in Eq. (1) is $C_{i}\alpha_{\mathrm{S}}/\pi$. The NLL resummation of the contributions from soft and collinear radiation is obtained (see, for instance, Eqs. (10), (12) and (26) in Ref. [6], or Eqs. (2.16) and (2.29) in Ref. [7]) by using the DL kernel of Eq. (1) and simply replacing the intensity of the soft-gluon coupling as follows

[TABLE]

where $\alpha_{\mathrm{S}}^{CMW}$ is the CMW coupling [1] and $\alpha_{\mathrm{S}}(\mu^{2})$ is the QCD running coupling at the renormalization scale $\mu$ in the ${\overline{\rm MS}}$ renormalization scheme. The value of the coefficient $K$ in Eq. (2) is ($n_{F}$ is the number of massless-quark flavours)

[TABLE]

as it turned out since early works on NLL resummation of several observables [12, 13, 14].

Two effects are embodied in the DL kernel $dw_{i}^{DL}$ through the replacement in Eq. (2). The QCD coupling $\alpha_{\mathrm{S}}$ is evaluated at the scale of the soft-gluon transverse momentum $q_{T}$ [15, 16]: this accounts for the resummation of the LL terms. The gluon coupling acquires a correction of ${\cal O}(\alpha_{\mathrm{S}}^{2})$ (which is controlled by the coefficient $K$ in the ${\overline{\rm MS}}$ renormalization scheme): this produces the resummation of NLL terms. Since the replacement takes place in the exponent of the Sudakov form factor, it is produced by the correlated radiation of soft partons (both two soft gluons and a soft $q{\bar{q}}$ pair), whereas the independent emission of soft gluons is taken into account through the exponentiation. We also note that ${{\mathscr{A}}}_{i}^{CMW}(\alpha_{\mathrm{S}})$ is an effective coupling at the inclusive level, since it is obtained by integrating over the momenta of the final-state correlated partons. The coupling ${{\mathscr{A}}}_{i}^{CMW}(\alpha_{\mathrm{S}})$ refers to radiation that is both soft and collinear. The Sudakov form factor includes other NLL terms due to soft wide-angle (i.e., non-collinear) radiation and hard (i.e., non-soft) collinear radiation: we postpone some comments on these terms.

Since the CMW Sudakov kernel refers to soft and collinear radiation, it can be viewed as obtained by considering the soft limit of multiple collinear radiation. In this respect it is natural to compare it with the DGLAP kernel [17] that controls the collinear evolution of the parton distribution functions (PDFs). In the soft limit, $z\rightarrow 1$, the flavour diagonal DGLAP kernel $P_{ii}(\alpha_{\mathrm{S}};z)$ ($1-z$ is the longitudinal-momentum fraction that is radiated in the final state) has the following behaviour [18]:

[TABLE]

where the dots on the right-hand side denote terms that are less singular than $(1-z)^{-1}$ (we have also neglected contact terms, proportional to $\delta(1-z)$, of virtual origin). The soft behaviour in Eq. (4) also applies to the collinear evolution of the parton fragmentation functions.

The perturbative function $A_{i}(\alpha_{\mathrm{S}})$ in Eq. (4) is usually called (light-like) cusp anomalous dimension, since it can also be related to the renormalization of cusp singularities of Wilson loops [19, 20]. In the context of our discussion, $A_{i}(\alpha_{\mathrm{S}})$ directly refers to the soft limit in Eq. (4), independently of any relations with Wilson loop renormalization. The perturbative expansion of $A_{i}(\alpha_{\mathrm{S}})$ reads

[TABLE]

where $\alpha_{\mathrm{S}}$ is the renormalized ${\overline{\rm MS}}$ coupling. The perturbative coefficients $A^{(1)}_{i},A^{(2)}_{i}$ [21, 22] and $A_{i}^{(3)}$ [17] are explicitly known. Using the ${\overline{\rm MS}}$ factorization scheme for PDFs and fragmentation functions, these coefficients are

[TABLE]

where $\zeta_{k}$ is the Riemann $\zeta$-function. The fourth-order coefficient $A^{(4)}_{i}$ is known in approximate numerical form [23, 24] (the calculation in full analytic form is under completion), and a first numerical estimate of $A_{q}^{(5)}$ has been presented recently [25]. By direct inspection of Eqs. (6)–(8) we note that the dependence on $i$ (the type of radiating parton) of the perturbative function $A_{i}(\alpha_{\mathrm{S}})$ is entirely specified up to ${\cal O}(\alpha_{\mathrm{S}}^{3})$ by the overall colour factor $C_{i}$. This overall dependence on $C_{i}$, which is customarily named as Casimir scaling relation, follows from the soft-parton origin of $A_{i}(\alpha_{\mathrm{S}})$ [26], and it is violated at higher perturbative orders [27], starting from ${\cal O}(\alpha_{\mathrm{S}}^{4})$.

From Eqs. (5)–(7) we see that, up to the second perturbative order, $A_{i}(\alpha_{\mathrm{S}})$ coincides with the CMW coupling ${{\mathscr{A}}}_{i}^{CMW}(\alpha_{\mathrm{S}})$ in Eqs. (2) and (3). One may be tempted to conclude that the cusp anomalous dimension provides a sensible definition of a physical (though effective) soft-gluon coupling beyond ${\cal O}(\alpha_{\mathrm{S}}^{2})$. The equivalence between $A_{i}(\alpha_{\mathrm{S}})$ and soft-gluon coupling, however, cannot hold in general. Indeed $A_{i}(\alpha_{\mathrm{S}})$ depends on the factorisation scheme of collinear singularities, while the physical coupling should not.

We add more comments on this point, since there are conceptual analogies (and differences) between the soft-collinear part of the Sudakov kernel and the soft limit of the DGLAP kernel. The DGLAP kernel is related to the probability of correlated emission of collinear partons with comparable values of transverse momenta (independent collinear emission is instead taken into account through the perturbative iteration of the kernel). To obtain the DGLAP kernel, the transverse momenta are integrated up to some value of the evolution (or factorization) scale. The transverse-momentum integral is collinear divergent in the low-momentum region: within the ${\overline{\rm MS}}$ factorization scheme, the divergences are handled by using dimensional regularization in $d=4-2\epsilon$ space time dimensions, and the DGLAP kernel is defined as the coefficient of the ensuing $1/\epsilon$ pole (see related comments after Eq. (22)). This is an unphysical procedure, although it is perfectly well defined for factorization purposes (a different factorization procedure would lead to a different DGLAP kernel). In contrast, the $q_{T}$ integration of the Sudakov kernel does not lead to collinear divergences since the low-$q_{T}$ region is ‘physically’ regularized by the definition of the measured observables. Nonetheless, the equality between the cusp anomalous dimension and the CMW coupling at ${\cal O}(\alpha_{\mathrm{S}}^{2})$ is not completely accidental, since at this perturbative order the ${\overline{\rm MS}}$ factorization procedure is equivalent to introduce a lower bound on the transverse momentum [1], which practically acts as the regularization procedure that can be implemented through the use of a collinear safe observable.

In the following we introduce the all-order definitions of two soft-gluon effective couplings, and we present some perturbative results. The results are obtained by regularizing ultraviolet and infrared divergences (which are encountered at intermediate stages of the calculations) through analytic continuation in $d=4-2\epsilon$ space time dimensions. Specifically, we use the customary scheme of conventional dimensional regularization (CDR). The QCD bare coupling $\alpha_{\mathrm{S}}^{u}$ and the renormalized running coupling $\alpha_{\mathrm{S}}(\mu_{R}^{2})$ in the ${\overline{\rm MS}}$ renormalization scheme are related by the following standard definition:

[TABLE]

where $\mu_{0}$ is the dimensional regularization scale, $\mu_{R}$ is the renormalization scale and $\gamma_{E}$ is the Euler number. The renormalization function $Z(\alpha_{\mathrm{S}};\epsilon)$ is

[TABLE]

where $\beta_{0}$ and $\beta_{1}$ are the first two perturbative coefficients of the QCD $\beta$-function $\beta(\alpha_{\mathrm{S}})$:

[TABLE]

[TABLE]

As we have already stated, an all-order definition of soft-gluon effective coupling has been given in Ref. [4]. We use the same starting point as in Ref. [4]. We consider a generic hard-scattering process that involves only two massless hard partons, which can be either a $q{\bar{q}}$ pair ($i=q$) or two gluons ($i=g$). We compute the probability for emitting a set of soft partons (soft gluons and soft $q{\bar{q}}$ pairs), and we consider the function $w_{i}(k;\epsilon)$ that gives the ‘probability’†††Note that this ‘probability’ is not positive definite since it refers to the correlation part of the total emission probability. of correlated emission (including the corresponding virtual corrections) of an arbitrary number of soft partons with total momentum $k$. This function is formally defined in Eq. (2.25) of Ref. [4], and it is called web function therein.

Contributions to $w_{i}(k;\epsilon)$ from virtual and real radiative corrections separately lead to ultraviolet and infrared divergences. However, the probability of correlated soft emission at fixed total momentum $k$ is a quantity that is infrared and collinear safe. Therefore, infrared singularities cancel in the computation of $w_{i}(k;\epsilon)$ and, after renormalization of $\alpha_{\mathrm{S}}$, the soft function $w_{i}(k;\epsilon)$ is finite in the physical four-dimensional limit $\epsilon\rightarrow 0$. For our subsequent purposes, we consider the general $d$-dimensional function $w_{i}(k;\epsilon)$, although it is well defined at $\epsilon=0$.

A relevant property of $w_{i}(k;\epsilon)$ is its invariance under longitudinal boosts along the direction of the momenta of the two hard partons in their centre–of–mass frame. It follows that $w_{i}(k;\epsilon)$ actually depends only on two kinematical variables: the transverse-momentum component $k_{T}$ of $k$ with respect to the direction of the radiating partons, and the transverse mass $m_{T}$ ($m_{T}^{2}=k_{T}^{2}+k^{2}$). We propose the definition of two different effective couplings, ${\widetilde{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}};\epsilon)$ and ${\widetilde{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}};\epsilon)$, which measure the intensity of inclusive soft-parton radiation. The definitions are

[TABLE]

where ${\widetilde{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}};\epsilon=0)$ corresponds‡‡‡The function $w_{i}(k;\epsilon)$ in Eqs. (13) and (14) and the web function in Eq. (2.25) of Ref. [4] are directly proportional, and the proportionality relation includes the overall factor $(k_{T}^{2})^{-\epsilon}$ that makes ${\widetilde{\mathscr{A}}}_{T,i}$ and ${\widetilde{\mathscr{A}}}_{0,i}$ dimensionless in any number $d$ of dimensions. to the soft coupling of Ref. [4].

The definitions in Eqs. (13) and (14) differ only in the kinematical variable that is kept fixed in the integration procedure over $k$: ${\widetilde{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}}(\mu^{2});\epsilon)$ is defined at fixed value $k_{T}=\mu$ of the transverse momentum, while ${\widetilde{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}}(\mu^{2});\epsilon)$ is defined at fixed value $m_{T}=\mu$ of the transverse mass. In the right-hand side of Eqs. (13) and (14), the factor $\mu^{2}$ is introduced for dimensional reasons (so that ${\widetilde{\mathscr{A}}}_{i}$ is dimensionless) and the factor $1/2$ takes into account the fact that the integration of $w_{i}(k;\epsilon)$ includes the angular regions where the soft momentum $k$ is collinear to the momentum of each of the two hard partons. In the definitions of Eqs. (13) and (14) the renormalization scale $\mu_{R}$ is set to the value $\mu_{R}=\mu$. Obviously, the soft couplings ${\widetilde{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}}(\mu^{2});\epsilon)$ and ${\widetilde{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}}(\mu^{2});\epsilon)$ are renormalization group invariant quantities, so that, at the perturbative level, they can equivalently be expressed in terms of the running coupling $\alpha_{\mathrm{S}}(\mu_{R}^{2})$ and the ratio $\mu^{2}/\mu_{R}^{2}$.

The integration over $k$ in Eqs. (13) and (14) is infrared and collinear safe, so that the limit $\epsilon\rightarrow 0$ is finite and well defined. Therefore, the soft-gluon effective couplings ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$ and ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ in the physical four-dimensional space time are simply

[TABLE]

Nonetheless we insist in using a $d$-dimensional definition of the soft-gluon coupling for a twofold (formal and practical) purpose. The formal aspects will be discussed below. At the practical level, the $d$-dimensional definition permits a direct application of the effective coupling in the context of hadron collisions, where Sudakov resummation can be sensitive to the PDFs of the colliding hadrons (and the related ${\overline{\rm MS}}$ factorization procedure in $d$ dimensions).

The coefficients of the perturbative expansion of ${\widetilde{\mathscr{A}}}_{i}$ and ${{\mathscr{A}}}_{i}$ are defined analogously to those in Eq. (5):

[TABLE]

The $\epsilon$-expansion at the $n$-th perturbative order is denoted as follows

[TABLE]

To make explicit the definition of the overall normalization of ${\widetilde{\mathscr{A}}}_{i}(\alpha_{\mathrm{S}};\epsilon)$ (and $w_{i}(k;\epsilon)$), we report the expression of the lowest-order contribution:

[TABLE]

where

[TABLE]

and $\Gamma(z)$ is the Euler $\Gamma$-function. We note that the two soft couplings ${\widetilde{\mathscr{A}}}_{T,i}^{(1)}$ and ${\widetilde{\mathscr{A}}}_{0,i}^{(1)}$ are exactly equal at the lowest perturbative order. This equality simply follows from the fact that the lowest-order contribution to $w_{i}(k;\epsilon)$ is proportional to $\delta(k^{2})=\delta(m_{T}^{2}-k_{T}^{2})$. We also note (see Eq. (19)) that the $\epsilon$ dependence of ${\widetilde{\mathscr{A}}}_{i}^{(1)}(\epsilon)$ starts at ${\cal O}(\epsilon^{2})$ (i.e., the coefficient ${\widetilde{\mathscr{A}}}_{i}^{(1;1)}$ at ${\cal O}(\epsilon)$ vanishes). This mild $\epsilon$ dependence is of entirely ‘kinematical’ origin (it arises from the $d$-dimensional phase space), since (due to helicity conservation) the dynamics of soft-gluon radiation does not produce any $\epsilon$ dependence at the lowest perturbative order.

We anticipate (see below) that, in the physical four-dimensional space time, both soft couplings in Eqs. (13) and (14) are equal to the CMW coupling ${{\mathscr{A}}}_{i}^{CMW}$ up to ${\cal O}(\alpha_{\mathrm{S}}^{2})$. Therefore, we have

[TABLE]

One of the main results of this paper is the following all-order relation between the cusp anomalous dimension $A_{i}(\alpha_{\mathrm{S}})$ and the soft-gluon couplings:

[TABLE]

This relation can be derived in differents ways. A procedure that we have used consists in considering threshold resummation [29, 14, 30] for the production of high-mass systems in hadron collisions. The threshold resummed cross section is related to the evolution of the PDFs in the soft limit (see Eq. (4)). We have applied both soft couplings in Eqs. (13) and (14) to the computation of the threshold resummed cross section and we have obtained the result in Eq. (21).

The relation in Eq. (21) can be rewritten in the following form:

[TABLE]

where ${\widetilde{\mathscr{A}}}_{i}$ is equivalently ${\widetilde{\mathscr{A}}}_{T,i}$ or ${\widetilde{\mathscr{A}}}_{0,i}$, and ${\cal P_{\epsilon}}$ is the projection operator [21] that extracts the $\epsilon$ poles (in ${\overline{\rm MS}}$ form) of the function of $\alpha_{\mathrm{S}}(\mu_{F}^{2})$ and $\epsilon$ in the curly bracket. The equivalence between Eqs. (21) and (22) can be proven by using some $d$-dimensional technicalities. We would like to point out that the relation between ${\widetilde{\mathscr{A}}}_{i}$ and $A_{i}$ as expressed in the form of Eq. (22) is in direct correspondence with our previous qualitative discussion about the relation between the Sudakov kernel and the DGLAP kernel. The soft coupling ${\widetilde{\mathscr{A}}}_{i}((\alpha_{\mathrm{S}}(q_{T}^{2});\epsilon)$ gives the intensity of the spectrum of correlated soft and collinear emission of partons with total transverse momentum $q_{T}$. In the right-hand side of Eq. (22), the $q_{T}$ spectrum is integrated over the region from $q_{T}=0$ up to some value of the factorization scale $\mu_{F}$. Following the ${\overline{\rm MS}}$ factorization procedure, the $\epsilon$ poles that arise from the $d$-dimensional regularization of the collinear singularities in the region around $q_{T}\simeq 0$ are then extracted to obtain (actually, to define) the intensity $A_{i}(\alpha_{\mathrm{S}}(\mu_{F}^{2}))$ of soft radiation in the DGLAP kernel (i.e., the cusp anomalous dimension in Eq. (4)).

Equation (21) relates§§§An equality between the soft anomalous dimension (which is related to the cusp anomalous dimension) and the $d$-dimensional rapidity anomalous dimension at the point $\epsilon=\beta(\alpha_{\mathrm{S}})$ is presented in Ref. [28]. the cusp anomalous dimension to the $d$-dimensional soft-gluon coupling at the conformal point $\epsilon=\beta(\alpha_{\mathrm{S}})$, where the $d$-dimensional QCD $\beta$-function $\beta(\alpha_{\mathrm{S}})-\epsilon$ vanishes. The relation (21) is not specific of QCD, and it also applies to other gauge theories. In particular, in the case of ${\cal N}=4$ maximally supersymmetric Yang–Mills theory we have $\beta(\alpha_{\mathrm{S}})=0$ and, therefore, the cusp anomalous dimension coincides with the physical (four-dimensional) soft-gluon coupling: ${{\mathscr{A}}}_{T}(\alpha_{\mathrm{S}})={{\mathscr{A}}}_{0}(\alpha_{\mathrm{S}})=A(\alpha_{\mathrm{S}})$.

According to Eq. (21), there is a non-trivial interplay between the perturbative dependence of the cusp anomalous dimension and the $d$-dimensional dependence of the soft-gluon coupling. In particular, since the $\epsilon$-dependence of ${\widetilde{\mathscr{A}}}_{i}^{(1)}(\epsilon)$ starts at ${\cal O}(\epsilon^{2})$ (see Eqs. (18) and (19)), Eq. (21) directly implies the equivalence up to ${\cal O}(\alpha_{\mathrm{S}}^{2})$ (see Eq. (20)) between the cusp anomalous dimension and the four-dimensional soft-gluon coupling (or the CMW coupling). As we have already recalled, this equivalence is not completely accidental [1] and, at the purely technical level, it can be viewed as a consequence of the mild ${\cal O}(\epsilon^{2})$ dependence in Eqs. (18) and (19).

The relation in Eq. (21) also states that the two $d$-dimensional soft couplings, ${\widetilde{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}};\epsilon)$ and ${\widetilde{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}};\epsilon)$, become equal by setting $\epsilon=\beta(\alpha_{\mathrm{S}})$. Starting from ${\cal O}(\alpha_{\mathrm{S}}^{2})$ (see Eqs. (23) and (25) below), the $\epsilon$-dependence of the two soft couplings is very different. In view of this, we find it remarkable that such a different $\epsilon$-dependence conspires to make the coupling equal at $\epsilon=\beta(\alpha_{\mathrm{S}})$. Incidentally, such a different $\epsilon$-dependence and the relation (21) imply that the two four-dimensional soft couplings, ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$ and ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$, inevitably differ starting from ${\cal O}(\alpha_{\mathrm{S}}^{3})$. Moreover, the difference ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})-{{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ is necessarily due to perturbative contributions that are proportional to the coefficients, $\beta_{0},\beta_{1}$ and so forth, of the QCD $\beta$-function.

In addition to be interesting for its intrinsic structure, the relation in Eq. (21) can be exploited for several different purposes. It can be used to crosscheck explicit perturbative computations of $A_{i}(\alpha_{\mathrm{S}})$ and ${\widetilde{\mathscr{A}}}_{i}(\alpha_{\mathrm{S}};\epsilon)$. Once one the the three functions $A_{i},{\widetilde{\mathscr{A}}}_{T,i}$ and ${\widetilde{\mathscr{A}}}_{0,i}$ is known at some perturbative order, Eq. (21) can exploited to extract information on the other two functions (in the following we explicitly make this use of Eq. (21)). The relation (21) can also be used to obtain the cusp anomalous dimension $A_{i}(\alpha_{\mathrm{S}})$ through the $d$-dimensional perturbative calculation of one of the two soft couplings ${\widetilde{\mathscr{A}}}_{i}(\alpha_{\mathrm{S}};\epsilon)$.

We have computed the soft function $w_{i}(k;\epsilon)$ at ${\cal O}(\alpha_{\mathrm{S}}^{2})$ by combining the one-loop correction to single soft-gluon radiation [31] with the $d$-dimensional integration of double soft-parton radiation at the tree level [32]. Then, using Eqs. (13) and (14), we have computed the soft-gluon effective couplings in $d$ dimensions at ${\cal O}(\alpha_{\mathrm{S}}^{2})$, and we obtain the following results [33]. In the case of ${\widetilde{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}};\epsilon)$ we find

[TABLE]

where

[TABLE]

and $\psi(1+z)=\frac{d\ln\Gamma(1+z)}{dz}$. In the case of ${\widetilde{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}};\epsilon)$ we find

[TABLE]

where

[TABLE]

and ${}_{3}F_{2}(\alpha,\beta,\gamma;\delta,\rho;z)$ is the generalized hypergeometric function of the variable $z$.

The $\epsilon$-expansion up to ${\cal O}(\epsilon^{2})$ of the second-order expressions in Eqs. (23) and (25) gives

[TABLE]

[TABLE]

where $A_{i}^{(2)}$ is given in Eq. (7). From these equations we see that the $\epsilon$ dependence of the two soft couplings ${\widetilde{\mathscr{A}}}_{T,i}^{(2)}(\epsilon)$ and ${\widetilde{\mathscr{A}}}_{0,i}^{(2)}(\epsilon)$ is already different at ${\cal O}(\epsilon)$. We also see that the limit $\epsilon\rightarrow 0$ of our explicit calculation at ${\cal O}(\alpha_{\mathrm{S}}^{2})$ leads to the equality in Eq. (20) between the two soft couplings and the CMW coupling. At the computational level the equality ${{\mathscr{A}}}_{T,i}^{(2)}={{\mathscr{A}}}_{0,i}^{(2)}$ originates as follows. Since ${\widetilde{\mathscr{A}}}_{T,i}^{(1)}(\epsilon)={\widetilde{\mathscr{A}}}_{0,i}^{(1)}(\epsilon)$, the value of ${\widetilde{\mathscr{A}}}_{i}^{(2)}(\epsilon)$ at $\epsilon=0$ is determined by the behaviour of the soft function $w_{i}(k;\epsilon)$ in the region where $k^{2}\simeq 0$. In this region we have $m_{T}^{2}\simeq k_{T}^{2}$ and, therefore, the difference between the right-hand side of Eqs. (13) and (14) (and, hence, between the two soft couplings) is not effective.

We now present our computation of the third-order coefficients ${{\mathscr{A}}}_{T,i}^{(3)}$ and ${{\mathscr{A}}}_{0,i}^{(3)}$ of both four-dimensional soft couplings. To this purpose we use Eq. (11) and we perturbatively expand Eq. (21) in terms of the coefficients ${\widetilde{\mathscr{A}}}_{i}^{(n;k)}$ that are defined in Eq. (17). We obtain

[TABLE]

This relation applies to both soft couplings ${\widetilde{\mathscr{A}}}_{T,i}$ and ${\widetilde{\mathscr{A}}}_{0,i}$ (we have omitted the corresponding subscripts $T$ and [math]), and we have also used ${\widetilde{\mathscr{A}}}_{i}^{(1;1)}=0$ (see Eqs. (18) and (19)). Since we have determined ${\widetilde{\mathscr{A}}}_{i}^{(1)}(\epsilon)$ and ${\widetilde{\mathscr{A}}}_{i}^{(2)}(\epsilon)$ to all orders in the $\epsilon$-expansion, the explicit values of the coefficients ${\widetilde{\mathscr{A}}}_{i}^{(1;2)}$ and ${\widetilde{\mathscr{A}}}_{i}^{(2;1)}$ can be directly read from Eqs. (18), (19), (S0.Ex4) and (S0.Ex5). Inserting these coefficients in Eq. (29) we can explicitly relate ${{\mathscr{A}}}_{i}^{(3)}$ to the coefficient $A_{i}^{(3)}$ (see Eq. (8)) of the cusp anomalous dimension. We obtain the following results:

[TABLE]

[TABLE]

Our result in Eq. (30) for the third-order coefficient of the soft-gluon coupling ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$ agrees with the corresponding result presented in Ref. [4] (see Eqs. (3.9) and (3.11b) therein).

Using the value of $A_{i}^{(3)}$ in Eq. (8), the results in Eqs. (30) and (31) explicitly relate the four-dimensional (physical) soft-gluon effective couplings ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$ and ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ with the ${\overline{\rm MS}}$ renormalized coupling $\alpha_{\mathrm{S}}$ up to ${\cal O}(\alpha_{\mathrm{S}}^{3})$. This relation generalizes the ${\cal O}(\alpha_{\mathrm{S}}^{2})$ CMW relation in Eq. (2) to the third order, and it can be used to construct the Sudakov kernel for soft-gluon resummation of infrared and collinear safe observables at NNLL accuracy [4].

In the case of the soft-gluon coupling ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ we have also computed its relation with the ${\overline{\rm MS}}$ coupling at ${\cal O}(\alpha_{\mathrm{S}}^{4})$. More precisely, we obtain an explicit relation between ${{\mathscr{A}}}_{0,i}^{(4)}$ and the corresponding coefficient $A_{i}^{(4)}$ of the cusp anomalous dimension. We find

[TABLE]

This fourth-order result can be used for applications to soft-gluon resummed calculations of infrared and collinear safe observables at the next-to-next-to-next-to-leading logarithmic (N3LL) accuracy.

Knowing the result in Eq. (S0.Ex6) and exploiting the relation in Eq. (21), we can also explicitly determine the third-order coefficient ${\widetilde{\mathscr{A}}}_{0,i}^{(3)}(\epsilon)$ of the $d$-dimensional soft coupling at ${\cal O}(\epsilon)$. To illustrate the procedure, we perturbatively expand Eq. (S0.Ex6) in terms of the coefficients ${\widetilde{\mathscr{A}}}_{i}^{(n;k)}$ of Eq. (17), and we obtain

[TABLE]

where we have used ${\widetilde{\mathscr{A}}}_{0,i}^{(1;1)}=0$. The explicit coefficients ${\widetilde{\mathscr{A}}}_{0,i}^{(1;2)}$ and ${\widetilde{\mathscr{A}}}_{0,i}^{(1;3)}$ at the first order and ${\widetilde{\mathscr{A}}}_{0,i}^{(2;1)}$ and ${\widetilde{\mathscr{A}}}_{0,i}^{(2;2)}$ at the second order can be read from Eqs. (18), (19) and Eq. (S0.Ex5), respectively. Therefore, by comparing Eqs. (S0.Ex6) and (33) we obtain

[TABLE]

with the explicit result

[TABLE]

We comment on our derivation of the result in Eq. (S0.Ex6). The soft-gluon effective coupling ${\widetilde{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}};\epsilon)$ is particularly suitable in the context of threshold resummation [29, 14, 30] for the production of colourless high-mass systems in hadron collisions. The threshold resummed cross section for these processes is presently known in explicit form up to N3LL accuracy [34]-[43]. We have applied ${\widetilde{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}};\epsilon)$ to threshold resummation and, exploiting the known N3LL results [42], we obtain Eq. (S0.Ex6).

The result in Eq. (S0.Ex6) relates the fourth-order perturbative term ${{\mathscr{A}}}_{0,i}^{(4)}$ of the soft coupling ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ to the corresponding term $A_{i}^{(4)}$ of the cusp anomalous dimension $A_{i}(\alpha_{\mathrm{S}})$. Since $A_{i}^{(4)}$ is not fully known in analytic form, we add some comments on the fourth-order results.

We have examined the colour structure of soft multiparton radiation from two hard partons at ${\cal O}(\alpha_{\mathrm{S}}^{4})$ and, consequently, we can obtain the general colour structure of the soft function $w_{i}(k;\epsilon)$ or, equivalently (due to Eqs. (13) and (14)), the colour structure of the soft coupling. We write this structure in the following form:

[TABLE]

where $N_{i}$ is the dimension of the colour representation of the hard parton $i$ ($N_{i}=N_{A}=N_{c}^{2}-1$ if $i=g$, and $N_{i}=N_{F}=N_{c}$ if $i=q,{\bar{q}}$), and $d_{xy}^{(4)}$ are the quartic Casimir invariants (we use the normalization of $d_{xy}^{(4)}$ as in Eqs. (2.6)–(2.10) of Ref. [24]). The entire dependence of ${\widetilde{\mathscr{A}}}_{i}^{(4)}(\epsilon)$ (for both couplings ${\widetilde{\mathscr{A}}}_{T,i}^{(4)}(\epsilon)$ and ${\widetilde{\mathscr{A}}}_{0,i}^{(4)}(\epsilon)$) on the colour of the hard parton $i$ is embodied in the Casimir dependent factors that we have explicitly written in the right-hand side of Eq. (36). The ‘quartic’ (${{\widetilde{\mathscr{A}}}}^{(4)}_{[4A]}(\epsilon)$ and $\;{{\widetilde{\mathscr{A}}}}^{(4)}_{[4F]}$) and ‘quadratic’ (${{\widetilde{\mathscr{A}}}}^{(4)}_{[2]}$) coefficients do not depend on the type of radiating parton $i$. In particular, ${{\widetilde{\mathscr{A}}}}^{(4)}_{[4A]}(\epsilon)$ and $\;{{\widetilde{\mathscr{A}}}}^{(4)}_{[4F]}$ are colour blind (they do not depend on $N_{c}$ and $n_{F}$). The coefficient ${{\widetilde{\mathscr{A}}}}^{(4)}_{[2]}$ still depends on $N_{c}$ and $n_{F}$, and this dependence involves all the colour structures that appear in the curly bracket of Eq. (S0.Ex6) plus an additional term with colour factor $C_{F}^{2}n_{F}$.

The presence in Eq. (36) of the quartic Casimir invariants violates Casimir scaling (i.e., the proportionality relation ${\widetilde{\mathscr{A}}}_{i}\propto C_{i}$). Nonetheless ${\widetilde{\mathscr{A}}}_{i}^{(4)}$ in Eq. (36) still fulfils a form of generalized Casimir scaling (in terms of three colour coefficients that depend on $i$) since ${{\widetilde{\mathscr{A}}}}^{(4)}_{[2]},{{\widetilde{\mathscr{A}}}}^{(4)}_{[4A]}$ and ${{\widetilde{\mathscr{A}}}}^{(4)}_{[4F]}$ do not depend on the hard parton $i$.

Setting $\epsilon=0$ in Eq. (36) and using Eq. (21), we obtain the colour structure of the four-dimensional soft coupling ${{\mathscr{A}}}_{0,i}^{(4)}$ (or, analogously¶¶¶The expression in Eq. (37) is equally valid for the soft coupling ${{\mathscr{A}}}_{T,i}^{(4)}$ through the replacement ${{\mathscr{A}}}^{(4)}_{0[2]}\rightarrow{{\mathscr{A}}}^{(4)}_{T[2]}$., ${{\mathscr{A}}}_{T,i}^{(4)}$) and of the cusp anomalous dimension $A_{i}^{(4)}$:

[TABLE]

[TABLE]

where, analogously to Eq. (36), the full dependence on the colour of the hard parton $i$ is entirely controlled by the Casimir dependent coefficients $C_{i},d_{Ai}^{(4)}/N_{i}$ and $d_{Fi}^{(4)}/N_{i}$.

We note that, to obtain Eqs. (37) and (38) from Eq. (36), we have exploited Eq. (21) and the property that the difference ${{\mathscr{A}}}_{0,i}^{(4)}-A_{i}^{(4)}$ fulfils Casimir scaling (see Eq. (33)), since the perturbative terms ${\widetilde{\mathscr{A}}}_{0,i}^{(n)}(\epsilon)$ with $n=1,2,3$ fulfil Casimir scaling. In particular, in Eqs. (37) and (38) we have set ${{\widetilde{\mathscr{A}}}}^{(4)}_{0[2]}(\epsilon=0)\equiv{{\mathscr{A}}}^{(4)}_{0[2]}$, then we have related the ‘quadratic’ coefficients of the soft coupling (${{\mathscr{A}}}^{(4)}_{0[2]})$) and of the cusp anomalous dimension ($A^{(4)}_{[2]}$) through Casimir scaling:

[TABLE]

and, finally, we have derived and implemented the following equalities

[TABLE]

between the ‘quartic’ coefficients of the soft coupling (${{\mathscr{A}}}^{(4)}_{0[4A]},{{\mathscr{A}}}^{(4)}_{0[4F]}$) and of the cusp anomalous dimension ($A^{(4)}_{[4A]},A^{(4)}_{[4F]}$). Our result in Eq. (S0.Ex6) is fully consistent with the Casimir scaling relation in Eq. (39).

We note that the generalized Casimir scaling of the soft coupling in Eq. (36) and the relation in Eq. (21) necessarily imply the same scaling for the cusp anomalous dimension in Eq. (38). The generalized Casimir scaling of the cusp anomalous dimension has been conjectured and verified to good numerical accuracy in Ref. [24]. We also note that at the fourth order the DGLAP kernel $P_{gg}(\alpha_{\mathrm{S}};z)$ includes a contribution with the quartic Casimir invariant $d_{FF}^{(4)}$, which is absent in $A_{g}^{(4)}$ of Eq. (38). Such contribution to $P_{gg}^{(4)}(\alpha_{\mathrm{S}};z)$ vanishes in the soft limit, consistently with the approximate numerical result of Ref. [24].

The fourth-order term $A_{i}^{(4)}$ of the cusp anomalous dimension is not yet known in full analytic form, although it is known with good numerical accuracy. The analytic results, which regard the coefficients of various colour factors, have been obtained by using different methods: the computation of the soft limit of the DGLAP kernel [44, 45, 23], the fourth-order evaluation of form factors [46, 47, 48], the cusp renormalization of Wilson loops [49] (as we have previously observed, the relation (21) leads to another method to compute $A_{i}(\alpha_{\mathrm{S}})$ through the evaluation of the $d$-dimensional soft coupling ${\widetilde{\mathscr{A}}}_{i}(\alpha_{\mathrm{S}};\epsilon)$). In particular, the ‘quartic’ coefficient $A^{(4)}_{[4F]}$ in Eqs. (37) and (38) has been computed very recently [47, 48]. The coefficients of the remaining colour factor contributions to $A_{i}^{(4)}$ have been evaluated in approximate numerical form [24].

The quantitative effect on the soft coupling ${{\mathscr{A}}}_{0,i}^{(4)}$ of the present numerical uncertainty of $A_{i}^{(4)}$ is very small, since the quantitative value of ${{\mathscr{A}}}_{0,i}^{(4)}$ turns out to be dominated by the contribution ${{\mathscr{A}}}_{0,i}^{(4)}-A_{i}^{(4)}$ that we have explicitly computed in Eq. (S0.Ex6). To see this, we write

[TABLE]

The term $\Delta_{i}$ depends on $n_{F}$. Using ${{\mathscr{A}}}_{0,i}^{(4)}-A_{i}^{(4)}$ from Eq. (S0.Ex6) and $A_{i}^{(4)}$ from Ref. [24] and setting $n_{F}=5$ (with $N_{c}=3$) we obtain

[TABLE]

where the numbers in brackets indicate the numerical uncertainty (due to $A_{i}^{(4)}$ [24]) of the preceding digit. Similar quantitative results are obtained for $n_{F}=3,4$. The term $\Delta_{i}$ turns out to contribute to ${{\mathscr{A}}}_{0,i}^{(4)}$ at the level of few percents, so that a small uncertainty on $A_{i}^{(4)}$ leads to a very small uncertainty on ${{\mathscr{A}}}_{0,i}^{(4)}$.

As observed in Ref. [24], due to the actual values of the ‘quartic’ coefficients $A^{(4)}_{[4A]}$ and $A^{(4)}_{[4F]}$ in Eq. (38), numerical Casimir scaling is completely broken in the fourth-order term $A_{i}^{(4)}$ of the cusp anomalous dimension. However, due to the smallness of $\Delta_{i}$, the soft coupling ${{\mathscr{A}}}_{0,i}^{(4)}$ still fulfils numerical Casimir scaling (${{\mathscr{A}}}_{0,i}^{(4)}\propto C_{i}$) modulo corrections at the few percent level.

We report the numerical value of the soft coupling ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ with $N_{c}=3$ up to ${\cal O}(\alpha_{\mathrm{S}}^{4})$. Using ${{\mathscr{A}}}_{0,i}^{(3)}$ from Eq. (31), ${{\mathscr{A}}}_{0,i}^{(4)}-A_{i}^{(4)}$ from Eq. (S0.Ex6), $A_{i}^{(4)}$ (with its numerical uncertainty) from Ref. [24] and setting $n_{F}=5$, we have

[TABLE]

The perturbative expansion in Eq. (43) can be compared with the corresponding perturbative expansion of the cusp anomalous dimension in Eq. (4.4) of Ref. [24]. From the comparison we can see that the third-order∥∥∥For comparison with the value $1.7157$ in Eq. (43), we note that the numerical value of the third-order coefficient (see Eq. (30)) of the soft coupling ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$ is 0.49121 . and fourth-order numerical coefficients in ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ are sizeably larger than those in $A_{i}(\alpha_{\mathrm{S}})$. Nonetheless the perturbative expansion of ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ is still numerically well behaved. We also see that the violation of Casimir scaling in the fourth-order term of ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ is numerically at the 4% level.

We add some general (though brief) comments on the soft-gluon effective coupling and Sudakov resummation.

The resummation procedure of logarithmic contributions of Sudakov type requires proper kinematical approximations of the phase space for multiparton final-state radiation. Such approximations are specific of the physical observables under consideration. As a consequence, the use of one or the other of the two soft-gluon couplings ${\widetilde{\mathscr{A}}}_{T,i}$ and ${\widetilde{\mathscr{A}}}_{0,i}$ can be more appropriate depending on the observables. The two soft couplings can alternatively (or equivalently) be used for the resummation treatment of different classes of observables. Some observables can also require a combined use of both soft couplings. In Ref. [4] the soft-gluon coupling ${{\mathscr{A}}}_{T,i}$ has been explicitly applied to the resummation of a wide class of observables. As we have previously mentioned, the soft-gluon coupling ${\widetilde{\mathscr{A}}}_{0,i}$ is particularly suitable in the context of threshold resummation and related observables, and its application to other classes of observables can be investigated.

The soft-gluon coupling ${\widetilde{\mathscr{A}}}_{i}$ controls the intensity of the spectrum of soft and collinear radiation in the Sudakov kernel. The Sudakov kernel has other dynamical components that, roughly speaking, are due to soft non-collinear (i.e., wide-angle) radiation and hard (i.e., non-soft) collinear radiation. Both components have to be included in a resummed calculation (see, e.g., Refs. [5, 6, 7] at NLL accuracy and Refs. [8, 4] at NNLL accuracy), and their inclusion has to be properly performed (i.e., properly matched) according to the soft coupling (either ${\widetilde{\mathscr{A}}}_{T,i}$ or ${\widetilde{\mathscr{A}}}_{0,i}$) that is specifically used in the soft-collinear component. However, we note that, at a given fixed perturbative order (say, $\alpha_{\mathrm{S}}^{n}$) in the Sudakov kernel, the soft-collinear component is logarithmically enhanced (by at least one power of log) with respect to the two other components. Therefore, the Sudakov kernel at NkLL accuracy requires the knowledge of the soft coupling ${{\mathscr{A}}}_{i}(\alpha_{\mathrm{S}})$ up to ${\cal O}(\alpha_{\mathrm{S}}^{k+1})$ and the computation of the other components up to ${\cal O}(\alpha_{\mathrm{S}}^{k})$ (i.e., one order lower than the soft coupling). For instance, to achieve NNLL accuracy in the Sudakov kernel, the third-order results in Eqs. (30) and (31) for the soft coupling have to be combined with the calculation at ${\cal O}(\alpha_{\mathrm{S}}^{2})$ of the other dynamical components.

A final comment regards the process dependence of the Sudakov kernel. The soft couplings in Eqs. (13) and (14) are computed by considering soft-parton radiation from two hard partons in a colour singlet configuration. Soft-gluon radiation in processes that involve several hard partons is definitely more complex than in the case of two hard-parton processes. This complex structure of soft-gluon radiation has to be taken properly into account. However, this does not affect the soft coupling ${\widetilde{\mathscr{A}}}_{i}$, since ${\widetilde{\mathscr{A}}}_{i}$ measures the intensity of radiation that is both soft and collinear to parton $i$. The complex structure of soft radiation in multiparton hard scattering only affects the soft wide-angle component of the Sudakov kernel (see, e.g., Refs. [5, 6, 7] at NLL accuracy).

We conclude the paper with a brief summary of its content. We have considered the all-order extension of the CMW effective coupling in the context of soft-gluon resummation beyond NLL accuracy. We have argued that there is no unique all-order extension, namely, no extension that shares all the universality (i.e., observable-independent) features of the CMW coupling at ${\cal O}(\alpha_{\mathrm{S}}^{2})$. Starting from the emission probability of an arbitrary number of soft partons, we have introduced the definition in $d=4-2\epsilon$ space-time dimensions of two effective couplings, ${\widetilde{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}};\epsilon)$ and ${\widetilde{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}};\epsilon)$, which measure the intensity of the inclusive spectrum for soft and collinear radiation from a massless hard parton $i$ ($i=q,{\bar{q}},g$). We have found that, to all perturbative orders, the two soft couplings are equal if they are evaluated at the $d$-dimensional point $\epsilon=\beta(\alpha_{\mathrm{S}})$, and they coincide with the (four-dimensional) cusp anomalous dimension $A_{i}(\alpha_{\mathrm{S}})$. The limit $\epsilon\rightarrow 0$ is smooth and it can be used to define the four-dimensional (‘physical’) couplings ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$ and ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$. The coupling ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$ has originally been defined in Ref. [4], and its explicit relation with $\alpha_{\mathrm{S}}$ up to ${\cal O}(\alpha_{\mathrm{S}}^{3})$ has been presented therein. We have computed both couplings, ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$ and ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$, up to ${\cal O}(\alpha_{\mathrm{S}}^{3})$ and, in the case of ${{\mathscr{A}}}_{T,i}(\alpha_{\mathrm{S}})$, our independent calculation confirm the result in Ref. [4]. In the case of ${\widetilde{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}};\epsilon)$ we are able to compute its third-order contribution up to ${\cal O}(\epsilon)$ and, in the four-dimensional limit, we obtain an explicit relation at ${\cal O}(\alpha_{\mathrm{S}}^{4})$ between ${{\mathscr{A}}}_{0,i}(\alpha_{\mathrm{S}})$ and the cusp anomalous dimension $A_{i}(\alpha_{\mathrm{S}})$. Moreover, we have presented the explicit $d$-dimensional results (e.g., to all orders in the $\epsilon$ expansion) for both soft couplings up to ${\cal O}(\alpha_{\mathrm{S}}^{2})$.

Acknowledgements. We would like to thank Bryan Webber for comments on the manuscript. This work is supported in part by the Swiss National Science Foundation (SNF) under contracts 200020_169041 and IZSAZ2_173357, by MINCyT under contract SUIZ/17/05, by Conicet and by ANPCyT.