A Complete Two-Loop, Five-Gluon Helicity Amplitude in Yang-Mills Theory

We compute the integrand of the full-colour, two-loop, five-gluon scattering amplitude in pure Yang-Mills theory with all helicities positive, using generalized unitarity cuts. Tree-level BCJ relations, satisfied by amplitudes appearing in the cuts, allow us to deduce all the necessary non-planar information for the full-colour amplitude from known planar data. We present our result in terms of irreducible numerators, with colour factors derived from the multi-peripheral colour decomposition. Finally, the leading soft divergences are checked to reproduce the expected infrared behaviour.


Introduction
Precision measurements at the Large Hadron Collider provide detailed information about the nature of the strong interaction and its role within the Standard Model. With new data already arriving from Run II, there is a growing need for higher-precision theoretical predictions for a variety of different observables. Recent years have seen considerable progress in this respect. Our ability to make predictions in perturbative QCD now covers most of the relevant 2 → 2 scattering processes at next-to-next-to-leading order (NNLO) precision [1][2][3][4][5][6][7][8][9][10][11][12][13], as well as the example of inclusive Higgs production at N 3 LO [14]. Despite this, processes with more than two particles in the final state remain beyond the reach of current NNLO methods.
Next-to-leading order (NLO) corrections to high multiplicity final states are by now commonplace in phenomenological studies. Such computations are possible thanks to automated techniques, which make use of integrand reduction [15], recursive techniques [16], (generalized) unitarity cuts [17][18][19][20] and the known basis of scalar integrals. Processes with up to five coloured partons [21,22] in the final state are feasible using on-shell methods that, by only working with the physical degrees of freedom, are efficient at controlling the complexity.
Two methods have already been explored in this direction. The first of these is maximal unitarity [41], which generalises the cutting techniques of Britto, Cachazo and Feng [19] and Forde [42] to compute the rational coefficients of the master integrals, incorporating information from IBPs. Maximal unitarity has been used to look at maximal cuts for a variety of high-multiplicity examples in four dimensions [43][44][45][46][47][48][49][50]. The second approach extends the integrand reduction program of Ossola, Papadopoulos and Pittau [15]. The initial steps in this direction [51,52] have now developed into a deeper understanding using the language of computational algebraic geometry [53][54][55][56][57]. The D-dimensional extension of this method has also been understood and applied in the context of the planar two-loop, five-gluon amplitude in QCD with all helicities positive (all plus) [58].
In the context of supersymmetric theories, computational methods based on an analysis of unitarity cuts have enabled a large number of high-loop computations. Other methods have also been developed, mainly in the context of these simplified theories. For example, the colour-kinematics duality of Bern, Carrasco and Johansson (BCJ) [59,60] has been successfully exploited to find the complete colour-dressed four-loop, four-gluon amplitude in N = 4 supersymmetric-Yang-Mills (sYM) [61]. The two-loop, five-gluon amplitude, computed in N = 4 sYM in ref. [62], has since been extended to the non-planar sector and cast into a complete set of numerators satisfying colour-kinematics duality [63]. The integrands in planar N = 4 sYM are known to all loop orders [64,65], and recent studies indicate that this simplicity may extend to the non-planar sector [66,67].
The observation that N = 4 sYM theory and the all-plus sector of QCD 1 are related by a dimension-shifting relation [68] suggests that the all-plus amplitude at two loops could be a useful testing ground for new techniques. Indeed, at two loops the planar sector of the all-plus amplitude was observed to be related to the N = 4 amplitude at the integrand level in a similar pattern to the one-loop story, although additional corrections to the N = 4 sector appeared in the form of one-loop squared (or butterfly) topologies [58]. This fact prompts the question as to how much the techniques applied in supersymmetric cases may help to simplify QCD applications.
In this article, we complete the computation of the two-loop, five-gluon, all-plus helicity amplitude including the non-planar sector. In order to deal with the increase in complexity of the full colour amplitude, we introduce a method to find compact colour decompositions that make full use of the underlying Kleiss-Kuijf (KK) relations [69] in a similar way to the previous treatment at tree level and one loop by Del Duca, Dixon and Maltoni (DDM) [70,71]. We then further exploit the on-shell construction of the irreducible numerators to show 1 The all-plus sector is equivalent to self-dual Yang-Mills at one loop. that all of the non-planar cuts can be obtained from the planar cuts. This is reminiscent of the colour-kinematics duality, and indeed we employ the BCJ relations [59] at tree level to relate planar and non-planar cuts.
Our paper is organised as follows. We first present the colour decomposition of the allplus, two-loop amplitude, exploiting the multi-peripheral decomposition of the underlying tree-level amplitudes. In the next section we describe how the complete kinematic structure can be constructed using knowledge from the planar sector and tree-level identities. After describing a worked example, we present compact results for the full integrand. We perform checks of the universal soft behaviour of the amplitude by evaluating the leading O( −2 ) poles of the integrals in the dimensional regularisation parameter . Finally, we draw some conclusions and discuss some future directions.

Review of irreducible numerators
We follow a multi-loop integrand reduction algorithm [51,53,54,58,72] which uses multivariate polynomial division to find an integrand representation of the two-loop amplitude. This section is intended as a brief overview of the approach; we encourage the reader to refer to the literature for more detailed information.
An integrand-reduced two-loop amplitude has the form where the sum over runs over graphs Γ, which are defined by a specific set of denominators D α (the set {α} labels the propagators in the graph Γ). Associated with each graph in the sum is a colour-dressed irreducible numerator∆ Γ ; these are functions of the external momenta p i , the loop momenta 1 and 2 , and also of the external colour indices a i . Each of these numerators has a colour decompositioñ where we must explicitly determine the permutation sum σ and the associated colour factors C Γ . We will present an algorithm to find a simple colour decomposition for our Yang-Mills amplitudes in the next section. An irreducible numerator ∆ Γ ({p i }, 1 , 2 ) can be written in terms of monomials of irreducible scalar products (ISPs). To determine a set of ISPs, we first choose a spanning set of momenta to expand the scalar products along the lines of the van Neerven-Vermaseren basis [73]. By re-expressing the propagators in terms of this spanning set of scalar products, we can see that many can be written as linear combinations of propagators and can therefore be removed and pushed down into simpler topologies. These scalar products are known as reducible scalar products (RSPs). The remaining scalar products are the ISPs; in general the propagators will be quadratic functions of them. To find a basis set of ISP monomials, these additional quadratic relations are removed using polynomial division with respect to a Gröbner basis of the relations. This technique has been developed in the public code -3 -

JHEP10(2015)064
BasisDet [53]. When working in d = 4 − 2 dimensions we also include three extradimensional ISPs Once a basis set of ISP monomials is identified, their rational coefficients are computed from the generalized unitarity cuts of the amplitude. As we take all the propagators contained in a particular graph on shell, the cut amplitude factorises into a product of treelevel amplitudes summed over internal helicity states. Following our schematic notation we can write this as where the trees A (0) are those associated with each vertex in the graph Γ. Making the distinction between the cut associated with a graph Γ and the irreducible numerator associated with the same graph is crucial for understanding this construction. The irreducible numerator contains only that information which is required on the cut associated with Γ, and which is not captured by irreducible numerators of graphs Γ that are "larger" than Γ, in the sense that the propagators contained in Γ are a proper superset of the propagators contained in Γ. In other words, by applying the cuts in a top-down approach we can isolate each topology systematically subtracting the higher-point singularities.
We will frequently specify the irreducible numerator associated with a graph Γ as ∆(Γ) for clarity; one should remember that the function ∆(Γ) depends on loop and external momenta. Furthermore, throughout this paper we will adopt an index notation for the graph labels which lists the number of propagators in each of the three two-loop branches 1 , 2 and 1 + 2 . In addition, we add extra labels to distinguish between topologies of this type. We follow the convention that the right branch is first index, the left branch the second and finally the central branch in the last entry. For example, the planar pentagonbox reads ∆ 431 = ∆( ), while double box with five legs in a non-planar ordering is written ∆ 331;5L 2 = ∆( ). A complete dictionary between this nomenclature and the graphs relevant for our two-loop, five-point calculation is given in table 1.
As shown in ref. [58], there is only ever a single branch to the set of solutions to the on-shell equation in d dimensions, which simplifies the inversion of the system in eq. (2.4) to find the coefficients of the ISP monomials in ∆ Γ , though at the cost of an increased number of monomials with respect to the four-dimensional case.
There are two important remarks about this construction. The integrand representation of eq. (2.1) is not unique, and there are different choices for both the set of ISPs and the set of monomials. Different sets of spanning vectors will result in different ISPs and the polynomial division requires a choice of monomial ordering.
In the following sections we will exploit two important consequences of this approach. Firstly, we will restrict the form of our irreducible numerators to ensure that the choice of ISPs and monomials satisfy the basic symmetries required by our colour decompositions.

JHEP10(2015)064
Secondly, we will make use of the factorisation of irreducible numerators into ordered treelevel amplitudes. These tree amplitudes satisfy a number of relationships among different orderings. This fact will allow us to determine all non-planar cuts of the two-loop all-plus amplitude from the planar irreducible numerators computed in ref. [58].

Colour decomposition
The main result of this work is the construction of the complete five-point, two-loop, all-plus amplitude in Yang-Mills theory. As we mentioned in section 2, it is necessary to choose a particular colour decomposition. This decomposition picks a set of colour tensors describing the colour structure of the amplitude. At the same time, it specifies an associated set of cut diagrams which must be computed. Each of these cut diagrams is, in turn, associated with a unique irreducible numerator. Thus the colour decomposition that we pick is of central importance, because it determines the set of irreducible numerators that we need to calculate.
In this section, we describe the general algorithm that we used for constructing an appropriate colour decomposition of the amplitude, before applying this algorithm to the specific case of the two-loop five-point amplitude.

Multi-peripheral colour decomposition
Our algorithm is applicable to the general case of an L-loop Yang-Mills amplitude. Following the generalized unitarity principle, we begin by writing the amplitude as a sum over all colour-dressed cuts. Diagrammatically, these cuts consist of vertices formed from colour-dressed tree amplitudes which are joined by on-shell propagators. At two loops the set of colour-dressed cuts can be classified by two basic topologies: the genuine two-loop topologies are shown in figure 2, and the one-loop squared (or butterfly) topologies in figure 3.
The central idea is to build the loop-level colour decomposition using knowledge of the underlying tree-level amplitudes. There are a variety of well-known presentations of these trees. We find it convenient to use the DDM form [70,71] for the tree amplitudes: wheref abc = i √ 2f abc are proportional to the standard structure constants in SU(N c ), the gauge group of our Yang-Mills theory. The main advantage of this form of the amplitude is that it contains (n − 2)! colour structures, as compared to (n − 1)! in the standard tracebased decomposition, for example. This fact helps to reduce the number of generated diagrams; in particular, an algorithm based directly on the trace decomposition of tree amplitudes generates a larger set of diagrams, some of which are rather obscure.
Each of the colour structures in the tree decomposition is a string of group theoretic structure constantsf abc . For an n-gluon amplitude the decomposition is constructed by . . . (a) ( 1 ) It is straightforward to build the loop colour structure from these DDM tree colour structures. The loop structure follows directly from the cut diagram: one simply inserts the DDM trees at the vertices; propagator lines connecting trees indicate that the ends of the DDM combs at either end of the propagator have the same colour index to be summed over. Notice, however, that we must pick two special lines in the DDM form of the tree amplitudes (corresponding to lines 1 and n in figure 1). These lines are on opposite ends of the DDM colour strings, so one can informally think of this choice as picking two lines and "stretching" the colour ordered tree between these two ends. We make canonical choices of which legs to pick as special, depending on the number of propagators that connect to the three-point amplitude. These choices are: • Two propagators In this case, it is natural to "stretch" using the two propagators as the special legs. Thus we build the colour structure by pasting a DDM multi-peripheral colour structure between the two propagators. We must sum over every ordering of the external legs. Pictorially, the operation is show in the upper insert (a) of figure 2. • Three propagators In the case of three propagators, we select two out of the three propagators to be the special lines in the DDM presentation. Notice that this choice hides some of the full symmetry of the diagram. In constructing the DDM tree, the propagators we have selected must be at the end of the multi-peripheral structure; we must sum over the positions of the other legs. The result is a sum of diagrams, as shown in the lower insert (b) of figure 2.
• Four propagators We again choose two propagators to "stretch" the cut amplitude into a DDM tree. At two loops, we only encounter this case in the butterfly topologies. We choose upper and lower propagators on the right side of the diagram as special; by symmetry, the result is the same as if we chose upper and lower propagators on the left of the diagram. The insert of figure 3 sketches out the procedure.
In this way, we build a set of colour structures. The kinematic structure associated with each colour structure is easily understood. Each time we insert a particular DDM colour trace, we also pick up a factor of the associated colour-ordered tree amplitude. Thus, the orientation of the legs in the kinematic diagram, associated to an irreducible numerator, is the same as in the colour structure; of course, the "stretching" procedure does not produce new propagators in the irreducible numerator.
One advantage of using the DDM presentation of the amplitude at tree level is that the KK relations are automatically satisfied. Our procedure recycles this property to loop level: we automatically generate a set of colour diagrams that is KK-independent. Along the way, we generate ordered diagrams for the kinematics. The same procedure works at L loops; the amplitude is expressed as where S i are the symmetry factors of the graphs, D i denote the products of the (inverse scalar) propagators, and the ∆ i are irreducible numerators for appropriate colour factors C i generated through our algorithm. Now let us see this procedure at work in the context of the five-point, two-loop amplitude, which is our main focus.

Five-point, two-loop amplitude
Now we describe the colour structure of the five-point, two-loop amplitude. We concentrate on the diagrams that do not vanish in the all-plus case according to ref. [58] and our calculations in section 4. A generic five-point two-loop amplitude can be constructed by a straightforward extension of the present discussion. Let us first write the amplitude and then explain its content in more detail. We label ∆ i and their colour factors by their diagrams directly in the formula: where the integration operator I acts on every ∆ i as The explicit symmetry factors compensate for the over-counts introduced by the overall sum over permutations of external legs. For convenience, we recapitulate all these ∆ i in table 1, where for each irreducible numerator ∆ i we also show its diagram, that of its colour factor, as well as the set of its non-equivalent permutations. The first three graphs in table 1, ∆ 431 , ∆ 332 and ∆ 422 , are the master diagrams corresponding to the maximal cuts. They are purely trivalent, thus their colour factors are unambiguously defined by their proper graphs.
The next three graphs, ∆ 331;M 1 , ∆ 232;M 1 and ∆ 322;M 1 , have a four-point vertex with two external and two internal edges. The two external legs automatically enter in the permutation sum with two possible orderings, hence multi-peripheral subgraphs are naturally obtained by fixing the internal lines, as in the insert (a) of figure 2. "Stretching" the four-point vertex by these lines gives a master graph for each colour factor.  The following two diagrams, ∆ 331;5L 1 and ∆ 331;5L 2 , share the same graph structure, up to the ordering of the four-point vertex. The apparent asymmetry introduced by our selecting these two diagrams, and omitting the graph with the external leg in the right loop, is an artefact of our colour decomposition. One could make other choices; the KK relations satisfied by the trees and their symmetries ensure that any choice is valid.
To expand the four-point vertex in the ninth graph, ∆ 322;5L 1 , we fixed the internal lines of the "diamond" subdiagram, hence its colour factor is C 332 , but with a minus sign due to one flipped vertex. The other permutation of the four-point vertex corresponds to the same topology and is present in the overall permutation sum with the right permuted colour diagram.

JHEP10(2015)064
The colour factor of the planar graph ∆ 430 follows in a straightforward manner from our algorithm (see figure 3), yielding C 431 as its colour factor. Its descendant ∆ 330;M 1 is more interesting, since it is the only graph in the all-plus case with two four-point vertices. They can be treated independently by linearity of colour decomposition. The external one is thus expanded in the same way as in ∆ 331;M 1 , producing C 430 as an intermediate step.
Expanding the internal four-point vertex gives C 431 again.
The last two graphs, ∆ 330;5L 1 and ∆ 330;5L 2 , share a five-point vertex. To explain their colour factors, let us consider the corresponding colour-dressed cut: .
We obtain the multi-peripheral decomposition of the five-point vertex by fixing the two right-hand loop edges and permuting the other three edges. The graphs in the second line can be vertically flipped to put leg 3 downstairs to match the presentation in table 1. Obviously, an equivalent decomposition could be achieved by fixing the loop edges on the left, which would change the orientation of leg 3 in the superficially non-planar graphs ∆ 330;5L 2 . The S 5 -summation in eq. (3.3) effectively symmetrises the colour structure over the two choices of multi-peripheral decompositions.
In the present paper we can avoid lower levels of the graph hierarchy thanks to the simplicity of the fully symmetric helicity configuration, but it already incorporates all the key elements of the general colour structure.

Kinematic structure
With our colour decomposition at hand, we turn our attention to the kinematic structure of the amplitude. We need to compute an irreducible numerator associated to each diagram in eq. (3.3); as Frellesvig, Zhang and one of the current authors have already computed all the planar irreducible numerators [58], our task is to determine the remaining nonplanar numerators. Of course, these numerators can be computed directly from their cuts. However, as we will see, it is easy to determine the complete set of non-planar irreducible numerators for this amplitude from its planar numerators and the knowledge of tree-level amplitude relations.

Non-planar from planar
The non-planar numerator ∆ 332 = ∆( ) can, of course, be obtained directly from its cut. However, we can avoid calculating this non-planar cut explicitly by relating it to a planar cut. We do so in two steps: first, we coalesce two (ordered) three-point trees into a limit of an ordered four-point tree; then we use the BCJ relations [59] satisfied by the ordered four-point tree to reorder the legs until the complete diagram becomes planar.

JHEP10(2015)064
In more detail, we use the following well-known relation, which is satisfied by on-shell amplitudes in the cuts: where s ij = (p i + p j ) 2 are the standard Mandelstam invariants. Since this identity is of central importance for us, we present a short proof. A four-point tree amplitude can be constructed from two three-point amplitudes using the BCFW recursion relation [74,75]: where hat signs on the right-hand side indicate that the amplitudes are evaluated on complex kinematics for some BCFW shift of external legs. The exact complex value of the shifted internal momentum ( 1+2) is defined by the on-shell condition where the precise expression for α depends on the particular BCFW shift. The key point is thatŝ 12 is a linear function of z, with the property that in the limit s 12 → 0, z → 0. In this limit eq. (4.2) becomes 1+2, 3, 4) .  This proves the identity (4.1). We proceed by applying our identity (4.1) to tree amplitudes inside the non-planar cut, rearranging the diagram until it becomes planar. It is simplest to begin with maximal diagrams, and then to continue to topologies with fewer propagators. We will work through the calculation of ∆ 332 , displayed in figure 4, as an example; we computed all non-planar irreducible numerators using the same technique.
The calculation starts at the level of the cuts:  where p i...j = p i + · · · + p j , and we understand that all internal helicities are implicitly summed over while all exposed propagators are cut. These cuts are decomposed into irreducible numerators as Cut 332 = ∆ 332 (12345, 1 , 2 ) , (4.7a) (4.7b) Using the fact that ( 1 − p 123 ) 2 = −( 2 − p 345 ) 2 on this cut, we see that A similar calculation for ∆ 422 leads to ∆ 422 (12345, 1 , 2 ) = ∆ 431 (12345, 1 , 2 ) . (4.9) So far the obtained non-planar numerators are valid only on their cuts, but they can be extended off shell. We may simply express the numerators in terms of a given set of ISPs and then define off-shell numerators unambiguously through these ISPs. The value of a given numerator depends on the choice of ISP basis off shell (in contrast to the situation on shell, of course). In this way, we determine a valid set of non-planar irreducible numerators. Notice that the ISP monomial choices made in the planar sector, such as the higher powers of µ ij preferred over high powers of ( i ·p j ), will then be easily translated to the non-planar numerators.
It is very useful to maintain the symmetries of the underlying graphs in this off-shell continuation. We achieve this by choosing an appropriate basis of ISPs on a graph-bygraph basis. One engineers the ISP basis such that the loop momentum-dependence of each irreducible numerator is captured by a set of ISPs, which map into one another under the graph symmetries without using any cut conditions. The symmetries of the maximal non-planar graphs, for example, are

JHEP10(2015)064
Graphs ISPs RSPs Table 2. The choices of ISPs and RSPs for the two non-planar masters, where the RSPs are chosen as the propagators of the respective graphs. Additionally, the higher-dimensional ISPs µ ij are shared by all topologies.
These symmetries motivate our choices of ISPs, given in table 2. For instance, the second symmetry of ∆ 332 in (4.10a) leads to a map of ISPs (4.11) After we express loop-momentum dependence in eqs. (4.8) and (4.9) in terms of the ISPs of table 2, using the fact that the RSPs (cut propagators) are zero on shell, we are left with appropriate off-shell irreducible numerators. These are listed in section 5. Note that the function ∆ 422 is not the same as the function ∆ 431 despite the on-shell equation (4.9): different ISPs are chosen to make different off-shell symmetries manifest. We obtained irreducible numerators for lower-level non-planar diagrams in the same way, using the BCJ relations on cuts and extending the results off shell. For the all-plus amplitude at hand we find that many lower-level irreducible numerators vanish. In other words, the higher-level numerators capture the unitarity cut structure of the full amplitude, which is given below.

The full-colour five-gluon all-plus amplitude
In this section we present a complete summary of all kinematic numerators contributing to the colour decomposition in eq. (3.3). We include both planar [58] and non-planar irreducible numerators computed using the technique described in section 4.1. The result is presented unrenormalized including the dependence on the spin dimension D s of the gluon, which is equal to 4 in the FDH scheme and 4 − 2 in CDR [37]. The dependence on the extra dimensional ISPs µ ij = can be collected into three general functions,  The remaining coefficients are expressed using the standard spinor-helicity formalism. In particular, we denote tr 5 = 4i µνρσ p µ 1 p ν 2 p ρ 3 p σ since tr ± = tr 1 2 (1 ± γ 5 )1234 . We also make use of "spurious" directions in order to find compact representations of the integrands, The full amplitude reads  -15 -

JHEP10(2015)064
We have found a representation of the full amplitude with no topologies with fewer than six propagators. We note that there are nonzero cuts at the integrand level, but the resulting integrals are scaleless and hence zero in dimensional regularisation. We have checked additional cuts at levels 2 and 3 to ensure that no nonzero topologies remain. To find an integrand with this property, the ISPs ( 1 ·ω 123 ) and ( 2 ·ω 543 ) in the numerators ∆ 330;5L 2 and ∆ 330;5L 1 were upgraded to include terms proportional to ( 1 − p 1 ) 2 and ( 2 − p 5 ) 2 .

Checking the soft divergences
Since the all-plus helicity configuration is zero at tree level, the universal infrared (IR) structure is the same as that of an ordinary one-loop amplitude. The poles of our twoloop amplitude should therefore be equivalent to those of the finite one-loop amplitude multiplied by the infrared pole operator including the sum over colour correlations [76] A (2) where we have used the • symbol to indicate that the colour matrices (in this purely gluonic case they will all be structure constants) should be inserted into the colour factors of the one-loop amplitude. The standard loop prefactor is given by 2 The one-loop amplitude to all orders in can be found in ref. [68].
There are a number of difficulties in checking eq. (6.1) in full since the five-point planar and non-planar integrals required are still unknown at this time. Resorting to numerical evaluation, as has been done in the planar case [58], is computationally prohibitive for two reasons. Firstly, the full colour expansion contains a large number of dimensionshifted integrals (∼ O(1000)) -an order of magnitude more than the leading colour terms. Secondly, there is no Euclidean region for the complete amplitude, and so contour deformation must be performed for many of these integrals, making them more complicated than the planar cases. This task is probably achievable with public tools like FIESTA [77] and SecDec [78].
There is, however, a much simpler method to check the leading soft singularities up to O( −1 ), which can be done analytically. In the leading soft limit, the colour correlations drop out of eq. (6.1): Clearly this is a weaker check than the full IR poles, but it does require non-trivial properties of the non-planar sector. The butterfly (one-loop squared) topologies, are all finite and JHEP10(2015)064 therefore not relevant for the IR properties. Scattering amplitudes in the soft or eikonal limit have many remarkable structures and universal properties. The interested reader may like to turn to ref. [79] for a recent introduction to the subject.

Evaluating the massless double box in the soft limit
The extra simplicity in our all-plus loop amplitude that sets it aside from most two-loop amplitudes is that the integrals contain at most a single soft divergence, rather than the maximum double soft divergence. We can therefore break our loop amplitudes up into sums of regions with soft singularities and evaluate the amplitude in the limit. In this limit the integral factorises into a product of two one-loop integrals and can be evaluated to extract the leading O( −2 ) divergence.
All of the poles of our amplitude are contained in the topologies proportional to the same dimension shifting numerator: which also has a simple behaviour, Taking the example of the two-loop double box, we find two soft regions by taking the limit of either loop. We find a soft singularity whenever we have two adjacent massless legs in one of the loop integrations. In each case, we factorise into an IR divergent triangle and a dimension-shifted box: [µ 2 11 ] . (6.6b) Recalling the one-loop integrals, and summing the two regions, we quickly arrive at the result:

Soft divergences of the five-point integrals
By following the method described in the previous section, we have derived the complete set of integrals required for the O( −2 ) part of the amplitude. All of the integrals have been checked numerically using the sector decomposition methods implemented in FIESTA [77] and SecDec [78]. Some of these integrals have been computed long ago in 4 − 2 dimensions and can be used to write the full integrals including finite terms via the dimensional reduction identities implemented in LiteRed [80] and IBP relations from FIRE5 [81]. 3 We have performed this task for the planar double box with an off-shell leg. Thus we arrived at the following soft limits for the integrals: Using these results, we have checked that eq. (6.3) does hold as expected for our amplitude (5.4).
In this paper we have explored the impact of tree-level amplitude relations in multi-loop integrand computations. There were two major aspects to our work. Firstly, we exploited the Kleiss-Kuijf relations to find a compact colour decomposition for the two-loop amplitude in terms of multi-peripheral colour factors in an analogous way to the tree-level and oneloop decompositions of Del Duca, Dixon and Maltoni [70,71].
Secondly, we applied the BCJ amplitude relations [59] to relate all non-planar generalized unitarity cuts to the previously computed planar ones. This allowed us to easily generate a compact representation the full colour two-loop, five-gluon, all-plus integrand building on previous planar work [58]. The soft infrared poles of the full amplitude were checked against the well-known universal pole structure.
We hope that the computational methods developed here will be of good use in the necessary extension to more general helicity configurations and other 2 → 3 scattering processes at two loops. They highlight some advantages of relating two-loop integrands to tree-level amplitudes via generalized unitarity cuts. As well as avoiding the large intermediate steps that make Feynman diagram computations at this loop order and multiplicity extremely computationally intensive, we are able to build known on-shell symmetries and relations into the amplitude by construction.
Another interesting aspect of the all-plus amplitude is the continuing connection to the previously known amplitudes in N = 4 sYM. Though the dimension shifting relation observed at one loop no longer holds, the integrands of our full all-plus amplitude and the expressions of Carrasco and Johansson [63] are related by the same dimension shifting operator seen in the planar case. For example, we find that ∆ xyz;T (1 + , 2 + , 3 + , 4 + , 5 + ) = where A and B are some functions of the external kinematics and loop momenta. The second term is proportional to (D s − 2) 2 and is a genuine contribution in QCD not related to N = 4. This additional numerator structure is enough to make the off-shell BCJ symmetries non-trivial to satisfy, even though the N = 4 integrand has been cast in such a form. It is an interesting question as to whether this would be possible for the amplitudes presented here and one that we intend to explore in the future.