Next Article in Journal
Universal Constants and Natural Systems of Units in a Spacetime of Arbitrary Dimension
Previous Article in Journal
Resonant Production of an Ultrarelativistic Electron–Positron Pair at the Gamma Quantum Scattering by a Field of the X-ray Pulsar
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Soft Anomalous Dimensions and Resummation in QCD

Department of Physics, Kennesaw State University, Kennesaw, GA 30144, USA
Universe 2020, 6(10), 165; https://doi.org/10.3390/universe6100165
Submission received: 23 August 2020 / Revised: 20 September 2020 / Accepted: 29 September 2020 / Published: 1 October 2020
(This article belongs to the Section High Energy Nuclear and Particle Physics)

Abstract

:
I discuss and review soft anomalous dimensions in QCD that describe soft-gluon threshold resummation for a wide range of hard-scattering processes. The factorization properties of the cross section in moment space and renormalization-group evolution are implemented to derive a general form for differential resummed cross sections. Detailed expressions are given for the soft anomalous dimensions at one, two, and three loops, including some new results, for a large number of partonic processes involving top quarks, electroweak bosons, Higgs bosons, and other particles in the standard model and beyond.

1. Introduction

This review discusses soft anomalous dimensions that control soft-gluon threshold resummation in QCD. Resummation follows from factorization properties of the cross section [1,2,3,4,5,6,7,8,9,10,11] and it provides a formalism for calculating contributions to higher-order corrections that are theoretically important and numerically significant. These soft-gluon contributions take the form of logarithms of a variable that is a measure of the available energy for additional radiation in a process.
Beyond leading logarithms, resummation is crucially dependent on the color exchange in the hard scattering. One-loop calculations of soft anomalous dimensions are necessary to achieve next-to-leading-logarithm (NLL) accuracy, while two-loop calculations are needed for next-to-next-to-leading-logarithm (NNLL) accuracy. The current state of the art are three-loop calculations which are needed for next-to-next-to-next-to-leading-logarithm (N 3 LL) accuracy.
Soft-gluon resummation is particularly relevant near the partonic threshold and it has been used for a large number of hard-scattering processes. In fact, soft-gluon resummation is often relevant even far from the threshold for many standard model (SM) and beyond the standard model (BSM) processes. There are a vast number of results in the literature that have been provided using various approaches, schemes, gauges, and definitions. Therefore, in addition to reviewing past results, it is useful to provide a comprehensive and unified treatment using a common formalism and notation for all these results. In this review all results are shown using the standard moment-space resummation formalism in QCD for, in general, differential cross sections in the MS ¯ scheme and in Feynman gauge. The expressions for the resummed cross section are shown in single-particle-inclusive (1PI) kinematics, but the soft anomalous dimensions are the same in other kinematics choices, such as pair-invariant-mass (PIM) kinematics.
This paper is meant as a focused review of theoretical work on resummation and as a compendium of results on soft anomalous dimensions and other quantities in the resummed expressions. It is not a review of phenomenological papers with applications of resummations; there are many hundreds of such papers, and a recent review for many processes involving top quarks can be found in reference [12]. Furthermore, this is a review on standard moment-space QCD resummation; it is not a review of related work using alternative approaches such as soft-collinear effective theory (SCET) or non-relativistic QCD (NRQCD). A comparative review of resummation approaches in QCD and SCET was given in reference [13]. A main goal of this paper is to provide common terminology and notation for the large number of results in the literature for a large variety of SM and BSM processes, in the hope that it will be useful in comparing past results and in future applications.
In Section 2 we provide the formalism for soft-gluon resummation based on factorization and renormalization-group evolution. Section 3 presents fixed-order expansions of resummed cross sections at next-to-leading order (NLO), next-to-next-to-leading order (NNLO), and next-to-next-to-next-to- leading order (N 3 LO). Section 4 has results for one, two, and three loops for the cusp anomalous dimension for the separate cases of two massive lines, or two massless lines, or one massive and one massless lines. The soft anomalous dimensions for many processes with trivial color structures are discussed in Section 5, for high- p T electroweak-boson production and related processes in Section 6, for single-top production and related processes in Section 7, for top-antitop production and related processes in Section 8, for jet production and related processes in Section 9, and for several 2 3 processes in Section 10. A concluding summary is given in Section 11.

2. Soft-Gluon Resummation

In this section we briefly review the moment-space QCD soft-gluon resummation formalism. We discuss the factorization and refactorization of the cross section, the renormalization-group evolution (RGE), the eikonal approximation, soft anomalous dimensions, and the resummed cross section. For simplicity we discuss 2 2 processes but also explain how this generalizes to 2 n processes [11]. For specificity we choose 1PI kinematics but also discuss modifications for PIM kinematics.

2.1. Factorization, RGE, and Resummation

The factorized form of the—in general—differential cross section d σ A B 12 in hadronic collisions for the process A B 12 is
d σ A B 12 = a , b d x a d x b ϕ a / A ( x a , μ F ) ϕ b / B ( x b , μ F ) d σ ^ a b 12 ( μ F , μ R ) ,
where μ F is the factorization scale, μ R is the renormalization scale, ϕ a / A ( ϕ b / B ) are parton distribution functions (pdf) for parton a (b) in hadron A (B), and d σ ^ a b 12 is the differential hard-scattering partonic cross section.
We consider partonic processes a b 12 with 4-momenta p a + p b p 1 + p 2 , and define the usual kinematical variables s = ( p a + p b ) 2 , t = ( p a p 1 ) 2 , and u = ( p b p 1 ) 2 . In 1PI kinematics we choose, without loss of generality, particle 1 as the observed particle. We also define the threshold variable s 4 = s + t + u m 1 2 m 2 2 where the masses m 1 and m 2 can be zero or finite. As we approach partonic threshold, with vanishing energy for additional radiation, we have s 4 0 . If an additional gluon with momentum p g is emitted in the final state, then we can equivalently write s 4 = ( p 2 + p g ) 2 m 2 2 , so as p g goes to 0 (i.e., we have a soft gluon), we again see that s 4 describes the extra energy in the soft emission and that s 4 0 . We note that we can extend our formulas to the general case of multi-particle final states, i.e., 2 n processes [11], by replacing m 2 2 with ( p 2 + + p n ) 2 in the expressions.
With the incoming partons a and b arising from hadrons (e.g., protons/antiprotons) A and B, we define the hadron-level variables S = ( p A + p B ) 2 , T = ( p A p 1 ) 2 , U = ( p B p 1 ) 2 , and S 4 = S + T + U m 1 2 m 2 2 . By writing p a = x a p A and p b = x b p B , where x a and x b are the fractions of the momenta carried by, respectively, partons a and b in hadrons A and B, we have the relations s = x a x b S , t = x a T + ( 1 x a ) m 1 2 , and u = x b U + ( 1 x b ) m 1 2 . Using the above relations, we find that
S 4 S = s 4 s ( 1 x a ) u m 2 2 s ( 1 x b ) t m 2 2 s + ( 1 x a ) ( 1 x b ) m 1 2 m 2 2 s .
The last term in the above equation can be ignored near threshold, in the limit x a 1 and x b 1 , since it involves the product ( 1 x a ) ( 1 x b ) .
We next discuss the factorization of the cross section in integral transform space [3,5,8]. We define Laplace transforms (shown with a tilde) of the partonic cross section as d σ ^ ˜ a b 12 ( N ) = 0 s ( d s 4 / s ) e N s 4 / s d σ ^ a b 12 ( s 4 ) , with transform variable N, and we also define the transforms of the pdf as ϕ ˜ ( N ) = 0 1 e N ( 1 x ) ϕ ( x ) d x . We note that under transforms the logarithms of s 4 in the perturbative series produce logarithms of N, and we will show that the latter exponentiate.
We then consider the parton–parton cross section d σ a b 12 , which has the same form as Equation (1) but with incoming partons instead of hadrons [3,4,5,6,7,8]:
d σ a b 12 ( S 4 ) = d x a d x b ϕ a / a ( x a ) ϕ b / b ( x b ) d σ ^ a b 12 ( s 4 ) ,
noting that the leading power as s 4 0 comes entirely from flavor-diagonal distributions ϕ a / a and ϕ b / b [5,6], and we define its transform as
d σ ˜ a b 12 ( N ) = 0 S d S 4 S e N S 4 / S d σ a b 12 ( S 4 ) .
Using Equation (2) (without the last term, which vanishes near threshold), we can rewrite the transform of the parton–parton cross section as
d σ ˜ a b 12 ( N ) = 0 1 d x a e N a ( 1 x a ) ϕ a / a ( x a ) 0 1 d x b e N b ( 1 x b ) ϕ b / b ( x b ) 0 s d s 4 s e N s 4 / s d σ ^ a b 12 ( s 4 ) = ϕ ˜ a / a ( N a ) ϕ ˜ b / b ( N b ) d σ ^ ˜ a b 12 ( N ) ,
where N a = N ( m 2 2 u ) / s and N b = N ( m 2 2 t ) / s in 1PI kinematics (while N a = N b = N in the corresponding formula in PIM kinematics).
Next, we introduce a refactorization of the cross section in terms of new functions H a b 12 , S a b 12 , ψ a / a , ψ b / b , J 1 , and J 2 [3,4,5,6,7,8]. The process-dependent hard function H a b 12 is purely short-distance, nonradiative, and infrared safe, and it comprises contributions from the amplitude and its complex conjugate. The soft function S a b 12 describes the emission of noncollinear soft gluons and is also process-dependent. The coupling of soft gluons to the partons in the hard scattering processes is described by eikonal (Wilson) lines as ordered exponentials of the gauge field. The hard and soft functions are in general Hermitian matrices in color-exchange space, and the refactorized cross section involves the trace of their product. The functions ψ a / a and ψ b / b differ from the pdf ϕ a / a and ϕ b / b , and they describe the dynamics of collinear emission from the incoming partons [1,3,4,5,6,7,8]. The functions J 1 and J 2 describe collinear emission from any final-state massless colored particles, and are absent otherwise. The refactorized form of the cross section [3,5,6,8] is then
d σ a b 12 = d w a d w b d w 1 d w 2 d w S ψ a / a ( w a ) ψ b / b ( w b ) J 1 ( w 1 ) J 2 ( w 2 ) × tr H a b 12 α s ( μ R ) S a b 12 w S s μ F δ S 4 S + w a ( u m 2 2 ) s + w b ( t m 2 2 ) s w S w 1 w 2
where the ws are dimensionless weights, with w a and w b for ψ a / a and ψ b / b , respectively, w 1 and w 2 for J 1 and J 2 , respectively, and w S for S a b 12 . The argument in the delta function of Equation (6) arises from the recasting of Equation (2) in terms of the new weights, as
S 4 S = ( 1 x a ) u m 2 2 s ( 1 x b ) t m 2 2 s + s 4 s = w a ( u m 2 2 ) s w b ( t m 2 2 ) s + w S + w 1 + w 2 .
We note that w a 1 x a and w b 1 x b because they refer to different functions.
After taking a transform of Equation (6), we have
d σ ˜ a b 12 ( N ) = 0 1 d w a e N a w a ψ a / a ( w a ) 0 1 d w b e N b w b ψ b / b ( w b ) 0 1 d w 1 e N w 1 J 1 ( w 1 ) 0 1 d w 2 e N w 2 J 2 ( w 2 ) × tr H a b 12 α s ( μ R ) 0 1 d w s e N w s S a b 12 w s s μ F = ψ ˜ a / a ( N a ) ψ ˜ b / b ( N b ) J ˜ 1 ( N ) J ˜ 2 ( N ) tr H a b 12 α s ( μ R ) S ˜ a b 12 s N μ F .
All N-dependence has now been absorbed into the functions S ˜ , ψ ˜ , and J ˜ .
By comparing Equations (5) and (8), we find an expression for the hard-scattering partonic cross section in transform space:
d σ ^ ˜ a b 12 ( N ) = ψ ˜ a ( N a ) ψ ˜ b ( N b ) J ˜ 1 ( N ) J ˜ 2 ( N ) ϕ ˜ a / a ( N a ) ϕ ˜ b / b ( N b ) tr H a b 12 α s ( μ R ) S ˜ a b 12 s N μ F .
We resum the N-dependence of the soft matrix S ˜ a b 12 via renormalization-group evolution [3,5],
S ˜ a b 12 b = Z a b 12 S ˜ a b 12 Z a b 12
where S ˜ a b 12 b is the bare quantity and Z a b 12 is a (in general, matrix of) renormalization constant(s). Thus, we have the renormalization group equation for S ˜ a b 12 ,
μ R d S ˜ a b 12 d μ R = μ R μ R + β ( g s , ϵ ) g s S ˜ a b 12 = Γ S a b 12 S ˜ a b 12 S ˜ a b 12 Γ S a b 12
where Γ S a b 12 is the soft anomalous dimension (matrix), g s 2 = 4 π α s , and β ( g s , ϵ ) = g s ϵ / 2 + β ( g s ) in 4 ϵ dimensions where β ( g s ) = μ R d g s / d μ R is the QCD beta function. We can also define the beta function in an alternative form in terms of α s as
β ( α s ) = d ln α s d ln μ R 2 = n = 0 β n α s 4 π n + 1 ,
with β 0 = ( 11 C A 2 n f ) / 3 [14,15], where C A = N c with N c the number of colors, and n f is the number of light quark flavors, β 1 = 34 C A 2 / 3 2 C F n f 10 C A n f / 3 [16,17,18], where C F = ( N c 2 1 ) / ( 2 N c ) , and
β 2 = 2857 54 C A 3 + C F 2 205 18 C F C A 1415 54 C A 2 n f + 11 9 C F + 79 54 C A n f 2
see [19,20]. Additionally, β 3 [21] is
β 3 = C A 4 150 , 653 486 44 9 ζ 3 + C A 3 n f 39 , 143 162 + 68 3 ζ 3 + C A 2 C F n f 7073 486 328 9 ζ 3 + C A C F 2 n f 2102 27 + 176 9 ζ 3 + 23 C F 3 n f + C A 2 n f 2 3965 162 + 56 9 ζ 3 + C F 2 n f 2 338 27 176 9 ζ 3 + C A C F n f 2 4288 243 + 112 9 ζ 3 + 53 243 C A n f 3 + 154 243 C F n f 3 + d A a b c d d A a b c d N A 80 9 + 704 3 ζ 3 + d F a b c d d A a b c d N A n f 512 9 1664 3 ζ 3 + d F a b c d d F a b c d N A n f 2 704 9 + 512 3 ζ 3
where ζ 3 = 1.202056903 , N A = N c 2 1 , d A a b c d d A a b c d / N A = N c 2 ( N c 2 + 36 ) / 24 , d F a b c d d A a b c d / N A = N c ( N c 2 + 6 ) / 48 , and d F a b c d d F a b c d / N A = ( N c 4 6 N c 2 + 18 ) / ( 96 N c 2 ) . The result for β 4 is given by a long expression in reference [22], and for QCD with N c = 3 it has the approximate numerical value β 4 537 , 148 186 , 162 n f + 17567.8 n f 2 231.28 n f 3 1.8425 n f 4 .

2.2. Eikonal Approximation and Soft Anomalous Dimensions

From Equation (11) we see that the evolution of the soft function is controlled by the soft anomalous dimension, Γ S a b 12 , which is also in general a matrix. In calculating Γ S a b 12 , we use the eikonal approximation, where the Feynman rules for diagrams with soft gluon emission simplify. For example in the emission of a soft gluon with four-momentum k μ from a quark with final four-momentum p μ , we have the simplification
u ¯ ( p ) ( i g s T c ) γ μ i ( p / + k / + m ) ( p + k ) 2 m 2 + i ϵ u ¯ ( p ) g s T c γ μ p / + m 2 p · k + i ϵ = u ¯ ( p ) g s T c v μ v · k + i ϵ
where u ¯ is a Dirac spinor, T c are the generators of SU(3), and p μ = ( s / 2 ) v μ , with v μ a four-velocity. Thus, a typical one-loop diagram involving eikonal lines i and j is of the form
g s 2 d n k ( 2 π ) n v i μ ( v i · k + i ϵ ) ( i ) g μ ν k 2 v j ν ( v j · k + i ϵ ) ,
which contributes a UV pole ( 1 / ϵ ) ( α s / π ) θ coth θ where θ = ln [ ( v i · v j + ( v i · v j ) 2 v i 2 v j 2 ) / v i 2 v j 2 ] is the cusp angle that we will discuss further in Section 4. For massless eikonal lines this expression for the UV pole simplifies to ( 1 / ϵ ) ( α s / π ) ln ( 2 v i · v j / v i 2 v j 2 ) .
The calculation of Γ S a b 12 requires determining the coefficients of the ultraviolet poles of relevant eikonal diagrams [3,5,7,23]. The counterterms for S ˜ a b 12 are the ultraviolet divergent coefficients times the basis color tensors. If there are m color tensors, then the counterterms are
S L a b 12 = I = 1 m c I a b 12 Z I L a b 12 ,
for the corrections to the color tensor c L a b 12 , where Z I L a b 12 denotes the I L matrix element of the renormalization matrix in Equation (10).
Then
Γ S a b 12 = d Z a b 12 d ln μ R Z a b 12 1 = β ( g s , ϵ ) Z a b 12 g s Z a b 12 1 ,
is the soft anomalous dimension matrix that controls the evolution of the soft function S a b 12 via Equation (11). In dimensional regularization Z a b 12 has 1 / ϵ poles. When expanding Z a b 12 in powers of the strong coupling, its I L matrix element is
Z I L a b 12 = δ I L + α s π Z I L a b 12 ( 1 ) + O ( α s 2 ) ,
and since Z a b 12 ( 1 ) has a 1 / ϵ pole while β ( g s , ϵ ) includes a g s ϵ / 2 term in dimensional regularization, we find that Γ S a b 12 is given at one loop simply by the negative of the residue of Z a b 12 .

2.3. Resummed cross Section

The resummed differential cross section in transform space is derived from the renormalization-group evolution of the N-dependent functions in Equation (9), i.e., ψ ˜ / ϕ ˜ , J ˜ , and of course S ˜ a b 12 . We find [3,5,8,12]
d σ ^ a b 12 resum ( N ) = exp i = a , b E i ( N i ) exp i = a , b 2 μ F s d μ μ γ i / i N ˜ i , α s ( μ ) exp j = f . s . q , g E j ( N ) × tr H a b 12 α s ( s ) P ¯ exp s s / N ˜ d μ μ Γ S a b 12 α s ( μ ) × S ˜ a b 12 α s ( s N ˜ ) P exp s s / N ˜ d μ μ Γ S a b 12 α s ( μ ) ,
where the symbols P ( P ¯ ) denote path-ordering in the same (reverse) sense as the integration variable μ , and N ˜ = N e γ E with γ E the Euler constant.
The first exponential in Equation (20) resums soft and collinear contributions from the incoming partons [1,2]. We have
E i ( N i ) = 0 1 d z z N i 1 1 1 z 1 ( 1 z ) 2 d λ λ A i α s ( λ s ) + D i α s ( ( 1 z ) 2 s ) ,
with A i = n = 1 ( α s / π ) n A i ( n ) . Here A i ( 1 ) = C i which is C F for a quark or antiquark and C A for a gluon, A i ( 2 ) = C i K / 2 where K = C A ( 67 / 18 ζ 2 ) 5 n f / 9 [24] with ζ 2 = π 2 / 6 , A i ( 3 ) is given by [25]
A i ( 3 ) = C i C A 2 245 96 67 36 ζ 2 + 11 24 ζ 3 + 11 8 ζ 4 + C F n f 55 96 + ζ 3 2 + C A n f 209 432 + 5 18 ζ 2 7 12 ζ 3 n f 2 108 ,
where ζ 4 = π 4 / 90 , and A i ( 4 ) is given by [26,27]
A i ( 4 ) = C i C A 3 42 , 139 10 , 368 5525 1296 ζ 2 + 1309 432 ζ 3 + 451 64 ζ 4 11 24 ζ 2 ζ 3 451 288 ζ 5 ζ 3 2 16 313 96 ζ 6 + C F 2 n f 143 576 + 37 48 ζ 3 5 4 ζ 5 + C F C A n f 17 , 033 10 , 368 + 55 96 ζ 2 + 29 18 ζ 3 11 16 ζ 4 1 2 ζ 2 ζ 3 + 5 8 ζ 5 + C A 2 n f 24 , 137 20 , 736 + 635 648 ζ 2 361 108 ζ 3 11 48 ζ 4 + 7 12 ζ 2 ζ 3 + 131 144 ζ 5 + C F n f 2 299 2592 5 18 ζ 3 + ζ 4 8 + C A n f 2 923 20 , 736 19 648 ζ 2 + 35 108 ζ 3 7 48 ζ 4 + n f 3 1 648 + ζ 3 108 + d i a b c d d A a b c d N R i ζ 2 2 + ζ 3 6 + 55 12 ζ 5 3 2 ζ 3 2 31 8 ζ 6 + d i a b c d d F a b c d N R i n f ζ 2 ζ 3 3 5 3 ζ 5
where ζ 5 = 1.036927755 , ζ 6 = π 6 / 945 , and where d i a b c d is d F a b c d for a quark or antiquark and d A a b c d for a gluon, while N R i is N F = N c for a quark or antiquark and N A = N c 2 1 for a gluon.
Additionally, D i = n = 1 ( α s / π ) n D i ( n ) , with D i ( 1 ) = 0 , D i ( 2 ) given by [4]
D i ( 2 ) = C i C A 101 54 + 11 6 ζ 2 + 7 4 ζ 3 + n f 7 27 ζ 2 3 ,
and D i ( 3 ) given by [28]
D i ( 3 ) = C i C A 2 297 , 029 46 , 656 + 6139 648 ζ 2 + 2509 216 ζ 3 187 48 ζ 4 11 12 ζ 2 ζ 3 3 ζ 5 + C A n f 31 , 313 23 , 328 1837 648 ζ 2 155 72 ζ 3 + 23 24 ζ 4 + C F n f 1711 1728 ζ 2 4 19 36 ζ 3 ζ 4 4 + n f 2 29 729 + 5 27 ζ 2 + 5 54 ζ 3 .
We note that our expression for the resummed cross section is for hadron–hadron scattering but can be easily adapted for hadron–lepton scattering (the sums over a, b in the first two exponentials reduce to one term) and lepton-lepton scattering (the sums vanish).
In the second exponential in Equation (20), involving the factorization scale, γ i / i is the moment-space anomalous dimension of the MS ¯ density ϕ i / i [29,30,31,32,33], γ i / i = A i ln N ˜ i + γ i with parton anomalous dimensions γ i = n = 1 ( α s / π ) n γ i ( n ) , where γ q ( 1 ) = 3 C F / 4 , γ g ( 1 ) = β 0 / 4 ,
γ q ( 2 ) = C F 2 3 32 3 4 ζ 2 + 3 2 ζ 3 + C F C A 17 96 + 11 12 ζ 2 3 4 ζ 3 C F n f 1 48 + ζ 2 6 ,
and
γ g ( 2 ) = C A 2 2 3 + 3 4 ζ 3 n f C F 8 + C A 6 .
Additionally, γ q ( 3 ) is given by [25]
γ q ( 3 ) = C F 3 29 128 + 9 32 ζ 2 + 17 16 ζ 3 + 9 4 ζ 4 1 2 ζ 2 ζ 3 15 4 ζ 5 + C F 2 n f 23 64 + 5 48 ζ 2 17 24 ζ 3 + 29 48 ζ 4 + C F n f 2 17 576 + 5 108 ζ 2 ζ 3 36 + C F 2 C A 151 256 205 96 ζ 2 + 211 48 ζ 3 247 96 ζ 4 + 1 4 ζ 2 ζ 3 + 15 8 ζ 5 + C F C A 2 1657 2304 + 281 108 ζ 2 97 36 ζ 3 5 64 ζ 4 + 5 8 ζ 5 + C F C A n f 5 16 167 216 ζ 2 + 25 72 ζ 3 + ζ 4 32 .
The third exponential in Equation (20) describes soft and collinear radiation from the final-state (f.s.) particles [1,2,6,8]. The exponent vanishes for colorless particles and for massive particles, i.e., E j = 0 in those cases. In the case of hadron production, where a quark or gluon hadronizes into the observed hadron, then this exponent is the same as the initial-state exponent [10], i.e., E 1 ( N ) = E 1 ( N ) plus the same term in the exponent for the scale dependence. For heavy jets there is a different expression [6]. For all other cases with final-state massless quarks or gluons—the majority of the cases studied—the exponent has the expression
E j ( N ) = 0 1 d z z N 1 1 1 z ( 1 z ) 2 1 z d λ λ A j α s λ s + B j α s ( ( 1 z ) s ) + D j α s ( ( 1 z ) 2 s ) ,
where B j = n = 1 ( α s / π ) n B j ( n ) with B q ( 1 ) = 3 C F / 4 and B g ( 1 ) = β 0 / 4 . Additionally, for B q ( 2 ) [4] and B g ( 2 ) (cf. [34]), we have
B q ( 2 ) = C F 2 3 32 + 3 4 ζ 2 3 2 ζ 3 + C F C A 57 32 11 12 ζ 2 + 3 4 ζ 3 + C F n f 5 16 + ζ 2 6 ,
and
B g ( 2 ) = C A 2 1025 432 3 4 ζ 3 + 79 108 C A n f + C F n f 8 5 108 n f 2 ,
while B q ( 3 ) and B g ( 3 ) are given by (cf. [34])
B q ( 3 ) = C F 3 29 128 9 32 ζ 2 17 16 ζ 3 9 4 ζ 4 + ζ 2 ζ 3 2 + 15 4 ζ 5 + C F 2 n f 77 128 17 32 ζ 2 + 7 12 ζ 3 + ζ 4 4 + C F C A n f 455 216 + 199 216 ζ 2 + 29 36 ζ 3 5 12 ζ 4 + C F C A 2 5599 864 2831 864 ζ 2 + 5 9 ζ 3 + 211 96 ζ 4 5 8 ζ 5 + C F 2 C A 23 32 + 287 64 ζ 2 89 24 ζ 3 17 8 ζ 4 ζ 2 ζ 3 4 15 8 ζ 5 + C F n f 2 127 864 11 216 ζ 2 ζ 3 9
and
B g ( 3 ) = C A 3 299 , 341 31 , 104 + 1307 576 ζ 2 523 144 ζ 3 275 96 ζ 4 + ζ 2 ζ 3 4 + 5 4 ζ 5 + C A 2 n f 41 , 453 10 , 368 39 32 ζ 2 + 59 36 ζ 3 + 25 48 ζ 4 + C A C F n f 191 128 11 12 ζ 3 C F 2 n f 64 + C A n f 2 557 1152 + 11 48 ζ 2 7 36 ζ 3 + C F n f 2 47 192 + ζ 3 6 + n f 3 25 1944 ζ 2 72 .
For the hard and soft functions we use the expansions H a b 12 = n = 0 ( α s d + n / π n ) H a b 12 ( n ) , where d is the power of α s in the leading-order cross section and it depends on the process, and S a b 12 = n = 0 ( α s / π ) n S a b 12 ( n ) . At lowest order, the soft matrix is given in terms of the color tensor basis by the expression S L I a b 12 ( 0 ) = tr c L * a b 12 c I a b 12 , while the hard matrix is real and symmetric, H L I a b 12 ( 0 ) = h L a b 12 ( 0 ) h I a b 12 ( 0 ) * with h L a b 12 ( 0 ) = ( S L K a b 12 ( 0 ) ) 1 tr ( c K a b 12 * M a b 12 ( 0 ) ) and h I a b 12 ( 0 ) * = tr ( M a b 12 ( 0 ) c K a b 12 ) ( S K I a b 12 ( 0 ) ) 1 where M a b 12 ( 0 ) is the lowest-order amplitude.
For the soft anomalous dimension we use the expansion Γ S a b 12 = n = 1 ( α s / π ) n Γ S a b 12 ( n ) . In the past, through the year 2008, most expressions for soft anomalous dimensions in the literature were given at one loop in axial gauge, while since 2009 most were given in Feynman gauge. The relation between the two is Γ S a b 12 ( 1 ) axial = Γ S a b 12 ( 1 ) Feyn . + ( 1 / 2 ) [ A a ( 1 ) + A b ( 1 ) + j = f . s . q , g A j ( 1 ) ] . Of course, the overall result for the resummed cross section is gauge-independent, and this is easily seen at one-loop from the fact that in axial gauge D i ( 1 ) axial = A i ( 1 ) , which precisely cancels the extra terms from change of gauge in the soft anomalous dimension.
We also note that in general we can disregard imaginary i π terms in the soft anomalous dimensions since such terms generally do not contribute. Of course, in cases with simple color flow when we only have functions—not matrices—it is clear from Equation (20) that any imaginary terms in Γ S a b 12 cancel out against imaginary terms in its Hermitian adjoint in the resummed cross section. Such terms also routinely cancel and do not contribute in fixed-order expansions even when we have processes that require matrices. Therefore, in the results in Section 5 through Section 10 we will drop any i π terms for simplicity.
In addition, we note that one can perform resummation in other schemes such as the DIS scheme (see, e.g., references [3,5,35] for the relation between the resummed expressions for the cross section in MS ¯ and DIS schemes). The expressions for the soft anomalous dimensions, however, remain the same.
When all external eikonal lines in the scattering process are massless, then the two-loop soft anomalous dimension matrix is proportional to the one-loop result [9]. Furthermore, three-parton correlations with massless eikonal lines do not contribute to the soft anomalous dimension at any order due to scaling symmetry constraints [36,37]. However, at three loops for the soft anomalous dimension in massless multi-leg scattering there are contributions from four-parton correlations [38].
When there are massive eikonal lines then the soft anomalous dimension at two or more loops is no longer proportional to the one-loop result [23]. Also, when two of the eikonal lines in the scattering process are massive, then the three-parton correlations no longer vanish; however, three-parton correlations still vanish when only one of the eikonal lines is massive.

3. NLO, NNLO, and N 3 LO Expansions of the Resummed Cross Section

In this section we expand the resummed cross section to NLO, NNLO, and N 3 LO, using Equations (20), (21), and (29) (see also [35,39]). Our expansions can be used for a variety of processes, but with the restrictions noted for the use of Equation (29) in the previous section.
The corrections take the form of plus distributions
D k ( s 4 ) = ln k ( s 4 / s ) s 4 + ,
defined by
0 s 4 max d s 4 f ( s 4 ) ln k ( s 4 / s ) s 4 + = 0 s 4 max d s 4 ln k ( s 4 / s ) s 4 f ( s 4 ) f ( 0 ) + 1 k + 1 ln k + 1 s 4 max s f ( 0 ) ,
where f is a smooth function, such as a pdf. Of course the above distributions, involving logarithms of s 4 / s , can readily be reexpressed in terms of logarithms of s 4 / M 2 for any hard scale M relevant to the process considered [35,39].
We note that the expressions can be simplified [35] for the cases of simple color structure where the soft anomalous dimensions are not matrices, as we will discuss in more detail below. We also note that we can extend the formulas for the fixed-order expansions to the general case of multi-particle final states, i.e., 2 n processes [11], by replacing m 2 2 by ( p 2 + + p n ) 2 in the formulas.

3.1. NLO Soft-Gluon Corrections

The NLO soft-gluon corrections are
d σ ^ a b 12 ( 1 ) = F a b 12 L O α s ( μ R ) π c 3 D 1 ( s 4 ) + c 2 D 0 ( s 4 ) + c 1 δ ( s 4 ) + α s d + 1 ( μ R ) π A a b 12 D 0 ( s 4 ) + T a b 12 δ ( s 4 ) ,
where F a b 12 L O = α s d tr H a b 12 ( 0 ) S a b 12 ( 0 ) denotes the leading-order (LO) coefficient,
c 3 = 2 ( A a ( 1 ) + A b ( 1 ) ) j = f . s . q , g A j ( 1 ) ,
and c 2 is given by c 2 = c 2 μ + T 2 , with
c 2 μ = ( A a ( 1 ) + A b ( 1 ) ) ln μ F 2 s
denoting the terms involving logarithms of the scale, and
T 2 = 2 A a ( 1 ) ln m 2 2 u s 2 A b ( 1 ) ln m 2 2 t s + D a ( 1 ) + D b ( 1 ) + j = f . s . q , g B j ( 1 ) + D j ( 1 )
denoting the scale-independent terms. Additionally,
A a b 12 = tr H a b 12 ( 0 ) Γ S a b 12 ( 1 ) S a b 12 ( 0 ) + H a b 12 ( 0 ) S a b 12 ( 0 ) Γ S a b 12 ( 1 ) .
For the cases with simple color structure where the soft anomalous dimension is not a matrix, we have α s d A a b 12 = 2 Re Γ S a b 12 ( 1 ) F a b 12 LO , and thus this term can be added to the term c 2 , thus simplifying the expressions at NLO and higher orders [35].
The δ ( s 4 ) terms involve a term c 1 that is proportional to the Born cross section, and a term T a b 12 that, in general, is not. We write c 1 = c 1 μ + T 1 , with
c 1 μ = A a ( 1 ) ln m 2 2 u s + A b ( 1 ) ln m 2 2 t s γ a ( 1 ) γ b ( 1 ) ln μ F 2 s + d β 0 4 ln μ R 2 s
denoting the terms involving logarithms of the scale. Additionally,
T 1 = A a ( 1 ) ln 2 m 2 2 u s + A b ( 1 ) ln 2 m 2 2 t s D a ( 1 ) ln m 2 2 u s D b ( 1 ) ln m 2 2 t s
and
T a b 12 = tr H a b 12 ( 0 ) S a b 12 ( 1 ) + H a b 12 ( 1 ) S a b 12 ( 0 ) .
We note that T a b 12 can also be determined via a comparison to a complete NLO calculation.

3.2. NNLO Soft-Gluon Corrections

The NNLO soft-gluon corrections are
d σ ^ a b 12 ( 2 ) = F a b 12 L O α s 2 ( μ R ) π 2 1 2 c 3 2 D 3 ( s 4 ) + 3 2 c 3 c 2 β 0 4 c 3 + β 0 8 j = f . s . q , g A j ( 1 ) D 2 ( s 4 ) + c 3 c 1 + c 2 2 ζ 2 c 3 2 β 0 2 T 2 + β 0 4 c 3 ln μ R 2 s + 2 ( A a ( 2 ) + A b ( 2 ) ) + j = f . s . q , g A j ( 2 ) + β 0 4 B j ( 1 ) D 1 ( s 4 ) + c 2 c 1 ζ 2 c 3 c 2 + ζ 3 c 3 2 β 0 2 T 1 + β 0 4 c 2 ln μ R 2 s 2 A a ( 2 ) ln m 2 2 u s 2 A b ( 2 ) ln m 2 2 t s + D a ( 2 ) + D b ( 2 ) + β 0 8 ( A a ( 1 ) + A b ( 1 ) ) ln 2 μ F 2 s ( A a ( 2 ) + A b ( 2 ) ) ln μ F 2 s + j = f . s . q , g B j ( 2 ) + D j ( 2 ) D 0 ( s 4 ) + α s d + 2 ( μ R ) π 2 3 2 c 3 A a b 12 D 2 ( s 4 ) + 2 c 2 β 0 2 A a b 12 + c 3 T a b 12 + F a b 12 D 1 ( s 4 ) + c 1 ζ 2 c 3 + β 0 4 ln μ R 2 s A a b 12 + c 2 T a b 12 + G a b 12 D 0 ( s 4 ) ,
where
F a b 12 = tr H a b 12 ( 0 ) Γ S a b 12 ( 1 ) 2 S a b 12 ( 0 ) + H a b 12 ( 0 ) S a b 12 ( 0 ) Γ S a b 12 ( 1 ) 2 + 2 H a b 12 ( 0 ) Γ S a b 12 ( 1 ) S a b 12 ( 0 ) Γ S a b 12 ( 1 )
and
G a b 12 = tr H a b 12 ( 1 ) Γ S a b 12 ( 1 ) S a b 12 ( 0 ) + H a b 12 ( 1 ) S a b 12 ( 0 ) Γ S a b 12 ( 1 ) + H a b 12 ( 0 ) Γ S a b 12 ( 1 ) S a b 12 ( 1 ) + H a b 12 ( 0 ) S a b 12 ( 1 ) Γ S a b 12 ( 1 ) + H a b 12 ( 0 ) Γ S a b 12 ( 2 ) S a b 12 ( 0 ) + H a b 12 ( 0 ) S a b 12 ( 0 ) Γ S a b 12 ( 2 ) .

3.3. N 3 LO Soft-Gluon Corrections

The N 3 LO soft-gluon corrections are
d σ ^ a b 12 ( 3 ) = F a b 12 L O α s 3 ( μ R ) π 3 1 8 c 3 3 D 5 ( s 4 ) + 5 8 c 3 2 c 2 5 24 c 3 2 β 0 + 5 48 c 3 β 0 j = f . s . q , g A j ( 1 ) D 4 ( s 4 ) + c 3 c 2 2 + 1 2 c 3 2 c 1 ζ 2 c 3 3 + β 0 2 12 c 3 5 6 β 0 c 3 c 2 + β 0 4 c 3 2 ln μ R 2 s β 0 2 c 3 A a ( 1 ) + A b ( 1 ) ln μ F 2 s + 2 c 3 A a ( 2 ) + A b ( 2 ) + j = f . s . q , g β 0 2 16 A j ( 1 ) c 3 A j ( 2 ) + β 0 6 c 2 A j ( 1 ) + c 3 β 0 4 B j ( 1 ) D 3 ( s 4 ) + 3 2 c 3 c 2 c 1 + 1 2 c 2 3 3 ζ 2 c 3 2 c 2 + 5 2 ζ 3 c 3 3 β 0 4 c 3 c 1 + 3 4 β 0 ζ 2 c 3 2 β 0 4 3 c 2 β 0 T 2 + β 0 8 c 3 6 c 2 β 0 ln μ R 2 s + ( 3 c 2 β 0 ) A a ( 2 ) + A b ( 2 ) A a ( 1 ) + A b ( 1 ) β 1 8 3 2 c 3 β 0 2 T 1 + 2 A a ( 2 ) ln m 2 2 u s + 2 A b ( 2 ) ln m 2 2 t s D a ( 2 ) D b ( 2 ) β 0 8 A a ( 1 ) + A b ( 1 ) ln 2 μ F 2 s + A a ( 2 ) + A b ( 2 ) ln μ F 2 s 3 2 c 3 j = f . s . q , g B j ( 2 ) D j ( 2 ) + β 0 4 ζ 2 A j ( 1 ) + 3 c 2 j = f . s . q , g A j ( 2 ) 2 + β 0 8 B j ( 1 ) + j = f . s . q , g A j ( 1 ) β 0 2 16 ln μ R 2 s + 3 β 0 4 A j ( 2 ) 3 β 0 2 16 B j ( 1 ) + 3 32 A j ( 1 ) β 1 + β 0 8 c 1 A j ( 1 ) D 2 ( s 4 ) + O D 1 ( s 4 ) + α s d + 3 ( μ R ) π 3 5 8 c 3 2 A a b 12 D 4 ( s 4 ) + 1 2 c 3 2 T a b 12 + 2 c 3 c 2 5 6 β 0 c 3 + β 0 6 j = f . s . q , g A j ( 1 ) A a b 12 + c 3 F a b 12 D 3 ( s 4 ) + 3 2 c 3 c 2 β 0 4 c 3 + β 0 8 j = f . s . q , g A j ( 1 ) T a b 12 + 3 2 c 2 F a b 12 + 3 2 c 3 G a b 12 + 1 2 K a b 12 + 3 2 c 3 c 1 + 3 2 c 2 2 3 ζ 2 c 3 2 3 4 β 0 ( c 2 + T 2 ) + 3 4 β 0 c 3 ln μ R 2 s + β 0 2 4 + 3 A a ( 2 ) + 3 A b ( 2 ) + 3 2 j = f . s . q , g ( A j ( 2 ) + β 0 4 B j ( 1 ) ) A a b 12 D 2 ( s 4 ) + O D 1 ( s 4 )
where
K a b 12 = tr H a b 12 ( 0 ) ( Γ S a b 12 ( 1 ) ) 3 S a b 12 ( 0 ) + H a b 12 ( 0 ) S a b 12 ( 0 ) ( Γ S a b 12 ( 1 ) ) 3 + 3 H a b 12 ( 0 ) ( Γ S a b 12 ( 1 ) ) 2 S a b 12 ( 0 ) Γ S a b 12 ( 1 ) + 3 H a b 12 ( 0 ) Γ S a b 12 ( 1 ) S a b 12 ( 0 ) ( Γ S a b 12 ( 1 ) ) 2 ,
and where we have omitted terms of order D 1 and D 0 (see also reference [39]).

4. Cusp Anomalous Dimension

The cusp anomalous dimension [23,40,41,42,43,44,45], involving two eikonal lines, is the simplest soft anomalous dimension, and it is an essential ingredient for all calculations of soft anomalous dimensions.
The cusp angle θ was introduced in Section 2.2, and it is given by
θ = ln p i · p j + ( p i · p j ) 2 p i 2 p j 2 p i 2 p j 2 = ln v i · v j + ( v i · v j ) 2 v i 2 v j 2 v i 2 v j 2
or equivalently θ = cosh 1 ( p i · p j / p i 2 p j 2 ) = cosh 1 ( v i · v j / v i 2 v j 2 ) , where p i μ = ( s / 2 ) v i μ and p j μ = ( s / 2 ) v j μ , where v i μ and v j μ are four-velocities.
The perturbative series for the cusp anomalous dimension in QCD is written as
Γ cusp = n = 1 α s π n Γ cusp ( n )
where α s is the strong coupling.
The cusp anomalous dimension at one loop [40] is given by
Γ cusp ( 1 ) = C F ( θ coth θ 1 ) .
The cusp anomalous dimension at two loops [23,43] is given by
Γ cusp ( 2 ) = K 2 Γ cusp ( 1 ) + C ( 2 )
where C ( 2 ) = C F C A C ( 2 ) , with [23]
C ( 2 ) = 1 2 + ζ 2 2 + θ 2 2 1 2 coth θ ζ 2 θ + θ 2 + θ 3 3 + Li 2 1 e 2 θ + 1 2 coth 2 θ ζ 3 + ζ 2 θ + θ 3 3 + θ Li 2 e 2 θ + Li 3 e 2 θ ,
and where
K 2 = C A 67 36 ζ 2 2 5 18 n f ;
i.e., K 2 = K / 2 = A i ( 2 ) / C i .
The cusp anomalous dimension at three loops [44,45] is given by
Γ cusp ( 3 ) = K 3 Γ cusp ( 1 ) + 2 K 2 Γ cusp ( 2 ) K 2 Γ cusp ( 1 ) + C ( 3 ) = K 3 Γ cusp ( 1 ) + 2 K 2 C ( 2 ) + C ( 3 ) ,
where
K 3 = C A 2 245 96 67 36 ζ 2 + 11 24 ζ 3 + 11 8 ζ 4 + C F n f 55 96 + ζ 3 2 + C A n f 209 432 + 5 ζ 2 18 7 ζ 3 12 n f 2 108 ,
i.e., K 3 = A i ( 3 ) / C i , and C ( 3 ) has a long expression which can be found in Equation (2.13) of reference [45].

4.1. Case with Two Massive Eikonal Lines

We now give explicit expressions for the cusp anomalous dimension at one, two, and three loops, when both eikonal lines represent massive quarks that have the same mass m (at the end of this subsection we will also consider two different masses). This is clearly applicable to the production of a heavy quark-antiquark pair. Here v i · v j = 1 + β 2 and v i 2 = v j 2 = 1 β 2 , where β = 1 4 m 2 / s is the speed of the massive quarks, and the cusp angle has the explicit form θ = ln [ ( 1 + β ) / ( 1 β ) ] ; also, we have β = tanh ( θ / 2 ) .
We note that coth θ = ( 1 + β 2 ) / ( 2 β ) , and we define
L β = ( 1 + β 2 ) 2 β ln 1 β 1 + β .
The cusp anomalous dimension at one loop is given by
Γ cusp β ( 1 ) = C F L β + 1
where we have added a superscript β to indicate that the cusp anomalous dimension is here given in terms of β , with both eikonal lines having mass m.
The cusp anomalous dimension at two loops is given by
Γ cusp β ( 2 ) = K 2 Γ cusp β ( 1 ) + C β ( 2 )
where C β ( 2 ) = C F C A C β ( 2 ) , with [23]
C β ( 2 ) = 1 2 + ζ 2 2 + 1 2 ln 2 1 β 1 + β + ( 1 + β 2 ) 4 β ζ 2 ln 1 β 1 + β ln 2 1 β 1 + β + 1 3 ln 3 1 β 1 + β Li 2 4 β ( 1 + β ) 2 + ( 1 + β 2 ) 2 8 β 2 ζ 3 ζ 2 ln 1 β 1 + β 1 3 ln 3 1 β 1 + β ln 1 β 1 + β Li 2 ( 1 β ) 2 ( 1 + β ) 2 + Li 3 ( 1 β ) 2 ( 1 + β ) 2 .
The cusp anomalous dimension at three loops is given by
Γ cusp β ( 3 ) = K 3 Γ cusp β ( 1 ) + 2 K 2 Γ cusp β ( 2 ) K 2 Γ cusp β ( 1 ) + C β ( 3 ) = K 3 Γ cusp β ( 1 ) + 2 K 2 C β ( 2 ) + C β ( 3 ) ,
where C β ( 3 ) = C F C A 2 C β ( 3 ) with
C β ( 3 ) = ζ 2 2 ζ 3 8 9 8 ζ 4 + ζ 2 2 ln 1 β 1 + β 1 4 ln 2 1 β 1 + β + 1 12 ln 3 1 β 1 + β 1 24 ln 4 1 β 1 + β + 1 4 ln 2 1 β 1 + β ln 4 β ( 1 + β ) 2 + 3 4 ln 1 β 1 + β Li 2 ( 1 β ) 2 ( 1 + β ) 2 5 8 Li 3 ( 1 β ) 2 ( 1 + β ) 2 + ( 1 + β 2 ) 2 β ζ 3 4 + 15 8 ζ 4 ζ 2 2 ζ 3 2 + 9 8 ζ 4 ln 1 β 1 + β + 1 4 + ζ 2 ln 2 1 β 1 + β 1 12 + ζ 2 3 ln 3 1 β 1 + β + 7 24 ln 4 1 β 1 + β 1 24 ln 5 1 β 1 + β + 1 2 ln 2 1 β 1 + β ln 4 β ( 1 + β ) 2 1 2 ln 3 1 β 1 + β ln 4 β ( 1 + β ) 2 3 4 ln 2 1 β 1 + β Li 2 ( 1 β ) 2 ( 1 + β ) 2 + 1 4 Li 2 4 β ( 1 + β ) 2 + 1 4 Li 3 ( 1 β ) 2 ( 1 + β ) 2 + 7 4 ln 1 β 1 + β Li 3 ( 1 β ) 2 ( 1 + β ) 2 + 1 2 Li 3 4 β ( 1 + β ) 2 15 8 Li 4 ( 1 β ) 2 ( 1 + β ) 2 + ( 1 + β 2 ) 2 4 β 2 ζ 2 ζ 3 2 19 8 ζ 4 + 3 2 ζ 5 3 2 ζ 3 15 8 ζ 4 ln 1 β 1 + β + ζ 3 4 ζ 2 ln 2 1 β 1 + β + 2 3 ζ 2 ln 3 1 β 1 + β 1 4 ln 4 1 β 1 + β + 11 120 ln 5 1 β 1 + β + ln 4 β ( 1 + β ) 2 ζ 3 + ζ 2 ln 1 β 1 + β ζ 2 ln 2 1 β 1 + β + 1 3 ln 3 1 β 1 + β 1 6 ln 4 1 β 1 + β ln 2 1 β 1 + β ln 2 4 β ( 1 + β ) 2 + ln 1 β 1 + β ln 3 4 β ( 1 + β ) 2 1 8 ln 4 4 β ( 1 + β ) 2 + ζ 2 2 ζ 2 ln 1 β 1 + β 2 ln 2 1 β 1 + β 1 12 ln 3 1 β 1 + β + ln 1 β 1 + β ln 4 β ( 1 + β ) 2 Li 2 ( 1 β ) 2 ( 1 + β ) 2 1 4 Li 2 2 ( 1 β ) 2 ( 1 + β ) 2 + 1 2 ln 2 4 β ( 1 + β ) 2 Li 2 4 β ( 1 + β ) 2 + 1 4 Li 2 2 4 β ( 1 + β ) 2 1 2 ln 2 4 β ( 1 β ) 2 Li 2 ( 1 β ) 2 4 β + ζ 2 2 + 3 2 ln 1 β 1 + β 1 4 ln 2 1 β 1 + β ln 4 β ( 1 + β ) 2 Li 3 ( 1 β ) 2 ( 1 + β ) 2 + ln 1 β 1 + β ln 4 β ( 1 + β ) 2 Li 3 4 β ( 1 + β ) 2 + 2 ln 1 β 1 + β ln 4 β ( 1 + β ) 2 Li 3 ( 1 β ) 2 4 β + 9 8 ln 1 β 1 + β Li 4 ( 1 β ) 2 ( 1 + β ) 2 + Li 4 4 β ( 1 + β ) 2 Li 4 ( 1 β ) 2 4 β 3 2 Li 5 ( 1 β ) 2 ( 1 + β ) 2 + 1 4 A ( β ) A ( 0 ) + B ( β ) B ( 0 )
with
A ( β ) = ( 1 + β 2 ) 3 8 β 3 3 ζ 5 4 ζ 4 ln 1 β 1 + β 3 ζ 3 ln 2 1 β 1 + β 4 3 ζ 2 ln 3 1 β 1 + β 1 5 ln 5 1 β 1 + β 2 3 ln 3 1 β 1 + β Li 2 ( 1 β ) 2 ( 1 + β ) 2 + ln 2 1 β 1 + β Li 3 ( 1 β ) 2 ( 1 + β ) 2 2 ln 1 β 1 + β Li 4 ( 1 β ) 2 ( 1 + β ) 2 + 3 Li 5 ( 1 β ) 2 ( 1 + β ) 2 + H 1 , 1 , 0 , 0 , 1 4 β ( 1 + β ) 2 + H 1 , 0 , 1 , 0 , 1 4 β ( 1 + β ) 2 ,
where A ( 0 ) = 2 and explicit expressions for the harmonic polylogarithms [46] H 1 , 1 , 0 , 0 , 1 and H 1 , 0 , 1 , 0 , 1 can be found in the Appendix of reference [45], and
B ( β ) = ( 1 β 2 ) 4 β 2 ζ 2 ζ 3 + 2 ζ 3 ln 2 1 β 1 + β + 3 2 ζ 4 1 6 ln 4 1 β 1 + β ln 2 β 1 β + 3 2 ζ 4 2 ζ 3 ln 1 β 1 + β 1 6 ln 4 1 β 1 + β ln 1 β 2 + 2 ζ 3 Li 2 1 + β 1 + β + Li 2 2 β 1 + β 2 3 ln 3 1 β 1 + β Li 2 1 β 1 + β Li 2 1 + β 1 + β + 2 ln 2 1 β 1 + β Li 3 1 β 1 + β Li 3 1 + β 1 + β 4 ln 1 β 1 + β Li 4 1 β 1 + β Li 4 1 + β 1 + β + 4 Li 5 1 β 1 + β 4 Li 5 1 + β 1 + β + 4 H 1 , 0 , 1 , 0 , 0 1 β 1 + β + H 1 , 0 , 1 , 0 , 0 1 β 1 + β H 1 , 0 , 1 , 0 , 0 1 β 1 + β H 1 , 0 , 1 , 0 , 0 1 β 1 + β ,
where B ( 0 ) = 3 ζ 3 / 2 and explicit expressions for the harmonic polylogarithms, H 1 , 0 , 1 , 0 , 0 , H 1 , 0 , 1 , 0 , 0 , H 1 , 0 , 1 , 0 , 0 , and H 1 , 0 , 1 , 0 , 0 can be found in the Appendix of reference [45].
The limit of the n-loop cusp anomalous dimension as β goes to 1, or equivalently as m 2 / s 0 , can be written as
Γ cusp β ( n ) β 1 K n C F ln m 2 s + R n
where we define K n = A i ( n ) / C i (and thus K 1 = 1 ), and the terms R n at one, two, and three loops are given, respectively, by R 1 = K 1 C F , R 2 = K 2 C F + ( 1 / 2 ) C F C A ( 1 ζ 3 ) , and
R 3 = K 3 C F + K 2 C F C A ( 1 ζ 3 ) + C F C A 2 1 2 + 3 4 ζ 2 ζ 3 4 3 4 ζ 2 ζ 3 + 9 8 ζ 5 .
Finally, we consider the case when the two lines have different masses, m 1 and m 2 . In that case, we define
L β 1 β 2 = ( 1 + β 1 β 2 ) 2 ( β 1 + β 2 ) ln ( 1 β 1 ) ( 1 β 2 ) ( 1 + β 1 ) ( 1 + β 2 ) ,
where β 1 = 1 4 m 1 2 s ( s + m 1 2 m 2 2 ) 2 1 / 2 and β 2 = 1 4 m 2 2 s ( s + m 2 2 m 1 2 ) 2 1 / 2 , and the cusp anomalous dimension at one loop is then given by
Γ cusp β 1 β 2 ( 1 ) = C F L β 1 β 2 + 1 .

4.2. Case with One Massless and One Massive Eikonal Line

We take the eikonal line i to represent a massive quark of mass m i , and the eikonal line j to represent a massless quark. To find the expressions for this case, we take the limit of the massive expression as the mass of eikonal line j goes to zero, taking into account the self energies. We note that in this limit θ = ln [ 2 p i · p j / ( m i p j 2 ) ] . At one loop, the heavy quark self-energy is - C F / 2 , so removing it adds C F / 2 , and then we add the massless contribution C F ln p j 2 / s . The overall change in the self-energy contributions relative to the fully massive case is thus C F / 2 + C F ln p j 2 / s , which equals R 1 / 2 + C F ln p j 2 / s , where R 1 was defined in the previous subsection. Thus, we find at one loop
Γ cusp m i ( 1 ) = C F ln 2 p i · p j m i s 1 2 ,
where we have added a superscript m i to indicate that only eikonal line i has a mass.
At two loops, the heavy-quark self-energy is ( K 2 / 2 ) C F + C F C A ( 1 ζ 2 ) / 4 , so again we remove it and then add the contribution for the massless eikonal line, which is C F K 2 ln p j 2 / s C F C A ( ζ 2 ζ 3 ) / 4 . The total additional contribution is thus ( K 2 / 2 ) C F C F C A ( 1 ζ 3 ) / 4 + C F K 2 ln p j 2 / s , which equals R 2 / 2 + C F K 2 ln p j 2 / s , where R 2 was defined in the previous subsection. Thus, we find the two-loop result
Γ cusp m i ( 2 ) = K 2 Γ cusp m i ( 1 ) + 1 4 C F C A ( 1 ζ 3 ) .
We note that this can be rewritten as Γ cusp m i ( 2 ) = K 2 C F ln ( 2 p i · p j / ( m i s ) ) + R 2 / 2 .
At three loops, again removing the heavy-quark self-energy and adding the massless contribution, we find a total additional contribution of R 3 / 2 + C F K 3 ln p j 2 / s , where R 3 was defined in the previous subsection in Equation (66). Thus, we find the three-loop result
Γ cusp m i ( 3 ) = K 3 Γ cusp m i ( 1 ) + 1 2 K 2 C F C A ( 1 ζ 3 ) + C F C A 2 1 4 + 3 8 ζ 2 ζ 3 8 3 8 ζ 2 ζ 3 + 9 16 ζ 5 .
We note that this can be rewritten as Γ cusp m i ( 3 ) = K 3 C F ln ( 2 p i · p j / ( m i s ) ) + R 3 / 2 .
We also note that in general the n-loop result can be written as Γ cusp m i ( n ) = K n C F ln ( 2 p i · p j / ( m i s ) ) + R n / 2 .

4.3. Case with Two Massless Eikonal Lines

In the case where both eikonal lines are massless, θ = ln ( 2 p i · p j / p i 2 p j 2 ) , and again removing the heavy-quark self-energies and adding the massless contributions, we find the lightlike cusp anomalous dimension as
Γ cusp massless = C F ln 2 p i · p j s n = 1 α s π n K n .

5. Γ S for Some Processes with Trivial Color Structure

The soft anomalous dimension vanishes to all orders for many processes with trivial color structure, with no colored particles in the final state. Some well-known processes of this type are the Drell–Yan processes q q ¯ γ * , q q ¯ Z ; W-boson production via q q ¯ W ± ; Higgs production via b b ¯ H and g g H ; production of electroweak boson pairs q q ¯ γ γ , q q ¯ Z Z , q q ¯ W + W , q q ¯ γ Z ; q q ¯ W ± γ ; q q ¯ W ± Z ; and charged Higgs production via b b ¯ H W + , b b ¯ H + H , g g H + H .
The soft anomalous dimension for deep-inelastic-scattering (DIS) l q l q with underlying process q γ * q is given at one, two, and three loops by
Γ S q γ * q ( 1 ) = C F ln ( t / s ) , Γ S q γ * q ( 2 ) = K 2 C F ln ( t / s ) , Γ S q γ * q ( 3 ) = K 3 C F ln ( t / s ) ,
respectively.

6. Γ S for Large- p T W Production and Related Processes

We next consider processes with a W-boson or a Z-boson or a Higgs boson produced at large- p T , as well as direct-photon production and related processes. For these processes there is only a single color tensor, coupling the two quarks (or quark and antiquark) to the gluon in an octet state. Thus the soft anomalous dimension is a simple function, not a matrix.
The soft anomalous dimensions for the processes q g W ± q and q g Z q and q g γ q and b g H b , are all identical. Using W production at large p T as the specific process, the soft anomalous dimension is given at one loop by [8,35,47]
Γ S q g W q ( 1 ) = C F ln u s + C A 2 ln t u ,
at two-loops by [47]
Γ S q g W q ( 2 ) = K 2 Γ S q g W q ( 1 ) ,
and at three loops by
Γ S q g W q ( 3 ) = K 3 Γ S q g W q ( 1 ) .
For q q ¯ W ± g or q q ¯ Z g or q q ¯ γ g or b b ¯ H g , the corresponding results are
Γ S q q ¯ W g ( 1 ) = C A 2 ln t u s 2 , Γ S q q ¯ W g ( 2 ) = K 2 Γ S q q ¯ W g ( 1 ) , Γ S q q ¯ W g ( 3 ) = K 3 Γ S q q ¯ W g ( 1 ) .
We also note that the soft anomalous dimensions for the reverse processes γ q q g and γ g q q ¯ are the same as for the corresponding processes above.

7. Γ S for Single-Top Production and Related Processes

We continue with various single-top production processes. They include s-channel, t-channel, and t W production, and various related FCNC single-top processes.

7.1. s-Channel Single-Top Production

Soft anomalous dimensions for s-channel single-top production were calculated at one loop in references [12,48,49], at two loops in [49,50], and at three loops in [50].
The partonic processes are q q ¯ t b ¯ . In this channel we have 2 2 processes at lowest order that involve a final-state top quark and a final-state massless quark. Thus, we have four colored particles involved in the scattering, one of which is massive.
The color structure of the hard scattering in s-channel single-top production is more complicated and thus the soft anomalous dimension is a 2 × 2 matrix in color space. We choose a singlet-octet s-channel color basis, c 1 q q ¯ t b ¯ = δ a b δ 12 and c 2 q q ¯ t b ¯ = T b a c T 12 c .
The four matrix elements of the s-channel soft anomalous dimension matrix, Γ S q q ¯ t b ¯ , are given at one loop by
Γ 11 q q ¯ t b ¯ ( 1 ) = C F ln s m t 2 m t s 1 2 , Γ 12 q q ¯ t b ¯ ( 1 ) = C F 2 N c ln t ( t m t 2 ) u ( u m t 2 ) , Γ 21 q q ¯ t b ¯ ( 1 ) = ln t ( t m t 2 ) u ( u m t 2 ) , Γ 22 q q ¯ t b ¯ ( 1 ) = C F C A 2 ln s m t 2 m t s 1 2 + 2 ln t ( t m t 2 ) u ( u m t 2 ) + C A 2 ln t ( t m t 2 ) m t s 3 / 2 1 2
where m t is the top-quark mass.
At two loops, we have
Γ 11 q q ¯ t b ¯ ( 2 ) = K 2 Γ 11 q q ¯ t b ¯ ( 1 ) + 1 4 C F C A ( 1 ζ 3 ) , Γ 12 q q ¯ t b ¯ ( 2 ) = K 2 Γ 12 q q ¯ t b ¯ ( 1 ) , Γ 21 q q ¯ t b ¯ ( 2 ) = K 2 Γ 21 q q ¯ t b ¯ ( 1 ) , Γ 22 q q ¯ t b ¯ ( 2 ) = K 2 Γ 22 q q ¯ t b ¯ ( 1 ) + 1 4 C F C A ( 1 ζ 3 ) .
At three loops, only the first element of the soft anomalous dimension matrix is needed to calculate the N 3 LO soft-gluon corrections, because only the first element of the leading-order hard matrix is nonzero. We have
Γ 11 q q ¯ t b ¯ ( 3 ) = K 3 Γ 11 q q ¯ t b ¯ ( 1 ) + 1 2 K 2 C F C A ( 1 ζ 3 ) + C F C A 2 1 4 + 3 8 ζ 2 ζ 3 8 3 8 ζ 2 ζ 3 + 9 16 ζ 5 .
Furthermore, up to unknown contributions from four-parton correlations, Γ 22 q q ¯ t b ¯ ( 3 ) should have the same form as Equation (80) (just replace the 11 subscripts by 22), and the two off-diagonal elements at three loops should have a similar form to Equation (79) (replace K 2 by K 3 ).

7.2. t-Channel Single-Top Production

Soft anomalous dimensions for t-channel single-top production were calculated at one loop in references [12,48,51], at two loops in [50,51], and at three loops in [50].
The partonic processes are b q t q . The color structure of the hard scattering in t-channel single-top production is again complicated, and the soft anomalous dimension is a 2 × 2 matrix in color space. We choose a singlet-octet t-channel color basis, c 1 b q t q = δ a 1 δ b 2 and c 2 b q t q = T 1 a c T 2 b c .
The four matrix elements of the t-channel soft anomalous dimension matrix Γ S b q t q are given at one loop by
Γ 11 b q t q ( 1 ) = C F ln t ( t m t 2 ) m t s 3 / 2 1 2 , Γ 12 b q t q ( 1 ) = C F 2 N c ln u ( u m t 2 ) s ( s m t 2 ) , Γ 21 b q t q ( 1 ) = ln u ( u m t 2 ) s ( s m t 2 ) , Γ 22 b q t q ( 1 ) = C F C A 2 ln t ( t m t 2 ) m t s 3 / 2 1 2 + 2 ln u ( u m t 2 ) s ( s m t 2 ) + C A 2 ln u ( u m t 2 ) m t s 3 / 2 1 2 .
At two loops, our calculation gives
Γ 11 b q t q ( 2 ) = K 2 Γ 11 b q t q ( 1 ) + 1 4 C F C A ( 1 ζ 3 ) , Γ 12 b q t q ( 2 ) = K 2 Γ 12 b q t q ( 1 ) , Γ 21 b q t q ( 2 ) = K 2 Γ 21 b q t q ( 1 ) , Γ 22 b q t q ( 2 ) = K 2 Γ 22 b q t q ( 1 ) + 1 4 C F C A ( 1 ζ 3 ) .
At three loops, only the first element of the soft anomalous dimension matrix is needed to calculate the N 3 LO soft-gluon corrections, as for the s-channel. We find
Γ 11 b q t q ( 3 ) = K 3 Γ 11 b q t q ( 1 ) + 1 2 K 2 C F C A ( 1 ζ 3 ) + C F C A 2 1 4 + 3 8 ζ 2 ζ 3 8 3 8 ζ 2 ζ 3 + 9 16 ζ 5 .
Furthermore, up to unknown contributions from four-parton correlations, Γ 22 b q t q ( 3 ) should have the same form as Equation (83) (just replace the 11 subscripts by 22), and the two off-diagonal elements at three loops should be of similar form to Equation (82) (replace K 2 by K 3 ).

7.3. b g t W and Related Processes

We next present the soft anomalous dimension for the associated production of a top quark with a W-boson via b g t W which is known at one-loop [35,48], two loops [52], and three loops [50]. The soft anomalous dimension for t W production is identical to that for other related processes in models of new physics, such as the associated production of a top quark with a charged Higgs boson via b g t H in two-Higgs-doublet models, and flavor-changing-neutral-current (FCNC) processes that proceed via anomalous top-quark couplings, such q g t Z (or t Z ) and q g t γ (see [12] for a review). In all these cases we have 2 2 processes at lowest order that involve a final-state top quark and a final-state boson.
The soft anomalous dimension at one-loop for b g t W (and all other processes in this set) is given by
Γ S b g t W ( 1 ) = C F ln m t 2 t m t s 1 2 + C A 2 ln u m t 2 t m t 2 .
The two-loop result is given by
Γ S b g t W ( 2 ) = K 2 Γ S b g t W ( 1 ) + 1 4 C F C A ( 1 ζ 3 ) .
The three-loop result is
Γ S b g t W ( 3 ) = K 3 Γ S b g t W ( 1 ) + 1 2 K 2 C F C A ( 1 ζ 3 ) + C F C A 2 1 4 + 3 8 ζ 2 ζ 3 8 3 8 ζ 2 ζ 3 + 9 16 ζ 5 .

7.4. FCNC u e t e

For the FCNC process u e t e , which proceeds via anomalous t-q- γ and t-q-Z couplings, we have [35,53,54]
Γ S u e t e ( 1 ) = C F ln m t 2 t m t s 1 2 .
The two-loop result is given by
Γ S u e t e ( 2 ) = K 2 Γ S u e t e ( 1 ) + 1 4 C F C A ( 1 ζ 3 )
and the three-loop result by
Γ S u e t e ( 3 ) = K 3 Γ S u e t e ( 1 ) + 1 2 K 2 C F C A ( 1 ζ 3 ) + C F C A 2 1 4 + 3 8 ζ 2 ζ 3 8 3 8 ζ 2 ζ 3 + 9 16 ζ 5 .

7.5. FCNC t g Production

We next consider t g production via an anomalous t-q-g coupling [55]. For the partonic process g u t g we choose the color basis c 1 g u t g = δ b 1 δ a 2 , c 2 g u t g = d a 2 c T 1 b c , c 3 g u t g = i f a 2 c T 1 b c . Then the one-loop soft anomalous dimension is
Γ S g u t g ( 1 ) = Γ 11 g u t g ( 1 ) 0 Γ 13 g u t g ( 1 ) 0 Γ 22 g u t g ( 1 ) Γ 23 g u t g ( 1 ) Γ 31 g u t g ( 1 ) Γ 32 g u t g ( 1 ) Γ 22 g u t g ( 1 )
where [55]
Γ 11 g u t g ( 1 ) = C F ln m t 2 u m t s 1 2 + C A ln u s , Γ 31 g u t g ( 1 ) = ln t ( t m t 2 ) s ( s m t 2 ) , Γ 32 g u t g ( 1 ) = ( N c 2 4 ) 4 N c ln t ( t m t 2 ) s ( s m t 2 ) , Γ 22 g u t g ( 1 ) = C F ln m t 2 u m t s 1 2 + C A 4 ln t u 2 ( s m t 2 ) ( t m t 2 ) ( u m t 2 ) 2 s 3 , Γ 23 g u t g ( 1 ) = C A 4 ln t ( t m t 2 ) s ( s m t 2 ) , Γ 13 g u t g ( 1 ) = 1 2 ln t ( t m t 2 ) s ( s m t 2 ) .
For the two-loop soft anomalous dimension matrix, as for s-channel and t-channel single-top production in the previous subsections, the two-loop matrix elements are given by K 2 times the corresponding one-loop elements, with an additional term C F C A ( 1 ζ 3 ) / 4 in the diagonal elements.

8. Γ S for t t ¯ Production and Related Processes

In this section we discuss the soft anomalous dimension matrices for top-antitop pair production, which of course are the same for bottom quark or charm quark production, and related processes such as DIS heavy-quark production, FCNC t t production, and squark and gluino production.

8.1. t t ¯ Production in Hadronic Collisions

The top-antitop pair production partonic processes at lowest order are q q ¯ t t ¯ and g g t t ¯ . Next, we present the one-loop and two-loop results for the soft anomalous matrices for these partonic processes [3,5,23,56,57,58] and a form for the three-loop results.

8.1.1. q q ¯ t t ¯

The soft anomalous dimension matrix Γ S q q ¯ t t ¯ for the process q q ¯ t t ¯ is a 2 × 2 matrix. We use a color tensor basis of s-channel singlet and octet exchange, c 1 q q ¯ t t ¯ = δ a b δ 12 , c 2 q q ¯ t t ¯ = T b a c T 12 c .
The four matrix elements of Γ S q q ¯ t t ¯ are given at one loop [3,5,58] by
Γ 11 q q ¯ t t ¯ ( 1 ) = = C F L β + 1 , Γ 12 q q ¯ t t ¯ ( 1 ) = C F N c ln t m t 2 u m t 2 , Γ 21 q q ¯ t t ¯ ( 1 ) = 2 ln t m t 2 u m t 2 , Γ 22 q q ¯ t t ¯ ( 1 ) = C F C A 2 L β 1 + 4 ln t m t 2 u m t 2 + C A 2 ln ( t m t 2 ) 2 s m t 2 1 ,
where L β is given by Equation (57). We note that the first element of the matrix is equal to Γ cusp β ( 1 ) , Equation (58).
At two loops we have [23,56,57,58]
Γ 11 q q ¯ t t ¯ ( 2 ) = Γ cusp β ( 2 ) , Γ 12 q q ¯ t t ¯ ( 2 ) = K 2 C A N 2 β Γ 12 q q ¯ t t ¯ ( 1 ) , Γ 21 q q ¯ t t ¯ ( 2 ) = K 2 + C A N 2 β Γ 21 q q ¯ t t ¯ ( 1 ) , Γ 22 q q ¯ t t ¯ ( 2 ) = K 2 Γ 22 q q ¯ t t ¯ ( 1 ) + C A C F C A 2 C β ( 2 ) + C A 2 4 ( 1 ζ 3 ) ,
where Γ cusp β ( 2 ) is given by Equation (59), C β ( 2 ) is given by Equation (60), and N 2 β is given by
N 2 β = 1 4 ln 2 1 β 1 + β + ( 1 + β 2 ) 8 β ζ 2 ln 2 1 β 1 + β Li 2 4 β ( 1 + β ) 2 .
At three loops, we expect a similar structure up to four-parton correlations; i.e., the first three matrix elements should have the same form as in Equation (93) (replace two-loop quantities by three-loop ones), while
Γ 22 q q ¯ t t ¯ ( 3 ) = K 3 Γ 22 q q ¯ t t ¯ ( 1 ) + C A C F C A 2 C A C β ( 3 ) + 2 K 2 C β ( 2 ) + K 2 2 C A 2 ( 1 ζ 3 ) + C A 3 1 4 + 3 8 ζ 2 ζ 3 8 3 8 ζ 2 ζ 3 + 9 16 ζ 5 + X 22 q q ¯ t t ¯ ( 3 ) 4 p ,
where C β ( 3 ) is given by Equation (62), and X 22 q q ¯ t t ¯ ( 3 ) 4 p denotes the unknown three-loop contributions from four-parton correlations.

8.1.2. g g t t ¯

The soft anomalous dimension matrix Γ S g g t t ¯ for the process g g t t ¯ in a color tensor basis c 1 g g t t ¯ = δ a b δ 12 , c 2 g g t t ¯ = d a b c T 12 c , c 3 g g t t ¯ = i f a b c T 12 c , is given by
Γ S g g t t ¯ = Γ 11 g g t t ¯ 0 Γ 13 g g t t ¯ 0 Γ 22 g g t t ¯ Γ 23 g g t t ¯ Γ 31 g g t t ¯ Γ 32 g g t t ¯ Γ 22 g g t t ¯ .
At one loop we have [5,58]
Γ 11 g g t t ¯ ( 1 ) = C F L β + 1 , Γ 13 g g t t ¯ ( 1 ) = ln t m t 2 u m t 2 , Γ 31 g g t t ¯ ( 1 ) = 2 ln t m t 2 u m t 2 , Γ 22 g g t t ¯ ( 1 ) = C F C A 2 L β 1 + C A 2 ln ( t m t 2 ) ( u m t 2 ) s m t 2 1 , Γ 23 g g t t ¯ ( 1 ) = C A 2 ln t m t 2 u m t 2 , Γ 32 g g t t ¯ ( 1 ) = ( N c 2 4 ) 2 N c ln t m t 2 u m t 2 .
At two loops we find [23,57,58]
Γ 11 g g t t ¯ ( 2 ) = Γ cusp β ( 2 ) , Γ 13 g g t t ¯ ( 2 ) = K 2 C A N 2 β Γ 13 g g t t ¯ ( 1 ) , Γ 31 g g t t ¯ ( 2 ) = K 2 + C A N 2 β Γ 31 g g t t ¯ ( 1 ) , Γ 22 g g t t ¯ ( 2 ) = K 2 Γ 22 g g t t ¯ ( 1 ) + C A C F C A 2 C β ( 2 ) + C A 2 4 ( 1 ζ 3 ) , Γ 23 g g t t ¯ ( 2 ) = K 2 Γ 23 g g t t ¯ ( 1 ) , Γ 32 g g t t ¯ ( 2 ) = K 2 Γ 32 g g t t ¯ ( 1 ) .
At three loops, we expect a similar structure but with additional four-parton correlations, as discussed in the previous subsection; i.e., replace two-loop quantities by three-loop ones, and Γ 22 g g t t ¯ ( 3 ) of the same general form as Equation (95).

8.2. DIS Heavy-Quark Production

For heavy-quark production in DIS (also known as electroproduction of heavy quarks), e p e Q Q ¯ , the underlying process is g γ * Q Q ¯ . As for direct photon production there is only a single color tensor, coupling the produced pair to the gluon in an octet state. The soft anomalous dimension is given at one loop [35,59] by
Γ S g γ * Q Q ¯ ( 1 ) = C F ( L β + 1 ) + C A 2 L β + ln ( t m Q 2 ) ( u m Q 2 ) m Q 2 s ,
where L β is given by Equation (57).

8.3. e + e t t ¯

The soft anomalous dimension for e + e t t ¯ is simply the cusp anomalous dimension for the case of both eikonal lines of mass m t [23,45] that was presented in Section 4.1.

8.4. FCNC q q t t

For the process u u t t [54], which proceeds via anomalous t-q- γ and t-q-Z couplings, we choose a color basis consisting of singlet exchange in the t and u channels, c 1 = δ a 1 δ b 2 and c 2 = δ a 2 δ b 1 . Then Γ S q q t t is a 2 × 2 soft anomalous dimension matrix, and its matrix elements at one loop [54] are given by
Γ 11 q q t t ( 1 ) = C F 2 ln m t 2 t m t s 1 + C F C A 2 2 ln m t 2 u m t s + L β , Γ 12 q q t t ( 1 ) = ln m t 2 t m t s + 1 2 L β , Γ 21 q q t t ( 1 ) = ln m t 2 u m t s + 1 2 L β , Γ 22 q q t t ( 1 ) = C F 2 ln m t 2 u m t s 1 + C F C A 2 2 ln m t 2 t m t s + L β ,
where L β is given by Equation (57). We, of course, note that this process is similar to t t ¯ production via the q q ¯ channel, but a different choice of color basis here leads to a different form for the results.

8.5. Squark and Gluino Production

For squark production via the process q q ¯ q ˜ q ¯ ˜ , the soft anomalous dimension is of the same form as for the top-production process q q ¯ t t ¯ in Section 8.1.1 (just replace the top-quark mass by the squark mass [35]), with a similar result for the channel q q q ˜ q ˜ (see also [60,61]).
For squark production via the process g g q ˜ q ¯ ˜ , the soft anomalous dimension is of the same form as for the top-production process g g t t ¯ in Section 8.1.2 (again, just replace the top-quark mass by the squark mass [35]). A modified form of this matrix describes gluino production via the process q q ¯ g ˜ g ˜ , now using the gluino mass [60]. An analogous result describes squark and gluino production via the process q g q ˜ g ˜ [61].
For gluino production via the process g g g ˜ g ˜ , the color structure is more complicated (the same as for g g g g in Section 9.5), and the soft anomalous dimension matrix is given in [60].

9. Γ S for Jet Production and Related Processes

In this section we present the soft anomalous dimension matrices for partonic processes involved in jet production [7]; these soft anomalous dimensions are also relevant to related processes such as hadron production.

9.1. q q ¯ q q ¯

We begin with the quark-antiquark annihilation processes, q q ¯ q q ¯ . There are three different types of quark-antiquark processes here, depending on the quark flavors: q j q ¯ j q j q ¯ j , q j q ¯ j q k q ¯ k , and q j q ¯ k q j q ¯ k .
In the t-channel singlet-octet color basis c 1 q q ¯ q q ¯ = δ a 1 δ b 2 , c 2 q q ¯ q q ¯ = T 1 a c T b 2 c , the one-loop soft anomalous dimension matrix is [7]
Γ S q q ¯ q q ¯ ( 1 ) = 2 C F ln ( t / s ) C F N c ln ( u / s ) 2 ln ( u / s ) 1 N c ln ( t s / u 2 ) .
At two loops, Γ S q q ¯ q q ¯ ( 2 ) = K 2 Γ S q q ¯ q q ¯ ( 1 ) .

9.2. q q q q and q ¯ q ¯ q ¯ q ¯

Next, we discuss quark–quark scattering processes, q q q q . There are two different types of quark–quark processes here, depending on the quark flavors: q j q j q j q j and q j q k q j q k .
In the t-channel octet-singlet color basis c 1 q q q q = T 1 a c T 2 b c , c 2 q q q q = δ a 1 δ b 2 , the one-loop soft anomalous dimension matrix is [7]
Γ S q q q q ( 1 ) = 1 N c ln ( t u / s 2 ) + 2 C F ln ( u / s ) 2 ln ( u / s ) C F N c ln ( u / s ) 2 C F ln ( t / s ) .
At two loops, Γ S q q q q ( 2 ) = K 2 Γ S q q q q ( 1 ) .
The same soft anomalous dimension matrix describes the process with antiquarks, q ¯ q ¯ q ¯ q ¯ .

9.3. q q ¯ g g and g g q q ¯

Next, we discuss the processes q q ¯ g g and g g q q ¯ .
For the process q q ¯ g g , in the s-channel color basis c 1 q q ¯ g g = δ a b δ 12 , c 2 q q ¯ g g = d 12 c T b a c , c 3 q q ¯ g g = i f 12 c T b a c , the one-loop soft anomalous dimension matrix is [7]
Γ S q q ¯ g g ( 1 ) = 0 0 ln ( u / t ) 0 C A 2 ln ( t u / s 2 ) C A 2 ln ( u / t ) 2 ln ( u / t ) N c 2 4 2 N c ln ( u / t ) C A 2 ln ( t u / s 2 ) .
At two loops, Γ S q q ¯ g g ( 2 ) = K 2 Γ S q q ¯ g g ( 1 ) .
This soft anomalous dimension matrix also describes the time-reversed process g g q ¯ q .

9.4. q g q g and q ¯ g q ¯ g

Here we discuss quark-gluon scattering, q g q g . In the t-channel color basis c 1 q g q g = δ a 1 δ b 2 , c 2 q g q g = d b 2 c T 1 a c , c 3 q g q g = i f b 2 c T 1 a c , the one-loop soft anomalous dimension matrix is [7]
Γ S q g q g ( 1 ) = C F + C A ln ( t / s ) 0 ln ( u / s ) 0 C F ln ( t / s ) + C A 2 ln ( u / s ) C A 2 ln ( u / s ) 2 ln ( u / s ) N c 2 4 2 N c ln ( u / s ) C F ln ( t / s ) + C A 2 ln ( u / s ) .
At two loops, Γ S q g q g ( 2 ) = K 2 Γ S q g q g ( 1 ) .
This soft anomalous dimension matrix also describes the process q ¯ g q ¯ g .

9.5. g g g g

Finally, we consider gluon-gluon scattering, g g g g . The color decomposition for this process is by far the most complicated among 2 2 processes. A complete color basis for the process g g g g is given by the eight color structures [7]
c 1 g g g g = = i 4 f a b c d 12 c d a b c f 12 c , c 2 g g g g = i 4 f a b c d 12 c + d a b c f 12 c , c 3 g g g g = = i 4 f a 1 c d b 2 c + d a 1 c f b 2 c , c 4 g g g g = P 1 ( a , b ; 1 , 2 ) = 1 N c 2 1 δ a 1 δ b 2 , c 5 g g g g = P 8 S ( a , b ; 1 , 2 ) = N c N c 2 4 d a 1 c d b 2 c , c 6 g g g g = P 8 A ( a , b ; 1 , 2 ) = 1 N c f a 1 c f b 2 c , c 7 g g g g = P 10 + 10 ¯ ( a , b ; 1 , 2 ) = 1 2 ( δ a b δ 12 δ a 2 δ b 1 ) 1 N c f a 1 c f b 2 c , c 8 g g g g = P 27 ( a , b ; 1 , 2 ) = 1 2 ( δ a b δ 12 + δ a 2 δ b 1 ) 1 N c 2 1 δ a 1 δ b 2 N c N c 2 4 d a 1 c d b 2 c ,
where we used the t-channel projectors P in the product 8 8 = 1 + 8 S + 8 A + 10 + 10 ¯ + 27 describing the color content of a set of two gluons.
The one-loop soft anomalous dimension matrix is [7]
Γ S g g g g ( 1 ) = Γ 3 × 3 ( 1 ) 0 3 × 5 0 5 × 3 Γ 5 × 5 ( 1 ) ,
with
Γ 3 × 3 ( 1 ) = N c ln ( t / s ) 0 0 0 N c ln ( u / s ) 0 0 0 N c ln ( t u / s 2 )
and
Γ 5 × 5 ( 1 ) = 2 N c ln ( t s ) 0 2 N c ln ( u s ) 0 0 0 N c 2 ln ( u t 2 s 3 ) N c 2 ln ( u s ) N c ln ( u s ) 0 2 N c N c 2 1 ln ( u s ) N c 2 ln ( u s ) N c 2 ln ( u t 2 s 3 ) 0 N c ( N c + 3 ) 2 ( N c + 1 ) ln ( u s ) 0 2 N c N c 2 4 ln ( u s ) 0 N c ln ( u s ) N c ( N c + 3 ) 2 ( N c + 2 ) ln ( u s ) 0 0 2 N c ln ( u s ) ( N c + 1 ) ( 2 N c ) N c ln ( u s ) ( N c + 1 ) ln ( u s ) 2 ln ( t s )
while at two loops we have Γ S g g g g ( 2 ) = K 2 Γ S g g g g ( 1 ) .

10. Γ S for Some 2 3 Processes

In this section we consider several processes that involve a three-particle final state at leading order.

10.1. t q H , t q Z , t q γ , t q W Production

We begin with processes with three-particle final states each involving a top quark produced in association with a Higgs boson, a photon, or a W or Z-boson [11].
We begin with the s-channel processes q ( p a ) + q ¯ ( p b ) t ( p 1 ) + b ¯ ( p 2 ) + H ( p 3 ) and q q ¯ t b ¯ Z , q q ¯ t b ¯ γ , q q ¯ t b ¯ W , q q ¯ t q ¯ W + . We define s = ( p a + p b ) 2 , t = ( p a p 1 ) 2 , and u = ( p b p 1 ) 2 , as before, and further define s = ( p 1 + p 2 ) 2 , t = ( p b p 2 ) 2 , and u = ( p a p 2 ) 2 . The soft anomalous dimension matrix is identical for all these processes. We choose q q ¯ t b ¯ H as the specific process, and we use the color basis c 1 q q ¯ t b ¯ H = δ a b δ 12 and c 2 q q ¯ t b ¯ H = T b a c T 12 c . Then, the four elements of the soft anomalous dimension matrix, Γ S q q ¯ t b ¯ H , are given at one loop by [11]
Γ 11 q q ¯ t b ¯ H ( 1 ) = C F ln s m t 2 m t s 1 2 , Γ 12 q q ¯ t b ¯ H ( 1 ) = C F 2 N c ln t ( t m t 2 ) u ( u m t 2 ) , Γ 21 q q ¯ t b ¯ H ( 1 ) = ln t ( t m t 2 ) u ( u m t 2 ) , Γ 22 q q ¯ t b ¯ H ( 1 ) = C F ln s m t 2 m t s 1 2 1 N c ln t ( t m t 2 ) u ( u m t 2 ) + N c 2 ln t ( t m t 2 ) s ( s m t 2 ) .
We note that this is very similar to s-channel single-top production. since in both cases we have two colored particles in the final state, the difference being an extra colorless boson in the case here. Thus, the soft anomalous dimension matrices are almost the same, the difference arising from the more complicated kinematics in t q H production—essentially, by replacing s with s , t with t , and u with u in selected places.
We continue with the t-channel processes b ( p a ) + q ( p b ) t ( p 1 ) + q ( p 2 ) + H ( p 3 ) and b q t q Z , b q t q γ , b q t q W , q q t q W + , which have the same soft anomalous dimension matrix. We define the kinematical variables as above. We choose b q t q H as the specific process, and we use the color basis c 1 b q t q H = δ a 1 δ b 2 and c 2 b q t q H = T 1 a c T 2 b c . Then, the four elements of the soft anomalous dimension matrix, Γ S b q t q H , are given at one loop by [11]
Γ 11 b q t q H ( 1 ) = C F ln t ( t m t 2 ) m t s 3 / 2 1 2 , Γ 12 b q t q H ( 1 ) = C F 2 N c ln u ( u m t 2 ) s ( s m t 2 ) , Γ 21 b q t q H ( 1 ) = ln u ( u m t 2 ) s ( s m t 2 ) , Γ 22 b q t q H ( 1 ) = C F ln t ( t m t 2 ) m t s 3 / 2 1 2 1 N c ln u ( u m t 2 ) s ( s m t 2 ) + N c 2 ln u ( u m t 2 ) t ( t m t 2 ) .
We note that this is very similar to t-channel single-top production, and the soft anomalous dimension matrices are almost the same, essentially differing by replacing s with s , t with t , and u with u in selected places.
At two loops, the soft anomalous dimension matrices for each of these s-channel or t-channel processes can be written compactly in terms of the corresponding one-loop results [11], in a way entirely analogous to the s-channel and t-channel single-top results in Section 7, i.e., as in Equations (79) and (82).

10.2. t t ¯ H , t t ¯ Z , t t ¯ γ , t t ¯ W Production

We next consider the processes q ( p a ) + q ¯ ( p b ) t ( p 1 ) + t ¯ ( p 2 ) + H ( p 3 ) and q q ¯ t t ¯ Z , q q ¯ t t ¯ γ , q q ¯ t t ¯ W ± , which have the same soft anomalous dimension matrix. We choose q q ¯ t t ¯ H as the specific process and use a color tensor basis of s-channel singlet and octet exchange, c 1 q q ¯ t t ¯ H = δ a b δ 12 , c 2 q q ¯ t t ¯ H = T b a c T 12 c . The four matrix elements of Γ S q q ¯ t t ¯ H are closely related to those for q q ¯ t t ¯ [3,5,58] that we presented in Section 8.1.1, and are given at one loop [62,63] by
Γ 11 q q ¯ t t ¯ H ( 1 ) = C F L β + 1 , Γ 12 q q ¯ t t ¯ H ( 1 ) = C F 2 N c Γ 21 q q ¯ t t ¯ H ( 1 ) , Γ 21 q q ¯ t t ¯ H ( 1 ) = ln ( t m t 2 ) ( t m t 2 ) ( u m t 2 ) ( u m t 2 ) , Γ 22 q q ¯ t t ¯ H ( 1 ) = C F C A 2 L β 1 + 2 ln ( t m t 2 ) ( t m t 2 ) ( u m t 2 ) ( u m t 2 ) + C A 2 ln ( t m t 2 ) ( t m t 2 ) s m t 2 1
where L β is of the form of Equation (57) but with β replaced by β = 1 4 m 2 / s .
The processes g g t t ¯ H , g g t t ¯ Z , g g t t ¯ γ , have the same soft anomalous dimension which is a 3 × 3 matrix of the form of Equation (96). The matrix elements are closely related to those for g g t t ¯ [5,58] that we presented in Section 8.1.2. We choose g g t t ¯ H as the specific process and use the color basis c 1 g g t t ¯ H = δ a b δ 12 , c 2 g g t t ¯ H = d a b c T 12 c , c 3 g g t t ¯ H = i f a b c T 12 c . At one loop we have [62,63]
Γ 11 g g t t ¯ H ( 1 ) = C F L β + 1 , Γ 13 g g t t ¯ H ( 1 ) = 1 2 ln ( t m t 2 ) ( t m t 2 ) ( u m t 2 ) ( u m t 2 ) , Γ 22 g g t t ¯ H ( 1 ) = C F C A 2 L β 1 + C A 2 1 2 ln ( t m t 2 ) ( t m t 2 ) ( u m t 2 ) ( u m t 2 ) s 2 m t 4 1 , Γ 31 g g t t ¯ H ( 1 ) = 2 Γ 13 g g t t ¯ H ( 1 ) , Γ 23 g g t t ¯ H ( 1 ) = C A 2 Γ 13 g g t t ¯ H ( 1 ) , Γ 32 g g t t ¯ H ( 1 ) = ( N c 2 4 ) 2 N c Γ 13 g g t t ¯ H ( 1 ) .

10.3. q q ¯ g , Q Q ¯ g , and g g g Final States

Soft anomalous dimension matrices at one loop for processes with three colored particles in the final state have appeared in references [64,65].
The soft anomalous dimension for the process q q ¯ q q ¯ g is a 4 × 4 matrix, for the process g g q q ¯ g it is an 11 × 11 matrix, and for the process g g g g g it is a 22 × 22 matrix, with details given in reference [64].
Results for related processes involving heavy quarks were given in reference [65]. The soft anomalous dimension for the process q q ¯ Q Q ¯ g is again a 4 × 4 matrix, and for the process g g Q Q ¯ g it is again an 11 × 11 matrix, with details given in reference [65].

11. Summary and Conclusions

Soft-gluon resummation provides a powerful method to calculate large and often dominant higher-order corrections in perturbative cross sections (see reference [12] for numerical results for many processes). Soft anomalous dimensions are essential in performing resummation beyond leading-logarithm accuracy, and in general, they are matrices in the space of color exchanges.
One-loop results for soft anomalous dimensions are available for virtually all 2 2 processes and many 2 3 processes. Two-loop results and even three-loop results are also known for many 2 2 processes and some 2 3 ones. We have reviewed these results using a consistent approach and consistent terminology in all cases. We have provided comprehensive and detailed expressions for a large number of 2 2 processes involving single-top and top-pair production, electroweak-boson and Higgs production, jet production, and other SM and BSM processes.
We have also provided results for soft anomalous dimensions for a number of 2 3 processes involving the production of single top quarks or top-antitop pairs in association with electroweak or Higgs bosons, and discussed processes with three final-state colored particles.
These results can be used and have been used for performing resummation and for calculating soft-gluon corrections at higher orders for a very large number of processes.

Funding

This research was funded by the National Science Foundation under grant number PHY 1820795.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Sterman, G. Summation of large corrections to short distance hadronic cross-sections. Nucl. Phys. B 1987, 281, 310–364. [Google Scholar] [CrossRef]
  2. Catani, S.; Trentadue, L. Resummation of the QCD perturbative series for hard processes. Nucl. Phys. B 1989, 327, 323–352. [Google Scholar] [CrossRef]
  3. Kidonakis, N.; Sterman, G. Subleading logarithms in QCD hard scattering. Phys. Lett. B 1996, 387, 867–874. [Google Scholar] [CrossRef]
  4. Contopanagos, H.; Laenen, E.; Sterman, G. Sudakov factorization and resummation. Nucl. Phys. B 1997, 484, 303–327. [Google Scholar] [CrossRef] [Green Version]
  5. Kidonakis, N.; Sterman, G. Resummation for QCD hard scattering. Nucl. Phys. B 1997, 505, 321–348. [Google Scholar] [CrossRef] [Green Version]
  6. Kidonakis, N.; Oderda, G.; Sterman, G. Threshold resummation for dijet cross sections. Nucl. Phys. B 1998, 525, 299–332. [Google Scholar] [CrossRef] [Green Version]
  7. Kidonakis, N.; Oderda, G.; Sterman, G. Evolution of color exchange in QCD hard scattering. Nucl. Phys. B 1998, 531, 365–402. [Google Scholar] [CrossRef] [Green Version]
  8. Laenen, E.; Oderda, G.; Sterman, G. Resummation of threshold corrections for single particle inclusive cross-sections. Phys. Lett. B 1998, 438, 173–183. [Google Scholar] [CrossRef] [Green Version]
  9. Aybat, S.M.; Dixon, L.J.; Sterman, G. The Two-loop anomalous dimension matrix for soft gluon exchange. Phys. Rev. Lett. 2006, 97, 072001. [Google Scholar] [CrossRef] [Green Version]
  10. Hinderer, P.; Ringer, F.; Sterman, G.; Vogelsang, W. Threshold resummation at NNLL for single-particle production in hadronic collisions. Phys. Rev. D 2019, 99, 054019. [Google Scholar] [CrossRef] [Green Version]
  11. Forslund, M.; Kidonakis, N. Resummation for 2→n processes in single-particle-inclusive kinematics. Phys. Rev. D 2020, 102, 034006. [Google Scholar] [CrossRef]
  12. Kidonakis, N. Soft-gluon corrections in top-quark production. Int. J. Mod. Phys. A 2018, 33, 1830021. [Google Scholar] [CrossRef]
  13. Kidonakis, N.; Pecjak, B.D. Top-quark production and QCD. Eur. Phys. J. C 2012, 72, 2084. [Google Scholar] [CrossRef] [Green Version]
  14. Gross, D.J.; Wilczek, F. Ultraviolet behavior of nonabelian gauge theories. Phys. Rev. Lett. 1973, 30, 1343–1346. [Google Scholar] [CrossRef] [Green Version]
  15. Politzer, H.D. Reliable perturbative results for strong interactions? Phys. Rev. Lett. 1973, 30, 1346–1349. [Google Scholar] [CrossRef] [Green Version]
  16. Caswell, W.E. Asymptotic behavior of nonabelian gauge theories to two-loop order. Phys. Rev. Lett. 1974, 33, 244–246. [Google Scholar] [CrossRef]
  17. Jones, D.R.T. Two-loop diagrams in Yang-Mills theory. Nucl. Phys. B 1974, 75, 531–538. [Google Scholar] [CrossRef]
  18. Egorian, E.; Tarasov, O.V. Two-loop renormalization of the QCD in an arbitrary gauge. Theor. Math. Phys. 1979, 41, 863–867. [Google Scholar]
  19. Tarasov, O.V.; Vladimirov, A.A.; Zharkov, A.Y. The Gell-Mann-Low function of QCD in the three-loop approximation. Phys. Lett. B 1980, 93, 429–432. [Google Scholar] [CrossRef]
  20. Larin, S.A.; Vermaseren, J.A.M. The three-loop QCD β-function and anomalous dimensions. Phys. Lett. B 1993, 303, 334–336. [Google Scholar] [CrossRef] [Green Version]
  21. Van Ritbergen, T.; Vermaseren, J.A.M.; Larin, S.A. The four-loop beta-function in Quantum Chromodynamics. Phys. Lett. B 1997, 400, 379–384. [Google Scholar] [CrossRef] [Green Version]
  22. Herzog, F.; Ruijl, B.; Ueda, T.; Vermaseren, J.A.M.; Vogt, A. The five-loop beta function of Yang-Mills theory with fermions. J. High Energy Phys. 2017, 2017, 090. [Google Scholar] [CrossRef] [Green Version]
  23. Kidonakis, N. Two-loop soft anomalous dimensions and next-to-next-to-leading-logarithm resummation for heavy quark production. Phys. Rev. Lett. 2009, 102, 232003. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kodaira, J.; Trentadue, L. Summing soft emission in QCD. Phys. Lett. B 1982, 112, 66–70. [Google Scholar] [CrossRef]
  25. Moch, S.; Vermaseren, J.A.M.; Vogt, A. The three-loop splitting functions in QCD: The non-singlet case. Nucl. Phys. B 2004, 688, 101–134. [Google Scholar] [CrossRef] [Green Version]
  26. Henn, J.M.; Korchemsky, G.P.; Mistlberger, B. The full four-loop cusp anomalous dimension in N=4 super Yang-Mills and QCD. J. High Energy Phys. 2020, 2004, 018. [Google Scholar] [CrossRef] [Green Version]
  27. Von Manteuffel, A.; Panzer, E.; Schabinger, R.M. Cusp and collinear anomalous dimensions in four-loop QCD from form factors. Phys. Rev. Lett. 2020, 124, 162001. [Google Scholar]
  28. Moch, S.; Vogt, A. Higher-order soft corrections to lepton pair and Higgs boson production. Phys. Lett. B 2005, 631, 48–57. [Google Scholar] [CrossRef] [Green Version]
  29. Floratos, E.G.; Ross, D.A.; Sachrajda, C.T. Higher-order effects in asymptotically free gauge theories: The anomalous dimensions of Wilson operators. Nucl. Phys. B 1977, 129, 66–88. [Google Scholar] [CrossRef]
  30. Floratos, E.G.; Ross, D.A.; Sachrajda, C.T. Higher-order effects in asymptotically free gauge theories: (II). flavour singlet Wilson operators and coefficient functions. Nucl. Phys. B 1979, 152, 493–520. [Google Scholar] [CrossRef] [Green Version]
  31. Gonzalez-Arroyo, A.; Lopez, C.; Yndurain, F.J. Second-order contributions to the structure functions in deep inelastic scattering (I). Theoretical calculations. Nucl. Phys. B 1979, 153, 161–186. [Google Scholar] [CrossRef]
  32. Curci, G.; Furmanski, W.; Petronzio, R. Evolution of parton densities beyond leading order: The nonsinglet case. Nucl. Phys. B 1980, 175, 27–92. [Google Scholar] [CrossRef] [Green Version]
  33. Furmanski, W.; Petronzio, R. Singlet parton densities beyond leading order. Phys. Lett. B 1980, 97, 437–442. [Google Scholar] [CrossRef] [Green Version]
  34. Moch, S.; Vermaseren, J.A.M.; Vogt, A. Higher-Order Corrections in Threshold Resummation. Nucl. Phys. B 2005, 726, 317–335. [Google Scholar] [CrossRef] [Green Version]
  35. Kidonakis, N. A unified approach to NNLO soft and virtual corrections in electroweak, Higgs, QCD, and SUSY processes. Int. J. Mod. Phys. A 2004, 19, 1793–1821. [Google Scholar] [CrossRef] [Green Version]
  36. Becher, T.; Neubert, M. Infrared singularities of scattering amplitudes in perturbative QCD. Phys. Rev. Lett. 2009, 102, 162001. [Google Scholar] [CrossRef] [Green Version]
  37. Gardi, E.; Magnea, L. Factorization constraints for soft anomalous dimensions in QCD scattering amplitudes. J. High Energy Phys. 2009, 2009, 079. [Google Scholar] [CrossRef]
  38. Almelid, O.; Duhr, C.; Gardi, E. Three-loop corrections to the soft anomalous dimension in multi-leg scattering. Phys. Rev. Lett. 2016, 117, 172002. [Google Scholar] [CrossRef]
  39. Kidonakis, N. Next-to-next-to-next-to-leading-order soft-gluon corrections in hard-scattering processes near threshold. Phys. Rev. D 2006, 73, 034001. [Google Scholar] [CrossRef] [Green Version]
  40. Polyakov, A.M. Gauge fields as rings of glue. Nucl. Phys. B 1980, 164, 171–188. [Google Scholar] [CrossRef]
  41. Brandt, R.A.; Neri, F.; Sato, M. Renormalization of loop functions for all loops. Phys. Rev. D 1981, 24, 879–902. [Google Scholar] [CrossRef]
  42. Ivanov, S.V.; Korchemsky, G.P.; Radyushkin, A.V. Infrared asymptotics of perturbative QCD. Contour gauges. Sov. J. Nucl. Phys. 1986, 44, 145–150. [Google Scholar]
  43. Korchemsky, G.P.; Radyushkin, A.V. Loop-space formalism and renormalization group for the infrared asymptotics of QCD. Phys. Lett. B 1986, 171, 459–467. [Google Scholar] [CrossRef]
  44. Grozin, A.; Henn, J.M.; Korchemsky, G.P.; Marquard, P. The three-loop cusp anomalous dimension in QCD. Phys. Rev. Lett. 2015, 114, 062006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Kidonakis, N. Three-loop cusp anomalous dimension and a conjecture for n loops. Int. J. Mod. Phys. A 2016, 31, 1650076. [Google Scholar] [CrossRef] [Green Version]
  46. Remiddi, E.; Vermaseren, J.A.M. Harmonic polylogarithms. Int. J. Mod. Phys. A 2000, 15, 725–754. [Google Scholar] [CrossRef] [Green Version]
  47. Kidonakis, N.; Gonsalves, R.J. Higher-order QCD corrections for the W-boson transverse momentum distribution. Phys. Rev. D 2013, 87, 014001. [Google Scholar] [CrossRef] [Green Version]
  48. Kidonakis, N. Single top quark production at the Fermilab Tevatron: Threshold resummation and finite-order soft gluon corrections. Phys. Rev. D 2006, 74, 114012. [Google Scholar] [CrossRef] [Green Version]
  49. Kidonakis, N. Next-to-next-to-leading logarithm resummation for s-channel single top quark production. Phys. Rev. D 2010, 81, 054028. [Google Scholar] [CrossRef] [Green Version]
  50. Kidonakis, N. Soft anomalous dimensions for single-top production at three loops. Phys. Rev. D 2019, 99, 074024. [Google Scholar] [CrossRef] [Green Version]
  51. Kidonakis, N. Next-to-next-to-leading-order collinear and soft gluon corrections for t-channel single top quark production. Phys. Rev. D 2011, 83, 091503. [Google Scholar] [CrossRef] [Green Version]
  52. Kidonakis, N. Two-loop soft anomalous dimensions for single top quark associated production with a W- or H-. Phys. Rev. D 2010, 82, 054018. [Google Scholar] [CrossRef] [Green Version]
  53. Belyaev, A.; Kidonakis, N. QCD corrections to flavor-changing-neutral-current single top production at the DESY HERA collider. Phys. Rev. D 2002, 65, 037501. [Google Scholar] [CrossRef] [Green Version]
  54. Kidonakis, N.; Belyaev, A. FCNC top quark production via anomalous tqV couplings beyond leading order. J. High Energy Phys. 2003, 2003, 004. [Google Scholar] [CrossRef] [Green Version]
  55. Kidonakis, N.; Martin, E. Soft-gluon corrections in FCNC top-quark production via anomalous gluon couplings. Phys. Rev. D 2014, 90, 054021. [Google Scholar] [CrossRef] [Green Version]
  56. Ferroglia, A.; Neubert, M.; Pecjak, B.D.; Yang, L.L. Two-loop divergences of scattering amplitudes with massive partons. Phys. Rev. Lett. 2009, 103, 201601. [Google Scholar] [CrossRef] [Green Version]
  57. Ferroglia, A.; Neubert, M.; Pecjak, B.D.; Yang, L.L. Two-loop divergences of massive scattering amplitudes in non-abelian gauge theories. J. High Energy Phys. 2009, 2009, 062. [Google Scholar] [CrossRef]
  58. Kidonakis, N. Next-to-next-to-leading soft-gluon corrections for the top quark cross section and transverse momentum distribution. Phys. Rev. D 2010, 82, 114030. [Google Scholar] [CrossRef] [Green Version]
  59. Laenen, E.; Moch, S. Soft gluon resummation for heavy quark electroproduction. Phys. Rev. D 1999, 59, 034027. [Google Scholar] [CrossRef] [Green Version]
  60. Kulesza, A.; Motyka, L. Threshold resummation for squark-antisquark and gluino-pair production at the LHC. Phys. Rev. Lett. 2009, 102, 111802. [Google Scholar] [CrossRef] [Green Version]
  61. Beenakker, W.; Brensing, S.; Kramer, M.; Kulesza, A.; Laenen, E.; Niessen, I. Soft-gluon resummation for squark and gluino hadroproduction. J. High Energy Phys. 2009, 2009, 041. [Google Scholar] [CrossRef] [Green Version]
  62. Li, H.T.; Li, C.S.; Li, S.A. Renormalization group improved predictions for tt¯W± production at hadron colliders. Phys. Rev. D 2014, 90, 094009. [Google Scholar] [CrossRef] [Green Version]
  63. Kulesza, A.; Motyka, L.; Stebel, T.; Theeuwes, V. Soft gluon resummation for associated tt¯H production at the LHC. J. High Energy Phys. 2016, 2016, 065. [Google Scholar] [CrossRef] [Green Version]
  64. Sjodahl, M. Color evolution of 2→3 processes. J. High Energy Phys. 2008, 2008, 083. [Google Scholar] [CrossRef] [Green Version]
  65. Szarek, E. Soft anomalous dimension matrices in heavy quark-antiquark hadroproduction in association with a gluon jet. Acta Phys. Polon. B 2018, 49, 1839–1864. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Kidonakis, N. Soft Anomalous Dimensions and Resummation in QCD. Universe 2020, 6, 165. https://doi.org/10.3390/universe6100165

AMA Style

Kidonakis N. Soft Anomalous Dimensions and Resummation in QCD. Universe. 2020; 6(10):165. https://doi.org/10.3390/universe6100165

Chicago/Turabian Style

Kidonakis, Nikolaos. 2020. "Soft Anomalous Dimensions and Resummation in QCD" Universe 6, no. 10: 165. https://doi.org/10.3390/universe6100165

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop