text:
# Two-loop QCD helicity amplitudes for massless quark-massless gauge boson scattering

## 1 Introduction

QCD perturbation theory has proven to be a reliable tool with which to generate theoretical predictions of hard scattering processes that can be compared with experimental results from electron-positron annihilation, electron-proton and hadron-hadron scattering. In such studies, the theoretical analysis of an observable at next-to-leading order (NLO) in $`\alpha _s`$ has become a standard candle, the best theoretical prediction to compare against experiment. These NLO calculations are also used to provide an estimate of the uncalculated theoretical uncertainties via the variation in the renormalisation (and factorisation) scale.

However, the large amounts of data soon to become available from for example, high-energy hadronic collider runs place a demand for more accurate theoretical predictions. To go beyond NLO to the next order in QCD perturbation theory, next-to-next-to-leading order (NNLO), has a direct impact on several aspects of the theoretical predictions. For example , the renormalisation (and factorization) scale uncertainty is significantly diminished, there is a better matching of the parton level jet algorithm used in theory and the hadron level jet algorithm used by experiment and the contributions attributed to non-perturbative power corections may also be reduced.

There are several ingredients that need to be considered to calculate the NNLO corrections for a specific $`22`$ scattering process,

* tree-level $`24`$ partonic amplitudes
* one-loop $`23`$ partonic amplitudes,
* two-loop $`22`$ partonic amplitudes.

The tree-level six-point amplitudes  and one-loop five point amplitudes  have been known for some time while there has been recent progress in computing two-loop amplitudes. Parton distribution functions from a global fit to available data with observables computed consistently at NNLO  coupled with a complete knowledge of the three-loop splitting functions  to generate the evolution of the parton distribution functions at this order are also required.

In this particular paper we are concerned with processes involving two massless quarks and two massless gauge bosons,

$$q\overline{q}gg,q\overline{q}g\gamma ,q\overline{q}\gamma \gamma $$

and the processes related by time reversal and crossing symmetry. In previous work , we computed both the self-interference of the one-loop amplitudes and the interference of the tree and two-loop amplitudes for these processes in conventional dimensional regularisation (CDR). In that work, the matrix-elements were averaged over all external helicities and summed over colours in the high energy limit where all quark masses can be neglected. Here, we extend these calculations to now compute the helicity amplitudes for these processes.

The one-loop helicity amplitudes for parton-parton scattering processes have been known for some time . More recently, there have been a few two-loop helicity amplitudes calculated within CDR with $`d=42ϵ`$ dimensions:

* for the $`gg\gamma \gamma `$ , $`\gamma \gamma \gamma \gamma `$  and $`gggg`$  scattering processes, the explicit contraction of the amplitude with the external polarization vectors was used.
* for the for the $`\gamma \gamma \gamma \gamma `$  and $`e^+e^{}3`$ jets  scattering processes, the individual components of the Lorentz-invariant amplitude decomposition were projected out.

After the explicit helicity structure is obtained via any of these methods, one needs to tackle a large number of two-loop integrals with different topologies and tensor structure. It is now well known (see for example  and references therein), that these integrals can be reduced  to a handful of master integrals using integration-by-parts  and Lorentz-invariance  identities. For the process at hand, the relevant master integrals are 10 massless four-point functions with all legs on-shell which have been calculated in the past few years using a wide variety of methods .

In this paper, we follow a similar appoach to extract the helicity structure of the amplitude as the one outlined in Refs. . We start with a general tensor that describes the two quark-two boson amplitude and use gauge invariance to find relations among the tensor coefficients. These coefficients can be extracted using a set of projectors that are valid at any order in perturbation theory. The projectors saturate the tensor structure and produce scalar integrals. The methodology of reducing these scalar integrals to master integrals is exactly the same as for the calculation of the spin summed matrix elements . Note that the tensor coefficients are evaluated in CDR.

Then, we provide the perturbative expansion for the one- and two-loop tensor coefficients and remove ultraviolet (UV) divergences at each order in $`\alpha _s`$ by renormalisation within the $`\overline{\mathrm{MS}}`$ scheme. The infrared (IR) divergent structure is shown to agree with the one obtained in Refs. . Note that in extracting the coefficients of the general tensor structure, conventional dimensional regularisation is preserved.

Finally, the helicity amplitudes are obtained from the general tensor using standard helicity techniques and within the t’Hooft-Veltman (HV) scheme where the external states are four-dimensional and the internal states are kept in 4-2$`ϵ`$-dimensions . The finite helicity amplitudes are the main new results presented in this paper and we give explicit analytic expressions valid for each helicity configuration and process in terms of logarithms and polylogarithms that are real in the physical domain.

During the writing of this paper, we became aware of an independent calculation of the two-loop helicity amplitudes for quark-gluon scattering . This paper  uses a rather different method to extract the helicity amplitudes in a variety of schemes. Where the two calculations overlap (in the HV scheme), we have compared the finite remainder terms for the two-loop $`q\overline{q}gg`$ processes and find complete agreement.

The paper is organised as follows. In Section 2 we introduce notation and discuss the general tensor structure of the two quark-two boson amplitude. Gauge invariance is used to reduce the number of independent tensor structures. We give projectors to isolate the various coefficients and relate the individual coefficients to the interference of amplitudes. Section 3 shows the perturbative expansion for the tensor coefficients computed within CDR and presents an analysis of their UV and IR divergent structure. The coefficients are renormalised in the $`\overline{\mathrm{MS}}`$ scheme to remove all UV divergences, while the IR poles are predicted using the results of Ref. . Expressions for the one-loop tensor coefficients to all orders in $`ϵ`$ appear in Section 3.2 in terms of one-loop bubble integral and the one-loop box integral in $`d=62ϵ`$. The relationship between the tensor coefficients and the helicity amplitudes in the HV scheme as described in Section 4. The UV and IR structure of the helicity amplitudes is exactly the same as for the tensor coefficients, and finite $`𝒪\left(ϵ^0\right)`$ remainders are obtained by subtracting the predicted pole structure from the explicit calculation both at one-loop and two-loop. The main results of this paper are the finite remainders of the two-loop helicity amplitudes. The expressions are rather lengthy and we confine ourselves to presenting the amplitudes for the $`q\overline{q}gg`$ process in Appendix A. For completeness, the finite one-loop remainders are collected in Appendix B. Results for the processes related by crossing symmetry and those processes involving photons can be obtained from the authors. Finally, in Section 6 we present a discussion of our results and conclude.

## 2 Notation

The generic process we consider in this paper is

$$0q(p_1,\lambda _1)+\overline{q}(p_2,\lambda _2)+b_3(p_3,\lambda _3)+b_4(p_4,\lambda _4),$$
(1)
where the bosons $`b_3`$ and $`b_4`$ are either gluons or photons. All particles are light-like with momenta satisfying,

$$0p_1^\mu +p_2^\mu +p_3^\mu +p_4^\mu ,p_i^2=0.$$
(2)
Physical processes are obtained by crossing particles into the initial state. The associated Mandelstam variables are given by

$$s_{12}=(p_1+p_2)^2,s_{23}=(p_2+p_3)^2,s_{13}=(p_1+p_3)^2,s_{12}+s_{23}+s_{13}=0.$$
(3)
We work in conventional dimensional regularisation treating all external quark and boson states in $`d`$ dimensions and renormalise the ultraviolet divergences in the $`\overline{\mathrm{MS}}`$ scheme. The bare coupling $`\alpha _0`$ is related to the running coupling $`\alpha _s\alpha _s(\mu ^2)`$ at renormalisation scale $`\mu `$, by

$$\alpha _0\mu _0^{2ϵ}S_ϵ=\alpha _s\mu ^{2ϵ}\left[1\frac{\beta _0}{ϵ}\left(\frac{\alpha _s}{2\pi }\right)+\left(\frac{\beta _0^2}{ϵ^2}\frac{\beta _1}{2ϵ}\right)\left(\frac{\alpha _s}{2\pi }\right)^2+𝒪\left(\alpha _s^3\right)\right],$$
(4)
where

$$S_ϵ=(4\pi )^ϵe^{ϵ\gamma },\gamma =0.5772\mathrm{}=\mathrm{Euler}\mathrm{constant}$$
(5)
is the typical phase-space volume factor in $`d=42ϵ`$ dimensions and $`\mu _0^2`$ is the mass parameter introduced in dimensional regularisation  to maintain a dimensionless coupling in the bare QCD Lagrangian density.

The first two coefficients of the QCD beta function, $`\beta _0`$ and $`\beta _1`$, for $`N_F`$ (massless) quark flavours, are

$$\beta _0=\frac{11C_A4T_RN_F}{6},\beta _1=\frac{17C_A^210C_AT_RN_F6C_FT_RN_F}{6},$$
(6)
where $`N`$ is the number of colours and

$$C_F=\frac{N^21}{2N},C_A=N,T_R=\frac{1}{2}.$$
(7)
### 2.1 The general tensor

The most general tensor structure for the amplitude, $`|`$, is

$`|`$ $`=`$ $`\overline{u}(p_1)/p_3u(p_2)(A_{11}ϵ_3.p_1ϵ_4.p_1+A_{12}ϵ_3.p_1ϵ_4.p_2)`$ (8)
$`+`$ $`\overline{u}(p_1)/p_3u(p_2)(A_{21}ϵ_3.p_2ϵ_4.p_1+A_{22}ϵ_3.p_2ϵ_4.p_2)`$
$`+`$ $`\overline{u}(p_1)/ϵ_4u(p_2)(B_1ϵ_3.p_1+B_2ϵ_3.p_2)`$
$`+`$ $`\overline{u}(p_1)/ϵ_3u(p_2)(C_1ϵ_4.p_1+C_2ϵ_4.p_2)`$
$`+`$ $`D_1\overline{u}(p_1)/ϵ_3/p_3/ϵ_4u(p_2)`$
$`+`$ $`D_2\overline{u}(p_1)/ϵ_4/p_3/ϵ_3u(p_2),`$
where the constraints $`ϵ_3p_3=ϵ_4p_4=0`$ have been applied. All coefficients depend on the identity of the gauge bosons and are functions of $`s_{13}`$ and $`s_{23}`$ (and implicitly $`s_{12}=s_{13}s_{23}`$) where $`s_{ij}=(p_i+p_j)^2`$. This tensor structure is a priori $`d`$-dimensional, since the dimensionality of the external states has not yet been specified. The amplitude $`|`$ obeys the Ward identity when one of the gauge boson polarization vectors is replaced with the corresponding boson momentum,

$$|(ϵ_3p_3)=|(ϵ_4p_4)=0.$$
(9)
These constraints yield relations amongst the 10 distinct tensor structures and applying these identities gives the gauge-invariant form of the tensor,

$`|=`$ $`A_{11}(s_{13},s_{23})T_{11}+A_{12}(s_{13},s_{23})T_{12}`$ (10)
$`+`$ $`A_{21}(s_{13},s_{23})T_{21}+A_{22}(s_{13},s_{23})T_{22}`$
$`+`$ $`B(s_{13},s_{23})T,`$
where $`A_{IJ}`$ and $`B`$ are five gauge-independent functions and the tensor structures $`T_{IJ}`$ and $`T`$ are given by

$`T_{1J}`$ $`=`$ $`\overline{u}(p_1)/p_3u(p_2)ϵ_3.p_1ϵ_4.p_J{\displaystyle \frac{s_{13}}{2}}\overline{u}(p_1)/ϵ_3u(p_2)ϵ_4.p_J{\displaystyle \frac{s_{J4}}{4}}\overline{u}(p_1)/ϵ_3/p_3/ϵ_4u(p_2),`$ (11)
$`T_{2J}`$ $`=`$ $`\overline{u}(p_1)/p_3u(p_2)ϵ_3.p_2ϵ_4.p_J{\displaystyle \frac{s_{23}}{2}}\overline{u}(p_1)/ϵ_3u(p_2)ϵ_4.p_J{\displaystyle \frac{s_{J4}}{4}}\overline{u}(p_1)/ϵ_4/p_3/ϵ_3u(p_2),`$ (12)
$`T`$ $`=`$ $`s_{23}(\overline{u}(p_1)/ϵ_4u(p_2)ϵ_3.p_1+{\displaystyle \frac{1}{2}}\overline{u}(p_1)/ϵ_3/p_3/ϵ_4u(p_2))`$ (13)
$``$ $`s_{13}(\overline{u}(p_1)/ϵ_4u(p_2)ϵ_3.p_2+{\displaystyle \frac{1}{2}}\overline{u}(p_1)/ϵ_4/p_3/ϵ_3u(p_2)),`$
with $`s_{14}=s_{23}`$. Each of the tensor structures satisfies the Ward identities. The tensor coefficients are further related by symmetry under the interchange of the momenta (and colour labels) of the two gauge bosons,

$`A_{21}(s_{13},s_{23})`$ $`=`$ $`A_{12}(s_{23},s_{13}),`$
$`A_{22}(s_{13},s_{23})`$ $`=`$ $`A_{11}(s_{23},s_{13}),`$
$`B(s_{13},s_{23})`$ $`=`$ $`B(s_{23},s_{13}).`$ (14)
### 2.2 Projectors for the tensor coefficients

The coefficients $`A_{IJ}`$ and $`B`$ may be easily extracted from a Feynman diagram calculation, using projectors such that

$$\underset{\mathrm{spins}}{}𝒫(X)|=X(s_{13},s_{23}).$$
(15)
The explicit forms for the projectors in $`d`$ space-time dimensions are,

$`𝒫(A_{11})`$ $`=`$ $`{\displaystyle \frac{d}{2s_{12}^2s_{13}s_{23}(d3)}}T_{11}^{}+{\displaystyle \frac{(d4)}{2s_{12}^2s_{13}s_{23}(d3)}}T_{22}^{}+{\displaystyle \frac{1}{2s_{12}^2s_{13}s_{23}(d3)}}T^{}`$
$`+{\displaystyle \frac{(2s_{13}(d4)s_{23})}{2s_{12}^2s_{13}^2s_{23}(d3)}}T_{12}^{}+{\displaystyle \frac{(2s_{23}(d4)s_{13})}{2s_{12}^2s_{13}s_{23}^2(d3)}}T_{21}^{},`$
$`𝒫(A_{12})`$ $`=`$ $`{\displaystyle \frac{(2s_{13}(d4)s_{23})}{2s_{12}^2s_{13}^2s_{23}(d3)}}T_{11}^{}+{\displaystyle \frac{(2s_{13}(d4)s_{23})}{2s_{12}^2s_{13}^2s_{23}(d3)}}T_{22}^{}{\displaystyle \frac{(2s_{13}(d4)s_{23})}{2s_{12}^2s_{13}^2s_{23}(d3)(d4)}}T^{}`$
$`+{\displaystyle \frac{(4(d3)s_{13}^2+4(d4)s_{13}s_{23}+d(d4)s_{23}^2)}{2s_{12}^2s_{13}^3s_{23}(d3)(d4)}}T_{12}^{}`$
$`+{\displaystyle \frac{(2(d4)s_{12}^2+(d2)(d6)s_{13}s_{23})}{2s_{12}^2s_{13}^2s_{23}^2(d3)(d4)}}T_{21}^{},`$
$`𝒫(B)`$ $`=`$ $`{\displaystyle \frac{1}{2s_{12}^2s_{13}s_{23}(d3)}}T_{11}^{}{\displaystyle \frac{1}{2s_{12}^2s_{13}s_{23}(d3)}}T_{22}^{}+{\displaystyle \frac{1}{2s_{12}^2s_{13}s_{23}(d4)}}T^{}`$
$`{\displaystyle \frac{(2s_{13}(d4)s_{23})}{2s_{12}^2s_{13}^2s_{23}(d3)(d4)}}T_{12}^{}+{\displaystyle \frac{(2s_{23}(d4)s_{13})}{2s_{12}^2s_{13}s_{23}^2(d3)(d4)}}T_{21}^{}.`$
The projectors for $`A_{21}`$ and $`A_{22}`$ are obtained by the permutations

$`𝒫(A_{22})`$ $`=`$ $`𝒫(A_{11})\{s_{13}s_{23},T_{11}^{}T_{22}^{},T_{12}^{}T_{21}^{},T^{}T^{}\},`$
$`𝒫(A_{21})`$ $`=`$ $`𝒫(A_{12})\{s_{13}s_{23},T_{11}^{}T_{22}^{},T_{12}^{}T_{21}^{},T^{}T^{}\}.`$
### 2.3 Reconstruction of the matrix elements

Starting from the tensor coefficients, the interference (or square) of amplitudes can be easily extracted in a variety of schemes. In the HV scheme, the external particles are treated in 4-dimensions and we find that,

$`{\displaystyle \underset{\mathrm{spins}}{}}|`$ $`=`$ $`+{\displaystyle \frac{s_{13}^3s_{23}}{2}}\left(2BA_{11}+A_{12}\right)\left(2BA_{11}+A_{12}\right)^{}`$ (18)
$`+{\displaystyle \frac{s_{13}s_{23}^3}{2}}\left(2BA_{21}+A_{22}\right)\left(2BA_{21}+A_{22}\right)^{}`$
$`+{\displaystyle \frac{s_{13}s_{23}s_{12}^2}{2}}\left(A_{11}A_{11}^{}+A_{22}A_{22}^{}\right).`$
Note that the tensor coefficients are in general matrices in colour space and the colour sum has still to be performed in evaluating, for example, $`A_{11}A_{11}^{}`$. In CDR, additional terms of order $`ϵ`$ and $`ϵ^2`$ are straightforwardly generated. The resulting formula is, however, rather lengthy and we do not give it here.

## 3 The perturbative expansion of the tensor coefficients

The three processes that we wish to consider in detail are,

$`0`$ $``$ $`q\overline{q}gg,`$ (19)
$`0`$ $``$ $`q\overline{q}g\gamma ,`$ (20)
$`0`$ $``$ $`q\overline{q}\gamma \gamma .`$ (21)
The amplitude for the process $`𝒫`$ (where $`𝒫=gg,g\gamma `$ and $`\gamma \gamma `$) has the perturbative expansion of the form,

$`|^𝒫`$ $`=`$ $`c^𝒫\left[|^{𝒫,(0)}+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)|^{𝒫,(1)}+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)^2|^{𝒫,(2)}+𝒪(\alpha _s^3)\right],`$ (22)
where the overall coupling factor is given by

$$c^𝒫=4\pi \alpha _s^q(\alpha e_q^2)^{1q},$$
(23)
with $`q=1,1/2,0`$ for $`𝒫=gg,g\gamma ,\gamma \gamma `$ respectively. $`|^{𝒫,(n)}`$ denotes the $`n`$-loop contribution to the amplitude. Denoting the tensor coefficients as $`X^𝒫`$ (for $`X=A_{IJ}`$ and $`B`$), we find the perturbative expansions of the form

$`X_{IJ}^𝒫`$ $`=`$ $`c^𝒫\left[X_{IJ}^{𝒫,(0)}+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)X_{IJ}^{𝒫,(1)}+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)^2X_{IJ}^{𝒫,(2)}+𝒪(\alpha _s^3)\right],`$ (24)
where the dependence on the kinematic scales is implicit. The $`n`$-loop perturbative coefficients are vectors in colour space and can be further decomposed as,

$`X^{gg,(n)}`$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{3}{}}}𝒞^{gg,[i]}X^{gg,(n),[i]},`$
$`X^{g\gamma ,(n)}`$ $`=`$ $`𝒞^{g\gamma ,[1]}X^{g\gamma ,(n),[1]},`$
$`𝒞^{\gamma \gamma ,[1]}`$ $`=`$ $`\delta _{i_1i_2}X^{\gamma \gamma ,(n),[1]},`$ (25)
where

$$𝒞^{gg,[1]}=\left(𝑻^{a_3}𝑻^{a_4}\right)_{i_2i_1},𝒞^{gg,[2]}=\left(𝑻^{a_4}𝑻^{a_3}\right)_{i_2i_1},𝒞^{gg,[3]}=𝜹^{a_3a_4}\delta _{i_1i_2},$$
(26)
and

$$𝒞^{g\gamma ,[1]}=\left(𝑻^{a_3}\right)_{i_2i_1},𝒞^{\gamma \gamma ,[1]}=\delta _{i_1i_2}.$$
(27)
Index $`i`$ denotes a colour label in the fundamental representation while index $`a`$ is in the adjoint. Note that the normalisation of the colour matrices is given by $`Tr(𝑻^{a_1}𝑻^{a_2})=\frac{1}{2}𝜹^{a_1a_2}`$.

Symmetry under interchange of the two gluons gives the additional relations,

$`A_{IJ}^{gg,(n),[2]}(s_{13},s_{23})`$ $`=`$ $`A_{JI}^{gg,(n),[1]}(s_{23},s_{13}),`$
$`A_{IJ}^{gg,(n),[3]}(s_{13},s_{23})`$ $`=`$ $`A_{JI}^{gg,(n),[3]}(s_{23},s_{13}),`$
$`B^{gg,(n),[2]}(s_{13},s_{23})`$ $`=`$ $`B^{gg,(n),[1]}(s_{23},s_{13}),`$
$`B^{gg,(n),[3]}(s_{13},s_{23})`$ $`=`$ $`B^{gg,(n),[3]}(s_{23},s_{13}).`$ (28)
For the processes involving one or more photons,

$`A_{IJ}^{g\gamma ,(n),[1]}(s_{13},s_{23})`$ $`=`$ $`A_{JI}^{g\gamma ,(n),[1]}(s_{23},s_{13}),`$
$`A_{IJ}^{\gamma \gamma ,(n),[1]}(s_{13},s_{23})`$ $`=`$ $`A_{JI}^{\gamma \gamma ,(n),[1]}(s_{23},s_{13}),`$
$`B^{g\gamma ,(n),[1]}(s_{13},s_{23})`$ $`=`$ $`B^{g\gamma ,(n),[1]}(s_{23},s_{13}),`$
$`B^{\gamma \gamma ,(n),[1]}(s_{13},s_{23})`$ $`=`$ $`B^{\gamma \gamma ,(n),[1]}(s_{23},s_{13}).`$ (29)
As in previous work, we use QGRAF  to generate the Feynman diagrams and then use MAPLE  and FORM  to manipulate the Dirac structures and traceology and then reduce the scalar loop integrals to master integrals.

### 3.1 Tree level results

At tree level, none of the $`T_{IJ}`$ structures are present so that,

$$A_{IJ}^{𝒫,(0),[i]}(s_{13},s_{23})=0,$$
(30)
while only the first two colour structures contribute to $`B^{gg,(0)}`$,

$`B^{gg,(0),[1]}(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{2}{s_{12}s_{23}}},`$
$`B^{gg,(0),[2]}(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{2}{s_{12}s_{13}}},`$
$`B^{gg,(0),[3]}(s_{13},s_{23})`$ $`=`$ $`0.`$ (31)
For the other processes,

$`B^{g\gamma ,(0),[1]}(s_{13},s_{23})=B^{\gamma \gamma ,(0),[1]}(s_{13},s_{23})`$ $`=`$ $`B^{gg,(0),[1]}(s_{13},s_{23})+B^{gg,(0),[2]}(s_{13},s_{23})`$ (32)
$`=`$ $`{\displaystyle \frac{2}{s_{13}s_{23}}}.`$
### 3.2 One-loop results

At one-loop, all of the tensor and colour structures contribute. For the two-gluon process, we find the following (unrenormalised) structure,

$`A_{11}^{gg,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_Ff_1(s_{13},s_{23})+C_Ag_1(s_{13},s_{23})+(2N_FC_A(d2))h(s_{12}),`$
$`A_{11}^{gg,(1),[2]}(s_{13},s_{23})`$ $`=`$ $`A_{11}^{gg,(1),[1]}(s_{23},s_{13}),`$
$`A_{11}^{gg,(1),[3]}(s_{13},s_{23})`$ $`=`$ $`0,`$
$`A_{12}^{gg,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_Ff_2(s_{13},s_{23})+C_Ag_2(s_{13},s_{23})+(2N_FC_A(d2))h(s_{12}),`$
$`A_{12}^{gg,(1),[2]}(s_{13},s_{23})`$ $`=`$ $`C_Ff_3(s_{13},s_{23})+C_Ag_3(s_{13},s_{23})(2N_FC_A(d2))h(s_{12}),`$
$`A_{12}^{gg,(1),[3]}(s_{13},s_{23})`$ $`=`$ $`h_1(s_{13},s_{23}),`$
$`B^{gg,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_Ff_4(s_{13},s_{23})+C_Ag_4(s_{13},s_{23}),`$
$`B^{gg,(1),[2]}(s_{13},s_{23})`$ $`=`$ $`B^{gg,(1),[1]}(s_{23},s_{13}),`$
$`B^{gg,(1),[3]}(s_{13},s_{23})`$ $`=`$ $`h_2(s_{13},s_{23}),`$ (33)
while for the gluon-photon and photon-photon processes we have,

$`A_{11}^{g\gamma ,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_F\left(f_1(s_{13},s_{23})f_1(s_{23},s_{13})\right)`$
$`+{\displaystyle \frac{C_A}{2}}\left(g_1(s_{13},s_{23})g_1(s_{23},s_{13})4h_3(s_{13},s_{23})\right),`$
$`A_{12}^{g\gamma ,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_F\left(f_2(s_{13},s_{23})+f_3(s_{13},s_{23})\right)`$
$`+{\displaystyle \frac{C_A}{2}}\left(g_2(s_{13},s_{23})+g_3(s_{13},s_{23})2h_1(s_{13},s_{23})\right),`$
$`B^{g\gamma ,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_F\left(f_4(s_{13},s_{23})+f_4(s_{23},s_{13})\right)`$
$`+{\displaystyle \frac{C_A}{2}}\left(g_4(s_{13},s_{23})+g_4(s_{23},s_{13})2h_2(s_{13},s_{23})2h_3(s_{13},s_{23})\right),`$
and,

$`A_{11}^{\gamma \gamma ,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_F\left(f_1(s_{13},s_{23})f_1(s_{23},s_{13})\right),`$
$`A_{12}^{\gamma \gamma ,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_F\left(f_2(s_{13},s_{23})+f_3(s_{13},s_{23})\right),`$
$`B^{\gamma \gamma ,(1),[1]}(s_{13},s_{23})`$ $`=`$ $`C_F\left(f_4(s_{13},s_{23})+f_4(s_{23},s_{13})\right).`$ (35)
The functions $`f_i`$, $`g_i`$, $`h_i`$ and $`h`$ can be written in terms of the one-loop box integral in $`d=62ϵ`$ dimensions, Box$`{}_{}{}^{6}(s_{ij},s_{ik})`$, and the one-loop bubble, Bub$`(s_{ij})`$. For all orders in $`ϵ=(4d)/2`$, we find that,

$$h(s_{12})=\frac{2(d4)}{s_{12}^2(d1)(d2)}\mathrm{Bub}(s_{12}),$$
(36)
while,

$`f_1(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{(d2)(d4)}{s_{13}}}\mathrm{Box}^6(s_{12},s_{23}){\displaystyle \frac{2(d4)}{s_{13}s_{23}}}\mathrm{Bub}(s_{23}){\displaystyle \frac{2(d4)}{s_{12}s_{13}}}\mathrm{Bub}(s_{12}),`$
$`f_2(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{s_{23}^2d(d2)s_{13}s_{23}(d2)(d14)2s_{13}^2(d10)}{s_{12}s_{13}^2}}\mathrm{Box}^6(s_{12},s_{23})`$
$`+{\displaystyle \frac{2(d12)s_{13}2s_{23}d}{s_{12}s_{13}^2}}\mathrm{Bub}(s_{23})+{\displaystyle \frac{2s_{23}d4(d8)s_{13}}{s_{12}s_{13}^2}}\mathrm{Bub}(s_{12}),`$
$`f_3(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{(d4)(d10)s_{23}^2+(3d^236d+92)s_{13}s_{23}2(d2)s_{13}^2}{s_{12}s_{13}s_{23}}}\mathrm{Box}^6(s_{12},s_{13})`$
$`+{\displaystyle \frac{4s_{13}4(d7)s_{23}}{s_{12}s_{13}s_{23}}}\mathrm{Bub}(s_{13})+{\displaystyle \frac{2(d10)s_{23}4s_{13}}{s_{12}s_{13}s_{23}}}\mathrm{Bub}(s_{12}),`$
$`f_4(s_{13},s_{23})`$ $`=`$ $`+{\displaystyle \frac{4(d3)s_{13}^2(d^215d+38)s_{23}s_{13}+(d2)s_{23}^2}{s_{12}s_{13}s_{23}}}\mathrm{Box}^6(s_{12},s_{23})`$
$`+{\displaystyle \frac{2(d7)s_{13}2s_{23}}{s_{12}s_{13}s_{23}}}\mathrm{Bub}(s_{23})+{\displaystyle \frac{2(d4)s_{23}+8(d3)s_{13}}{(d4)s_{12}s_{13}s_{23}}}\mathrm{Bub}(s_{12}),`$
$`g_1(s_{13},s_{23})`$ $`=`$ $`f_1(s_{13},s_{23}),`$
$`g_2(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{s_{23}(d2)(4s_{13}+ds_{23})}{s_{12}s_{13}^2}}\mathrm{Box}^6(s_{12},s_{23})+{\displaystyle \frac{2s_{23}d+4s_{13}}{s_{12}s_{13}^2}}\mathrm{Bub}(s_{23})`$
$`+{\displaystyle \frac{2(d6)s_{13}2ds_{23}}{s_{12}s_{13}^2}}\mathrm{Bub}(s_{12}),`$
$`g_3(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{(d2)(2s_{13}(d6)s_{23})}{s_{12}s_{23}}}\mathrm{Box}^6(s_{12},s_{13})+{\displaystyle \frac{4}{s_{12}s_{23}}}\mathrm{Bub}(s_{12})`$
$`+{\displaystyle \frac{2(d4)s_{23}4s_{13}}{s_{12}s_{13}s_{23}}}\mathrm{Bub}(s_{13}),`$
$`g_4(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{(d^216d+36)s_{13}2(d2)s_{23}}{2s_{12}s_{13}}}\mathrm{Box}^6(s_{12},s_{23})`$
$`{\displaystyle \frac{(d^216d+40)s_{13}2(d4)s_{23}}{(d4)s_{12}s_{13}s_{23}}}\mathrm{Bub}(s_{23}){\displaystyle \frac{2}{s_{12}s_{13}}}\mathrm{Bub}(s_{12}),`$
$`h_1(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{(d10)((d2)s_{23}+2s_{13})}{2s_{12}s_{13}}}\mathrm{Box}^6(s_{12},s_{23})`$
$`+{\displaystyle \frac{(d4)(d10)(2s_{13}+s_{23})}{2s_{12}s_{13}}}\mathrm{Box}^6(s_{12},s_{13})`$
$`{\displaystyle \frac{(d10)(2s_{13}(d4)s_{23})}{2s_{12}s_{13}}}\mathrm{Box}^6(s_{13},s_{23}),`$
$`h_2(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{(8(d3)s_{13}(d4)(d10)s_{23})}{4s_{12}s_{23}}}(\mathrm{Box}^6(s_{12},s_{23})+\mathrm{Box}^6(s_{13},s_{23}))`$
$`+{\displaystyle \frac{(8(d3)s_{23}(d4)(d10)s_{13})}{4s_{12}s_{13}}}(\mathrm{Box}^6(s_{12},s_{13})+\mathrm{Box}^6(s_{13},s_{23}))`$
$`+{\displaystyle \frac{4+3(d4)^2}{4s_{12}}}\mathrm{Box}^6(s_{13},s_{23})`$
$`{\displaystyle \frac{4(d3)}{(d4)}}\left({\displaystyle \frac{1}{s_{12}s_{13}}}\mathrm{Bub}(s_{23})+{\displaystyle \frac{1}{s_{12}s_{23}}}\mathrm{Bub}(s_{13})+{\displaystyle \frac{1}{s_{23}s_{13}}}\mathrm{Bub}(s_{12})\right),`$
$`h_3(s_{13},s_{23})`$ $`=`$ $`{\displaystyle \frac{1}{4}}f_1(s_{12},s_{23}).`$ (37)
### 3.3 Ultraviolet renormalisation

The renormalisation of the matrix element is carried out by replacing the bare coupling $`\alpha _0`$ with the renormalised coupling $`\alpha _s\alpha _s(\mu ^2)`$, evaluated at the renormalisation scale $`\mu ^2`$ according to Eq. (4).

We denote the $`i`$-loop contribution to the unrenormalised tensor coefficients for process $`𝒫`$ by $`A_{IJ}^{𝒫,(i),\mathrm{un}}`$ and $`B^{𝒫,(i),\mathrm{un}}`$, using the same normalization as for the decomposition of the tensor coefficients (22); the dependence on $`(s_{13},s_{23})`$ is always understood implicitly. The renormalised coefficients are then obtained as

$`A_{IJ}^{𝒫,(0)}`$ $`=`$ $`0,`$
$`A_{IJ}^{𝒫,(1)}`$ $`=`$ $`S_ϵ^1A_{IJ}^{𝒫,(1),\mathrm{un}},`$
$`A_{IJ}^{𝒫,(2)}`$ $`=`$ $`S_ϵ^2A_{IJ}^{𝒫,(2),\mathrm{un}}{\displaystyle \frac{(q+1)\beta _0}{ϵ}}S_ϵ^1A_{IJ}^{𝒫,(1),\mathrm{un}},`$ (38)
and

$`B^{𝒫,(0)}`$ $`=`$ $`B^{𝒫,(0),\mathrm{un}},`$
$`B^{𝒫,(1)}`$ $`=`$ $`S_ϵ^1B^{𝒫,(1),\mathrm{un}}{\displaystyle \frac{q\beta _0}{ϵ}}B^{𝒫,(0),\mathrm{un}},`$
$`B^{𝒫,(2)}`$ $`=`$ $`S_ϵ^2B^{𝒫,(2),\mathrm{un}}{\displaystyle \frac{(q+1)\beta _0}{ϵ}}S_ϵ^1B^{𝒫,(1),\mathrm{un}}{\displaystyle \frac{q}{2}}\left({\displaystyle \frac{\beta _1}{ϵ}}{\displaystyle \frac{(q+1)\beta _0^2}{ϵ^2}}\right)B^{𝒫,(0),\mathrm{un}},`$
where $`q`$ is the number of powers of the strong coupling for process $`𝒫`$ as described in Eq. (23).

### 3.4 Infrared behaviour of the tensor coefficients

After performing ultraviolet renormalisation, the amplitudes still contain singularities which are of infrared origin and will be analytically cancelled by those occurring in radiative processes of the same order. Catani  has shown how to organize the infrared pole structure of the one- and two-loop contributions renormalised in the $`\overline{\mathrm{MS}}`$ scheme in terms of the tree and renormalised one-loop amplitudes. The same procedure applies to the tensor coefficients.

In particular, the infrared behaviour of the one-loop coefficients is given by

$`A_{IJ}^{𝒫,(1)}`$ $`=`$ $`A_{IJ}^{𝒫,(1),\mathrm{finite}},`$
$`B^{𝒫,(1)}`$ $`=`$ $`𝑰^{𝒫,(1)}(ϵ)B^{𝒫,(0)}+B^{𝒫,(1),\mathrm{finite}},`$ (40)
while the two-loop singularity structure is

$`A_{IJ}^{𝒫,(2)}`$ $`=`$ $`𝑰^{𝒫,(1)}(ϵ)A_{IJ}^{𝒫,(1)}+A_{IJ}^{𝒫,(2),\mathrm{finite}},`$
$`B^{𝒫,(2)}`$ $`=`$ $`𝑰^{𝒫,(2)}(ϵ)B^{𝒫,(0)}+𝑰^{𝒫,(1)}(ϵ)B^{(1)}+B^{(2),\mathrm{finite}},`$ (41)
where $`𝑰^{𝒫,(2)}(ϵ)`$ is given by,

$`𝑰^{𝒫,(2)}(ϵ)`$ $`=`$ $`{\displaystyle \frac{1}{2}}𝑰^{𝒫,(1)}(ϵ)𝑰^{𝒫,(1)}(ϵ){\displaystyle \frac{\beta _0}{ϵ}}𝑰^{𝒫,(1)}(ϵ)+e^{ϵ\gamma }{\displaystyle \frac{\mathrm{\Gamma }(12ϵ)}{\mathrm{\Gamma }(1ϵ)}}\left({\displaystyle \frac{\beta _0}{ϵ}}+K\right)𝑰^{𝒫,(1)}(2ϵ)`$ (42)
$`+𝑯^{𝒫,(2)}(ϵ)`$
and the constant $`K`$ is

$$K=\left(\frac{67}{18}\frac{\pi ^2}{6}\right)C_A\frac{10}{9}T_RN_F.$$
(43)
In QCD the $`𝑰^{𝒫,(1)}(ϵ)`$ operator is given by

$`𝑰^{𝒫,(1)}(ϵ)`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{e^{ϵ\gamma }}{\mathrm{\Gamma }(1ϵ)}}{\displaystyle \underset{i}{}}\nu _i^{sing}(ϵ){\displaystyle \underset{ji}{}}𝐓_i𝐓_j\left({\displaystyle \frac{\mu ^2}{2p_ip_j}}\right)^ϵ`$ (44)
where the sum runs over the pairs of external coloured particles involved in process $`𝒫`$ and the singular function is

$$\nu _i^{sing}(ϵ)=\frac{1}{ϵ^2}+\gamma _i\frac{1}{ϵ},\gamma _{q,\overline{q}}=\frac{3}{2},\gamma _g=\frac{\beta _0}{C_A}$$
(45)
The colour charge algebra is defined as

$$𝐓_i𝐓_j=\{\begin{array}{cc}𝐓_j𝐓_i\hfill & \text{if }ij\text{,}\hfill \\ 𝐓_i^2=C_i\hfill & \text{otherwise,}\hfill \end{array}$$

with the Casimir $`C_i=C_F`$ ($`C_i=C_A`$) if parton $`i`$ is a quark (gluon) and where the colour charge $`𝐓_i`$ is $`t_{cb}^a`$ ($`t_{cb}^a`$) for a quark (anti-quark) and $`if_{cab}`$ for a gluon.

Then, given the colour basis for the two-quark two-gluon process presented in Eq. (3), we have the following infrared-singularity operator

$`𝑰^{gg,(1)}(ϵ)={\displaystyle \frac{e^{ϵ\gamma }}{\mathrm{\Gamma }(1ϵ)}}\left(\begin{array}{ccc}𝑨(ϵ,s_{12},s_{13},s_{23})& 0& 𝑩(ϵ,s_{23},s_{13},s_{12})\\ 0& 𝑨(ϵ,s_{12},s_{23},s_{13})& 𝑩(ϵ,s_{13},s_{23},s_{12})\\ 𝑩(ϵ,s_{12},s_{13},s_{23})& 𝑩(ϵ,s_{12},s_{23},s_{13})& 𝑪(ϵ,s_{12})\end{array}\right)`$ (49)
where

$`𝑨(ϵ,s_{12},s_{23},s_{13})`$ $`=`$ $`\left({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{2ϵ}}\right)\left[N\left({\displaystyle \frac{\mu ^2}{s_{13}}}\right)^ϵ+C_F\left({\displaystyle \frac{\mu ^2}{s_{12}}}\right)^ϵ\right]`$ (51)
$`+N\left({\displaystyle \frac{\beta _0}{2Nϵ}}{\displaystyle \frac{3}{4ϵ}}\right)\left[\left({\displaystyle \frac{\mu ^2}{s_{12}}}\right)^ϵ+\left({\displaystyle \frac{\mu ^2}{s_{13}}}\right)^ϵ\right],`$
$`𝑩(ϵ,s_{12},s_{23},s_{13})`$ $`=`$ $`\left({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{4ϵ}}+{\displaystyle \frac{\beta _0}{2Nϵ}}\right)\left[\left({\displaystyle \frac{\mu ^2}{s_{23}}}\right)^ϵ\left({\displaystyle \frac{\mu ^2}{s_{12}}}\right)^ϵ\right],`$ (52)
$`𝑪(ϵ,s_{12})`$ $`=`$ $`\left[\left({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{2ϵ}}\right)C_F+\left({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{\beta _0}{Nϵ}}\right)N\right]\left({\displaystyle \frac{\mu ^2}{s_{12}}}\right)^ϵ.`$ (53)
Similarly, for the other processes we have

$`𝑰^{g\gamma ,(1)}(ϵ)`$ $`=`$ $`{\displaystyle \frac{e^{ϵ\gamma }}{2\mathrm{\Gamma }(1ϵ)}}\{{\displaystyle \frac{1}{N}}({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{2ϵ}})({\displaystyle \frac{\mu ^2}{s_{12}}})^ϵ`$ (54)
$`N({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{4ϵ}}+{\displaystyle \frac{\beta _0}{2Nϵ}})[({\displaystyle \frac{\mu ^2}{s_{23}}})^ϵ+({\displaystyle \frac{\mu ^2}{s_{13}}})^ϵ]\}`$
and

$`𝑰^{\gamma \gamma ,(1)}(ϵ)`$ $`=`$ $`C_F{\displaystyle \frac{e^{ϵ\gamma }}{\mathrm{\Gamma }(1ϵ)}}\left({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{2ϵ}}\right)\left({\displaystyle \frac{\mu ^2}{s_{12}}}\right)^ϵ.`$ (55)
In Eq. (3.4) the function $`𝑯^{𝒫,(2)}`$ contains poles of $`𝒪\left(1/ϵ\right)`$ and is process and renormalisation scheme dependent. The specific structure of this term was not given in Ref. . However, there it was pointed out that in the case of a QCD process that involves a $`q\overline{q}`$ pair, the function $`𝑯^{𝒫,(2)}`$ is related to the quark electromagnetic form factor. Now, we know that the form factor is not merely related, but is the backbone of the singular structure of multi-loop QCD amplitudes .

It was only after the direct Feynman diagram evaluation of two-loop matrix elements , that the colour uncorrelated part of the $`𝑯^{𝒫,(2)}`$ function for all $`22`$ partonic processes was fully identified. More precisely, it was found that each external coloured leg in the partonic process contributes independently to the matrix element given by

$`_0|𝑯^{𝒫,(2)}|_0={\displaystyle \frac{e^{ϵ\gamma }}{4ϵ\mathrm{\Gamma }(1ϵ)}}H^{𝒫,(2)}_0|_0`$ (56)
with

$`H^{𝒫,(2)}=n_qH_q^{(2)}+n_gH_g^{(2)}`$ (57)
where $`n_g`$ is the number of external gluons, $`n_q`$ is the number of external quarks and anti-quarks for the process $`𝒫`$ and

$`H_q^{(2)}=H_{\overline{q}}^{(2)}`$ $`=`$ $`\left({\displaystyle \frac{\pi ^2}{2}}6\zeta _3{\displaystyle \frac{3}{8}}\right)C_F^2+\left({\displaystyle \frac{13}{2}}\zeta _3+{\displaystyle \frac{245}{216}}{\displaystyle \frac{23}{48}}\pi ^2\right)C_AC_F`$ (58)
$`+\left({\displaystyle \frac{25}{54}}+{\displaystyle \frac{\pi ^2}{12}}\right)T_RN_FC_F,`$
$`H_g^{(2)}`$ $`=`$ $`{\displaystyle \frac{20}{27}}T_R^2N_F^2+T_RC_FN_F\left({\displaystyle \frac{\pi ^2}{36}}+{\displaystyle \frac{58}{27}}\right)T_RN_FC_A`$ (59)
$`+\left({\displaystyle \frac{\zeta _3}{2}}+{\displaystyle \frac{5}{12}}+{\displaystyle \frac{11}{144}}\pi ^2\right)C_A^2.`$
Already in Ref. , Catani pointed out the possibility of $`𝑯^{𝒫,(2)}`$ having additional non-trivial colour correlations which could be explored by calculating the hard scattering amplitude, rather than the interference of two-loop and tree amplitudes.<sup>1</sup><sup>1</sup>1The sum over colours causes the colour correlations to vanish Such colour correlations were found for the four gluon interaction, which is the only partonic process (involving four coloured particles) for which the two-loop helicity amplitudes  exist. According to Ref. , the colour correlations form a non-trivial matrix in colour space and is the same for all helicities.

For the processes we study, we find that the colour correlations are indeed helicity independent and are given by

$`𝑯^{𝒫,(2)}(ϵ)={\displaystyle \frac{e^{ϵ\gamma }}{4ϵ\mathrm{\Gamma }(1ϵ)}}\left({\displaystyle \frac{\mu ^2}{s_{12}}}\right)^{2ϵ}\left(H^{𝒫,(2)}\mathrm{𝟏}+\widehat{𝐇}^{𝒫,(2)}\right)`$ (60)
with $`H^{𝒫,(2)}`$ as given in Eq. (57) and

$`\widehat{𝐇}^{gg,(2)}`$ $`=`$ $`4if_{abc}𝑻^{a[1]}𝑻^{b[2]}𝑻^{c[3]}\mathrm{ln}\left({\displaystyle \frac{s_{12}}{s_{23}}}\right)\mathrm{ln}\left({\displaystyle \frac{s_{23}}{s_{13}}}\right)\mathrm{ln}\left({\displaystyle \frac{s_{13}}{s_{12}}}\right),`$ (61)
$`\widehat{𝐇}^{g\gamma ,(2)}`$ $`=`$ $`0,`$ (62)
$`\widehat{𝐇}^{\gamma \gamma ,(2)}`$ $`=`$ $`0.`$ (63)
These results are in complete agreement with the results found in Ref.  and show that the colour correlations are also independent of the external particles. This also agrees with the analysis of Ref.  where it was shown that the contributions to $`\widehat{𝐇}^{𝒫,(2)}`$ are of soft nature and directly related to soft anomalous dimension matrices obtained by studying the evolution of colour exchange in QCD hard scattering (see  and references therein). It is worth mentioning that these type of colour correlations between three different partons also appeared in the analysis of higher-order contributions to the soft-gluon current in Ref. .

Note that the matrix $`f_{abc}𝑻^{a[1]}𝑻^{b[2]}𝑻^{c[3]}`$ collects all colour weights that result from the analysis of triple gluon vertices connecting any three different external partons, projected into the colour basis given in Eq. (3), so that for the two quark-two gluon process,

$`if_{abc}𝑻^{a[1]}𝑻^{b[2]}𝑻^{c[3]}={\displaystyle \frac{1}{8}}\left(\begin{array}{ccc}2& 2& 4N\\ 2& 2& 4N\\ N& N& 0\end{array}\right).`$ (67)
Given the tree level results presented in Section 3.1 and the divergent structure presented in Eq.(3.4), we can see that in fact the third column of this matrix is not needed at $`𝒪\left(1/ϵ\right)`$, since the $`𝜹^{a_3a_4}\delta _{i_1i_2}`$ colour structure is not present at tree level.

## 4 Helicity amplitudes

The helicity amplitudes $`|_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ can be obtained from the general $`d`$-dimensional tensor of Eq. (10) using standard four-dimensional helicity techniques corresponding to working in the ’t Hooft-Veltman scheme. Conventionally, the two helicity states of a four-dimensional light-like spinor $`\psi (p)`$ are denoted by,

$$\psi _\pm (p)=\frac{1}{2}(1\pm \gamma _5)\psi (p),$$
(68)
with the further notation,

$$|p\pm =\psi _\pm (p),p\pm |=\overline{\psi _\pm (p)}.$$
(69)
Particles may thus be crossed to the initial state by reversing the sign of the helicity. The basic quantity is the spinor product,

$$pq=p|q+,[pq]=[p+|q],$$
(70)
such that

$$pq[qp]=2p.q.$$
(71)
In this notation, the polarization vector of the out-going light-like gauge boson of momentum $`p`$ is defined with respect to a reference momentum $`k`$ as,

$$ϵ^\pm (p,k)=\pm \frac{p\pm |\gamma _\mu |k\pm }{\sqrt{2}k|p\pm }.$$
(72)
Fixing the helicity of the (final-state) quark to be positive, $`\lambda _1=+`$ we find that the helicity amplitudes, upto a helicity dependent phase, are given by

$`|_{+++}`$ $``$ $`{\displaystyle \frac{\sqrt{s_{12}^2s_{13}s_{23}}}{2}}A_{11}(s_{13},s_{23}),`$
$`|_{++}`$ $``$ $`{\displaystyle \frac{\sqrt{s_{13}^3s_{23}}}{2}}\left(2B(s_{13},s_{23})+A_{12}(s_{13},s_{23})A_{11}(s_{13},s_{23})\right),`$
$`|_{++}`$ $``$ $`{\displaystyle \frac{\sqrt{s_{13}s_{23}^3}}{2}}\left(2B(s_{13},s_{23})+A_{22}(s_{13},s_{23})A_{21}(s_{13},s_{23})\right),`$
$`|_+`$ $``$ $`{\displaystyle \frac{\sqrt{s_{12}^2s_{13}s_{23}}}{2}}A_{22}(s_{13},s_{23}).`$ (73)
By inspection of Eqs. (3) and (3) we see that, as expected, there are two independent helicity amplitudes and that the other two helicity amplitudes can be obtained by Bose symmetry and,

$`|_{++}`$ $`=`$ $`|_{++}(s_{13}s_{23}),`$
$`|_+`$ $`=`$ $`|_{+++}(s_{13}s_{23}).`$ (74)
Helicity amplitudes where the quark helicity is negative are obtained through parity,

$$|_{+\lambda _3\lambda _4}=\left(|_{+\lambda _3\lambda _4}\right)^{}.$$
(75)
By inspection, we see that “squaring” and summing over the helicities generates the interference of amplitudes given in Eq. (18). Thus, the helicity amplitudes can be obtained directly from a Feynman diagram calculation by projecting out the individual tensor coefficients and taking the appropriate linear combination.

Just as for the individual tensor coefficients, the helicity amplitudes for each process $`𝒫`$ can be perturbatively decomposed as,

$`|^𝒫_{\lambda _1\lambda _2\lambda _3\lambda _4}`$
$`=`$ $`c^𝒫\left[|^{𝒫,(0)}_{\lambda _1\lambda _2\lambda _3\lambda _4}+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)|^{𝒫,(1)}_{\lambda _1\lambda _2\lambda _3\lambda _4}+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)^2|^{𝒫,(2)}_{\lambda _1\lambda _2\lambda _3\lambda _4}+𝒪(\alpha _s^3)\right],`$
where the overall coupling factor is given by Eq. (23). As in Eqs. (3.3) and (3.3), the renormalised helicity amplitudes are obtained from the unrenormalised amplitudes by

$`|^{𝒫,(0)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ $`=`$ $`|^{𝒫,(0),\mathrm{un}}_{\lambda _1\lambda _2\lambda _3\lambda _4},`$
$`|^{𝒫,(1)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ $`=`$ $`S_ϵ^1|^{𝒫,(1),\mathrm{un}}_{\lambda _1\lambda _2\lambda _3\lambda _4}{\displaystyle \frac{q\beta _0}{ϵ}}|^{𝒫,(0),\mathrm{un}}_{\lambda _1\lambda _2\lambda _3\lambda _4},`$
$`|^{𝒫,(2)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ $`=`$ $`S_ϵ^2|^{𝒫,(2),\mathrm{un}}_{\lambda _1\lambda _2\lambda _3\lambda _4}{\displaystyle \frac{(q+1)\beta _0}{ϵ}}S_ϵ^1|^{𝒫,(1),\mathrm{un}}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$
$`{\displaystyle \frac{q}{2}}\left({\displaystyle \frac{\beta _1}{ϵ}}{\displaystyle \frac{(q+1)\beta _0^2}{ϵ^2}}\right)|^{𝒫,(0),\mathrm{un}}_{\lambda _1\lambda _2\lambda _3\lambda _4},`$
while the infrared singularity structure is given by,

$`|^{𝒫,(1)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ $`=`$ $`𝑰^{𝒫,(1)}(ϵ)|^{𝒫,(0)}_{\lambda _1\lambda _2\lambda _3\lambda _4}+|^{𝒫,(1),\mathrm{finite}}_{\lambda _1\lambda _2\lambda _3\lambda _4},`$ (78)
$`|^{𝒫,(2)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ $`=`$ $`𝑰^{𝒫,(2)}(ϵ)|^{𝒫,(0)}_{\lambda _1\lambda _2\lambda _3\lambda _4}+𝑰^{𝒫,(1)}(ϵ)|^{𝒫,(1)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ (79)
$`+|^{𝒫,(2),\mathrm{finite}}_{\lambda _1\lambda _2\lambda _3\lambda _4},`$
where the operators $`𝑰^{𝒫,(i)}`$ are the process dependent matrices in colour space given in Section 3.4.

Similarly, the $`n`$-loop helicity amplitudes are vectors in colour space and can be further decomposed as in Eq. (3)

$`|^{gg,(n)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{3}{}}}𝒞^{gg,[i]}|^{gg,(n),[i]}_{\lambda _1\lambda _2\lambda _3\lambda _4},`$
$`|^{g\gamma ,(n)}`$ $`=`$ $`𝒞^{g\gamma ,[1]}|^{g\gamma ,(n),[1]}_{\lambda _1\lambda _2\lambda _3\lambda _4},`$
$`|^{\gamma \gamma ,(n)}`$ $`=`$ $`𝒞^{\gamma \gamma ,[1]}|^{\gamma \gamma ,(n),[1]}_{\lambda _1\lambda _2\lambda _3\lambda _4}.`$ (80)
### 4.1 Helicity amplitudes for physical processes

The physically relevant scattering amplitudes we wish to describe are the quark-antiquark annihilation processes,

$`q(p_2,+)+\overline{q}(p_1,)`$ $``$ $`b_3(p_3,+)+b_4(p_4,+),`$ (81)
$`q(p_2,+)+\overline{q}(p_1,)`$ $``$ $`b_3(p_3,+)+b_4(p_4,),`$ (82)
$`q(p_2,+)+\overline{q}(p_1,)`$ $``$ $`b_3(p_3,)+b_4(p_4,+),`$ (83)
$`q(p_2,+)+\overline{q}(p_1,)`$ $``$ $`b_3(p_3,)+b_4(p_4,),`$ (84)
and the quark-boson scattering processes,

$`q(p_2,+)+b_3(p_1,+)`$ $``$ $`q(p_3,+)+b_4(p_4,+),`$ (85)
$`q(p_2,+)+b_3(p_1,+)`$ $``$ $`q(p_3,+)+b_4(p_4,),`$ (86)
$`q(p_2,+)+b_3(p_1,)`$ $``$ $`q(p_3,+)+b_4(p_4,+),`$ (87)
$`q(p_2,+)+b_3(p_1,)`$ $``$ $`q(p_3,+)+b_4(p_4,),`$ (88)
$`q(p_2,+)+b_4(p_1,+)`$ $``$ $`b_3(p_3,+)+q(p_4,+),`$ (89)
$`q(p_2,+)+b_4(p_1,+)`$ $``$ $`b_3(p_3,)+q(p_4,+),`$ (90)
$`q(p_2,+)+b_4(p_1,)`$ $``$ $`b_3(p_3,+)+q(p_4,+),`$ (91)
$`q(p_2,+)+b_4(p_1,)`$ $``$ $`b_3(p_3,)+q(p_4,+).`$ (92)
where the bosons $`b_3`$ or $`b_4`$ are gluons or photons. Processes where the initial state quark has negative helicity are obtained by a parity transformation, while processes with two bosons in the initial state are given by time reversal. There are additional symmetries when the gauge bosons are identical and further symmetries between the various colour components for the two gluon process.

For convenience, we denote the processes of Eqs. (81)–(84) as the $`s`$-channel, the processes of Eqs. (85)–(88) as the $`u`$-channel and the processes of Eqs. (89)–(92) as the $`t`$-channel. $`s`$, $`t`$ and $`u`$ are the usual Mandelstam variables, $`s=s_{12}>0`$, $`t=s_{23}<0`$, $`u=s_{13}<0`$ and $`s+t+u=0`$.

To present the helicity amplitudes for the various processes, it is convenient to organise the amplitude where particles 1 and 2 are in the initial state in terms of a spinor factor $`𝒮`$ and the colour factor $`𝒞`$ so that,

$$|_c^{𝒫,(n)}_{\lambda _2\lambda _1\lambda _3\lambda _4}=\underset{i}{}𝒮_{c\lambda _2\lambda _1\lambda _3\lambda _4}^{𝒫,[i]}\times 𝒞_c^{𝒫,[i]}\times |\widehat{}_c^{𝒫,(n),[i]}_{\lambda _2\lambda _1\lambda _3\lambda _4},$$
(93)
where $`c`$ (= $`s`$, $`t`$, $`u`$) denotes the channel and the sum runs over the number of colour structures: three for the two-gluon processes and one otherwise.

Explicitly, we have

$`𝒞_s^{gg,[1]}`$ $`=`$ $`\left(𝑻^{a_3}𝑻^{a_4}\right)_{i_1i_2},`$
$`𝒞_s^{gg,[2]}`$ $`=`$ $`\left(𝑻^{a_4}𝑻^{a_3}\right)_{i_1i_2},`$
$`𝒞_s^{gg,[3]}`$ $`=`$ $`𝜹^{a_3a_4}\delta _{i_1i_2},`$ (94)
while for the processes with photons we have,

$`𝒞_s^{g\gamma ,[1]}`$ $`=`$ $`\left(𝑻^{a_3}\right)_{i_1i_2},`$
$`𝒞_s^{\gamma \gamma ,[1]}`$ $`=`$ $`\delta _{i_1i_2},`$ (95)
The $`u`$\- and $`t`$-channel colour factors are obtained by crossing symmetry,

$$𝒞_u^{𝒫,[i]}=𝒞_s^{𝒫,[i]}(13),𝒞_t^{𝒫,[i]}=𝒞_s^{𝒫,[i]}(14).$$
(96)
There is considerable freedom in the definition of the helicity dependent spinor factors. For the $`s`$-channel, for process (82) we choose,

$`𝒮_{s++}^{gg,[1]}`$ $`=`$ $`i{\displaystyle \frac{24^314}{12233441}},`$
$`𝒮_{s++}^{gg,[2]}`$ $`=`$ $`i{\displaystyle \frac{24^314}{12244331}},`$
$`𝒮_{s++}^{gg,[3]}=𝒮_{++}^{g\gamma }=𝒮_{++}^{\gamma \gamma }`$ $`=`$ $`i{\displaystyle \frac{24^314}{13142324}},`$ (97)
while for process (81) we choose,

$$𝒮_{s+++}^{𝒫,[i]}=i\frac{23[13]}{34^2}.$$
(98)
The factors for the other $`s`$-channel processes, (83,84) are obtained by Bose symmetry (exchanging labels 3 and 4). Note that for the two-gluon process, this symmetry also exchanges the first two colour structures, $`[1]`$ and $`[2]`$. The spinor prefactors for the $`u`$\- and $`t`$ channels are also obtained by crossing symmetry,

$$𝒮_{u\lambda _2\lambda _3\lambda _1\lambda _4}^{𝒫,[i]}=𝒮_{s\lambda _2\lambda _1\lambda _3\lambda _4}^{𝒫,[i]}(p_1p_3),𝒮_{t\lambda _2\lambda _4\lambda _3\lambda _1}^{𝒫,[i]}=𝒮_{s\lambda _2\lambda _1\lambda _3\lambda _4}^{𝒫,[i]}(p_1p_4).$$
(99)
### 4.2 Tree level helicity amplitudes

At tree-level, inserting the results for the tensor coefficients for the various processes given in Section 3.1, the only non-vanishing helicity amplitudes in the $`s`$-channel are,

$`|\widehat{}_s^{gg,(0),[1]}_{++}=|\widehat{}_s^{gg,(0),[1]}_{++}`$ $`=`$ $`2,`$
$`|\widehat{}_s^{gg,(0),[2]}_{++}=|\widehat{}_s^{gg,(0),[2]}_{++}`$ $`=`$ $`2,`$
$`|\widehat{}_s^{g\gamma ,(0)}_{++}=|\widehat{}^{g\gamma ,(0)}_{++}`$ $`=`$ $`2,`$
$`|\widehat{}_s^{\gamma \gamma ,(0)}_{++}=|\widehat{}^{\gamma \gamma ,(0)}_{++}`$ $`=`$ $`2.`$ (100)
The appropriate spinor factors are obtained from Eq. (4.1) while the colour factors are given in Eqs. (4.1) and (4.1). Tree helicity amplitudes for the other channels are obtained by the appropriate crossing symmetry.

### 4.3 One-loop helicity amplitudes

The one-loop helicity amplitudes are straightforwardly obtained by inserting the all-orders results for the tensor coefficients given in Section 3.2 into Eq. (4), expanding the one-loop bubble and box integrals around $`ϵ=0`$, and renormalising according to Eq. (4). Explicit expansions in $`ϵ`$ for the $`d=42ϵ`$ one-loop bubble and $`d=62ϵ`$ one-loop box graphs can be found, for example, in the Appendix of Ref. .

The finite remainders of the one-loop amplitudes defined through Eq. (78) for the two-gluon process can be decomposed according to the number of colours and massless quark flavours,

$`|\widehat{}_c^{gg,(1),[i],finite}`$ $`=`$ $`NA_c^{gg,(1),[i]}+{\displaystyle \frac{1}{N}}B_c^{gg,(1),[i]}+N_FC_c^{gg,(1),[i]}\beta _0S|\widehat{}_c^{gg,(0),[i]},`$
for $`i=1,2`$ and

$$|\widehat{}_c^{gg,(1),[3],finite}=D_c^{gg,(1),[3]}+\frac{N_F}{N}E_c^{gg,(1),[3]}$$
(102)
for $`i=3`$ and where

$$S=\mathrm{log}\left(\frac{s}{\mu ^2}\right)$$
(103)
and $`\beta _0`$ is defined in Eq. (6). For the other two processes we have,

$`|\widehat{}_c^{g\gamma ,(1),[1],finite}`$ $`=`$ $`NA_c^{g\gamma ,(1),[1]}+{\displaystyle \frac{1}{N}}B_c^{g\gamma ,(1),[1]}+N_FC_c^{g\gamma ,(1),[1]}{\displaystyle \frac{1}{2}}\beta _0S|\widehat{}_c^{g\gamma ,(0),[1]},`$
$`|\widehat{}_c^{\gamma \gamma ,(1),[1],finite}`$ $`=`$ $`NA_c^{\gamma \gamma ,(1),[1]}+{\displaystyle \frac{1}{N}}B_c^{\gamma \gamma ,(1),[1]}+N_FC_c^{\gamma \gamma ,(1),[1]}.`$
For clarity, the dependence on the helicities has been suppressed. Explicit expressions for $`A,\mathrm{},E`$ for $`q\overline{q}gg`$ in the physical region, $`s=s_{12}>0`$, $`t=s_{23}<0`$ and $`u=s_{13}<0`$, are given in Appendix B. Explicit results for $`q\overline{q}gg`$ scattering in the physical region, $`s=s_{12}>0`$, $`t=s_{23}<0`$ and $`u=s_{13}<0`$ are given in Appendix B. Formulae for the other channels and processes are available from the authors as FORM outputs. We have checked that the finite one-loop helicity amplitudes presented here are in agreement with the helicity amplitudes for $`q\overline{q}gg`$ (and crossed processes) given for the HV scheme in Ref. .

### 4.4 Two-loop helicity amplitudes

The main results of this paper are the two-loop amplitudes for processes (81)–(92). As in the tree and one-loop cases, the helicity amplitudes can be directly extracted from the appropriate linear combination, Eq. (4), of unrenormalised two-loop tensor coefficients obtained by direct evaluation of the projectors given in Section 2.2 acting on two-loop graphs. Renormalisation is achieved via Eq. (4) and the two-loop master integrals expanded around $`ϵ=0`$. The finite remainder of the two-loop amplitudes are defined through Eq. (79) and, for each process can be decomposed according to the number of colours and massless quark flavours,

$`|\widehat{}_c^{gg,(2),[i],finite}`$ $`=`$ $`N^2A_c^{gg,(2),[i]}+B_c^{gg,(2),[i]}+{\displaystyle \frac{1}{N^2}}C_c^{gg,(2),[i]}+NN_FD_c^{gg,(2),[i]}`$ (106)
$`+{\displaystyle \frac{N_F}{N}}E_c^{gg,(2),[i]}+N_F^2F_c^{gg,(2),[i]}+{\displaystyle \frac{N_F^2}{N^2}}G_c^{gg,(2),[i]}`$
$`2\beta _0S|\widehat{}_c^{gg,(1),[i],finite}(\beta _1S+\beta _0^2S^2)|\widehat{}_c^{gg,(0),[i]},`$
for $`i=1,2`$ and

$`|\widehat{}_c^{gg,(2),[3],finite}`$ $`=`$ $`NH_c^{gg,(2),[3]}+{\displaystyle \frac{1}{N}}I_c^{gg,(2),[3]}+N_FJ_c^{gg,(2),[3]}`$
$`+{\displaystyle \frac{N_F}{N^2}}K_c^{gg,(2),[3]}+{\displaystyle \frac{N_F^2}{N}}L_c^{gg,(2),[3]}`$
$`2\beta _0S|\widehat{}_c^{gg,(1),[3],finite}(\beta _1S+\beta _0^2S^2)|\widehat{}_c^{gg,(0),[3]}.`$
For the other processes we have,

$`|\widehat{}_c^{g\gamma ,(2),[1],finite}`$ $`=`$ $`N^2A_c^{g\gamma ,(2),[1]}+B_c^{g\gamma ,(2),[1]}+{\displaystyle \frac{1}{N^2}}C_c^{g\gamma ,(2),[1]}+NN_FD_c^{g\gamma ,(2),[1]}`$ (108)
$`+{\displaystyle \frac{N_F}{N}}E_c^{g\gamma ,(2),[1]}+N_F^2F_c^{g\gamma ,(2),[1]}+\left({\displaystyle \underset{q}{}}Q_q\right)\left({\displaystyle \frac{4}{N}}N\right)G_c^{g\gamma ,(2),[1]}`$
$`{\displaystyle \frac{3}{2}}\beta _0S|\widehat{}_c^{g\gamma ,(1),[1],finite}\left({\displaystyle \frac{1}{2}}\beta _1S+{\displaystyle \frac{3}{8}}\beta _0^2S^2\right)|\widehat{}_c^{g\gamma ,(0),[1]},`$
$`|\widehat{}_c^{\gamma \gamma ,(2),[1],finite}`$ $`=`$ $`N^2A_c^{\gamma \gamma ,(2),[1]}+B_c^{\gamma \gamma ,(2),[1]}+{\displaystyle \frac{1}{N^2}}C_c^{\gamma \gamma ,(2),[1]}+NN_FD_c^{\gamma \gamma ,(2),[1]}`$
$`+{\displaystyle \frac{N_F}{N}}E_c^{\gamma \gamma ,(2),[1]}+\left({\displaystyle \underset{q}{}}Q_q^2\right)C_FF_c^{\gamma \gamma ,(2),[1]}\beta _0S|\widehat{}_c^{\gamma \gamma ,(1),[1],finite},`$
where $`Q_q`$ is the electric charge of a massless quark of flavour $`q`$.

As in the previous section, the dependence on the helicities has been suppressed. Explicit expressions for $`A,\mathrm{},L`$ for $`q\overline{q}gg`$ process in the physical region, $`s=s_{12}>0`$, $`t=s_{23}<0`$ and $`u=s_{13}<0`$, are given in Appendix A. Formulae for the other channels and processes are available from the authors as FORM outputs. We have checked that the finite two-loop helicity amplitudes presented here are in agreement with the helicity amplitudes for $`q\overline{q}gg`$ (and crossed processes) given for the HV scheme in Ref. .

## 5 Comparison with previous results

In recovering the square, or interference, of amplitudes in a particular channel $`c`$, the contraction of the $`n`$-loop colour vector $`|^{𝒫,(n)}`$ with a conjugate $`m`$-loop colour vector $`^{𝒫,(m)}|`$ obeys the rule

$`^{𝒫,(m)}|^{𝒫,(n)}`$ (110)
$`=`$ $`{\displaystyle \underset{\mathrm{helicities}}{}}{\displaystyle \underset{\mathrm{colours}}{}}|_c^{𝒫,(n)}_{\lambda _1\lambda _2\lambda _3\lambda _4}^{}|_c^{𝒫,(n)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$
where $`|_c^{𝒫,(n)}_{\lambda _1\lambda _2\lambda _3\lambda _4}`$ contains all of the spinor and colour information for a particular helicity. Expanding the colour and spinor factors $`𝒞`$ and $`𝒮`$, and dropping the explicit dependence on the helicities, we see that,

$`^{𝒫,(m)}|^{𝒫,(n)}`$ (111)
$`=`$ $`{\displaystyle \underset{\mathrm{helicities}}{}}{\displaystyle \underset{\mathrm{colours}}{}}{\displaystyle \underset{i,j}{}}𝒞_c^{𝒫,[i]}𝒞_c^{𝒫,[j]}𝒮_c^{𝒫,[i]}𝒮_c^{𝒫,[j]}|\widehat{}_c^{𝒫,(n),[i]}^{}|\widehat{}_c^{𝒫,(n),[j]}`$
$`=`$ $`{\displaystyle \underset{\mathrm{helicities}}{}}{\displaystyle \underset{i,j}{}}𝒞𝒞_{ij}^𝒫𝒮𝒮_{ij}^𝒫|\widehat{}_c^{𝒫,(n),[i]}^{}|\widehat{}_c^{𝒫,(n),[j]}.`$
The symmetric colour sum matrix $`𝒞𝒞^𝒫=_{\mathrm{colours}}𝒞_c^{𝒫,[i]}𝒞_c^{𝒫,[j]}`$ is given by

$$𝒞𝒞^{gg}=\frac{N^21}{2}\left(\begin{array}{ccc}\frac{N^21}{2N}& \frac{1}{2N}& 1\\ \frac{1}{2N}& \frac{N^21}{2N}& 1\\ 1& 1& 2N\end{array}\right),$$
(112)
and,

$$𝒞𝒞^{g\gamma }=\frac{N^21}{2},𝒞𝒞^{\gamma \gamma }=N.$$
(113)
In the $`s`$-channel, the process and helicity dependent spinor matrices $`𝒮𝒮_c^𝒫=𝒮_c^{𝒫,[i]}𝒮_c^{𝒫,[j]}`$ are given by,

$$𝒮𝒮_{s++}^{gg}=\frac{u}{s^2t}\left(\begin{array}{ccc}u^2& ut& us\\ tu& t^2& ts\\ su& ts& s^2\end{array}\right),𝒮𝒮_{s++}^{g\gamma }=𝒮𝒮_{s++}^{\gamma \gamma }=\frac{u}{t},$$
(114)
and

$$𝒮𝒮_{s+++}^{gg}=\frac{ut}{s^2}\left(\begin{array}{ccc}1& 1& 1\\ 1& 1& 1\\ 1& 1& 1\end{array}\right),𝒮𝒮_{s+++}^{g\gamma }=𝒮𝒮_{s+++}^{\gamma \gamma }=\frac{ut}{s^2}.$$
(115)
Similar matrices for the other helicities are obtained by Bose symmetry, while corresponding matrices for the $`u`$ and $`t`$-channels are obtained by crossing symmetry,

$$𝒮𝒮_u^𝒫=𝒮𝒮_s^𝒫(st),𝒮𝒮_t^𝒫=𝒮𝒮_s^𝒫(su).$$
(116)
The explicit minus sign accounts for the exchange of a single fermion between the initial and final states.

Eqs. (110) and (111) can be used to recreate the full interference of tree and two-loop amplitudes. However, we notice that we can immediately separate the singularities from the finite parts using Eq. (79). The finite remainder, denoted by $`_{inite}^{2\times 0}(s,t,u)`$ in Refs.  is thus given by,

$$2Re\underset{\mathrm{helicities}}{}\underset{i,j}{}𝒞𝒞_{ij}^𝒫𝒮𝒮_{ij}^𝒫|\widehat{}_c^{𝒫,(0),[i]}^{}|\widehat{}_c^{𝒫,(2),[j],finite}.$$
(117)
Note that by using helicities, we are implicitly treating the external particle states in 4-dimensions - the HV scheme. However, the singular contributions to the renormalised amplitudes defined in Eqs. (78), (79) are given relative to the tree amplitude and therefore any additional terms that would be produced in CDR are automatically removed. Treating the external states differently for both the $`n`$-loop and tree amplitudes does not alter the finite contribution. Direct application of Eq. (110) for the real part of the tree and two-loop interference ($`n=0`$ and $`m=2`$) in the HV scheme therefore generates the same finite parts given in section 3.2 of Ref.  for quark-gluon scattering and in Section 4 and Appendix A of Ref.  obtained using CDR.

Similarly, the square of one-loop graphs is obtained using $`n=1`$ and $`m=1`$ in Eq. (110). In principle, the one-loop amplitude should be expanded through to $`𝒪(ϵ^2)`$. However, as noted in , the finite self interference only contains logarithms (up to the fourth power), but does not contain triple and quartic polylogs. These terms naturally arise when the $`𝒪(ϵ)`$ and $`𝒪(ϵ^2)`$ terms in the expansion of the one-loop box are multiplied by singularities from the conjugate amplitude. Such terms appear in the Catani pole structure, $`^{(1)}|𝑰^{(1)}|^{(0)}`$ where $`|^{(1)}`$ does need to be expanded through to $`𝒪(ϵ^2)`$. The finite parts of the renormalised one-loop self interference remaining after Catani’s prediction for the singularities has been removed should therefore be identical to those obtained by summing the finite parts using Eq. (110). In this instance, the finite remainder denoted by $`_{inite}^{1\times 1}(s,t,u)`$ in  is obtained from

$$\underset{\mathrm{helicities}}{}\underset{i,j}{}𝒞𝒞_{ij}^𝒫𝒮𝒮_{ij}^𝒫|\widehat{}_c^{𝒫,(1),[i],finite}^{}|\widehat{}_c^{𝒫,(1),[j],finite}.$$
(118)
We have checked that this procedure correctly reproduces the coefficients given in Section 5 and Appendix B of Ref.  for the processes involving one or more photons. Unfortunately, the results of Ref.  are not presented in quite the same way and a direct comparison is more difficult.<sup>2</sup><sup>2</sup>2In fact there is a misprint in the renormalisation term, $`_𝒜`$, of Eq. 4.6 of Ref. . However, we have checked the finite remainders after the Catani poles have been subtracted, both directly from the square of the one-loop amplitude (the CDR scheme) and the one-loop square reconstructed from the helicity amplitudes (the HV scheme) with the HV results of Ref.  and again find agreement.

## 6 Conclusions

In this paper we have presented the analytic expressions for the one- and two-loop virtual helicity amplitudes for quark-gauge boson scattering, $`q\overline{q}gg`$, $`q\overline{q}g\gamma `$ and $`q\overline{q}\gamma \gamma `$ and those processes related by crossing symmetry. These amplitudes were obtained by the construction of $`d`$-dimensional projection operators that extract the coefficients of the most general tensor representation for the amplitude, order by order in perturbation theory. Once the renormalised tensor coefficients are known, the renormalised helicity amplitudes can be constructed straightforwardly as a linear combination of the coefficients.

At the level of the diagrammatic calculation, we applied the projection operators directly to the Feynman diagrams that contribute at each order in $`\alpha _s`$. The projectors saturate the tensor structure of the Feynman diagram and yield a set of scalar integrals that can be related to a basis set of master integrals with the application of widely used reduction algorithms. In fact, we are using exactly the same tools as we did in , except in that case, the projector was merely the conjugate tree amplitude. Because the projectors exist in $`d`$-dimensions, conventional dimensional regularisation is preserved and there is no ambiguity in dealing with $`\gamma _5`$. Once the general tensor coefficients are determined, the external particles can be treated in 4-dimensions (the ’t Hooft-Veltman scheme) and standard helicity methods used to construct the helicity amplitudes.

By summing over helicities and colours, the full interference of tree and two-loop graphs can be reconstructed. In previous work , the finite part has been separated from infrared singular parts using the Catani formalism . In this procedure, the singular operators multiply tree and one-loop amplitudes. Changing the scheme for the external particles therefore changes the overall pole contribution. However, the finite remainder left after subtracting the poles is invariant under changing from the conventional dimensional regularisation scheme to the ’t Hooft-Veltman scheme. Therefore, the finite parts of the helicity amplitudes presented here precisely reproduce the finite parts of the interference of tree and two-loop graphs given in . Similarly, the finite parts of the self interference of one-loop amplitudes can be reconstructed purely from the finite parts of the one-loop helicity amplitudes. In addition, the finite remainders of the two-loop amplitudes for quark-gluon scattering in the HV scheme agree with the independent calculation of Ref. .

For the case of quark-gluon scattering, the presence of four coloured external particles gives rise to additional colour correlations proportional to $`1/ϵ`$. While the existence of such additional non-trivial colour correlations was pointed out in Ref.  and are expected on general grounds , the precise form of the colour structure of $`𝑯^{𝒫,(2)}`$ was not predicted. So far, the results for gluon-gluon scattering  as well as those for quark-gluon scattering presented in this paper and in Ref. , support the form given in Eqs. (60) and (61). These colour correlations vanish when the interference of tree and two-loop amplitudes is constructed.

The more general goal is to evaluate the NNLO QCD corrections to the associated scattering processes which should yield a significant improvement in theoretical stability, see for example . However, the reality is that the phenomenology is still some way off and there is still a significant amount of work remaining to be done. The two-loop corrections to the $`22`$ partonic process presented here, need to be combined with infrared singular contributions from both the one-loop $`23`$ partonic process, when one of the partons becomes unresolved  and the tree-level $`24`$ partonic process when two partons become unresolved . These different contributions need to be put together in such a way to make the explicit singularities analytically cancel, so that the four-dimensional limit can be safely taken, and then by encoding the finite remainders in a stable numerical implementation. There has been recent progress in identifying appropriate counterterms  and we are optimistic that the technical problems will soon be overcome enabling predictions of both dijet, jet and photon and diphoton production at the Tevatron and the LHC at an accuracy of $`𝒪(5\%)`$ which is expected to be competitive with the experimental measurements.

## Acknowledgements

We thank Babis Anastasiou for helpful discussions and Lance Dixon, Zvi Bern and Abilio De Freitas for communicating their results to us at an early stage and help in comparing our results. M.E. Tejeda-Yeomans thanks the IPPP for hospitality.

## Appendix A Finite two-loop contributions

In this appendix we give explicit formulae for the coefficients $`A,\mathrm{},L`$ for the $`q\overline{q}gg`$ process. Expressions for quark-gluon scattering as well as the processes involving photons are available from the authors as a FORM output file.

As usual, the polylogarithms $`\mathrm{Li}_n(w)`$ are defined by

$`\mathrm{Li}_n(w)`$ $`=`$ $`{\displaystyle _0^w}{\displaystyle \frac{dt}{t}}\mathrm{Li}_{n1}(t)\mathrm{for}n=2,3,4`$
$`\mathrm{Li}_2(w)`$ $`=`$ $`{\displaystyle _0^w}{\displaystyle \frac{dt}{t}}\mathrm{log}(1t).`$ (119)
Using the standard polylogarithm identities , we retain the polylogarithms with arguments $`x`$, $`1x`$ and $`(x1)/x`$, where

$$x=\frac{t}{s},y=\frac{u}{s}=1x,z=\frac{u}{t}=\frac{x1}{x}.$$
(120)
For convenience, we also introduce the following logarithms

$$X=\mathrm{log}\left(\frac{t}{s}\right),Y=\mathrm{log}\left(\frac{u}{s}\right),S=\mathrm{log}\left(\frac{s}{\mu ^2}\right),$$
(121)
where $`\mu `$ is the renormalisation scale.

### A.1 $`q(p_2,+)+\overline{q}(p_1,)g(p_3,+)+g(p_4,+)`$

$`A_{s+++}^{gg,(2),[1]}`$ $`=`$ $`\left({\displaystyle \frac{25}{24}}X^2{\displaystyle \frac{35}{18}}X{\displaystyle \frac{175}{54}}\right)+\left({\displaystyle \frac{17}{8}}X^2{\displaystyle \frac{5}{2}}X\right){\displaystyle \frac{t}{u}}{\displaystyle \frac{5}{4}}X^2{\displaystyle \frac{t^2}{u^2}}`$
$`+\left({\displaystyle \frac{1}{8}}X^2+{\displaystyle \frac{11}{8}}X{\displaystyle \frac{475}{144}}\right){\displaystyle \frac{s}{t}}`$
$`+\left(\left({\displaystyle \frac{25}{12}}X{\displaystyle \frac{19}{6}}\right)+\left({\displaystyle \frac{17}{4}}X{\displaystyle \frac{5}{2}}\right){\displaystyle \frac{t}{u}}{\displaystyle \frac{5}{2}}X{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{4}}X{\displaystyle \frac{11}{24}}\right){\displaystyle \frac{s}{t}}\right)i\pi `$
$`B_{s+++}^{gg,(2),[1]}`$ $`=`$ $`\left({\displaystyle \frac{1}{4}}X^2{\displaystyle \frac{7}{4}}XY+{\displaystyle \frac{1}{2}}X+{\displaystyle \frac{7}{8}}Y^2{\displaystyle \frac{3}{2}}Y+{\displaystyle \frac{7}{8}}\pi ^2{\displaystyle \frac{2}{3}}\right)+\left({\displaystyle \frac{9}{8}}X^2X\right){\displaystyle \frac{t}{u}}`$
$`{\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{3}{4}}X^2{\displaystyle \frac{3}{2}}XY+{\displaystyle \frac{3}{4}}Y^2+{\displaystyle \frac{3}{4}}\pi ^2\right){\displaystyle \frac{t}{s}}+\left({\displaystyle \frac{7}{4}}X+{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{3}{4}}Y{\displaystyle \frac{139}{36}}\right){\displaystyle \frac{s}{t}}`$
$`+{\displaystyle \frac{1}{8}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left(\left({\displaystyle \frac{5}{4}}X1\right)+\left({\displaystyle \frac{9}{4}}X1\right){\displaystyle \frac{t}{u}}X{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{2}}Y{\displaystyle \frac{5}{6}}\right){\displaystyle \frac{s}{t}}+{\displaystyle \frac{1}{4}}Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`C_{s+++}^{gg,(2),[1]}`$ $`=`$ $`\left({\displaystyle \frac{7}{8}}X^2{\displaystyle \frac{7}{4}}XY+2X+{\displaystyle \frac{7}{8}}Y^2{\displaystyle \frac{3}{2}}Y+{\displaystyle \frac{7}{8}}\pi ^2\right)+\left({\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{1}{2}}X\right){\displaystyle \frac{t}{u}}`$
$`+{\displaystyle \frac{1}{4}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{3}{4}}X^2{\displaystyle \frac{3}{2}}XY+{\displaystyle \frac{3}{4}}Y^2+{\displaystyle \frac{3}{4}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`+\left({\displaystyle \frac{1}{8}}X^2+{\displaystyle \frac{3}{8}}X+{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{3}{4}}Y{\displaystyle \frac{17}{16}}\right){\displaystyle \frac{s}{t}}+{\displaystyle \frac{1}{8}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left({\displaystyle \frac{1}{2}}+\left({\displaystyle \frac{1}{2}}X+{\displaystyle \frac{1}{2}}\right){\displaystyle \frac{t}{u}}+{\displaystyle \frac{1}{2}}X{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{4}}X+{\displaystyle \frac{1}{2}}Y{\displaystyle \frac{3}{8}}\right){\displaystyle \frac{s}{t}}+{\displaystyle \frac{1}{4}}Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`D_{s+++}^{gg,(2),[1]}`$ $`=`$ $`\left({\displaystyle \frac{5}{12}}X^2+{\displaystyle \frac{7}{18}}X+{\displaystyle \frac{397}{108}}\right)+\left({\displaystyle \frac{3}{4}}X^2+X\right){\displaystyle \frac{t}{u}}+{\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{4}}X+{\displaystyle \frac{37}{36}}\right){\displaystyle \frac{s}{t}}`$
$`+\left(\left({\displaystyle \frac{5}{6}}X+{\displaystyle \frac{11}{6}}\right)+\left({\displaystyle \frac{3}{2}}X+1\right){\displaystyle \frac{t}{u}}+X{\displaystyle \frac{t^2}{u^2}}+{\displaystyle \frac{1}{12}}{\displaystyle \frac{s}{t}}\right)i\pi `$
$`E_{s+++}^{gg,(2),[1]}`$ $`=`$ $`\left({\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{3}{2}}X+{\displaystyle \frac{1}{4}}Y^23Y+{\displaystyle \frac{1}{4}}\pi ^2+{\displaystyle \frac{11}{12}}\right)+\left(X^2{\displaystyle \frac{3}{2}}X\right){\displaystyle \frac{t}{u}}`$
$`{\displaystyle \frac{3}{4}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{2}}X^2XY+{\displaystyle \frac{1}{2}}Y^2+{\displaystyle \frac{1}{2}}\pi ^2\right){\displaystyle \frac{t}{s}}+\left({\displaystyle \frac{1}{4}}X+{\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{3}{2}}Y+{\displaystyle \frac{19}{36}}\right){\displaystyle \frac{s}{t}}`$
$`+{\displaystyle \frac{3}{4}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left(\left({\displaystyle \frac{1}{2}}X{\displaystyle \frac{3}{2}}\right)+\left(2X{\displaystyle \frac{3}{2}}\right){\displaystyle \frac{t}{u}}{\displaystyle \frac{3}{2}}X{\displaystyle \frac{t^2}{u^2}}+\left(Y{\displaystyle \frac{17}{12}}\right){\displaystyle \frac{s}{t}}+{\displaystyle \frac{3}{2}}Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`F_{s+++}^{gg,(2),[1]}`$ $`=`$ $`{\displaystyle \frac{1}{6}}i\pi +\left({\displaystyle \frac{1}{18}}X{\displaystyle \frac{5}{27}}\right)`$
$`G_{s+++}^{gg,(2),[1]}`$ $`=`$ $`0`$ (128)
$`A_{s+++}^{gg,(2),[2]}`$ $`=`$ $`\left({\displaystyle \frac{1}{24}}Y^2+{\displaystyle \frac{59}{72}}Y{\displaystyle \frac{25}{432}}\right)+\left({\displaystyle \frac{1}{8}}Y^2+{\displaystyle \frac{11}{8}}Y{\displaystyle \frac{475}{144}}\right){\displaystyle \frac{t}{u}}+\left({\displaystyle \frac{3}{8}}Y^2{\displaystyle \frac{5}{2}}Y\right){\displaystyle \frac{s}{t}}`$
$`+{\displaystyle \frac{5}{4}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left(\left({\displaystyle \frac{1}{12}}Y+{\displaystyle \frac{5}{24}}\right)+\left({\displaystyle \frac{1}{4}}Y{\displaystyle \frac{11}{24}}\right){\displaystyle \frac{t}{u}}+\left({\displaystyle \frac{3}{4}}Y{\displaystyle \frac{5}{2}}\right){\displaystyle \frac{s}{t}}+{\displaystyle \frac{5}{2}}Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`B_{s+++}^{gg,(2),[2]}`$ $`=`$ $`\left({\displaystyle \frac{1}{4}}XY+{\displaystyle \frac{3}{4}}X{\displaystyle \frac{1}{8}}Y^2+{\displaystyle \frac{1}{4}}Y{\displaystyle \frac{1}{8}}\pi ^2{\displaystyle \frac{115}{36}}\right)+\left({\displaystyle \frac{3}{4}}X+{\displaystyle \frac{7}{4}}Y{\displaystyle \frac{139}{36}}\right){\displaystyle \frac{t}{u}}`$
$`{\displaystyle \frac{1}{8}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{3}{4}}X^2{\displaystyle \frac{3}{2}}XY+{\displaystyle \frac{3}{4}}Y^2+{\displaystyle \frac{3}{4}}\pi ^2\right){\displaystyle \frac{t}{s}}+\left({\displaystyle \frac{1}{8}}Y^2Y\right){\displaystyle \frac{s}{t}}+{\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left(\left({\displaystyle \frac{1}{4}}X{\displaystyle \frac{5}{6}}\right){\displaystyle \frac{5}{6}}{\displaystyle \frac{t}{u}}{\displaystyle \frac{1}{4}}X{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{4}}Y1\right){\displaystyle \frac{s}{t}}+Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`C_{s+++}^{gg,(2),[2]}`$ $`=`$ $`\left({\displaystyle \frac{1}{4}}XY+{\displaystyle \frac{3}{4}}X{\displaystyle \frac{9}{8}}Y{\displaystyle \frac{1}{8}}\pi ^2{\displaystyle \frac{17}{16}}\right)+\left({\displaystyle \frac{3}{4}}X+{\displaystyle \frac{1}{8}}Y^2+{\displaystyle \frac{3}{8}}Y{\displaystyle \frac{17}{16}}\right){\displaystyle \frac{t}{u}}`$
$`{\displaystyle \frac{1}{8}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{3}{4}}X^2{\displaystyle \frac{3}{2}}XY+{\displaystyle \frac{3}{4}}Y^2+{\displaystyle \frac{3}{4}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`+\left({\displaystyle \frac{1}{4}}Y^2+{\displaystyle \frac{1}{2}}Y\right){\displaystyle \frac{s}{t}}{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left(\left({\displaystyle \frac{1}{4}}X+{\displaystyle \frac{1}{4}}Y{\displaystyle \frac{3}{8}}\right)+\left({\displaystyle \frac{1}{4}}Y{\displaystyle \frac{3}{8}}\right){\displaystyle \frac{t}{u}}{\displaystyle \frac{1}{4}}X{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{2}}Y+{\displaystyle \frac{1}{2}}\right){\displaystyle \frac{s}{t}}{\displaystyle \frac{1}{2}}Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`D_{s+++}^{gg,(2),[2]}`$ $`=`$ $`\left({\displaystyle \frac{1}{6}}Y^2+{\displaystyle \frac{13}{36}}Y{\displaystyle \frac{143}{54}}\right)+\left({\displaystyle \frac{1}{4}}Y+{\displaystyle \frac{37}{36}}\right){\displaystyle \frac{t}{u}}+\left({\displaystyle \frac{1}{4}}Y^2+Y\right){\displaystyle \frac{s}{t}}{\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left(\left({\displaystyle \frac{1}{3}}Y{\displaystyle \frac{3}{4}}\right)+{\displaystyle \frac{1}{12}}{\displaystyle \frac{t}{u}}+\left({\displaystyle \frac{1}{2}}Y+1\right){\displaystyle \frac{s}{t}}Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`E_{s+++}^{gg,(2),[2]}`$ $`=`$ $`\left({\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{3}{2}}X+{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{13}{4}}Y+{\displaystyle \frac{1}{4}}\pi ^2{\displaystyle \frac{7}{18}}\right)+\left(X^2{\displaystyle \frac{3}{2}}X{\displaystyle \frac{1}{4}}Y+{\displaystyle \frac{19}{36}}\right){\displaystyle \frac{t}{u}}`$
$`{\displaystyle \frac{3}{4}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{2}}X^2XY+{\displaystyle \frac{1}{2}}Y^2+{\displaystyle \frac{1}{2}}\pi ^2\right){\displaystyle \frac{t}{s}}+\left({\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{3}{2}}Y\right){\displaystyle \frac{s}{t}}+{\displaystyle \frac{3}{4}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left(\left({\displaystyle \frac{1}{2}}X{\displaystyle \frac{17}{12}}\right)+\left(2X{\displaystyle \frac{17}{12}}\right){\displaystyle \frac{t}{u}}{\displaystyle \frac{3}{2}}X{\displaystyle \frac{t^2}{u^2}}+\left(Y{\displaystyle \frac{3}{2}}\right){\displaystyle \frac{s}{t}}+{\displaystyle \frac{3}{2}}Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`F_{s+++}^{gg,(2),[2]}`$ $`=`$ $`{\displaystyle \frac{1}{6}}i\pi +\left({\displaystyle \frac{1}{18}}Y+{\displaystyle \frac{5}{27}}\right)`$
$`G_{s+++}^{gg,(2),[2]}`$ $`=`$ $`0`$ (135)
$`H_{s+++}^{gg,(2),[3]}`$ $`=`$ $`\left({\displaystyle \frac{1}{48}}X^2{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{11}{18}}X+{\displaystyle \frac{5}{24}}Y^2{\displaystyle \frac{16}{9}}Y+{\displaystyle \frac{1}{4}}\pi ^2\right)`$
$`+\left({\displaystyle \frac{11}{16}}X^2{\displaystyle \frac{1}{4}}XY{\displaystyle \frac{7}{6}}X+{\displaystyle \frac{1}{8}}\pi ^2\right){\displaystyle \frac{t}{u}}{\displaystyle \frac{5}{8}}X^2{\displaystyle \frac{t^2}{u^2}}+{\displaystyle \frac{5}{8}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left({\displaystyle \frac{1}{4}}X^2{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{1}{4}}Y^2+{\displaystyle \frac{1}{4}}\pi ^2\right){\displaystyle \frac{t}{s}}+\left({\displaystyle \frac{1}{4}}XY+{\displaystyle \frac{9}{16}}Y^2{\displaystyle \frac{7}{6}}Y+{\displaystyle \frac{1}{8}}\pi ^2\right){\displaystyle \frac{s}{t}}`$
$`+(({\displaystyle \frac{13}{24}}X{\displaystyle \frac{1}{12}}Y{\displaystyle \frac{7}{6}})+({\displaystyle \frac{13}{8}}X{\displaystyle \frac{1}{4}}Y{\displaystyle \frac{7}{6}}){\displaystyle \frac{t}{u}}{\displaystyle \frac{5}{4}}X{\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{1}{4}}X+{\displaystyle \frac{7}{8}}Y{\displaystyle \frac{7}{6}}){\displaystyle \frac{s}{t}}+{\displaystyle \frac{5}{4}}Y{\displaystyle \frac{s^2}{t^2}})i\pi `$
$`I_{s+++}^{gg,(2),[3]}`$ $`=`$ $`\left({\displaystyle \frac{1}{16}}X^2+{\displaystyle \frac{1}{4}}XY{\displaystyle \frac{11}{12}}X{\displaystyle \frac{1}{4}}Y^2+Y{\displaystyle \frac{1}{8}}\pi ^2\right)`$
$`+\left({\displaystyle \frac{3}{16}}X^2{\displaystyle \frac{1}{4}}XY+{\displaystyle \frac{1}{12}}X+{\displaystyle \frac{1}{8}}\pi ^2\right){\displaystyle \frac{t}{u}}`$
$`+\left({\displaystyle \frac{1}{2}}X^2+XY{\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{1}{2}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`+\left({\displaystyle \frac{1}{4}}XY+{\displaystyle \frac{3}{16}}Y^2+{\displaystyle \frac{1}{12}}Y+{\displaystyle \frac{1}{8}}\pi ^2\right){\displaystyle \frac{s}{t}}`$
$`+\left(\left({\displaystyle \frac{1}{8}}X{\displaystyle \frac{1}{4}}Y+{\displaystyle \frac{1}{12}}\right)+\left({\displaystyle \frac{1}{8}}X{\displaystyle \frac{1}{4}}Y+{\displaystyle \frac{1}{12}}\right){\displaystyle \frac{t}{u}}+\left({\displaystyle \frac{1}{4}}X+{\displaystyle \frac{1}{8}}Y+{\displaystyle \frac{1}{12}}\right){\displaystyle \frac{s}{t}}\right)i\pi `$
$`J_{s+++}^{gg,(2),[3]}`$ $`=`$ $`\left({\displaystyle \frac{1}{24}}X^2+{\displaystyle \frac{1}{2}}XY{\displaystyle \frac{55}{72}}X{\displaystyle \frac{1}{3}}Y^2+{\displaystyle \frac{47}{36}}Y{\displaystyle \frac{1}{4}}\pi ^2\right)+\left({\displaystyle \frac{3}{8}}X^2+{\displaystyle \frac{13}{24}}X\right){\displaystyle \frac{t}{u}}`$
$`+{\displaystyle \frac{1}{4}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{2}}X^2+XY{\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{1}{2}}\pi ^2\right){\displaystyle \frac{t}{s}}+\left({\displaystyle \frac{1}{8}}Y^2+{\displaystyle \frac{13}{24}}Y\right){\displaystyle \frac{s}{t}}`$
$`{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+(({\displaystyle \frac{5}{12}}X{\displaystyle \frac{1}{6}}Y+{\displaystyle \frac{13}{24}})+({\displaystyle \frac{3}{4}}X+{\displaystyle \frac{13}{24}}){\displaystyle \frac{t}{u}}+{\displaystyle \frac{1}{2}}X{\displaystyle \frac{t^2}{u^2}}+({\displaystyle \frac{1}{4}}Y+{\displaystyle \frac{13}{24}}){\displaystyle \frac{s}{t}}`$
$`{\displaystyle \frac{1}{2}}Y{\displaystyle \frac{s^2}{t^2}})i\pi `$
$`K_{s+++}^{gg,(2),[3]}`$ $`=`$ $`\left({\displaystyle \frac{1}{4}}XY{\displaystyle \frac{17}{24}}X{\displaystyle \frac{1}{8}}Y^2+{\displaystyle \frac{3}{2}}Y{\displaystyle \frac{1}{8}}\pi ^2\right)+\left({\displaystyle \frac{1}{2}}X^2+{\displaystyle \frac{19}{24}}X\right){\displaystyle \frac{t}{u}}`$
$`+{\displaystyle \frac{3}{8}}X^2{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{1}{2}}XY{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{1}{4}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`+\left({\displaystyle \frac{1}{4}}Y^2+{\displaystyle \frac{19}{24}}Y\right){\displaystyle \frac{s}{t}}{\displaystyle \frac{3}{8}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+\left(\left({\displaystyle \frac{1}{4}}X+{\displaystyle \frac{19}{24}}\right)+\left(X+{\displaystyle \frac{19}{24}}\right){\displaystyle \frac{t}{u}}+{\displaystyle \frac{3}{4}}X{\displaystyle \frac{t^2}{u^2}}+\left({\displaystyle \frac{1}{2}}Y+{\displaystyle \frac{19}{24}}\right){\displaystyle \frac{s}{t}}{\displaystyle \frac{3}{4}}Y{\displaystyle \frac{s^2}{t^2}}\right)i\pi `$
$`L_{s+++}^{gg,(2),[3]}`$ $`=`$ $`\left({\displaystyle \frac{1}{36}}X{\displaystyle \frac{1}{36}}Y\right)`$ (140)
### A.2 $`q(p_2,+)+\overline{q}(p_1,)g(p_3,+)+g(p_4,)`$

$`A_{s++}^{gg,(2),[1]}`$ $`=`$ $`({\displaystyle \frac{1}{2}}\mathrm{Li}_4(x)+{\displaystyle \frac{35}{12}}\mathrm{Li}_3(x){\displaystyle \frac{1}{4}}\mathrm{Li}_2(x)X^2{\displaystyle \frac{35}{12}}\mathrm{Li}_2(x)X+{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)\pi ^2{\displaystyle \frac{1}{16}}X^4`$
$`{\displaystyle \frac{1}{6}}X^3Y+{\displaystyle \frac{77}{72}}X^3{\displaystyle \frac{35}{24}}X^2Y+{\displaystyle \frac{1}{3}}X^2\pi ^2{\displaystyle \frac{505}{144}}X^2+{\displaystyle \frac{1}{3}}XY\pi ^2+2X\zeta _3`$
$`+{\displaystyle \frac{143}{72}}X\pi ^2{\displaystyle \frac{2273}{432}}X{\displaystyle \frac{241}{24}}\zeta _3{\displaystyle \frac{19}{288}}\pi ^4+{\displaystyle \frac{5}{4}}\pi ^2+{\displaystyle \frac{36077}{1728}})`$
$`+({\displaystyle \frac{1}{4}}\mathrm{Li}_3(x)+{\displaystyle \frac{1}{4}}\mathrm{Li}_2(x)X{\displaystyle \frac{7}{24}}X^3+{\displaystyle \frac{1}{8}}X^2Y+{\displaystyle \frac{31}{16}}X^2+{\displaystyle \frac{2}{3}}X\pi ^2+{\displaystyle \frac{203}{144}}X`$
$`{\displaystyle \frac{1}{2}}\zeta _3+{\displaystyle \frac{17}{36}}\pi ^2+{\displaystyle \frac{475}{144}}){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{1}{2}}\mathrm{Li}_3(x){\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{1}{3}}X^3{\displaystyle \frac{1}{4}}X^2Y+{\displaystyle \frac{11}{16}}X^2{\displaystyle \frac{5}{6}}X\pi ^2+{\displaystyle \frac{475}{72}}X`$
$`+\zeta _3{\displaystyle \frac{23}{18}}\pi ^2){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{3}{2}}\mathrm{Li}_4(x)\mathrm{Li}_3(x)X+{\displaystyle \frac{1}{4}}\mathrm{Li}_2(x)X^2+{\displaystyle \frac{7}{48}}X^4{\displaystyle \frac{43}{36}}X^3{\displaystyle \frac{2}{3}}X^2\pi ^2`$
$`+{\displaystyle \frac{475}{144}}X^2+X\zeta _3{\displaystyle \frac{23}{18}}X\pi ^2{\displaystyle \frac{1}{60}}\pi ^4){\displaystyle \frac{t^3}{u^3}}`$
$`+(({\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{35}{12}}\mathrm{Li}_2(x){\displaystyle \frac{1}{4}}X^3{\displaystyle \frac{1}{2}}X^2Y+{\displaystyle \frac{11}{8}}X^2{\displaystyle \frac{35}{12}}XY`$
$`+{\displaystyle \frac{1}{6}}X\pi ^2{\displaystyle \frac{7}{24}}X+2\zeta _3+{\displaystyle \frac{11}{24}}\pi ^2+{\displaystyle \frac{823}{432}})`$
$`+\left({\displaystyle \frac{1}{4}}\mathrm{Li}_2(x){\displaystyle \frac{7}{8}}X^2+{\displaystyle \frac{1}{4}}XY+{\displaystyle \frac{49}{24}}X+{\displaystyle \frac{1}{12}}\pi ^2+{\displaystyle \frac{467}{144}}\right){\displaystyle \frac{t}{u}}`$
$`+\left({\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)+X^2{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{121}{24}}X{\displaystyle \frac{1}{6}}\pi ^2+{\displaystyle \frac{475}{72}}\right){\displaystyle \frac{t^2}{u^2}}`$
$`+(\mathrm{Li}_3(x)+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{7}{12}}X^3{\displaystyle \frac{7}{4}}X^2{\displaystyle \frac{1}{6}}X\pi ^2+{\displaystyle \frac{475}{72}}X+\zeta _3){\displaystyle \frac{t^3}{u^3}})i\pi `$
$`B_{s++}^{gg,(2),[1]}`$ $`=`$ $`(2\mathrm{Li}_4(x)7\mathrm{Li}_4(y)+\mathrm{Li}_4(z)+\mathrm{Li}_3(x)X\mathrm{Li}_3(x)Y+{\displaystyle \frac{9}{2}}\mathrm{Li}_3(x)+2\mathrm{Li}_3(y)X`$
$`+4\mathrm{Li}_3(y)Y+6\mathrm{Li}_3(y){\displaystyle \frac{5}{2}}\mathrm{Li}_2(x)X^2{\displaystyle \frac{9}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{1}{6}}\mathrm{Li}_2(x)\pi ^24\mathrm{Li}_2(y)XY`$
$`6\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{6}}X^4{\displaystyle \frac{11}{6}}X^3Y{\displaystyle \frac{25}{36}}X^3{\displaystyle \frac{3}{2}}X^2Y^2+{\displaystyle \frac{13}{12}}X^2Y+{\displaystyle \frac{3}{4}}X^2\pi ^2`$
$`{\displaystyle \frac{95}{36}}X^2+{\displaystyle \frac{5}{6}}XY^3{\displaystyle \frac{13}{2}}XY^2{\displaystyle \frac{1}{2}}XY\pi ^2+{\displaystyle \frac{25}{9}}XY2X\zeta _3{\displaystyle \frac{40}{9}}X\pi ^2+{\displaystyle \frac{97}{12}}X`$
$`{\displaystyle \frac{1}{12}}Y^4+{\displaystyle \frac{5}{3}}Y^3{\displaystyle \frac{1}{4}}Y^2\pi ^2{\displaystyle \frac{191}{72}}Y^2+Y\zeta _3+{\displaystyle \frac{7}{6}}Y\pi ^2+{\displaystyle \frac{3}{2}}Y+{\displaystyle \frac{11}{72}}\zeta _3`$
$`+{\displaystyle \frac{73}{1440}}\pi ^4{\displaystyle \frac{329}{144}}\pi ^2{\displaystyle \frac{14135}{1296}})`$
$`+({\displaystyle \frac{9}{2}}\mathrm{Li}_4(x)15\mathrm{Li}_4(y)11\mathrm{Li}_4(z)+7\mathrm{Li}_3(x)X7\mathrm{Li}_3(x)Y+{\displaystyle \frac{13}{2}}\mathrm{Li}_3(x)`$
$`+14\mathrm{Li}_3(y)X+6\mathrm{Li}_3(y){\displaystyle \frac{15}{4}}\mathrm{Li}_2(x)X^2{\displaystyle \frac{13}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{17}{6}}\mathrm{Li}_2(x)\pi ^2`$
$`10\mathrm{Li}_2(y)XY+2\mathrm{Li}_2(y)Y^26\mathrm{Li}_2(y)Y{\displaystyle \frac{11}{48}}X^4{\displaystyle \frac{31}{12}}X^35X^2Y^2+{\displaystyle \frac{1}{4}}X^2Y`$
$`{\displaystyle \frac{5}{6}}X^2\pi ^2+{\displaystyle \frac{305}{72}}X^2+{\displaystyle \frac{3}{2}}XY^3{\displaystyle \frac{29}{6}}XY^2{\displaystyle \frac{5}{6}}XY\pi ^2{\displaystyle \frac{37}{36}}XY13X\zeta _3{\displaystyle \frac{25}{3}}X\pi ^2`$
$`+{\displaystyle \frac{383}{36}}X{\displaystyle \frac{5}{4}}Y^2\pi ^2+7Y\zeta _3+{\displaystyle \frac{3}{2}}Y\pi ^2+{\displaystyle \frac{3}{4}}Y10\zeta _3+{\displaystyle \frac{19}{45}}\pi ^4{\displaystyle \frac{407}{72}}\pi ^2+{\displaystyle \frac{139}{36}}){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{9}{2}}\mathrm{Li}_4(x)12\mathrm{Li}_4(y)12\mathrm{Li}_4(z)+12\mathrm{Li}_3(x)X6\mathrm{Li}_3(x)Y+{\displaystyle \frac{7}{2}}\mathrm{Li}_3(x)`$
$`+12\mathrm{Li}_3(y)X{\displaystyle \frac{15}{4}}\mathrm{Li}_2(x)X^2{\displaystyle \frac{7}{2}}\mathrm{Li}_2(x)X3\mathrm{Li}_2(x)\pi ^26\mathrm{Li}_2(y)XY{\displaystyle \frac{5}{16}}X^4`$
$`+X^3Y{\displaystyle \frac{17}{12}}X^33X^2Y^2{\displaystyle \frac{5}{2}}X^2YX^2\pi ^2+{\displaystyle \frac{173}{24}}X^22XY\pi ^2+{\displaystyle \frac{3}{2}}XY`$
$`12X\zeta _3{\displaystyle \frac{77}{12}}X\pi ^2+{\displaystyle \frac{305}{36}}X+6Y\zeta _3+{\displaystyle \frac{1}{3}}Y\pi ^2{\displaystyle \frac{11}{2}}\zeta _3+{\displaystyle \frac{11}{20}}\pi ^4{\displaystyle \frac{103}{36}}\pi ^2){\displaystyle \frac{t^2}{u^2}}`$
$`+(2\mathrm{Li}_4(x)+\mathrm{Li}_3(x)X+{\displaystyle \frac{3}{2}}\mathrm{Li}_3(x){\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)\pi ^2+{\displaystyle \frac{1}{8}}X^4`$
$`{\displaystyle \frac{43}{36}}X^3{\displaystyle \frac{1}{4}}X^2\pi ^2+{\displaystyle \frac{359}{72}}X^2X\zeta _3{\displaystyle \frac{91}{36}}X\pi ^2{\displaystyle \frac{3}{2}}\zeta _3+{\displaystyle \frac{7}{90}}\pi ^4){\displaystyle \frac{t^3}{u^3}}`$
$`+\left({\displaystyle \frac{3}{8}}X^2+{\displaystyle \frac{3}{4}}XY{\displaystyle \frac{3}{8}}Y^2{\displaystyle \frac{3}{8}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`+((6\mathrm{Li}_3(y)\mathrm{Li}_2(x)X+{\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)4\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{3}}X^3{\displaystyle \frac{7}{2}}X^2Y+{\displaystyle \frac{37}{12}}X^2`$
$`+2XY^2{\displaystyle \frac{29}{6}}XY+{\displaystyle \frac{2}{3}}X\pi ^2+3X{\displaystyle \frac{1}{6}}Y^3+{\displaystyle \frac{3}{2}}Y^2{\displaystyle \frac{5}{6}}Y\pi ^2{\displaystyle \frac{91}{36}}Y\zeta _3+\pi ^2{\displaystyle \frac{13}{4}})`$
$`+(14\mathrm{Li}_3(y)+{\displaystyle \frac{5}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)6\mathrm{Li}_2(y)Y+{\displaystyle \frac{5}{12}}X^32X^2Y+{\displaystyle \frac{5}{4}}X^2`$
$`+3XY^2{\displaystyle \frac{19}{6}}XY{\displaystyle \frac{2}{3}}X\pi ^2+{\displaystyle \frac{299}{18}}X+{\displaystyle \frac{1}{6}}Y^3{\displaystyle \frac{11}{6}}Y^2{\displaystyle \frac{5}{6}}Y\pi ^2{\displaystyle \frac{37}{36}}Y6\zeta _3`$
$`+{\displaystyle \frac{4}{3}}\pi ^2+{\displaystyle \frac{119}{9}}){\displaystyle \frac{t}{u}}`$
$`+(6\mathrm{Li}_3(x)+12\mathrm{Li}_3(y){\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{7}{2}}\mathrm{Li}_2(x)6\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{4}}X^3+{\displaystyle \frac{1}{2}}X^2`$
$`5XY{\displaystyle \frac{3}{2}}X\pi ^2+{\displaystyle \frac{235}{12}}X+{\displaystyle \frac{3}{2}}Y6\zeta _3+{\displaystyle \frac{7}{12}}\pi ^2+{\displaystyle \frac{305}{36}}){\displaystyle \frac{t^2}{u^2}}`$
$`+(\mathrm{Li}_3(x){\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)+{\displaystyle \frac{1}{3}}X^3X^2{\displaystyle \frac{1}{6}}X\pi ^2+{\displaystyle \frac{359}{36}}X\zeta _3+{\displaystyle \frac{1}{4}}\pi ^2){\displaystyle \frac{t^3}{u^3}})i\pi `$
$`C_{s++}^{gg,(2),[1]}`$ $`=`$ $`({\displaystyle \frac{5}{2}}\mathrm{Li}_4(x)+\mathrm{Li}_4(y)+\mathrm{Li}_4(z)+\mathrm{Li}_3(x)X+\mathrm{Li}_3(x)Y{\displaystyle \frac{3}{4}}\mathrm{Li}_3(x)+3\mathrm{Li}_3(y)`$
$`{\displaystyle \frac{1}{4}}\mathrm{Li}_2(x)X^2+{\displaystyle \frac{3}{4}}\mathrm{Li}_2(x)X{\displaystyle \frac{1}{6}}\mathrm{Li}_2(x)\pi ^23\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{48}}X^4{\displaystyle \frac{1}{6}}X^3Y{\displaystyle \frac{3}{8}}X^3`$
$`+{\displaystyle \frac{1}{4}}X^2Y^2+{\displaystyle \frac{3}{8}}X^2Y+{\displaystyle \frac{1}{3}}X^2\pi ^2+{\displaystyle \frac{5}{16}}X^2{\displaystyle \frac{3}{2}}XY^2{\displaystyle \frac{1}{6}}XY\pi ^2X\zeta _3+X\pi ^2+{\displaystyle \frac{39}{16}}X`$
$`{\displaystyle \frac{11}{8}}Y^2Y\zeta _3+{\displaystyle \frac{3}{2}}Y+{\displaystyle \frac{9}{4}}\zeta _3+{\displaystyle \frac{31}{360}}\pi ^4+{\displaystyle \frac{5}{16}}\pi ^2{\displaystyle \frac{187}{64}})`$
$`+({\displaystyle \frac{7}{2}}\mathrm{Li}_4(x)+7\mathrm{Li}_4(y)+7\mathrm{Li}_4(z)\mathrm{Li}_3(x)X+7\mathrm{Li}_3(x)Y+{\displaystyle \frac{3}{4}}\mathrm{Li}_3(x)+9\mathrm{Li}_3(y)`$
$`{\displaystyle \frac{3}{4}}\mathrm{Li}_2(x)X^2{\displaystyle \frac{3}{4}}\mathrm{Li}_2(x)X+{\displaystyle \frac{1}{6}}\mathrm{Li}_2(x)\pi ^29\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{16}}X^4{\displaystyle \frac{7}{6}}X^3Y{\displaystyle \frac{31}{24}}X^3`$
$`+{\displaystyle \frac{7}{4}}X^2Y^2+{\displaystyle \frac{5}{8}}X^2Y+{\displaystyle \frac{2}{3}}X^2\pi ^2+{\displaystyle \frac{27}{16}}X^2{\displaystyle \frac{9}{2}}XY^2{\displaystyle \frac{7}{6}}XY\pi ^2{\displaystyle \frac{5}{4}}XY+X\zeta _3`$
$`+2X\pi ^2+{\displaystyle \frac{97}{16}}X{\displaystyle \frac{3}{2}}Y^27Y\zeta _3{\displaystyle \frac{2}{3}}Y\pi ^2+{\displaystyle \frac{3}{4}}Y{\displaystyle \frac{1}{2}}\zeta _3+{\displaystyle \frac{1}{90}}\pi ^4+{\displaystyle \frac{65}{24}}\pi ^2+{\displaystyle \frac{17}{16}}){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{3}{2}}\mathrm{Li}_4(x)+12\mathrm{Li}_4(y)+12\mathrm{Li}_4(z)6\mathrm{Li}_3(x)X+12\mathrm{Li}_3(x)Y+6\mathrm{Li}_3(y)`$
$`{\displaystyle \frac{3}{4}}\mathrm{Li}_2(x)X^2+\mathrm{Li}_2(x)\pi ^26\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{16}}X^42X^3Y{\displaystyle \frac{11}{12}}X^3+3X^2Y^2+{\displaystyle \frac{1}{4}}X^2Y`$
$`+{\displaystyle \frac{1}{4}}X^2\pi ^2+{\displaystyle \frac{35}{16}}X^23XY^22XY\pi ^2{\displaystyle \frac{3}{2}}XY+6X\zeta _3+{\displaystyle \frac{9}{4}}X\pi ^2+{\displaystyle \frac{23}{8}}X`$
$`12Y\zeta _3{\displaystyle \frac{2}{3}}Y\pi ^2+{\displaystyle \frac{5}{2}}\zeta _3{\displaystyle \frac{11}{60}}\pi ^4+{\displaystyle \frac{23}{12}}\pi ^2){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{5}{2}}\mathrm{Li}_4(x)+6\mathrm{Li}_4(y)+6\mathrm{Li}_4(z)4\mathrm{Li}_3(x)X+6\mathrm{Li}_3(x)Y{\displaystyle \frac{3}{2}}\mathrm{Li}_3(x)`$
$`{\displaystyle \frac{1}{4}}\mathrm{Li}_2(x)X^2+{\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)\pi ^2{\displaystyle \frac{1}{48}}X^4X^3Y+{\displaystyle \frac{3}{2}}X^2Y^2{\displaystyle \frac{1}{12}}X^2\pi ^2`$
$`+{\displaystyle \frac{11}{16}}X^2XY\pi ^2+4X\zeta _3+{\displaystyle \frac{5}{4}}X\pi ^26Y\zeta _3+{\displaystyle \frac{3}{2}}\zeta _3{\displaystyle \frac{5}{36}}\pi ^4){\displaystyle \frac{t^3}{u^3}}`$
$`+\left({\displaystyle \frac{3}{8}}X^2+{\displaystyle \frac{3}{4}}XY{\displaystyle \frac{3}{8}}Y^2{\displaystyle \frac{3}{8}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`+((2\mathrm{Li}_3(x){\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{15}{4}}\mathrm{Li}_2(x){\displaystyle \frac{1}{4}}X^3{\displaystyle \frac{9}{8}}X^2+{\displaystyle \frac{3}{4}}XY+{\displaystyle \frac{5}{8}}X`$
$`{\displaystyle \frac{11}{4}}Y2\zeta _3{\displaystyle \frac{1}{4}}\pi ^2+{\displaystyle \frac{63}{16}})`$
$`+(6\mathrm{Li}_3(x){\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{33}{4}}\mathrm{Li}_2(x){\displaystyle \frac{3}{4}}X^3{\displaystyle \frac{23}{8}}X^2+{\displaystyle \frac{5}{4}}XY+{\displaystyle \frac{17}{8}}X`$
$`{\displaystyle \frac{17}{4}}Y6\zeta _3{\displaystyle \frac{3}{4}}\pi ^2+{\displaystyle \frac{109}{16}}){\displaystyle \frac{t}{u}}`$
$`+(6\mathrm{Li}_3(x){\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)X+6\mathrm{Li}_2(x){\displaystyle \frac{3}{4}}X^3{\displaystyle \frac{5}{2}}X^2+{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{23}{8}}X`$
$`{\displaystyle \frac{3}{2}}Y6\zeta _3{\displaystyle \frac{3}{4}}\pi ^2+{\displaystyle \frac{23}{8}}){\displaystyle \frac{t^2}{u^2}}`$
$`+\left(2\mathrm{Li}_3(x){\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{3}{2}}\mathrm{Li}_2(x){\displaystyle \frac{1}{4}}X^3{\displaystyle \frac{3}{4}}X^2+{\displaystyle \frac{11}{8}}X2\zeta _3{\displaystyle \frac{1}{4}}\pi ^2\right){\displaystyle \frac{t^3}{u^3}}`$
$`)i\pi `$
$`D_{s++}^{gg,(2),[1]}`$ $`=`$ $`({\displaystyle \frac{2}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)X{\displaystyle \frac{7}{36}}X^3+{\displaystyle \frac{1}{3}}X^2Y+{\displaystyle \frac{43}{36}}X^2{\displaystyle \frac{13}{36}}X\pi ^2{\displaystyle \frac{5}{9}}X+{\displaystyle \frac{17}{12}}\zeta _3`$
$`{\displaystyle \frac{77}{108}}\pi ^2{\displaystyle \frac{1307}{432}})`$
$`+\left({\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{19}{18}}X{\displaystyle \frac{5}{36}}\pi ^2{\displaystyle \frac{37}{36}}\right){\displaystyle \frac{t}{u}}+\left({\displaystyle \frac{3}{4}}X^2{\displaystyle \frac{37}{18}}X+{\displaystyle \frac{5}{18}}\pi ^2\right){\displaystyle \frac{t^2}{u^2}}`$
$`+\left({\displaystyle \frac{7}{36}}X^3{\displaystyle \frac{37}{36}}X^2+{\displaystyle \frac{5}{18}}X\pi ^2\right){\displaystyle \frac{t^3}{u^3}}`$
$`+(({\displaystyle \frac{2}{3}}\mathrm{Li}_2(x){\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{2}{3}}XY{\displaystyle \frac{1}{18}}X{\displaystyle \frac{1}{12}}\pi ^2{\displaystyle \frac{13}{18}})`$
$`+({\displaystyle \frac{2}{3}}X{\displaystyle \frac{25}{18}}){\displaystyle \frac{t}{u}}+({\displaystyle \frac{13}{6}}X{\displaystyle \frac{37}{18}}){\displaystyle \frac{t^2}{u^2}}+({\displaystyle \frac{1}{4}}X^2{\displaystyle \frac{37}{18}}X){\displaystyle \frac{t^3}{u^3}})i\pi `$
$`E_{s++}^{gg,(2),[1]}`$ $`=`$ $`({\displaystyle \frac{2}{3}}\mathrm{Li}_3(y)+{\displaystyle \frac{2}{3}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{7}{36}}X^3{\displaystyle \frac{1}{3}}X^2Y+{\displaystyle \frac{2}{9}}X^2+{\displaystyle \frac{1}{2}}XY^2{\displaystyle \frac{5}{9}}XY+{\displaystyle \frac{4}{9}}X\pi ^2`$
$`{\displaystyle \frac{11}{6}}X{\displaystyle \frac{1}{6}}Y^3{\displaystyle \frac{17}{18}}Y^2+{\displaystyle \frac{7}{18}}Y\pi ^2+3Y+{\displaystyle \frac{23}{36}}\zeta _3+{\displaystyle \frac{37}{72}}\pi ^2+{\displaystyle \frac{3401}{1296}})`$
$`+(4\mathrm{Li}_4(x)+4\mathrm{Li}_4(y)+4\mathrm{Li}_4(z)+4\mathrm{Li}_3(x)Y+8\mathrm{Li}_3(y){\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)\pi ^28\mathrm{Li}_2(y)Y`$
$`+{\displaystyle \frac{1}{6}}X^4{\displaystyle \frac{2}{3}}X^3Y+{\displaystyle \frac{1}{2}}X^3+X^2Y^2{\displaystyle \frac{1}{4}}X^2Y+{\displaystyle \frac{1}{3}}X^2\pi ^2+{\displaystyle \frac{1}{9}}X^2{\displaystyle \frac{11}{3}}XY^2{\displaystyle \frac{2}{3}}XY\pi ^2`$
$`+{\displaystyle \frac{8}{9}}XY+X\pi ^2+{\displaystyle \frac{7}{9}}X2Y^24Y\zeta _3+{\displaystyle \frac{7}{6}}Y\pi ^2+{\displaystyle \frac{5}{3}}Y8\zeta _3+{\displaystyle \frac{7}{45}}\pi ^4{\displaystyle \frac{7}{12}}\pi ^2{\displaystyle \frac{19}{36}}){\displaystyle \frac{t}{u}}`$
$`+(12\mathrm{Li}_4(x)+12\mathrm{Li}_4(y)+12\mathrm{Li}_4(z)+12\mathrm{Li}_3(x)Y+8\mathrm{Li}_3(y)2\mathrm{Li}_2(x)\pi ^2`$
$`8\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{2}}X^42X^3Y+{\displaystyle \frac{1}{3}}X^3+3X^2Y^2+{\displaystyle \frac{1}{4}}X^2Y+X^2\pi ^2+{\displaystyle \frac{11}{12}}X^24XY^2`$
$`2XY\pi ^2{\displaystyle \frac{2}{3}}XY+{\displaystyle \frac{5}{6}}X\pi ^2+{\displaystyle \frac{11}{18}}X12Y\zeta _3+{\displaystyle \frac{4}{3}}Y\pi ^28\zeta _3+{\displaystyle \frac{7}{15}}\pi ^4{\displaystyle \frac{1}{2}}\pi ^2){\displaystyle \frac{t^2}{u^2}}`$
$`+(8\mathrm{Li}_4(x)+8\mathrm{Li}_4(y)+8\mathrm{Li}_4(z)+8\mathrm{Li}_3(x)Y{\displaystyle \frac{2}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)X`$
$`{\displaystyle \frac{4}{3}}\mathrm{Li}_2(x)\pi ^2+{\displaystyle \frac{1}{3}}X^4{\displaystyle \frac{4}{3}}X^3Y+{\displaystyle \frac{7}{36}}X^3+2X^2Y^2+{\displaystyle \frac{2}{3}}X^2\pi ^2{\displaystyle \frac{1}{36}}X^2{\displaystyle \frac{4}{3}}XY\pi ^2`$
$`+{\displaystyle \frac{5}{6}}X\pi ^28Y\zeta _3+{\displaystyle \frac{2}{3}}\zeta _3+{\displaystyle \frac{14}{45}}\pi ^4){\displaystyle \frac{t^3}{u^3}}{\displaystyle \frac{3}{2}}Y^2{\displaystyle \frac{s}{t}}`$
$`+(({\displaystyle \frac{2}{3}}\mathrm{Li}_2(x){\displaystyle \frac{1}{12}}X^2{\displaystyle \frac{1}{3}}XY{\displaystyle \frac{10}{9}}X{\displaystyle \frac{1}{3}}Y^2{\displaystyle \frac{22}{9}}Y+3)`$
$`+\left(4\mathrm{Li}_3(x)+8\mathrm{Li}_2(x)+{\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{1}{6}}XY{\displaystyle \frac{5}{9}}X+{\displaystyle \frac{1}{3}}Y^2{\displaystyle \frac{28}{9}}Y4\zeta _3+{\displaystyle \frac{19}{9}}\right){\displaystyle \frac{t}{u}}`$
$`+\left(12\mathrm{Li}_3(x)+8\mathrm{Li}_2(x)+{\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{1}{2}}X{\displaystyle \frac{2}{3}}Y12\zeta _3+{\displaystyle \frac{11}{18}}\right){\displaystyle \frac{t^2}{u^2}}`$
$`+(8\mathrm{Li}_3(x)+{\displaystyle \frac{2}{3}}\mathrm{Li}_2(x){\displaystyle \frac{1}{12}}X^2{\displaystyle \frac{1}{18}}X8\zeta _3{\displaystyle \frac{1}{9}}\pi ^2){\displaystyle \frac{t^3}{u^3}}3Y{\displaystyle \frac{s}{t}})i\pi `$
$`F_{s++}^{gg,(2),[1]}`$ $`=`$ $`\left({\displaystyle \frac{1}{18}}X+{\displaystyle \frac{5}{27}}\right)i\pi +\left({\displaystyle \frac{1}{12}}X^2+{\displaystyle \frac{5}{27}}X+{\displaystyle \frac{2}{27}}\pi ^2\right)`$
$`G_{s++}^{gg,(2),[1]}`$ $`=`$ $`\left({\displaystyle \frac{1}{36}}XY{\displaystyle \frac{1}{36}}Y^2\right)+\left({\displaystyle \frac{1}{36}}X^2{\displaystyle \frac{1}{36}}XY\right){\displaystyle \frac{t}{u}}`$
$`+\left(\left({\displaystyle \frac{1}{36}}X{\displaystyle \frac{1}{36}}Y\right)+\left({\displaystyle \frac{1}{36}}X{\displaystyle \frac{1}{36}}Y\right){\displaystyle \frac{t}{u}}\right)i\pi `$
$`A_{s++}^{gg,(2),[2]}`$ $`=`$ $`(2\mathrm{Li}_4(x)+\mathrm{Li}_4(y)+2\mathrm{Li}_4(z)+2\mathrm{Li}_3(x)Y+\mathrm{Li}_3(y)Y+{\displaystyle \frac{41}{3}}\mathrm{Li}_3(y){\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)\pi ^2`$
$`{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y^2{\displaystyle \frac{41}{3}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{12}}X^4{\displaystyle \frac{1}{3}}X^3Y+{\displaystyle \frac{1}{2}}X^2Y^2+{\displaystyle \frac{1}{6}}X^2\pi ^2{\displaystyle \frac{1}{6}}XY^3`$
$`{\displaystyle \frac{41}{6}}XY^2{\displaystyle \frac{1}{3}}XY\pi ^2{\displaystyle \frac{5}{24}}Y^4+{\displaystyle \frac{26}{9}}Y^3+Y^2\pi ^2{\displaystyle \frac{145}{18}}Y^2Y\zeta _3+{\displaystyle \frac{247}{72}}Y\pi ^2`$
$`{\displaystyle \frac{2}{27}}Y{\displaystyle \frac{445}{24}}\zeta _3+{\displaystyle \frac{121}{1440}}\pi ^4{\displaystyle \frac{13}{6}}\pi ^2+{\displaystyle \frac{30593}{1728}})`$
$`+(12\mathrm{Li}_4(x)+12\mathrm{Li}_4(y)+12\mathrm{Li}_4(z)+12\mathrm{Li}_3(x)Y+10\mathrm{Li}_3(y)2\mathrm{Li}_2(x)\pi ^2`$
$`10\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{2}}X^42X^3Y+3X^2Y^2+X^2\pi ^25XY^22XY\pi ^212Y\zeta _3`$
$`+{\displaystyle \frac{5}{3}}Y\pi ^210\zeta _3+{\displaystyle \frac{7}{15}}\pi ^4{\displaystyle \frac{5}{3}}\pi ^2+{\displaystyle \frac{1}{8}}){\displaystyle \frac{t}{u}}`$
$`+(10\mathrm{Li}_4(x)+10\mathrm{Li}_4(y)+10\mathrm{Li}_4(z)+10\mathrm{Li}_3(x)Y{\displaystyle \frac{5}{3}}\mathrm{Li}_2(x)\pi ^2+{\displaystyle \frac{5}{12}}X^4`$
$`{\displaystyle \frac{5}{3}}X^3Y+{\displaystyle \frac{5}{2}}X^2Y^2+{\displaystyle \frac{5}{6}}X^2\pi ^2{\displaystyle \frac{5}{3}}XY\pi ^210Y\zeta _3+{\displaystyle \frac{7}{18}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{3}{2}}\mathrm{Li}_4(y)+\mathrm{Li}_3(y)Y+3\mathrm{Li}_3(y){\displaystyle \frac{1}{4}}\mathrm{Li}_2(y)Y^23\mathrm{Li}_2(y)Y{\displaystyle \frac{3}{2}}XY^2`$
$`{\displaystyle \frac{7}{48}}Y^4+{\displaystyle \frac{13}{9}}Y^3+{\displaystyle \frac{2}{3}}Y^2\pi ^2+{\displaystyle \frac{713}{144}}Y^2Y\zeta _3+{\displaystyle \frac{23}{18}}Y\pi ^25Y3\zeta _3+{\displaystyle \frac{1}{60}}\pi ^4){\displaystyle \frac{s}{t}}`$
$`+{\displaystyle \frac{5}{2}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+((2\mathrm{Li}_3(x)+\mathrm{Li}_3(y)+{\displaystyle \frac{41}{3}}\mathrm{Li}_2(x)\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{2}}XY^2{\displaystyle \frac{5}{6}}Y^3+5Y^2+{\displaystyle \frac{1}{3}}Y\pi ^2`$
$`{\displaystyle \frac{35}{9}}Y\zeta _3{\displaystyle \frac{29}{72}}\pi ^2+{\displaystyle \frac{142}{27}})`$
$`+\left(12\mathrm{Li}_3(x)+10\mathrm{Li}_2(x)12\zeta _3\right){\displaystyle \frac{t}{u}}+5Y{\displaystyle \frac{s^2}{t^2}}+\left(10\mathrm{Li}_3(x)10\zeta _3\right){\displaystyle \frac{t^2}{u^2}}`$
$`+(\mathrm{Li}_3(y)+3\mathrm{Li}_2(x){\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y{\displaystyle \frac{7}{12}}Y^3+{\displaystyle \frac{5}{2}}Y^2+{\displaystyle \frac{1}{6}}Y\pi ^2+{\displaystyle \frac{713}{72}}Y\zeta _35){\displaystyle \frac{s}{t}})i\pi `$
$`B_{s++}^{gg,(2),[2]}`$ $`=`$ $`(6\mathrm{Li}_4(x)+15\mathrm{Li}_4(y)+20\mathrm{Li}_4(z)8\mathrm{Li}_3(x)X+14\mathrm{Li}_3(x)Y+\mathrm{Li}_3(x)`$
$`12\mathrm{Li}_3(y)X+5\mathrm{Li}_3(y)Y+7\mathrm{Li}_3(y)+\mathrm{Li}_2(x)X^2\mathrm{Li}_2(x)X+{\displaystyle \frac{7}{3}}\mathrm{Li}_2(x)\pi ^2`$
$`+6\mathrm{Li}_2(y)XY{\displaystyle \frac{5}{2}}\mathrm{Li}_2(y)Y^27\mathrm{Li}_2(y)Y+{\displaystyle \frac{3}{4}}X^43X^3Y+{\displaystyle \frac{1}{6}}X^3+{\displaystyle \frac{11}{2}}X^2Y^2`$
$`+{\displaystyle \frac{1}{3}}X^2Y+{\displaystyle \frac{5}{12}}X^2\pi ^2{\displaystyle \frac{41}{18}}X^2{\displaystyle \frac{5}{6}}XY^3{\displaystyle \frac{31}{6}}XY^2{\displaystyle \frac{4}{3}}XY\pi ^2+{\displaystyle \frac{19}{18}}XY+12X\zeta _3`$
$`{\displaystyle \frac{5}{6}}X\pi ^2+{\displaystyle \frac{3}{2}}X{\displaystyle \frac{1}{24}}Y^4+{\displaystyle \frac{5}{3}}Y^3+{\displaystyle \frac{5}{6}}Y^2\pi ^2{\displaystyle \frac{145}{36}}Y^214Y\zeta _3{\displaystyle \frac{5}{6}}Y\pi ^2+{\displaystyle \frac{71}{12}}Y`$
$`+{\displaystyle \frac{731}{72}}\zeta _3{\displaystyle \frac{383}{1440}}\pi ^4+{\displaystyle \frac{343}{144}}\pi ^2{\displaystyle \frac{18815}{1296}})`$
$`+(3\mathrm{Li}_4(x)+27\mathrm{Li}_4(y)+27\mathrm{Li}_4(z)18\mathrm{Li}_3(x)X+21\mathrm{Li}_3(x)Y{\displaystyle \frac{5}{2}}\mathrm{Li}_3(x)`$
$`12\mathrm{Li}_3(y)X+7\mathrm{Li}_3(y)+3\mathrm{Li}_2(x)X^2+{\displaystyle \frac{5}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{7}{2}}\mathrm{Li}_2(x)\pi ^2+6\mathrm{Li}_2(y)XY`$
$`7\mathrm{Li}_2(y)Y+{\displaystyle \frac{7}{8}}X^4{\displaystyle \frac{7}{2}}X^3Y{\displaystyle \frac{1}{2}}X^3+{\displaystyle \frac{27}{4}}X^2Y^2+{\displaystyle \frac{9}{4}}X^2YX^2\pi ^2{\displaystyle \frac{13}{8}}X^2`$
$`{\displaystyle \frac{7}{2}}XY^2{\displaystyle \frac{1}{2}}XY\pi ^23XY+18X\zeta _3{\displaystyle \frac{5}{12}}X\pi ^2+{\displaystyle \frac{3}{4}}X21Y\zeta _3{\displaystyle \frac{5}{6}}Y\pi ^2`$
$`+{\displaystyle \frac{15}{2}}\zeta _3{\displaystyle \frac{37}{60}}\pi ^4+{\displaystyle \frac{5}{4}}\pi ^2+{\displaystyle \frac{1}{4}}){\displaystyle \frac{t}{u}}`$
$`+(5\mathrm{Li}_4(x)+7\mathrm{Li}_4(y)+7\mathrm{Li}_4(z)6\mathrm{Li}_3(x)X+7\mathrm{Li}_3(x)Y{\displaystyle \frac{7}{2}}\mathrm{Li}_3(x)+{\displaystyle \frac{7}{2}}\mathrm{Li}_2(x)X`$
$`+{\displaystyle \frac{5}{6}}\mathrm{Li}_2(x)\pi ^2+{\displaystyle \frac{1}{24}}X^4{\displaystyle \frac{7}{6}}X^3Y+X^3+{\displaystyle \frac{7}{4}}X^2Y^2+{\displaystyle \frac{1}{4}}X^2Y{\displaystyle \frac{5}{12}}X^2\pi ^2{\displaystyle \frac{11}{8}}X^2`$
$`{\displaystyle \frac{7}{6}}XY\pi ^2+6X\zeta _3+{\displaystyle \frac{5}{12}}X\pi ^27Y\zeta _3+{\displaystyle \frac{7}{2}}\zeta _3{\displaystyle \frac{7}{36}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+\left({\displaystyle \frac{3}{8}}X^2{\displaystyle \frac{3}{4}}XY+{\displaystyle \frac{3}{8}}Y^2+{\displaystyle \frac{3}{8}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`+(3\mathrm{Li}_4(x)7\mathrm{Li}_4(y)+\mathrm{Li}_4(z)2\mathrm{Li}_3(x)Y+6\mathrm{Li}_3(x)+\mathrm{Li}_3(y)X+4\mathrm{Li}_3(y)Y`$
$`+3\mathrm{Li}_3(y)2\mathrm{Li}_2(x)X^26\mathrm{Li}_2(x)X{\displaystyle \frac{1}{6}}\mathrm{Li}_2(x)\pi ^24\mathrm{Li}_2(y)XY3\mathrm{Li}_2(y)Y`$
$`+{\displaystyle \frac{1}{24}}X^4{\displaystyle \frac{5}{3}}X^3Y{\displaystyle \frac{7}{4}}X^2Y^2+{\displaystyle \frac{1}{3}}X^2Y+{\displaystyle \frac{4}{3}}X^2\pi ^2+{\displaystyle \frac{5}{6}}XY^35XY^2{\displaystyle \frac{5}{6}}XY\pi ^2`$
$`+{\displaystyle \frac{25}{9}}XYX\zeta _3{\displaystyle \frac{5}{3}}X\pi ^2{\displaystyle \frac{1}{6}}Y^4+{\displaystyle \frac{28}{9}}Y^3+{\displaystyle \frac{1}{6}}Y^2\pi ^2{\displaystyle \frac{199}{72}}Y^2`$
$`+2Y\zeta _3+{\displaystyle \frac{35}{18}}Y\pi ^22Y3\zeta _3+{\displaystyle \frac{7}{72}}\pi ^4){\displaystyle \frac{s}{t}}+Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+((6\mathrm{Li}_3(x)7\mathrm{Li}_3(y)4\mathrm{Li}_2(x)X+6\mathrm{Li}_2(x)+\mathrm{Li}_2(y)Y+{\displaystyle \frac{2}{3}}X^3X^2Y+{\displaystyle \frac{4}{3}}X^2`$
$`{\displaystyle \frac{3}{2}}XY^2{\displaystyle \frac{8}{3}}XY+{\displaystyle \frac{1}{2}}X\pi ^2{\displaystyle \frac{7}{2}}X{\displaystyle \frac{1}{6}}Y^3+{\displaystyle \frac{10}{3}}Y^27Y2\zeta _3{\displaystyle \frac{5}{6}}\pi ^2{\displaystyle \frac{29}{4}})`$
$`+(3\mathrm{Li}_3(x)12\mathrm{Li}_3(y)+{\displaystyle \frac{19}{2}}\mathrm{Li}_2(x)+6\mathrm{Li}_2(y)Y+X^3{\displaystyle \frac{1}{2}}X^2+{\displaystyle \frac{9}{2}}XY+{\displaystyle \frac{1}{2}}X\pi ^2`$
$`{\displaystyle \frac{25}{4}}X3Y3\zeta _3{\displaystyle \frac{17}{12}}\pi ^2+{\displaystyle \frac{3}{4}}){\displaystyle \frac{t}{u}}`$
$`+\left(\mathrm{Li}_3(x)+{\displaystyle \frac{7}{2}}\mathrm{Li}_2(x)+{\displaystyle \frac{3}{2}}X^2+{\displaystyle \frac{1}{2}}XY{\displaystyle \frac{11}{4}}X\zeta _3{\displaystyle \frac{7}{12}}\pi ^2\right){\displaystyle \frac{t^2}{u^2}}`$
$`+(2\mathrm{Li}_3(x)+5\mathrm{Li}_3(y)3\mathrm{Li}_2(x)4\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{6}}X^33X^2Y+{\displaystyle \frac{10}{3}}X^2+2XY^2`$
$`{\displaystyle \frac{19}{3}}XY+{\displaystyle \frac{5}{6}}X\pi ^2+{\displaystyle \frac{25}{9}}X{\displaystyle \frac{1}{2}}Y^3+4Y^2{\displaystyle \frac{2}{3}}Y\pi ^2{\displaystyle \frac{11}{4}}Y+\zeta _3+2\pi ^22){\displaystyle \frac{s}{t}}`$
$`+2Y{\displaystyle \frac{s^2}{t^2}})i\pi `$
$`C_{s++}^{gg,(2),[2]}`$ $`=`$ $`(2\mathrm{Li}_4(x)+2\mathrm{Li}_4(y)+2\mathrm{Li}_4(z)+2\mathrm{Li}_3(x)Y+\mathrm{Li}_3(x)+\mathrm{Li}_3(y)\mathrm{Li}_2(x)X`$
$`{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)\pi ^2\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{12}}X^4{\displaystyle \frac{1}{3}}X^3Y+{\displaystyle \frac{1}{2}}X^2Y^2{\displaystyle \frac{1}{2}}X^2Y+{\displaystyle \frac{1}{6}}X^2\pi ^2+{\displaystyle \frac{1}{2}}X^2`$
$`{\displaystyle \frac{1}{2}}XY^2{\displaystyle \frac{1}{3}}XY\pi ^2{\displaystyle \frac{3}{2}}XY+{\displaystyle \frac{1}{6}}X\pi ^2+{\displaystyle \frac{3}{2}}X+{\displaystyle \frac{3}{4}}Y^22Y\zeta _3+{\displaystyle \frac{1}{6}}Y\pi ^2{\displaystyle \frac{3}{4}}Y`$
$`+{\displaystyle \frac{7}{4}}\zeta _3+{\displaystyle \frac{13}{120}}\pi ^4{\displaystyle \frac{5}{48}}\pi ^2{\displaystyle \frac{247}{64}})`$
$`+(3\mathrm{Li}_4(x)+3\mathrm{Li}_4(y)+3\mathrm{Li}_4(z)+3\mathrm{Li}_3(x)Y+{\displaystyle \frac{1}{2}}\mathrm{Li}_3(x)+\mathrm{Li}_3(y){\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X`$
$`{\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)\pi ^2\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{8}}X^4{\displaystyle \frac{1}{2}}X^3Y+{\displaystyle \frac{3}{4}}X^2Y^2{\displaystyle \frac{1}{4}}X^2Y+{\displaystyle \frac{1}{4}}X^2\pi ^2{\displaystyle \frac{1}{8}}X^2`$
$`{\displaystyle \frac{1}{2}}XY^2{\displaystyle \frac{1}{2}}XY\pi ^2+{\displaystyle \frac{1}{12}}X\pi ^2+{\displaystyle \frac{3}{4}}X3Y\zeta _3+{\displaystyle \frac{1}{6}}Y\pi ^2{\displaystyle \frac{3}{2}}\zeta _3+{\displaystyle \frac{7}{60}}\pi ^4{\displaystyle \frac{1}{4}}\pi ^2+{\displaystyle \frac{1}{8}}){\displaystyle \frac{t}{u}}`$
$`+(\mathrm{Li}_4(x)+\mathrm{Li}_4(y)+\mathrm{Li}_4(z)+\mathrm{Li}_3(x)Y{\displaystyle \frac{1}{2}}\mathrm{Li}_3(x)+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{1}{6}}\mathrm{Li}_2(x)\pi ^2`$
$`+{\displaystyle \frac{1}{24}}X^4{\displaystyle \frac{1}{6}}X^3Y+{\displaystyle \frac{1}{4}}X^2Y^2+{\displaystyle \frac{1}{4}}X^2Y+{\displaystyle \frac{1}{12}}X^2\pi ^2+{\displaystyle \frac{1}{8}}X^2{\displaystyle \frac{1}{6}}XY\pi ^2{\displaystyle \frac{1}{12}}X\pi ^2`$
$`Y\zeta _3+{\displaystyle \frac{1}{2}}\zeta _3+{\displaystyle \frac{7}{180}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+\left({\displaystyle \frac{3}{8}}X^2{\displaystyle \frac{3}{4}}XY+{\displaystyle \frac{3}{8}}Y^2+{\displaystyle \frac{3}{8}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`+(\mathrm{Li}_4(x)+{\displaystyle \frac{5}{2}}\mathrm{Li}_4(y)+\mathrm{Li}_4(z)\mathrm{Li}_3(y)X\mathrm{Li}_3(y)Y3\mathrm{Li}_3(y){\displaystyle \frac{1}{6}}\mathrm{Li}_2(x)\pi ^2`$
$`+{\displaystyle \frac{1}{4}}\mathrm{Li}_2(y)Y^2+3\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{24}}X^4{\displaystyle \frac{1}{6}}X^3Y+{\displaystyle \frac{1}{12}}X^2\pi ^2+{\displaystyle \frac{3}{2}}XY^2{\displaystyle \frac{1}{6}}XY\pi ^2`$
$`+X\zeta _3+{\displaystyle \frac{1}{16}}Y^4{\displaystyle \frac{1}{4}}Y^2\pi ^2{\displaystyle \frac{15}{16}}Y^2+Y\zeta _3{\displaystyle \frac{1}{2}}Y\pi ^2+Y+3\zeta _3{\displaystyle \frac{1}{120}}\pi ^4){\displaystyle \frac{s}{t}}`$
$`{\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{s^2}{t^2}}+((2\mathrm{Li}_3(x)XY{\displaystyle \frac{1}{2}}X2\zeta _3+{\displaystyle \frac{1}{6}}\pi ^2+{\displaystyle \frac{3}{4}})`$
$`+\left(3\mathrm{Li}_3(x)+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(x){\displaystyle \frac{1}{2}}XY{\displaystyle \frac{1}{4}}X3\zeta _3+{\displaystyle \frac{1}{12}}\pi ^2+{\displaystyle \frac{3}{4}}\right){\displaystyle \frac{t}{u}}`$
$`+\left(\mathrm{Li}_3(x)+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)+{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{1}{4}}X\zeta _3{\displaystyle \frac{1}{12}}\pi ^2\right){\displaystyle \frac{t^2}{u^2}}`$
$`+\left(2\mathrm{Li}_3(y)3\mathrm{Li}_2(x)+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{4}}Y^3{\displaystyle \frac{15}{8}}Y+2\zeta _3+1\right){\displaystyle \frac{s}{t}}`$
$`Y{\displaystyle \frac{s^2}{t^2}})i\pi `$
$`D_{s++}^{gg,(2),[2]}`$ $`=`$ $`(2\mathrm{Li}_4(x)2\mathrm{Li}_4(y)2\mathrm{Li}_4(z)2\mathrm{Li}_3(x)Y{\displaystyle \frac{14}{3}}\mathrm{Li}_3(y)+{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)\pi ^2+{\displaystyle \frac{14}{3}}\mathrm{Li}_2(y)Y`$
$`{\displaystyle \frac{1}{12}}X^4+{\displaystyle \frac{1}{3}}X^3Y{\displaystyle \frac{1}{2}}X^2Y^2{\displaystyle \frac{1}{6}}X^2\pi ^2+{\displaystyle \frac{7}{3}}XY^2+{\displaystyle \frac{1}{3}}XY\pi ^2{\displaystyle \frac{7}{18}}Y^3+{\displaystyle \frac{71}{36}}Y^2`$
$`+2Y\zeta _3{\displaystyle \frac{47}{36}}Y\pi ^2{\displaystyle \frac{14}{9}}Y+{\displaystyle \frac{65}{12}}\zeta _3{\displaystyle \frac{7}{90}}\pi ^4+{\displaystyle \frac{10}{27}}\pi ^2{\displaystyle \frac{863}{432}})`$
$`+(6\mathrm{Li}_4(x)6\mathrm{Li}_4(y)6\mathrm{Li}_4(z)6\mathrm{Li}_3(x)Y4\mathrm{Li}_3(y)+\mathrm{Li}_2(x)\pi ^2+4\mathrm{Li}_2(y)Y`$
$`{\displaystyle \frac{1}{4}}X^4+X^3Y{\displaystyle \frac{3}{2}}X^2Y^2{\displaystyle \frac{1}{2}}X^2\pi ^2+2XY^2+XY\pi ^2+6Y\zeta _3{\displaystyle \frac{2}{3}}Y\pi ^2+4\zeta _3`$
$`{\displaystyle \frac{7}{30}}\pi ^4+{\displaystyle \frac{2}{3}}\pi ^2){\displaystyle \frac{t}{u}}`$
$`+(4\mathrm{Li}_4(x)4\mathrm{Li}_4(y)4\mathrm{Li}_4(z)4\mathrm{Li}_3(x)Y+{\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)\pi ^2{\displaystyle \frac{1}{6}}X^4+{\displaystyle \frac{2}{3}}X^3Y`$
$`X^2Y^2{\displaystyle \frac{1}{3}}X^2\pi ^2+{\displaystyle \frac{2}{3}}XY\pi ^2+4Y\zeta _3{\displaystyle \frac{7}{45}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+\left({\displaystyle \frac{7}{36}}Y^3{\displaystyle \frac{71}{36}}Y^2{\displaystyle \frac{5}{18}}Y\pi ^2+2Y\right){\displaystyle \frac{s}{t}}`$
$`Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+((2\mathrm{Li}_3(x){\displaystyle \frac{14}{3}}\mathrm{Li}_2(x){\displaystyle \frac{1}{2}}Y^2+{\displaystyle \frac{1}{2}}Y+2\zeta _3+{\displaystyle \frac{1}{36}}\pi ^2{\displaystyle \frac{25}{18}})`$
$`+\left(6\mathrm{Li}_3(x)4\mathrm{Li}_2(x)+6\zeta _3\right){\displaystyle \frac{t}{u}}+\left(4\mathrm{Li}_3(x)+4\zeta _3\right){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{71}{18}}Y+2){\displaystyle \frac{s}{t}}2Y{\displaystyle \frac{s^2}{t^2}})i\pi `$
$`E_{s++}^{gg,(2),[2]}`$ $`=`$ $`(4\mathrm{Li}_4(x)+4\mathrm{Li}_4(y)+4\mathrm{Li}_4(z)+4\mathrm{Li}_3(x)Y+{\displaystyle \frac{26}{3}}\mathrm{Li}_3(y){\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)\pi ^2{\displaystyle \frac{26}{3}}\mathrm{Li}_2(y)Y`$
$`+{\displaystyle \frac{1}{6}}X^4{\displaystyle \frac{2}{3}}X^3Y+{\displaystyle \frac{1}{12}}X^3+X^2Y^2{\displaystyle \frac{1}{12}}X^2Y+{\displaystyle \frac{1}{3}}X^2\pi ^2+{\displaystyle \frac{14}{9}}X^2{\displaystyle \frac{23}{6}}XY^2`$
$`{\displaystyle \frac{2}{3}}XY\pi ^2{\displaystyle \frac{7}{9}}XY+3X{\displaystyle \frac{1}{6}}Y^3+{\displaystyle \frac{1}{18}}Y^24Y\zeta _3+{\displaystyle \frac{7}{9}}Y\pi ^22Y{\displaystyle \frac{313}{36}}\zeta _3`$
$`+{\displaystyle \frac{7}{45}}\pi ^4{\displaystyle \frac{85}{72}}\pi ^2+{\displaystyle \frac{4085}{1296}})`$
$`+(12\mathrm{Li}_4(x)+12\mathrm{Li}_4(y)+12\mathrm{Li}_4(z)+12\mathrm{Li}_3(x)Y+8\mathrm{Li}_3(y)2\mathrm{Li}_2(x)\pi ^2`$
$`8\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{2}}X^42X^3Y+{\displaystyle \frac{1}{4}}X^3+3X^2Y^2{\displaystyle \frac{1}{4}}X^2Y+X^2\pi ^2+2X^24XY^2`$
$`2XY\pi ^2{\displaystyle \frac{2}{3}}XY+{\displaystyle \frac{5}{3}}X12Y\zeta _3+{\displaystyle \frac{4}{3}}Y\pi ^28\zeta _3+{\displaystyle \frac{7}{15}}\pi ^4{\displaystyle \frac{7}{9}}\pi ^2){\displaystyle \frac{t}{u}}`$
$`+(8\mathrm{Li}_4(x)+8\mathrm{Li}_4(y)+8\mathrm{Li}_4(z)+8\mathrm{Li}_3(x)Y{\displaystyle \frac{2}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)X`$
$`{\displaystyle \frac{4}{3}}\mathrm{Li}_2(x)\pi ^2+{\displaystyle \frac{1}{3}}X^4{\displaystyle \frac{4}{3}}X^3Y+2X^2Y^2+{\displaystyle \frac{2}{3}}X^2\pi ^2+{\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{4}{3}}XY\pi ^2+{\displaystyle \frac{5}{9}}X\pi ^2`$
$`8Y\zeta _3+{\displaystyle \frac{2}{3}}\zeta _3+{\displaystyle \frac{14}{45}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{2}{3}}\mathrm{Li}_3(y){\displaystyle \frac{2}{3}}\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{3}}X^2Y+{\displaystyle \frac{1}{2}}XY^2{\displaystyle \frac{5}{9}}XY+{\displaystyle \frac{1}{6}}X\pi ^2{\displaystyle \frac{13}{36}}Y^3`$
$`+{\displaystyle \frac{49}{12}}Y^2Y\pi ^23Y{\displaystyle \frac{2}{3}}\zeta _3){\displaystyle \frac{s}{t}}`$
$`+{\displaystyle \frac{3}{2}}Y^2{\displaystyle \frac{s^2}{t^2}}`$
$`+((4\mathrm{Li}_3(x)+{\displaystyle \frac{26}{3}}\mathrm{Li}_2(x)+{\displaystyle \frac{1}{6}}X^2+{\displaystyle \frac{5}{6}}XY+{\displaystyle \frac{7}{3}}X{\displaystyle \frac{2}{3}}Y4\zeta _3+{\displaystyle \frac{19}{6}})`$
$`+\left(12\mathrm{Li}_3(x)+8\mathrm{Li}_2(x)+{\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{10}{3}}X{\displaystyle \frac{2}{3}}Y12\zeta _3+{\displaystyle \frac{5}{3}}\right){\displaystyle \frac{t}{u}}`$
$`+\left(8\mathrm{Li}_3(x)+{\displaystyle \frac{2}{3}}\mathrm{Li}_2(x){\displaystyle \frac{1}{3}}X^2+X8\zeta _3{\displaystyle \frac{1}{9}}\pi ^2\right){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{2}{3}}\mathrm{Li}_2(x){\displaystyle \frac{1}{3}}X^2+XY{\displaystyle \frac{5}{9}}X+{\displaystyle \frac{1}{12}}Y^2+{\displaystyle \frac{137}{18}}Y3){\displaystyle \frac{s}{t}}+3Y{\displaystyle \frac{s^2}{t^2}})i\pi `$
$`F_{s++}^{gg,(2),[2]}`$ $`=`$ $`\left({\displaystyle \frac{1}{18}}Y+{\displaystyle \frac{5}{27}}\right)i\pi +\left({\displaystyle \frac{1}{12}}Y^2+{\displaystyle \frac{5}{27}}Y+{\displaystyle \frac{2}{27}}\pi ^2\right)`$
$`G_{s++}^{gg,(2),[2]}`$ $`=`$ $`\left({\displaystyle \frac{1}{36}}X^2+{\displaystyle \frac{1}{18}}XY{\displaystyle \frac{1}{36}}Y^2\right)+\left({\displaystyle \frac{1}{36}}X{\displaystyle \frac{1}{36}}Y\right)i\pi {\displaystyle \frac{s}{t}}+\left({\displaystyle \frac{1}{36}}XY{\displaystyle \frac{1}{36}}Y^2\right){\displaystyle \frac{s}{t}}`$
$`H_{s++}^{gg,(2),[3]}`$ $`=`$ $`(4\mathrm{Li}_4(x){\displaystyle \frac{13}{2}}\mathrm{Li}_4(y)+4\mathrm{Li}_4(z)+\mathrm{Li}_3(x)X+{\displaystyle \frac{7}{2}}\mathrm{Li}_3(x)Y{\displaystyle \frac{13}{6}}\mathrm{Li}_3(x){\displaystyle \frac{5}{2}}\mathrm{Li}_3(y)X`$
$`+{\displaystyle \frac{15}{2}}\mathrm{Li}_3(y)Y+{\displaystyle \frac{9}{4}}\mathrm{Li}_3(y)\mathrm{Li}_2(x)X^2+{\displaystyle \frac{13}{6}}\mathrm{Li}_2(x)X+{\displaystyle \frac{11}{6}}\mathrm{Li}_2(x)\pi ^2+2\mathrm{Li}_2(y)XY`$
$`{\displaystyle \frac{9}{4}}\mathrm{Li}_2(y)Y^2{\displaystyle \frac{9}{4}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{6}}X^4X^3Y+{\displaystyle \frac{7}{4}}X^2Y^2+{\displaystyle \frac{4}{3}}X^2Y{\displaystyle \frac{1}{6}}X^2\pi ^2{\displaystyle \frac{1}{4}}XY^3`$
$`{\displaystyle \frac{23}{12}}XY^2+{\displaystyle \frac{11}{6}}XY+{\displaystyle \frac{5}{2}}X\zeta _3{\displaystyle \frac{91}{24}}X\pi ^2{\displaystyle \frac{1}{16}}Y^4+{\displaystyle \frac{11}{18}}Y^3{\displaystyle \frac{1}{4}}Y^2\pi ^2+{\displaystyle \frac{161}{72}}Y^23Y\zeta _3`$
$`+{\displaystyle \frac{43}{48}}Y\pi ^2{\displaystyle \frac{211}{108}}Y{\displaystyle \frac{9}{4}}\zeta _3+{\displaystyle \frac{43}{120}}\pi ^4{\displaystyle \frac{391}{72}}\pi ^2)`$
$`+({\displaystyle \frac{25}{4}}\mathrm{Li}_4(x){\displaystyle \frac{7}{2}}\mathrm{Li}_3(x)X\mathrm{Li}_3(x)Y{\displaystyle \frac{13}{3}}\mathrm{Li}_3(x)2\mathrm{Li}_3(y)X{\displaystyle \frac{29}{2}}\mathrm{Li}_3(y)`$
$`+{\displaystyle \frac{3}{8}}\mathrm{Li}_2(x)X^2+{\displaystyle \frac{13}{3}}\mathrm{Li}_2(x)X+{\displaystyle \frac{5}{3}}\mathrm{Li}_2(x)\pi ^2+\mathrm{Li}_2(y)XY+{\displaystyle \frac{29}{2}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{3}{32}}X^4`$
$`+{\displaystyle \frac{1}{6}}X^3Y{\displaystyle \frac{11}{12}}X^3+{\displaystyle \frac{7}{6}}X^2Y{\displaystyle \frac{13}{8}}X^2\pi ^2+{\displaystyle \frac{25}{24}}X^2+{\displaystyle \frac{29}{4}}XY^2+{\displaystyle \frac{2}{3}}XY\pi ^2+{\displaystyle \frac{17}{4}}XY`$
$`+{\displaystyle \frac{7}{2}}X\zeta _3{\displaystyle \frac{23}{9}}X\pi ^2{\displaystyle \frac{5}{6}}X+Y\zeta _3{\displaystyle \frac{25}{12}}Y\pi ^2+{\displaystyle \frac{46}{3}}\zeta _3{\displaystyle \frac{25}{72}}\pi ^4{\displaystyle \frac{77}{72}}\pi ^2){\displaystyle \frac{t}{u}}`$
$`+(11\mathrm{Li}_4(x){\displaystyle \frac{33}{2}}\mathrm{Li}_4(y){\displaystyle \frac{33}{2}}\mathrm{Li}_4(z)+{\displaystyle \frac{7}{2}}\mathrm{Li}_3(x)X{\displaystyle \frac{31}{2}}\mathrm{Li}_3(x)Y+{\displaystyle \frac{31}{12}}\mathrm{Li}_3(x)`$
$`+2\mathrm{Li}_3(y)X{\displaystyle \frac{3}{4}}\mathrm{Li}_2(x)X^2{\displaystyle \frac{31}{12}}\mathrm{Li}_2(x)X+{\displaystyle \frac{19}{12}}\mathrm{Li}_2(x)\pi ^2\mathrm{Li}_2(y)XY{\displaystyle \frac{11}{16}}X^4`$
$`+{\displaystyle \frac{31}{12}}X^3Y{\displaystyle \frac{1}{12}}X^3{\displaystyle \frac{33}{8}}X^2Y^2+{\displaystyle \frac{1}{3}}X^2Y{\displaystyle \frac{7}{8}}X^2\pi ^2+{\displaystyle \frac{15}{8}}X^2+{\displaystyle \frac{25}{12}}XY\pi ^2`$
$`{\displaystyle \frac{7}{2}}X\zeta _3{\displaystyle \frac{191}{72}}X\pi ^2+{\displaystyle \frac{31}{2}}Y\zeta _3{\displaystyle \frac{31}{12}}\zeta _3{\displaystyle \frac{139}{360}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{5}{2}}\mathrm{Li}_4(x){\displaystyle \frac{5}{2}}\mathrm{Li}_4(y)+8\mathrm{Li}_4(z){\displaystyle \frac{13}{2}}\mathrm{Li}_3(x)X+6\mathrm{Li}_3(x)Y{\displaystyle \frac{31}{6}}\mathrm{Li}_3(x)`$
$`6\mathrm{Li}_3(y)X+{\displaystyle \frac{13}{2}}\mathrm{Li}_3(y)Y+{\displaystyle \frac{31}{6}}\mathrm{Li}_3(y)+{\displaystyle \frac{5}{4}}\mathrm{Li}_2(x)X^2+{\displaystyle \frac{31}{6}}\mathrm{Li}_2(x)X+{\displaystyle \frac{11}{3}}\mathrm{Li}_2(x)\pi ^2`$
$`+4\mathrm{Li}_2(y)XY{\displaystyle \frac{5}{4}}\mathrm{Li}_2(y)Y^2{\displaystyle \frac{31}{6}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{19}{48}}X^4{\displaystyle \frac{17}{12}}X^3Y{\displaystyle \frac{11}{18}}X^3+3X^2Y^2`$
$`+{\displaystyle \frac{7}{4}}X^2Y+{\displaystyle \frac{5}{12}}X^2\pi ^2+{\displaystyle \frac{43}{18}}X^2+{\displaystyle \frac{1}{12}}XY^3{\displaystyle \frac{7}{4}}XY^2+{\displaystyle \frac{5}{6}}XY\pi ^2+{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{11}{2}}X\zeta _3`$
$`{\displaystyle \frac{23}{16}}X\pi ^2+{\displaystyle \frac{19}{27}}X{\displaystyle \frac{1}{16}}Y^4+{\displaystyle \frac{11}{18}}Y^3+{\displaystyle \frac{1}{4}}Y^2\pi ^2{\displaystyle \frac{26}{9}}Y^2{\displaystyle \frac{11}{2}}Y\zeta _3+{\displaystyle \frac{23}{16}}Y\pi ^2`$
$`{\displaystyle \frac{19}{27}}Y{\displaystyle \frac{41}{180}}\pi ^4{\displaystyle \frac{1}{4}}\pi ^2){\displaystyle \frac{t}{s}}+{\displaystyle \frac{3}{2}}Y^2{\displaystyle \frac{s}{t}}`$
$`+(({\displaystyle \frac{9}{2}}\mathrm{Li}_3(x)+5\mathrm{Li}_3(y)4\mathrm{Li}_2(x)X+{\displaystyle \frac{53}{12}}\mathrm{Li}_2(x){\displaystyle \frac{5}{2}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{6}}X^3{\displaystyle \frac{5}{4}}X^2Y`$
$`+{\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{19}{4}}XY{\displaystyle \frac{1}{2}}X\pi ^2+{\displaystyle \frac{11}{6}}X{\displaystyle \frac{1}{4}}Y^3{\displaystyle \frac{19}{24}}Y^2{\displaystyle \frac{3}{2}}Y\pi ^2+{\displaystyle \frac{149}{12}}Y{\displaystyle \frac{1}{2}}\zeta _3`$
$`{\displaystyle \frac{187}{144}}\pi ^2{\displaystyle \frac{211}{108}})`$
$`+({\displaystyle \frac{9}{2}}\mathrm{Li}_3(x)2\mathrm{Li}_3(y){\displaystyle \frac{1}{4}}\mathrm{Li}_2(x)X{\displaystyle \frac{61}{6}}\mathrm{Li}_2(x)+\mathrm{Li}_2(y)Y+{\displaystyle \frac{7}{8}}X^3X^2+{\displaystyle \frac{7}{3}}XY`$
$`{\displaystyle \frac{1}{3}}X\pi ^2+{\displaystyle \frac{19}{3}}X+{\displaystyle \frac{17}{4}}Y+{\displaystyle \frac{9}{2}}\zeta _3{\displaystyle \frac{5}{9}}\pi ^2{\displaystyle \frac{5}{6}}){\displaystyle \frac{t}{u}}`$
$`+(12\mathrm{Li}_3(x)+2\mathrm{Li}_3(y){\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{31}{12}}\mathrm{Li}_2(x)\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{4}}X^3+{\displaystyle \frac{11}{8}}X^2`$
$`+{\displaystyle \frac{2}{3}}XY{\displaystyle \frac{1}{6}}X\pi ^2+{\displaystyle \frac{15}{4}}X+12\zeta _3+{\displaystyle \frac{31}{72}}\pi ^2){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{1}{2}}\mathrm{Li}_3(x)+{\displaystyle \frac{1}{2}}\mathrm{Li}_3(y){\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{31}{3}}\mathrm{Li}_2(x)+{\displaystyle \frac{3}{2}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{5}{12}}X^3`$
$`{\displaystyle \frac{5}{4}}X^2Y{\displaystyle \frac{5}{6}}X^2+{\displaystyle \frac{5}{4}}XY^2+{\displaystyle \frac{31}{6}}XY+X\pi ^2{\displaystyle \frac{5}{6}}X{\displaystyle \frac{5}{12}}Y^3+{\displaystyle \frac{5}{6}}Y^2Y\pi ^2`$
$`+{\displaystyle \frac{5}{6}}Y{\displaystyle \frac{31}{36}}\pi ^2){\displaystyle \frac{t}{s}}+3Y{\displaystyle \frac{s}{t}})i\pi `$
$`I_{s++}^{gg,(2),[3]}`$ $`=`$ $`(2\mathrm{Li}_4(x){\displaystyle \frac{5}{2}}\mathrm{Li}_4(y)+\mathrm{Li}_3(x)X+{\displaystyle \frac{1}{2}}\mathrm{Li}_3(x)Y{\displaystyle \frac{3}{2}}\mathrm{Li}_3(x)+{\displaystyle \frac{1}{2}}\mathrm{Li}_3(y)X+{\displaystyle \frac{3}{2}}\mathrm{Li}_3(y)Y`$
$`+{\displaystyle \frac{9}{4}}\mathrm{Li}_3(y)+{\displaystyle \frac{3}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{1}{6}}\mathrm{Li}_2(x)\pi ^2{\displaystyle \frac{1}{4}}\mathrm{Li}_2(y)Y^2{\displaystyle \frac{9}{4}}\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{12}}X^3Y`$
$`+{\displaystyle \frac{3}{8}}X^2Y^2{\displaystyle \frac{5}{12}}X^2Y+{\displaystyle \frac{5}{8}}X^2\pi ^2{\displaystyle \frac{1}{6}}XY^3{\displaystyle \frac{7}{6}}XY^2{\displaystyle \frac{1}{3}}XY\pi ^2+{\displaystyle \frac{1}{2}}XY`$
$`{\displaystyle \frac{1}{2}}X\zeta _3+{\displaystyle \frac{47}{24}}X\pi ^2+{\displaystyle \frac{1}{48}}Y^4+{\displaystyle \frac{13}{24}}Y^2\pi ^2{\displaystyle \frac{9}{8}}Y^2{\displaystyle \frac{1}{2}}Y\zeta _3+{\displaystyle \frac{17}{24}}Y\pi ^2+{\displaystyle \frac{47}{12}}Y{\displaystyle \frac{9}{4}}\zeta _3`$
$`+{\displaystyle \frac{1}{36}}\pi ^4+{\displaystyle \frac{3}{8}}\pi ^2)`$
$`+({\displaystyle \frac{23}{4}}\mathrm{Li}_4(x)3\mathrm{Li}_4(y)3\mathrm{Li}_4(z)+{\displaystyle \frac{11}{2}}\mathrm{Li}_3(x)X\mathrm{Li}_3(x)Y4\mathrm{Li}_3(x)+4\mathrm{Li}_3(y)X`$
$`+{\displaystyle \frac{1}{2}}\mathrm{Li}_3(y){\displaystyle \frac{9}{8}}\mathrm{Li}_2(x)X^2+4\mathrm{Li}_2(x)X{\displaystyle \frac{11}{6}}\mathrm{Li}_2(x)\pi ^22\mathrm{Li}_2(y)XY{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y`$
$`{\displaystyle \frac{5}{32}}X^4+{\displaystyle \frac{1}{6}}X^3Y{\displaystyle \frac{3}{4}}X^2Y^2{\displaystyle \frac{1}{12}}X^2Y+X^2\pi ^2+{\displaystyle \frac{13}{8}}X^2{\displaystyle \frac{1}{4}}XY^2{\displaystyle \frac{5}{6}}XY\pi ^2`$
$`{\displaystyle \frac{7}{12}}XY{\displaystyle \frac{11}{2}}X\zeta _3+{\displaystyle \frac{10}{3}}X\pi ^2+{\displaystyle \frac{1}{2}}X+Y\zeta _3+{\displaystyle \frac{5}{12}}Y\pi ^2+{\displaystyle \frac{133}{360}}\pi ^4+{\displaystyle \frac{11}{24}}\pi ^2){\displaystyle \frac{t}{u}}`$
$`+(4\mathrm{Li}_4(x){\displaystyle \frac{3}{2}}\mathrm{Li}_4(y){\displaystyle \frac{3}{2}}\mathrm{Li}_4(z)+{\displaystyle \frac{7}{2}}\mathrm{Li}_3(x)X{\displaystyle \frac{1}{2}}\mathrm{Li}_3(x)Y{\displaystyle \frac{9}{4}}\mathrm{Li}_3(x)+2\mathrm{Li}_3(y)X`$
$`{\displaystyle \frac{3}{4}}\mathrm{Li}_2(x)X^2+{\displaystyle \frac{9}{4}}\mathrm{Li}_2(x)X{\displaystyle \frac{11}{12}}\mathrm{Li}_2(x)\pi ^2\mathrm{Li}_2(y)XY{\displaystyle \frac{1}{16}}X^4+{\displaystyle \frac{1}{12}}X^3Y{\displaystyle \frac{1}{12}}X^3`$
$`{\displaystyle \frac{3}{8}}X^2Y^2+{\displaystyle \frac{1}{3}}X^2Y+{\displaystyle \frac{3}{8}}X^2\pi ^2+{\displaystyle \frac{1}{8}}X^2{\displaystyle \frac{5}{12}}XY\pi ^2{\displaystyle \frac{7}{2}}X\zeta _3+{\displaystyle \frac{11}{8}}X\pi ^2`$
$`+{\displaystyle \frac{1}{2}}Y\zeta _3+{\displaystyle \frac{9}{4}}\zeta _3+{\displaystyle \frac{71}{360}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{1}{2}}\mathrm{Li}_4(x){\displaystyle \frac{1}{2}}\mathrm{Li}_4(y){\displaystyle \frac{1}{2}}\mathrm{Li}_3(x)X+{\displaystyle \frac{1}{2}}\mathrm{Li}_3(y)Y+{\displaystyle \frac{1}{4}}\mathrm{Li}_2(x)X^2{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)\pi ^2`$
$`{\displaystyle \frac{1}{4}}\mathrm{Li}_2(y)Y^2{\displaystyle \frac{1}{48}}X^4+{\displaystyle \frac{1}{12}}X^3Y+{\displaystyle \frac{1}{12}}X^2\pi ^2+{\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{1}{12}}XY^3{\displaystyle \frac{1}{6}}XY\pi ^2`$
$`+XY{\displaystyle \frac{23}{6}}X+{\displaystyle \frac{1}{48}}Y^4{\displaystyle \frac{1}{12}}Y^2\pi ^2{\displaystyle \frac{3}{2}}Y^2+{\displaystyle \frac{23}{6}}Y+{\displaystyle \frac{1}{36}}\pi ^4{\displaystyle \frac{1}{2}}\pi ^2){\displaystyle \frac{t}{s}}`$
$`+{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{s}{t}}+\left({\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{1}{2}}XY{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{1}{4}}\pi ^2\right){\displaystyle \frac{t^2}{s^2}}`$
$`+(({\displaystyle \frac{3}{2}}\mathrm{Li}_3(x)+2\mathrm{Li}_3(y)+{\displaystyle \frac{15}{4}}\mathrm{Li}_2(x){\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{4}}X^3+{\displaystyle \frac{1}{4}}X^2Y{\displaystyle \frac{7}{6}}X^2`$
$`{\displaystyle \frac{11}{12}}XY+{\displaystyle \frac{1}{4}}X\pi ^2+{\displaystyle \frac{1}{2}}X{\displaystyle \frac{1}{6}}Y^3{\displaystyle \frac{1}{24}}Y^2+{\displaystyle \frac{1}{4}}Y\pi ^2{\displaystyle \frac{7}{4}}Y\zeta _3{\displaystyle \frac{1}{8}}\pi ^2+{\displaystyle \frac{47}{12}})`$
$`+({\displaystyle \frac{9}{2}}\mathrm{Li}_3(x)+4\mathrm{Li}_3(y){\displaystyle \frac{1}{4}}\mathrm{Li}_2(x)X+{\displaystyle \frac{9}{2}}\mathrm{Li}_2(x)2\mathrm{Li}_2(y)Y{\displaystyle \frac{5}{8}}X^3{\displaystyle \frac{25}{12}}X^2{\displaystyle \frac{1}{6}}XY`$
$`{\displaystyle \frac{1}{12}}X\pi ^2+{\displaystyle \frac{8}{3}}X{\displaystyle \frac{7}{12}}Y{\displaystyle \frac{9}{2}}\zeta _3{\displaystyle \frac{1}{2}}\pi ^2+{\displaystyle \frac{1}{2}}){\displaystyle \frac{t}{u}}`$
$`+(3\mathrm{Li}_3(x)+2\mathrm{Li}_3(y){\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X+{\displaystyle \frac{9}{4}}\mathrm{Li}_2(x)\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{4}}X^3{\displaystyle \frac{25}{24}}X^2+{\displaystyle \frac{2}{3}}XY`$
$`{\displaystyle \frac{1}{6}}X\pi ^2+{\displaystyle \frac{1}{4}}X3\zeta _3{\displaystyle \frac{3}{8}}\pi ^2){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{1}{2}}\mathrm{Li}_3(x)+{\displaystyle \frac{1}{2}}\mathrm{Li}_3(y)+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(x)X{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{12}}X^3+{\displaystyle \frac{1}{4}}X^2Y`$
$`{\displaystyle \frac{1}{4}}XY^2+2X+{\displaystyle \frac{1}{12}}Y^32Y){\displaystyle \frac{t}{s}}+{\displaystyle \frac{1}{2}}Y{\displaystyle \frac{s}{t}})i\pi `$
$`J_{s++}^{gg,(2),[3]}`$ $`=`$ $`({\displaystyle \frac{2}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{1}{3}}\mathrm{Li}_3(y){\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)X{\displaystyle \frac{1}{3}}\mathrm{Li}_2(y)Y{\displaystyle \frac{5}{24}}X^2Y+{\displaystyle \frac{1}{24}}XY^2`$
$`{\displaystyle \frac{29}{72}}XY+{\displaystyle \frac{5}{12}}X\pi ^2{\displaystyle \frac{1}{9}}Y^3{\displaystyle \frac{13}{24}}Y^2{\displaystyle \frac{17}{72}}Y\pi ^2+{\displaystyle \frac{19}{72}}Y{\displaystyle \frac{1}{3}}\zeta _3+{\displaystyle \frac{17}{12}}\pi ^2)`$
$`+(\mathrm{Li}_4(x)+\mathrm{Li}_4(y)+\mathrm{Li}_4(z)+\mathrm{Li}_3(x)Y+2\mathrm{Li}_3(y){\displaystyle \frac{1}{6}}\mathrm{Li}_2(x)\pi ^22\mathrm{Li}_2(y)Y`$
$`+{\displaystyle \frac{1}{24}}X^4{\displaystyle \frac{1}{6}}X^3Y+{\displaystyle \frac{1}{6}}X^3+{\displaystyle \frac{1}{4}}X^2Y^2{\displaystyle \frac{1}{24}}X^2Y+{\displaystyle \frac{1}{12}}X^2\pi ^2{\displaystyle \frac{5}{12}}X^2XY^2`$
$`{\displaystyle \frac{1}{6}}XY\pi ^2+{\displaystyle \frac{1}{12}}XY+{\displaystyle \frac{3}{4}}X\pi ^2+{\displaystyle \frac{1}{2}}XY\zeta _3+{\displaystyle \frac{1}{3}}Y\pi ^22\zeta _3+{\displaystyle \frac{7}{180}}\pi ^4{\displaystyle \frac{5}{12}}\pi ^2){\displaystyle \frac{t}{u}}`$
$`+(2\mathrm{Li}_4(x)+2\mathrm{Li}_4(y)+2\mathrm{Li}_4(z)+2\mathrm{Li}_3(x)Y{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)\pi ^2+{\displaystyle \frac{1}{12}}X^4{\displaystyle \frac{1}{3}}X^3Y`$
$`+{\displaystyle \frac{1}{2}}X^2Y^2+{\displaystyle \frac{1}{24}}X^2Y+{\displaystyle \frac{1}{6}}X^2\pi ^2+{\displaystyle \frac{1}{4}}X^2{\displaystyle \frac{1}{3}}XY\pi ^2{\displaystyle \frac{1}{12}}X\pi ^22Y\zeta _3+{\displaystyle \frac{7}{90}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{1}{3}}\mathrm{Li}_3(x){\displaystyle \frac{1}{3}}\mathrm{Li}_3(y){\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)X+{\displaystyle \frac{1}{3}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{9}}X^3{\displaystyle \frac{1}{4}}X^2Y{\displaystyle \frac{29}{24}}X^2`$
$`+{\displaystyle \frac{1}{4}}XY^2+XY+{\displaystyle \frac{47}{72}}X\pi ^2+{\displaystyle \frac{5}{18}}X{\displaystyle \frac{1}{9}}Y^3+{\displaystyle \frac{5}{24}}Y^2{\displaystyle \frac{47}{72}}Y\pi ^2{\displaystyle \frac{5}{18}}Y{\displaystyle \frac{1}{2}}\pi ^2){\displaystyle \frac{t}{s}}`$
$`{\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{s}{t}}`$
$`+\left({\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{1}{2}}XY{\displaystyle \frac{1}{4}}Y^2{\displaystyle \frac{1}{4}}\pi ^2\right){\displaystyle \frac{t^2}{s^2}}`$
$`+(({\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)+{\displaystyle \frac{1}{8}}X^2{\displaystyle \frac{2}{3}}XY{\displaystyle \frac{29}{72}}X+{\displaystyle \frac{5}{24}}Y^2{\displaystyle \frac{77}{24}}Y+{\displaystyle \frac{17}{72}}\pi ^2+{\displaystyle \frac{19}{72}})`$
$`+\left(\mathrm{Li}_3(x)+2\mathrm{Li}_2(x){\displaystyle \frac{1}{24}}X^2{\displaystyle \frac{1}{12}}XY{\displaystyle \frac{3}{4}}X+{\displaystyle \frac{1}{12}}Y\zeta _3+{\displaystyle \frac{1}{2}}\right){\displaystyle \frac{t}{u}}`$
$`+\left(2\mathrm{Li}_3(x)+{\displaystyle \frac{1}{24}}X^2+{\displaystyle \frac{1}{12}}XY+{\displaystyle \frac{1}{2}}X2\zeta _3\right){\displaystyle \frac{t^2}{u^2}}`$
$`+\left({\displaystyle \frac{2}{3}}\mathrm{Li}_2(x){\displaystyle \frac{1}{12}}X^2{\displaystyle \frac{1}{3}}XY+{\displaystyle \frac{11}{36}}X+{\displaystyle \frac{1}{12}}Y^2{\displaystyle \frac{11}{36}}Y+{\displaystyle \frac{1}{18}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`Y{\displaystyle \frac{s}{t}})i\pi `$
$`K_{s++}^{gg,(2),[3]}`$ $`=`$ $`({\displaystyle \frac{1}{3}}\mathrm{Li}_3(y)+{\displaystyle \frac{1}{3}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{24}}X^2Y+{\displaystyle \frac{1}{24}}XY^2+{\displaystyle \frac{1}{8}}XY{\displaystyle \frac{1}{12}}X\pi ^2{\displaystyle \frac{3}{2}}Y^2+{\displaystyle \frac{7}{36}}Y\pi ^2`$
$`+{\displaystyle \frac{29}{24}}Y+{\displaystyle \frac{1}{3}}\zeta _3{\displaystyle \frac{1}{8}}\pi ^2)`$
$`+(2\mathrm{Li}_4(x)+2\mathrm{Li}_4(y)+2\mathrm{Li}_4(z)+2\mathrm{Li}_3(x)Y+{\displaystyle \frac{1}{3}}\mathrm{Li}_3(x)+4\mathrm{Li}_3(y){\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)X`$
$`{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)\pi ^24\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{12}}X^4{\displaystyle \frac{1}{3}}X^3Y+{\displaystyle \frac{1}{2}}X^2Y^2+{\displaystyle \frac{1}{12}}X^2Y+{\displaystyle \frac{1}{6}}X^2\pi ^2`$
$`+{\displaystyle \frac{3}{4}}X^22XY^2{\displaystyle \frac{1}{3}}XY\pi ^2{\displaystyle \frac{1}{4}}XY{\displaystyle \frac{4}{9}}X\pi ^2+{\displaystyle \frac{5}{6}}X2Y\zeta _3+{\displaystyle \frac{2}{3}}Y\pi ^2`$
$`{\displaystyle \frac{13}{3}}\zeta _3+{\displaystyle \frac{7}{90}}\pi ^4{\displaystyle \frac{17}{36}}\pi ^2){\displaystyle \frac{t}{u}}`$
$`+(4\mathrm{Li}_4(x)+4\mathrm{Li}_4(y)+4\mathrm{Li}_4(z)+4\mathrm{Li}_3(x)Y{\displaystyle \frac{1}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)X{\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)\pi ^2`$
$`+{\displaystyle \frac{1}{6}}X^4{\displaystyle \frac{2}{3}}X^3Y+X^2Y^2+{\displaystyle \frac{1}{24}}X^2Y+{\displaystyle \frac{1}{3}}X^2\pi ^2+{\displaystyle \frac{1}{4}}X^2{\displaystyle \frac{2}{3}}XY\pi ^2+{\displaystyle \frac{7}{36}}X\pi ^2`$
$`4Y\zeta _3+{\displaystyle \frac{1}{3}}\zeta _3+{\displaystyle \frac{7}{45}}\pi ^4){\displaystyle \frac{t^2}{u^2}}`$
$`+({\displaystyle \frac{1}{3}}\mathrm{Li}_3(x){\displaystyle \frac{1}{3}}\mathrm{Li}_3(y){\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)X+{\displaystyle \frac{1}{3}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{1}{4}}X^2{\displaystyle \frac{1}{2}}XY{\displaystyle \frac{5}{18}}X\pi ^2{\displaystyle \frac{1}{3}}X`$
$`+{\displaystyle \frac{1}{4}}Y^2+{\displaystyle \frac{5}{18}}Y\pi ^2+{\displaystyle \frac{1}{3}}Y+{\displaystyle \frac{1}{4}}\pi ^2){\displaystyle \frac{t}{s}}`$
$`{\displaystyle \frac{3}{4}}Y^2{\displaystyle \frac{s}{t}}`$
$`+(({\displaystyle \frac{1}{3}}\mathrm{Li}_2(x)+{\displaystyle \frac{1}{24}}X^2{\displaystyle \frac{1}{6}}XY+{\displaystyle \frac{1}{8}}X{\displaystyle \frac{1}{8}}Y^2{\displaystyle \frac{23}{8}}Y+{\displaystyle \frac{29}{24}})`$
$`+\left(2\mathrm{Li}_3(x)+{\displaystyle \frac{11}{3}}\mathrm{Li}_2(x)+{\displaystyle \frac{1}{4}}X^2+{\displaystyle \frac{1}{6}}XY+{\displaystyle \frac{5}{4}}X{\displaystyle \frac{1}{4}}Y2\zeta _3+{\displaystyle \frac{1}{18}}\pi ^2+{\displaystyle \frac{5}{6}}\right){\displaystyle \frac{t}{u}}`$
$`+\left(4\mathrm{Li}_3(x)+{\displaystyle \frac{1}{3}}\mathrm{Li}_2(x){\displaystyle \frac{1}{8}}X^2+{\displaystyle \frac{1}{12}}XY+{\displaystyle \frac{1}{2}}X4\zeta _3{\displaystyle \frac{1}{18}}\pi ^2\right){\displaystyle \frac{t^2}{u^2}}`$
$`+\left({\displaystyle \frac{2}{3}}\mathrm{Li}_2(x)+{\displaystyle \frac{1}{6}}X^2{\displaystyle \frac{1}{3}}XY{\displaystyle \frac{1}{6}}Y^2+{\displaystyle \frac{1}{18}}\pi ^2\right){\displaystyle \frac{t}{s}}`$
$`{\displaystyle \frac{3}{2}}Y{\displaystyle \frac{s}{t}})i\pi `$
$`L_{s++}^{gg,(2),[3]}`$ $`=`$ $`\left({\displaystyle \frac{1}{72}}XY{\displaystyle \frac{1}{36}}Y^2+{\displaystyle \frac{5}{54}}Y{\displaystyle \frac{5}{72}}\pi ^2\right)+\left({\displaystyle \frac{1}{36}}X^2{\displaystyle \frac{5}{54}}X{\displaystyle \frac{1}{36}}Y^2+{\displaystyle \frac{5}{54}}Y\right){\displaystyle \frac{t}{s}}`$
$`+\left(\left({\displaystyle \frac{1}{72}}X+{\displaystyle \frac{1}{24}}Y+{\displaystyle \frac{5}{54}}\right)+\left({\displaystyle \frac{1}{18}}X+{\displaystyle \frac{1}{18}}Y\right){\displaystyle \frac{t}{s}}\right)i\pi `$
### A.3 $`q(p_2,+)+\overline{q}(p_1,)g(p_3,)+g(p_4,+)`$

The finite remainders are obtained by Bose symmetry from the $`q(p_2,+)+\overline{q}(p_1,)g(p_3,+)+g(p_4,)`$ process.

$`A_{s++}^{gg,(2),[1]}`$ $`=`$ $`A_{s++}^{gg,(2),[2]}(tu),`$
$`A_{s++}^{gg,(2),[2]}`$ $`=`$ $`A_{s++}^{gg,(2),[1]}(tu),`$
$`A_{s++}^{gg,(2),[3]}`$ $`=`$ $`A_{s++}^{gg,(2),[3]}(tu),`$ (160)
etc.

### A.4 $`q(p_2,+)+\overline{q}(p_1,)g(p_3,)+g(p_4,)`$

The finite remainders are obtained by Bose symmetry from the $`q(p_2,+)+\overline{q}(p_1,)g(p_3,+)+g(p_4,)`$ process.

$`A_{s+}^{gg,(2),[1]}`$ $`=`$ $`A_{s+++}^{gg,(2),[2]}(tu),`$
$`A_{s+}^{gg,(2),[2]}`$ $`=`$ $`A_{s+++}^{gg,(2),[1]}(tu),`$
$`A_{s+}^{gg,(2),[3]}`$ $`=`$ $`A_{s+++}^{gg,(2),[3]}(tu),`$ (161)
etc.

## Appendix B Finite one-loop contributions

In this appendix we list the coefficients $`A,\mathrm{},E`$ for the finite one-loop amplitudes defined in Section 4.3 for $`q\overline{q}gg`$. Expressions for quark-gluon scattering as well as the processes involving photons are available from the authors as a FORM output file.

### B.1 $`q(p_2,+)+\overline{q}(p_1,)g(p_3,+)+g(p_4,+)`$

$`A_{s+++}^{gg,(1),[1]}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{s}{t}}{\displaystyle \frac{1}{3}}`$
$`B_{s+++}^{gg,(1),[1]}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{s}{t}}`$
$`C_{s+++}^{gg,(1),[1]}`$ $`=`$ $`{\displaystyle \frac{1}{3}}`$
$`A_{s+++}^{gg,(1),[2]}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{t}{u}}{\displaystyle \frac{1}{6}}`$
$`B_{s+++}^{gg,(1),[2]}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{t}{u}}{\displaystyle \frac{1}{2}}`$
$`C_{s+++}^{gg,(1),[2]}`$ $`=`$ $`{\displaystyle \frac{1}{3}}`$
$`D_{s+++}^{gg,(1),[3]}=E_{s+++}^{gg,(1),[3]}`$ $`=`$ $`0`$ (162)
### B.2 $`q(p_2,+)+\overline{q}(p_1,)g(p_3,+)+g(p_4,)`$

$`A_{s++}^{gg,(1),[1]}`$ $`=`$ $`\left({\displaystyle \frac{1}{2}}X^2+{\displaystyle \frac{11}{6}}X+{\displaystyle \frac{7}{2}}\right)+\left({\displaystyle \frac{1}{2}}X+{\displaystyle \frac{1}{2}}\right){\displaystyle \frac{t}{u}}+X{\displaystyle \frac{t^2}{u^2}}+{\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{t^3}{u^3}}`$ (163)
$`+\left(\left(X{\displaystyle \frac{11}{6}}\right){\displaystyle \frac{1}{2}}{\displaystyle \frac{t}{u}}+{\displaystyle \frac{t^2}{u^2}}+X{\displaystyle \frac{t^3}{u^3}}\right)i\pi `$
$`B_{s++}^{gg,(1),[1]}`$ $`=`$ $`\left({\displaystyle \frac{1}{2}}X^2+{\displaystyle \frac{3}{2}}X{\displaystyle \frac{7}{2}}\right)+\left({\displaystyle \frac{3}{2}}X^2+{\displaystyle \frac{5}{2}}X+{\displaystyle \frac{1}{2}}\right){\displaystyle \frac{t}{u}}+\left({\displaystyle \frac{3}{2}}X^2+X\right){\displaystyle \frac{t^2}{u^2}}+{\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{t^3}{u^3}}`$ (164)
$`+\left(\left(X+{\displaystyle \frac{3}{2}}\right)+\left(3X+{\displaystyle \frac{5}{2}}\right){\displaystyle \frac{t}{u}}+\left(3X+1\right){\displaystyle \frac{t^2}{u^2}}+X{\displaystyle \frac{t^3}{u^3}}\right)i\pi `$
$`C_{s++}^{gg,(1),[1]}`$ $`=`$ $`{\displaystyle \frac{1}{3}}X+{\displaystyle \frac{1}{3}}i\pi `$ (165)
$`A_{s++}^{gg,(1),[2]}`$ $`=`$ $`\left(Y^2+{\displaystyle \frac{10}{3}}Y+3\right){\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{s}{t}}+\left(\left(2Y{\displaystyle \frac{1}{3}}\right)Y{\displaystyle \frac{s}{t}}\right)i\pi `$ (166)
$`B_{s++}^{gg,(1),[2]}`$ $`=`$ $`4{\displaystyle \frac{1}{2}}Y^2{\displaystyle \frac{s}{t}}Yi\pi {\displaystyle \frac{s}{t}}`$ (167)
$`C_{s++}^{gg,(1),[2]}`$ $`=`$ $`{\displaystyle \frac{1}{3}}Y+{\displaystyle \frac{1}{3}}i\pi `$ (168)
$`D_{s++}^{gg,(1),[3]}`$ $`=`$ $`\left(XY{\displaystyle \frac{1}{2}}Y^2+{\displaystyle \frac{5}{3}}Y{\displaystyle \frac{1}{2}}\pi ^2\right)+{\displaystyle \frac{3}{4}}X^2{\displaystyle \frac{t}{u}}`$ (169)
$`+\left({\displaystyle \frac{1}{2}}X^2{\displaystyle \frac{5}{3}}X{\displaystyle \frac{1}{2}}Y^2+{\displaystyle \frac{5}{3}}Y\right){\displaystyle \frac{t}{s}}`$
$`+\left(\left(X+{\displaystyle \frac{5}{3}}\right)+{\displaystyle \frac{3}{2}}X{\displaystyle \frac{t}{u}}+\left(XY\right){\displaystyle \frac{t}{s}}\right)i\pi `$
$`E_{s++}^{gg,(1),[3]}`$ $`=`$ $`{\displaystyle \frac{1}{6}}Y{\displaystyle \frac{1}{6}}i\pi +\left({\displaystyle \frac{1}{6}}X{\displaystyle \frac{1}{6}}Y\right){\displaystyle \frac{t}{s}}`$ (170)
### B.3 $`q(p_2,+)+\overline{q}(p_1,)g(p_3,)+g(p_4,+)`$

The finite remainders are obtained by Bose symmetry from the $`q(p_2,+)+\overline{q}(p_1,)g(p_3,+)+g(p_4,)`$ process.

$`A_{s++}^{gg,(1),[1]}`$ $`=`$ $`A_{s++}^{gg,(1),[2]}(tu),`$
$`A_{s++}^{gg,(1),[2]}`$ $`=`$ $`A_{s++}^{gg,(1),[1]}(tu),`$
$`A_{s++}^{gg,(1),[3]}`$ $`=`$ $`A_{s++}^{gg,(1),[3]}(tu),`$ (171)
etc.

### B.4 $`q(p_2,+)+\overline{q}(p_1,)g(p_3,)+g(p_4,)`$

The finite remainders are obtained by Bose symmetry from the $`q(p_2,+)+\overline{q}(p_1,)g(p_3,+)+g(p_4,)`$ process.

$`A_{s+}^{gg,(1),[1]}`$ $`=`$ $`A_{s+++}^{gg,(1),[2]}(tu),`$
$`A_{s+}^{gg,(1),[2]}`$ $`=`$ $`A_{s+++}^{gg,(1),[1]}(tu),`$
$`A_{s+}^{gg,(1),[3]}`$ $`=`$ $`A_{s+++}^{gg,(1),[3]}(tu),`$ (172)
etc.