Conformal invariance of the one-loop all-plus helicity scattering amplitudes

The massless QCD Lagrangian is conformally invariant and, as a consequence, so are the tree-level scattering amplitudes. However, the implications of this powerful symmetry at loop level are only beginning to be explored systematically. Even for finite loop amplitudes, the way conformal symmetry manifests itself may be subtle, e.g. in the form of anomalous conformal Ward identities. As they are finite and rational, the one-loop all-plus and single-minus amplitudes are a natural first step towards understanding the conformal properties of Yang-Mills theory at loop level. Remarkably, we find that the one-loop all-plus amplitudes are conformally invariant, whereas the single-minus are not. Moreover, we present a formula for the one-loop all-plus amplitudes where the symmetry is manifest term by term. Surprisingly, each term transforms covariantly under directional dual conformal variations. We prove the formula directly using recursive techniques, and check that it has the correct physical factorisations.


Introduction
Particle collisions take place at the Large Hadron Collider at extremely high energies. As a consequence, the masses of the scattered particles can often be neglected, in which case (the Lagrangian of) the Standard Model becomes conformally invariant at classical level.
The implications of conformal symmetry have been extensively explored in position space, where it found many successful applications in the calculation of anomalous dimensions and correlation functions, or in the search for consistency conditions constraining the space of conformal field theories through the operator product expansion. This is a very active field of research. However, much less is known about the consequences of conformal symmetry for scattering amplitudes, which live in on-shell momentum space.
Scattering amplitudes play a crucial role in quantum field theory. They connect theory and experiment, as they are the most basic ingredient in the computation of cross sections. Moreover, being the gauge-invariant building blocks for cross sections, their study often gives precious insights into the structure of the underlying theory. For instance, the dual (super)conformal symmetry of planar N = 4 super Yang-Mills is nowhere to be seen in the Lagrangian, and emerges only at the level of the scattering amplitudes [1].
However, studying the consequences of conformal symmetry for scattering amplitudes is challenging for several reasons. Primarily, scattering amplitudes in quantum field theory are typically divergent at quantum level. The necessity of regulating infrared and ultraviolet JHEP02(2020)019 divergences forces one to introduce a mass scale, which obscures the conformal symmetry of the Lagrangian. This is a major issue, but there is hope. Infrared and ultraviolet divergences factorise separately in a well understood way [2]. It is possible to subtract them, and define a finite "hard part" of the amplitude, in which the regulator can be removed. It is therefore not unreasonable to speculate that it may be possible to find a particular subtraction and renormalisation scheme in which this hard function exhibits, in some way, the underlying conformal properties. Another approach, following [3,4] and references therein, is to study the theory at a conformal fixed point.
One might bypass the fundamental problem of divergences by studying finite loop amplitudes, only to find that the idea of conformal symmetry being broken at loop level only by the dimensionful regulators is actually naïve. Indeed, even for finite loop integrals that are exactly conformally invariant off-shell, the way conformal symmetry is implemented on-shell can be very subtle. In refs. [5][6][7][8][9][10], it was shown that on-shell effects can lead to conformal Ward identities with a non-zero anomaly term.
On top of all this, there is also a practical complication. The generator of special conformal transformations, which is a first-order differential operator in coordinate space, becomes second-order in momentum space. This means that, even to answer the most basic questions -for instance, what is the most general conformal invariant for a given particle configuration? -requires the solution of a system of second-order partial differential equations. As the number of variables increases, this soon becomes a formidable task. Hence the importance of finding, by any means, explicit examples of conformal invariants, to understand and reproduce the mechanism underlying their symmetry.
Despite all the issues, the study of the implications of conformal symmetry for scattering amplitudes is not only of deep theoretical interest, but is also rewarding from the practical point of view. Although a systematic investigation has been started only recently, we can already quote several applications. The anomalous conformal Ward identities of refs. [8][9][10][11], for instance, have been used to compute non-trivial loop integrals, including complicated two-loop five-particle examples. Conformal symmetry can be used also in relation to the rational functions appearing in a scattering amplitude. The on-shell diagrams [12] built by assembling conformally invariant vertices are conformally invariant in their turn. Knowing the conformal objects for a given particle configuration therefore gives a handle on the leading singularities [13] of the corresponding scattering amplitudes, which are related to rational factors appearing in them.
For example, the rational factors in the polylogarithmic part of the two-loop all-plus amplitudes have been computed for arbitrary number of gluons in the planar limit using four-dimensional unitarity [14]. In this approach, the one-loop all-plus amplitudes have no cuts and may be regarded as vertices. The rational factors are then given by onshell diagrams which involve tree-level amplitudes and the one-loop all-plus amplitudes as vertices. The conformal invariance of the one-loop amplitudes, which we prove in this paper, straightforwardly implies that of these rational factors. It is natural to expect the same to hold for the non-planar contributions. Indeed, in the recently computed full two-loop fivegluon all-plus amplitude [15,16], the polylogarithmic part contains the same conformally invariant rational factors which appear in our formula for the one-loop amplitude.

JHEP02(2020)019
The on-shell effects breaking (super)conformal symmetry for finite (super)integrals have been studied in scalar and Yukawa theories, and matter supermultiplet [8][9][10][11]. We aim to make the first step towards understanding the conformal properties of Yang-Mills theory at loop level. In any supersymmetric theory, the gluon scattering amplitudes with no or one negative helicity gluon vanish. Given this, they are the natural first objects to study in a general Yang-Mills theory. Since they vanish at tree-level, they are finite and rational at one loop. The one-loop all-plus amplitudes meet the naïve expectation of being conformally invariant. The single-minus amplitudes do not. 1 In this work, we prove that the one-loop all-plus helicity amplitudes are conformally invariant for any number of gluons. We do this by deriving a new formula in which the symmetry is manifest term by term. We present and discuss the new expression in section 2. Its manifest conformal symmetry is proven in section 3. Remarkably, the new formula has non-trivial properties under dual conformal symmetry as well, which we highlight in section 4. Section 5 is then devoted to the analysis of the analytic properties: we show that the new formula has the correct leading behaviour in the soft and collinear limits, and that it is free of spurious poles. In section 6 we prove our formula by combining two complementary BCFW recursions. We draw our conclusions and discuss the future research directions in section 7.

Manifestly conformal result for the one-loop all-plus amplitude in QCD
The n-gluon all-plus amplitudes vanish at tree level, and are given by finite and rational functions at one loop [19][20][21]. In this section we present a new, manifestly conformally invariant formula for the one-loop all-plus amplitudes.
We start by defining our conventions. We expand the n-gluon all-plus helicity scattering amplitude in powers of the coupling constant g as where p i and a i are the momentum and adjoint colour index of the i-th gluon. In order to make the colour dependence explicit, we adopt the following decomposition into colourordered partial amplitudes [22], where S n is the set of all permutations of n objects, Z n is the subset of cyclic permutations leaving the single trace invariant, and S n;c is the subset which leaves the corresponding double trace invariant. The T a i are the generators of the SU(N c ) gauge group in the fundamental representation, normalised such that Tr T a T b = δ ab .

JHEP02(2020)019
The double-trace components A (1,0) n in eq. (2.2) are determined by the single-trace ones A (1,1) n through U(1) decoupling relations [22]. We will therefore focus on the latter. The single-trace components A (1,1) n are colour ordered, namely they receive contributions only from planar graphs where the cyclic ordering of the external legs is fixed to match that of the generators in the corresponding trace. Because of this, they are called planar (partial) amplitudes, and their singularities can occur only in channels of adjacent momenta. Note that the complete amplitude is invariant under permutations of the external legs. From this, thanks to the cyclic symmetry of the traces, the planar partial amplitudes A (1,1) n inherit symmetry under cyclic permutations of the external legs.
In QCD with n f quarks and gauge group SU(N c ), the planar one-loop all-plus amplitudes in four dimensions take a remarkably compact form [19][20][21], where M is an overall constant factor, and we defined the Parke-Taylor denominator as 3) was conjectured in refs. [19,21] based on the behaviour as external momenta become collinear [23][24][25]. The conjecture was then proven in ref. [20] by off-shell recursion relations [24]. In this paper we derive a different expression, 2 A (1,1) with the dual coordinates defined as and 1|x 1k x km |1 = λ α 1 (x 1k ) αα (x km )α β λ 1 β . As we will show in section 3, the summands C kmn in our formula (2.6) are individually conformally invariant. This implies in a very neat way that the one-loop all-plus amplitude is conformally invariant for any number of external gluons. This is not obvious at all in the original formula (2.3), where the summands are not separately conformally invariant.

JHEP02(2020)019
Moreover, it is intriguing to note that our formula (2.6) is reminiscent of similar formulas for NMHV super-amplitudes in N = 4 super Yang-Mills [1]. In particular, the C kmn show some hints of dual conformal symmetry, which can be highlighted by rewriting them in another way, The dual conformal properties of the summands C kmn will be discussed in section 4.
The vanishing of the all-plus and single-minus tree-level amplitudes implies not only that the one-loop all-plus amplitude is finite and rational, but also that it does not contain multi-particle poles. Only the singularities where two (colour-)adjacent momenta become collinear are allowed. While this is clear in eq. (2.3), we will show in section 5 that the same holds for eq. (2.6) as well.
Finally, note that the cyclic symmetry is not manifest in eq. (2.6). The gluon with momentum p 1 , for instance, appears to have a somewhat special role in the summands (2.7), which is tightly linked to the dual conformal symmetry properties discussed in section 4. This is just an artefact of the specific representation we chose, and different choices are of course possible.

Proof of manifest conformal symmetry
The tree-level gluon scattering amplitudes are conformally invariant. This means that, on top of enjoying Poincaré symmetry, they are invariant under dilatations and special conformal transformations. While the invariance under dilatations is guaranteed by the absence of dimensionful parameters, that under special conformal transformations is more interesting. It means that the tree-level gluon amplitudes are annihilated by the secondorder generator [28] The all-plus amplitudes vanish at tree-level, and are therefore finite at one loop. The absence of divergences, and thus of dimensionful scales to treat them, suggests that conformal symmetry might survive at one loop in this special case. Recent studies [8][9][10][11] have however warned us against very subtle conformal symmetry breaking mechanisms even for finite loop integrals. Indeed, the single-minus amplitudes, which are finite at one-loop just like the all-plus ones, are not conformally invariant.
In this section we show that the all-plus amplitudes are conformally invariant at one loop. We begin by spelling out some low-multiplicity cases. Then we will prove analytically that our formula (2.6) is conformally invariant term by term for any number of external gluons.

JHEP02(2020)019
Let us begin with the simplest case: the four-gluon amplitude. From eq. (2.3), for n = 4, we have Using momentum conservation and Schouten identities, this can be rewritten in many equivalent ways. In particular, This expression, although exhibiting spurious double poles, has the advantage of displaying in a spectacularly manifest way its conformal invariance, In fact, no calculation is needed to show it. If we recall the form of the conformal generator (3.1), it contains a second-order mixed partial derivative with respect to λ i andλ i for all particles i = 1, . . . , 4. The expression (3.3) contains only λ orλ spinors for each particle, and it is therefore trivially annihilated.
Let us now look at the case n = 5. Five-particle amplitudes have a special feature: they depend on a pseudo-scalar invariant, 5 = 4i µνρσ p µ 1 p ν 2 p ρ 3 p σ 4 , which vanishes in all collinear limits. While the constraints on the spurious poles and the collinear limits fix entirely the form of the one-loop all-plus amplitudes for n > 5 [19,21], the dependence on 5 at n = 5 remains elusive. Quite remarkably, conformal symmetry fixes it as well, once the rest of the amplitude is given. The five-gluon amplitude is usually presented in the following way [18], There is no linear combination of the terms appearing in eq. (3.5) that is conformally invariant, other than the one which constitutes the amplitude. As noticed in ref. [15], the five-gluon all-plus amplitude can be rewritten in a manifestly conformally invariant way, Each term in this expression is trivially annihilated by the conformal generator (3.1), in the same way the mechanism worked at n = 4.
Comparing the collinear properties of the original (3.5) and of the new (3.6) formula reveals a tension between manifest conformal symmetry and manifest analytic properties. While eq. (3.5) obviously contains poles only where two colour-adjacent momenta become collinear, eq. (3.6) exhibits spurious poles. What is more, although eq. (3.6) can be rewritten in several manifestly conformally invariant ways, none of them is free of spurious poles. In fact, there are 5 2 = 10 conformal objects of the form [ij] 2 /( ab bc ca ) with i, j, a, b, c all different, and each one of them contains necessarily at least one spurious pole. Very JHEP02(2020)019 interestingly, the residues at the spurious poles vanish in a beautiful way thanks to cancellations between different conformal invariants. This can be shown for arbitrary multiplicity, as we discuss in section 5.
Some readers might recognise that the appearance of spurious poles is a typical feature of BCFW representations. Indeed, we will show in section 6 that our formula (2.6) is constructed from a suitably chosen BCFW recursion.
So far, the mechanism allowing for conformal symmetry has been trivial. More interesting conformally invariant objects appear starting from n = 6. In this case, our formula (2.6) gives  In the first line, we recognise the generalisation of the terms already seen in the n = 4, 5 cases, together with new objects in the second line. The conformal invariance of the latter is less obvious, and was proven in ref. [30].
We now proceed to prove analytically that each summand C kmn in eq. (2.6) is individually conformally invariant for generic multiplicity n. In order to do this, it is convenient to write it as This rewriting in fact exposes a separation in the dependence on the particles. Particles {1, m, m + 1, . . . , n} appear only in a holomorphic way. Their contribution to the generator k αα (3.1) therefore vanishes trivially. The remaining particles are then divided into two sets, {2, . . . , k − 1} and {k, . . . , m − 1}. Since C (1) kmnβσ depends only on the former set, and C (2)βσ kmn only on the latter, the action of the generator of k αα (3.1) splits. Therefore, to show the invariance of eq. (3.8) it is sufficient to prove that k αα C (1)βσ kmn = 0 and k αα C (2)βσ kmn = 0. We show here in some detail the first computation: where the indicesβ andσ are symmetrised, and the labels of the spinors in the parentheses are defined modulus k − 1, namely λ α k mod k−1 = λ α 1 . Next, we apply a Schouten identity and perform standard manipulations to draw out a telescoping sum, In the second line of eq. (3.11) we recognise the telescoping sum with f a,α = 1a a,a+1 λ a+1,α mod k−1 . Since f 1,α = 0 and f k−1,α = −λ 1α , the sum exactly cancels out the first line of eq. (3.11), thus giving which completes the proof. Note that the conformal invariance of the planar amplitudes A (1,1) n implies that of the subleading-colour ones A (1,0) n . Although the reverse is not true, it is still very interesting that the remarkably compact expressions for the subleading-colour components presented in ref. [16] were proven to be conformally invariant for arbitrary n in ref. [30].

Hints for dual conformal symmetry
The numerators of the summands in our formula (2.7) have an intriguing reminiscence of the NMHV super-amplitudes in N = 4 super Yang-Mills. The latter enjoy an additional symmetry, dual (super)conformal symmetry [1]. It is therefore natural to ask whether the one-loop all-plus amplitudes have traces of this symmetry as well.
Dual conformal symmetry is a dynamic symmetry, namely it is not present in the Lagrangian, and emerges only at the level of the scattering amplitudes. Its presence can be most naturally revealed by viewing the scattering amplitudes as functions of the dual coordinates x i , related to the momenta by p i = x i − x i+1 with x n+1 ≡ x 1 . Then, dual conformal symmetry is the ordinary conformal symmetry in dual space.
We find that the one-loop four-gluon all-plus amplitude is indeed a dual conformal invariant, namely it is annihilated by the infinitesimal generator of dual special conformal transformations [1,31], At higher multiplicity, the full all-plus amplitude does not retain any sign of dual conformal symmetry as a whole. Quite remarkably, however, the individual summands C kmn are actually covariant under dual special conformal transformations projected along a specific direction.
The easiest way to see that the one-loop four-gluon all-plus amplitude is dual conformally invariant is to look at the infinitesimal variation generated by K αα (4.1) projected along some direction b, For convenience of the reader we list the infinitesimal transformations of the spinors and of the dual variables: Using these equations it is easy to verify that the one-loop four-gluon all-plus amplitude given by eq. (3.2) or (3.3) enjoys dual conformal symmetry. In fact, it is a Yangian invariant [31]. At n = 5 several spinor brackets have non-adjacent indices, therefore the dual conformal properties are not obvious at first sight. Looking at the general formula (2.6), we see that λ 1 breaks the dual conformal covariance. For instance, 1|x 1k x km |m is dual conformally covariant, but projecting by |1 rather than |m breaks the symmetry, (4.8) We recognise on the left-hand side of eq. (4.8) the factor in the numerator of the summands C kmn (2.7). The same conclusion holds for the full C kmn . The fact that λ 1 alone breaks the symmetry, however, suggests to project the dual conformal variation along a specific direction, in our case p 1 . 4 The existence of such a subgroup of dual conformal symmetry, called directional dual conformal symmetry, has been recently unveiled and used in the computation of certain non-planar loop integrals [32][33][34] (see also ref. [35]). The projection by p 1 in fact removes all the terms which would break the dual conformal covariance, as can be clearly seen by taking b ∝ p 1 in eq. (4.8). This implies JHEP02(2020)019 that the numerator of C kmn (2.7) transforms in a covariant way under dual conformal variations along the direction b = p 1 , for some infinitesimal parameter . Remarkably, the same holds for the denominator, so that each summand C kmn in our formula (2.6) is directionally dual conformally covariant, It is clear, however, that each term in our formula (2.6) transforms with different dual conformal weights, so that the entire formula does not have a symmetry that can be clearly stated. If the whole formula had been covariant along some direction, then -because of the cyclic symmetry of the all-plus planar amplitude -it would have been covariant along any direction. This is what happens in the n = 4 case: the sum contains only one summand, which indeed is exactly dual conformal. At higher multiplicity, the presence of a preferred direction, p 1 , is a mere artefact of the representation of the summands we chose, but has no meaning for the complete amplitude. Nevertheless, it is remarkable that the one-loop all-plus amplitude can be written in such a form that each term, separately, exhibits dual conformal covariance properties.

Analytic structure of the amplitude
There is a tension between manifest conformal symmetry and manifest analytic properties. As discussed in section 2, our formula (2.6) is manifestly conformally invariant, but its analytic properties are obscured. In this section we verify that it has the correct leading behaviour in the soft and collinear limits, and that the spurious poles where non-adjacent momenta become collinear have vanishing residues.
Note that the collinear behaviour and the absence of spurious poles fix entirely the one-loop all-plus amplitudes for n > 5 [19,21]. The fact that our formula (2.6) matches the known n = 4, 5 results and has the correct analytic properties is a proof of its validity.

Soft limit
Let us first consider the soft limit, namely the limit in which all components of a gluon momentum, say p j , go to zero. The leading divergence of the one-loop all-plus planar amplitude factorises into a universal tree-level soft factor and the one-loop amplitude with gluon j removed [36,37], where Note that, for a generic one-loop amplitude, eq. (5.1) would contain also a term with the one-loop soft factor and the tree-level amplitude, but the latter vanishes in the all-plus helicity case.

JHEP02(2020)019
Let us now investigate the leading soft behaviour of our formula (2.6). For this purpose, it is convenient to start from eq. (2.9) for the summands C kmn . The terms with k = j, m = j + 1 in the sum vanish as p j → 0. The others give We see that the single sums in eq. (5.3) extend the second double sum to m = j +1, and the third to k = j + 1, respectively. The resulting expression is therefore that of the one-loop (n − 1)-gluon all-plus amplitude A (1,1) n−1 (1 + , . . . ,ĵ + , . . . , n + ), in perfect agreement with the expected leading soft behaviour (5.1).

Collinear limit
We now check that our expression (2.6) exhibits the correct behaviour as two adjacent momenta, say p j and p j+1 , become collinear, i.e.
for some factor z, with P = p j + p j+1 . At one loop, the leading behaviour of the planar n-gluon all-plus amplitude in the collinear limit (5.4) has the form [23-25, 38, 39] A (1,1) n (. . . , j + , (j + 1) + , . . .) n−1 . . . , P + , . . . , (5.5) where Split (0) − is the tree-level g − → g + g + splitting function, Just like in the soft limit case, the behaviour of the one-loop all-plus amplitude is particularly simple, because the all-plus and single-minus amplitudes vanish at tree-level. In general, there would be a contribution from the one-loop splitting functions and the treelevel amplitudes as well.

Absence of spurious poles
Finally, we know that the planar amplitude does not contain poles where two non-adjacent momenta become collinear. Let us show that eq. (2.6) is free of such poles. Looking at the expression of the summands given by eq. (2.9), it is clear that the only spurious poles are of the type 1/ 1j . The terms of the sum in eq. (2.6) that appear to contain this pole are the ones with k, k − 1 = j or m, m − 1 = j. The term C j,j+1,n does not actually contain any 1/ 1 j pole, as becomes apparent on rewriting starting from eq. (2.9). The remaining terms contributing to the pole are then Through Schouten identities it is however easy to see that so that no 1/ 1 j pole is left.

Proof of the new formula via BCFW recursion
On-shell recursive techniques have proven very successful for the computation of scattering amplitudes. In this section we prove our all-n formula (2.6) for the one-loop all-plus amplitudes through the Britto-Cachazo-Feng-Witten (BCFW) recursion [41]. We start with a briew review of this technique. If an n-particle scattering amplitude A n is finite and rational, its form can be reconstructed from the knowledge of its poles and residues. The idea at the basis of the BCFW recursion is to perform such analysis by shifting by a complex parameter z the external momenta, and studying the amplitude evaluated at the deformed kinematics A n (z) as a function of z.

JHEP02(2020)019
The momenta, say p a and p b , are shifted in such a way that on-shellness and momentum conservation are preserved. This can be achieved for instance by shifting the corresponding spinors asλ Then, if the shifted amplitude A n (z) is a finite rational function with poles at points z i , we can apply Cauchy's theorem to A n (z)/z with a contour encircling infinity. As a result, the original unshifted amplitude A n = A n (z = 0) is rewritten in terms of the residues of the shifted one, where the sum runs over the finite-distance poles z i , i.e. |z i | < ∞, and C ∞ n is the residue at z = ∞, i.e.
Tree-level amplitudes are rational functions of the spinor products, and hence of z. Furthermore, they have only simple poles, corresponding to multi-particle and collinear poles. On such poles a tree-level amplitude A (0) n factorises into the product of lowermultiplicity tree amplitudes, so that eq. (6.2) becomes where σ = ±1 is the helicity of the intermediate state with momentum P i , and i labels all possible factorisations such that the shifted legs a and b appear on opposite sides of the pole. Finally, at tree level it is always possible to find a choice for the shifted momenta (6.1) such that there is no pole at infinity. In such case, eq. (6.4) gives a completely algorithmic recursion relation for the amplitude. At loop level the situation is more complicated. First of all, loop amplitudes are in general not rational. The all-plus and single-minus amplitudes however vanish at tree-level, so that they are finite and rational at one loop.
The second complication is due to the potential presence of double poles, introduced by the one-loop all-plus and all-minus splitting functions regulating the collinear singularities [39]. The latter are however removed by vanishing tree-level amplitudes, so that the one-loop all-plus amplitude has only simple poles. Note however that the one-loop single-minus amplitude, instead, does contain double poles.
Thanks to these simplifications, the BCFW recursion relation for the one-loop allplus amplitudes is identical to the one for tree-level amplitudes given by eq. (6.4), with an additional sum over the two ways of assigning the loop to the pair of lower-multiplicity amplitudes, Finally, there is the issue of the residue at infinity, which produces the surface term C ∞ n in the recursions (6.4) and (6.5). As can be understood from the known result (2.3), there

JHEP02(2020)019
is no shift of the form (6.1) for which the amplitude would vanish as z → ∞. In ref. [42] the authors propose a more refined shift of three spinors which removes the pole at infinity. However, that shift breaks the manifest conformal invariance of the separate terms in the recursion. The ordinary BCFW shift of the form (6.1), on the other hand, is conformally invariant: any term constructed via BCFW recursion with shift (6.1) from conformally invariant lower-multiplicity amplitudes is in its turn conformally invariant. This remarkable property becomes transparent in twistor space [43].
Since our goal is to expose the conformal invariance of the one-loop all-plus amplitude, we prefer to tame the pole at infinity by combining two complementary BCFW shifts [44,45]. The basic idea is that the residue at infinity for one shift may be given by finite-distance residues of a second, auxiliary shift. By combining two suitably chosen BCFW shifts, therefore, it may be possible to reconstruct the whole amplitude, including the surface term.
We now present a proof by induction of our formula (2.6). The proof is based on the recursion relation (6.5) coming from the shift which we denote by 2, 1] z . In section 6.1 we show that the amplitude has a non-vanishing residue at z → ∞ under the deformation 2, 1] z , so that the recursion relation contains a surface term. We find that this surface term is contained in the finite-distance residues of a second, auxiliary BCFW shift, λ n−1 −→λ n−1 = λ n−1 + wλ 1 ,λ 1 −→λ 1 =λ 1 − wλ n−1 , (6.7) labeled by n − 1, 1] w . We then use the latter shift to determine by recursion the surface term of the original recursion based on 2, 1] z . Finally, we solve the recursion for the complete one-loop all-plus amplitude in section 6.2.

Term at infinity
Let us consider the behaviour at z → ∞ of our formula (2.6) under the deformation 2, 1] z (6.6). The numerator in eq. (2.6) is independent of z thanks to the projection by λ 1 . Analysing the denominator, we see that one type of term leads to a non-vanishing behaviour at infinity. The ensuing residue at infinity is

JHEP02(2020)019
Note that, after some manipulation and the use of telescoping series, this can be shown to match the residue at infinity obtained from the original formula (2.3). 5 In this section we show that the surface term C ∞ n (6.8) can be obtained via the auxiliary BCFW recursion n − 1, 1] w (6.7), without assuming the knowledge of the full amplitude.
The surface term given by eq. (6.8) vanishes at infinity under certain other shifts, in particular by deforming momenta p 1 and p n−1 . It is thus natural to expect that we can find C ∞ n from C ∞ n−1 using the BCFW shift n − 1, 1] w (6.7). In general there is some overlap between the finite-distance residues arising from two different shifts, and -provided that the shifts complement each other -the full solution can be reconstructed by careful comparison. Doing this with the two shifts 2, 1] z (6.6) and n − 1, 1] w (6.7) for n = 5, 6, 7, a clear pattern emerges: the n-point surface term C ∞ n is given by the residues at the finite-distance poles of the (n − 1)-surface term C ∞ n−1 under the shift n − 1, 1] w (6.7).
This translates into a homogeneous recursion relation for the surface term C ∞ n , based on the BCFW shift n − 1, 1] w (6.7), which is completely decoupled from the rest of the amplitude. This recursion receives contributions from the two diagrams shown in figure 1, and explicitly reads where the hats denote deformation according to eq. (6.7), P ij = p i + p j , and A MHV is the three-point anti-MHV amplitude, . (6.10) It is understood that the two terms on the right-hand side of eq. (6.9) are evaluated at the values of the deformation parameter w for whichP 2 n−1,n = 0 andP 2 n−2,n−1 = 0, respectively. We assume, by induction hypothesis, that the (n − 1)-point surface term C ∞ n−1 is given by eq. (6.8). Then, the contribution from the diagram in figure 1a, in the first line of eq. (6.9), can be written as  which is exactly C ∞ n as given by eq. (6.8). We have therefore determined that the surface term in the recursion (6.5) is given by eq. (6.8). We can now move on to solve the recursion for the full one-loop all-plus amplitude.

Proof by induction of the new formula
We can now show that our formula (2.6) is the solution of the BCFW recursion relation for the one-loop all-plus amplitude with shift 1, 2] z (6.6).
The n-gluon one-loop all-plus amplitude factorises into the (n − 1)-gluon one, on the left, and the three-gluon anti-MHV amplitude, on the right.
We assume, by induction hypothesis, that the (n−1)-gluon one-loop all-plus amplitude in eq. (6.14) is given by our formula (2.6). We spell it out for convenience of the reader, C kmn , (6.18) where we have recognised in the double sum the C kmn summands as given by eq. (2.9). Finally, we note that the surface term given by eq. (6.8) can be rewritten as Substituting eq. (6.19) into eq. (6.18) produces our formula (2.6) and completes the proof.

JHEP02(2020)019 7 Conclusion and outlook
We showed that the all-plus amplitudes in QCD are conformally invariant at one loop.
The new formula, given by eq. (2.6), together with our proof of the conformal invariance of each summand, makes this manifest for arbitrary number of gluons. This constitutes the main result of our paper. Remarkably, the four-gluon case exhibits also dual conformal symmetry. Together, conformal and dual conformal symmetry make the four-gluon one-loop all-plus amplitude a Yangian invariant. Traces of this propagate also to higher multiplicity: the separate terms of eq. (2.6) transform in a covariant way under directional dual conformal symmetry. Our result has consequences at two loops. The rational factors in the polylogarithmic part of the two-loop all-plus amplitudes are known in the planar limit for any number of gluons [14]. They are given by on-shell diagrams involving tree-level amplitudes and oneloop all-plus amplitudes as vertices. Therefore, the conformal invariance of the one-loop all-plus amplitudes implies that of these two-loop rational factors. We expect the same to hold for the non-planar contributions as well, as can already be seen in the recently computed full two-loop five-gluon all-plus amplitude [15,16].
This work is one step towards unraveling the consequences of conformal symmetry for scattering amplitudes. The latter is in general obscured at quantum level due to ultraviolet and infrared divergences, and due to the on-shell effects revealed by refs. [8][9][10]28]. Nevertheless, we have shown that the all-plus amplitudes are conformally invariant at one loop. Interestingly, we find that the one-loop single-minus amplitudes are not conformally invariant, although finite and rational like the all-plus case. Uncovering the underlying conformal properties of the single-minus amplitudes will be an important milestone.
More generally, we would like to understand how conformal symmetry is implemented in the presence of ultraviolet and infrared divergences. We consider this avenue to be of great interest for future studies.