Positive helicity Einstein-Yang-Mills amplitudes from the Double Copy

All positive helicity four-point gluon-graviton amplitudes in Einstein-Yang-Mills theory coupled to a dilaton and axion field are computed at the leading one-loop order using colour-kinematics duality. In particular, all relevant contributions in the gravitational and gauge coupling are established. This extends a previous generalized unitarity based computation beyond the leading terms in the gravitational coupling $\kappa$. The resulting purely rational expressions take very compact forms. The previously seen vanishing of the single-graviton-three-gluon amplitude at leading order in $\kappa$ is seen to be lifted at order $\kappa^{3}$.

This holds true for an amplitude in an arbitrary matter coupled Einstein-Yang-Mills (EYM) with n gluons and m gravitons and persists to all-loops in supersymmetric theories. In nonsupersymmetric theories, in particular in the "pure" Yang-Mills (YM) and Einstein gravity examples, the leading contribution to these amplitudes is at the one-loop order. The resulting one-loop pure gluon or pure graviton positive helicity amplitudes turn out to be remarkably simple rational functions of the kinematic invariants and spinor products. This is a mandatory property in order to have vanishing unitarity cuts in four dimensions for these amplitudes. In fact, these pure gluon or graviton one-loop amplitudes are identical to the ones one finds in selfdual YM or self-dual gravity, respectively. Here all-multiplicity expressions exists for the pure gluon and pure graviton amplitudes of uniform helicities [2][3][4].
Computing scattering amplitudes involving graviton via Feynman diagrams, as they follow from an expansion of the metric tensor about a flat Minkowski background, is a daunting task, because of the sheer complexity and the infinite number of vertices involved. An important insight that arose from the study of these amplitudes, however, is that they are much simpler than expected and display an intimate connection to amplitudes in YM theory. This relation is known as the double copy, as gravitational amplitudes may be obtained by taking the product of two gauge theory quantities.
The earliest such connection is due to Kawai, Lewellen and Tye (KLT) who expressed gravitational tree-level amplitudes as sums over products of gauge theory amplitudes weighted by Mandelstam invariants, originating from a string theory analysis [5]. A decade ago Bern, Carrasco and Johansson (BCJ) [6] showed that these relations may be understood elegantly in terms of a specific diagrammatic expansion of the amplitudes: If one represents the gauge theory amplitudes in a trivalent fashion such that all kinematic numerators are arranged to obey a Jacobi-like relation mirroring the property of the colour degrees of freedom, then the numerator of a gravitational amplitude follows by simply squaring the gauge theory kinematic numerator. This duality -known as the colour kinematic duality (CKD) -was proven at tree-level [7] and conjectured to hold at loop-level as well [8]. Plenty of examples have been presented which confirm this conjecture and furnish the state of the art technique to compute highest loop orders in supergravities and beyond, see e.g. [9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26]. Generally it is difficult to find numerators that satisfy CKD at looplevel, therefore it has been shown in [27] that a modified double copy can be used to construct gravity integrands where contact terms are included due to the violation of kinematical Jacobi identities. Another one-loop generalization of the KLT-formula has been presented in [28,29].
In this paper we focus on the positive helicity sector of amplitudes for scattering processes involving gravitons interacting with gluons described by Einstein-Yang-Mills theory coupled to a dilaton and axion field (EYM), 1 as it arises as the low energy limit of bosonic string theory. The evaluation of tree-level S-matrix elements in pure EYM has made rapid progress in the recent years and may be considered as completely solved. The key insight was to relate colour-ordered amplitudes of pure EYM to YM theory via an expansion of the form M tree n,m (1, · · · , n; h 1 , · · · , h m ) = α∈Perm(2,··· ,n−1;h 1 ,··· ,hm) which has been proposed by [30] for one graviton and n gluons based on a string computation and proven by field theory techniques [31]. Higher graviton extensions of this result were presented in [31][32][33]. The generalization to the entire single-colour-trace and multi-colour-trace sector has been carried out in [34] and [35] by giving an algorithm to construct the coefficients N (1, α, n).
The colour ordered amplitudes in YM theory on the other hand have been determined in [36], which implies the complete solution of the problem of EYM amplitudes at tree-level. An extension to arbitrary matter-single graviton amplitudes was recently established in [37]. The next step has been to tackle the loop amplitudes in (pure) EYM. In [38] all four point amplitudes at one-loop with one particle species circulating in the loop have been calculated for up to one negative helicity state. To make the calculation tractable the authors used the twoparticle cut method in 4 − 2 dimensions. Restricting to one particle species circulating the loop allowed them to use supersymmetric Ward-Takahashi identities to replace the virtual particles in the loop by complex scalars, since the all-plus and all-plus-but-one-minus helicity amplitudes vanish in supersymmetric and supergravitational theories [1]. However, these identities do not hold for a mixed propagation of gravitons and gluons in the loop and therefore this technique is not applicable in the most general setting.
In this work we shall partly fill this gap by computing in four dimensions all positive helicity amplitudes of EYM at four points to all orders in κ at one-loop precision. 2 We now use the double copy method to generate the integrands of EYM by tensoring the two gauge theory integrands of YM and Yang-Mills coupled to a biadjoint scalar (YM+φ 3 ) as was established in a series of papers [18-20, 33, 39] starting from matter coupled N = 2 supergravities. Focusing on the all-plus sector the calculation is greatly simplified since in YM theory in this case the four point amplitude only receives contributions from the box integrals, which implies that this side of the double copy trivially satisfies the CKD. Therefore in the YM+φ 3 sector also a non duality respecting representation may be used to perform the double copy. The later we generate through usual Feynman diagrammatics.
We can report on rather compact formulae, the summary of our results reads 3 : The set of all reduced Feynman propagators 1/(p 2 − m 2 ) associated to the j th graph are denoted by the inverse of the product α j D α j .
• The numerators n j account for the remaining kinematical dependence of the amplitude. Note that factors of ±i in the Feynman propagators are also absorbed in n j .
The second sum runs over all distinct, nonisomorphic, trivalent graphs Γ and the first one over all |S m | = m! permutations of the external legs. Any overcounting of the j th diagram is removed by the symmetry factor S j . Note that by using the identity 1 = D α j /D α j any graph in a diagrammatic expansion can be made trivalent formally.
Representing the amplitude in the form given in (4) reveals the parallel treatment of colour degrees of freedom c j and kinematical degrees of freedom n j and is especially powerful if one arranges the kinematical numerators in such a way that they obey the same algebraic relations as the corresponding colour factors It has been conjectured by Bern, Carrasco and Johanson (BCJ) [6,8] and shown at tree-level (L = 0) in refs. [7,[41][42][43][44][45] that it is always possible to arrange all the numerators n i of a diagram in such a way that they obey (5).
It is a striking feature of this colour-kinematics duality (CKD) that integrands for gravity amplitudes can be easily constructed from integrands of gauge theories if at least one set of the gauge theory numerators n j orñ j satisfies (5): The set Γ g includes all pairs of trivalent graphs from both gauge theories with numerators n j and n j , respectively. This feature of obtaining gravity amplitude integrands has been proven at treelevel for pure gravity using BCFW recursion relations [46] in four dimensions [7]. Furthermore, gauge invariance of both gauge theory amplitudes implies invariance of the corresponding double copied amplitude under linearized diffeomorphisms, i.e. (6) is an amplitude of some gravity theory [33]. The unitarity method can extend this feature to loop-level by reducing amplitudes containing loops to tree-level amplitudes and demanding that CKD holds for all cuts. However, a general proof is still missing [47] 6 . This double copy (DC) construction can be viewed as a generalization of the Kawai-Lewellen-Tye relations [5] and a loop generalization thereof. Since gauge theory integrands are much easier to calculate than gravity integrands, the DC prescription gives a powerful tool to build gravity integrands. In particular, it is enough to construct a certain number of master integrands, because all other integrands are determined by relation (5). Diagrammatically this can be sketched as Applying the same reasoning the bubble graphs may be reexpressed to triangle graphs. Therefore the master numerators for four points at one-loop are given by the box diagrams since all other diagrams can be obtained by applying the kinematical Jacobi identity (5). 7

The double copy procedure for Einstein-Yang-Mills
In this section we outline the methods we use to obtain the amplitudes presented in (3). The strategy is to write the amplitude in the form (6), for which we have to determine the numerators of pure Yang-Mills (YM) and Yang-Mills coupled to a biadjoint scalar (YM+φ 3 ). Afterwards we perform the integrals by decomposing the tensor integral into scalar Feynman integrals and evaluating these in four dimensions.
We shall start our discussion by defining YM theory and YM+φ 3 theory through the Lagrangians with the definitions We have chosen the normalization Tr T a T b = δ ab for the Hermitian (fundamental) generators T a of the gauge group SU(N ), which implies that the structure constant is given by For the global gauge group with structure constants F ABC it is simply demanded that they shall obey the Jacobi identity and that they are antisymmetric in all of their indices A, B, C, . . .. 8 Dimensional analysis reveals that the coupling constant λ is of mass dimensions one. 7 Here we do not consider the bubble-on-external-leg or tadpole diagrams that vanish in dimensional regularization after integration. 8 After double copying we promote the global gauge group to the local gauge group SU(N ) of EYM.

6
The EYM Lagrangian emerging from the double copy of (8) and (9) reads explicitly [49] where G is the determinant of the metric G µν , the scalar curvature is encoded in the Ricci tensor R andF a µν = i 2 √ −G µν σ F σ,a represents the dual field strength tensor. The scalars ϕ and χ are the dilaton and axion, respectively. 9 The authors of ref. [49] confirmed at tree level that the DC procedure of the theories (8) and (9) give the same amplitudes as the one derived from (11). Clearly the dilaton as well as the axion interact with the gauge fields. To obtain the pure EYM result one has to remove the contributions originating from the two scalars. In [15] it has been shown that for N = 0 supergravity the axion and dilaton can be subtracted by a DC of fields in the fundamental representations with opposite statistics.
After discussing the relevant Lagrangians of the theories, we start with the DC procedure. The first step is to determine both numerators n YM and n YM+φ 3 which can be extracted from the corresponding one-loop integrands of pure YM and YM+φ 3 . It turns out that the four point YM numerator n YM 1 + 2 + 3 + 4 + at one-loop for an all-plus helicity amplitude is very simple because it is exclusively given by graphs with box topology. 10 Let us quickly review this property. The corresponding YM amplitude reads in the colour trace basis where S 4 is the permutation group which is quotiented by the subgroups Z 4 and S 4;3 that leave the single and double traces invariant, respectively. N is the degree of the Lie group SU(N ). In this case the primitive amplitude A 4;1 is related to the partial amplitude A 4;3 by A 4;3 = 6A 4;1 [52]. In [53][54][55][56] A 4;1 has been calculated to all orders in . Its Veltman-Passarino [57] reduced form reads where the scalar-box integral I 4 [µ 4 ; S, T ] is defined in the appendix A.1 and [ij] , ij are the helicity spinor brackets of the external momenta p i , p j . The helicity spinor representation of these read pα α = |i]α i| α (see [40,58,59] for reviews). µ is the fictitious mass of the propagating complex scalar field in the loop which needs to be integrated over [54] in order to emulate 4 − 2 dimensions. 11 Thus the corresponding numerator of the partial amplitude is Using the conventions for the colour-ordered Feynman rules from [40,58,59] one can transform the fundamental generators T a of the single and double trace colour structure into adjoint generators f abc . It turns out that the entire colour structure of (12) is encoded precisely in the chain of structure constants c a 1 a 2 a 3 a 4 , c a 1 a 2 a 4 a 3 and c a 1 a 4 a 2 a 3 , which are given by Feynman graphs with box topology, i.e. c a 1 a 2 a 3 a 4 : This enables us to write the amplitude (12) in the form (4) with the YM numerators n YM 1 + 2 + 3 + 4 + . The exact relation between both numerators reads 12 Furthermore the numerators obey the following symmetry properties These properties imply that for all 24 positions of the external legs the numerator reads as (15). Once the symmetries of the numerator are known it is trivial to prove that they obey colourkinematic duality (CKD). CKD demands the relation (5) to hold which in this case translates into It has been discussed that n YM triangle = 0. Besides, we have shown that all numerators with box topology are the same n YM i + j + k + l + = n YM 1 + 2 + 3 + 4 + . These facts imply that (18) is trivially obeyed. An alternative way to obtain BCJ numerators using light-cone Feynman rules for the one-loop rational pure YM amplitudes was introduced in [60].
Thus, in restricting ourselves to the all-plus sector, only the colour structures (14) appear for the YM amplitudes and their numerators satisfy CKD automatically. Therefore to obtain integrands of EYM we have to collect all diagrams of YM+φ 3 theory which have the same colour structures as (14) and multiply the two numerators divided by the stripped propagators. More details of the DC construction are given in A.2.
After we have analyzed the first gauge theory let us turn to the second one, i.e. YM+φ 3 . The corresponding amplitude of YM+φ 3 at one-loop reads It follows from the Lagrangian (9) that the φ 3 interaction in YM+φ 3 is proportional to the coupling constants gλ. Thus the exponent n ≤ 4 indicates how often the φ 3 -interaction appears in the Feynman diagrammatic decomposition of the amplitude. n YM+φ 3 j shall be computed by Feynman diagrams generated from the Feynman rules of the theory given in appendix A.2.
Once we know both gauge theory numerators we can determine the gravity amplitude via In [18,33] it has been deduced how the global colour structure constant F ABC of YM+φ 3 maps into the local one f abc of EYM. Furthermore, it follows from (19) and (20) that the mapping of the coupling constants from the gauge theories to EYM is given by Thus the first step to obtain the full result in four dimensions for the amplitudes presented in section 3 is to evaluate the contributing numerators n YM+φ 3 j by Feynman diagrams. After the gravity numerators of the amplitude are determined we shall decompose the tensor integral into scalar integrals by using the technique of Veltman-Passarino reduction [57]. This feature is nicely implemented in the Mathematica package FeynCalc [61,62], which is used in the following calculations. The final scalar Veltman-Passarino functions still depend on the fictitious mass squared µ 2 . The next step is then to evaluate the remaining integral over µ 2 . However, as the authors of [54] have shown, we can express the integrals containing µ 2 as higher dimensional loop integrals, which gives an easy analytic way to determine the four-dimensional amplitude. The relevant formulae to this reduction are given in appendix A.1.

Amplitudes
In this section, all positive helicity amplitudes in EYM shall be calculated using the DC construction discussed in section 2. New results are obtained for 1 + 2 + 3 + 4 ++ and 1 + 2 + 3 ++ 4 ++ at order κ 3 and κ 4 , respectively. Furthermore, we also present the κ-corrections at order κ 2 and κ 4 for the four-gluon amplitude 1 + 2 + 3 + 4 + . Note that a DC expression of these integrands for arbitrary helicity configurations has been given in [18] by using the results of [12], however, the integrated amplitude has not been published. In the last part of this section we present the result for 1 ++ 2 ++ 3 ++ 4 ++ . The relevant integrands of YM+φ 3 which shall be used in the DC prescription are collected in appendix A.3.

Amplitudes
We shall start the evaluation by calculating gravitional corrections to the four-gluon amplitude. At first we compute the κ 2 correction to 1 + 2 + 3 + 4 + | κ 2 with the DC method introduced in the previous section. Therefore we have to calculate in YM+φ 3 all Feynman graphs which are proportional to λ 2 according to the DC dictionary presented in (21). In section 2.2 we also pointed out that only the box diagrams are non-vanishing for the all-plus YM amplitude. Hence, only the diagrams which carry the same colour structure as the box diagrams have to be determined in YM+φ 3 . A careful analysis shows that only the graph topologies shown in figure 1 contribute. The integrands in figure 1 have a very simple form and can be obtained by the Feynman rules given in appendix A.2. For example the integrand of the first graph in figure 1 reads: The factors in the denominator D i = Q 2 i + i = q 2 i − µ 2 + i are the denominators of the Feynman propagators. These are defined in appendix A.1 eq. (30). The colour structure reads c abcd = f a ab f b bc f c cd f d da and the global SU(N ) group information is encoded in F ABE F E CD , which will be mapped into the adjoint gauge group generators of EYM in the DC.
The amplitude representation (19) only contains cubic graphs, however, we see that in figure 1 Feynman graphs with quartic vertices also appear, e.g. the first graph in the second line of figure 1 We use the Jacobi-identity c abdc ≡ c abcd + c triangle of (17) in our computations which implies that we have to add the numerators of the first two terms. However, we also know that all numerators in YM theory associated to colour structures different from (14) vanish. This implies that the parts of the integrand which effectively contribute to EYM are § The next step is to insert D 3 D 3 which makes the graph trivalent formally such that the gravity integrand can be obtained by (20). Figure 1: Graphs of these topologies are the only ones that have to be considered in YM+φ 3 at order λ 2 g 4 . The other graphs do not contain the colour structures (14). Curly lines represent propagating gluons and dashed lines represent scalar fields. The internal momenta Q i are d-dimensional.
After collecting all the non-vanishing contributions and building up the complete integrand, we Veltman-Passarino reduce the amplitude with the Mathematica package FeynCalc [61,62]. Using the mapping conventions given in (21) the full amplitude at order κ 2 reads where perm indicates the permutations of the legs 2 and 3 as well as 2 and 4. The kinematical dependence is encoded in the spinor brackets and the Mandelstam variables S = 12 [21], T = 23 [32] and U = 13 [31]. The next correction term 1 + 2 + 3 + 4 + | κ 4 can be obtained by the same technique as for 1 + 2 + 3 + 4 + | κ 2 g 2 YM . All the graphs which contribute are depicted in figure 2. The numerators are fairly simple and are listed in the appendix A.3.1. A straight forward calculation gives the following integrated amplitude where again perm indicates the permutations of legs 2 and 3 as well as 2 and 4 and N g = δ a a = N 2 − 1 is the number of adjoint generators of the Lie algebra.

Amplitudes:
It has been explicitly shown that 1 + 2 + 3 + 4 ++ | κg 3 YM vanishes in four dimensions [38] using both generalized unitarity and the DC. Therefore the DC calculation is not reproduced in this paper. Figure 2: At order g 4 the box, triangle and bubble topologies are non-vanishing after double copying.
We therefore move on to the order κ 3 contribution to this amplitude. It can be proven that for the last graph in figure 3 all box colour structures vanish after applying Jacobi identites.
A typical integrand of figure 3 is of the form where four-dimensional spinor helicity variables are used to represent the polarization vector. Compared to the previous two amplitudes we have now a gauge choice that is encoded in the reference vector r 4 which can be chosen arbitrarily but not such that it is proportional to p 4 . We have done the calculation with the choices r 4 ∈ {p 1 , p 2 , p 3 }. Since the amplitude has to be invariant under different gauges, this represents a powerful crosscheck for the final result. After all integrands of YM+φ 3 have been determined one can construct the gravity integrand using (20). Evaluating and reducing this expression yields the simple result 3.3 Amplitudes: The leading order in κ 2 for the amplitude 1 + 2 + 3 ++ 4 ++ | κ 2 g 2 YM has been determined in [38] using the unitarity based two cut method. We shall see that with the DC prescription (20) Figure 3: To obtain 1 + 2 + 3 + 4 ++ | κ 3 g YM we have to analyze this type of graphs in YM+φ 3 . Figure 4: This pair of graph topologies represent graphs which appear at leading order in λ 2 g 2 for 1 + 2 + 3 ++ 4 ++ .
is much easier obtained. Only the graph topologies drawn in figure 4 have to be evaluated on the YM+φ 3 side. Besides, the calculation can even further be simplified by choosing the reference momenta r i for the gluon polarization vectors at the legs three and four to be the same such that the integrands containing the quartic vertex give zero identically. 13 Thus only the first type of graph from figure 4 has to be determined. The resulting amplitude is given by This result agrees with the expression given in [38] if one inserts the pair of coupling constants (g YM κ/2) 2 and the colour structure N Tr T a T b = N δ ab = f a ab f b ba following from the decomposition of a one-loop amplitude into partial amplitudes. Since the result stated in [38] is the only partial amplitude which contributes, both expression coincide.
The κ 4 contribution can be determined by the graphs given in figure 5. Applying the same steps as before we arrive at the result Here N g represents again the dimension of the adjoint representation of the gauge group. This result has been calculated for the following choices of reference momenta r 3 ∈ {p 1 , p 2 , p 4 }, r 4 ∈ {p 1 , p 2 , p 3 } yielding identical results. 13 This follows from + 3 · + 4 ∼ r 3 r 4 .

13
3 ' 4 2 0 ) 5 6 7 1 ( & Figure 5: At g 4 these graph topologies shall contribute to the gravity amplitude The remaining all-plus amplitude of EYM only contains gravitons as asymptotic states. In principle all the diagrams given in figure 6 contribute to the integrand. It can be seen that the graphs in the first row arise from pure YM theory. It has been discussed in section 2.2 that the calculation of these diagrams simplifies by choosing an appropriate gauge such that the graphs with box topology are the only non-vanishing ones, i.e. the only integrand which survives is exactly given by (15). For example this has been demonstrated in [56] using unitarity cuts. To simplify their analysis, the propagating gluon has been replaced by a complex scalar which is possible due to the supersymmetric Ward-Takahashi identities. Since a complex scalar has two degrees of freedom we can represent this relation diagrammatically by This immediately shows that the graphs in the first and second row are intimately related. Hence we can write the amplitude as the sum of three scalar boxes in 4−2 dimensions, which are defined 8 9 @ A B C Figure 6: These topologies in YM+φ 3 have to be evaluated to obtain in equation (29) of the appendix A.1: In four dimensions the result simply reduces to [34] 12 34 According to [15] we can remove the dilaton and axion by subtracting twice the contribution generated by the adjoint scalar circulating in the loop. For the all-plus amplitude the scalar part is given by [34] 12 34 which implies that the pure EYM result reads [34] 12 34 We note that the pure gravity part of (28) also agrees with [4,63]. 14 14 Note that [12]

Conclusions
In this paper we have calculated the one-loop corrections to EYM for all four point amplitudes with positive helicity configuration. Using the DC of YM and YM+φ 3 as well as the fact that diagrams with box topology are the only ones on the YM side that have a non-vanishing contribution, we have evaluated all the diagrams of YM+φ 3 which contain the colour structure of graphs with box topology. Furthermore, we have checked that our results obtained for 1 + , 2 + , 3 ++ , 4 ++ | κ 2 g 2 YM and the pure gravitational part of 1 ++ , 2 ++ , 3 ++ , 4 ++ | κ 4 coincide with known results in the literature.
Generally, it would be very interesting to be able to obtain the pure EYM results from the amplitudes calculated in this paper. An analysis of possible interaction terms generated by the Lagrangian (11) shows that the axion and dilaton cannot contribute at one loop for a four point amplitude if the number of asymptotic gluon states matches the power of the coupling constant g. The reason is that in this case only the gauge fields can propagate in the loop, because all fields couple at least quadratically to it. Similarly, we were able to seperate the axion and dilaton contribution to the four graviton amplitude from the gluon and graviton contributions. However, all these cases have in common that only one particle type is circulating in the loop. All the remaining amplitudes contain a mixed type of particles circulating in the loop, which makes a subtraction of the axion and dilaton rather difficult.
Therefore it would be very interesting to remove the axion and dilaton contribution from all the amplitudes to obtain pure EYM amplitudes. It might be possible to subtract these fields by a similar method used in [15].
Moreover, one can compute the remaining amplitudes with arbitrary helicity configurations at four points. These can be obtained by using the DC method in the same way as it has been done here because in [12] the colour-kinematics representation for the YM amplitudes at one-loop has been presented for arbitrary helicity configurations. However, the calculation is much more tedious since the numerators given in [12] are non-zero for all possible graph topologies. However, this is left for future work. with and the zero-mass box one can derive exact expressions in arbitrary dimensions. The reason is that these expression are exact to all orders in such that one obtains higher dimensional integrals by simply shifting the value of : which follows from the expression (29). Together with the formula (31), all scalar Veltman-Passarino functions can be expressed in terms of Mandelstam variables S = 12 [21], T = 14 [41], and U = 13 [31]. We shall use the following explicit expressions in four dimensions:

A.2 Feynman rules
We shall give a short review of how we can determine the actual gravity amplitude from the DC of the two theories (8). It turns out that in this case the gravity theory is uniquely determined by its spectrum and its cubic coupling, i.e. its three point amplitude. Scrutinizing the spectrum of the gravity theory (YM + φ 3 ) ⊗ DC YM one can deduce that it contains a graviton, an antisymmetric two-form and the dilaton. However, the latter two can be removed by introducing ghost fields in the double copy of the loop amplitudes. Thus the next step is to normalize the three point vertex in L YM+φ 3 appropriately to obtain EYM. This has been done in [18,33] by explicitly calculating the three point amplitude of N = 2 Einstein-Yang-Mills and then truncating this to the bosonic sector.

A.3 Integrands
All integrands are presented in four dimensions. We shall use the notation of the colour factors from (14). The inverse of the stripped propagators is given by Integrands which are proportional to the colour ordering c adbc shall be written with permuted external momenta p i , i.e.Q 0 = l,Q 1 = p 2 ,Q 2 = l + p 2 + p 3 ,Q 3 = Q 3 . The corresponding stripped propagators are given by (D j ) −1 = (Q 2 j + i ) −1 . We only write the parts of the integrand in YM+φ 3 which contribute to the gravity amplitude in EYM, hence all terms that are not proportional to the colour structure (14) are not displayed.
To simplify the expression for the amplitudes we shall make use of the explicit value of the quadratic Casimir of the adjoint representation c A which for SU(N ) is given by These are the two types of integrands that appear in the computation of 1 A a , 2 B b , 3 C c , 4 D d λ 2 g 4 . The remaining integrands can be obtained by permuting the external legs. In total twelve non-equal box integrands and six non-equal triangle integrands contribute.
For the computation of 1 A a , 2 B b , 3 C c , 4 D d g 4 the following diagrams have to be evaluated. N g = N 2 − 1 is the number of adjoint generators of SU(N ). By permuting the external legs the remaining non-equal integrands can be obtained (six box graphes, twelve triangles and six bubbles).
These four types of integrands can appear in general for the computation of 1 A a , 2 B b , 3 C c , 4 + d g 4 λ . However, once the double copy is performed the last type of graphs are vanishing. For the first type of box graphes and the triangle graphes we have two additional permutations whereas for the other box graphs we have two.