Multiple soft radiation at one-loop order and the emission of a soft quark-antiquark pair

We consider the radiation of two or more soft partons in QCD hard-scattering at one-loop order. The corresponding scattering amplitude is singular, and the singular behaviour is controlled by a process-independent soft current. Using regularization in $d=4 - 2\epsilon$ space-time dimensions, we explicitly evaluate the ultraviolet and infrared divergent ($\epsilon$-pole) terms of the one-loop soft current for emission of an arbitrary number of soft partons in a generic hard-scattering process. Then we consider the specific case of soft quark-antiquark ($q{\bar q}$) emission and we compute the one-loop current by including the finite terms. We find that the one-loop soft-$q{\bar q}$ current exhibits a new type of transverse-momentum singularity, which has a quantum (absorptive) origin and a purely non-abelian character. At the squared amplitude (cross section) level, this transverse-momentum singularity produces contributions to multijet production processes in hadron collisions. The one-loop squared current also leads to charge asymmetry terms, which are a distinctive features of soft-$q{\bar q}$ radiation. We also extend these results to the cases of QED and mixed QCD$\times$QED radiative corrections for soft fermion-antifermion emission.


Introduction
The physics program carried out at the Large Hadron Collider (LHC) has already produced an impressive amount of high-precision data, and similar data will be obtained in the next runs of the LHC. Theoretical predictions are thus demanded to achieve a corresponding high precision.
In the context of the perturbative evaluation of QCD radiative corrections, the present high-precision frontier is represented by computations at the next-to-next-to-next-to-leading order (N 3 LO) in the QCD coupling α S . Some N 3 LO results for LHC processes are already available (see, e.g., related references in Ref. [1]). In the case of observables that are highly sensitive to multiple radiation of soft and collinear partons, the fixed-order QCD predictions have to be supplemented with the all-order resummed calculations of classes of large logarithmic contributions. In few specific cases (see, e.g., Refs. [2,3,4,5,6,7,8,9,10,11,12,13,14]) resummed QCD calculations have reached the next-to-next-to-next-to-leading logarithmic (N 3 LL) accuracy.
An important feature of QCD scattering amplitudes is the presence of singularities in soft and collinear regions of the phase space, and the corresponding presence of infrared (IR) divergences in virtual contributions at the loop level. The soft and collinear singularities have a process-independent structure, and they are controlled by universal factorization formulae and corresponding soft/collinear factors. As briefly recalled below, these factorization properties are relevant for both fixed-order and resummed QCD calculations.
In the computation of physical observables for hard-scattering processes, phase space soft/collinear singularities and virtual IR divergences cancel between themselves. However, much technical effort is required to achieve and implement the cancellation, and the effort highly increases by increasing the perturbative order. Soft/collinear factorization formulae can be used to organize and greatly simplify the cancellation mechanism of the IR divergences in fixed-oder calculations.
In the evaluation of observables close to the exclusive boundary of the phase space, real and virtual radiative corrections in the scattering amplitudes are kinematically strongly unbalanced. As a consequence, the cancellation mechanism of the IR divergences leaves residual effects in the form of large logarithmic contributions. Soft/collinear factorization formulae and the corresponding singular factors are the basic ingredients for the explicit computation and resummation of these large logarithmic contributions.
The singular factors at O(α S ) and O(α 2 S ) for soft and collinear factorization of scattering amplitudes are known since long time. The explicit knowledge of soft/collinear factorization at O(α S ) has been essential to devise fully general (process-independent and observableindependent) methods to carry out next-to-leading order (NLO) QCD calculations (see, e.g., Refs. [15,16]). Similarly, the knowledge of soft/collinear factorization formulae at O(α 2 S ) [17,18,19,20,21,22,23,24,25,26,27] is exploited to develop methods (see, e.g., the review in Ref. [1]) at the next-to-next-to-leading order (NNLO). Soft/collinear factorization up to O(α 2 S ) contributes to resummed calculations up to next-to-next-toleading logarithmic (NNLL) accuracy (see, e.g., Refs. [28,29]). Soft and collinear factorization at O(α 3 S ) can be used in the context of N 3 LO calculations and of resummed calculations at N 3 LL accuracy. The process-independent singular factors for the various collinear limits at O(α 3 S ) are presented in Refs. [30,31,32,33,34,35,36,37,38,26]. Soft factorization of scattering amplitudes at O(α 3 S ) requires the study of various tree-level and loop contributions. Triple soft-gluon radiation at the tree level is studied in Ref. [39]. Double soft emission at one loop level has been considered recently in Ref. [40]. Single soft-gluon radiation at two loop order is examined in detail in Refs. [37,41,42,43]. This paper is devoted to a study of soft-parton emission at O(α 3 S ) and beyond this order. More precisely, we consider the singular behaviour of scattering amplitudes in the limit in which two or more external partons are soft. The singularity is controlled in factorized form by a current for soft multiparton radiation from hard partons. At one-loop order the soft current contains IR and ultraviolet (UV) divergent contributions that we explicitly evaluate for the emission of an arbitrary number of soft partons. In the particular case of emission of a soft quark-antiquark (qq) pair, we explicitly compute also the finite contributions to the one-loop current. We comment on the related results of Ref. [40] in the paper.
The outline of the paper is as follows. In Sect. 2 we introduce our notation, and we recall the soft factorization formula for scattering amplitudes and the known results on the tree-level currents for emission of a single soft gluon and of soft-qq pair. We use analytic continuation in d = 4 − 2ǫ space-time dimensions to regularize IR and UV divergences in loop contributions. In Sect. 3 we discuss general features of the current for multiple soft radiation at the loop level. In particular, we present in explicit form the result of the IR and UV divergent (ǫ-pole) terms of the one-loop soft current. In Sect. 4 we consider the emission of a soft-qq pair and we compute the corresponding one-loop current by including the finite (i.e., O(ǫ 0 )) terms. We comment on general features of our result that is valid for generic multiparton scattering processes in arbitrary kinematical configurations. Section 5 is devoted to consider soft-qq radiation at the squared amplitude level. We first recall the results for the squared current at the tree level, and then we explicitly compute the one-loop squared current. We discuss the structure of the charge asymmetry contributions, which are a distinctive feature of soft-qq radiation at the loop level. In Sect. 5.3 we present simplified expressions for processes with two and three hard partons. In Sect. 6 we generalize our QCD results for soft qq emission to the cases of QED and mixed QCD×QED radiative corrections for soft fermion-antifermion emission. A brief summary of our results is presented in Sect. 7.

Soft factorization
We consider the amplitude (the S-matrix element) M of a generic scattering process whose external particles (the external legs of M) are QCD partons (quarks, antiquarks and gluons) and, possibly, additional non-QCD particles (i.e., partons with no colour charge such as leptons, Higgs and electroweak vector bosons and so forth). We use the notation M(p 1 , p 2 , . . . , p n ), where p i (i = 1, . . . , n) is the momentum of the QCD parton A i (A i = g, q orq ). Unless otherwise specified, the dependence of M on the momenta (and quantum numbers) of additional colourless particles is not explicitly denoted.
The external QCD partons are on-shell with physical spin polarizations (thus, M in-cludes the corresponding spin wave functions), and we always define the external momenta p i 's as outgoing momenta. Note, however, that we do no restrict our treatment to processes with physical partons in the final state. In particular, the time-component (i.e. the 'energy') p 0 i of the momentum vector p ν i (ν = 0, 1, . . . , d − 1) in d space-time dimensions is not positive definite. Different types of physical processes are described by considering different kinematical regions of the parton momenta and by simply applying crossing symmetry to the wave functions and quantum numbers of the external partons of the same matrix element M(p 1 , p 2 , . . . , p n ). According to our definition of the momenta, if p i has positive energy, M(. . . , p i , . . . ) describes a physical process that produces the parton A i in the final state; if p i has negative energy, M(. . . , p i , . . . ) describes a physical process produced by the collision of the antiparton A i in the initial state.
It is convenient to directly work in colour (and spin) space, and to use the notation of Ref. [16] (see also Ref. [44]). We treat the colour and spin structures by formally introducing an orthonormal basis {|c 1 , c 2 , . . . , c n ⊗ |s 1 , s 2 , . . . , s n } in colour + spin space. The scattering amplitude can be written as Thus |M(p 1 , p 2 , . . . , p n ) is a vector in colour + spin (helicity) space.
As previously stated, we define the external momenta p i 's as outgoing momenta. The colour indices {c 1 , c 2 , . . . c n } are consistently defined as outgoing colour indices: c i is the colour index of the parton A i with outgoing momentum p i (if p i has negative energy, c i is the colour index of the physical parton A i that collides in the initial state). An analogous comment applies to spin indices.
The amplitude M can be evaluated in QCD perturbation theory as a power series expansion (i.e., loop expansion) in the QCD coupling g S (or, equivalently, in the strong coupling α S = g 2 S /(4π)). We write where M (0) is the tree-level scattering amplitude, M (1) is the one-loop contribution, M (2) is the two-loop contribution, and so forth. More generally, M (0) is not necessarily a tree amplitude, but rather the lowest-order amplitude for that given process. Thus, M (L) (L = 1, 2 . . . ) is the corresponding L-loop correction. For instance, in the cases of the diphoton production process gg → γγ or the Higgs boson (H) production process gg → H, the corresponding amplitude M (0) involves a quark loop (a massive-quark loop in the case of gg → H). Note that in Eq. (3) we have not made explicit the dependence on powers of g S . Thus, M (0) includes an integer power of g S as overall factor, and M (1) includes an extra factor of g 2 S (i.e., M (1) /M (0) ∝ g 2 S ).
Physical processes take place in four-dimensional space-time. The four-dimensional evaluation of the L-loop amplitude M (L) leads to UV and IR divergences that have to be properly regularized. We regularize both kind of divergences by performing the analytic continuation of the loop momenta and phase-space in d = 4−2ǫ space-time dimensions. We postpone comments on different variants of dimensional regularization. The dimensionalregularization scale is denoted by µ. After regularization, the UV and IR divergences appears as ǫ-poles of the Laurent series expansion in powers of ǫ around ǫ = 0. Throughout the paper we formally consider expressions for arbitrary values of d = 4−2ǫ (equivalently, in terms of ǫ expansions, the expressions are valid to all orders in ǫ before they are eventually truncated at some order in ǫ). We always consider unrenormalized amplitudes, and g S denotes the bare (unrenormalized) coupling constant.
We are interested in studying the behaviour of M in the kinematical configuration where one or more of the momenta of the external massless partons (gluons or massless quark and antiquarks) become soft. In this kinematical configuration, M becomes singular. To make the notation more explicit, the soft momenta are denoted by q ν k , while the other parton momenta are still denoted by p ν i . The behaviour of M(. . . , q k , . . . , p i , . . . ) in this multiparton soft region is formally specified by performing an overall rescaling of all soft momenta as q k → λq k (the rescaling parameter λ is the same for each soft momentum q k ) and by considering the limit λ → 0. In this limit, if the set of soft partons has m (m ≥ 1) momenta q k 's (k = 1, . . . , m), the amplitude M is singular and it behaves as The power-like behaviour (λ) −m that we have specified in the right-hand side of Eq. (4) determines the dominant singular terms of M in the multiple soft region. The logarithmic corrections ln r λ (r = 0, 1, 2, . . . ) arise from scaling violation, since the naïve (power-like) scaling behaviour is violated by the effects of the UV and IR divergences of the scattering amplitude at the loop level (see Sect. 3). The dots on the right-hand side of Eq. (4) denote the subdominant singular behaviour of M. The relative suppression factor between subdominant and dominant terms is (at least) of O( √ λ).
The computation of physical observables eventually requires the phase-space integration of the squared amplitude |M| 2 . We note that, after phase-space integration over the soft momenta, the dominant singular behaviour of |M| 2 produces logarithmic soft (IR) divergences (i.e., ǫ poles), whereas the subdominant singular behaviour does not lead to soft divergences. In this paper we are interested in the dominant singular behaviour of Eq. (4).
In the case of tree-level scattering amplitudes [47,48,22], the factorization formula (5) can be simply derived by considering soft-parton radiation from the hard-parton (the partons with momenta p 1 , . . . , p n ) external legs of the amplitude and by directly applying the eikonal approximation for emission vertices and propagators. At the one-loop level, the factorization structure of Eq. (5) was worked out in Refs. [45,21,23]. In particular, as discussed in detail in Ref. [23], the one-loop soft current can still be computed by using the eikonal approximation for soft-parton radiation from the external hard partons, and this discussion generalizes to two-loop and higher-loop orders. Owing to its origin by eikonal radiation from the hard partons, the soft current can also straightforwardly be expressed in equivalent form as matrix element of Wilson line operators [46].
The soft current J (q 1 , . . . , q m ) depends on the soft partons, specifically on their momenta and their quantum numbers (flavour, spin, colour), and it also depends on the hard partons (on their momenta and their quantum numbers), though we use a customary notation in which the dependence of J on p 1 , . . . , p n is not explicitly denoted in its argument. The current J is an operator (a 'rectangular' matrix) that acts from the (lower-dimensional) colour+spin space of the hard partons to the (higher-dimensional) colour+spin space of the soft and hard partons. We remark on the fact that the soft current J is simply proportional to the unit operator in the spin subspace of the hard partons, since soft radiation is insensitive to the spin of the hard radiating partons.
In spite of its dependence on hard partons, the soft current J is completely universal, namely, it does not depend on the specific scattering amplitude M and on its corresponding specific physical process. The universality of J also implies that it is directly applicable in contexts that do not directly refer to the soft behaviour of scattering amplitudes. For instance, J (or, more specifically, J † J ) is precisely the integrand of any specific soft function (see, e.g., Ref. [28] and references therein) that can be introduced through soft-collinear effective theory (SCET) [49] methods.
The colour-space factorization formula (5) does not require any specifications about the detailed colour structure of the scattering amplitudes in its left-hand and right-hand sides. Scattering amplitudes can be decomposed in a form that factorizes the QCD colour from colourless kinematical coefficients, which are colour-ordered subamplitudes (see, e.g., Ref. [50]). Colour-ordered subamplitudes fulfil soft factorization formulae that are analogous to Eq. (5) in terms of colour-stripped (though colour-ordered) soft factors (see, e.g., Refs. [48,21]). The factorization properties of colour-order subamplitudes and the corresponding soft factors can be directly and explicitly derived from Eq. (5). To this purpose it is sufficient to insert the colour decomposition of M and the explicit colour structure of J in Eq. (5). Therefore, the colour-space factorization of Eq. (5) and soft factorization of colour-ordered subamplitudes are equivalent formulations. The advantage of Eq. (5) is that it leads to a more compact formulation, without the necessity of introducing the explicit colour decomposition of M, whose actual form depends on the specific partonic content of the amplitude (e.g., on the number of gluons and quark-antiquark pairs * ) and on the loop order. Moreover, the colour space formulation can simplify the direct computation of the soft limit of squared amplitudes (see, e.g., Sect. 5 ).
The soft current J in Eq. (5) can be evaluated in QCD perturbation theory, and it can * In particular, if the set of soft partons includes one or more quark-antiquark pairs, the scattering amplitudes in the left-hand and right-hand sides of Eq. (5) have different numbers of quark-antiquark pairs and, therefore, they have different colour decompositions. be expressed in terms of a loop expansion that is completely analogous to that in Eq. (3). We write where J (0) is the tree-level current, J (1) is the one-loop current, and so forth. Analogously to Eq. (3), the loop label L in J (L) refers to the unrenormalized current. Inserting the expansions (3) and (6) in Eq. (5) we obtain factorization formulae that are valid order-byorder in the number of loops. The soft factorization formula for tree-level (lowest-order) amplitudes is where the symbol '≃' means † that we are neglecting subdominant terms in the soft limit (i.e., the terms denoted by dots in the right-hand side of Eqs. (4) and (5)). The soft factorization formula for one-loop amplitudes is The tree-level current for the emission of a single soft gluon of momentum q ν is well known [47]: where the notation i ∈ H means that the sum extends over all hard partons (with momenta p i ) in M, ε ν (q) is the spin polarization vector of the soft gluon, and T i is the colour charge of the hard parton i.
The spin index s and the colour index a (a = 1, . . . , N 2 c − 1, for SU(N c ) QCD with N c colours) of the soft gluon can be specified by acting onto Eq. (9) in colour+spin space as in Eq. (2).
where T a i denotes the generators of SU(N c ) of the representation of the parton i. We have c i | T a i |c ′ i = (T a i ) c i c ′ i , where (T a ) cb ≡ if cab (colour-charge matrix in the adjoint representation) if the parton i is a gluon and (T a ) αβ ≡ t a αβ (colour-charge matrix in the fundamental representation, with α, β = 1, . . . , N c ) if the parton i is a quark ((T a ) αβ ≡ t a αβ = −t a βα if the parton i is an antiquark). We normalize the colour matrices such as [t a , t b ] = if abc and Tr(t a t b ) = T R δ ab with T R = 1/2. The colour-charge algebra gives We note that the colour charge operators T i fulfil some relevant properties related to colour conservation. For instance, we have [16] i∈M which follows from the fact that the scattering amplitude M is a colour-singlet state (the notation i ∈ M means that the sum in Eq. (10) extends over all external partons of M). † The symbol '≃' is used throughout the paper with the same meaning as in Eq. (7).
A property analogous to that in Eq. (10) is fulfilled by the soft multiparton current J . We have or, equivalently, where the notation k ∈ S (i ∈ H) means that the sum extends over all the soft (hard) partons in J . The relations in Eqs. (11) and (12) express the property of colour flow conservation and follow from the fact that the total colour charge is conserved in the radiation process of soft partons by hard partons. Note from Eq. (12) that the total charge of the soft partons acts on J as a rotation of its hard-parton charges. We remark that Eqs. (10)- (12) are valid to all orders in the loop expansion or, equivalently, order-by-order in QCD perturbation theory. It is straightforward to explicitly check that the soft singlegluon current in Eq. (9) fulfils the colour flow conservation property in Eq. (11).
Since scattering amplitudes are colour-singlet states, the structure of soft factorization in Eq. (5) implies that the explicit expression of the universal soft current J necessarily involves some degrees of arbitrariness. Different expressions for J are indeed permitted, provided the difference is proportional to an operator that is proportional to the total colour charge of the hard partons. Owing to Eq. (10), this degree of arbitrariness is physically harmless (it does not affect the soft behaviour of the scattering amplitude) and the ensuing different expressions of J are fully equivalent (although they are not exactly equal at the formal level, before acting onto colour-singlet states). The property of colour-flow conservation in Eq. (11) does not remove this degree of arbitrariness.
For subsequent use (and similarly to Ref. [39]) we introduce the notation where the subscript CS in the symbol cs = means that the equality between the colour operators O and O ′ (e.g., soft currents or their corresponding squared currents) is valid if these operators act (either on the left or on the right) onto colour-singlet states. The notation in Eq. (13) permits to directly relate (and equate) expressions that simply differ by contributions that are due to the physically harmless arbitrariness of the soft current J .
The soft factor for radiation of two soft gluons from tree-level colour-ordered subamplitudes with external gluons and with external gluons and an additional quark-antiquark pair was computed in Ref. [48]. The tree-level current J (0) (q 1 , q 2 ) for emission of two soft gluons in a generic scattering amplitude was given in Ref. [22]. The tree-level current for emission of three soft gluons was computed in Ref. [39].
The tree-level current for emission of a single soft quark (or antiquark) vanishes. This result is equivalent to say that the dominant singular behaviour in Eq. (4) is absent in the soft single-quark limit (the radiation of a single soft quark only produces a subdominant behaviour of O(1/ √ λ) in the right-hand side of Eq. (4)).
The tree-level current for emission of a soft qq pair was computed in Ref. [22], where the result was explicitly reported at the level of squared amplitudes (i.e., the result refers to J † J ). The corresponding result for the qq current is where we have introduced the fermionic current j ν (1, 2), The soft quark and antiquark have momenta q ν 1 and q ν 2 , respectively, and u(q) and v(q) are the customary Dirac spinors. The spin indices (s 1 and s 2 ) and the colour indices (α 1 and α 2 ) of the quark and antiquark are embodied in the colour+spin space notation of Eq. (14). Analogously to Eq. (9), we can consider ( α 1 , α 2 | ⊗ s 1 , s 2 | ) J(q 1 , q 2 ) ≡ J α 1 ,α 2 s 1 ,s 2 (q 1 , q 2 ) and we have ( 3 One-loop current for multiple soft emission: UV and IR divergences The soft singular behaviour of one-loop amplitudes (see Eq. (8)) is controlled by J (0) and by an additional new ingredient, the one-loop soft current J (1) .
The one-loop soft limit for emission of a single soft gluon was worked out independently by two groups [21,23], finding results that are in agreement. The analysis of Ref. [21] is based on the study of colour-ordered subamplitudes, while Ref. [23] considers generic scattering amplitudes. The results of Refs. [21,23] are valid for the case of massless hard partons. The generalization of the results of Ref. [23] to include massive hard partons (such as heavy quarks) was carried out in Ref. [25]. In the remaining part of this paper we limit ourselves to consider scattering amplitudes with massless hard partons.
The result of the one-loop current for single gluon emission is [23] (we explicitly write J (1) a ≡ a| J (1) , where a is the colour index of the soft gluon): where 'x − i0' denotes the customary Feynman prescription for analytic continuation in different kinematical regions (x > 0 and x < 0) and c Γ is the typical volume factor of d-dimensional one-loop integrals: We remark that Eq. (16) gives the complete result to all orders in the ǫ expansion around ǫ = 0 (equivalently, the result in arbitrary d = 4 − 2ǫ space-time dimensions).
We comment on some features of Eq. (16). The one-loop current is proportional to the structure constants f abc of the gauge group and, therefore, it is purely non-abelian. This is in agreement with the absence of one-loop corrections to the soft current for single softphoton emission in massless QED [51]. The current in Eq. (16) involves non-abelian colour correlations, if abc T b i T c j , with two hard partons. Its kinematical structure has a rational dependence on p i · ε(q)/p i · q (which is analogous to that in the tree-level current of Eq. (9)) that is only modified through logarithmic corrections by the one-loop interactions. The logarithmic corrections are due to the ǫ expansion of the last factor in the right-hand side of Eq. (16), and they are proportional to powers of ln q 2 ⊥ ij (modulo branch-cut effects), where q ⊥ ij , has a simple kinematical interpretation since it is the transverse component of the gluon momentum q with respect to the longitudinal direction singled out by the momenta p i and p j (in a reference frame in which p i and p j are back-to-back) of the colour-correlated partons.
The overall scaling behaviour of J (1) (λq) (with λ > 0) in the limit λ → 0 is proportional to (λ 2 ) −ǫ /λ = (1/λ)mod(ǫ r ln r λ), and it is in agreement with Eq. (4). In particular, we explicitly see that the ln λ-enhancement is produced by the use of dimensional regularization to avoid the IR and UV divergences in the one-loop contribution to J . Performing the ǫ expansion of Eq. (16), this IR and UV behaviour produces double (1/ǫ 2 ) and single (1/ǫ) poles near ǫ = 0.
The two-loop current for single soft-gluon emission was computed in Ref. [37] up to including contributions of O(ǫ 0 ) for the simplest case of scattering amplitudes with only two hard partons. Subsequently this result was extended up to O(ǫ 2 ) [41,42] and to all orders in ǫ [42]. The two-loop result of Refs. [37,41,42] has a structure that is very similar to the one-loop current in Eq. (16). More involved structures, in terms of both colour correlations and kinematical dependence, do appear in the general case of scattering amplitudes with three or more hard partons, and the corresponding two-loop current for single soft-gluon emission was considered and explicitly computed in Ref. [43], by including the finite contributions up to O(ǫ 0 ).
We now discuss multiple soft radiation at one-loop order. The structure of the loop-level current J for multiple soft radiation is expected to be definitely more complex (in terms of both colour and kinematical dependence) than the single soft-gluon current in Eq. (16). The presence of two or more soft partons and the ensuing dependence on their momenta increases the number of relevant kinematical invariants, which drive an increased complexity of colour and kinematical correlations (especially at high orders in the ǫ expansion). In the remaining part of this Section we deal with general properties of the soft current J with m ≥ 2 soft partons. In particular, we consider the UV and IR divergences of J and we discuss their regularization scheme dependence.
The one-loop current for multiple soft emission has (analogously to Eq. (16)) double (1/ǫ 2 ) and single (1/ǫ) pole contributions due to the presence of IR and UV divergences in the four-dimensional case (ǫ = 0). At L-loop order, the current J (L) has poles of the type 1/ǫ k with 2L ≥ k ≥ 1. These ǫ-pole contributions are directly related to the corresponding contributions to the multiparton scattering amplitudes [44,52,53,54,55,56]. The ǫ-pole contributions to the one-loop current J (1) have a general structure, whose explicit form can be directly derived from the known universal structure of the IR and UV divergences of oneloop scattering amplitudes [57,16,44]. Starting from the results in Refs. [57,16,44], the procedure to derive the ǫ-pole contributions to J (1) is completely analogous to that used in Refs. [33,26] for the study of the multiparton collinear limit of scattering amplitudes (see, in particular, Eqs. (104)-(109) in the arXiv version of Ref. [26] and replace the collinear splitting matrix Sp (1) with the soft current J (1) ). Moreover that procedure can be extended to higher-loop orders and it leads to a compact representation of the ǫ-pole contributions to J at arbitrary perturbative orders (see the analogous procedure in Sect. 6.1 and, in particular, Eq. (137) in the arXiv version of Ref. [26] and replace Sp with J ). Owing to the complete analogy with the collinear limit studied in Ref. [26], we limit ourselves to present the final results for the soft limit.
The general all-order representation of the ǫ-pole contributions to J is where J [no ǫ−poles] is obtained from J by properly subtracting its ǫ-pole part order by order in the loop expansion. Therefore at each perturbative order J [no ǫ−poles] is finite in the limit ǫ → 0 order by order in the ǫ expansion around ǫ = 0. Note that this statement refers to the loop expansion of J [no ǫ−poles] with respect to renormalized QCD coupling (the use of the renormalized QCD coupling removes ǫ poles of UV origin, which cannot be absorbed in the V factors of Eq. (19)). Obviously, at the tree level J and J [no ǫ−poles] coincides (J (0) = J (0) [no ǫ−poles] ). The colour space operator V(q 1 , . . . , q m , p 1 , . . . , p n ) is the processindependent factor [44,52,53,54,55,56] that controls the IR ǫ-pole contributions to the scattering amplitude |M(q 1 , . . . , q m , p 1 , . . . , p n ) in the left-hand side of the factorization formula (5) (at the tree level, V (0) = 1). The operator V(p 1 , . . . , p n ) to the right of J [no ǫ−poles] in Eq. (19) is the restriction of V to the scattering amplitude |M(p 1 , . . . , p n ) in the left-hand side of Eq. (5) (i.e., the amplitude with the external soft legs removed). Since the V factors in Eq. (19) only depends on the colour, flavour and momentum of the soft and hard partons, their perturbative knowledge determines in a recursive manner (i.e., order by order in the loop expansion) the ǫ-pole part of J at a given order in terms of J at lower perturbative orders. At the one-loop level, using Eq. (19) and the known expression [57,16,44] of the one-loop term V (1) of the operator V, we obtain the following expression for the ǫ-pole contributions to the soft multiparton current J (1) : where C k (T 2 k = C k ) is the Casimir coefficient of the parton k and, analogously, the coefficient γ k depends on the flavour of the parton k and, explicitly, we have where N f is the number of flavours of massless quarks. The coefficient b 0 is the first perturbative coefficient of the QCD β function, Note that, in our normalization, we have b 0 = γ g .
The various ǫ-pole terms in Eq. (20) have different origins. The single-pole term that is proportional to b 0 is of UV origin; it can be removed by renormalizing the soft current J (we recall that we are considering unrenormalized scattering amplitudes and, correspondingly, unrenormalized soft currents). The other ǫ-pole terms are of IR origin. The double-pole terms, which are proportional to the Casimir coefficients C k (C k = C F and C A for quarks and gluons, respectively), originate from one-loop contributions in which the loop momentum is nearly on-shell, very soft and parallel to the momentum of one of the soft partons involved in the current. The single-pole terms with γ k coefficients are produced by contributions in which the loop momentum is not soft, though it is nearly on-shell and parallel to the momentum of one of the external soft partons of the current. The single-pole terms with logarithmic dependence on soft-parton and hard-partons subenergies (q k · q l , p i · q k , p i · p j ) originate from configurations in which the loop momentum is very soft and at wide angle with respect to the direction of the external-leg (soft and hard) partons. Specifically, the radiative part of these terms (i.e., the real part of the logarithms) is due to a nearly onshell virtual gluon in the loop, while the absorptive part (i.e., the imaginary part of the logarithms) is due to the exchange of an off-shell Coulomb-type gluon.
The expression in the right-hand side of Eq. (20) is valid for an arbitrary number m of soft partons in the current (and for an arbitrary number of hard partons in the scattering amplitude). This expression is given in terms of explicit coefficients and of the tree-level current J (0) for the corresponding parton configuration. Once J (0) (and, in particular, its colour structure) is explicitly known, Eq. (20) can be directly applied to determine the explicit ǫ-pole contributions to the one-loop current J (1) .
In particular, in the case of a single (m = 1) soft parton, using Eq. (9) it is straightforward to check that Eq. (20) gives the ǫ-pole terms of the one-loop result J (1) in Eq. (16). In this respect, we note that the expression in the right-hand side of Eq. (20) does not identically (in its precise algebraic form) correspond to the ǫ-pole terms in Eq. (16): the difference is due to terms of O(1/ǫ 2 ) and O(1/ǫ) that are proportional to the total colour charge ( i∈H T i ) of the hard partons. As previously discussed (see Eq. (13) and related comments) the presence of such terms in J is physically harmless.
The tree-level currents J (0) for emission of two (m = 2) soft partons (either two gluons or a qq pair) are also explicitly known [22]. Therefore Eq. (20) can also be straightforwardly applied to explicitly obtain the ǫ-pole terms of the one-loop current J (1) for double softparton emission. The case of a soft quark and antiquark is discussed in detail in Sect. 4.
The case of two soft gluons is studied in Ref. [40]. We have checked that the ǫ-pole terms of the one-loop double-gluon current computed in Ref. [40] agree with the corresponding result that is obtained by using Eq. (20) and the colour conservation relation (10). To be precise about the absolute normalization of the one-loop current, we think that the factor e −ǫγ E has to be removed from the expansion parameter in Eq. (3.2) of Ref. [40].
We comment on the behaviour of the one-loop current J (1) (q 1 , . . . , q m ) with respect to the overall rescaling q k → λq k of all the momenta of the soft partons. To avoid the effects of branch-cut contributions from crossing different kinematical regions of soft and hard momenta, with limit ourselves to considering the case with λ > 0. According to Eq. (4) (see also the discussion below it) and Eq. (5), the tree-level current J (0) behaves as and the expected one-loop behaviour is The behaviour as in Eqs. (23) and (24) is indeed observed in the tree-level results of Eqs. (9) and (14) and in the one-loop soft single-parton current of Eq. (16). Using Eq. (23) and applying the λ rescaling to the explicit expression in the right-hand side of Eq. (20), we obtain the result . . , q m )+O(ǫ 0 ) (note that we neglect harmeless contributions proportional to the total colour charge of the hard partons). This result is perfectly consistent with Eq. (24), since Eq. (20) only embodies the correct ǫ-pole contributions to J (1) . In particular, in Eq. (20) these contributions are embodied in a 'minimal' form by systematically neglecting terms of O(ǫ n ) (n ≥ 0), with the sole exception of terms that arise from the ǫ-expansion of the overall factor c Γ ((4π) 2 c Γ = 1+O(ǫ)). The ǫpole contributions to J (1) can be expressed in alternative forms with respect to Eq. (20). In particular, the right-hand side of Eq. (20) can be supplemented with terms of O(ǫ n ) (n ≥ 0) in a manner that restores the behaviour in Eq. (24) to all orders in the ǫ expansion.
An alternative explicit form of the ǫ-pole contributions to J (1) for the soft multiparton (m ≥ 2) limit is as follows where the total soft momentum is denoted by q 1...m , and we have introduced the logarithmic functions ℓ ik and L ij of hard and soft momenta: It can be explicitly checked that the two expressions in Eqs. (20) and (25) only differ by terms of O(ǫ 0 ) and higher orders in ǫ while acting onto colour singlet quantities. We note that the logarithmic functions ℓ ik and L ij in Eqs. (27) and (28) are invariant under the overall rescaling q k → λq k (λ > 0) of the soft momenta. Therefore, the explicit expression in the right-hand side of Eq. (25) exactly fulfils the scaling behaviour in Eq. (24). We also note that both expressions of Eqs. (20) and (25) fulfil the colour flow conservation property in Eq. (11).
Throughout the paper we use the dimensional regularization procedure to deal with UV and IR divergences and, therefore, the momenta (and their associated phase space) of the virtual particles inside loops are analytically continued to d = 4 − 2ǫ space-time dimensions [58,59,60]. Different variants of dimensional regularization can be used, and each variant defines a specific regularization scheme (RS). The RSs that are mostly used are conventional dimensional regularization (CDR) [59,60], the 't Hooft-Veltman (HV) [58] scheme, dimensional reduction (DR) [61] and the four-dimensional helicity (4DH) scheme [62]. The momenta of the external-leg particles in the scattering amplitude can be either d-dimensional (CDR and DR schemes) or four-dimensional (HV and 4DH schemes). The number of spin polarization (helicity) states of the gluon also depends on the RS: externalleg gluons can have either d − 2 = 2 − 2ǫ polarizations (CDR) or 2 polarizations (HV, DR, 4DH), and virtual gluons can have either d − 2 = 2 − 2ǫ polarizations (CDR, HV) or 2 polarizations (DR, 4DH). Scattering amplitudes and, consequently, soft currents (as defined by the soft limit) depend on the RS. As for the RS dependence on external-leg particles, throughout the paper we formally express soft (tree-level and one-loop) currents in terms of external-leg momenta (p i , q k ) and corresponding polarization wave functions (ε(q k ), u(q k ), v(q k )): these expressions are formally RS invariant, although momenta and wave functions implicitly embody an RS dependence (which can be regarded as a dependence of O(ǫ)). At the one-loop level, soft currents (and scattering amplitudes) have a residual RS dependence that can be explicitly parametrized by the number of polarization states h g of virtual gluons. We write h g = 2 − 2ǫδ R and, therefore, we have (this is the same notation as used, e.g., in Refs. [21,20]) To formally express the explicit δ R dependence of the one-loop soft current J (1) we then define J (1) where both terms in the right-hand side are expressed through the same formal external-leg variables (momenta and wave functions), which embody an implicit dependence on the RS, and J (1) δ R =1 is obtained by setting δ R = 1 in the explicit expression of J (1) . Roughly speaking (e.g., modulo the implicit RS dependence due to the number of polarizations of the external partons), J RS in Eq. (30) represents the difference of J (1) between a given RS and the CDR (or HV) scheme.
One-loop scattering amplitudes have an explicit RS dependence on δ R . Considering the ǫ expansion up to including terms of O(ǫ 0 ), the dependence on δ R can be written in factorized form through the tree-level scattering amplitude and universal (process independent) coefficients [63]. Using these scattering amplitude results, we can obtain the explicit δ R dependence of the one-loop current up to the same order in the ǫ expansion. In particular, the δ R dependence of the one-loop scattering amplitude can be controlled through an ensuing δ R dependence of the one-loop expression V (1) [63,44,64] of the operator V in Eq. (19) and, therefore, we can explicitly compute the right-hand side of Eq. (30) where the coefficient γ k depends on the flavour of the soft parton k and it has an IR origin, while the coefficient b 0 has an UV origin. The explicit IR coefficients [63] and the UV coefficient [65] are Note that, analogously to the structure of its ǫ-pole contributions (see Eqs. (20) and (25)), the δ R dependence dependence of J (1) has a factorized structure in terms of its corresponding tree-level current J (0) . Since γ g = b 0 , in the case of a single soft gluon (m = 1), Eq. (31) agrees with the explicit result in Eq. (16). Incidentally, we recall [23] that the result in Eq. (16) (to all orders in the ǫ expansion) is valid in any RS, and thus the expression of the single soft gluon current in Eq. (16) is basically RS invariant (it does not depend on δ R , and the RS dependence is formally encoded in the corresponding RS dependence on ε(q) and the dimensionality of the external momenta p i , p j , q). As shown in Eq. (31), in the soft multiparton (m ≥ 2) case J  (25)) through the product J (0) · 1/ǫ 2 (J (0) conceptually embodies an RS dependence at O(ǫ) through its external-leg momenta and polarization vectors).
Throughout the paper we explicitly consider unrenormalized scattering amplitudes and currents. However, UV renormalization commutes with the soft limit and, therefore, the renormalization procedure can be straightforwardly applied to all the explicit expressions presented herein. In particular, since we are considering amplitudes and currents with (on-shell) massless hard partons, the renormalization procedure simply amounts to replace the bare coupling g S (or α S = g 2 S /(4π)) with its expression in terms of the renormalized running coupling g S (µ R ) (or α S (µ R )) at the renormalization scale µ R . In this respect, we recall that also the coupling renormalization is affected by RS subtleties. For instance, renormalizing the coupling by subtraction of the sole UV ǫ-poles (e.g., the term proportional to b 0 in Eq. (20)) in a given RS does not lead to an RS invariant definition of the renormalized coupling g S (µ R ): an additional finite renormalization shift of g S (µ R ) (whose size depends on the RS dependent coefficient δ R b 0 in Eq. (32)) [65] is necessary to achieve an RS independent definition of g S (µ R ).
As discussed and presented in this Section, the ǫ-pole contributions (and also the RS dependent contributions at O(ǫ 0 )) to the one-loop current J (1) for the general case of m (m ≥ 2) soft partons and an arbitrary number of hard partons are completely determined by Eqs. (20) or (25) (and Eq. (31)), and they are explicitly known as soon as the corresponding tree-level current J (0) is known. The determination of J (1) at O(ǫ 0 ) and, possibly, at higher orders in ǫ requires detailed one-loop computations and they have a high complexity. To have a rough idea of the computational complexity, we can simply observe that J (1) can (in principle) be determined by performing the soft limit of one-loop amplitudes according to Eq. (8). To apply Eq. (8) we have to consider amplitudes with m+n external legs, and the number n of non-soft external legs cannot be 'too small', otherwise the amplitude on the right-hand side of Eq. (8) vanishes. For example, the amplitude should have at least two external hard QCD partons (because of colour conservation) and one additional colourless external leg (because of momentum conservation): the soft limit of such an amplitude with m + 3 external legs leads to the current J (1) in the simplest case with two hard QCD partons. To get information on the colour-correlation structure of J (1) in the general case of several hard partons, the amplitude should have at least n = 4 hard QCD partons in its external legs (owing to colour-conservation relations, the cases with n = 2 and 3 hard partons lead to simplified colour structures; see, e.g., Sect. 5.3). In summary, the amplitudes to be considered should have m + n external legs with n ≥ 4: even in the case of double soft-parton radiation (m = 2), this implies (at least) six external legs. As is well known, one-loop computations of these multileg scattering amplitudes are definitely complex to be carried out in analytic form (which is necessary to perform the soft limit). The computation of the one-loop current J (1) can be highly simplified by using general methods (e.g., the method of Ref. [23])) that do not require a full direct computation of scattering amplitudes. However, despite some relevant simplification, even these methods have to deal with multileg one-loop Feynman integrals, whose evaluation is definitely complex, especially at high orders in the ǫ expansion.

Soft qq emission: the one-loop current
In this Section we present and discuss the results of our explicit computation of the QCD one-loop current for soft qq radiation. In Sect. 6 we also generalize the results to the cases of QED and mixed QCD×QED radiative corrections.
The tree-level current J (0) for emission of a soft-qq pair in a scattering amplitude with an arbitrary number of hard partons is given in Eq. (14). To evaluate the one-loop contribution J (1) we use the general (process-independent) method of Ref. [23] (the same method is used in the computations of Refs. [25] and [41]). The computational procedure is completely analogous to that in Ref. [23] (though it is extended from the case of a single soft gluon to the case of a soft-qq pair) and we do not repeat all the details. We have to evaluate a set of one-loop Feynman diagrams (as example, in Fig. 1 we show two contributing Feynman diagrams) in which the external-leg hard partons are coupled to virtual gluons by using the eikonal approximation (for both vertices and propagators), while the other vertices and propagators are computed by using the customary QCD Feynman rules. We perform the calculation by using both the Feynman gauge and the axial gauge n · A = 0, with an auxiliary light-like (n 2 = 0) gauge vector n µ . Combining all the contributing Feynman diagrams, the dependence on the gauge vector cancels at the integrand level (i.e., before performing the integration over the loop momentum) and the total axial-gauge integrand coincides with the Feynman gauge integrand: this provides us with an explicit check of the gauge invariance of the procedure and of the calculation.
As usual in the context of dimensional regularization, scaleless one-loop integrals vanish. Eventually we have to compute several (non-vanishing) tensor, vector and scalar one-loop Feynman integrals. Tensor and vector integrals are expressed in terms of scalar integrals by using customary techniques [66]. One-loop integrals with five external legs (pentagon integrals; see, e.g. the Feynman diagram in Fig. 1(a)) are expressed [67] in terms of one-loop integrals with four external legs (box integrals) plus remaining pentagon integrals in 6 − 2ǫ space-time dimensions, which only contribute at O(ǫ) (and higher orders in ǫ). We do not explicitly evaluate these contributions at O(ǫ). We eventually express the complete result in terms of a minimal set of basic one-loop scalar integrals. The set involves customary two-point and three-point (with at least one on-shell leg) Feynman integrals and some soft box integrals (box integrals with eikonal propagators). Part of these soft box integrals was already computed in Ref. [23] and the additional integrals are analogous to those encountered in Ref. [68]. We have performed an independent calculation of these soft box integrals and we find agreement with the results reported in Ref. [68] (see 'soft box 2' and 'soft box 4' in Sect. 4.2 of Ref. [68]). Our final result for the one-loop current J (1) is reported below. To present our results, we first define the tree-level and one-loop rescaled currentsĴ andĴ (1) as follows whereĴ (1) is written in terms of two components,Ĵ (1, div) andĴ (1, fin) . The rescaled current J (0) can be read from comparing Eqs. (14) and (33). The componentĴ (1, div) of Eq. (35) embodies the ǫ-pole contributions to J (1) , whileĴ (1, fin) includes all the remaining UV/IR finite contributions at O(ǫ 0 ) and higher orders in ǫ. The explicit expressions ofĴ (1, div) and where we have used the logarithmic functions of Eqs. (27) and (28) and we have introduced the shorthand notation ℓ ik (q 12 ) ≡ ℓ ik (with k = 1, 2) and L ij (q 12 ) ≡ L ij (i.e., we omit the explicit dependence on the argument q 12 ). The kinematical variable q 2 12⊥ ij that is used in Eq. (37) is We remark that the results in Eqs. (36) and (37) are valid in arbitrary kinematical regions, since the time component ('energy') of the outgoing momenta {q 1 , q 2 , p i , p j } of the soft and hard partons can have an arbitrary sign. According to the notation in Eq. (14) the colour indices α 1 and α 2 of the soft quark and antiquark are specified by considering α 1 , α 2 | J (1) (q 1 , q 2 ) ≡ J (1) α 1 ,α 2 (q 1 , q 2 ), and this leads to α 1 , α 2 | t a t b = t a t b α 1 α 2 in Eqs. (36) and (37).
The result in Eq. (36), which follows from our direct computation of J (1) , agrees with the ǫ-pole contributions that can straightforwardly be obtained by applying the general results in Eqs. (20) or (25) to the specific case of a soft qq pair (note that this agreement is valid modulo harmless terms that are proportional to the total colour charge i∈H T i of the hard partons). The expression in Eq. (36) has a term that is directly proportional to J (0) and additional terms that involve colour (and kinematical) correlations of the soft qq pair with two hard partons. These colour correlations are produced by the colour matrix factor t a t b T a i T b j . We remark that these correlations are not purely non-abelian, but they also include a component that is still present in the abelian limit of commuting colour matrices (this feature has to be contrasted with the one-loop single soft-gluon case of Eq. (16), in which correlations are purely non-abelian). In particular, this also implies that the soft current for lepton-antilepton radiation in massless QED has non-vanishing QED radiative corrections at one-loop level (see Sect. 6). The kinematical coefficients of these colour-correlation terms are proportional to the momentum function whose real part is the logarithm of a conformally invariant cross ratio, namely, We note that (analogously to the treatment in Sect. 3) in the computation ofĴ (1) (q 1 , q 2 ) we have dressed gluon propagators with one-loop vacuum polarization effects that are due only to massless partons. In particular, the terms proportional to N fĴ (0) in the right-hand side of Eqs. (36) and (37) are due to the vacuum polarization of N f massless quarks. Vacuum polarization effects of massive quarks can straightforwardly be included inĴ (1) (q 1 , q 2 ), and they produce corresponding (mass-dependent) contributions that are proportional tô We comment on the structure ofĴ (1, fin) . We have explicitly computed it up to O(ǫ 0 ) and the result is presented in Eq. (37). The expression ofĴ (1, fin) at O(ǫ 0 ) is quite compact and remarkably much simpler than expected. In particular, although it involves momentum functions of trascendentality equal to two, they are only powers of logarithmic functions with no additional dependence on dilogarithms Li 2 . Dilogarithms do appear in the computation of individual Feynman diagrams and loop integrals at O(ǫ 0 ), but they cancel in the complete result forĴ (1, fin) . The finite componentĴ (1, fin) includes a term that is proportional to J (0) and additional correlation terms with two hard partons whose colour structure is exactly analogous to that in Eq. (36) (and it embodies both abelian and non-abelian components). We have explicitly checked that no different colour-correlation structures occur at any higher orders in the ǫ expansion. The term that is proportional to J (0) explicitly depends on the RS parameter δ R : this dependence exactly agrees with that of the general result in Eq. (31).
We also comment on the kinematical dependence of the colour correlation terms. At the tree level the soft-qq current J (0) has a kinematical structure with a rational dependence on j(1, 2) · p i /p i · q 12 (see Eqs. (14) and (15)). In particular, this dependence leads to a collinear singularity if q 2 12 → 0 (i.e., if the momenta of the soft quark and antiquark are parallel). Exactly the same rational dependence (though possibly modified by logarithmic factors) occurs in the one-loop contributionsĴ (1, div) andĴ (1, fin) . However, by inspection of Eq. (37) we see that the one-loop interaction at O(ǫ 0 ) also produces a different type of kinematical dependence as given by the factor j(1, 2) · p i q 2 12 /(p i · q 12 q 2 12⊥ij ). This rational factor has no collinear singularity at q 2 12 → 0, but it potentially leads to a singularity in the limit q 2 12⊥ij → 0. This is a 'transverse-momentum singularity', since the kinematical variable q 2 12⊥ij in Eq. (38) is the transverse component of the momentum q 12 of the soft qq pair with respect to the momenta p i and p j of the colour-correlated hard partons in a reference frame in which p i and p j are back-to-back.
The transverse-momentum singularity in the current is partly screened by the logarithmic function L ij , and we have which shows that the current has a one-loop singularity of absorptive origin. Considering the physically most relevant kinematical region in which the soft quark and antiquark are produced in the final state (q 0 which shows that the transverse-momentum singularity is present in the scattering amplitude of a physical process in which the final-state soft qq pair is produced by the collision of the hard partons i and j in the initial state. We remark that this singularity has a pure quantum mechanics (loop) origin, and it occurs in the limit q 2 12⊥ij → 0 even if the transverse momenta q 1⊥ij and q 2⊥ij (q 2 k⊥ij = 2p i · q k p j · q k /p i · p j , k = 1, 2) of the soft quark and antiquark are separately large (i.e., they are separately non-vanishing) and q 2 12 is large. We also note that, setting q 12⊥ij = 0 at fixed non-vanishing values of q 2 12 and q 1⊥ij (or q 2⊥ij ), we have ℓ i1 − ℓ j2 = − ℓ j1 − ℓ i2 . Therefore, in the limit q 2 12⊥ij → 0 the factor ℓ i1 − ℓ j2 is (approximately) antisymmetric with respect to the exchange p i ↔ p j and this implies that we can perform the following replacement in Eq. (37): and it follows that, in the kinematical region with q 0 1 > 0 and q 0 2 > 0, the transversemomentum singularity has a purely non-abelian character (see the factor f abc in the righthand side of the relation (43)).
As we have just discussed, the singularity of the soft-qq current in the limit q 12⊥ij → 0 originates from one-loop interactions of the two soft partons. Therefore, we expect the presence of the transverse-momentum singularity also in the case of double soft-gluon emission at one-loop level. The one-loop double-gluon current computed in Ref. [40] indeed shows singular terms at q 12⊥ij → 0.
We have also computed the soft-qq one-loop current J (1) by explicitly evaluating its dependence on the RS parameter δ R to all orders in ǫ. Using the notation of Eq. (30) and considering the rescaled currents in Eqs. (33) and (35), we find the result Note that the δ R dependence at one-loop order is completely factorized with respect to J (0) . We also note that this factorized structure and the explicit expression of the ǫ-dependent factor in Eq. (44) are exactly equal to the corresponding RS dependence of the splitting function, Split(g → q(q 1 )q(q 2 )), of one-loop scattering amplitudes for radiation of a qq pair in the collinear limit [21,20,27].
The one-loop current J (1) for soft-qq emission has been independently computed in Ref. [40], and the corresponding result is presented in Sect. 3.3 therein. We first note that the one-loop result of Ref. [40] differs from our result already at the level of ǫ-pole contributions. However, we also note that we can remove such difference by adjusting the relative size of the four contributions in the right-hand side of Eq. (3.20) of Ref. [40]. More precisely, we modify the size of M s.l. by applying the replacement − 4 Nc + Nc 2 → −C F + Nc 2 = 1 2Nc to its colour coefficient (see the line 10 of Eq. (3.21)). We have contacted the author of Ref. [40] and he agreed with this correction. Performing such replacement, we have explicitly checked that the expression of J (1) in Eq.(3.20) of Ref. [40] agrees with our result (modulo the overall normalization of the one-loop current, which is not clearly specified in Ref. [40]) for both the ǫ-pole terms and the finite contributions at O(ǫ 0 ). However, we note that this check and comparison involve some 'limitations'. The explicit result of Ref. [40] only refers to the 'time-like' region, namely to the kinematical region in which the soft partons and all the hard partons are physically produced in the final state. Moreover, the result of Ref. [40] is specified for fixed (four dimensional) helicities of the soft quark and antiquark, and the comparison with our result requires the repeated use of the Schouten identity (which, precisely speaking, is valid only in four space-time dimensions) for the product of helicity spinors.

Soft qq radiation: squared amplitudes and current
Using the colour+spin space notation of Sect. 2, the squared amplitude |M| 2 (summed over the colours and spins of its external legs) is written as follows Accordingly, the square of the soft-emission factorization formula in Eq. (5) gives |M(q 1 , . . . , q m , p 1 , . . . , p n )| 2 ≃ M(p 1 , . . . , p n )| |J(q 1 , . . . , q m ) | 2 |M(p 1 , . . . , p n ) , (46) where, analogously to Eqs. (7) and (8), the symbol ≃ means that we have neglected contributions that are subdominant in the soft multiparton limit (i.e., the contributions that are denoted by the dots on the right-hand side of Eq. (5)). In the right-hand side of Eq.
The squared current |J | 2 is a colour operator that depends on the colour charges (and momenta) of the hard partons in M(p 1 , . . . , p n ). These colour charges produce colour correlations and, therefore, the right-hand side of Eq. (46) is not proportional to |M(p 1 , . . . , p n )| 2 in the case of a generic scattering amplitude † . As remarked on in Sect. 2, J is simply proportional to the unit operator in the spin subspace of the hard partons. Therefore, we note that the squared current |J | 2 of Eq. (47) still applies to spin-polarized hard-scattering processes, namely, to processes in which the spin polarizations of the hard partons are fixed (rather than summed over). Obviously, Eqs. (45)- (47) can also be properly generalized to the case in which the spin polarizations of one or more soft partons are fixed.
In the following part of this Section, we only consider soft-qq radiation and the corresponding soft current J (q 1 , q 2 ) (see Eq. (14) and Sect. 4). We define the loop expansion of the squared current as follows where |Ĵ| 2 (0ℓ) and |Ĵ| 2 (1ℓ) are the tree-level (0 loop) and one-loop rescaled contributions to |J | 2 , respectively.

The tree-level squared current
The tree-level squared current in Eq. (48) is whereĴ (0) is the rescaled current in Eqs. (14) and (33). The computation of the right-hand side of Eq. (49) is straightforward and the explicit result was first presented in Sect. 3.2 of Ref. [22]. We have where the momentum-dependent function I ij (q 1 , q 2 ) is (see Eq. (96) in Ref. [22]) Using colour charge conservation (see Eq. (10)), the tree-level squared current |Ĵ| 2 (0ℓ) can be recast in the following different form where the soft function w ij is w ij (q 1 , q 2 ) = I ii (q 1 , q 2 ) + I jj (q 1 , q 2 ) − 2 I ij (q 1 , q 2 ) .
The expressions in the right-hand side of Eqs. (50) and (52) are not identical at the algebraic level, but they are fully equivalent by acting onto scattering amplitudes (or, generically, colour-singlet states). The expression in Eq. (52) has a more straightforward physical interpretation, since the function w ij (q 1 , q 2 ) is directly related (see Sect. 5.3.1) to the intensity of soft-qq radiation from two hard partons, i and j, in a colour-singlet configuration.
The tree-level squared current in Eqs. (50) or (52) produces two-particle correlations between the hard partons. Their colour structure has the form of dipole contributions T i ·T j . We note that the momentum-dependent functions I ij (q 1 , q 2 ) and w ij (q 1 , q 2 ) are symmetric with respect to the exchange q 1 ↔ q 2 (they are also symmetric with respect to p i ↔ p j ). In contrast, our result for the one-loop squared current (see Sect. 5.2) produces both twoparticle and three-particle correlations and, moreover, it involves also an antisymmetric dependence on the momenta q 1 and q 2 .

The one-loop squared current
The one-loop squared current in Eq. (48) is where 'h.c' denotes the hermitian-conjugate contribution, and the rescaled currentsĴ The explicit computation of Eq. (54) produces some contributions that involve the fully-symmetric colour tensor d abc , with indices {a, b, c} in the adjoint representation of SU(N c ). The presence of d abc is a distinctive feature of (squared) currents for radiation of soft quarks and antiquarks.
Using d abc we also define the d-conjugated (quadratic) charge operator D i of the parton i as follows Performing the SU(N c ) colour algebra, we explicitly find i =q : where we have used Note that the tensor d abc is odd under charge conjugation. This fact is responsible for the opposite overall sign between the d-charge D i and the colour charge T i of quarks and antiquarks (see Eqs. (57) and (58)). Analogously, in the gluon case the d-charge b|D a i |c in Eq. (59) is symmetric with respect to b ↔ c, while the colour charge b|T a i |c is antisymmetric with respect to b ↔ c.
The explicit expression of the one-loop squared current |Ĵ(q 1 , q 2 )| 2 (1ℓ) is obtained by insertingĴ (0) (see Eqs. (14) and (33)) andĴ (1) (see Eqs. (35)- (37)) in the right-hand side of Eq. (54), and by performing the sum over the colours and spins of the soft quark and antiquark. We find the following result: ijk denotes the fact that these functions are symmetric under the exchange q 1 ↔ q 2 of the momenta of the soft quark and antiquark: Analogously, the superscript [A] in w ijk highlights the fact that these functions are antisymmetric under the exchange q 1 ↔ q 2 : ijk (q 2 , q 1 ) .
Therefore, w [A] and F [A] produce a quark-antiquark charge asymmetry in the one-loop squared current. We note that the charge-asymmetry contributions appear in Eq. (61) with the associated colour factors D i ·T j = d abc T a i T b i T c j and d abc T a i T b j T c k that have a linear dependence on the colour tensor d abc (which is odd under charge conjugation). The chargeasymmetry contributions to |J (q 1 , q 2 )| 2 have a quantum origin and are characteristic of the radiation of soft quark-antiquark pairs (the squared current |J(q 1 , . . . , q m )| 2 for radiation of m soft gluons is instead fully symmetric with respect to the soft-gluon momenta q 1 , . . . , q m ).
We present the explicit result of the ǫ expansion of the functions w [S] , w [A] , F [S] and F [A] up to O(ǫ 0 ). More precisely, we limit ourselves to presenting the expressions of these functions in the kinematical region where q 0 1 > 0 and q 0 2 > 0 (i.e., the soft quark and antiquark are produced in the physical final state), which is the most relevant physical region ‡ . In this region, the squared current depends (see Eqs. (36) and (37)) on the logarithms ℓ i1 ± ℓ j2 (which are purely real, independently of whether the momenta p i and p j are physically incoming or outgoing) and on the real part L ijR and discontinuity Θ (in) ij of the logarithm L ij . We have (see Eqs. (27) and (28)) where The function w ij has the following expression in the region where q 0 1 > 0 and q 0 2 > 0: This function (which is symmetric under the exchange i ↔ j) controls the size of the one-loop radiative corrections to the tree-level colour dipole correlations T i · T j .
We note that w ij also depends on colour coefficients, while F ijk only depends on parton momenta.

The function F [S]
ijk is associated with non-abelian three-particle correlations with colour charge factor f abc T a i T b j T c k . In the region where q 0 1 > 0 and q 0 2 > 0, we have the explicit result The charge-asymmetry contributions to Eq. (61) can be expressed through the function F ijk . In the region where q 0 1 > 0 and q 0 2 > 0 we have At arbitrary orders in the ǫ expansion, the two-particle correlation function w ijk as follows In contrast to w ij , we note that w ij is antisymmetric under the exchange i ↔ j of the hard-parton momenta. In particular, this antisymmetry of w [A] ij implies that in the sum over i and j of Eq. (61) we can replace D i · T j by its antisymmetric component, namely, ij has the following expression:

By inspection of Eqs. (66)-(70) we note that only the function F [S]
ijk exhibits a discontinuity with respect to the momenta of the hard partons (see Θ (in) ij in Eqs. (65) and (67)). The discontinuity contributes in the kinematical region where two hard-parton momenta i and j have negative time component (p 0 i < 0 and p 0 j < 0), namely, the partons i and j collide in the physical initial state. This discontinuity term of the squared current in Eq. (61) originates as interference between a one-loop absorptive (imaginary) contribution and the antihermitian colour factor if abc T a i T b j T c k (we recall that i, j and k refer to three distinct partons). Actually, the entire term proportional to f abc T a i T b j T c k in Eq. (61) has this origin § as absorptive/colour interference (the absorptive term being related to the kinematical region where q 0 1 > 0 and q 0 2 > 0).
As discussed in Sect. 4 (see Eqs. (41)- (43) and accompanying comments) the oneloop current of soft-qq emission has a transverse-momentum singularity at q 12⊥ij → 0. This singularity has a non-abelian character and an absorptive origin. At the level of the one-loop squared current, this singularity does appear in the function F  (70) we see the term (q 2 12⊥ij ) −1 L ijR → (q 2 12 ) −1 , which is not singular at q 12⊥ij → 0). Therefore, the transverse-momentum singularity at q 12⊥ij → 0 contributes through colour correlation f abc T a i T b j T c k to one-loop squared amplitudes for the class of processes with initial-state colliding partons i and j and two or more final-state hard partons (as recalled below in Eq. (71), the colour correlation vanishes if there is only one final-state hard parton). This class of processes includes, for instance, dijet (or heavy-quark pair) production in hadron-hadron collisions and the transverse-momentum singularity is directly related to the transverse momentum of the dijet system (heavy-quark pair). Interestingly, we note that this is the same class of processes that is sensistive to effects due to the violation of strict collinear factorization [26]. However, we remark on the fact that the transversemomentum singularity at q 12⊥ij → 0 and violation of strict collinear factorization are § The one-loop squared current for single soft-gluon radiation [23] has three-particle correlations of the type f abc T a i T b j T c k , which have an analogous origin as absorptive/colour interference. independent phenomena (e.g., the singularity at q 12⊥ij → 0 is not due to violation of strict factorization in the one-loop collinear limit of three partons, such as the soft quark and antiquark and a hard parton i or j).
Regarding three-particle correlations of the type f abc T a i T b j T c k with three distinct hard partons, we also recall two general features. As first noticed in Ref. [23], such colour correlations vanish by acting onto scattering amplitudes with only three hard partons (plus additional colourless external particles). Indeed, we have [23] where |i j k denotes a generic colour singlet state of three distinct hard partons i, j and k (the result in Eq. (71) simply follows from the colour conservation relation (T i + T j + T k )|i j k = 0). As pointed out in Refs. [69,70], such colour correlations vanish by considering their expectation value onto pure QCD amplitudes at the tree level. Namely, we have [69,71] where M (0) is a generic scattering amplitude with only quark and gluons external lines (and no additional colourless external particles) as obtained by tree-level QCD interactions. Therefore, the three-particle correlations f abc T a i T b j T c k contributes to the one-loop squared current in Eq. (61) only (see Eq. (71)) for processes with four or more hard partons and only (see Eq. (72)) through the introduction of either QCD loop corrections or electroweak interactions (see Refs. [23,70]) in the hard-parton scattering amplitude M(p 1 , . . . , p n ).
We present some general comments on the charge-asymmetry contributions to the oneloop squared current of Eq. (61). Such contributions produce non-vanishing effects only for specific classes of scattering amplitudes (see the discussion below and in Sect. 5.3) and quantities that are not invariant under charge conjugation. Obviously, due to their antisymmetry under the exchange q 1 ↔ q 2 , the charge-asymmetry contributions in Eq. (61) give vanishing effects after phase-space symmetric integration over the momenta q 1 and q 2 of the soft quark and antiquark. At the cross section level, the charge-asymmetry contributions can give non-vanishing effects to quantities in which the soft quark (or antiquark) is triggered, either directly (as it can be done for bottom or charm quark) or indirectly (e.g., through its fragmentation), in the final state. For instance, we recall that the Altarelli-Parisi splitting functions for collinear evolution of parton densities and fragmentation functions have a quark-antiquark charge asymmetry [33,72,73,74], which starts at O(α 3 S ) (the same perturbative order of the soft-qq one-loop squared current) and which does not vanish in the soft limit.
Considering M(p 1 , . . . , p n )| |J(q 1 , q 2 ) | 2 |M(p 1 , . . . , p n ) , the charge-asymmetry contributions vanish if M(p 1 , . . . , p n ) is a pure multigluon scattering amplitude, namely, if it has only gluon external lines (with no additional external qq pairs or colourless particles). This is a general consequence of the fact that the original (i.e., before performing the soft-qq limit) squared amplitude |M(q 1 , q 2 , p 1 , . . . , p n )| 2 is charge-conjugation invariant, since its external legs are gluons and a single qq pair (one cannot distinguish between the quark and the antiquark at the squared amplitude level). At the purely technical level, it turns out (as it can be verified) that the colour charge operators D i · T j and d abc T a i T b j T c k in Eq. (61) have vanishing expectation value onto pure multigluon amplitudes.
We also note that the three-particle correlations of the type d abc T a i T b j T c k in Eq. (61) contribute only for processes with four or more hard partons. Indeed, in the case of only three hard partons we have where the three distinct hard partons (i, j and k) in the colour singlet state |i j k are either three gluons or a gluon and a qq pair. The proof of Eq. (73) is given in Sect. 5.3.2 (see Eqs. (85), (89) and related comments).

Processes with two and three hard partons
The soft-emission factorization formula (46) for squared amplitudes embodies colour correlations produced by the squared current |J | 2 . In the case of scattering amplitudes with two or three hard partons (plus, necessarily, additional colourless external particles) the colour correlations have a simplified structure. A related discussion and some general results for multiple soft-gluon radiation can be found in Ref. [39]. The radiation of soft-qq pairs produces additional colour correlations from charge-asymmetry contributions: their main features are discussed in this Section.

Processes with two hard partons
We consider a generic scattering amplitude M BC (q 1 , q 2 , p B , p C ) whose external legs are two hard partons (denoted as B and C), a soft qq pair and additional colourless particles (which are never explicitly denoted). The two hard partons can be either a qq pair (note that we specify B = q and C =q) or two gluons ({BC} = {gg}). The corresponding scattering amplitude |M BC (p B , p C ) without the soft-qq pair is a colour singlet state. There is only one colour singlet configuration of the two hard partons, and the corresponding one-dimensional colour space is generated by a single colour state vector that we denote as |BC .
Since the soft-qq squared current |J (q 1 , q 2 )| 2 conserves the colour charge of the hard partons, the state |J| 2 |BC is also proportional to |BC . We write where |J | 2 BC is a c-number (it is the eigenvalue of the operator |J | 2 onto the colour state |BC ). Therefore, the soft-factorization formula (46) has the following factorized c-number form: with no residual correlation effects in colour space (the dependence on SU(N c ) colour coefficients is embodied in the c-number factors |J | 2 BC and |M BC | 2 ). In this respect, the structure of Eq. (75) is similar to that of soft-photon factorization formulae in QED. We recall that a c-number factorization formula analogous to Eq. (75) is equally valid for multiple soft-gluon radiation from two hard partons [39].
We note that Eqs. (74) and (75) are valid at arbitrary loop orders in the perturbative expansion of both the squared amplitude and the squared current. Therefore, by considering Eq. (74) and the loop expansion in Eq. (48), we can limit ourselves to evaluate the eigenvalues (c-numbers) |Ĵ| 2 (0ℓ)BC and |Ĵ| 2 (1ℓ)BC , which are the tree-level and one-loop contributions to |J | 2 BC .
The tree-level squared current |Ĵ| 2 (0ℓ) in Eqs. (48) and (52) depends on the colour dipole factor T B · T C and, by simply using charge conservation (T C |BC = −T B |BC ), we have . This leads to the tree-level result first presented in Ref. [22]: where w ij (q 1 , q 2 ) is given in Eq. (53) and C B is the quadratic Casimir coefficient of the hard parton (either The one-loop squared current in Eq. (61) depends on the colour dipole T B · T C (as at the tree level) and on charge-asymmetry colour correlations.
By using Eqs. (56)-(59), we have the following colour algebra results: i =q : We note that the operators Combining all the contributions in Eq. (61), we find the following final result: where the functions w ij are given in Eqs. (66) and (70), respectively In the case of soft-qq radiation from the hard partons {BC} = {qq} (see Eq. (80)), we do find charge-asymmetry contributions in |Ĵ (q 1 , q 2 )| 2 (1ℓ)BC . We recall that the function w ij (q 1 , q 2 ) is antisymmetric with respect to the separate exchanges q 1 ↔ q 2 and i ↔ j. Therefore, in Eq. (80) the asymmetry in the momenta of the soft-qq pair is correlated with a corresponding asymmetry in the momenta p B and p C of the hard q andq. In particular, |Ĵ (q 1 , q 2 )| 2 (1ℓ)BC is invariant under the overall exchange of fermions and antifermions (i.e. {q 1 , p B } ↔ {q 2 , p C }), consistently with charge-conjugation invariance.
In the case of soft-qq radiation from two hard gluons, the one-loop result in Eq. (81) shows no charge-asymmetry effects. We state that this feature persists at arbitrary orders in the QCD loop expansion. The absence of charge-asymmetry effects follows from the fact that the c-number squared current |J(q 1 , q 2 )| 2 BC for {BC} = {gg} is entirely controlled by QCD interactions, with absolutely no dependence (both explicitly and implicitly) on the production mechanism of the two hard gluons. Therefore, such squared current is chargeconjugation invariant (similarly to the squared amplitude for the process gg → qq) and one cannot distinguish between the soft quark and antiquark.

Processes with three hard partons
Before considering the explicit evaluation of the soft-qq squared current |J(q 1 , q 2 )| 2 for processes with three hard partons, we recall and derive some general algebraic relations for the action of colour charge correlations operators onto a generic colour singlet state |ijk formed by three distinct partons i, j and k (i = j, j = k, k = i) in arbitrary representations of the gauge group SU(N c ). We consider the correlations operators that appear in |J (q 1 , q 2 )| 2 up to one-loop level, namely, As is well known, the action of dipole factors onto |ijk can be evaluated in terms of quadratic Casimir coefficients T 2 i = C i (see the Appendix A of Ref. [16]). We have and related permutations of i, j, k. In particular, any generic colour singlet state |ijk is an eigenstate of T i · T j or, equivalently, the action of T i · T j onto |ijk is always proportional to the unit operator in colour space. The result in Eq. (82) simply follows from the charge conservation relation (T i + T j + T k )|ijk = 0, which also leads to the result in Eq. (71) for the operator f abc T a i T b j T c k .
Considering charge-asymmetry correlations and using colour conservation ( T a k |ijk = −(T a i + T a j )|ijk ), we have the following relations and related permutations of i, j, k. We note that we are dealing with seven colour correlations operators (six two-particle correlations of the type D i · T j , and the three-particle correlation d abc T a i T b j T c k ) whose action onto |ijk is 'non-trivial', while the action of the three operators D i · T i is directly worked out in c-number form in terms of cubic Casimir coefficients (see Eqs. (77)-(79)). Since colour conservation leads to the six linear relations (exploiting permutations) in Eqs. (83) and (84), all the non-trivial colour correlations can be expressed in terms of a single correlation operators. To explicitly show this, we derive the following relations. The three-particle correlation is directly related to cubic Casimir coefficients as follows The three symmetric (with respect to i ↔ j) two-particle correlations are equal and directly related to cubic Casimir coefficients as follows The three antisymmetric (with respect to i ↔ j) two-particle correlations fulfil two independent linear relations (which are related through i ↔ j) as follows The derivation of Eqs. (85)-(88) from Eqs. (83) and (84) is relatively straightforward. For instance, Eq. (85) is derived by first summing Eq. (84) and its two independent permutations to obtain 3d abc T a In the specific cases in which i, j, k are either three gluons or a gluon and a qq pair, we can use the explicit results for D i · T i in Eqs. (77)-(79) and, consequently, Eq. (85) gives d abc T a i T b j T c k |ijk = 0 (this proves Eq. (73)) and from Eq. (86) we obtain and related permutations of i, j, k.
Regarding the vanishing value of the correlations f abc T a i T b j T c k and d abc T a i T b j T c k in Eqs. (71) and (73) We summarize our general discussion on colour correlations for processes with three hard partons i, j, k in arbitrary colour representations of SU(N c ). The charge-symmetric component of |J (q 1 , q 2 )| 2 up to one-loop order is proportional to the unit operator in colour space, and it can be expressed in c-number form, in terms of quadratic Casimir coefficients (see Eqs. (71) and (82)). The charge-asymmetry component of |J (q 1 , q 2 )| 2 at one-loop order can eventually be expressed (see Eqs. (85)-(88)) in terms of cubic Casimir coefficients (cnumbers) and a single operator (e.g., D i · T j ) whose action onto the colour singlet state |ijk has to be explicitly computed (the result depends on the specific state |ijk ).
We come to explicitly discuss soft-qq radiation from scattering amplitudes with three hard partons in the specific cases that are relevant within perturbative QCD. We consider a generic scattering amplitude M ABC (q 1 , q 2 , p A , p B , p C ) whose external legs are colourless particles (which are not explicitly denoted), a soft-qq pair and three hard partons (denoted as A, B, C) that can be either a gluon and a qq pair ({ABC} = {gqq}) or three gluons ({ABC} = {ggg}). The corresponding scattering amplitude |M ABC (p A , p B , p C ) without the soft-qq pair is a colour singlet state formed by the three hard partons A, B and C. We consider the cases {ABC} = {gqq} and {ABC} = {ggg} in turn.

g qq case
We specifically set A = g, B = q and C =q.
There is only one colour singlet configuration of the three hard partons, gqq, and the corresponding one-dimensional colour space is generated by a single colour state vector that we denote as |ABC . Therefore, we are in a situation in which we can apply the same reasoning of Sect. 5.3.1 (see Eqs. (74) and (75) and the accompanying discussion). The state |ABC is an eigenstate of the soft-qq squared current |J (q 1 , q 2 )| 2 , and the soft-factorization formula (46) has the following factorized form: The tree-level squared current |Ĵ | 2 (0ℓ) in Eqs. (48) and (52) involves colour dipole correlations. Using Eq. (82), dipole correlations can be expressed in terms of quadratic Casimir coefficients and this leads to the tree-level result first presented in Ref. [22]: where w ij (q 1 , q 2 ) is given in Eq. (53). Note that the result in Eq. (92) is symmetric under the exchange p B ↔ p C of the momenta of the hard quark and antiquark.
The one-loop squared current |Ĵ(q 1 , q 2 )| 2 (1ℓ) in Eq. (61) has contributions with and without charge asymmetry. Owing to Eq. (71), the charge-symmetric contributions only involve colour dipole correlations, as at the tree level. As discussed in Eqs. (85)-(88), the chargeasymmetry contributions require the explicit evaluation of a single correlation operator of the type D i · T j . We consider the operator D B · T C , whose action onto |ABC can be related to the action of the dipole operator T B · T C . Indeed, we have  88)). We find the following result for the eigenvalue |Ĵ (q 1 , q 2 )| 2 (1ℓ)ABC of the one-loop squared current for soft-qq radiation: BC (q 1 , q 2 ) where the functions w ij (q 1 , q 2 ) and w ij (q 1 , q 2 ) are given in Eqs. (66) and (70), respectively. We note that the charge symmetric contribution in Eq. (94) is symmetric under the exchange p B ↔ p C of the hard quark and antiquark. The charge-asymmetry contribution in Eq. (94) is instead antisymmetric under the exchange p B ↔ p C , in complete analogy with the corresponding contribution for soft-qq radiation from two hard partons (see Eq. (80)).

g g g case
We now consider the case in which the three hard partons A, B and C are gluons. The colour singlet space spanned by the three hard gluons is two-dimensional. It is convenient to choose the basis formed by the orthogonal colour state vectors |(ABC) f and |(ABC) d that are defined as follows where a, b, c are the colour indices of the three gluons. We note that the two states in Eq. (95) have different charge conjugation. The scattering amplitude |M ABC (p A , p B , p C ) is, in general, a linear combination of the colour antisymmetric state |(ABC) f and the colour symmetric state |(ABC) d , and we write As examples of the scattering amplitude |M ABC (p A , p B , p C ) , we can mention the three scattering processes H → ggg, γ → ggg and Z → ggg. In the Higgs boson process H → ggg (see, e.g., Ref. [75]) the amplitude component M d of Eq. (96) vanishes, while in the photon process γ → ggg (see, e.g., Ref. [76]) we have M f = 0. In the case of the Z boson process Z → ggg both components M f and M d are not vanishing (see, e.g., Ref. [76]).
We also note that all these scattering amplitudes are produced through QCD interactions involving quark loops (within the Standard Model, gluons have tree-level interactions only with quarks and, consequently, M ABC vanishes at the tree level).
We have previously discussed the case of the three hard partons {ABC} = {gqq}, which generate a one-dimensional colour singlet space. The fact that the colour singlet space is two-dimensional for {ABC} = {ggg} is an essential difference. In particular, in the case {ABC} = {ggg} the all-order soft-factorization formula (46) for squared amplitudes cannot be recast in the factorized c-number form of Eq. (91). The action of the squared current |J | 2 onto |M ABC of Eq. (96) is colour conserving, but it can produce colour correlations between the two colour singlet states |(ABC) f and |(ABC) d of the three hard gluons. In general, the squared soft current |J | 2 can be represented as a 2 × 2 correlation matrix that acts onto the two-dimensional space generated by |(ABC) f and |(ABC) d . The all-order structure of this correlation matrix is discussed in Ref. [39] for the case of multiple softgluon radiation. In the following we explicitly consider soft-qq radiation at the tree level and one-loop order.
The tree-level squared current |Ĵ | 2 (0ℓ) in Eqs. (48) and (52) only involves colour dipole correlations, whose action onto both |(ABC) f and |(ABC) d is proportional to the unit matrix in colour space (see Eq. (82)). Therefore, the contribution of |Ĵ | 2 (0ℓ) to the factorization formula (46) can be expressed in factorized c-number form and, using Eq. (82), we have (see also Ref. [22]) where w {ABC} (q 1 , q 2 ) = w AB (q 1 , q 2 ) + w BC (q 1 , q 2 ) + w CA (q 1 , q 2 ) , and w ij (q 1 , q 2 ) is given in Eq. (53). Since w ij is symmetric under the exchange p i ↔ p j , we note that the function w {ABC} has a completely symmetric dependence on the gluon momenta p A , p B , p C (as required by Bose symmetry). We also note that Eq. (97) is valid at arbitrary orders in the loop expansion of the amplitude M ABC (p A , p B , p C ).
The action of the one-loop squared current |Ĵ (q 1 , q 2 )| 2 (1ℓ) in Eqs. (48) and (61) onto |M ABC involves charge symmetric and charge-asymmetry contributions. As summarized in the discussion below Eq. (89), the charge symmetric contributions are proportional to the unit matrix in colour space, while the charge-asymmetry contributions can be expressed in terms of a single colour correlation operator. Specifically, by using Eq. (82)) and Eqs. (85)-(88), we explicitly find [ABC] (q 1 , q 2 ) |ABC , (99) where w and the functions w ij (q 1 , q 2 ) and w ij (q 1 , q 2 ) are given in Eqs. (66) and (70), respectively. We note that the charge symmetric contribution to Eq. (99) depends on the function w [S] {ABC} that has a fully symmetric dependence on the hard-gluon momenta p A , p B , p C . The charge-asymmetry function w [ABC] is instead antisymmetric under the exchange of two gluon momenta (e.g., p A ↔ p B ).
The charge-asymmetry operator D B · T A in the right-hand side of Eq. (99) acts differently onto the two colour states |(ABC) f and |(ABC) d of Eq. (96). By explicitly performing the SU(N c ) colour algebra, we find the following result: and we note that the operator D B · T A produces 'pure' transitions between the colour symmetric and colour antisymmetric states |(ABC) f and |(ABC) d , which have different charge conjugation.
Using Eqs. (96), (99) and (102), we obtain the final result for the contribution of the one-loop soft-qq squared current to squared amplitudes with three hard gluons. We find [ABC] (q 1 , q 2 ) which is not simply proportional to |M ABC | 2 (unlike the corresponding result in Eq. (91) for {ABC} = {gqq})). In contrast with the case of scattering amplitudes with two hard gluons (see Eq. (81)), we note that the expression in Eq. (103) involves a charge-asymmetry contribution that is not vanishing, provided the hard-scattering amplitude includes nonvanishing components M f and M d (i.e., M ABC has no definite charge conjugation). Such feature of M ABC depends on the specific production mechanism of the three hard gluons. The functions w [ABC] and (M † d M f + h.c.) are separately antisymmetric under the exchange of two gluon momenta and, consequently, their product is symmetric. Therefore, the righthand side of Eq. (103) (including its charge-asymmetry contribution) is fully symmetric under permutations of the three hard gluons, as expected and required by Bose symmetry.

Soft fermion-antifermion radiation in QED and mixed QCD×QED
Our results in Sects. 4 and 5 for soft-qq emission can be generalized to consider the emission of a soft fermion-antifermion (ff ) pair through QED (photon) interactions and mixed QCD×QED (gluon and photon) interactions. Before presenting the results, we precisely specify our framework.
The soft fermions can be either massless quarks (f = q) or electrically-charged massless leptons (f = ℓ). We consider generic scattering amplitudes, M, whose external particles are massless quarks and gluons, massless leptons and, additionally, particles that carry no colour charge and no electric charge (i.e., photons, Higgs and Z bosons in the context of Standard Model). The external particles (i.e., their momenta and quantum numbers) of M are treated as outgoing particles (as already specified in Sect. 2 for the pure QCD case). The internal legs of M can include massless (photons, gluons) and massive (e.g., heavy quarks and/or W ± bosons) particles. If an external ff pair becomes soft, the scattering amplitude M is singular and the singular behaviour is due to the production of the soft-ff pair through QCD (gluon) and QED (photon) interactions. We formally treat QCD, QED and mixed QCD×QED interactions on equal footing. Therefore, the scattering amplitude M has a generalized perturbative (loop) expansion in powers of two unrenormalized couplings: the QCD coupling g S and the QED coupling g (g 2 /(4π) = α is the fine structure constant at the unrenormalized level). Regarding the RS of the UV and IR divergences, photons and charged leptons are treated in the same way (see Sect. 3) as gluons and massless quarks, respectively.
The factor ∆ f in the right-hand side of Eq. (106) is a colour operator that depends on the type of soft fermion f . If f = ℓ, we simply have ∆ f = 1. If f = q, ∆ f is the projection operator onto the colour singlet state of the ff pair, namely, by using the colour space notation of Sect. 2 we have α 1 , α 2 | ∆ f = δ α 1 α 2 .
The rescaled currentĴ (1) (2γ) in Eq. (107) is entirely due to QED interactions, and it has the following expression: whereĴ (0) (1γ) is given in Eq. (106). In Eq. (108) the soft fermion f can be either a quark or a lepton and, similarly, the charged hard particles i, j ∈ H can include quarks and leptons. The factor N ch. in the right-hand side of Eq. (108) is analogous to the factor T R N f of the QCD expressions in the right-hand side of Eqs. (36) and (37). The coefficient N ch. depends on the squared electric charges of the massless ‡ quarks and leptons in the theory, and we have The one-loop termĴ (1) (1γ) in Eq. (107) is due to mixed QCD×QED interactions. It has ‡ Analogously to the QCD case, in the one-loop currentĴ (1) (2γ) (q 1 , q 2 ) we have not considered and included vacuum polarization effects due to massive particles (charged leptons, quarks and W ± ).
The one-loop term in the squared current of Eq. (111) includes all possible contributions that are proportional to the powers (g 2 S ) 3−n (g 2 ) n with 0 ≤ n ≤ 3. We write it in the following form: two-particle and three-particle charge-asymmetry correlations (i.e., the summed index k can also be equal to either i or j). In the QCD case, we also noticed that three-particle correlations do not contribute to the squared of the soft-ff current for emission from three hard partons (see Eqs. (71) and (73)). A corresponding observation does not apply to the one-loop contribution in Eq. (115). For example, we can consider soft-ff emission from the hard-scattering processūd → W − →ν e e − (the charges of the outgoing hard particles are {+2/3, +1/3, −1}) and we see that the product e k e i e j of three distinct charges in Eq. (115) does not vanish.
where the charge symmetric function w ij (q 1 , q 2 ) and the charge-asymmetry function F ijk (q 1 , q 2 ) are given in Eqs. (53) and (68), respectively. We note that the one-loop term in Eq. (117) is not vanishing only if the soft fermion is a quark. Similarly to Eq. (115), the summed index k in Eq. (117) can also be equal to either i or j.
By direct inspection of the tree-level and one-loop results in Eqs. (113) , (115) , (117) and (118), we see that the charge symmetric (charge-asymmetry) contributions are proportional to even (odd) powers of e f , as expected from charge conjugation symmetry.

Summary
We have considered the radiation of two or more soft partons in QCD hard scattering. In this soft limit the scattering amplitude is singular, and the singular behaviour is controlled in factorized form by a multiparton soft current, which has a process-independent structure. At loop level, the scattering amplitudes and the soft current have UV and IR divergences, which we regularize in the form of ǫ poles by analytic continuation in d = 4 − 2ǫ space-time dimensions.
We have discussed the general structure of the ǫ-pole divergences of the multiparton soft current. We have considered the soft current at one-loop order and we have presented the explicit form of its ǫ-pole (divergent) contributions. We have also discussed the RS dependence of the one-loop soft current.
In the remaining part of the paper we have considered the specific case of soft qq radiation, by presenting a detailed study at one-loop order. Considering arbitrary kinematical regions of the soft-parton and hard-parton momenta, we have explicitly computed the oneloop current by including the finite terms at O(ǫ 0 ). We find a relatively simple expression, which, for instance, includes powers of logarithmic functions but no dilog functions.
We find that the one-loop current produces a new type of singularity if the soft-qq pair is radiated with a vanishing transverse momentum with respect to the direction of two colliding hard partons in the initial state. This new transverse-momentum singularity has a quantum (more precisely, absorptive) origin and a purely non-abelian character. Owing to its dynamical origin, the transverse-momentum singularity can appear also in the one-loop current for double soft-gluon emission.
We have computed the one-loop contribution of the squared current for soft-qq emission and the ensuing colour correlations for squared amplitudes of generic multiparton hardscattering processes. We have also explicitly considered the specific cases of processes with two or three hard partons, in which the colour correlation structure can be partly simplified.
We find that, despite its absorptive origin, the new one-loop transverse-momentum singularity contributes to squared amplitudes (and, hence, cross sections) of scattering processes with two initial-state colliding partons (hadrons) and two or more hard partons (jets) in the final state.
At variance with the case of multiple soft-gluon radiation, the emission of soft fermions and antifermions lead to charge asymmetry effects. We have discussed in details the charge asymmetry contributions of the one-loop squared current for soft qq radiation.
We have finally generalized our QCD study of soft qq emission to the study of QED and mixed QCD×QED radiative corrections in the context of soft fermion-antifermion radiation. We have presented the corresponding one-loop results for the soft current and its square. agreement number 754496 and by the Generalitat Valenciana (Spain) through the plan GenT program (CIDEGENT/2020/011).