Top-quark loop corrections in Z + jet and Z + 2 jet production

: The sophistication of current predictions for Z +jet production at hadron colliders necessitates a re-evaluation of any approximations inherent in the theoretical calculations. In this paper we address one such issue, the inclusion of mass e(cid:11)ects in top-quark loops. We ameliorate an existing calculation of Z + 1 jet and Z + 2 jet production by presenting exact analytic formulae for amplitudes containing top-quark loops that enter at next-to-leading order in QCD. Although approximations based on an expansion in powers of 1 =m 2 t can lead to poor high-energy behavior, an exact treatment of top-quark loops demonstrates that their e(cid:11)ect is small and has limited phenomenological interest.


JHEP01(2017)020
The production of a Z-boson in association with jets is of considerable importance as a tool for understanding the Standard Model (SM). The Z+jet process has been proposed as a probe of both parton distribution functions and the high-energy running of the strong coupling, α s . The production of more than one jet is especially important in the environment of the LHC, where typical jet reconstruction algorithms routinely result in multiple jets. This means that Z+jets processes represent significant backgrounds in many searches for New Physics, notably when the Z-boson decays to neutrinos so that it is a source of large missing transverse energy. Therefore they must be predicted precisely within the SM. In order to obtain the level of theoretical precision required to match the small experimental uncertainties [1,2], it is imperative to perform perturbative calculations of Z+jet processes beyond the leading order. The dominant source of corrections arises from QCD, with next-to-leading order (NLO) calculations available for processes involving up to four jets [3][4][5][6][7]. The expected experimental precision of measurements of the Z + 1 jet state has motivated the calculation of this process to the next order (NNLO) [8][9][10][11], so that experimental and theoretical uncertainties in this case are commensurate over a relatively large kinematic range. At this level of theoretical accuracy it is also necessary to have control over corrections arising from the electroweak sector. These effects are known for up to two jets in the final state [12][13][14].
With these results in hand it is important to revisit assumptions and approximations inherent in some of the calculations performed so far. One such approximation relates to the inclusion of the effect of the top quark in one-loop virtual corrections to these processes. Since the mass of the top quark introduces a new scale into the problem, including its effect results in a significantly more complex analytic calculation than the usual case in which all quarks are considered massless. In their classic 1997 paper [15], Bern, Dixon and Kosower (BDK) gave results for such contributions to the Z + 1 and Z + 2 jet processes by performing a large mass expansion in the top-quark mass. Although this approximation was appropriate in the last century, and in particular for e + e − annihilation at LEP energies, it may no longer be appropriate at the LHC and higher energy machines where scales above the top quark mass are probed.
In this paper we shall compute a class of one-loop corrections to Z + 1 and Z + 2 jet processes, specifically considering the effects of fermion loops in which the full dependence on the top-quark mass is retained. The one-loop results for these processes can be obtained with a number of numerical programs, such as MadLoop/Madgraph aMC@NLO [16,17], GoSam [18,19] and OpenLoops [20], since the presence of a massive particle does not complicate the computation in the same way. Although the results of these programs are therefore sufficient for many practical purposes, the existence of analytic results enables a more efficient and numerically-stable evaluation of the one-loop amplitudes. This has been important for their use in other applications, such as the evaluation of NNLO corrections to Z+jet production [9] and matching jet substructure observables at next-to-next-to-leading logarithmic accuracy [21]. JHEP01(2017)020 Figure 1. Examples of fermion loop diagrams contributing to Z + 1 jet production. The only non-zero contribution enters through the axial coupling of the Z-boson to third-generation quarks.
The amplitudes that we have computed may be useful for other crossed, or related, processes and are provided in the appendix. The phenomenological impact of these calculations is assessed for the Z +1 jet case in section 2 and for the Z +2 jet process in section 3.

Top-loop effects in Z + 1 jet production
In the case of Z + 1 jet production, top-quark loop contributions only enter through diagrams such as the ones shown in figure 1. Furry's theorem means that diagrams containing a vector coupling of the Z-boson to the quark loop vanish, so that only the axial coupling contributes. In fact, since we consider all quarks other than the top quark to be massless, due to the opposite weak isospin of up-and down-type quarks, the only contribution from these diagrams comes from the third generation. In the original BDK treatment of these diagrams [15], these contributions are computed in the limit that m t → ∞, with the leading term in a 1/m 2 t expansion retained. We have recomputed these contributions retaining the full top-quark mass dependence; the analytic form of the amplitudes representing this contribution is given in appendix A.
This expansion can be extended to include higher-order terms but in the high-energy regime this can lead to problems since the expansion is properly of the form s/m 2 t , where s becomes large. This is illustrated in figure 2 (left), which shows results obtained using the CT14.NN pdf set [22] with both renormalization and factorization scales equal to H T /2, where H T is the scalar sum of the transverse momenta of all leptons and partons. The leading term in the expansion (as presented in BDK) agrees very well with the exact result over the range shown. Including further terms in the 1/m 2 t expansion spoils this agreement. Although the exact treatment and the 1/m 6 t approximation agree up to jet transverse momenta around 1 TeV, beyond that the approximation is no longer under control and results in a wildly different prediction for the spectrum. The lower panel shows the ratio of the approximation with the leading term to the exact result. The two differ by around 0.7% for a jet with 3 TeV transverse momentum. Since the number of events in this region is negligible this is not a significant difference. We conclude that, although the exact result should be preferred, there is no observable impact on the phenomenology of this process when using only the leading term in the 1/m 2 t expansion. At a 100 TeV collider the differences are more significant, as shown in figure 2 (right). Even for 10 TeV jets, which would be abundant at such a collider, the effect of the approximate top-quark loop is a few percent. Since this is at the same level as the NNLO corrections, it is important that the exact result be available and taken into account. JHEP01(2017)020 Figure 2. The jet p T spectrum for Z+jet production at NLO, computed for the 13 TeV LHC (left) and a 100 TeV pp collider (right). The calculation uses a scale µ r = µ f = H T /2 and no cuts are applied apart from p T (jet) > 25 GeV. The red (solid) histogram corresponds to the exact result while the blue (dot-dash) and magenta (dash) histograms represent the large-mass expansion up to 1/m 2 t and 1/m 6 t respectively, as detailed in the text.
Since the thrust of this paper is to examine the impact of the top quark loops themselves, not to perform an exhaustive analysis of the Z+jet process, we have not considered the uncertainty associated with the choice of scales used in figure 2. Indeed, the effects discussed here may be much smaller than other extant theoretical shortcomings. The residual scale dependence in NNLO calculations of the Z + 1 jet process [8][9][10][11] is typically at the level of a few percent in the bulk, but much larger at high jet transverse momentum. The effect of virtual electroweak corrections is also known to be significant in this region [12][13][14] and approximations for combining the effect of QCD and electroweak corrections gives rise to a further uncertainty.

Top-loop effects in Z + 2 jet production
The Z + 2 jet process is sensitive to a much wider range of virtual corrections that involve a closed loop of top quarks. This is partly due to the fact that the process is represented by two separate parton-level reactions (and all appropriate crossings): The labels q and Q refer to two (possibly distinct) flavours of massless quark. We will refer to these by the abbreviated forms, qqggZ and qqQQZ processes. In the original BDK paper, all of the top-quark loop contributions have been included using the 1/m 2 t expansion.
JHEP01(2017)020  In this work we have computed all of the corrections retaining the full dependence on the top-quark mass. The addition of the mass complicates the analytic form of the amplitudes but we have still obtained relatively compact expressions. This is achieved through the use of analytic unitarity methods for computing one-loop box and triangle coefficients [23][24][25] and by recycling BDK results for the massless case whenever possible. Full details of our calculation, including explicit expressions for all amplitudes, are presented in appendices B and C.
The top-quark loop contributions can be categorized according to the manner in which the Z-boson couples to the partons: 1. Contributions where the Z boson couples to the light quark line. These correspond to vacuum polarization contributions to the qqQQZ process shown in figure 3 and to the loop corrections to the qqggZ process depicted in figure 4. These amplitudes are described in detail in appendices B.3 and C.3.

2.
A vector coupling of the Z-boson to a closed loop of top quarks, occurring in diagrams such as the one shown in figure 5(c). These are only present in the qqggZ process and are described in appendix C.4.
3. An axial coupling of the Z-boson to a closed loop of quarks, as shown in figures 5 and 6 for the qqggZ and qqQQZ processes, respectively. This contribution vanishes for all but the third generation of quarks, whose effect is captured here. For the qqQQZ process these corrections are discussed in appendix B.4 while the corresponding contributions to the qqggZ process are detailed in appendices C.5 and C.6. JHEP01(2017)020  We will now examine the effect of each of these contributions separately, both in the 1/m 2 t approximation used in the BDK form of the amplitudes and with the improved treatment provided by the exact expressions presented here. Our calculation is performed by incorporating our newly-calculated amplitudes in the Monte Carlo program MCFM [26][27][28], which already includes a complete calculation of Z + 2 jet production at NLO that makes use of the BDK loop amplitudes. The expressions for the amplitudes with the exact top-mass dependence are written in terms of the scalar integrals described in appendix D, that are evaluated numerically using the ff [29,30] and QCDLoop [31,32] libraries.
For all of the results in this section we will consider the production of an on-shell Zboson that decays to an electron-positron pair, with no cuts applied to the leptons. This is a choice made for the presentation of our results, and not an intrinsic limitation of MCFM. We use the same pdf set and scale choices as in section 2.

Results: 100 TeV collider
Since we expect the problems associated with the 1/m 2 t expansion used in the original BDK expressions to be exacerbated at high energies, we first present results for a putative 100 TeV proton-proton collider. We define jets using the anti-k T clustering algorithm with a jet separation R = 0.5 and demand that they satisfy, p T (jet) > 500 GeV , y(jet) < 4 .  relatively small deviations in this range. In contrast, approximating the contributions that involve the Z-boson coupling to light quarks in the same way leads to substantial errors at jet transverse momenta of about 3 TeV and higher. The NLO rate is over-estimated by a factor of four for a 10 TeV jet. Using the exact result for the top-quark loops yields a prediction that is essentially unchanged from the one in which they are not included at all.

Results: LHC at √ s = 14 TeV
We now turn to results of more immediate interest, namely predictions for the LHC operating at √ s = 14 TeV. We adjust the jet cuts accordingly and now demand, A comparison of our calculations under these cuts is shown in figure 8. Note that, in comparison to the previous figure, the lower panel has a much smaller scale since we consider transverse momenta for the jet that are much lower. In addition, having observed that the effect of the diagrams in which the Z-boson couples to a top-quark loop is small, in this case we simply show the sum of the contributions from the vector and axial couplings of the Z-boson to top quarks. As expected, at the energies that are accessible at the LHC JHEP01(2017)020 the error made when using the 1/m 2 t approximation is much less severe. Even at a jet transverse momentum of 1 TeV it only results in a 4% deviation from the result with no top-loops included. The effect of the approximation on the cross-section for both jets above 50 GeV is an enhancement of a mere 0.05%.

Conclusions
In this paper we have reviewed the importance of top-quark loops in NLO corrections to Z + 1 jet and Z + 2 jet production. To do so, we have computed the effect of these loops with an exact treatment of the top-quark mass and given analytic forms for all the relevant amplitudes. We find that the effect of these loops is very small and not important for phenomenology at the LHC. For a putative 100 TeV proton-proton collider the effects are more significant and, for the Z +1 jet case, lead to a few percent change in the prediction for jets with transverse momentum of 10 TeV. Attempting to include the effect of these loops by using an expansion in powers of 1/m 2 t leads to the theoretical prediction being overestimated due to poor high-energy behaviour. While this may be at a level that is tolerable at the LHC, it can lead to results at 100 TeV that are incorrect by factors of two or more.

JHEP01(2017)020
Boson Feynman rule Coupling Table 1. Feynman rules and couplings of a photon and a Z to a fermion-antifermion pair. For massless fermions it is convenient to use the left-and right-handed couplings, rather than the vector and axial couplings, so both are shown. Q f is the charge of the fermion in units of the positron electric charge.

Acknowledgments
RKE would like to thank the Fermilab theory group for hospitality during the preparation of this paper. The research of JMC is supported by the US DOE under contract DE-AC02-07CH11359.

A Five point amplitude
In this appendix we consider the five-point amplitudes that enter the calculation of Z +1 jet production. Specifically, we consider the process:

A.1 Tree graphs
We write the tree-level amplitude as, where we have omitted the labels of the electron-positron pair, (5 and 4 respectively). We further define e is the QED coupling, g the QCD coupling, Q q is the charge of quark q in units of e, (the positron charge), and the ratio of Z and photon propagators is given by where M Z and Γ Z are the mass and width of the Z. The definition of the Z/γ * couplings is given in table 1. Colour matrices are normalized such that For the tree amplitude , the result is, We note that this matrix element has the same sign as BDK, eq. (D1), as do all of the results in this section. The i j and [i j] are the normal spinor products for massless vectors, such that i j [j i] = s ij . For details of their definition see refs. [33,34].

A.2 Fermion loop corrections to the tree level amplitude
The one-loop colour decomposition is given by The results for the functions A 5;1 and A 5;2 are given in BDK [15]. Since they do not involve fermion loops we do not repeat them here. The function A 5;3 contains the terms where a Z couples to a loop of quarks via the axial coupling, as shown in figure 1, where our conventions for the overall coupling factors are given in table 1. If we consider all quarks except the top to be massless then there is a net contribution only for the third generation, because of the opposite weak isospin of the up-and down-type quarks.
The result for the leading order interfered with the NLO and summed over colours is given in terms of partial amplitudes as, Note that the axial part A f 5;3 depends on the flavour of the quark (f ) and we have to sum over the contributions of the top and bottom loops. We choose to follow the conventions of the original BDK presentation and write the subleading contributions with permuted momentum labels. In this scheme we further define In this case the terms of order and higher in c Γ are not needed because the amplitude is finite. The result for this amplitude is, including both the top and bottom contributions, where m f is the mass of the quark running in the triangular loop. We shall take the bottom quark to be massless, m b = 0.

JHEP01(2017)020
The function f is the axial triangle function that depends on m f , for which results have been given in ref. [15] and are detailed in appendix F. For a massive quark, such as the top quark, in the special case where one of the legs of the triangle is light-like, we have, (cf. eq. (F.13)) B 0 and C 0 are the scalar integral functions defined in appendix D. In the limit m → ∞ we get The result for a massless quark is, where L 0 , L 1 are the cut-completed functions, Summing over the third generation isodoublet in the limit m b = 0 we get, Keeping only the leading term in the m → ∞ limit given in eq. (A.13), this agrees with BDK, eq. (D.11).
B Six point amplitude, A(1 q , 2Q, 3 Q , 4q, 5ē, 6 e ) We now consider processes with one more parton in the final state, starting with processes containing four quarks. 1

B.1 Tree graphs
The general decomposition for the tree process requires that we include the two terms corresponding to the Z/γ * attaching to one or the other of the quark lines,

JHEP01(2017)020
The result for the tree process is, For the case of identical quark flavours (q = Q) see ref. [36].

B.3 Top loops -vacuum polarization contribution
The one-loop contribution to the unrenormalized vacuum polarization is given by, with c Γ given in eq. (A.10). The contribution of a top quark loop to π(p 2 ) is where T R = 1 2 . Renormalization is effected by performing subtraction at zero momentum transfer (p 2 = 0), so that the effect of the top quark decouples at large momentum transfer. In this scheme both the running of the coupling and the evolution of the parton distributions remain in the five flavour scheme. We find In this subtraction scheme, the renormalized contribution coming from the diagrams shown in figure 3 is, Performing the large mass expansion, in the limit m → ∞ we get, This result agrees with BDK, eq. (12.2).

B.4 Top loops -axial vector coupling contribution
The contribution of the top and bottom quarks to the diagrams shown in figure 6 is, The axial triangle function f is presented in appendix F. In particular, the reduction of f to scalar integrals for the case at hand is given in eq. (F.10).

C.1 Tree graphs
Following ref. [15], the colour decomposition of the tree-level contribution to A 6 is The independent results for helicities of the gluons in the tree amplitude are, cf. BDK eqs. (8.4), (8.9) and (8.15). The remaining helicity combination may be obtained by combining the operations of parity (interchanging i j and [i j]) and charge conjugation (exchanging identities of external fermions and anti-fermions). Thus we have, where the operation flip 1 is defined in eq. (C.3) (and also BDK eq. (6.7)).

JHEP01(2017)020
C.2 General structure at one-loop The one-loop colour decomposition is given by [15] A 1−loop where Q i is the electric charge (in units of the positron charge) of the ith quark and n f is the number of light quark flavours. The partial amplitudes A 6;1 and A 6;3 represent contributions where the Z couples to the fermion line as shown in figure 4. The partial amplitudes A v 6;4 , A ax 6;4 and A ax 6;5 represent the contributions from a photon or Z coupling to a fermion loop through a vector or axial-vector coupling. The full results with massless partons running in the loop have been given in BDK. The addition of this paper is to insert the full top quark mass dependence of A 6;1 , A v 6;4 , A ax 6;4 and A ax 6;5 . The partial amplitudes were further decomposed in the original BDK paper into primitive amplitudes as follows: (1 q , 4q, 2, 3) . (C.10) We must therefore provide new expressions, containing the full top quark mass dependence, for the following quantities: • A t 6 (1 q , 2, 3, 4q), in which the Z boson couples to the light quark line.

JHEP01(2017)020
• A v 6;4 (1 q , 4q; 2, 3), in which the Z boson is radiated from a quark loop through the vector coupling. In our approach it is not useful to perform an additional decomposition into A vs 6 and A vf 6 .
For the quantities A v 6;4 , A ax 6 and A ax,sl 6 we will follow the conventions of the original BDK paper and not present expressions for the momentum labelling as in eq. (C.10), but instead do so for the configuration (1 q , 2q; 3, 4).
The colour sum for e + e − →qqgg in terms of partial amplitudes is, (1 q , 2, 3, 4q) C.3 Result for A t 6 (1 q , 2, 3, 4q) The aim of this section is to calculate the full mass dependence of the quantity A t 6 , which is part of A 6;1 that is defined in eq. (C.10). The relevant diagrams do not contribute to A 6;3 . The minus sign for the fermion loop is included in A t 6 . For this case the only non-zero amplitudes occur when the gluons have the same helicity.

(C.16)
We will not present explicit results for the term A v,t 6;4 since they can be simply related to previously published results for the process gg → ZZ [37]. This exploits the fact that A v,t 6;4 only receives contributions from box (not triangle) diagrams, so that replacing a single Z → ¯ current by a g * → qq one is trivial. We have, This is in accord with the procedure for extracting the vector-vector contribution given in eq. (24) of ref. [37], up to an expected change in the overall factor and a sign to match the conventions of BDK.
C.5 Result for A ax,sl (1 q , 2q, 3 g , 4 g ) The sub-leading colour piece receives contributions from the diagram of the type shown in figure 5a. The full result for the third generation isodoublet is A ax,sl This agrees with BDK eq. (11.12).

JHEP01(2017)020
C.6 Result for A ax (1 q , 2q, 3 g , 4 g ) The most complicated case in which to account for the top-quark mass is the calculation of the leading-colour contribution from a loop of massive fermions with an axial vector coupling to the Z-boson. For a complete isodoublet of massless quarks there is no net contribution of this type since the diagrams precisely cancel between the isospin partners. For the (t, b) isodoublet this is no longer the case once a non-zero mass for the top quark is assumed. The contribution of this isodoublet has been presented, retaining only the leading 1/m 2 t terms in an expansion of the top-quark diagrams, in the paper of BDK. The result for the massless diagrams can be extracted from their eqs. (11.3)-(11.4) and eqs. (11.8)-(11.12), simply by discarding the terms proportional to 1/m 2 t . Our base amplitude can be written as follows, This is an expansion in terms of the scalar box (D x|y|z 0 ), triangle (C x|y 0 ) and bubble (B x 0 ) integrals, defined explicitly in appendix D, as well as a left-over rational part (R). The box and triangle coefficients in that expansion have a further mass expansion, while the bubble coefficients and rational part are independent of the mass m. We use this feature to simplify the presentation of our results by replacing the expansion of eq. (C.21) by the more compact form, The function A ax 6,BDK collects the bubble and rational terms as well as the contribution from the triangle coefficient c 12|34 , all of which may be extracted from the previous calculation of BDK. In the paper of BDK, the bubble coefficients have been re-organized to perform cut completion, leading to more compact expressions. It is thus more efficient to use this compact form as our point of departure in presenting the results. Note though that, in our case, the relevant completed functions will be replaced by combinations of scalar bubble integrals that involve the internal top-quark mass. This is a consequence of the fact that the bubble coefficients are unchanged in the massive case, but the integrals themselves are changed.
Apart from the contribution that can be extracted from the results of BDK, eq. (C.24) also enumerates all of the remaining box and triangle integral coefficients that must be JHEP01(2017)020 specified to complete the amplitudes. Although it appears that we should specify three box integral coefficients and six triangle coefficients this is not the case. A number of relations between the various coefficients can be used to minimize the number of independent expressions that must be given explicitly. The simplest relations are those that just correspond to a relabelling of momenta, for example, An additional simplification is due to the structure of the infrared divergences that are present when m = 0, which requires that the box and triangle coefficients are related. Explicitly, we make use of the identities, 27) together with the partner relation that can be obtained by exchanging labels 3 and 4. The coefficients of the m 2 term in the triangle expansion of eq. (C.23), i.e. c (2) x|y , are related to the rational part, R [25]. 2 We exploit this relation in order to determine the coefficient c (2) 12|34 which would normally require much simplification in an explicit analytic calculation, Explicit results for the remaining independent coefficients will be given below.

C.6.1 Box coefficients
For the box coefficients d (i) it is sufficient to consider only 3 + 4 + and 3 + 4 − helicity combinations. The remaining helicities are obtained from these ones according to, where flip 5 is defined in eq. (C.5).

C.6.2 Triangle coefficients
In general there are six possible kinematic configurations of triangle integrals that may contribute to this partial amplitude. These are: where the third leg is clear from momentum conservation. A summary of the method for determining each of these coefficients is shown in table 2. Note that, since the box integral coefficients d 4|123 , which will be specified in sections C.6.2 and C.6.2 respectively below. c 12|3 coefficients. For the triangle coefficients c 12|3 it is sufficient to consider only 3 + 4 + and 3 + 4 − helicity combinations. The remaining helicities are obtained from these ones according to, where flip 5 is defined in eq. (C.5). As indicated in table 2, the mass-dependent terms in the coefficient vanish: c  c 4|123 coefficients. There are no simple symmetry relations between the c 4|123 coefficients of different helicities. We must therefore specify them all.
The coefficients that appear in the m → 0 limit are simply obtained by using the infrared relation, eq. (C.27). The c (2) coefficients are more complicated:

C.6.3 BDK contribution
The final contribution in eq. (C.24) that must be specified is A ax 6,BDK . This consists of terms representing the contributions of the bubble integrals and rational terms, as well as the mass-independent coefficient of the triangle c 12|34 . Although the bubble and triangle coefficients are the same as in the original BDK paper, the integrals that they multiply are of course the ones with non-zero masses in the loop. Our recasting therefore necessitates the introduction of the following functions related to scalar bubble integrals, such that x = p 2 x , y = p 2 y . In the limit m → 0 these reduce to the standard BDK functions, . The overall sign of our expressions is also opposite to the one of BDK, due to the fact that our result describes the amplitudes for a top quark (τ f 3 = +1/2) rather than the (massless) bottom quark (τ f 3 = −1/2) in BDK. As in the -25 -JHEP01(2017)020 massless case, we only need consider three helicity combinations. The final one is related by, With the preliminaries understood we can make use of BDK eq. (11.3) to write the amplitude with two gluons of positive helicity as, The amplitudes with gluons of opposite helicity are not related by a symmetry, but do share a common structure. We note also that the recasting of BDK eqs. (11.9) and (11.10) also requires the following replacements to be made in the BDK formulae, The auxiliary common quantity is adapted from BDK eq. (11.9) and is given by, Finally, we note that the determination of the triangle coefficient c 12|34 using eq. (C.28) requires knowledge of the rational part of the amplitudes. We do not list these explicitly here since they may be simply obtained from the expressions for A ax 6,BDK through the relation, The scalar integrals themselves are defined as follows, where the denominator function is For the purposes of this paper we take the masses in the propagators to be real. Near four dimensions we use d = 4 − 2 (and for clarity the small imaginary part which fixes the analytic continuations is specified by +i ε). µ is a scale introduced so that the integrals preserve their natural dimensions, despite excursions away from d = 4. We have removed the overall constant which occurs in d-dimensional integrals

E Numerical values of coefficients
The test momenta are, in the notation p = (E, p x , p y , p z ) (in GeV), with p 1 + p 2 + p 3 + p 4 + p 5 + p 6 = 0 .
The results for the various contributions to the A ax 6 partial amplitudes are shown in tables 3-7. We show results for the non-zero box and triangle coefficients as well as the remaining contribution A ax 6,BDK that includes both bubbles and rational terms. The coefficient c (0) 12|34 is not shown explicitly for the amplitudes with opposite gluon helicity even though it is non-zero, since its effect is also included in A ax 6,BDK .

F Axial triangle
The amplitude for a Z coupling to two gluons is denoted by T µνρ AB . We calculate the triangle shown in figure 9, where all momenta are outgoing q 1 + q 2 + q 3 = 0 and q 2 i = 0. The result for the two triangle diagrams shown in figure 9, (including the minus sign for a fermion loop) is, T µνρ AB (q 1 , q 2 ) = i g 2 e 16π 2 δ AB 2v f A Γ µνρ , (F.1)