Asymptotic Charges and Coherent States in QCD

We study the connection between asymptotic symmetries in non-Abelian gauge theories and the generalised coherent states following from the application to QCD of the Faddeev-Kulish approach to asymptotic dynamics. We compute the large gauge transformation properties of the soft evolution operators and use this to define the quantum corrected, non-linear contribution to the asymptotic charges. We then compute the one-loop, leading IR-divergent, correction to matrix elements of the charges inserted between dressed scattering states and show that the results depend on a particular order of soft limits. For one choice of ordering we find that the conservation law for the asymptotic charges is not corrected, while for a second we find a correction proportional to the one-loop soft current.


Introduction
The study of the infrared behaviour of gauge theory scattering amplitudes has a long history and in the case of QED has essentially been understood since the work of Bloch and Nordsieck [1], though it has been refined over the years [2][3][4][5]. The standard approach involves the computation of amplitudes which are formally singular -they in fact vanish after exponentiation of the perturbative divergences. One then focuses on inclusive quantities involving arbitrary numbers of real soft-photons which cancel the IR-divergences from virtual photons in loops. An alternative approach is to directly formulate infrared-finite Smatrix elements by choosing appropriate asymptotic states. For QED this approach where the asymptotic states are not eigenstates of the photon number operator but rather have the form of coherent states, was, starting from the work of Chung [6] 1 , developed by a number of authors, e.g. [8][9][10][11]. Faddeev and Kulish, building on the work of Dollard for the Coulomb problem in non-relativistic quantum mechanics [12], related the structure of these coherent states to the form of the large-time Hamiltonian [13,94]. The approach was partially extended to the much more complicated case of non-Abelian gauge theory [14][15][16][17][18], [19][20][21][22][23][24] and more recently to perturbative gravity [25].
The observation that asymptotic conservation laws, which follow from Noether's second theorem [26] for large gauge transformations, pave the way to understanding the infrared dynamics of gauge theories [27,28] (see [29] for a review and more complete references), has lead to a renewed interest in the study of coherent state operators and soft-dressing more generally [30][31][32][33][34][35][36][37][38]. It has already been demonstrated that for QED [32] and perturbative gravity [34] the coherent states relevant to the construction of an infrared finite S-matrix follow from the symmetry of asymptotic charges. The existence of an infinite dimensional symmetry group has lead to the interpretation of the QED vacuum as being degenerate and that scattering processes are accompanied by a shift in the vacua. Infrared divergences due to massless particles which result in the vanishing of S-matrix elements are thus connected with the "wrong" choice of the in-and out-vacua, and a cure can be sought in a systematic way by considering the asymptotic charges. An analogous statement can be made for perturbative gravity, using BMS supertranslation charges to find suitable asymptotic states for an infrared-finite gravity S-matrix.
The corresponding understanding of QCD infrared dynamics is significantly less complete. Compared to QED, the natural complication that arises is that gauge bosons selfinteract in a non-trivial way. The persistence of these non-trivial self-interactions at early and late times is central to the behaviour of the infrared properties of QCD [75]. It is our aim here to understand the connection between asymptotic symmetries in non-Abelian gauge theories and generalized IR leading coherent states. We will follow the approach of Catani, Ciafaloni and Marchesini, [19][20][21][22][23] (see [39] for an introduction) which uses energy ordering in each interaction to systematically organise the divergences due to soft gluons. Non-Abelian gauge theories of course also have collinear divergences which have been treated in the coherent operator approach [40], see also [41] and the recent interesting work [42,43]. We will focus on the leading IR-divergence, though the method can be extended in principle to include sub-leading divergences. To be slightly more precise, when we extend the tree-level considerations to one-loop we will focus on the IR double pole terms arising from simultaneously soft and collinear regions of parameter space.
In the following a central goal is the definition of asymptotic charges at loop level. As we discuss below we start from the linear part of the charge that arises in the classical construction and which can be understood to act on Fock, or bare, asymptotic states at null-infinity. However, as the S-matrix is ill-defined for such states we must introduce the soft-evolution operators generated by the QCD Hamiltonian which is used to define dressed states. We then consider the transformation properties of these soft-evolution operators under transformations generated by the asymptotic charges and we show how the softevolution operators can be viewed as generating the non-linear part of the charge from the linear part. This implies the existence of an infinite sequence of quantum corrections to the non-linear charge. This charge can then be used to define a Ward identity for scattering amplitudes. The main result of this paper is the computation, at leading IR divergence, of the one-loop corrections to the Ward identity of the asymptotic charge which is related to the soft gluon theorem. We demonstrate explicitly that the contribution from the soft gluon dressing factor cancels the contribution from the one-loop soft factor in QCD ( [44,45]), showing that the asymptotic charge produces a vacuum that is orthogonal to all scattering states built on the original vacuum as in [32]. We will see that this result depends on the precise prescription for the order of soft limits in the definition of the charge and dressing factor which is analogous to the order of limits discussed in [46,47]. There are two order-of-limit choices we must make, for the first, in the definition of the charge we take a order of limits used in [46], which corresponds to the standard choice in dimensionallyregularized soft limits. This choice is motivated, as in [46], by its use in problems such as the computation of physical cross-sections. For the second choice, in the definition of the dressed states, we compute the result using both orderings for comparison. As we discuss in Sec. 6 it is however natural to choose the prescription which preserves, where possible, the symmetries of the theory. We will see that there is indeed a particular prescription that preserves the asymptotic symmetries to the order we compute, as we find that with this choice the Ward identity receives no corrections at one-loop and leading IR-divergence.
Preliminaries Perturbative computations are relevant in QCD as they are related to experimental observables due to two important properties: asymptotic freedom and factorisation. Confinement naturally sets a scale Λ QCD such that for partons to be well-defined objects we require the existence of a lower cutoff on the lowest momentum transfer in a given process. Perturbative quantities are then related to physical cross-sections by convolution with non-pertubative but universal objects, e.g. parton distribution functions. In this work we will initially consider asymptotic states formed by acting on the Fock vacuum with creation/annihilation operators for the hard incoming or outgoing partons (these can be either massless gluons or massive quarks) which are labelled by momenta p , colour indices α (corresponding to the fundamental representation for quarks and the adjoint for gluons) and helicity indices σ as appropriate.
We will consider n-particle IR divergent S-matrix elements between such asymptotic states which give the usual perturbative scattering amplitudes. The behaviour of such amplitudes as individual gluons become soft and the relation to the asymptotic charge plays a key role in our work. We will properly introduce and define the linearized asymptotic charge, Q lin , in subsequent sections however for now it is sufficient to state that it involves soft gluon creation/annihilation operators a a † σ (ω q )/a a σ (ω q ) with vanishing gluon energy ω q 0 and has matrix elements of the form with a careful interpretation of the limit. At tree level it is known that the limit in the definition of the charge can be understood as These terms can be computed using the tree-level soft-gluon theorem with the soft-current given by where the soft gluon of momentum q, colour a and helicity σ is taken to be outgoing. In this expression, and similar expressions below, as the limit does not strictly exist the notation lim ωq→0 should be understood as referring to the leading term in an expansion in small ω q . If the gluon was incoming, there would be an overall minus sign and the corresponding polarisation vector would beε σ . At loop level the issue of the soft limit is more subtle: after computing the perturbative terms defining the matrix elements one may attempt to take the soft limit at the level of the integrands before performing loop integrations or alternatively one may keep ω q finite and take the limit only after performing the loop integrations. It is known, from the case of subleading IR behaviour of graviton amplitudes [46,47,93], that the order has important consequences for the interpretation of quantum corrections to the Ward identities for asymptotic symmetries. In this work for the definition of matrix elements of the asymptotic charge we always take the soft limit after the computation of matrix elements. That is we will define a regularised charge Q lin (ω q ) and then define out|Q lin (0)S|in := lim ωq→0 out|Q lin (ω q )S|in (1.7) so that for the insertion of the soft gluon operator we use the usual soft gluon theorem. The regularisation of the charge will be related to that of the soft-limits of amplitudes. When we discuss the soft-evolution operators we will follow [22] and use an energy cut-off. However for explicit computations we will make use of the one-loop soft-limits which were computed using dimensional regularisation, with parameterˆ = (d − 4)/2, in [45] (see also [44, 82-85, 95, 96] for an earlier work in the case where the quarks are massless) and can be written as where the first term on the right-hand side is the iterated tree result involving the tree-level soft current (1.6) while the second term is due to the one-loop soft current which is, to leading divergence, where C A is the adjoint quadratic Casimir. It is important to note that this leading double pole is due to both collinear and soft divergences at leading logarithmic accuracy and as we will see the coherent state will deal with both of them at this order, as explained in [20,23]. However before computing such matrix elements we first understand in more detail the definition of the asymptotic charge and its expression in terms of free-field operators. In order to do this we must review both the classical construction of the charges, the softevolution operators in QCD and the transformation properties of the evolution operators under large gauge transformations.

Asymptotic Charges for QCD
We will be interested in the asymptotic charges of Yang-Mills theory related to large gauge transformations, which are those that are non-vanishing on the boundary of space-time. Quite generally, Noether's second theorem relates the existence of a local symmetry to a two-form κ νµ which can be integrated over a codimension-two sphere to define a charge. That is, the local symmetry implies that there exists a conserved current where the function (x) parametrizes local symmetry transformations and the current S µ ( ) vanishes on-shell If is constant, one recovers the usual conserved current that couples to the gauge field. If we consider the integral of j µ ( ) over a manifold Σ this charge is non-vanishing only if the function has support at the manifold boundary ∂Σ, otherwise it is trivially zero. This is equivalent to the statement that gauge symmetries are not really "symmetries" in the proper sense. Let us review the construction of classical non-Abelian asymptotic charges at I ± for QCD in four dimensional Minkowski space using the usual conformal compactification (for details see [27,29] and see Appendix A for the notations used in this paper). The two-form, κ, is given by the field strength associated to the gauge field A µ = A a µ T a Figure 1. Penrose diagram picture of Minkowski compactification which obeys the equations of motion where j M ν is the matter colour current. The classical charge for non-Abelian large gauge transformations with parameter is then given as, see for example [48,49] and also [50]. The relevant surfaces for our asymptotic charges are future and past null infinity, usually denoted I + and I − (see Fig. 1). The natural coordinates for discussing I + are the so-called retarded Bondi coordinates such that the flat-space metric becomes (2.8) The boundary at r → ∞, I + , is thus parametrized by the coordinates (u, z,z). The advanced coordinates, (r, v = t + r, z,z), are most convenient for I − and where necessary we will apply the antipodal matching conditions as in [28]. Furthermore, as it is convenient for making contact with perturbative computations, relevant for the construction of the coherent states, we will focus on Lorenz gauge in this paper with the gauge fixing condition The same charge can be found in non-covariant gauges such as retarded radial or radiation gauge. In general the asymptotic symmetries may depend on the choice of gauge fixing, which can also be seen at the path integral level (see [49]), but for the leading asymptotic charges one finds the same result for both covariant and physical gauge choices. Finally we will be interested in large gauge transformations parametrised by a function = a (z,z)T a which labels the transformations at infinity, for example at I + , In retarded Bondi coordinates on the celestial sphere S 2 defined by the u → −∞ region of I + , i.e. I + − , the charge (2.6) becomes where F ru = lim r→∞ r 2 F ru . This charge is usually denoted Q + to distinguish it from the corresponding charge Q − on I − . In order to avoid a proliferation of superscripts we will mostly neglect this index and hopefully it is clear from context to which charge we are referring.
The charge can now be rewritten, assuming the field strengths vanish at I + + (for details see appendix B), as and (2.14) Here we have adopted the standard splitting of the asymptotic charge into a piece linear in the gauge field (usually called soft) and a piece non-linear in the gauge fields (usually called hard). In the construction of asymptotic charges it is usually argued that it is the linearised theory which is used to define the two-form. This is the case in the framework of covariant charges [48,51] and in a related context was noted by Abbott and Deser [50] for Yang-Mills theory. In the current context this would correspond to linearising the theory around the background A a z = A ā z = 0 and would thus discard the non-linear gluon piece of the asymptotic charge. However, it was argued that one should include the non-gluon piece inside the hard part of the charge [52]. As we will explicitly discuss later, we start from the linear charge (2.13) acting on asymptotic states however such a hard contribution in our case is indeed present as it will be generated by the time evolution given by the Møller operator. At I − there is an analogous construction where the linearised asymptotic charge is given by in advanced Bondi coordinates. We will be interested in understanding the action of these charges on scattering states and in order to hew closely to the standard S-matrix formulation we will start with the usual equal-time commutators on a space-like hypersurface and by taking the t → ±∞ limit of the usual free field expansion The sum is over helicities with polarisation vectors ε σ µ and the equal time commutation relations are given in Appendix A. We can then evaluate the fields arbitrarily close to I + by using the saddle point approximation in the r → ∞ limit [29,53,77,78] such that, where q 0 = ω q andx = x/r. The asymptotic charges can now be expressed in terms of the creation and annihilation operators and for the linear, or soft, charge we have which provides the connection with the action discussed in the introduction eq.(1.3). The non-linear part of the charge is more complicated to interpret in the quantum theory as it requires normal ordering. For the gluon part we have where we have normal ordered the expression and dropped the resulting constant. As the energy integral is over only positive values, the delta-function in the second term has no non-zero support in the integration region and it is therefore vanishing. Defining the gluon number density operator we can write the first term in the particularly simple form There is one important subtlety in relating the linear charge to large gauge transformations. Using the equal-time commutators and the expression (2.18) one can show that Thus the gauge transformation of the asymptotic field A a z (u, z,z) is given by This will give the incorrect large gauge transformations using the usual definition of the delta-function but it can be remedied by inserting a factor of two for the zero-mode contribution. This is related to the fact noted in the context of the Abelian theory [31,53,54] that the radiative phase space at I + defined as Γ + = {F a uz , F a uz } will not give the usual linearised large gauge transformation. It is easy to check that using the standard Poisson brackets of the non-Abelian gauge theory The problem is that the zero modes at the boundary correspond to a single real scalar field and not to a complex one: in the zero frequency limit the two helicities are identified and there is a miscounting of the degrees of freedom. This problem can be solved as in [53,54] by imposing additional constraints at the boundaries of I + . Alternatively one can add the factor of two for zero modes as was done by the authors of [31] used in the definition of the QED coherent state operator. We will follow a similar procedure however in our case we modify the charge, adding a factor of two for zero-modes, and leave the coherent state, and the amplitudes, unchanged.
Splitting of the charge contributions according to the energy scale E A key point is that once we pick a scale E we split the Fock space of free particles into the hard and the soft parts H = H E s ⊗ H E h . Indeed we could make a splitting of the energy integral to separate the two contributions where we have introduced λ as an infrared cutoff. As mentioned previously, we will in fact take the linearised charge in the asymptotic region as our starting point and see how the non-linearity emerges from commuting with the evolution operators. As we will see this corresponds to effectively having where we have a non-linear hard term for both gluons and quarks. This is similar to the approach taken in [52].

Asymptotic Hamiltonian and Soft Evolution Operators
The starting point for the coherent state approach to IR divergences [13] is the choice of an appropriate asymptotic Hamiltonian describing the parton dynamics in the far future and far past. We will review, following the arguments of [22], how one can carry out the Faddeev-Kulish approach in the non-Abelian case at the leading order in the IR divergences. We start from the splitting of the standard QCD interaction Hamiltonian in the interaction representation into soft and hard parts: This is done by introducing at each interaction vertex the energy transfer ν where ω i are the energies of the interacting particles with η i = +1 (resp. −1) for incoming (resp. outgoing) particles. We define the soft part of the Hamiltonian as containing only energies below a scale ν < E and we also introduce a lower cut-off λ < ν. The lower energy bound λ is not only required by our use of perturbative QCD but also to have a good definition of the FK states [13] however in the following we will sometimes leave it implicit. It is important to note that the region ν < E contains both soft and collinear subregions. For example [24] considered a cubic vertex with incoming gluon energy ω 1 = |q 1 | and outgoing energies ω i = |q i | for i = 2, 3. In particular ω 3 = ω 2 1 + ω 2 2 − 2ω 1 ω 2 cos θ 12 , where θ 12 is the angle between q 1 and q 2 . The condition (3.2) then defines a hyperbola in ω 2 -cos θ 12 plane which defines regions of small ω 2 and cos θ 12 ∼ 1. Alternatively one can define the asymptotic region by imposing both an angular cutoff θ 12 < Θ and an energy cutoff ω 2 < M . In this way either form of the cutoffs can be used to treat the soft and collinear regions for any of the particles. The soft Møller operators are defined as the standard time-ordered product but in order to isolate the leading IR singularities, which are of the form (g 2 YM log E λ ) k at k-loops, it is useful to transform to frequency space and so write the Hamiltonian as [20] Using this expression we can write the soft-evolution operator as Infrared singularities come from vanishing energy denominators, and in particular leading logarithms come from the region specified by the strong ordering where we can thus approximate To leading order (double pole accuracy) in the IR divergences, for which we can neglect the distinction between Ω E + and Ω E − due to the i0 prescription at leading order, the Møller operator is then given by the frequency-ordered exponential Furthermore, it is possible to show that at this order in the IR: • In each three-gluon vertex we may assume there is always one gluon which is much softer than the others and so we can use the eikonal approximation, • Quarks interact only via eikonal vertices and pair production is neglected because the process g → qq is IR finite, • Four-gluon vertices and ghost contributions can be neglected.
With these assumptions it is possible to write the soft interaction Hamiltonian as a sum of two terms, one depending on quarks and a soft gluon and a second purely gluonic cubic term These can be combined so that where ρ now contains a contribution not only from fermionic matter ρ f but also from the hard gluons with energies ω p > ω q . It is important to note that while the density operator involves gluons which are harder than the soft gluon at that vertex it does not commute with all soft gluon operators. This is fundamentally different from the QED case and corresponds to the non-linear nature of the gauge symmetry. Using the Hamiltonian H E s (t) in (3.11) the soft Møller operator becomes where Π a µ (q) = a a µ (q) − a a † µ (q) is the displacement operator and J a q µ = g YM ωq dp ρ a (p) p µ p · q (3.13) and the dot product denotes the contraction of both Lorentz and gauge indices where appropriate. The exponential is interpreted as being ordered in the soft gluon energies with smaller energies to the right. Here we have switched from the frequency ordering to an energy ordering. For a complete equivalence to the frequency ordering we should also include an ordering on angles and a cutoff on the angular region. However as we are only interested in the leading IR region -which is both soft and collinear-it is sufficient to only consider the energy cutoff. Where this leads to expressions appearing singular an angular cutoff can be mentally added. For our explicit computations in subsequent sections we in fact use dimensional regularisation which simultaneous treats both soft and collinear singularities but the frequency/energy ordering is useful for the physical picture it provides. Finally, the soft evolution operator can be used to define an IR-finite S-matrix S E by removing the IR singularities due to initial and final state interactions as shown schematically in Fig. 2.

Large Gauge Transformations for Soft Evolution Operators
We are interested in understanding the transformation properties of the S-matrix under large gauge transformations and the corresponding Ward identity for amplitudes. As the Figure 3. The Ward identity for the asymptotic charge. The red clouds represent the parton dressing factors comprising soft gluons.
S-matrix relates states in the far future and far past, where the theory is taken to be free, we will use the linearised charge, (2.13), when we compute matrix elements of the commutator where we use the appropriate linearised charge Q ± for incoming and outgoing states. The key to our approach is to use the soft-evolution operator to relate the free theory to the interacting theory and so relate the linearised charge to the non-linear contributions. As the soft-evolution operator can be defined in the quantum theory, this gives a method of defining the correct quantum non-linear corrections to the charge. This can be done by analysing the gauge transformation properties of the soft-evolution operator. Before carrying out the computation, we can see how such non-linear terms imply a Ward identity for the S-matrix: If we had that the evolution-operator transformation was given by a relation of the form for some non-linear charge Q h , then we would have the relation Here we have used dressed states, e.g.
in computing matrix elements. If the linearised charge produces states orthogonal to scattering states constructed on the original vacuum, this becomes This relation, graphically represented in Fig. 3, is what we refer to as the Ward identity. Our goal is to study this relation in the context of QCD using the soft-evolution operators and capturing the leading IR singularities. To this end we consider the explicit expression for the linearised large gauge transformations using the definition of the linearised charge (2.13) with the commutators (2.21). It is convenient to first integrate by parts for the z andz variables and start with the expression In principle we must carefully account for any boundary terms that occur in this step, however at the end of the computation we will undo this integration by parts and so remove the boundary terms again. Using the definition of the energy ordering, we calculate As the charge only involves zero modes, it acts only on the last term in each element of the sum, i.e. E λ dq 1 J 1 · Π 1 , and then only for ω 1 = λ ∼ 0. That is, we are assuming that the charge acts on the mode with energy ω 1 = λ ∼ 0 but not on ω 2 > ω 1 . By straightforward computation one finds where N ( ) = λ dp ρ a (p)N a ( , p) with This expression for N a depends on having performed an integration by parts and so ignores any potential boundary terms. For example it can be seen that such a term vanishes for the case a being a constant. One can undo the integration by parts, and so recover the dropped boundary terms, and write where we use the notationq µ = q µ /ω q and which for the case of constant becomes N a ( , p) = g YM a dωδ(ω). In these expressions we must interpret the δ-functions as taking the soft limit in the appropriate fashion. We have thus found that which is almost what we would have expected from the classical analysis of the non-linear charge. However as , we must correctly order the terms to put the expression in the form (4.2). Performing the non-trivial commutators to one-loop, O(g 2 YM ), we find where (4.14) Here we see that there are corrections to the tree-level expression for Q h which arise due to the non-Abelian nature of the theory. Indeed we can see that there will be further nonlinear corrections at each higher loop orders. These corrections do vanish in the case of a constant , which is to say that the total colour charge does not receive any corrections. Moreover the correction involves a soft-gluon operator which we might expect to have vanishing contribution when computing matrix elements of the IR-finite S-matrix. We now turn to the computation of exactly such matrix elements.

Ward Identities for Dressed S-matrix
To compute matrix elements between dressed states, we make use of the non-trivial fact that, to leading order in the soft divergence, the dressing of external states factorises in colour space [19,20,22] where the coherent-state operator U p i E α i β i (Π) is a functional of the soft gluons only. The coherent-state operator for a parton in the gauge group representation with generators t a αβ is defined by the energy-ordered integral where the dressed gluon field is similarly defined by and U pE ab is the adjoint coherent-state operator. These non-linear equations can be solved iteratively so that to O(g 2 YM ) we have This dressing factor captures the leading-order effects of soft-gluon radiation of each of the hard partons. It includes all-order effects arising from arbitrary numbers of gluons being radiated, as well as loop effects which arise from normal ordering each of the terms.
Tree-level Ward Identity While the calculation is essentially identical to the Abelian case, to fix conventions we will start with the tree-level result. The relevant terms for computing matrix elements of the commutator of the linearised charge with the S-matrix are: where the subscripts on the colour generators indicate the parton leg upon which they act and for convenience we denote, for example, t a α j β j b † β j = t a j b † α j . The two contributions to the matrix element can be represented graphically as in Fig. 4. The first term, corresponding to diagram (4.a), corresponds to the usual absorption or emission of the soft gluon in the linearised charge. We compute these terms by using the soft-gluon theorem, see (1.2) and (1.6), and due to the sign difference in the soft limits of incoming and outgoing gluons the two terms of the commutator add rather than cancel. The next terms, corresponding to diagram (4.b) in Fig. 4, arise from contracting oscillators in the linearised charge with those in the coherent states. These terms cancel due to the sign from the commutator. Thus we find where we have introduced the eigenvalue of the non-linear charge which, as expected, is N a ( , p) as defined in (4.10). As in the QED case the charge acts on the Fock vacuum to produce a state orthogonal to all scattering states built on the original vacuum. This can be seen by keeping only those terms where the charge corresponds to an emitted gluon, and for convenience considering only incoming hard partons. While we now only have half the terms from the commutator, we can still have contributions of the form represented in Fig. 4, however due to a cancellation between diagrams (4.a) and (4.b), one finds at tree-level. (6.a) (6.b) (6.c) Figure 6. The IR divergent contributions at one-loop involving one external parton.
Finite One-loop S-matrix Now let us proceed to a one-loop calculation, and again for now considering only in-particles, we compute the one-loop matrix elements The first term is the usual IR-divergent S-matrix element which arises from diagrams such as (5.a) in Fig. 5, while the remaining terms are the compensating IR-divergent terms from the coherent state which are graphically represented in diagrams (5.b) and (5.c). Here we will focus on those diagrams which involve two external partons and neglect one-parton contributions such as those in Fig. 6. These one-parton contributions are subleading in the IR divergences and so are not needed, but can in fact be reconstructed by insisting on gauge invariance. The second term (5.b) gives an O(g 2 YM ) contribution by using the tree-level soft limit of the S-matrix. In particular it can be written as where we have again dropped the one-parton, i.e. j = k, terms. The next term, corresponding to (5.c), involves contracting two gluons from the dressing factors of two partons. We thus need σε σ µ ε σ ν + terms with two oscillators (5.11) so that we find which is the same as the (5.b) contribution up to an overall factor. We can rewrite the product of polarisation vectors using (A.3) which, if we include the one-parton contributions and impose total colour conservation to remove the c µ dependent terms, becomes σε σ µ ε σ ν = η µν . Hence we find the S-matrix elements in terms of the IR-divergent loop integral I jk = dq p j · p k p j · q p k · q . (5.14) While we have previously considered an energy cut-off to regularise divergences, in order to compare with known results it is useful to instead use dimensional regularisation. Keeping the leading divergence in the parameterˆ = d 2 − 2, the loop integral is given by which encodes the double-pole singularity due to both soft and collinear divergences at leading logarithmic accuracy. To one-loop order and to our accuracy the amplitude is known to be, [45], and hence we see that, as expected, the singular parts cancel in the S-matrix elements. This can be repeated for generic outgoing states and again one finds that the one-loop leading singularities cancel as expected. An analogous approach to removing the singularities would be the Wilson-line dressing which reproduces the same leading IR-divergences [55] (see also [18, 71-73, 98, 99]). One-loop Ward Identity We now want to compute the one-loop correction to the Ward identity using the dressed states. We thus consider the terms t es js λ dq s p js ·Π es (qs) p js ·qs where as before the subscripts on the generators denote the hard parton upon which they act. This expression simplifies significantly as the contributions where the gluon in the charge is contracted with a gluon in the dressing of the partons cancel between the two terms in the commutator. Thus we only need to keep the contributions where the charge contracts with the S-matrix as these terms add. The non-vanishing contributions are graphically represented in Fig. 7 where it can be seen that these contributions closely parallel those of Fig. 5. The result can be straightforwardly computed using the previous tree-level and one-loop results. For example diagram (7.b) gives where Q h,a is defined in (5.7) and I kr is the integral (5.14). The diagram (7.c) gives the same contribution but with a coefficient of − 1 2 , so that I (7.c) = − 1 2 I (7.b) . These terms cancel against contributions such as those in diagram (7.a), which can be computed using the soft limit of one-loop amplitudes given in (1.8). As can be seen in, for example, (5.19) there are contributions with j = k = r that involve three external hard-parton legs. These contributions straightforwardly cancel against the appropriate terms in the first term of the one-loop soft limit involving the one-loop amplitude, that is the iterated tree term, using (5.16) for the leading singular terms in the amplitude. Slightly more subtle are the contributions which involve only two external legs. In particular for the diagram (7.b) there is an ambiguity regarding the order in which one takes the soft limits corresponding to the gluon in the charge and the lowest-energy gluon in the parton dressing. The two orderings are shown schematically in Fig. 8. Diagram (8.b.1) corresponds to first taking the soft-limit for the charge before computing the loop integral involving the soft gluon in the dressing factor. In diagram (8.b.2) we instead take the soft limit for the dressing-factor gluon and then for the gluon in the charge. In particular this gives an ordering of the colour generators which combines with the other two-parton contributions from (9.c) to cancel against the two-parton contributions from (7.a) which come from the soft limit of the one-loop amplitude. The final result is that using the ordering (8.b.2) we find that the commutator vanishes at one-loop to leading singularity There are in principle contributions which only involve a single external parton momentum which we have not included as they are subleading in the IR expansion, however these can again be found by using gauge invariance. Thus we find that, using this prescription, there are no corrections to the Ward identity. If we had chosen the alternative prescription (8.b.1) we would find a correction which is related to the one-loop soft current and which can be written as Figure 9. Non-ambiguous two-parton contributions to one-loop Ward identity. In this case we see that there appears to be a one-loop effect which in principle could be interpreted as a correction to the hard charge. Related to this we must also establish that the linearised soft charge acting on the vacuum produces states that are orthogonal to the dressed scattering states at one-loop, that is we must compute In order to do this, one needs to compute all O(g 3 YM ) contributions, a calculation which closely parallels the computation of the one-loop soft current [44]. There are additional three-parton contributions, see Fig. 10, which combine, using essentially the manipulations from the tree-level calculation, to zero. That is the contributions from Fig. 10 cancel those of Fig. 7. The more non-trivial calculation involves the two-parton contributions (see Fig. 9 and Fig. 11) and there are contributions with the same ordering ambiguity as in Fig. 8. Once again if we make the choice corresponding to (8.b.2) we find that there is a cancellation of the two-parton contributions so that to leading order in the IR divergences. Alternatively if we choose the ordering (8.b.1) of Fig. 8, we find that all the diagrams that arise from dressing the tree-level amplitude cancel amongst themselves. The diagram corresponding to dressing the one-loop amplitude, (10.a) of Fig. 10, cancels the contribution due to the iterated tree term in the soft-limit of the loop Q lin Q lin Q lin (11.a) (11.b) (11.c) Figure 11. The additional two parton one-loop contributions to the orthogonality condition.
amplitude (1.8), so that the term arising from the one-loop soft current remains. Thus we find which is the same as the one-loop contribution to the Ward identity for that choice.
While for simplicity we focused on the case with only incoming hard particles, it is possible to generalize our results to generic in-and out-states by essentially using crossing symmetry. In order to check diagrammatically that crossing is satisfied for dressed states, one can argue that the non-Abelian gluon clouds at leading logarithmic order weakly commute with the S-matrix. The steps are essentially the same as in [34], with the important difference that here we have non-Abelian generators inside the clouds. However even in the non-Abelian theory the coherent state-operators have the property We have explicitly checked that for the ordering prescription (8.b.2) there are no corrections to the Ward identity for generic states, while for ordering (8.b.1) the correction term now has a sum over both outgoing and incoming partons with a corresponding sign.

Discussion and Conclusions
In this work we have computed matrix elements of the commutator between asymptotic states corresponding to hard partons dressed with soft-gluon coherent operators. At tree-level the result following from the soft-theorem is essentially the same as in the Abelian case In the Abelian case the right-hand side of this relation corresponds to the contribution from the hard part of the asymptotic charge and so this expression is equivalent to the conservation law for the full charge. The non-Abelian theory is quite different as the correct definition of the non-linear charge is unclear and these difficulties are likely only exacerbated by including quantum corrections. We choose a pragmatic approach of defining our quantum non-linear hard charge by means of the soft evolution operators and then computing matrix elements to check whether the resulting Ward identities continue to hold. We have seen in Sec. 4 that the hard charge does receive modifications at loop order however they involve soft gluon operators which can be shown to not contribute at one-loop order if we use IR-finite S-matrix elements with a particular ordering prescription.
In the computation of matrix elements of the commutator we noted that there are several order-of-limits issues. First, we have assumed a "weak" definition of the soft charge insertion, where the soft limit for the gluon operator in the charge is taken after the evaluation of the matrix element. This makes the connection between the asymptotic charge and the standard soft limit of amplitudes most direct. A similar issue arose in the case of sub-leading soft-theorems for graviton amplitudes and our approach is analogous to that of [46]. For the case of the dressed states there remains a second ordering ambiguity, which can be seen in diagrams (8.b.1) and (8.b.2) in Fig. 8, as one must choose whether to take the soft-limit in the charge before-ordering O 1 -or after-ordering O 2 -the contribution from the coherent state operator.
These two orderings imply: • Ordering O 1 : The orthogonality relation is broken by pure one loop effects (i.e. no tree-level iterated terms) and the soft-charge no longer produces states orthogonal to scattering states. The Ward identity receives the exact same correction and one may attempt to correct the hard charge at one loop in fashion similar to [56] to preserve the conservation law.
• Ordering O 2 : Using this definition the orthogonality relation is still valid at one-loop and leading singularity and moreover the Ward identity is preserved In this case the interpretation is clearer and we can see that loop effects at this order do not affect the asymptotic symmetries of the S-matrix.
It is quite natural to choose a prescription that preserves the symmetries of the theory where it does not lead to inconsistencies which argues in favor of the second choice. That this second choice is to be preferred on symmetry grounds can be further seen as follows. In this work we have used the dressing following from the evolution operator (3.12). However one could, following [13], modify the dressing so that J a q,f µ = g YM ωq dp ρ a (p) where f µ is a q-dependent vector and whose value could be determined by insisting on physical asymptotic states. Relatedly one can choose a modified vacuum which is an eigenstate of the linearised charge if we also assume that the Fock vacuum is an eigenstate of the linearised charge with eigenvalue 0. Strictly speaking this corresponds to the identification of the Fock vacuum with the dressed vacuum |0 ≡ 0 0 but this is in fact what we have used as we have set Ω E |0 = |0 in our regularised, perturbative calculations. In this sense we may interpret the failure of the orthogonality relation in the case of ordering O 1 as being related to the coherent states no longer being eigenstates of the linearised charge at one-loop and so requiring a one-loop modification. Thus if we make the physically sensible choice that the asymptotic states remain eigenstates of the linearised asymptotic charges we are forced to use ordering O 2 and we thus find that there are no corrections to the Ward identity. While our explicit computations have been at one-loop, the leading IR divergence coherent state construction is valid to all-loop orders and so with the correct ordering prescription there should be no leading IR divergent quantum corrections to the Ward identity at any loop order. More non-trivially it should also be possible to repeat the analysis at subleading order in the IR divergences. The construction of the dressed states is more complicated since single-parton coherent states can no longer be defined and one needs non-factorisable coherent states which take into account two-parton correlations [20]. At sub-leading order the issue of collinear IR-divergences will also become more involved. In this work we have included collinear effects to the order we consider, that is in the softregion, which gives rise to the double pole inˆ , and we see that the asymptotic charges continue to be conserved with the appropriate ordering prescription. At subleading order this will need to be reconsidered as there are also collinear singularities in the non-soft region [88,90,102]. Such collinear divergences play an important role in understanding treelevel asymptotic symmetries where they are related to the algebra of charges, [27,57,58]. Though at loop-level they are less well studied, the coherent state approach can treat such divergences, see [40,41] and the recent work [42,43] which proposed using the soft-collinear effective theory Hamiltonian to generate the asymptotic time evolution. It can in fact be shown quite generally [59] that the dressed S-matrix is completely IR finite and as the proof makes no use of an explicit expression of the asymptotic Hamiltonian it includes any collinear divergences.
The Faddeev-Kulish dressing of asymptotic particles with soft gauge boson clouds is related to the dressing of hard external particle fields with Wilson lines. This method was first introduced by Mandelstam [60] in order to achieve a gauge-invariant formulation of QED. Nowadays, it is, amongst other applications, an important tool to efficiently compute the infrared-divergence structure of scattering amplitudes. See e.g. [55] for a discussion of the virtual soft gluon contribution to multi-leg amplitudes in non-Abelian gauge theories obtained in this formalism. It is also known that the Wilson line approach to soft radiation produces the Faddeev-Kulish dressing when taking the Wilson lines along time-like paths [61]. Furthermore, it was recently used in [62] to study soft photon hair of black holes. It would be interesting to further study the connection between these formulations and asymptotic symmetries in the context of QCD.
Finally in this work we have focused on single insertions of the asymptotic charge but it would be of interest to consider the generalisation to multiple insertions to study the algebra of asymptotic charges. He et al in [57], see also [63,64], showed that the tree-level double-soft limit of two positive helicity gluons could be rewritten as the level zero Kac-Moody algebra. This can in fact even be used to construct a stress-energy tensor for gluons by the Sugawara method and to derive a KZ-like equation for MHV amplitudes [58]. By considering more general tree-level double-soft limits [65,66] the algebra of currents can be extended to include mixed-helicity gluons and sub-leading, in the soft expansion, currents [58]. The coherent state approach may make it possible to understand quantum corrections to the algebra of such currents.

• The free mode expansion of the gluon field is
A a µ (x) = dq ε σ µ (q)a a σ (q)e iq·x + ε σ µ (q)a a † σ (q)e −iq·x , and we use the commutator and are normalised such that Tr(t a t b ) = 1 2 δ ab . Quarks transform in the fundamental representation t a = T a and gluons in the adjoint (t a ) bc = (T a A ) bc = −if a bc .
• The covariant derivative is defined as

B Asymptotic Expansions
In this appendix we review the construction of the asymptotic charges for non-Abelian gauge theory. We follow [27,53,57,67,68] and the review [29]. In order to understand the asymptotic symmetries, one must impose fall-off conditions compatible with the equations of motion and which allow for relevant solutions. We consider Lorenz gauge, which in our coordinates is −∂ u (r 2 A r ) − ∂ r (r 2 A u − r 2 A r ) + γ zz (∂ z Az + ∂zA z ) = 0 (B.1) and impose the conditions at null infinity (r, t → ∞, u = t − r = const) A z = A z (u, z,z) + O(r −1 ) , Az = Az(u, z,z) + O(r −1 ) The corresponding condition on the field strength component is