format: text

text:
# Two-loop QCD corrections to massless quark-gluon scattering11footnote 1Work supported in part by the UK Particle Physics and Astronomy Research Council and by the EU Fourth Framework Programme ‘Training and Mobility of Researchers’, Network ‘Quantum Chromodynamics and the Deep Structure of Elementary Particles’, contract FMRX-CT98-0194 (DG 12 - MIHT). C.A. acknowledges the financial support of the Greek government and M.E.T. acknowledges financial support from CONACyT and the CVCP. We thank the British Council and German Academic Exchange Service for support under ARC project 1050.

## 1 Introduction

In hadronic collisions, jet cross sections are computed as a convolution of hard partonic cross sections with parton-distribution functions, followed by the fragmentation of the final-state partons into hadrons. The theoretical prediction can be improved by including higher order corrections which have the effect of reducing the unphysical renormalisation- and factorisation-scale dependences and by improving the matching of the parton-level theoretical jet algorithm with the hadron-level experimental jet algorithm. At present jet production is described at next-to-leading order (NLO) and several numerical programs are available  which have been extensively used to compare with data from the TEVATRON and CERN S$`p\overline{p}`$S.

Improving the theoretical prediction to next-to-next-to-leading order (NNLO) requires several ingredients. First, the parton-density functions are needed to NNLO accuracy which in turn requires knowledge of the three-loop splitting functions. At present, the even moments of the splitting functions are known for the flavour singlet and non-singlet structure functions $`F_2`$ and $`F_L`$ up to $`N=12`$ while the odd moments up to $`N=13`$ are known for $`F_3`$ . The numerically small $`N_F^2`$ non-singlet contribution is also known . These moments are sufficient to parameterise the splitting functions in $`x`$-space  and NNLO global analyses  are starting to appear. Second, the hard scattering matrix elements can be computed at NNLO. In the high energy limit, the quarks can be assumed to be massless and at this order, there are contributions from $`24`$ tree-level diagrams , from $`23`$ one-loop-level diagrams  and two-loop $`22`$ diagrams.

The evaluation of the two-loop diagrams has been a challenge for the past few years due to the presence of two-loop planar and crossed boxes. In the massless parton limit and in dimensional regularisation, analytic expressions for these basic scalar integrals have now been provided by Smirnov  and Tausk  as series in $`ϵ=(4D)/2`$, where $`D`$ is the space-time dimension, together with constructive procedures for reducing tensor integral to a basis set of known scalar (master) integrals . This makes the calculation of the two-loop amplitudes for $`22`$ QCD scattering processes possible. Bern, Dixon and Kosower  were the first to address such scattering processes and provided analytic expressions for the maximal-helicity-violating two-loop amplitude for $`gggg`$. Subsequently, Bern, Dixon and Ghinculov  completed the two-loop calculation of physical $`22`$ scattering amplitudes for the QED processes $`e^+e^{}\mu ^+\mu ^{}`$ and $`e^+e^{}e^{}e^+`$. More recently, we have provided expressions relevant for unlike- and like-quark scattering in the massless limit in Refs.  respectively. The corresponding matrix elements for the self-interference of one-loop $`22`$ quark scattering processes are given in .

In this paper, we address the $`𝒪\left(\alpha _s^4\right)`$ one- and two-loop corrections to the QCD process

$$q+\overline{q}g+g,$$
(1)
together with the time-reversed and crossed processes

$`q+g`$ $``$ $`q+g,`$ (2)
$`g+\overline{q}`$ $``$ $`g+\overline{q},`$ (3)
$`g+g`$ $``$ $`q+\overline{q}.`$ (4)
As is in Refs. , we use the $`\overline{\mathrm{MS}}`$ renormalisation scheme to remove the ultraviolet singularities and conventional dimensional regularisation, where all external particles are treated in $`D`$ dimensions. We provide expressions for both the interference of tree-level and two-loop graphs as well as the self-interference of one-loop amplitudes. In each case, we find that the infrared pole structure agrees with that obtained using Catani’s general factorisation formulae . The finite remainders are the main new results presented in this paper and we give explicit analytic expressions valid for each of the processes of Eqs. (1)–(4) in terms of logarithms and polylogarithms that are real in the physical domain. For simplicity, we decompose our results according to the powers of the number of colours $`N`$ and the number of light-quark flavours $`N_F`$ .

Our paper is organised as follows. We first establish our notation in Sec. 2. Analytic expressions for the interference of the two-loop and tree-level amplitudes are given in Sec. 3, while formulae describing the self-interference of one-loop graphs are given in Sec. 4. In Sec. 3.1 we adopt the notation used in Ref. , to isolate the infrared singularity structure of the two-loop amplitudes in the $`\overline{\mathrm{MS}}`$ scheme and we demonstrate that the anticipated singularity structure agrees with our explicit calculation. The finite $`𝒪\left(ϵ^0\right)`$ remainder of the two-loop graphs is one of the main results of our paper and expressions appropriate for the $`q\overline{q}gg`$, $`qgqg`$ ($`g\overline{q}g\overline{q}`$) and $`ggq\overline{q}`$ scattering processes are given in Sec. 3.2 in terms of logarithms and polylogarithms that have no imaginary parts. Expressions for the self-interference of one-loop graphs are given in Sec. 4 in terms of the one-loop bubble integral in $`D=42ϵ`$ and the one-loop box integral in $`D=62ϵ`$. Analytic formulae connecting these integrals in the various kinematic regions are given in Appendix A. As for the two-loop contributions, the singularity structure agrees with that obtained using Catani’s formalism. Finally we conclude with a brief summary of the results in Sec. 5.

## 2 Notation

For calculational purposes, the process we consider is

$$q(p_1)+\overline{q}(p_2)+g(p_3)+g(p_4)0,$$
(5)
where the particles are all incoming and carry light-like momenta, satisfying

$$p_1^\mu +p_2^\mu +p_3^\mu +p_4^\mu =0,p_i^2=0.$$

The associated Mandelstam variables are given by

$$s=(p_1+p_2)^2,t=(p_2+p_3)^2,u=(p_1+p_3)^2,s+t+u=0.$$
(6)
We work in conventional dimensional regularisation treating all external quark and gluon 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 _0S_ϵ=\alpha _s\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],$$
(7)
where

$$S_ϵ=(4\pi )^ϵe^{ϵ\gamma },\gamma =0.5772\mathrm{}=\mathrm{Euler}\mathrm{constant}$$
(8)
is the typical phase-space volume factor in $`D=42ϵ`$ dimensions. 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},$$
(9)
where $`N`$ is the number of colours and

$$C_F=\frac{N^21}{2N},C_A=N,T_R=\frac{1}{2}.$$
(10)
The renormalised four point amplitude in the $`\overline{\mathrm{MS}}`$ scheme is thus

$`|`$ $`=`$ $`4\pi \alpha _s\left[|^{(0)}+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)|^{(1)}+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)^2|^{(2)}+𝒪\left(\alpha _s^3\right)\right],`$
where the $`|^{(i)}`$ represents the colour-space vector describing the $`i`$-loop amplitude. The dependence on both renormalisation scale $`\mu `$ and renormalisation scheme is implicit.

We denote the squared amplitude summed over spins and colours by

$$|=|(q+\overline{q}g+g)|^2=𝒞(s,t,u).$$
(12)
which is symmetric under the exchange of $`t`$ and $`u`$.

The squared matrix elements for the crossed processes are obtained by exchanging the Mandelstam variables and introducing a minus sign for each quark change between initial and final states

$`{\displaystyle |(g+gq+\overline{q})|^2}`$ $`=`$ $`𝒞(s,t,u),`$ (13)
$`{\displaystyle |(q+gq+g)|^2}`$ $`=`$ $`𝒞(u,t,s),`$ (14)
$`{\displaystyle |(g+\overline{q}g+\overline{q})|^2}`$ $`=`$ $`𝒞(u,t,s).`$ (15)
The function $`𝒞`$ can be expanded perturbatively to yield

$`𝒞(s,t,u)`$ $`=`$ $`16\pi ^2\alpha _s^2\left[𝒞^4(s,t,u)+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)𝒞^6(s,t,u)+\left({\displaystyle \frac{\alpha _s}{2\pi }}\right)^2𝒞^8(s,t,u)+𝒪\left(\alpha _s^3\right)\right],`$
where

$`𝒞^4(s,t,u)`$ $`=`$ $`^{(0)}|^{(0)}=`$ (17)
$`=`$ $`2{\displaystyle \frac{N^21}{N}}(1ϵ)\left({\displaystyle \frac{N^21}{ut}}2{\displaystyle \frac{N^2}{s^2}}\right)\left(t^2+u^2ϵs^2\right),`$
$`𝒞^6(s,t,u)`$ $`=`$ $`\left(^{(0)}|^{(1)}+^{(1)}|^{(0)}\right),`$ (18)
$`𝒞^8(s,t,u)`$ $`=`$ $`\left(^{(1)}|^{(1)}+^{(0)}|^{(2)}+^{(2)}|^{(0)}\right).`$ (19)
Expressions for $`𝒞^6`$ are given in Ref.  using dimensional regularisation to isolate the infrared and ultraviolet singularities.

In the following sections, we present expressions for the infrared singular and finite contributions to $`𝒞^8`$ and the crossed processes. For convenience, we divide $`𝒞^8(s,t,u)`$ into two pieces

* the pure two-loop contributions

$$𝒞^{8(2\times 0)}(s,t,u)=^{(0)}|^{(2)}+^{(2)}|^{(0)},$$
(20)
described in Sec. 3 and
* the self-interference of the one-loop amplitude

$$𝒞^{8(1\times 1)}(s,t,u)=^{(1)}|^{(1)},$$
(21)
described in Sec. 4.

As in Refs. , we use QGRAF  to produce the two-loop Feynman diagrams to construct $`|^{(2)}`$. We then project by $`^{(0)}|`$ and perform the summation over colours and spins. It should be noted that when summing over the gluon polarisations, we ensure that the polarisations states are transversal (i.e. physical) by using an axial gauge

$$\underset{\mathrm{spins}}{}ϵ_i^\mu ϵ_i^\nu =g^{\mu \nu }+\frac{n_i^\mu p_i^\nu +n_i^\nu p_i^\mu }{n_ip_i}$$
(22)
where $`p_i`$ is the momentum of gluon $`i`$ and $`n_i`$ is an arbitrary light-like 4-vector. For simplicity, we choose $`n_3^\mu =p_4^\mu `$ and $`n_4^\mu =p_3^\mu `$. Finally, the trace over the Dirac matrices is carried out in $`D`$ dimensions using conventional dimensional regularisation. It is then straightforward to identify the scalar and tensor integrals present and replace them with combinations of the basis set of master integrals using the tensor reduction of two-loop integrals described in , based on integration-by-parts  and Lorentz invariance  identities. The final result is a combination of master integrals in $`D=42ϵ`$ for which the expansions around $`ϵ=0`$ are given in .

## 3 Two-loop contribution

We further decompose the two-loop contributions as a sum of two terms

$$𝒞^{8(2\times 0)}(s,t,u)=𝒫oles(s,t,u)+inite(s,t,u).$$
(23)
$`𝒫oles`$ contains infrared singularities that will be analytically canceled by the infrared singularities occurring in radiative processes of the same order (ultraviolet divergences are removed by renormalisation). $`inite`$ is the remainder which is finite as $`ϵ0`$.

### 3.1 Infrared Pole Structure

Following the procedure outlined in Ref. , we can write the infrared pole structure, renormalised in the $`\overline{\mathrm{MS}}`$ scheme as

$`𝒫oles(s,t,u)=Re\{`$ (24)
$`[{\displaystyle \frac{V^2}{4N}}(𝑨(ϵ,s,t,u)^2𝑩(ϵ,s,t,u)𝑫(ϵ,s,t,u)+2𝑨(2ϵ,s,t,u)4{\displaystyle \frac{\beta _0}{ϵ}}𝑨(ϵ,s,t,u))`$
$`+{\displaystyle \frac{V}{4}}(𝑨(ϵ,s,t,u)𝑩(ϵ,s,t,u)𝑩(ϵ,s,t,u)𝑪(ϵ,s,t,u)+2𝑩(2ϵ,s,t,u)`$
$`4{\displaystyle \frac{\beta _0}{ϵ}}𝑩(ϵ,s,t,u))+{\displaystyle \frac{V}{4N}}𝑩(ϵ,s,t,u)𝑫(ϵ,s,u,t)]𝒯_1(s,t,u)`$
$`+[{\displaystyle \frac{V}{4N}}(𝑨(ϵ,s,t,u)^2+𝑩(ϵ,s,t,u)𝑫(ϵ,s,t,u)+4{\displaystyle \frac{\beta _0}{ϵ}}𝑨(ϵ,s,t,u)2𝑨(2ϵ,s,t,u))`$
$`+{\displaystyle \frac{V}{4}}(𝑨(ϵ,s,t,u)𝑩(ϵ,s,t,u)𝑩(ϵ,s,t,u)𝑪(ϵ,s,t,u)4{\displaystyle \frac{\beta _0}{ϵ}}𝑩(ϵ,s,t,u)`$
$`+2𝑩(2ϵ,s,t,u)){\displaystyle \frac{V^2}{4N}}𝑩(ϵ,s,t,u)𝑫(ϵ,s,u,t)]𝒯_2(s,t,u)`$
$`+\left[{\displaystyle \frac{V^2}{4N}}𝑨(ϵ,s,t,u)+{\displaystyle \frac{V}{4}}𝑩(ϵ,s,t,u)\right]_1(s,t,u)`$
$`+\left[{\displaystyle \frac{V}{4N}}𝑨(ϵ,s,t,u)+{\displaystyle \frac{V}{4}}𝑩(ϵ,s,t,u)\right]_2(s,t,u)`$
$`+\left[{\displaystyle \frac{V^2}{4N}}𝑫(ϵ,s,t,u)+{\displaystyle \frac{V}{4}}𝑪(ϵ,s,t,u){\displaystyle \frac{V}{4N}}𝑫(ϵ,s,u,t)\right]_3(s,t,u)+_2`$
$`+(ut)\},`$
where

$`𝑨(ϵ,s,t,u)`$ $`=`$ $`\{({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{2ϵ}})[N({\displaystyle \frac{\mu ^2}{u}})^ϵ+C_F({\displaystyle \frac{\mu ^2}{s}})^ϵ]`$ (25)
$`+N({\displaystyle \frac{\beta _0}{2Nϵ}}{\displaystyle \frac{3}{4ϵ}})[({\displaystyle \frac{\mu ^2}{s}})^ϵ+({\displaystyle \frac{\mu ^2}{u}})^ϵ]\}{\displaystyle \frac{e^{ϵ\gamma }}{\mathrm{\Gamma }(1ϵ)}}`$
$`𝑩(ϵ,s,t,u)`$ $`=`$ $`\left({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{4ϵ}}+{\displaystyle \frac{\beta _0}{2Nϵ}}\right)\left[\left({\displaystyle \frac{\mu ^2}{t}}\right)^ϵ\left({\displaystyle \frac{\mu ^2}{s}}\right)^ϵ\right]{\displaystyle \frac{e^{ϵ\gamma }}{\mathrm{\Gamma }(1ϵ)}}`$ (26)
$`𝑪(ϵ,s,t,u)`$ $`=`$ $`\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}}\right)^ϵ{\displaystyle \frac{e^{ϵ\gamma }}{\mathrm{\Gamma }(1ϵ)}}`$ (27)
$`𝑫(ϵ,s,t,u)`$ $`=`$ $`\left({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{4ϵ}}+{\displaystyle \frac{\beta _0}{2Nϵ}}\right)\left[\left({\displaystyle \frac{\mu ^2}{t}}\right)^ϵ\left({\displaystyle \frac{\mu ^2}{u}}\right)^ϵ\right]{\displaystyle \frac{e^{ϵ\gamma }}{\mathrm{\Gamma }(1ϵ)}}`$ (28)
and

$`V`$ $`=`$ $`N^21,`$ (29)
$``$ $`=`$ $`e^{ϵ\gamma }{\displaystyle \frac{\mathrm{\Gamma }(12ϵ)}{\mathrm{\Gamma }(1ϵ)}}\left({\displaystyle \frac{\beta _0}{ϵ}}+K\right),`$ (30)
$`K`$ $`=`$ $`\left({\displaystyle \frac{67}{18}}{\displaystyle \frac{\pi ^2}{6}}\right)C_A{\displaystyle \frac{10}{9}}T_RN_F.`$ (31)
The tree-type structures, $`𝒯_1`$ and $`𝒯_2`$, are given by

$`𝒯_1(s,t,u)`$ $`=`$ $`{\displaystyle \frac{t}{u}}𝒯_2(s,t,u),`$ (32)
$`𝒯_2(s,t,u)`$ $`=`$ $`8(1ϵ){\displaystyle \frac{\left(t^2+u^2ϵs^2\right)}{s^2}},`$ (33)
while the interference of tree with one-loop structures are represented by $`_1`$, $`_2`$ and $`_3`$

$`_1(s,t,u)`$ $`=`$ $`{\displaystyle \frac{N^2+1}{N}}f_2(s,t,u){\displaystyle \frac{t}{Nu}}f_1(s,u,t)+f_3(s,t,u)`$ (34)
$`6\beta _0{\displaystyle \frac{(1ϵ)}{32ϵ}}𝒯_1(s,t,u)\mathrm{Bub}(s)`$
$`2ϵ(12ϵ)[{\displaystyle \frac{N}{ϵ^2}}\mathrm{Bub}(u)+{\displaystyle \frac{N}{2}}({\displaystyle \frac{1}{ϵ^2}}{\displaystyle \frac{2\beta _0}{Nϵ}}+{\displaystyle \frac{3}{2ϵ}})\mathrm{Bub}(s)`$
$`{\displaystyle \frac{1}{2N}}({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{3}{2ϵ}})\mathrm{Bub}(s)]𝒯_1(s,t,u)`$
$`_2(s,t,u)`$ $`=`$ $`{\displaystyle \frac{t}{u}}_1(s,u,t)`$ (35)
$`_3(s,t,u)`$ $`=`$ $`16(12ϵ){\displaystyle \frac{1}{u}}\left[t^2+u^2+\left(ut2t^22u^2\right)ϵ+\left(t^2+u^2+3ut\right)ϵ^2\right]\mathrm{Box}^6(t,u)`$ (36)
$`+2\left({\displaystyle \frac{1}{ϵ}}2\right)\left[\mathrm{Bub}(t)\mathrm{Bub}(s)\right]𝒯_1(s,t,u)+f_1(s,t,u)`$
$`+2\left({\displaystyle \frac{1}{ϵ}}2\right)\left[\mathrm{Bub}(u)\mathrm{Bub}(s)\right]𝒯_2(s,t,u)+{\displaystyle \frac{t}{u}}f_1(s,u,t),`$
where the infrared-finite functions $`f_1,f_2`$ and $`f_3`$ are

$`f_1(s,t,u)`$ $`=`$ $`8(12ϵ){\displaystyle \frac{1}{s^2}}[u(2u^2+5t^2+3tu)+ϵ(3t^34u^33tu^2)`$ (37)
$`ϵ^2s(4t^2+2u^2+5tu)+s^2tϵ^3]\mathrm{Box}^6(s,t)`$
$`f_2(s,t,u)`$ $`=`$ $`8(12ϵ){\displaystyle \frac{t}{us^2}}[(tu)(t^2+2u^2+tu)+ϵ(2t^3+ut^2+4tu^2+5u^3)`$ (38)
$`s^3ϵ^2]\mathrm{Box}^6(s,u)`$
$`f_3(s,t,u)`$ $`=`$ $`4(1ϵ){\displaystyle \frac{t}{su}}\{{\displaystyle \frac{V}{N}}[2stϵ(2s3u)3sϵ^2]`$
$`+4N(uuϵ+sϵ^2)\left\}\right[\mathrm{Bub}(u)\mathrm{Bub}(s)]`$
$`+{\displaystyle \frac{4t}{s^2u(1ϵ)(32ϵ)}}\{N[18u^2+15t^23t(st)`$
$`ϵ(78u^236t(st)+st)+ϵ^2(80u^2+10s(st)69st)]`$
$`{\displaystyle \frac{1}{N}}[24u^2+3tu21t^2+ϵ(85u^243t(st)+3st)`$
$`ϵ^2(112u^2+6s(st)109st)]+\beta _0[20s^240tu`$
$`2ϵ(38s^23us62tu)+4ϵ^2(27s^226tu)]\}ϵ\mathrm{Bub}(s)+𝒪(ϵ^3).`$
These expressions are valid in all kinematic regions. However, to evaluate the pole structure in a particular region, the one-loop bubble graph $`\mathrm{Bub}`$ and the one-loop box integral in $`D=62ϵ`$ dimensions, $`\mathrm{Box}^6`$, must be expanded as a series in $`ϵ`$. This analytic expansion is given in Appendix A.

The function $`_2`$, that appears in Eq. (24), exhibits only a single pole in $`ϵ`$ and is given by

$$_2^{(0)}|𝑯^{(2)}(ϵ)|^{(0)}=\frac{e^{ϵ\gamma }}{4ϵ\mathrm{\Gamma }(1ϵ)}H^{(2)}^{(0)}|^{(0)}$$

where the constant $`H^{(2)}`$ is

$`H^{(2)}`$ $`=`$ $`\left(\zeta _3+{\displaystyle \frac{5}{6}}+{\displaystyle \frac{11}{72}}\pi ^2\right)C_A^2+\left(13\zeta _3+{\displaystyle \frac{245}{108}}{\displaystyle \frac{23}{24}}\pi ^2\right)C_AC_F+{\displaystyle \frac{10}{27}}N_F^2`$
$`+\left({\displaystyle \frac{\pi ^2}{36}}{\displaystyle \frac{58}{27}}\right)C_AN_F+\left({\displaystyle \frac{\pi ^2}{12}}+{\displaystyle \frac{29}{54}}\right)N_FC_F+\left({\displaystyle \frac{3}{4}}+\pi ^212\zeta _3\right)C_F^2,`$
and $`\zeta _n`$ is the Riemann Zeta function with $`\zeta _2=\pi ^2/6`$ and $`\zeta _3=1.202056\mathrm{}`$ We note that $`H^{(2)}`$ is renormalisation-scheme dependent and Eq. (3.1) is valid in the $`\overline{\mathrm{MS}}`$ scheme. We also note that Eq. (3.1) differs from the corresponding expressions found in the singularity structure of two-loop quark-quark scattering in all but the $`C_F^2`$ coefficient. This is due to the presence of infrared emissions from gluons which modify the terms involving either $`C_A`$ or $`N_F`$.

It can be easily noted that the leading infrared singularity in Eq. (24) is $`𝒪\left(1/ϵ^4\right)`$. It is a very stringent check on the reliability of our calculation that the pole structure obtained by computing the Feynman diagrams directly and introducing series expansions in $`ϵ`$ for the scalar master integrals agrees with Eq. (24) through to $`𝒪\left(1/ϵ\right)`$. We therefore construct the finite remainder by subtracting Eq. (24) from the full result.

### 3.2 Finite contributions

The finite two-loop contribution to $`𝒞^8(s,t,u)`$ is defined as

$$inite(s,t,u)=2Re\left[^{(0)}|^{(2)\mathrm{fin}}\right].$$
(41)
In hadronic collisions, all parton scattering processes (Eqs. (1)–(4)) contribute simultaneously. We therefore need to evaluate $`inite(s,t,u)`$ for the $`q\overline{q}gg`$ and $`ggq\overline{q}`$ process (which we denote as the $`s`$-channel since, although the tree-level process contains graphs in all three channels, the squared tree matrix elements are proportional to $`1/s^2`$) and $`inite(u,t,s)`$ for the QCD Compton processes $`qgqg`$ and $`g\overline{q}g\overline{q}`$ (which we label as the $`u`$-channel).

Of course, the analytic expressions for the various processes are related by crossing symmetry. However, the master crossed boxes have cuts in all three channels yielding complex parts in all physical regions. The analytic continuation is therefore rather involved and prone to error. We therefore choose to give expressions describing $`𝒞^8(s,t,u)`$ and $`𝒞^8(u,t,s)`$ which are directly valid in the physical region $`s>0`$ and $`u,t<0`$, and are given in terms of logarithms and polylogarithms that have no imaginary parts.

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).`$ (42)
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}.$$
(43)
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),U=\mathrm{log}\left(\frac{u}{\mu ^2}\right),$$
(44)
where $`\mu `$ is the renormalisation scale.

For each channel, we choose to present our results by grouping terms according to the power of the number of colours $`N`$ and the number of light quarks $`N_F`$, so that in the generic $`c`$-channel we write

$`inite_c(s,t,u)=\left(N^21\right)`$
$`\times \left(N^3A_c+NB_c+{\displaystyle \frac{1}{N}}C_c+{\displaystyle \frac{1}{N^3}}D_c+N_FN^2E_c+N_FF_c+{\displaystyle \frac{N_F}{N^2}}G_c+N_F^2NH_c+{\displaystyle \frac{N_F^2}{N}}I_c\right).`$
#### 3.2.1 $`𝓕𝒊𝒏𝒊𝒕𝒆\mathbf{(}𝒔\mathbf{,}𝒕\mathbf{,}𝒖\mathbf{)}`$: the $`𝒔`$-channel process

We first give expressions valid for the annihilation processes, $`q\overline{q}gg`$ and $`gg\overline{q}q`$. We find that

$`A_s`$ $`=`$ $`\{({\displaystyle \frac{36077}{864}}{\displaystyle \frac{109}{144}}\pi ^4+{\displaystyle \frac{241}{12}}\zeta _3{\displaystyle \frac{5}{2}}\pi ^2+{\displaystyle \frac{2327}{216}}S+20\mathrm{Li}_4(x)+{\displaystyle \frac{11}{72}}\pi ^2S9\zeta _3S`$ (46)
$`+{\displaystyle \frac{121}{9}}S^221\mathrm{Li}_4(y)20\mathrm{Li}_4(z){\displaystyle \frac{35}{6}}\mathrm{Li}_3(y)+{\displaystyle \frac{1}{3}}XY^3+{\displaystyle \frac{35}{12}}Y^2X+{\displaystyle \frac{10}{3}}YX\pi ^2`$
$`+4\pi ^2\mathrm{Li}_2(x)+16Y\zeta _3+{\displaystyle \frac{1}{2}}Y^2\mathrm{Li}_2(y)+{\displaystyle \frac{35}{6}}Y\mathrm{Li}_2(y)+{\displaystyle \frac{121}{9}}YS5X^2Y^2`$
$`+{\displaystyle \frac{505}{72}}Y^25X^2{\displaystyle \frac{77}{36}}Y^3+{\displaystyle \frac{1}{8}}Y^4{\displaystyle \frac{5}{6}}X^4+{\displaystyle \frac{10}{3}}YX^3{\displaystyle \frac{2}{3}}\pi ^2Y^2{\displaystyle \frac{5}{3}}\pi ^2X^2`$
$`{\displaystyle \frac{143}{36}}Y\pi ^2{\displaystyle \frac{11}{3}}Y^2S20Y\mathrm{Li}_3(x)+{\displaystyle \frac{2273}{216}}Y){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{11}{90}}\pi ^4+{\displaystyle \frac{9}{2}}\zeta _3+{\displaystyle \frac{169}{36}}\pi ^2+{\displaystyle \frac{457}{144}}{\displaystyle \frac{11}{6}}S+{\displaystyle \frac{33}{4}}X\mathrm{Li}_2(x)8Y\zeta _3+{\displaystyle \frac{33}{8}}X^2Y`$
$`{\displaystyle \frac{33}{4}}\mathrm{Li}_3(x){\displaystyle \frac{2417}{144}}X+{\displaystyle \frac{2}{3}}Y\pi ^2+8Y\mathrm{Li}_3(x){\displaystyle \frac{17}{24}}X^3+{\displaystyle \frac{141}{16}}Y^2+{\displaystyle \frac{55}{6}}YS+X^2Y^2)`$
$`+({\displaystyle \frac{29}{72}}\pi ^4{\displaystyle \frac{205}{6}}\zeta _32\pi ^2{\displaystyle \frac{1535}{108}}S+{\displaystyle \frac{30377}{432}}22\mathrm{Li}_4(x){\displaystyle \frac{11}{36}}\pi ^2S+18\zeta _3S`$
$`{\displaystyle \frac{242}{9}}S^2+18\mathrm{Li}_4(y)+16\mathrm{Li}_4(z)+{\displaystyle \frac{5}{12}}\mathrm{Li}_3(y)+{\displaystyle \frac{57}{4}}\mathrm{Li}_3(x){\displaystyle \frac{2}{3}}XY^3{\displaystyle \frac{5}{24}}Y^2X`$
$`{\displaystyle \frac{8}{3}}YX\pi ^2+{\displaystyle \frac{977}{144}}X4\pi ^2\mathrm{Li}_2(x)Y^2\mathrm{Li}_2(y){\displaystyle \frac{5}{12}}Y\mathrm{Li}_2(y){\displaystyle \frac{55}{6}}XS`$
$`+{\displaystyle \frac{22}{3}}SX^28\mathrm{Li}_3(y)X{\displaystyle \frac{715}{18}}YS+2X^2Y^2{\displaystyle \frac{57}{4}}X\mathrm{Li}_2(x){\displaystyle \frac{57}{8}}X^2Y`$
$`+4X\mathrm{Li}_3(x)+4X\zeta _3X^2\mathrm{Li}_2(x)+{\displaystyle \frac{40}{9}}X\pi ^2{\displaystyle \frac{3847}{144}}Y^2{\displaystyle \frac{793}{144}}X^2+{\displaystyle \frac{401}{72}}Y^3`$
$`{\displaystyle \frac{1}{4}}Y^4+{\displaystyle \frac{1}{12}}X^4+{\displaystyle \frac{431}{72}}X^3{\displaystyle \frac{8}{3}}YX^3+{\displaystyle \frac{4}{3}}\pi ^2Y^2+4\pi ^2X^2+{\displaystyle \frac{47}{18}}Y\pi ^2+{\displaystyle \frac{22}{3}}Y^2S`$
$`+8Y\mathrm{Li}_3(x){\displaystyle \frac{3059}{432}}Y){\displaystyle \frac{t^2}{s^2}}\}+\{ut\}`$
$`B_s`$ $`=`$ $`\{({\displaystyle \frac{317}{180}}\pi ^4{\displaystyle \frac{733}{18}}\zeta _3{\displaystyle \frac{1055}{72}}\pi ^2{\displaystyle \frac{7543}{216}}S+2Y\mathrm{Li}_3(y)+8\mathrm{Li}_4(x)20YX\mathrm{Li}_2(y)`$ (47)
$`{\displaystyle \frac{44}{3}}YXS+{\displaystyle \frac{65}{9}}\pi ^2S+10\zeta _3S{\displaystyle \frac{121}{9}}S^279\mathrm{Li}_4(y)74\mathrm{Li}_4(z){\displaystyle \frac{71}{6}}\mathrm{Li}_3(y)`$
$`+{\displaystyle \frac{25}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{164771}{2592}}+{\displaystyle \frac{5}{3}}XY^3+{\displaystyle \frac{1}{4}}Y^2X+{\displaystyle \frac{19}{3}}YX\pi ^2{\displaystyle \frac{22}{27}}X{\displaystyle \frac{19}{3}}\pi ^2\mathrm{Li}_2(x)`$
$`+74Y\zeta _3+{\displaystyle \frac{1}{2}}Y^2\mathrm{Li}_2(y)+{\displaystyle \frac{71}{6}}Y\mathrm{Li}_2(y)+{\displaystyle \frac{16}{9}}XY{\displaystyle \frac{220}{9}}XS+24\mathrm{Li}_3(y)X`$
$`{\displaystyle \frac{22}{9}}YS{\displaystyle \frac{49}{2}}X^2Y^2{\displaystyle \frac{25}{3}}X\mathrm{Li}_2(x)+{\displaystyle \frac{55}{6}}X^2Y+44X\mathrm{Li}_3(x)36X\zeta _3`$
$`11X^2\mathrm{Li}_2(x){\displaystyle \frac{551}{36}}X\pi ^2{\displaystyle \frac{125}{72}}Y^2+{\displaystyle \frac{583}{36}}X^2+{\displaystyle \frac{127}{36}}Y^3{\displaystyle \frac{3}{8}}Y^4{\displaystyle \frac{25}{12}}X^4{\displaystyle \frac{17}{2}}X^3`$
$`+9YX^3{\displaystyle \frac{8}{3}}\pi ^2Y^2{\displaystyle \frac{16}{3}}\pi ^2X^2{\displaystyle \frac{83}{36}}Y\pi ^2+{\displaystyle \frac{22}{3}}Y^2S70Y\mathrm{Li}_3(x){\displaystyle \frac{5765}{216}}Y){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{229}{240}}\pi ^4+{\displaystyle \frac{118}{3}}\zeta _3{\displaystyle \frac{799}{72}}\pi ^2{\displaystyle \frac{11}{6}}S+{\displaystyle \frac{139}{36}}18Y\mathrm{Li}_3(y)+{\displaystyle \frac{65}{2}}\mathrm{Li}_4(y)`$
$`{\displaystyle \frac{299}{6}}\mathrm{Li}_3(x)+{\displaystyle \frac{2}{3}}XY^3{\displaystyle \frac{5}{6}}YX\pi ^2+{\displaystyle \frac{14}{9}}X+19Y\zeta _3+{\displaystyle \frac{7}{4}}Y^2\mathrm{Li}_2(y)+{\displaystyle \frac{73}{4}}XY`$
$`{\displaystyle \frac{55}{6}}YS{\displaystyle \frac{13}{8}}X^2Y^2+{\displaystyle \frac{299}{6}}X\mathrm{Li}_2(x)+{\displaystyle \frac{62}{3}}X^2Y{\displaystyle \frac{311}{36}}Y^2{\displaystyle \frac{17}{48}}Y^4+{\displaystyle \frac{73}{18}}X^3`$
$`{\displaystyle \frac{47}{12}}\pi ^2Y^2{\displaystyle \frac{55}{6}}Y\pi ^2{\displaystyle \frac{11}{6}}Y^2S9Y\mathrm{Li}_3(x))`$
$`+({\displaystyle \frac{17}{60}}\pi ^4+{\displaystyle \frac{443}{18}}\zeta _3+{\displaystyle \frac{16}{9}}\pi ^2+{\displaystyle \frac{1502}{27}}S+20Y\mathrm{Li}_3(y)+{\displaystyle \frac{7}{2}}\mathrm{Li}_4(x)+8YX\mathrm{Li}_2(y)`$
$`{\displaystyle \frac{1}{12}}\pi ^2S2\zeta _3S{\displaystyle \frac{79}{2}}\mathrm{Li}_4(y)22\mathrm{Li}_4(z)+{\displaystyle \frac{71}{6}}\mathrm{Li}_3(y){\displaystyle \frac{71}{6}}\mathrm{Li}_3(x){\displaystyle \frac{16}{3}}XY^3`$
$`+{\displaystyle \frac{32}{3}}Y^2X{\displaystyle \frac{8}{3}}YX\pi ^2+{\displaystyle \frac{359}{36}}X+7\pi ^2\mathrm{Li}_2(x)Y\zeta _3{\displaystyle \frac{45}{4}}Y^2\mathrm{Li}_2(y){\displaystyle \frac{71}{6}}Y\mathrm{Li}_2(y)`$
$`+{\displaystyle \frac{391}{18}}XY+{\displaystyle \frac{11}{6}}XS+{\displaystyle \frac{11}{6}}SX^2+7\mathrm{Li}_3(y)X{\displaystyle \frac{11}{6}}YS+{\displaystyle \frac{13}{4}}X^2Y^2+{\displaystyle \frac{71}{6}}X\mathrm{Li}_2(x)`$
$`24X^2Y+16X\mathrm{Li}_3(x)+X\zeta _3{\displaystyle \frac{27}{4}}X^2\mathrm{Li}_2(x)+{\displaystyle \frac{103}{18}}X\pi ^2{\displaystyle \frac{203}{36}}Y^2{\displaystyle \frac{277}{12}}X^2`$
$`{\displaystyle \frac{101}{18}}Y^3+{\displaystyle \frac{31}{48}}Y^4{\displaystyle \frac{33}{16}}X^4+{\displaystyle \frac{341}{18}}X^3+3YX^3+{\displaystyle \frac{11}{4}}\pi ^2Y^2+{\displaystyle \frac{11}{4}}\pi ^2X^2{\displaystyle \frac{103}{18}}Y\pi ^2`$
$`{\displaystyle \frac{11}{6}}Y^2S7Y\mathrm{Li}_3(x){\displaystyle \frac{19139}{324}}+{\displaystyle \frac{601}{36}}Y){\displaystyle \frac{t^2}{s^2}}`$
$`X^2{\displaystyle \frac{t^4}{u^2s^2}}\}+\{ut\}`$
$`C_s`$ $`=`$ $`\{({\displaystyle \frac{41393}{2592}}+{\displaystyle \frac{301}{720}}\pi ^4+{\displaystyle \frac{173}{36}}\zeta _3{\displaystyle \frac{133}{36}}\pi ^2+{\displaystyle \frac{5297}{216}}S+4Y\mathrm{Li}_3(y)20\mathrm{Li}_4(x)`$ (48)
$`4YX\mathrm{Li}_2(y){\displaystyle \frac{13}{24}}\pi ^2S+5\zeta _3S+\mathrm{Li}_4(y)+8\mathrm{Li}_4(z){\displaystyle \frac{15}{2}}\mathrm{Li}_3(y)3\mathrm{Li}_3(x)`$
$`{\displaystyle \frac{2}{3}}XY^3{\displaystyle \frac{11}{12}}Y^2X{\displaystyle \frac{14}{3}}YX\pi ^2+{\displaystyle \frac{59}{3}}X{\displaystyle \frac{8}{3}}\pi ^2\mathrm{Li}_2(x)12Y\zeta _3{\displaystyle \frac{1}{2}}Y^2\mathrm{Li}_2(y)`$
$`+{\displaystyle \frac{15}{2}}Y\mathrm{Li}_2(y)+2XY{\displaystyle \frac{22}{3}}SX^2+10\mathrm{Li}_3(y)X11YS+{\displaystyle \frac{7}{2}}X^2Y^2+3X\mathrm{Li}_2(x)`$
$`{\displaystyle \frac{10}{3}}X^2Y+10X\mathrm{Li}_3(x)6X\zeta _34X^2\mathrm{Li}_2(x)+{\displaystyle \frac{49}{18}}X\pi ^2{\displaystyle \frac{425}{72}}Y^2+{\displaystyle \frac{179}{18}}X^2`$
$`{\displaystyle \frac{23}{36}}Y^3+{\displaystyle \frac{3}{8}}Y^4+{\displaystyle \frac{1}{3}}X^4{\displaystyle \frac{65}{18}}X^3{\displaystyle \frac{8}{3}}YX^3+{\displaystyle \frac{5}{6}}\pi ^2Y^2+{\displaystyle \frac{3}{2}}\pi ^2X^2+{\displaystyle \frac{5}{18}}Y\pi ^2`$
$`{\displaystyle \frac{11}{3}}Y^2S+14Y\mathrm{Li}_3(x)+{\displaystyle \frac{271}{24}}Y){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{271}{240}}\pi ^4+{\displaystyle \frac{3}{2}}\zeta _3{\displaystyle \frac{127}{72}}\pi ^2+{\displaystyle \frac{13}{4}}X\mathrm{Li}_2(x)+26Y\mathrm{Li}_3(y){\displaystyle \frac{21}{4}}Y^2\mathrm{Li}_2(y)25Y\zeta _3`$
$`{\displaystyle \frac{119}{24}}X^2Y{\displaystyle \frac{13}{4}}\mathrm{Li}_3(x)+{\displaystyle \frac{3133}{144}}X{\displaystyle \frac{4}{3}}XY^3+{\displaystyle \frac{7}{4}}\pi ^2Y^2+{\displaystyle \frac{55}{18}}Y\pi ^2{\displaystyle \frac{55}{12}}XY`$
$`{\displaystyle \frac{83}{2}}\mathrm{Li}_4(y){\displaystyle \frac{7}{6}}YX\pi ^2{\displaystyle \frac{22}{3}}Y^2S+15Y\mathrm{Li}_3(x)+{\displaystyle \frac{23}{16}}{\displaystyle \frac{209}{72}}X^3+{\displaystyle \frac{919}{144}}Y^2`$
$`+{\displaystyle \frac{9}{16}}Y^4{\displaystyle \frac{22}{3}}YS+{\displaystyle \frac{23}{8}}X^2Y^2)`$
$`+({\displaystyle \frac{13}{360}}\pi ^4+15\zeta _3+{\displaystyle \frac{1}{3}}\pi ^2{\displaystyle \frac{3}{4}}S+{\displaystyle \frac{11}{2}}\mathrm{Li}_4(x)+\pi ^2S12\zeta _3S{\displaystyle \frac{11}{2}}\mathrm{Li}_4(y)6\mathrm{Li}_4(z)`$
$`+{\displaystyle \frac{11}{4}}\mathrm{Li}_3(y){\displaystyle \frac{11}{4}}\mathrm{Li}_3(x){\displaystyle \frac{255}{16}}{\displaystyle \frac{25}{8}}Y^2X+{\displaystyle \frac{2}{3}}YX\pi ^2{\displaystyle \frac{21}{16}}X+{\displaystyle \frac{1}{3}}\pi ^2\mathrm{Li}_2(x)`$
$`+3Y\zeta _3{\displaystyle \frac{1}{4}}Y^2\mathrm{Li}_2(y){\displaystyle \frac{11}{4}}Y\mathrm{Li}_2(y){\displaystyle \frac{11}{2}}XY+3\mathrm{Li}_3(y)X{\displaystyle \frac{3}{4}}X^2Y^2`$
$`+{\displaystyle \frac{11}{4}}X\mathrm{Li}_2(x)+{\displaystyle \frac{25}{8}}X^2Y3X\zeta _3+{\displaystyle \frac{1}{4}}X^2\mathrm{Li}_2(x){\displaystyle \frac{7}{6}}X\pi ^2+{\displaystyle \frac{83}{16}}Y^2{\displaystyle \frac{11}{16}}X^2`$
$`+{\displaystyle \frac{13}{24}}Y^3{\displaystyle \frac{5}{48}}Y^4{\displaystyle \frac{7}{48}}X^4{\displaystyle \frac{13}{24}}X^3+YX^3+{\displaystyle \frac{1}{12}}\pi ^2Y^2{\displaystyle \frac{7}{12}}\pi ^2X^2+{\displaystyle \frac{7}{6}}Y\pi ^2`$
$`3Y\mathrm{Li}_3(x)+{\displaystyle \frac{21}{16}}Y){\displaystyle \frac{t^2}{s^2}}+X^2{\displaystyle \frac{t^4}{u^2s^2}}\}+\{ut\}`$
$`D_s`$ $`=`$ $`\{({\displaystyle \frac{1}{90}}\pi ^4+{\displaystyle \frac{5}{2}}\zeta _3+{\displaystyle \frac{5}{8}}\pi ^2{\displaystyle \frac{3}{8}}S+2Y\mathrm{Li}_3(y)+{\displaystyle \frac{1}{2}}\pi ^2S6\zeta _3S+11\mathrm{Li}_4(y)`$ (49)
$`+14\mathrm{Li}_4(z){\displaystyle \frac{3}{2}}\mathrm{Li}_3(y)+8\mathrm{Li}_3(x){\displaystyle \frac{187}{32}}+{\displaystyle \frac{3}{4}}Y^2X{\displaystyle \frac{7}{3}}YX\pi ^2+X+\pi ^2\mathrm{Li}_2(x)`$
$`18Y\zeta _3{\displaystyle \frac{1}{2}}Y^2\mathrm{Li}_2(y)+{\displaystyle \frac{3}{2}}Y\mathrm{Li}_2(y)+2\mathrm{Li}_3(y)X+4X^2Y^28X\mathrm{Li}_2(x)`$
$`{\displaystyle \frac{11}{2}}X^2Y6X\mathrm{Li}_3(x)+4X\zeta _3X^2\mathrm{Li}_2(x)+{\displaystyle \frac{10}{3}}X\pi ^2+{\displaystyle \frac{5}{8}}Y^2{\displaystyle \frac{9}{4}}X^2{\displaystyle \frac{3}{4}}Y^3`$
$`{\displaystyle \frac{1}{8}}Y^4{\displaystyle \frac{1}{12}}X^4{\displaystyle \frac{7}{3}}YX^3+{\displaystyle \frac{1}{2}}\pi ^2Y^2+{\displaystyle \frac{1}{2}}\pi ^2X^2+2Y\pi ^2+16Y\mathrm{Li}_3(x)+{\displaystyle \frac{39}{8}}Y){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{1}{4}}{\displaystyle \frac{4}{45}}\pi ^42\zeta _3+{\displaystyle \frac{43}{12}}\pi ^219X\mathrm{Li}_2(x)4Y\mathrm{Li}_3(y)Y^2\mathrm{Li}_2(y)12Y\zeta _3`$
$`{\displaystyle \frac{15}{2}}X^2Y+19\mathrm{Li}_3(x)+{\displaystyle \frac{29}{4}}X{\displaystyle \frac{1}{3}}\pi ^2Y^2+{\displaystyle \frac{4}{3}}Y\pi ^2{\displaystyle \frac{7}{2}}XY+10\mathrm{Li}_4(y)`$
$`{\displaystyle \frac{4}{3}}YX\pi ^2+16Y\mathrm{Li}_3(x){\displaystyle \frac{11}{6}}X^3+4Y^2{\displaystyle \frac{7}{12}}Y^4+2X^2Y^2)`$
$`+\left(3X+{\displaystyle \frac{7}{2}}X^2+3Y+{\displaystyle \frac{3}{2}}\pi ^2+{\displaystyle \frac{1}{2}}Y^23XY\right){\displaystyle \frac{t^2}{s^2}}`$
$`X^2{\displaystyle \frac{t^4}{u^2s^2}}\}+\{ut\}`$
$`E_s`$ $`=`$ $`\{({\displaystyle \frac{14}{45}}\pi ^4{\displaystyle \frac{17}{6}}\zeta _3+{\displaystyle \frac{77}{54}}\pi ^2+{\displaystyle \frac{185}{54}}S8\mathrm{Li}_4(x){\displaystyle \frac{1}{36}}\pi ^2S{\displaystyle \frac{44}{9}}S^2+8\mathrm{Li}_4(y)`$ (50)
$`+8\mathrm{Li}_4(z)+{\displaystyle \frac{4}{3}}\mathrm{Li}_3(y){\displaystyle \frac{2}{3}}Y^2X{\displaystyle \frac{4}{3}}YX\pi ^2{\displaystyle \frac{4}{3}}\pi ^2\mathrm{Li}_2(x)8Y\zeta _3{\displaystyle \frac{4}{3}}Y\mathrm{Li}_2(y)`$
$`{\displaystyle \frac{44}{9}}YS+2X^2Y^2{\displaystyle \frac{43}{18}}Y^2+2X^2+{\displaystyle \frac{7}{18}}Y^3+{\displaystyle \frac{1}{3}}X^4{\displaystyle \frac{4}{3}}YX^3+{\displaystyle \frac{2}{3}}\pi ^2X^2+{\displaystyle \frac{13}{18}}Y\pi ^2`$
$`+{\displaystyle \frac{2}{3}}Y^2S+8Y\mathrm{Li}_3(x)+{\displaystyle \frac{10}{9}}Y+{\displaystyle \frac{1307}{216}}){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{37}{36}}{\displaystyle \frac{11}{360}}\pi ^44\zeta _3{\displaystyle \frac{49}{36}}\pi ^2+{\displaystyle \frac{1}{3}}S4X\mathrm{Li}_2(x)+2Y\zeta _32X^2Y+4\mathrm{Li}_3(x)`$
$`+{\displaystyle \frac{91}{18}}X+{\displaystyle \frac{2}{3}}Y\pi ^22Y\mathrm{Li}_3(x)2Y^2{\displaystyle \frac{5}{3}}YS{\displaystyle \frac{1}{4}}X^2Y^2)`$
$`+({\displaystyle \frac{17}{180}}\pi ^4+{\displaystyle \frac{17}{3}}\zeta _3{\displaystyle \frac{32}{27}}\pi ^2{\displaystyle \frac{221}{27}}S+4\mathrm{Li}_4(x)+{\displaystyle \frac{1}{18}}\pi ^2S+{\displaystyle \frac{88}{9}}S^24\mathrm{Li}_4(y)`$
$`4\mathrm{Li}_4(z)+{\displaystyle \frac{4}{3}}\mathrm{Li}_3(y)4\mathrm{Li}_3(x){\displaystyle \frac{2}{3}}Y^2X+{\displaystyle \frac{2}{3}}YX\pi ^2{\displaystyle \frac{19}{18}}X+{\displaystyle \frac{2}{3}}\pi ^2\mathrm{Li}_2(x)`$
$`+2Y\zeta _3{\displaystyle \frac{4}{3}}Y\mathrm{Li}_2(y){\displaystyle \frac{863}{108}}+{\displaystyle \frac{5}{3}}XS{\displaystyle \frac{4}{3}}SX^2+2\mathrm{Li}_3(y)X+{\displaystyle \frac{109}{9}}YS`$
$`{\displaystyle \frac{1}{2}}X^2Y^2+4X\mathrm{Li}_2(x)+2X^2Y2X\zeta _3{\displaystyle \frac{16}{9}}X\pi ^2+{\displaystyle \frac{70}{9}}Y^2+{\displaystyle \frac{1}{9}}X^2{\displaystyle \frac{7}{9}}Y^3`$
$`{\displaystyle \frac{1}{6}}X^4{\displaystyle \frac{7}{9}}X^3+{\displaystyle \frac{2}{3}}YX^3{\displaystyle \frac{1}{3}}\pi ^2X^2{\displaystyle \frac{7}{9}}Y\pi ^2{\displaystyle \frac{4}{3}}Y^2S2Y\mathrm{Li}_3(x){\displaystyle \frac{31}{6}}Y){\displaystyle \frac{t^2}{s^2}}\}`$
$`+\{ut\}`$
$`F_s`$ $`=`$ $`\{({\displaystyle \frac{28}{45}}\pi ^4{\displaystyle \frac{10}{9}}\zeta _3+{\displaystyle \frac{347}{108}}\pi ^2+{\displaystyle \frac{7}{9}}S+16\mathrm{Li}_4(x)+{\displaystyle \frac{8}{3}}YXS{\displaystyle \frac{11}{9}}\pi ^2S+{\displaystyle \frac{44}{9}}S^2`$ (51)
$`16\mathrm{Li}_4(y)16\mathrm{Li}_4(z)+{\displaystyle \frac{4}{3}}\mathrm{Li}_3(y)+4\mathrm{Li}_3(x)+{\displaystyle \frac{1}{2}}Y^2X+{\displaystyle \frac{8}{3}}YX\pi ^2{\displaystyle \frac{161}{18}}X`$
$`+{\displaystyle \frac{8}{3}}\pi ^2\mathrm{Li}_2(x)+16Y\zeta _3{\displaystyle \frac{4}{3}}Y\mathrm{Li}_2(y){\displaystyle \frac{1}{2}}XY{\displaystyle \frac{3661}{324}}+{\displaystyle \frac{62}{9}}XS+{\displaystyle \frac{26}{9}}YS`$
$`4X^2Y^24X\mathrm{Li}_2(x){\displaystyle \frac{2}{3}}X^2Y{\displaystyle \frac{37}{18}}X\pi ^2+{\displaystyle \frac{35}{18}}Y^2{\displaystyle \frac{79}{18}}X^2{\displaystyle \frac{7}{9}}Y^3{\displaystyle \frac{2}{3}}X^4`$
$`+{\displaystyle \frac{2}{3}}X^3+{\displaystyle \frac{8}{3}}YX^3{\displaystyle \frac{4}{3}}\pi ^2X^2+{\displaystyle \frac{1}{18}}Y\pi ^2{\displaystyle \frac{4}{3}}Y^2S16Y\mathrm{Li}_3(x)+{\displaystyle \frac{23}{9}}Y){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{11}{90}}\pi ^4+{\displaystyle \frac{52}{3}}\zeta _3+{\displaystyle \frac{35}{12}}\pi ^2+{\displaystyle \frac{1}{3}}S+{\displaystyle \frac{52}{3}}X\mathrm{Li}_2(x)+8Y\zeta _3+{\displaystyle \frac{47}{6}}X^2Y{\displaystyle \frac{52}{3}}\mathrm{Li}_3(x)`$
$`{\displaystyle \frac{50}{9}}X+{\displaystyle \frac{1}{18}}Y\pi ^2{\displaystyle \frac{1}{3}}XY+{\displaystyle \frac{1}{3}}Y^2S8Y\mathrm{Li}_3(x){\displaystyle \frac{11}{36}}X^3+{\displaystyle \frac{49}{18}}Y^2{\displaystyle \frac{19}{36}}`$
$`+{\displaystyle \frac{5}{3}}YSX^2Y^2)`$
$`+({\displaystyle \frac{263}{27}}S{\displaystyle \frac{1}{3}}XS+{\displaystyle \frac{3}{2}}X\pi ^2{\displaystyle \frac{1}{3}}SX^2+{\displaystyle \frac{5}{2}}X^2Y+{\displaystyle \frac{1}{3}}YS{\displaystyle \frac{1}{6}}\pi ^2S`$
$`+{\displaystyle \frac{1}{3}}Y^2S+{\displaystyle \frac{7}{9}}X+{\displaystyle \frac{4}{3}}X\mathrm{Li}_2(x)+{\displaystyle \frac{211}{18}}X^2{\displaystyle \frac{4}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{4085}{324}}+{\displaystyle \frac{4}{3}}\mathrm{Li}_3(y){\displaystyle \frac{1}{9}}\zeta _3`$
$`{\displaystyle \frac{31}{9}}Y{\displaystyle \frac{4}{3}}Y\mathrm{Li}_2(y){\displaystyle \frac{3}{2}}Y\pi ^2{\displaystyle \frac{73}{36}}X^3{\displaystyle \frac{7}{6}}Y^2X+{\displaystyle \frac{25}{36}}Y^3+{\displaystyle \frac{25}{18}}\pi ^2{\displaystyle \frac{41}{18}}Y^2`$
$`{\displaystyle \frac{58}{9}}XY){\displaystyle \frac{t^2}{s^2}}3X^2{\displaystyle \frac{t^4}{u^2s^2}}\}+\{ut\}`$
$`G_s`$ $`=`$ $`\{({\displaystyle \frac{28}{15}}\pi ^4+{\displaystyle \frac{143}{18}}\zeta _3+{\displaystyle \frac{19}{36}}\pi ^2{\displaystyle \frac{227}{54}}S48\mathrm{Li}_4(x){\displaystyle \frac{1}{12}}\pi ^2S+48\mathrm{Li}_4(y)+48\mathrm{Li}_4(z)`$ (52)
$`8\mathrm{Li}_3(x)+{\displaystyle \frac{3401}{648}}+{\displaystyle \frac{1}{6}}Y^2X8YX\pi ^2+{\displaystyle \frac{101}{6}}X8\pi ^2\mathrm{Li}_2(x)48Y\zeta _3`$
$`+{\displaystyle \frac{1}{2}}XY+{\displaystyle \frac{4}{3}}SX^2+2YS+12X^2Y^2+8X\mathrm{Li}_2(x)+{\displaystyle \frac{1}{3}}X^2Y+{\displaystyle \frac{64}{9}}X\pi ^2`$
$`+{\displaystyle \frac{4}{9}}Y^2+{\displaystyle \frac{89}{9}}X^2+{\displaystyle \frac{7}{18}}Y^3+2X^4+{\displaystyle \frac{7}{9}}X^38YX^3+4\pi ^2X^2+{\displaystyle \frac{2}{9}}Y\pi ^2+{\displaystyle \frac{2}{3}}Y^2S`$
$`+48Y\mathrm{Li}_3(x){\displaystyle \frac{11}{3}}Y){\displaystyle \frac{t}{u}}`$
$`+({\displaystyle \frac{11}{30}}\pi ^452\zeta _3{\displaystyle \frac{107}{18}}\pi ^252X\mathrm{Li}_2(x)24Y\zeta _3{\displaystyle \frac{71}{3}}X^2Y+52\mathrm{Li}_3(x)+{\displaystyle \frac{71}{9}}X`$
$`+{\displaystyle \frac{46}{9}}Y\pi ^2{\displaystyle \frac{2}{3}}XY+{\displaystyle \frac{4}{3}}Y^2S+24Y\mathrm{Li}_3(x)+{\displaystyle \frac{7}{9}}X^3{\displaystyle \frac{37}{9}}Y^2+{\displaystyle \frac{4}{3}}YS+3X^2Y^2)`$
$`+({\displaystyle \frac{10}{3}}X\pi ^28X4X\mathrm{Li}_2(x)18X^2+4\mathrm{Li}_3(x)4\mathrm{Li}_3(y)+8Y+4Y\mathrm{Li}_2(y)`$
$`+{\displaystyle \frac{10}{3}}Y\pi ^2+9Y^2){\displaystyle \frac{t^2}{s^2}}+9X^2{\displaystyle \frac{t^4}{u^2s^2}}\}+\{ut\}`$
$`H_s`$ $`=`$ $`\{({\displaystyle \frac{4}{27}}\pi ^2+{\displaystyle \frac{4}{9}}S^2{\displaystyle \frac{20}{27}}S{\displaystyle \frac{10}{27}}Y+{\displaystyle \frac{4}{9}}YS+{\displaystyle \frac{1}{6}}Y^2){\displaystyle \frac{t}{u}}`$ (53)
$`+({\displaystyle \frac{8}{27}}\pi ^2+{\displaystyle \frac{40}{27}}S{\displaystyle \frac{8}{9}}S^2+{\displaystyle \frac{20}{27}}Y{\displaystyle \frac{8}{9}}YS{\displaystyle \frac{1}{3}}Y^2){\displaystyle \frac{t^2}{s^2}}\}+\{ut\}`$
$`I_s`$ $`=`$ $`\{({\displaystyle \frac{1}{9}}XY+{\displaystyle \frac{10}{27}}Y{\displaystyle \frac{1}{6}}Y^2+{\displaystyle \frac{20}{27}}S{\displaystyle \frac{4}{9}}XS{\displaystyle \frac{4}{9}}S^2{\displaystyle \frac{1}{18}}X^2{\displaystyle \frac{7}{54}}\pi ^2`$ (54)
$`{\displaystyle \frac{4}{9}}YS+{\displaystyle \frac{10}{27}}X){\displaystyle \frac{t}{u}}{\displaystyle \frac{2}{9}}X(XY){\displaystyle \frac{t^2}{s^2}}\}+\{ut\}`$
#### 3.2.2 $`𝓕𝒊𝒏𝒊𝒕𝒆\mathbf{(}𝒖\mathbf{,}𝒕\mathbf{,}𝒔\mathbf{)}`$: the $`𝒖`$-channel process

Similarly, for the Compton scattering processes, $`qgqg`$ and $`g\overline{q}g\overline{q}`$, we find that the coefficients for the finite part in Eq. (3.2) are given by

$`A_u`$ $`=`$ $`({\displaystyle \frac{2417}{144}}X+{\displaystyle \frac{11}{45}}\pi ^4{\displaystyle \frac{19}{144}}\pi ^2+{\displaystyle \frac{3}{4}}\zeta _3{\displaystyle \frac{11}{3}}U+{\displaystyle \frac{33}{4}}\mathrm{Li}_3(x)+{\displaystyle \frac{457}{72}}+8\mathrm{Li}_3(x)Y`$ (55)
$`{\displaystyle \frac{55}{3}}UY{\displaystyle \frac{10}{3}}\pi ^2Y+8\zeta _3Y{\displaystyle \frac{2}{3}}\pi ^2Y^22X^2Y{\displaystyle \frac{33}{4}}X\mathrm{Li}_2(x)8X\zeta _3`$
$`{\displaystyle \frac{35}{24}}X\pi ^2+{\displaystyle \frac{4}{3}}XY^3+8\mathrm{Li}_3(y)X+{\displaystyle \frac{55}{6}}UX+{\displaystyle \frac{141}{16}}X^2{\displaystyle \frac{17}{24}}X^3{\displaystyle \frac{141}{8}}YX`$
$`+2Y^2X+2X^2Y^2{\displaystyle \frac{8}{3}}\pi ^2XY{\displaystyle \frac{2}{3}}Y^4{\displaystyle \frac{4}{3}}Y^3+{\displaystyle \frac{141}{8}}Y^2+{\displaystyle \frac{2417}{72}}Y)`$
$`+({\displaystyle \frac{121}{9}}U^2{\displaystyle \frac{36077}{864}}9U\zeta _3{\displaystyle \frac{121}{72}}U\pi ^2+{\displaystyle \frac{2273}{432}}X{\displaystyle \frac{181}{720}}\pi ^4{\displaystyle \frac{215}{144}}\pi ^2+{\displaystyle \frac{103}{6}}\zeta _3`$
$`+{\displaystyle \frac{2327}{216}}U+{\displaystyle \frac{1}{2}}\mathrm{Li}_4(z)+{\displaystyle \frac{1}{2}}\mathrm{Li}_4(y)+{\displaystyle \frac{35}{12}}\mathrm{Li}_3(x)10\mathrm{Li}_3(x)Y{\displaystyle \frac{121}{9}}UY{\displaystyle \frac{29}{18}}\pi ^2Y`$
$`6\zeta _3Y+{\displaystyle \frac{7}{12}}\pi ^2Y^2+{\displaystyle \frac{7}{4}}X^2Y{\displaystyle \frac{35}{12}}X\mathrm{Li}_2(x)+8X\zeta _3{\displaystyle \frac{1}{4}}X^2\mathrm{Li}_2(x){\displaystyle \frac{55}{36}}X\pi ^2`$
$`{\displaystyle \frac{11}{6}}UX^2{\displaystyle \frac{13}{6}}XY^3+{\displaystyle \frac{1}{12}}\pi ^2X^2{\displaystyle \frac{1}{4}}\pi ^2\mathrm{Li}_2(x){\displaystyle \frac{11}{3}}UY^210\mathrm{Li}_3(y)X`$
$`+{\displaystyle \frac{121}{18}}UX+{\displaystyle \frac{1}{12}}X^4+{\displaystyle \frac{145}{144}}X^2{\displaystyle \frac{77}{72}}X^3{\displaystyle \frac{1}{2}}X^3Y{\displaystyle \frac{145}{72}}YX{\displaystyle \frac{7}{4}}Y^2X{\displaystyle \frac{17}{8}}X^2Y^2`$
$`+{\displaystyle \frac{37}{12}}\pi ^2XY+{\displaystyle \frac{11}{3}}UXY{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)XY+{\displaystyle \frac{13}{12}}Y^4+{\displaystyle \frac{7}{6}}Y^3+{\displaystyle \frac{145}{72}}Y^2`$
$`{\displaystyle \frac{2273}{216}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]`$
$`+({\displaystyle \frac{11}{6}}U\pi ^2{\displaystyle \frac{2273}{432}}X{\displaystyle \frac{35}{36}}\pi ^4{\displaystyle \frac{865}{144}}\pi ^2+{\displaystyle \frac{35}{12}}\zeta _3{\displaystyle \frac{41}{2}}\mathrm{Li}_4(z){\displaystyle \frac{35}{6}}\mathrm{Li}_3(y)+20\mathrm{Li}_4(x)`$
$`+{\displaystyle \frac{41}{2}}\mathrm{Li}_4(y){\displaystyle \frac{35}{12}}\mathrm{Li}_3(x)10\mathrm{Li}_3(x)Y{\displaystyle \frac{7}{4}}\pi ^2Y20\mathrm{Li}_3(y)Y+10\zeta _3Y{\displaystyle \frac{13}{4}}\pi ^2Y^2`$
$`{\displaystyle \frac{7}{4}}X^2Y+{\displaystyle \frac{35}{12}}X\mathrm{Li}_2(x)8X\zeta _3+{\displaystyle \frac{1}{4}}X^2\mathrm{Li}_2(x)+{\displaystyle \frac{55}{36}}X\pi ^2+{\displaystyle \frac{11}{6}}UX^2{\displaystyle \frac{3}{2}}XY^3`$
$`{\displaystyle \frac{7}{4}}\pi ^2X^2+{\displaystyle \frac{35}{6}}\mathrm{Li}_2(y)Y{\displaystyle \frac{15}{4}}\pi ^2\mathrm{Li}_2(x)+10\mathrm{Li}_3(y)X{\displaystyle \frac{121}{18}}UX{\displaystyle \frac{11}{12}}X^4{\displaystyle \frac{865}{144}}X^2`$
$`+{\displaystyle \frac{77}{72}}X^3{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y^2+{\displaystyle \frac{23}{6}}X^3Y+{\displaystyle \frac{865}{72}}YX+{\displaystyle \frac{14}{3}}Y^2X{\displaystyle \frac{23}{8}}X^2Y^2+{\displaystyle \frac{1}{4}}\pi ^2XY`$
$`{\displaystyle \frac{11}{3}}UXY+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)XY\left)\right[{\displaystyle \frac{t^2s^2}{st}}]`$
$`+({\displaystyle \frac{242}{9}}U^2+18U\zeta _3+{\displaystyle \frac{253}{36}}U\pi ^2{\displaystyle \frac{4}{27}}X{\displaystyle \frac{7}{40}}\pi ^4+{\displaystyle \frac{35}{9}}\pi ^2{\displaystyle \frac{161}{6}}\zeta _3{\displaystyle \frac{1535}{108}}U`$
$`+2\mathrm{Li}_4(z)+2\mathrm{Li}_4(y){\displaystyle \frac{22}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{30377}{432}}+2\mathrm{Li}_3(x)Y+{\displaystyle \frac{440}{9}}UY+{\displaystyle \frac{155}{18}}\pi ^2Y`$
$`6\zeta _3Y+{\displaystyle \frac{5}{3}}\pi ^2Y^22X\mathrm{Li}_3(x){\displaystyle \frac{41}{3}}X^2Y+{\displaystyle \frac{22}{3}}X\mathrm{Li}_2(x)+4X\zeta _3+X^2\mathrm{Li}_2(x)`$
$`+{\displaystyle \frac{223}{36}}X\pi ^2+{\displaystyle \frac{22}{3}}UX^2+2XY^3{\displaystyle \frac{1}{3}}\pi ^2X^2+\pi ^2\mathrm{Li}_2(x)+{\displaystyle \frac{44}{3}}UY^22\mathrm{Li}_3(y)X`$
$`{\displaystyle \frac{220}{9}}UX{\displaystyle \frac{1}{3}}X^4{\displaystyle \frac{145}{9}}X^2+{\displaystyle \frac{52}{9}}X^3+{\displaystyle \frac{5}{3}}X^3Y+{\displaystyle \frac{290}{9}}YX+{\displaystyle \frac{41}{3}}Y^2X{\displaystyle \frac{3}{2}}X^2Y^2`$
$`+\pi ^2XY{\displaystyle \frac{44}{3}}UXY+2\mathrm{Li}_2(y)XYY^4{\displaystyle \frac{82}{9}}Y^3{\displaystyle \frac{290}{9}}Y^2+{\displaystyle \frac{8}{27}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{u^2}}]`$
$`+({\displaystyle \frac{2995}{432}}X+{\displaystyle \frac{38}{45}}\pi ^4+{\displaystyle \frac{509}{48}}\pi ^2+{\displaystyle \frac{83}{12}}\zeta _3+20\mathrm{Li}_4(z){\displaystyle \frac{83}{6}}\mathrm{Li}_3(y)16\mathrm{Li}_4(x)20\mathrm{Li}_4(y)`$
$`{\displaystyle \frac{83}{12}}\mathrm{Li}_3(x)+10\mathrm{Li}_3(x)Y+{\displaystyle \frac{17}{6}}\pi ^2Y+20\mathrm{Li}_3(y)Y10\zeta _3Y+2\pi ^2Y^22X\mathrm{Li}_3(x)`$
$`+{\displaystyle \frac{17}{6}}X^2Y+{\displaystyle \frac{83}{12}}X\mathrm{Li}_2(x)+4X\zeta _3+{\displaystyle \frac{37}{24}}X\pi ^2+2XY^3+{\displaystyle \frac{4}{3}}\pi ^2X^2+{\displaystyle \frac{83}{6}}\mathrm{Li}_2(y)Y`$
$`+4\pi ^2\mathrm{Li}_2(x)10\mathrm{Li}_3(y)X+{\displaystyle \frac{275}{18}}UX+{\displaystyle \frac{2}{3}}X^4+{\displaystyle \frac{509}{48}}X^2+{\displaystyle \frac{5}{24}}X^33X^3Y`$
$`{\displaystyle \frac{509}{24}}YX+{\displaystyle \frac{49}{12}}Y^2X+2X^2Y^2\left)\right[{\displaystyle \frac{t^2s^2}{u^2}}]`$
$`B_u`$ $`=`$ $`({\displaystyle \frac{11}{6}}U\pi ^2+{\displaystyle \frac{14}{9}}X+{\displaystyle \frac{119}{180}}\pi ^4+24\pi ^2+{\displaystyle \frac{173}{6}}\zeta _3{\displaystyle \frac{11}{3}}U{\displaystyle \frac{65}{2}}\mathrm{Li}_4(z){\displaystyle \frac{65}{2}}\mathrm{Li}_4(y)`$ (56)
$`+{\displaystyle \frac{299}{6}}\mathrm{Li}_3(x)27\mathrm{Li}_3(x)Y+{\displaystyle \frac{55}{3}}UY{\displaystyle \frac{29}{6}}\pi ^2Y11\zeta _3Y+{\displaystyle \frac{1}{4}}\pi ^2Y^2`$
$`+18X\mathrm{Li}_3(x){\displaystyle \frac{197}{6}}X^2Y{\displaystyle \frac{299}{6}}X\mathrm{Li}_2(x)+X\zeta _3{\displaystyle \frac{7}{4}}X^2\mathrm{Li}_2(x){\displaystyle \frac{11}{6}}X\pi ^2`$
$`{\displaystyle \frac{11}{6}}UX^2{\displaystyle \frac{23}{6}}XY^3{\displaystyle \frac{5}{4}}\pi ^2X^2{\displaystyle \frac{7}{4}}\pi ^2\mathrm{Li}_2(x){\displaystyle \frac{11}{3}}UY^2+9\mathrm{Li}_3(y)X{\displaystyle \frac{55}{6}}UX`$
$`{\displaystyle \frac{41}{24}}X^4{\displaystyle \frac{311}{36}}X^2+{\displaystyle \frac{73}{18}}X^3+{\displaystyle \frac{37}{6}}X^3Y{\displaystyle \frac{173}{9}}YX+{\displaystyle \frac{73}{3}}Y^2X{\displaystyle \frac{55}{8}}X^2Y^2`$
$`+{\displaystyle \frac{27}{4}}\pi ^2XY+{\displaystyle \frac{11}{3}}UXY{\displaystyle \frac{7}{2}}\mathrm{Li}_2(y)XY+{\displaystyle \frac{139}{18}}+{\displaystyle \frac{23}{12}}Y^4{\displaystyle \frac{146}{9}}Y^3+{\displaystyle \frac{173}{9}}Y^2`$
$`{\displaystyle \frac{28}{9}}Y)`$
$`+({\displaystyle \frac{164771}{2592}}{\displaystyle \frac{121}{9}}U^2+10U\zeta _3{\displaystyle \frac{34}{9}}U\pi ^2{\displaystyle \frac{5941}{432}}X{\displaystyle \frac{4}{3}}\pi ^4+{\displaystyle \frac{1123}{144}}\pi ^2{\displaystyle \frac{1529}{36}}\zeta _3`$
$`{\displaystyle \frac{7543}{216}}U+{\displaystyle \frac{71}{2}}\mathrm{Li}_4(z)+{\displaystyle \frac{71}{2}}\mathrm{Li}_4(y)+{\displaystyle \frac{7}{4}}\mathrm{Li}_3(x)+{\displaystyle \frac{242}{9}}UY{\displaystyle \frac{347}{36}}\pi ^2Y38\zeta _3Y`$
$`{\displaystyle \frac{3}{4}}\pi ^2Y^223X\mathrm{Li}_3(x)+{\displaystyle \frac{11}{4}}X^2Y{\displaystyle \frac{7}{4}}X\mathrm{Li}_2(x)+42X\zeta _3+{\displaystyle \frac{21}{4}}X^2\mathrm{Li}_2(x)`$
$`{\displaystyle \frac{91}{72}}X\pi ^2+{\displaystyle \frac{11}{3}}UX^2{\displaystyle \frac{23}{6}}XY^3+{\displaystyle \frac{5}{12}}\pi ^2X^2+{\displaystyle \frac{21}{4}}\pi ^2\mathrm{Li}_2(x){\displaystyle \frac{22}{3}}UY^246\mathrm{Li}_3(y)X`$
$`{\displaystyle \frac{121}{9}}UX+{\displaystyle \frac{43}{24}}X^4+{\displaystyle \frac{347}{48}}X^2{\displaystyle \frac{179}{72}}X^3{\displaystyle \frac{19}{3}}X^3Y{\displaystyle \frac{1169}{72}}YX+{\displaystyle \frac{59}{12}}Y^2X`$
$`+{\displaystyle \frac{37}{8}}X^2Y^2+{\displaystyle \frac{31}{4}}\pi ^2XY+{\displaystyle \frac{22}{3}}UXY+{\displaystyle \frac{21}{2}}\mathrm{Li}_2(y)XY+{\displaystyle \frac{23}{12}}Y^4{\displaystyle \frac{59}{18}}Y^3`$
$`+{\displaystyle \frac{1169}{72}}Y^2+{\displaystyle \frac{5941}{216}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]`$
$`+({\displaystyle \frac{11}{3}}U\pi ^2+{\displaystyle \frac{207}{16}}X{\displaystyle \frac{319}{90}}\pi ^4+{\displaystyle \frac{1291}{144}}\pi ^2+{\displaystyle \frac{121}{12}}\zeta _3{\displaystyle \frac{87}{2}}\mathrm{Li}_4(z){\displaystyle \frac{121}{6}}\mathrm{Li}_3(y)`$
$`+74\mathrm{Li}_4(x)+{\displaystyle \frac{87}{2}}\mathrm{Li}_4(y){\displaystyle \frac{121}{12}}\mathrm{Li}_3(x)26\mathrm{Li}_3(x)Y+{\displaystyle \frac{163}{12}}\pi ^2Y52\mathrm{Li}_3(y)Y`$
$`+26\zeta _3Y{\displaystyle \frac{19}{4}}\pi ^2Y^221X\mathrm{Li}_3(x)+{\displaystyle \frac{163}{12}}X^2Y+{\displaystyle \frac{121}{12}}X\mathrm{Li}_2(x)34X\zeta _3`$
$`+{\displaystyle \frac{23}{4}}X^2\mathrm{Li}_2(x)+{\displaystyle \frac{43}{24}}X\pi ^2{\displaystyle \frac{11}{3}}UX^2+{\displaystyle \frac{1}{6}}XY^3{\displaystyle \frac{17}{4}}\pi ^2X^2+{\displaystyle \frac{121}{6}}\mathrm{Li}_2(y)Y`$
$`{\displaystyle \frac{95}{12}}\pi ^2\mathrm{Li}_2(x)+26\mathrm{Li}_3(y)X11UX{\displaystyle \frac{9}{8}}X^4+{\displaystyle \frac{1291}{144}}X^2{\displaystyle \frac{433}{72}}X^3+{\displaystyle \frac{17}{2}}\mathrm{Li}_2(y)Y^2`$
$`+7X^3Y{\displaystyle \frac{1291}{72}}YX{\displaystyle \frac{7}{2}}Y^2X{\displaystyle \frac{89}{8}}X^2Y^2+{\displaystyle \frac{25}{4}}\pi ^2XY+{\displaystyle \frac{22}{3}}UXY`$
$`{\displaystyle \frac{17}{2}}\mathrm{Li}_2(y)XY\left)\right[{\displaystyle \frac{t^2s^2}{st}}]`$
$`+(2U\zeta _3{\displaystyle \frac{1}{12}}U\pi ^2+{\displaystyle \frac{40}{3}}X{\displaystyle \frac{883}{360}}\pi ^4+{\displaystyle \frac{329}{36}}\pi ^2+{\displaystyle \frac{443}{18}}\zeta _3+{\displaystyle \frac{1502}{27}}U{\displaystyle \frac{19139}{324}}`$
$`+18\mathrm{Li}_4(z)+18\mathrm{Li}_4(y)+18\mathrm{Li}_3(x)Y+{\displaystyle \frac{40}{3}}\pi ^2Y18\zeta _3Y+4\pi ^2Y^218X\mathrm{Li}_3(x)`$
$`{\displaystyle \frac{40}{3}}X^2Y+18X\zeta _3+9X^2\mathrm{Li}_2(x)+{\displaystyle \frac{20}{3}}X\pi ^2+{\displaystyle \frac{11}{3}}\pi ^2X^2+9\pi ^2\mathrm{Li}_2(x)18\mathrm{Li}_3(y)X`$
$`+{\displaystyle \frac{1}{2}}X^4{\displaystyle \frac{517}{36}}X^2+{\displaystyle \frac{20}{3}}X^3+X^3Y+7YX+{\displaystyle \frac{19}{2}}X^2Y^2+\pi ^2XY`$
$`+18\mathrm{Li}_2(y)XY7Y^2{\displaystyle \frac{80}{3}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{u^2}}]`$
$`+({\displaystyle \frac{11}{6}}U\pi ^2{\displaystyle \frac{121}{36}}X{\displaystyle \frac{119}{180}}\pi ^4{\displaystyle \frac{157}{18}}\pi ^2{\displaystyle \frac{71}{6}}\zeta _3{\displaystyle \frac{43}{2}}\mathrm{Li}_4(z)+{\displaystyle \frac{71}{3}}\mathrm{Li}_3(y)`$
$`+22\mathrm{Li}_4(x)+{\displaystyle \frac{43}{2}}\mathrm{Li}_4(y)+{\displaystyle \frac{71}{6}}\mathrm{Li}_3(x)9\mathrm{Li}_3(x)Y{\displaystyle \frac{39}{2}}\pi ^2Y18\mathrm{Li}_3(y)Y+9\zeta _3Y`$
$`{\displaystyle \frac{17}{4}}\pi ^2Y^2+2X\mathrm{Li}_3(x){\displaystyle \frac{39}{2}}X^2Y{\displaystyle \frac{71}{6}}X\mathrm{Li}_2(x)X\zeta _3{\displaystyle \frac{9}{4}}X^2\mathrm{Li}_2(x)`$
$`{\displaystyle \frac{137}{18}}X\pi ^2+{\displaystyle \frac{11}{6}}UX^2{\displaystyle \frac{5}{2}}XY^3{\displaystyle \frac{11}{4}}\pi ^2X^2{\displaystyle \frac{71}{3}}\mathrm{Li}_2(y)Y{\displaystyle \frac{5}{4}}\pi ^2\mathrm{Li}_2(x)`$
$`+9\mathrm{Li}_3(y)X+{\displaystyle \frac{11}{6}}UX{\displaystyle \frac{43}{24}}X^4{\displaystyle \frac{157}{18}}X^2+{\displaystyle \frac{221}{18}}X^3{\displaystyle \frac{7}{2}}\mathrm{Li}_2(y)Y^2+{\displaystyle \frac{29}{6}}X^3Y`$
$`+{\displaystyle \frac{157}{9}}YX+{\displaystyle \frac{23}{3}}Y^2X{\displaystyle \frac{13}{8}}X^2Y^2{\displaystyle \frac{37}{12}}\pi ^2XY{\displaystyle \frac{11}{3}}UXY`$
$`+{\displaystyle \frac{7}{2}}\mathrm{Li}_2(y)XY\left)\right[{\displaystyle \frac{t^2s^2}{u^2}}]Y^2{\displaystyle \frac{s^4}{t^2u^2}}+(Y^2+2YX\pi ^2X^2){\displaystyle \frac{t^4}{s^2u^2}}`$
$`C_u`$ $`=`$ $`({\displaystyle \frac{22}{3}}U\pi ^2+{\displaystyle \frac{3133}{144}}X{\displaystyle \frac{203}{180}}\pi ^4+{\displaystyle \frac{401}{48}}\pi ^2{\displaystyle \frac{1}{4}}\zeta _3+{\displaystyle \frac{83}{2}}\mathrm{Li}_4(z)+{\displaystyle \frac{23}{8}}+{\displaystyle \frac{83}{2}}\mathrm{Li}_4(y)`$ (57)
$`+{\displaystyle \frac{13}{4}}\mathrm{Li}_3(x)+41\mathrm{Li}_3(x)Y+{\displaystyle \frac{44}{3}}UY+{\displaystyle \frac{41}{9}}\pi ^2Y+9\zeta _3Y{\displaystyle \frac{13}{12}}\pi ^2Y^226X\mathrm{Li}_3(x)`$
$`+{\displaystyle \frac{41}{3}}X^2Y{\displaystyle \frac{13}{4}}X\mathrm{Li}_2(x)+X\zeta _3+{\displaystyle \frac{21}{4}}X^2\mathrm{Li}_2(x){\displaystyle \frac{299}{72}}X\pi ^2{\displaystyle \frac{22}{3}}UX^2+{\displaystyle \frac{3}{2}}XY^3`$
$`+{\displaystyle \frac{25}{12}}\pi ^2X^2+{\displaystyle \frac{21}{4}}\pi ^2\mathrm{Li}_2(x){\displaystyle \frac{44}{3}}UY^211\mathrm{Li}_3(y)X{\displaystyle \frac{22}{3}}UX+{\displaystyle \frac{55}{24}}X^4+{\displaystyle \frac{919}{144}}X^2`$
$`{\displaystyle \frac{209}{72}}X^3{\displaystyle \frac{47}{6}}X^3Y{\displaystyle \frac{259}{72}}YX{\displaystyle \frac{161}{6}}Y^2X+{\displaystyle \frac{125}{8}}X^2Y^2{\displaystyle \frac{33}{4}}\pi ^2XY`$
$`+{\displaystyle \frac{44}{3}}UXY+{\displaystyle \frac{21}{2}}\mathrm{Li}_2(y)XY{\displaystyle \frac{3}{4}}Y^4+{\displaystyle \frac{161}{9}}Y^3+{\displaystyle \frac{259}{72}}Y^2{\displaystyle \frac{3133}{72}}Y)`$
$`+(5U\zeta _3{\displaystyle \frac{145}{24}}U\pi ^2+{\displaystyle \frac{743}{48}}X{\displaystyle \frac{113}{240}}\pi ^4+{\displaystyle \frac{1631}{144}}\pi ^2{\displaystyle \frac{4}{9}}\zeta _3+{\displaystyle \frac{5297}{216}}U+{\displaystyle \frac{19}{2}}\mathrm{Li}_4(z)`$
$`+{\displaystyle \frac{19}{2}}\mathrm{Li}_4(y){\displaystyle \frac{41393}{2592}}+{\displaystyle \frac{21}{4}}\mathrm{Li}_3(x)+19\mathrm{Li}_3(x)Y+11UY+{\displaystyle \frac{5}{2}}\pi ^2Y\zeta _3Y`$
$`+{\displaystyle \frac{5}{12}}\pi ^2Y^27X\mathrm{Li}_3(x)+{\displaystyle \frac{17}{2}}X^2Y{\displaystyle \frac{21}{4}}X\mathrm{Li}_2(x)2X\zeta _3+{\displaystyle \frac{9}{4}}X^2\mathrm{Li}_2(x)`$
$`{\displaystyle \frac{27}{8}}X\pi ^2{\displaystyle \frac{11}{2}}UX^2+{\displaystyle \frac{5}{6}}XY^3+{\displaystyle \frac{3}{4}}\pi ^2X^2+{\displaystyle \frac{9}{4}}\pi ^2\mathrm{Li}_2(x)11UY^2+5\mathrm{Li}_3(y)X`$
$`{\displaystyle \frac{11}{2}}UX+{\displaystyle \frac{7}{12}}X^4+{\displaystyle \frac{97}{48}}X^2{\displaystyle \frac{17}{8}}X^3{\displaystyle \frac{4}{3}}X^3Y{\displaystyle \frac{145}{24}}YX18Y^2X+{\displaystyle \frac{57}{8}}X^2Y^2`$
$`{\displaystyle \frac{61}{12}}\pi ^2XY+11UXY+{\displaystyle \frac{9}{2}}\mathrm{Li}_2(y)XY{\displaystyle \frac{5}{12}}Y^4+12Y^3+{\displaystyle \frac{145}{24}}Y^2`$
$`{\displaystyle \frac{743}{24}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]`$
$`+({\displaystyle \frac{11}{6}}U\pi ^2+{\displaystyle \frac{67}{16}}X+{\displaystyle \frac{14}{45}}\pi ^4+{\displaystyle \frac{1141}{144}}\pi ^2+{\displaystyle \frac{9}{4}}\zeta _3+{\displaystyle \frac{21}{2}}\mathrm{Li}_4(z){\displaystyle \frac{9}{2}}\mathrm{Li}_3(y)8\mathrm{Li}_4(x)`$
$`{\displaystyle \frac{21}{2}}\mathrm{Li}_4(y){\displaystyle \frac{9}{4}}\mathrm{Li}_3(x)+5\mathrm{Li}_3(x)Y+{\displaystyle \frac{17}{3}}\pi ^2Y+10\mathrm{Li}_3(y)Y5\zeta _3Y+{\displaystyle \frac{3}{4}}\pi ^2Y^2`$
$`3X\mathrm{Li}_3(x)+{\displaystyle \frac{17}{3}}X^2Y+{\displaystyle \frac{9}{4}}X\mathrm{Li}_2(x)+6X\zeta _3+{\displaystyle \frac{7}{4}}X^2\mathrm{Li}_2(x)+{\displaystyle \frac{19}{72}}X\pi ^2{\displaystyle \frac{11}{6}}UX^2`$
$`+{\displaystyle \frac{3}{2}}XY^3+{\displaystyle \frac{3}{4}}\pi ^2X^2+{\displaystyle \frac{9}{2}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{5}{12}}\pi ^2\mathrm{Li}_2(x)5\mathrm{Li}_3(y)X+{\displaystyle \frac{11}{2}}UX+{\displaystyle \frac{1}{4}}X^4`$
$`+{\displaystyle \frac{1141}{144}}X^2{\displaystyle \frac{107}{72}}X^3+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y^2{\displaystyle \frac{2}{3}}X^3Y{\displaystyle \frac{1141}{72}}YX{\displaystyle \frac{41}{12}}Y^2X+{\displaystyle \frac{3}{8}}X^2Y^2`$
$`{\displaystyle \frac{1}{4}}\pi ^2XY+{\displaystyle \frac{11}{3}}UXY{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)XY\left)\right[{\displaystyle \frac{t^2s^2}{st}}]`$
$`+({\displaystyle \frac{3}{4}}U{\displaystyle \frac{255}{16}}+15\zeta _3+{\displaystyle \frac{11}{90}}\pi ^4{\displaystyle \frac{35}{12}}\pi ^2+{\displaystyle \frac{9}{4}}X^2+YXY^2+U\pi ^2`$
$`12U\zeta _3\left)\right[{\displaystyle \frac{t^2+s^2}{u^2}}]`$
$`+({\displaystyle \frac{21}{16}}X{\displaystyle \frac{47}{180}}\pi ^4{\displaystyle \frac{47}{16}}\pi ^2{\displaystyle \frac{11}{4}}\zeta _3{\displaystyle \frac{11}{2}}\mathrm{Li}_4(z)+{\displaystyle \frac{11}{2}}\mathrm{Li}_3(y)+6\mathrm{Li}_4(x)+{\displaystyle \frac{11}{2}}\mathrm{Li}_4(y)`$
$`+{\displaystyle \frac{11}{4}}\mathrm{Li}_3(x)3\mathrm{Li}_3(x)Y{\displaystyle \frac{3}{2}}\pi ^2Y6\mathrm{Li}_3(y)Y+3\zeta _3Y{\displaystyle \frac{1}{4}}\pi ^2Y^2{\displaystyle \frac{3}{2}}X^2Y`$
$`{\displaystyle \frac{11}{4}}X\mathrm{Li}_2(x)3X\zeta _3{\displaystyle \frac{1}{4}}X^2\mathrm{Li}_2(x)+{\displaystyle \frac{17}{24}}X\pi ^2{\displaystyle \frac{1}{2}}XY^3{\displaystyle \frac{5}{12}}\pi ^2X^2{\displaystyle \frac{11}{2}}\mathrm{Li}_2(y)Y`$
$`{\displaystyle \frac{7}{12}}\pi ^2\mathrm{Li}_2(x)+3\mathrm{Li}_3(y)X{\displaystyle \frac{1}{8}}X^4{\displaystyle \frac{47}{16}}X^2{\displaystyle \frac{13}{24}}X^3+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y^2+{\displaystyle \frac{1}{2}}X^3Y`$
$`+{\displaystyle \frac{47}{8}}YX{\displaystyle \frac{5}{4}}Y^2X{\displaystyle \frac{5}{8}}X^2Y^2+{\displaystyle \frac{1}{4}}\pi ^2XY{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)XY\left)\right[{\displaystyle \frac{t^2s^2}{u^2}}]`$
$`+Y^2{\displaystyle \frac{s^4}{t^2u^2}}+\left(Y^22YX+\pi ^2+X^2\right){\displaystyle \frac{t^4}{s^2u^2}}`$
$`D_u`$ $`=`$ $`({\displaystyle \frac{1}{2}}+{\displaystyle \frac{29}{4}}X+{\displaystyle \frac{3}{5}}\pi ^4+{\displaystyle \frac{25}{6}}\pi ^2+15\zeta _310\mathrm{Li}_4(z)10\mathrm{Li}_4(y)19\mathrm{Li}_3(x)`$ (58)
$`+12\mathrm{Li}_3(x)Y+{\displaystyle \frac{7}{3}}\pi ^2Y+12\zeta _3Y+{\displaystyle \frac{1}{3}}\pi ^2Y^2+4X\mathrm{Li}_3(x)+13X^2Y+19X\mathrm{Li}_2(x)`$
$`16X\zeta _3+X^2\mathrm{Li}_2(x){\displaystyle \frac{1}{6}}X\pi ^2+2XY^3{\displaystyle \frac{2}{3}}\pi ^2X^2+\pi ^2\mathrm{Li}_2(x)+20\mathrm{Li}_3(y)X`$
$`X^4+4X^2{\displaystyle \frac{11}{6}}X^3+4X^3YYX9Y^2X+{\displaystyle \frac{5}{2}}X^2Y^2{\displaystyle \frac{7}{3}}\pi ^2XY`$
$`+2\mathrm{Li}_2(y)XYY^4+6Y^3+Y^2{\displaystyle \frac{29}{2}}Y)`$
$`+(6U\zeta _3+{\displaystyle \frac{1}{2}}U\pi ^2+{\displaystyle \frac{47}{16}}X+{\displaystyle \frac{7}{18}}\pi ^4{\displaystyle \frac{3}{16}}\pi ^2+{\displaystyle \frac{23}{4}}\zeta _3{\displaystyle \frac{3}{8}}U{\displaystyle \frac{11}{2}}\mathrm{Li}_4(z)`$
$`{\displaystyle \frac{11}{2}}\mathrm{Li}_4(y){\displaystyle \frac{13}{4}}\mathrm{Li}_3(x){\displaystyle \frac{187}{32}}+7\mathrm{Li}_3(x)Y+{\displaystyle \frac{7}{6}}\pi ^2Y+7\zeta _3Y+{\displaystyle \frac{5}{12}}\pi ^2Y^2`$
$`+2X\mathrm{Li}_3(x)+{\displaystyle \frac{7}{2}}X^2Y+{\displaystyle \frac{13}{4}}X\mathrm{Li}_2(x)9X\zeta _3+{\displaystyle \frac{3}{4}}X^2\mathrm{Li}_2(x)+{\displaystyle \frac{1}{24}}X\pi ^2+{\displaystyle \frac{7}{6}}XY^3`$
$`{\displaystyle \frac{5}{12}}\pi ^2X^2+{\displaystyle \frac{3}{4}}\pi ^2\mathrm{Li}_2(x)+11\mathrm{Li}_3(y)X{\displaystyle \frac{5}{8}}X^4{\displaystyle \frac{13}{16}}X^2{\displaystyle \frac{3}{8}}X^3+{\displaystyle \frac{5}{2}}X^3Y+{\displaystyle \frac{13}{8}}YX`$
$`5Y^2X+{\displaystyle \frac{11}{8}}X^2Y^2{\displaystyle \frac{13}{12}}\pi ^2XY+{\displaystyle \frac{3}{2}}\mathrm{Li}_2(y)XY{\displaystyle \frac{7}{12}}Y^4+{\displaystyle \frac{10}{3}}Y^3{\displaystyle \frac{13}{8}}Y^2`$
$`{\displaystyle \frac{47}{8}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]`$
$`+({\displaystyle \frac{31}{16}}X+{\displaystyle \frac{103}{180}}\pi ^4{\displaystyle \frac{23}{16}}\pi ^2+{\displaystyle \frac{19}{4}}\zeta _3+{\displaystyle \frac{11}{2}}\mathrm{Li}_4(z){\displaystyle \frac{19}{2}}\mathrm{Li}_3(y)14\mathrm{Li}_4(x)`$
$`{\displaystyle \frac{11}{2}}\mathrm{Li}_4(y){\displaystyle \frac{19}{4}}\mathrm{Li}_3(x)+3\mathrm{Li}_3(x)Y+2\pi ^2Y+6\mathrm{Li}_3(y)Y3\zeta _3Y+{\displaystyle \frac{1}{4}}\pi ^2Y^2`$
$`+4X\mathrm{Li}_3(x)+2X^2Y+{\displaystyle \frac{19}{4}}X\mathrm{Li}_2(x)+7X\zeta _3+{\displaystyle \frac{1}{4}}X^2\mathrm{Li}_2(x)+{\displaystyle \frac{7}{24}}X\pi ^2+{\displaystyle \frac{1}{2}}XY^3`$
$`+{\displaystyle \frac{3}{4}}\pi ^2X^2+{\displaystyle \frac{19}{2}}\mathrm{Li}_2(y)Y{\displaystyle \frac{3}{4}}\pi ^2\mathrm{Li}_2(x)3\mathrm{Li}_3(y)X{\displaystyle \frac{1}{24}}X^4{\displaystyle \frac{23}{16}}X^2+{\displaystyle \frac{3}{8}}X^3`$
$`{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)Y^2+{\displaystyle \frac{1}{6}}X^3Y+{\displaystyle \frac{23}{8}}YX+{\displaystyle \frac{11}{4}}Y^2X+{\displaystyle \frac{5}{8}}X^2Y^2{\displaystyle \frac{19}{12}}\pi ^2XY`$
$`+{\displaystyle \frac{1}{2}}\mathrm{Li}_2(y)XY\left)\right[{\displaystyle \frac{t^2s^2}{st}}]+(2X^2YX+{\displaystyle \frac{1}{2}}\pi ^2+Y^2\left)\right[{\displaystyle \frac{t^2+s^2}{u^2}}]`$
$`+\left({\displaystyle \frac{3}{2}}X^23YX+{\displaystyle \frac{3}{2}}\pi ^23X\right)\left[{\displaystyle \frac{t^2s^2}{u^2}}\right]Y^2{\displaystyle \frac{s^4}{t^2u^2}}`$
$`+\left(Y^2+2YX\pi ^2X^2\right){\displaystyle \frac{t^4}{s^2u^2}}`$
$`E_u`$ $`=`$ $`({\displaystyle \frac{91}{18}}X{\displaystyle \frac{11}{180}}\pi ^4{\displaystyle \frac{25}{18}}\pi ^24\zeta _3+{\displaystyle \frac{2}{3}}U4\mathrm{Li}_3(x)2\mathrm{Li}_3(x)Y+{\displaystyle \frac{10}{3}}UY`$ (59)
$`+{\displaystyle \frac{2}{3}}\pi ^2Y2\zeta _3Y+{\displaystyle \frac{1}{6}}\pi ^2Y^2{\displaystyle \frac{37}{18}}+2X^2Y+4X\mathrm{Li}_2(x)+2X\zeta _3+{\displaystyle \frac{2}{3}}X\pi ^2`$
$`{\displaystyle \frac{1}{3}}XY^32\mathrm{Li}_3(y)X{\displaystyle \frac{5}{3}}UX2X^2+4YX2Y^2X{\displaystyle \frac{1}{2}}X^2Y^2+{\displaystyle \frac{2}{3}}\pi ^2XY`$
$`+{\displaystyle \frac{1}{6}}Y^4+{\displaystyle \frac{4}{3}}Y^34Y^2{\displaystyle \frac{91}{9}}Y)`$
$`+({\displaystyle \frac{44}{9}}U^2+{\displaystyle \frac{11}{36}}U\pi ^2+{\displaystyle \frac{5}{9}}X+{\displaystyle \frac{11}{90}}\pi ^4+{\displaystyle \frac{133}{108}}\pi ^2{\displaystyle \frac{13}{6}}\zeta _3+{\displaystyle \frac{185}{54}}U{\displaystyle \frac{2}{3}}\mathrm{Li}_3(x)`$
$`+4\mathrm{Li}_3(x)Y+{\displaystyle \frac{44}{9}}UY+{\displaystyle \frac{13}{36}}\pi ^2Y+4\zeta _3Y{\displaystyle \frac{1}{3}}\pi ^2Y^2{\displaystyle \frac{1}{4}}X^2Y+{\displaystyle \frac{2}{3}}X\mathrm{Li}_2(x)`$
$`4X\zeta _3+{\displaystyle \frac{5}{18}}X\pi ^2+{\displaystyle \frac{1}{3}}UX^2+{\displaystyle \frac{2}{3}}XY^3+{\displaystyle \frac{2}{3}}UY^2+4\mathrm{Li}_3(y)X{\displaystyle \frac{22}{9}}UX`$
$`{\displaystyle \frac{7}{36}}X^2+{\displaystyle \frac{7}{36}}X^3+{\displaystyle \frac{7}{18}}YX+{\displaystyle \frac{1}{4}}Y^2X+X^2Y^2{\displaystyle \frac{4}{3}}\pi ^2XY{\displaystyle \frac{2}{3}}UXY`$
$`{\displaystyle \frac{1}{3}}Y^4{\displaystyle \frac{1}{6}}Y^3{\displaystyle \frac{7}{18}}Y^2+{\displaystyle \frac{1307}{216}}{\displaystyle \frac{10}{9}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]`$
$`+({\displaystyle \frac{1}{3}}U\pi ^2{\displaystyle \frac{5}{9}}X+{\displaystyle \frac{11}{30}}\pi ^4+{\displaystyle \frac{79}{36}}\pi ^2{\displaystyle \frac{2}{3}}\zeta _3+8\mathrm{Li}_4(z)+{\displaystyle \frac{4}{3}}\mathrm{Li}_3(y)8\mathrm{Li}_4(x)`$
$`8\mathrm{Li}_4(y)+{\displaystyle \frac{2}{3}}\mathrm{Li}_3(x)+4\mathrm{Li}_3(x)Y+{\displaystyle \frac{1}{4}}\pi ^2Y+8\mathrm{Li}_3(y)Y4\zeta _3Y+\pi ^2Y^2`$
$`+{\displaystyle \frac{1}{4}}X^2Y{\displaystyle \frac{2}{3}}X\mathrm{Li}_2(x)+4X\zeta _3{\displaystyle \frac{5}{18}}X\pi ^2{\displaystyle \frac{1}{3}}UX^2+{\displaystyle \frac{2}{3}}XY^3+{\displaystyle \frac{2}{3}}\pi ^2X^2`$
$`{\displaystyle \frac{4}{3}}\mathrm{Li}_2(y)Y+{\displaystyle \frac{4}{3}}\pi ^2\mathrm{Li}_2(x)4\mathrm{Li}_3(y)X+{\displaystyle \frac{22}{9}}UX+{\displaystyle \frac{1}{3}}X^4+{\displaystyle \frac{79}{36}}X^2{\displaystyle \frac{7}{36}}X^3`$
$`{\displaystyle \frac{4}{3}}X^3Y{\displaystyle \frac{79}{18}}YX{\displaystyle \frac{11}{12}}Y^2X+X^2Y^2+{\displaystyle \frac{2}{3}}UXY\left)\right[{\displaystyle \frac{t^2s^2}{st}}]`$
$`+({\displaystyle \frac{88}{9}}U^2{\displaystyle \frac{28}{9}}X{\displaystyle \frac{221}{27}}U+{\displaystyle \frac{4}{3}}\mathrm{Li}_3(x)+{\displaystyle \frac{13}{3}}\zeta _3{\displaystyle \frac{67}{54}}\pi ^2{\displaystyle \frac{124}{9}}UY{\displaystyle \frac{17}{18}}X\pi ^2`$
$`+{\displaystyle \frac{5}{3}}X^2Y{\displaystyle \frac{4}{3}}X\mathrm{Li}_2(x){\displaystyle \frac{4}{3}}UX^2{\displaystyle \frac{8}{3}}UY^2+{\displaystyle \frac{62}{9}}UX+{\displaystyle \frac{71}{18}}X^2{\displaystyle \frac{7}{9}}X^3{\displaystyle \frac{71}{9}}YX`$
$`{\displaystyle \frac{5}{3}}Y^2X+{\displaystyle \frac{8}{3}}UXY+{\displaystyle \frac{10}{9}}Y^3+{\displaystyle \frac{71}{9}}Y^2{\displaystyle \frac{10}{9}}\pi ^2Y{\displaystyle \frac{863}{108}}{\displaystyle \frac{23}{18}}U\pi ^2`$
$`+{\displaystyle \frac{56}{9}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{u^2}}]`$
$`+({\displaystyle \frac{37}{18}}X{\displaystyle \frac{11}{60}}\pi ^4{\displaystyle \frac{23}{6}}\pi ^2{\displaystyle \frac{8}{3}}\zeta _34\mathrm{Li}_4(z)+{\displaystyle \frac{16}{3}}\mathrm{Li}_3(y)+4\mathrm{Li}_4(x)+4\mathrm{Li}_4(y)`$
$`+{\displaystyle \frac{8}{3}}\mathrm{Li}_3(x)2\mathrm{Li}_3(x)Y{\displaystyle \frac{4}{3}}\pi ^2Y4\mathrm{Li}_3(y)Y+2\zeta _3Y{\displaystyle \frac{1}{2}}\pi ^2Y^2{\displaystyle \frac{4}{3}}X^2Y`$
$`{\displaystyle \frac{8}{3}}X\mathrm{Li}_2(x)2X\zeta _3{\displaystyle \frac{1}{2}}X\pi ^2{\displaystyle \frac{1}{3}}XY^3{\displaystyle \frac{1}{3}}\pi ^2X^2{\displaystyle \frac{16}{3}}\mathrm{Li}_2(y)Y{\displaystyle \frac{2}{3}}\pi ^2\mathrm{Li}_2(x)`$
$`+2\mathrm{Li}_3(y)X{\displaystyle \frac{47}{9}}UX{\displaystyle \frac{1}{6}}X^4{\displaystyle \frac{23}{6}}X^2+{\displaystyle \frac{2}{3}}X^3Y+{\displaystyle \frac{23}{3}}YX{\displaystyle \frac{4}{3}}Y^2X`$
$`{\displaystyle \frac{1}{2}}X^2Y^2\left)\right[{\displaystyle \frac{t^2s^2}{u^2}}]`$
$`F_u`$ $`=`$ $`({\displaystyle \frac{1}{3}}U\pi ^2{\displaystyle \frac{50}{9}}X{\displaystyle \frac{11}{45}}\pi ^4+{\displaystyle \frac{41}{9}}\pi ^2+{\displaystyle \frac{52}{3}}\zeta _3+{\displaystyle \frac{2}{3}}U{\displaystyle \frac{19}{18}}+{\displaystyle \frac{52}{3}}\mathrm{Li}_3(x)8\mathrm{Li}_3(x)Y`$ (60)
$`{\displaystyle \frac{10}{3}}UY{\displaystyle \frac{85}{36}}\pi ^2Y8\zeta _3Y+{\displaystyle \frac{2}{3}}\pi ^2Y^2{\displaystyle \frac{83}{12}}X^2Y{\displaystyle \frac{52}{3}}X\mathrm{Li}_2(x)+8X\zeta _3`$
$`{\displaystyle \frac{115}{36}}X\pi ^2+{\displaystyle \frac{1}{3}}UX^2{\displaystyle \frac{4}{3}}XY^3+{\displaystyle \frac{2}{3}}UY^28\mathrm{Li}_3(y)X+{\displaystyle \frac{5}{3}}UX+{\displaystyle \frac{49}{18}}X^2`$
$`{\displaystyle \frac{11}{36}}X^3{\displaystyle \frac{43}{9}}YX+{\displaystyle \frac{21}{4}}Y^2X2X^2Y^2+{\displaystyle \frac{8}{3}}\pi ^2XY{\displaystyle \frac{2}{3}}UXY+{\displaystyle \frac{2}{3}}Y^4{\displaystyle \frac{7}{2}}Y^3`$
$`+{\displaystyle \frac{43}{9}}Y^2+{\displaystyle \frac{100}{9}}Y)`$
$`+({\displaystyle \frac{44}{9}}U^2{\displaystyle \frac{3661}{324}}+{\displaystyle \frac{7}{9}}U\pi ^2{\displaystyle \frac{115}{36}}X{\displaystyle \frac{11}{45}}\pi ^4{\displaystyle \frac{367}{108}}\pi ^2+{\displaystyle \frac{14}{9}}\zeta _3+{\displaystyle \frac{7}{9}}U{\displaystyle \frac{8}{3}}\mathrm{Li}_3(x)`$
$`8\mathrm{Li}_3(x)Y{\displaystyle \frac{88}{9}}UY{\displaystyle \frac{1}{12}}\pi ^2Y8\zeta _3Y+{\displaystyle \frac{2}{3}}\pi ^2Y^2+{\displaystyle \frac{1}{4}}X^2Y+{\displaystyle \frac{8}{3}}X\mathrm{Li}_2(x)`$
$`+8X\zeta _3{\displaystyle \frac{2}{3}}UX^2{\displaystyle \frac{4}{3}}XY^3+{\displaystyle \frac{4}{3}}UY^28\mathrm{Li}_3(y)X+{\displaystyle \frac{44}{9}}UX{\displaystyle \frac{11}{9}}X^2{\displaystyle \frac{1}{18}}X^3`$
$`+{\displaystyle \frac{53}{18}}YX+{\displaystyle \frac{9}{4}}Y^2X2X^2Y^2+{\displaystyle \frac{8}{3}}\pi ^2XY{\displaystyle \frac{4}{3}}UXY+{\displaystyle \frac{2}{3}}Y^4{\displaystyle \frac{3}{2}}Y^3{\displaystyle \frac{53}{18}}Y^2`$
$`+{\displaystyle \frac{115}{18}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]`$
$`+({\displaystyle \frac{2}{3}}U\pi ^2{\displaystyle \frac{23}{4}}X{\displaystyle \frac{11}{15}}\pi ^4{\displaystyle \frac{19}{6}}\pi ^2+{\displaystyle \frac{4}{3}}\zeta _316\mathrm{Li}_4(z){\displaystyle \frac{8}{3}}\mathrm{Li}_3(y)+16\mathrm{Li}_4(x)`$
$`+16\mathrm{Li}_4(y){\displaystyle \frac{4}{3}}\mathrm{Li}_3(x)8\mathrm{Li}_3(x)Y{\displaystyle \frac{19}{12}}\pi ^2Y16\mathrm{Li}_3(y)Y+8\zeta _3Y2\pi ^2Y^2`$
$`{\displaystyle \frac{19}{12}}X^2Y+{\displaystyle \frac{4}{3}}X\mathrm{Li}_2(x)8X\zeta _3{\displaystyle \frac{1}{18}}X\pi ^2+{\displaystyle \frac{2}{3}}UX^2{\displaystyle \frac{4}{3}}XY^3{\displaystyle \frac{4}{3}}\pi ^2X^2`$
$`+{\displaystyle \frac{8}{3}}\mathrm{Li}_2(y)Y{\displaystyle \frac{8}{3}}\pi ^2\mathrm{Li}_2(x)+8\mathrm{Li}_3(y)X+2UX{\displaystyle \frac{2}{3}}X^4{\displaystyle \frac{19}{6}}X^2+{\displaystyle \frac{13}{18}}X^3+{\displaystyle \frac{8}{3}}X^3Y`$
$`+{\displaystyle \frac{19}{3}}YX+{\displaystyle \frac{35}{12}}Y^2X2X^2Y^2{\displaystyle \frac{4}{3}}UXY\left)\right[{\displaystyle \frac{t^2s^2}{st}}]`$
$`+({\displaystyle \frac{4}{3}}X{\displaystyle \frac{263}{27}}U{\displaystyle \frac{1}{9}}\zeta _3{\displaystyle \frac{1}{3}}\pi ^2{\displaystyle \frac{2}{3}}X\pi ^2+{\displaystyle \frac{4}{3}}X^2Y+{\displaystyle \frac{85}{18}}X^2{\displaystyle \frac{2}{3}}X^33YX`$
$`+3Y^2{\displaystyle \frac{4}{3}}\pi ^2Y+{\displaystyle \frac{4085}{324}}{\displaystyle \frac{1}{6}}U\pi ^2+{\displaystyle \frac{8}{3}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{u^2}}]`$
$`+({\displaystyle \frac{19}{9}}X+{\displaystyle \frac{8}{3}}\mathrm{Li}_3(y)+{\displaystyle \frac{4}{3}}\mathrm{Li}_3(x){\displaystyle \frac{4}{3}}\zeta _3+7\pi ^2+{\displaystyle \frac{5}{12}}X\pi ^2+{\displaystyle \frac{9}{4}}X^2Y{\displaystyle \frac{4}{3}}X\mathrm{Li}_2(x)`$
$`{\displaystyle \frac{1}{3}}UX^2{\displaystyle \frac{8}{3}}\mathrm{Li}_2(y)Y{\displaystyle \frac{1}{3}}UX+7X^2{\displaystyle \frac{49}{36}}X^314YX{\displaystyle \frac{43}{12}}Y^2X`$
$`+{\displaystyle \frac{2}{3}}UXY+{\displaystyle \frac{9}{4}}\pi ^2Y{\displaystyle \frac{1}{3}}U\pi ^2\left)\right[{\displaystyle \frac{t^2s^2}{u^2}}]3Y^2{\displaystyle \frac{s^4}{t^2u^2}}`$
$`+\left(3\pi ^23X^23Y^2+6YX\right){\displaystyle \frac{t^4}{s^2u^2}}`$
$`G_u`$ $`=`$ $`({\displaystyle \frac{4}{3}}U\pi ^2+{\displaystyle \frac{71}{9}}X+{\displaystyle \frac{11}{15}}\pi ^420\pi ^252\zeta _352\mathrm{Li}_3(x)+24\mathrm{Li}_3(x)Y{\displaystyle \frac{8}{3}}UY`$ (61)
$`+{\displaystyle \frac{64}{9}}\pi ^2Y+24\zeta _3Y2\pi ^2Y^2+{\displaystyle \frac{64}{3}}X^2Y+52X\mathrm{Li}_2(x)24X\zeta _3+{\displaystyle \frac{85}{9}}X\pi ^2`$
$`+{\displaystyle \frac{4}{3}}UX^2+4XY^3+{\displaystyle \frac{8}{3}}UY^2+24\mathrm{Li}_3(y)X+{\displaystyle \frac{4}{3}}UX{\displaystyle \frac{37}{9}}X^2+{\displaystyle \frac{7}{9}}X^3+{\displaystyle \frac{86}{9}}YX`$
$`{\displaystyle \frac{50}{3}}Y^2X+6X^2Y^28\pi ^2XY{\displaystyle \frac{8}{3}}UXY2Y^4+{\displaystyle \frac{100}{9}}Y^3{\displaystyle \frac{86}{9}}Y^2{\displaystyle \frac{142}{9}}Y)`$
$`+({\displaystyle \frac{11}{12}}U\pi ^2+{\displaystyle \frac{79}{12}}X+{\displaystyle \frac{11}{15}}\pi ^4+{\displaystyle \frac{151}{36}}\pi ^2+{\displaystyle \frac{71}{18}}\zeta _3+{\displaystyle \frac{3401}{648}}{\displaystyle \frac{227}{54}}U+4\mathrm{Li}_3(x)`$
$`+24\mathrm{Li}_3(x)Y2UY+{\displaystyle \frac{5}{3}}\pi ^2Y+24\zeta _3Y2\pi ^2Y^22X^2Y4X\mathrm{Li}_2(x)`$
$`24X\zeta _3{\displaystyle \frac{1}{12}}X\pi ^2+UX^2+4XY^3+2UY^2+24\mathrm{Li}_3(y)X+UX+{\displaystyle \frac{31}{6}}X^2`$
$`+{\displaystyle \frac{7}{12}}X^3{\displaystyle \frac{65}{6}}YX{\displaystyle \frac{3}{2}}Y^2X+6X^2Y^28\pi ^2XY2UXY2Y^4+Y^3`$
$`+{\displaystyle \frac{65}{6}}Y^2{\displaystyle \frac{79}{6}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]`$
$`+({\displaystyle \frac{1}{3}}U\pi ^2+{\displaystyle \frac{41}{4}}X+{\displaystyle \frac{11}{5}}\pi ^4+{\displaystyle \frac{85}{18}}\pi ^24\zeta _3+48\mathrm{Li}_4(z)+8\mathrm{Li}_3(y)48\mathrm{Li}_4(x)`$
$`48\mathrm{Li}_4(y)+4\mathrm{Li}_3(x)+24\mathrm{Li}_3(x)Y{\displaystyle \frac{2}{3}}\pi ^2Y+48\mathrm{Li}_3(y)Y24\zeta _3Y+6\pi ^2Y^2`$
$`{\displaystyle \frac{2}{3}}X^2Y4X\mathrm{Li}_2(x)+24X\zeta _3{\displaystyle \frac{17}{36}}X\pi ^2+{\displaystyle \frac{1}{3}}UX^2+4XY^3+4\pi ^2X^2`$
$`8\mathrm{Li}_2(y)Y+8\pi ^2\mathrm{Li}_2(x)24\mathrm{Li}_3(y)XUX+2X^4+{\displaystyle \frac{85}{18}}X^2+{\displaystyle \frac{7}{36}}X^38X^3Y`$
$`{\displaystyle \frac{85}{9}}YX{\displaystyle \frac{10}{3}}Y^2X+6X^2Y^2{\displaystyle \frac{2}{3}}UXY\left)\right[{\displaystyle \frac{t^2s^2}{st}}]`$
$`+\left({\displaystyle \frac{9}{2}}X^2+9YX{\displaystyle \frac{9}{2}}\pi ^29Y^2\right)\left[{\displaystyle \frac{t^2+s^2}{u^2}}\right]`$
$`+(4X\mathrm{Li}_2(x)4\mathrm{Li}_3(x)8\mathrm{Li}_3(y)+4\zeta _3{\displaystyle \frac{27}{2}}X^2+27YX+4Y^2X+{\displaystyle \frac{2}{3}}X\pi ^2`$
$`+8\mathrm{Li}_2(y)Y{\displaystyle \frac{27}{2}}\pi ^28X\left)\right[{\displaystyle \frac{t^2s^2}{u^2}}]+9Y^2{\displaystyle \frac{s^4}{t^2u^2}}`$
$`+\left(9X^2+9\pi ^218YX+9Y^2\right){\displaystyle \frac{t^4}{s^2u^2}}`$
$`H_u`$ $`=`$ $`({\displaystyle \frac{1}{12}}X^2{\displaystyle \frac{5}{27}}X+{\displaystyle \frac{1}{6}}Y^2{\displaystyle \frac{1}{6}}YX{\displaystyle \frac{4}{9}}UY{\displaystyle \frac{20}{27}}U{\displaystyle \frac{7}{108}}\pi ^2+{\displaystyle \frac{2}{9}}UX+{\displaystyle \frac{4}{9}}U^2`$ (62)
$`+{\displaystyle \frac{10}{27}}Y\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]`$
$`+\left({\displaystyle \frac{1}{12}}\pi ^2{\displaystyle \frac{1}{12}}X^2+{\displaystyle \frac{5}{27}}X+{\displaystyle \frac{1}{6}}YX{\displaystyle \frac{2}{9}}UX\right)\left[{\displaystyle \frac{t^2s^2}{st}}\right]`$
$`+({\displaystyle \frac{8}{9}}U^2+{\displaystyle \frac{40}{27}}U{\displaystyle \frac{1}{6}}X^2{\displaystyle \frac{20}{27}}Y+{\displaystyle \frac{1}{3}}YX+{\displaystyle \frac{7}{54}}\pi ^2{\displaystyle \frac{1}{3}}Y^2+{\displaystyle \frac{10}{27}}X`$
$`+{\displaystyle \frac{8}{9}}UY{\displaystyle \frac{4}{9}}UX\left)\right[{\displaystyle \frac{t^2+s^2}{u^2}}]`$
$`+\left({\displaystyle \frac{1}{6}}X^2{\displaystyle \frac{1}{3}}YX+{\displaystyle \frac{1}{6}}\pi ^2{\displaystyle \frac{10}{27}}X+{\displaystyle \frac{4}{9}}UX\right)\left[{\displaystyle \frac{t^2s^2}{u^2}}\right]`$
$`I_u`$ $`=`$ $`({\displaystyle \frac{8}{9}}UY+{\displaystyle \frac{1}{3}}YX{\displaystyle \frac{4}{9}}UX+{\displaystyle \frac{10}{27}}X+{\displaystyle \frac{5}{54}}\pi ^2{\displaystyle \frac{1}{9}}X^2+{\displaystyle \frac{20}{27}}U{\displaystyle \frac{20}{27}}Y`$ (63)
$`{\displaystyle \frac{1}{3}}Y^2{\displaystyle \frac{4}{9}}U^2\left)\right[{\displaystyle \frac{t^2+s^2}{st}}]+({\displaystyle \frac{1}{18}}\pi ^2+{\displaystyle \frac{1}{18}}X^2{\displaystyle \frac{1}{9}}YX\left)\right[{\displaystyle \frac{t^2s^2}{st}}]`$
$`{\displaystyle \frac{1}{9}}\left(X^2\pi ^2\right)\left[{\displaystyle \frac{t^2+s^2}{u^2}}\right]+\left({\displaystyle \frac{2}{9}}YX{\displaystyle \frac{1}{9}}\pi ^2{\displaystyle \frac{1}{9}}X^2\right)\left[{\displaystyle \frac{t^2s^2}{u^2}}\right]`$
## 4 One-loop self-interference

In this section, we present the self-interference of the one-loop amplitude of Eq. (21) in terms of the one-loop box in $`D=62ϵ`$ dimensions and the one-loop bubble in $`D=42ϵ`$ dimensions. The $`ϵ`$ expansion of these integrals is given in the Appendix and can be inserted directly into the following expressions. We can write

$$𝒞^{8(1\times 1)}(s,t,u)=𝒫_𝒮(s,t,u)𝒯_2(s,t,u)s^2tu+2Re\left\{𝒫_{}(s,t,u)\right\}+_𝒯(s,t,u),$$
(64)
where the self-interference of the singular terms of the amplitude is

$`𝒫_𝒮(s,t,u)`$ $`=`$ $`{\displaystyle \frac{VN}{4}}\left|𝒫_𝒞(s,t,u)\right|^2`$
$`+{\displaystyle \frac{V}{4}}2Re\left\{𝒫_𝒞(s,t,u)^{}[𝒫_𝒜(s,t,u)+_𝒜(s,t,u)+(tu)]\right\}`$
$`+{\displaystyle \frac{V^2}{4N}}\left\{\right|𝒫_𝒜(s,t,u)+_𝒜(s,t,u)|^2+(tu)\}`$
$`{\displaystyle \frac{V}{4N}}\mathrm{\hspace{0.17em}2}Re\left\{\left[𝒫_𝒜(s,t,u)^{}+_𝒜(s,t,u)^{}\right]\left[𝒫_𝒜(s,u,t)+_𝒜(s,u,t)\right]\right\},`$
and the interference of the singular terms with the one-loop amplitude is

$`𝒫_{}(s,t,u)`$ $`=`$ $`{\displaystyle \frac{s}{2}}\{{\displaystyle \frac{VN}{4}}𝒫_𝒞(s,t,u)^{}u_3(s,t,u)`$ (66)
$`+{\displaystyle \frac{V}{4}}[(𝒫_𝒜(s,t,u)^{}+_𝒜(s,t,u)^{})u_3(s,t,u)`$
$`+𝒫_𝒞(s,t,u)^{}u_2(s,t,u)+(tu)]`$
$`+{\displaystyle \frac{V^2}{4N}}[(𝒫_𝒜(s,t,u)^{}+_𝒜(s,t,u)^{})u_2(s,t,u)+(tu)]`$
$`{\displaystyle \frac{V}{4N}}[(𝒫_𝒜(s,u,t)^{}+_𝒜(s,u,t)^{})u_2(s,t,u)+(tu)]\},`$
where the functions $`_2`$ and $`_3`$ are defined in Eqs. (35) and (36) respectively. In Eqs. (4) and (66)

$`𝒫_𝒜(s,t,u)`$ $`=`$ $`{\displaystyle \frac{12ϵ}{2st}}\{\beta _0[\mathrm{Bub}\left(s\right)+\mathrm{Bub}\left(t\right)]+[{\displaystyle \frac{V}{N}}({\displaystyle \frac{1}{ϵ}}+{\displaystyle \frac{3}{2}}){\displaystyle \frac{3N}{2}}]\mathrm{Bub}\left(s\right)`$ (67)
$`+N({\displaystyle \frac{2}{ϵ}}+{\displaystyle \frac{3}{2}})\mathrm{Bub}\left(t\right)\},`$
$`𝒫_𝒞(s,t,u)`$ $`=`$ $`{\displaystyle \frac{12ϵ}{ϵ}}\left[{\displaystyle \frac{1}{tu}}\mathrm{Bub}\left(s\right)+{\displaystyle \frac{1}{su}}\mathrm{Bub}\left(t\right)+{\displaystyle \frac{1}{st}}\mathrm{Bub}\left(u\right)\right]`$ (68)
are infrared divergent while

$`_𝒜(s,t,u)={\displaystyle \frac{1}{2st}}\beta _0\left\{{\displaystyle \frac{(2ϵ)(38ϵ+8ϵ^2)}{(1ϵ)(32ϵ)}}\mathrm{Bub}\left(s\right){\displaystyle \frac{2}{ϵ}}\right\},`$ (69)
originates from ultraviolet terms which are rendered finite after renormalisation in the $`\overline{\mathrm{MS}}`$ scheme.

Finally, the finite contribution is given by

$`_𝒯(s,t,u)`$ $`=`$ $`{\displaystyle \frac{VN}{4}}_𝒜(s,t,u)+{\displaystyle \frac{V}{2}}\mathrm{\hspace{0.17em}2}Re\left\{_{}(s,t,u)\right\}+{\displaystyle \frac{V^2}{4N}}_𝒟(s,t,u){\displaystyle \frac{V}{4N}}_𝒞(s,t,u)`$ (70)
$`+(tu),`$
where

$`_𝒜(s,t,u)`$ $`=`$ $`\left[{\displaystyle \frac{2u\left(8t^4+4s^4+12st^3+13s^2t^2+4s^3t\right)}{s^2t}}\right]\left|\mathrm{Box}^6(s,t)\right|^2`$ (71)
$`+\left[{\displaystyle \frac{4u\left(2s^2+ts+4t^2\right)}{t}}\right]2Re\left\{\mathrm{Box}^6(s,t)^{}\mathrm{Box}^6(t,u)\right\}`$
$`+\left[{\displaystyle \frac{4tu\left(5s^2+4ts+4t^2\right)}{s^2}}\right]\mathrm{Box}^6(s,t)^{}\mathrm{Box}^6(s,u)`$
$`+\left[{\displaystyle \frac{4s^2\left(t^2+u^2\right)}{tu}}\right]\left|\mathrm{Box}^6(t,u)\right|^2,`$
$`_{}(s,t,u)`$ $`=`$ $`\beta _0_{}(t,s)^{}\{\left[{\displaystyle \frac{2\left(t^2+u^2\right)}{t}}\right]\mathrm{Box}^6(t,u)+\left[{\displaystyle \frac{u^2\left(2s^2+ts+4t^2\right)}{s^2t}}\right]\mathrm{Box}^6(s,t)`$ (72)
$`+\left[{\displaystyle \frac{t\left(4t^2+5s^2+7ts\right)}{s^2}}\right]\mathrm{Box}^6(s,u)\}`$
$`+\left[{\displaystyle \frac{(ut)(2s^4+3s^3t7sut^2+4t^4)N}{s^2t}}+{\displaystyle \frac{2s^4+2t^4+ts^34sut^2}{stN}}\right]\left|\mathrm{Box}^6(s,t)\right|^2`$
$`+\left[{\displaystyle \frac{(ut)t(4t^25su)N}{s^2}}+{\displaystyle \frac{t(2t^2+3ts+5s^2)}{sN}}\right]\mathrm{Box}^6(s,u)^{}\mathrm{Box}^6(s,t)`$
$`+\left[{\displaystyle \frac{2(ut)(2t^2su)N}{t}}+{\displaystyle \frac{2s(t^2+s^2)}{tN}}\right]\mathrm{Box}^6(t,u)^{}\mathrm{Box}^6(s,t)`$
$`+\left[(6t+s)N+{\displaystyle \frac{s(t+3s)}{tN}}\right]\mathrm{Box}^6(t,u)^{}_{}(t,s)`$
$`+\left[{\displaystyle \frac{(6u^2+s^2)N}{t}}+{\displaystyle \frac{2(5s^28tu)}{Nt}}\right]ϵ\mathrm{Box}^6(t,u)^{}\mathrm{Bub}\left(s\right)`$
$`+\left[{\displaystyle \frac{(12t^3+17t^2s+13ts^2+2s^3)N}{2s^2}}+{\displaystyle \frac{(t+3s)(2t^2+ts+2s^2)}{2stN}}\right]\mathrm{Box}^6(s,t)^{}_{}(t,s)`$
$`+\left[{\displaystyle \frac{(12t^2+11ts+9s^2)tN}{2s^2}}+{\displaystyle \frac{(t+3s)t}{sN}}\right]\mathrm{Box}^6(s,u)^{}_{}(t,s)`$
$`+\left[{\displaystyle \frac{(12t^2+33ts+23s^2)tN}{2s^2}}+{\displaystyle \frac{2(8t^2+14ts+11s^2)t}{s^2N}}\right]ϵ\mathrm{Box}^6(s,u)^{}\mathrm{Bub}\left(s\right)`$
$`+[{\displaystyle \frac{(12t^4+39t^3s+38t^2s^2+31ts^3+14s^4)N}{2s^2t}}`$
$`+{\displaystyle \frac{16t^4+36t^3s+34t^2s^2+21ts^3+10s^4}{s^2tN}}\left]ϵ\mathrm{Box}^6\right(s,t\left)^{}\mathrm{Bub}\right(s),`$
$`_𝒞(s,t,u)`$ $`=`$ $`2\beta _0^2\left\{\left[{\displaystyle \frac{t^2+u^2}{s^2}}\right]_{}(t,s)^{}_{}(u,s){\displaystyle \frac{2ut}{s^2}}ϵ^2\left|\mathrm{Bub}\left(s\right)\right|^2\right\}`$ (73)
$`+\beta _0\left\{[{\displaystyle \frac{2(t^2+8ut+u^2)N}{s^2}}+{\displaystyle \frac{2}{N}}]ϵ^2\right|\mathrm{Bub}\left(s\right)|^2`$
$`+\left[{\displaystyle \frac{(6t^2+s^2)N}{s^2}}+{\displaystyle \frac{2(5s^28tu)}{s^2N}}\right]2ϵRe\left[\mathrm{Bub}\left(s\right)^{}_{}(t,s)\right]`$
$`+\left[{\displaystyle \frac{2(ut)(2t^2su)N}{s^2}}+{\displaystyle \frac{2(t^2+s^2)}{sN}}\right]\mathrm{\hspace{0.17em}2}Re\left[\mathrm{Box}^6(s,t)^{}_{}(u,s)\right]`$
$`+[{\displaystyle \frac{(5s^212tu)N}{s^2}}+{\displaystyle \frac{5}{N}}]_{}(t,s)^{}_{}(u,s)\}`$
$`+\left[{\displaystyle \frac{4(s^2ut)(tu)^2N^2}{s^2}}+{\displaystyle \frac{4s^2}{N^2}}\right]\mathrm{Box}^6(s,t)^{}\mathrm{Box}^6(s,u)`$
$`+\left[{\displaystyle \frac{3(s^23tu)N^2}{s^2}}+{\displaystyle \frac{15}{2}}\right]_{}(t,s)^{}_{}(u,s)`$
$`+\left[{\displaystyle \frac{(ut)(6t^2+8ts+5s^2)N^2}{s^2}}+{\displaystyle \frac{5t^2+2ts+3s^2}{s}}+{\displaystyle \frac{2st}{N^2}}\right]\mathrm{\hspace{0.17em}2}Re\left\{\mathrm{Box}^6(s,t)^{}_{}(u,s)\right\}`$
$`+[{\displaystyle \frac{(tu)(6t^23ts+s^2)N^2}{s^2}}{\displaystyle \frac{32t^3+35t^2s+22ts^2+11s^3}{s^2}}`$
$`+{\displaystyle \frac{2(4t^2+5s^2)}{sN^2}}]2ϵRe\{\mathrm{Box}^6(s,t)^{}\mathrm{Bub}\left(s\right)\}`$
$`+\left[{\displaystyle \frac{3(2s^2+5ts6t^2)N^2}{2s^2}}+{\displaystyle \frac{48t^2+5ts3s^2}{2s^2}}+{\displaystyle \frac{4(t+3s)}{sN^2}}\right]2ϵRe\left\{_{}(t,s)^{}\mathrm{Bub}\left(s\right)\right\}`$
$`+\left[{\displaystyle \frac{(7s^236tu)N^2}{2s^2}}{\displaystyle \frac{43s^248tu}{s^2}}+{\displaystyle \frac{47s^264tu}{s^2N^2}}\right]ϵ^2\left|\mathrm{Bub}\left(s\right)\right|^2`$
and

$`_𝒟(s,t,u)`$ $`=`$ $`\beta _0^2\left\{{\displaystyle \frac{4ut}{s^2}}ϵ^2\left|\mathrm{Bub}\left(s\right)\right|^2+{\displaystyle \frac{2u(t^2+u^2)}{s^2t}}\left|_{}(t,s)\right|^2\right\}`$ (74)
$`+\beta _0\left\{[{\displaystyle \frac{2u(6t+s)N}{s^2}}+{\displaystyle \frac{2u(t+3s)}{stN}}]\right|_{}(t,s)|^2`$
$`+\left[{\displaystyle \frac{u(6t^2+12ts+7s^2)N}{s^2t}}+{\displaystyle \frac{2u(5s^28tu)}{s^2tN}}\right]2ϵRe\left[\mathrm{Bub}\left(s\right)^{}_{}(t,s)\right]`$
$`+\left[{\displaystyle \frac{2u(ut)(2t^2su)N}{s^2t}}+{\displaystyle \frac{2u(t^2+s^2)}{stN}}\right]\mathrm{\hspace{0.17em}2}Re\left[\mathrm{Box}^6(s,t)^{}_{}(t,s)\right]`$
$`+[{\displaystyle \frac{4u(s3t)N}{s^2}}+{\displaystyle \frac{4u}{sN}}]ϵ^2\left|\mathrm{Bub}\left(s\right)|^2\right\}`$
$`+[{\displaystyle \frac{(tu)(6t^4+21t^3s+25t^2s^2+19ts^3+7s^4)N^2}{s^2tu}}`$
$`{\displaystyle \frac{32t^5+99t^4s+130t^3s^2+100t^2s^3+58ts^4+17s^5}{s^2tu}}`$
$`+{\displaystyle \frac{2(4t^48t^2su+8ts^3+5s^4)}{stuN^2}}]2ϵRe\{\mathrm{Box}^6(s,t)^{}\mathrm{Bub}\left(s\right)\}`$
$`+[{\displaystyle \frac{2(tu)^2(2t^44t^2su+2ts^3+s^4)N^2}{s^2tu}}+{\displaystyle \frac{4(ut)(t^42t^2sus^3u)}{stu}}`$
$`+{\displaystyle \frac{2(s^4+t^2u^2)}{tuN^2}}\left]\right|\mathrm{Box}^6(s,t)|^2`$
$`+[{\displaystyle \frac{(ut)(6t^3+10t^2s+7ts^2+s^3)N^2}{s^2u}}+{\displaystyle \frac{t^43t^3s6t^2s^29ts^33s^4}{stu}}`$
$`+{\displaystyle \frac{s^2(t+3s)}{tuN^2}}]\mathrm{\hspace{0.17em}2}Re\{\mathrm{Box}^6(s,t)^{}_{}(t,s)\}`$
$`+[{\displaystyle \frac{(36t^4+114t^3s+191t^2s^2+160ts^3+51s^4)N^2}{2s^2tu}}`$
$`{\displaystyle \frac{2(24t^4+99t^3s+156t^2s^2+117ts^3+34s^4)}{s^2tu}}`$
$`+{\displaystyle \frac{64t^4+192t^3s+241t^2s^2+162ts^3+51s^4}{s^2tuN^2}}\left]ϵ^2\right|\mathrm{Bub}\left(s\right)|^2`$
$`+[{\displaystyle \frac{(18t^3+57t^2s+42ts^2+7s^3)N^2}{2s^2u}}+{\displaystyle \frac{48t^4+101t^3s+49t^2s^233ts^321s^4}{2s^2tu}}`$
$`+{\displaystyle \frac{(t+3s)(4u^2+s^2)}{stuN^2}}]2ϵRe\{\mathrm{Bub}\left(s\right)^{}_{}(t,s)\}`$
$`+\left[{\displaystyle \frac{t(9t^2+12ts+5s^2)N^2}{s^2u}}+{\displaystyle \frac{(t+3s)(s+3t)}{su}}+{\displaystyle \frac{(t+3s)^2}{2tuN^2}}\right]\left|_{}(t,s)\right|^2,`$
with

$$_{}(t,s)=\mathrm{Bub}\left(t\right)(1+ϵ)\mathrm{Bub}\left(s\right).$$
(75)
Only the first term in the expansions of the box integrals are required in Eqs. (71)–(74), where we have systematically discarded contributions of $`𝒪\left(ϵ\right)`$. The bubble integrals must be expanded through to $`𝒪\left(ϵ^0\right)`$.

In this section, we have isolated the infrared divergences of the one-loop amplitude into the terms $`𝒫_𝒜`$ and $`𝒫_𝒞`$. Both functions depend on one-loop bubble integrals and diverge as $`1/ϵ^2`$. The separation of the singular terms becomes evident when we complete our basis of one-loop master integrals with the finite box in $`D=62ϵ`$. In the context of one-loop integrals, this form for the divergences arises naturally. The singular behaviour of the one-loop amplitude can also be predicted applying the formalism of  which yields

$$𝒫_{𝒜}^{}{}_{,C}{}^{}(s,t,u)=\frac{1}{st}𝑨(ϵ,s,u,t)$$
(76)
and

$$𝒫_{𝒞}^{}{}_{,C}{}^{}(s,t,u)=\frac{1}{su}𝑩(ϵ,s,t,u)\frac{1}{st}𝑩(ϵ,s,u,t),$$
(77)
where $`𝑨(ϵ,s,t,u)`$ and $`𝑩(ϵ,s,t,u)`$ are defined in Eqs. (25) and (26). The singular structure of $`𝒫_{𝒜}^{}{}_{,C}{}^{}`$ and $`𝒫_{𝒞}^{}{}_{,C}{}^{}`$ around $`ϵ=0`$ precisely matches that given in Eq. (67) and (68) and the two formalisms differ only in the finite remainder. In order to make the agreement explicit one could replace $`𝒫_𝒜`$($`𝒫_𝒞`$) with $`𝒫_{𝒜}^{}{}_{,C}{}^{}`$($`𝒫_{𝒞}^{}{}_{,C}{}^{}`$) in Eqs. (4) and (66) with appropriate modifications due to the finite differences $`𝒫_𝒜𝒫_{𝒜}^{}{}_{,C}{}^{}`$ and $`𝒫_𝒞𝒫_{𝒞}^{}{}_{,C}{}^{}`$ in Eqs. (66) and (70).

## 5 Summary

In this paper we presented the $`𝒪\left(\alpha _s^4\right)`$ QCD corrections to the $`22`$ scattering processes $`q\overline{q}gg`$, $`gg\overline{q}q`$ and the associated crossed processes $`qgqg`$ and $`\overline{q}g\overline{q}g`$ in the high energy limit, where the quark masses can be ignored. We computed renormalised analytic expressions for the interference of the tree-level diagrams with the two-loop ones and for the self-interference of one-loop graphs in the $`\overline{\mathrm{MS}}`$ scheme. Throughout we employed conventional dimensional regularisation.

The renormalised matrix elements are infrared divergent and contain poles down to $`𝒪\left(1/ϵ^4\right)`$. The singularity structure of one- and two-loop diagrams has been thoroughly studied by Catani  who provided a procedure for predicting the infrared behaviour of renormalised amplitudes. The anticipated pole structure agrees exactly with that obtained by direct Feynman diagram evaluation. In fact Catani’s method does not determine the $`1/ϵ`$ poles exactly, but expects that the remaining unpredicted $`1/ϵ`$ poles are non-logarithmic and proportional to constants (colour factors, $`\pi ^2`$ and $`\zeta _3`$). We find that this is indeed the case, and the constant $`H^{(2)}`$ is given in Eq. (3.1). This provides a very strong check on the reliability of our results. Similarly, the infrared divergent structure of the squared one-loop diagrams we found by direct evaluation agrees with the expected pole structure.

The pole structure of the two-loop contribution is given in Eq. (24) while expressions for the finite parts in the $`s`$-channel and $`u`$-channels according to the colour decomposition of Eq. (3.2) are given in Secs. 3.2.1 and 3.2.2 respectively. Similarly, the infrared divergent one-loop contributions along with the remaining finite parts are detailed in Sec. 4. The expressions for the two-loop pole structure and the singular and finite parts of the self-interference of one-loop graphs are analytic and are given in terms of the one-loop bubble graph and one-loop box graph in $`D=62ϵ`$ dimensions. To evaluate these formulae in the appropriate physical region requires the insertion of the series expansion of the one-loop graphs around $`ϵ=0`$. These expansions are given in Appendix A.

In this paper, we have concentrated only on QCD processes. However, the QED processes $`e^+e^{}\gamma \gamma `$ and Compton scattering $`e^{}\gamma e^{}\gamma `$ as well as $`\gamma \gamma e^{}e^+`$ are also of interest. We note that the expressions given here are related to those for the corresponding QED processes in the limit where the electron mass can be ignored by the alteration of various colour factors. We expect to address this in a separate article.

The results presented here, together with those previously computed for quark-quark scattering  are necessary ingredients for the next-to-next-to-leading order predictions for jet cross sections in hadron-hadron collisions. On their own, they are insufficient to make physical predictions and much work remains to be done. First, at the matrix element level, similar expressions to those presented here for gluon-gluon scattering are needed. Given the recent progress in the field, we anticipate that this problem will soon be solved. Second, a systematic procedure for analytically canceling the infrared divergences between the tree-level $`24`$, the one-loop $`23`$ and the $`22`$ processes needs to be established for semi-inclusive jet cross sections. Again, recent progress in determining the singular limits of tree-level matrix elements when two particles are unresolved  and the soft and collinear limits of one-loop amplitudes , together with the analytic cancellation of the infrared singularities in the somewhat simpler case of $`e^+e^{}\mathrm{photon}+\mathrm{jet}`$ at next-to-leading order , suggest that the technical problems will soon be solved for generic $`22`$ scattering processes. There are additional problems due to initial state radiation. However, the recent steps taken towards the determination of the three-loop splitting functions  are also promising.

Third, a numerical implementation of the various contributions must be developed. The next-to-leading order programs for three jet production that have already been written provide a first step in this direction . We are confident that the problem of the numerical cancellation of residual infrared divergences will soon be addressed thereby enabling the construction of numerical programs to provide next-to-next-to-leading order QCD estimates of jet production in hadron collisions.

## Acknowledgements

C.A. acknowledges the financial support of the Greek Government and M.E.T. acknowledges financial support from CONACyT and the CVCP. We gratefully acknowledge the support of the British Council and German Academic Exchange Service under ARC project 1050. This work was supported in part by the EU Fourth Framework Programme ‘Training and Mobility of Researchers’, Network ‘Quantum Chromodynamics and the Deep Structure of Elementary Particles’, contract FMRX-CT98-0194 (DG-12-MIHT), in part by the U.S. Department of Energy under Grant No. DE-FG02-95ER40896 and in part by the University of Wisconsin Research Committee with funds granted by the Wisconsin Alumni Research Foundation.

## Appendix A One-loop master integrals

In this appendix, we list the expansions for the one-loop box integrals in $`D=62ϵ`$. We remain in the physical region $`s>0`$, $`u,t<0`$, and write coefficients in terms of logarithms and polylogarithms that are real in this domain. More precisely, we use the notation of Eqs. (43) and (44) to define the arguments of the logarithms and polylogarithms. The polylogarithms are defined as in Eq. (42).

We find that the box integrals have the expansion,

$`\mathrm{Box}^6(u,t)`$ $`=`$ $`{\displaystyle \frac{e^{ϵ\gamma }\mathrm{\Gamma }\left(1+ϵ\right)\mathrm{\Gamma }\left(1ϵ\right)^2}{2s\mathrm{\Gamma }\left(12ϵ\right)\left(12ϵ\right)}}\left({\displaystyle \frac{\mu ^2}{s}}\right)^ϵ\{{\displaystyle \frac{1}{2}}[(XY)^2+\pi ^2]`$ (78)
$`+2ϵ\left[\mathrm{Li}_3(x)X\mathrm{Li}_2(x){\displaystyle \frac{1}{3}}X^3{\displaystyle \frac{\pi ^2}{2}}X\right]`$
$`2ϵ^2[\mathrm{Li}_4(x)+Y\mathrm{Li}_3(x){\displaystyle \frac{1}{2}}X^2\mathrm{Li}_2(x){\displaystyle \frac{1}{8}}X^4{\displaystyle \frac{1}{6}}X^3Y+{\displaystyle \frac{1}{4}}X^2Y^2`$
$`{\displaystyle \frac{\pi ^2}{4}}X^2{\displaystyle \frac{\pi ^2}{3}}XY{\displaystyle \frac{\pi ^4}{45}}]+(ut)\}+𝒪(ϵ^3),`$
and

$`\mathrm{Box}^6(s,t)`$ $`=`$ $`{\displaystyle \frac{e^{ϵ\gamma }\mathrm{\Gamma }(1+ϵ)\mathrm{\Gamma }(1ϵ)^2}{2u\mathrm{\Gamma }(12ϵ)(12ϵ)}}({\displaystyle \frac{\mu ^2}{u}})^ϵ\{(X^2+2i\pi X)`$ (79)
$`+ϵ[(2\mathrm{L}\mathrm{i}_3(x)+2X\mathrm{Li}_2(x){\displaystyle \frac{2}{3}}X^3+2YX^2\pi ^2X+2\zeta _3)`$
$`+i\pi (2\mathrm{L}\mathrm{i}_2(x)+4YXX^2{\displaystyle \frac{\pi ^2}{3}})]`$
$`+ϵ^2[(2\mathrm{L}\mathrm{i}_4(z)+2\mathrm{L}\mathrm{i}_4(y)2Y\mathrm{Li}_3(x)2X\mathrm{Li}_3(y)+(2XYX^2\pi ^2)\mathrm{Li}_2(x)`$
$`+{\displaystyle \frac{1}{3}}X^4{\displaystyle \frac{5}{3}}X^3Y+{\displaystyle \frac{3}{2}}X^2Y^2+{\displaystyle \frac{2}{3}}\pi ^2X^22\pi ^2XY+2Y\zeta _3+{\displaystyle \frac{1}{6}}\pi ^4)`$
$`+i\pi (2\mathrm{L}\mathrm{i}_3(x)2\mathrm{L}\mathrm{i}_3(y)+2Y\mathrm{Li}_2(x)+{\displaystyle \frac{1}{3}}X^32X^2Y+3XY^2`$
$`{\displaystyle \frac{\pi ^2}{3}}Y+2\zeta _3)]\}+𝒪(ϵ^3).`$
$`\mathrm{Box}^6(s,u)`$ is obtained from Eq. (79) by exchanging $`u`$ and $`t`$.

Finally, the one-loop bubble integral in $`D=42ϵ`$ dimensions is given by

$$\mathrm{Bub}(s)=\frac{e^{ϵ\gamma }\mathrm{\Gamma }\left(1+ϵ\right)\mathrm{\Gamma }\left(1ϵ\right)^2}{\mathrm{\Gamma }\left(22ϵ\right)ϵ}\left(\frac{\mu ^2}{s}\right)^ϵ.$$
(80)