---
author:
- 'Taushif Ahmed,'
- 'Maguni Mahakhud,'
- 'Prakash Mathews,'
- Narayan Rana
- 'and V. Ravindran'
title: 'Two-Loop QCD Corrections to Higgs $\rightarrow b + \bar{b} + g$ Amplitude'
---

Introduction
============

\[sec:intro\] The tests of the Standard Model (SM) have been going on for several decades in various experiments and most of its predictions have been tested in an unprecedented accuracy. The recent discovery of Higgs boson by ATLAS [@atlashiggs] and CMS [@cmshiggs] collaborations at the Large Hadron Collider (LHC) puts the SM on firm footing. The Higgs boson results from Higgs mechanism that provides a framework for electroweak symmetry breaking. Elementary particles such as leptons, quarks, gauge bosons and Higgs boson acquire masses through the Higgs mechanism. The mass of the Higgs boson being a parameter of the theory can not be predicted by the SM and hence its discovery provides a valuable information on this. Results from Higgs searches at LEP [@higgslep] and Tevatron [@higgstev] were crucial ingredients to the recent discovery in narrowing down the search regions for the LHC collaborations. The direct searches at the LEP excluded Higgs of mass below 114.4 GeV and the precision electroweak measurements [@lepprecis] hinted for Higgs boson in the mass less than 152 GeV at $95\%$ confidence level (CL). Tevatron on the other hand excluded Higgs of mass in the range $162-166$ GeV at $95\%$ CL.

The dominant production mechanism for the Higgs production at the LHC is gluon gluon fusion through top quark loop. The subdominant ones come from vector boson fusion, associated production of Higgs with vector bosons and top anti-top pairs and bottom anti-bottom annihilation. The inclusive production cross section for the Higgs production is known to an unprecedented accuracy due to many breakthroughs in the computation of amplitudes, loop and phase space integrals. For gluon-gluon [@gghNNLO], vector boson fusion processes [@bolzoni], and associated production with vector bosons [@Han:1991ia], the inclusive rates are known to NNLO accuracy in QCD. There are also studies related to the Higgs production in association with bottom quarks which were also motivated to study Higgs boson in certain SUSY models, namely MSSM. The coupling of bottom quarks become large in the large $\tan\beta$ region, where $\tan\beta$ is the vacuum expectation values of up and down type Higgs fields in the Higgs sector of MSSM. Such large couplings can enhance gluon fusion as well as bottom quark fusion subprocesses. Fully inclusive cross section for Higgs production in association with bottom quark to NNLO level accuracy is also known in the variable flavour scheme (VFS) [@vfs], while it is known only up to NLO level in the fixed flavour scheme (FFS) [@ffs]. In the VFS, one assumes the initial state bottom quarks inside the proton. They are there as a result of emission of collinear bottom anti-bottom states from the gluons intrinsically present inside the proton. They being collinear give large logs which need to be resummed. The resummed contribution is the source for non-vanishing bottom and anti-bottom parton distribution functions inside the proton in the VFS scheme.

The differential distributions for Higgs production and its decay to pair of photons [@Anastasiou:2005qj] or massive vector bosons [@Anastasiou:2007mz; @Grazzini:2008tf] have also been known at NNLO level in QCD in the infinite top quark mass limit. Such exclusive observables allow direct comparison of theoretical predictions with experimental results which include kinematical cuts on the final state particles. In particular, observables with jet vetos enhance the significance of the signal considerably allowing us to study the properties of Higgs boson and its coupling to other SM particles. NNLO QCD prediction [@Boughezal:2013uia] for production of Higgs with one jet through effective gluon-gluon-higgs vertex in the infinite top quark mass limit is available, thanks to various ingredients that are computed to the required accuracy by different groups [@Hggg; @Gehrmann:2011aa]. As the experimental accuracy improves, it will be important to include other subdominant production mechanisms. In this article, we provide the relevant one and two loop amplitudes for the process $H \rightarrow b + \overline b + g $ which is analytically continued also to obtain the production of Higgs boson with one jet in bottom anti-bottom annihilation, i.e., $b+\overline b \rightarrow H+g$, where Higgs couples to bottom quark through Yukawa coupling denoted by $\lambda$. We use VFS scheme throughout. This will be an important supplement to the Higgs boson with one jet at NNLO level as it includes the bottom quark effects in VFS scheme.

Beyond leading order in perturbation theory, one encounters large number of Feynman amplitudes with rich Lorentz and gauge structures. In addition, the loop integrals become increasingly complicated due to their multiple kinematic dependence. Generation of diagrams, simplification of Lorentz, Dirac and color indices can be done symbolically. Using integration by parts (IBP) and Lorentz invariant (LI) identities the large number of loop integrals can be reduced in a rather straight forward way to few master integrals (MI). The two loop MIs for four legs processes where all fields but one external leg are massless were solved by Gehrmann and Remiddi [@Gehrmann:2000zt] using an elegant method of differential equations.

In this article we present one and two loop QCD amplitudes for the process $H \rightarrow b + \overline b + g $ treating both bottom and other four light quarks massless. We do not include top quark in our analysis. To obtain infrared safe observables, we require, in addition to these two loop amplitudes, one loop corrected $H\rightarrow b+\overline b + 2~{\rm partons }$ and tree level $H \rightarrow b+\overline b + 3~{\rm partons}$ amplitudes. Note that they are individually infrared singular due to the presence of massless partons in the amplitudes. There exist several equally efficient frameworks which use these infrared sensitive contributions to combine them to obtain infrared safe observables. They go by the names sector decomposition [@secdec], $q_T$-subtraction [@qtsub] and antenna subtraction [@antsub] methods. More recently the method developed by Czakon using sector decomposition and FKS [@fks] phase space slicing, was applied to obtain top quark pair production [@topnnlo] at NNLO level and NNLO QED corrections [@nnloqed] to $Z\rightarrow e^+ e^-$. Antenna subtraction was used to obtain NNLO QCD corrections to di-jet production at the LHC. The NNLO corrections to Higgs plus one jet resulting from only gluon-gluon-Higgs effective interaction are obtained recently in [@Gehrmann:2011aa] making best use of the subtraction methods in an efficient way. The amplitudes presented in this article will constitute contributions coming from bottom-antibottom-higgs interactions to Higgs plus one jet observable at NNLO level. We have presented the amplitudes in the form suitable for easier implementation to study infra-red safe hadron level observables involving Higgs plus one jet at NNLO in QCD.

In the next section, we discuss the Lagrangian that describes coupling of Higgs boson with bottom quark, explain how the projector technique can be used to obtain the amplitudes and describe the renormalization and factorization properties of the amplitudes. Section \[sec:calc\] is dedicated to the computational details. Final results in compact form are given in Section \[sec:result\] and corresponding coefficients are given in the Appendix. In section \[sec:conc\], we conclude with our findings.

Theory
======

The interaction part of the action involving bottom quarks and Higgs boson is given by $$\begin{aligned}
S^b_{I} = - \lambda \int d^4 x \, \phi(x) \overline \psi_b(x) \psi_b(x)\end{aligned}$$ where, $\psi_b(x)$ denotes the bottom quark field and $\phi(x)$ the scalar field. $\lambda$ is the Yukawa coupling given by $\sqrt{2} m_b/\upsilon$, with the bottom quark mass $m_b$ and the vacuum expectation value $\upsilon \approx 246$ GeV. For the pseudoscalar Higgs of MSSM, we need to replace $\lambda \phi(x) \overline \psi_b(x) \psi_b(x)$ by $\tilde \lambda \tilde \phi(x) \overline \psi_b(x) \gamma_5 \psi_b(x)$ in the above equation. The MSSM couplings are $$\tilde{\lambda} = \left\{
  \begin{array}{ll}
    -  \frac{\sqrt{2} m_b \sin\alpha}{\upsilon \cos\beta}  \,,& \qquad \tilde{\phi} = h\,,\\
    \phantom{-}  \frac{\sqrt{2} m_b \cos\alpha}{\upsilon \cos\beta}  \,,& \qquad \tilde{\phi}=H\,,\\
    \phantom{-}  \frac{\sqrt{2} m_b \tan\beta}{\upsilon } \,, & \qquad \tilde{\phi}=A\,
  \end{array}
  \right.$$ respectively. The angle $\alpha$ is the measure of mixing of weak and mass eigenstates of neutral Higgs bosons. In the VFS scheme, except in the Yukawa coupling, $m_b$ is taken to be zero like other light quarks in the theory. The number of active flavours is taken to be $n_f=5$. We work in Feynman gauge throughout.

Notation and kinematics {#subsec:defs}
-----------------------

We consider the decay of Higgs boson to a bottom quark, anti-bottom quark and a gluon $$H(q) \longrightarrow b(p_1) + \bar{b} (p_2) + g(p_3) \, .$$ The associated Mandelstam variables are defined as $$s \equiv (p_1 + p_2)^2, \hspace{1cm} t \equiv (p_2 + p_3)^2, \hspace{1cm} u \equiv (p_1 + p_3)^2$$ which satisfy $$s >0,~ t > 0, ~ u > 0, ~~ s + t + u = M_H^2 \equiv Q^2 > 0$$ where, $M_H$ is the mass of the Higgs boson. We also define the following dimensionless invariants which appear in harmonic polylogarithms (HPL) [@remiddi] and 2dHPL [@Gehrmann:2000zt] as $$x \equiv s/Q^2, \hspace{1cm} y \equiv u/Q^2, \hspace{1cm} z \equiv t/Q^2$$ satisfying $$0 < x < 1,~ 0 < y < 1,~ 0 < z < 1,~ ~\text{and}~ x + y + z = 1.$$

### Analytical continuation {#sec:ancont .unnumbered}

In order to compute the Higgs + 1 jet production at hadron colliders, the decay amplitudes must be analytically continued to the appropriate kinematical regions. The corresponding processes are $$\begin{aligned}
  \nonumber
  1. ~~~~ \overline{b} (-p_1) + b(-p_2) &\rightarrow g(p_3) + H (p_4) 
  \\ \nonumber
  2. ~~~~ b (-p_2) + g(-p_3) &\rightarrow b (p_1) + H (p_4)
  \\ \label{subprocess}
  3. ~~~~ \overline{b}(-p_1) + g(-p_3) &\rightarrow \overline{b} (p_2) + H (p_4)\end{aligned}$$

For the process 1, $Q^2 = M_H^2 > 0$, $s > 0 ,~ t < 0 $ and $u < 0$. Hence we introduce the dimensionless parameters $u_1$ and $v_1 $ with the following definitions $$u_1 \equiv - \frac{u}{s}, \hspace{1cm} v_1 \equiv \frac{Q^2}{s}$$ such that $0 < u_1 < 1$ and $0< v_1 < 1$.

Similarly, for the process 2, $Q^2 = M_H^2 > 0$, $s < 0 ,~ t > 0 $ and $u < 0$ and the dimensionless parameters are $u_2$ and $v_2 $ with the following definitions $$u_2 \equiv - \frac{u}{t}, \hspace{1cm} v_2 \equiv \frac{Q^2}{t}$$ such that $0 < u_2 < 1$ and $0< v_2 < 1$. The last one is trivially related to the second one.

The general structure of the amplitude {#sec:gentensor}
--------------------------------------

In this section, we describe how the amplitude for $H \rightarrow b+\overline b + g$ can be obtained using projector technique. Since the amplitude contains one external gluon, it can be expressed as $$| \M \rangle = {\cal S}_{\mu} (b,\bar{b};g) \varepsilon^{\mu}$$ where, $\varepsilon^{\mu}$ is the gluon polarization vector.

We observe the amplitude has the following general structure in terms of the coefficients $A',A''$ and $A_2$: $${\cal S}_{\mu} (b,\bar{b};g) = \bar{u} (p_1) \Big\{
A'~ p_{1 \mu} + A''~ p_{2 \mu} + A_2~\slashed p_3 \gamma_{\mu} 
\Big\} v (p_2)$$ where, we have used $p_3.\varepsilon = 0$. QCD Ward identity gives $$A'~ p_{1}. p_3 + A''~ p_{2} . p_3 = 0 
 ~~ \Rightarrow ~~
 A' = - A''~ \frac{p_{2} . p_3}{p_{1}. p_3} \equiv A_1~ p_2.p_3 \, .$$ Hence, the amplitude takes the following form: $$\begin{aligned}
{\cal S}_{\mu} (b,\bar{b};g) ~\varepsilon^{\mu} &= \bar{u} (p_1) \Big\{
A_1~( p_2.p_3 ~ p_{1 \mu} - p_1 . p_3 ~ p_{2 \mu} ) 
+ A_2~\slashed p_3 \gamma_{\mu} 
\Big\} v(p_2) ~\varepsilon^{\mu}
\nonumber\\ \label{tencoeff}
&\equiv A_1~ {\rm T_1} + A_2~ {\rm T_2} \; .\end{aligned}$$ The coefficients $A_m ~ (m = 1, 2)$ can be obtained from the amplitude $| \M \rangle$ using appropriate projectors ${\cal P}(A_m)$ $$A_m=\sum_{\rm spins} {\cal P} (A_m) {\cal S}_{\mu} (b,\bar{b};g) ~\varepsilon^{\mu}  =\sum_{\rm spins} {\cal P} (A_m) | \M \rangle$$ where, in $d$ space-time dimensions, the projectors are found to be $$\begin{aligned}
 \label{projectors}
{\cal P} (A_1) &= \frac{2 (d-2)}{s^2 \, t \, u \, (d-3)} {\rm T_1}^{\dagger}
                + \frac{1}{s \, t \, u \, (d-3)} {\rm T_2}^{\dagger} \, ,
\nonumber\\
{\cal P} (A_2) &= \frac{1}{s \, t \, u \, (d-3)} {\rm T_1}^{\dagger}  
                + \frac{1}{2 \, t \, u \, (d-3)} {\rm T_2}^{\dagger} \, .       \end{aligned}$$ Expanding the coefficients $A_m$ in powers of strong coupling constant $a_s = g_s^2/16 \pi^2$, we obtain $$\label{coeffexpnd}
A_m = \frac{\lambda}{\mu_R^\epsilon} ~4\pi \sqrt{a_s} T^a_{ij} \Big\{
A_m^{(0)} + a_s A_m^{(1)}
+ a_s^2 A_m^{(2)} + {\cal O}(a_s^3)
\Big\}$$ where, $T^a$ are the Gell-Mann matrices, $a$ is adjoint and $i$, $j$ are fundamental indices of SU(3) and $\mu_R$ is the renormalization scale. These coefficients $A_m^{(l)}$ completely specify the amplitude order by order in perturbation theory.

As described in section \[sec:ancont\], for Higgs + 1 jet production, the above amplitudes have to be suitably crossed and the coefficients $A_{m}$ will be expressed in terms of corresponding $u_i$ and $v_i$.

Ultraviolet renormalization
---------------------------

The Feynman amplitudes for the process $H\rightarrow b+\overline b+g$ beyond leading order develop ultraviolet divergences in QCD. We have used dimensional regularization to regulate them taking space-time dimension to be $d=4+\epsilon$. The scale $\mu_0$ is introduced to scale the mass dimension of the dimension-full strong coupling constant in $d$ dimensions. If we denote the dimensionless strong coupling constant by $\hat g_s$ in $d$ dimensions, then the unrenormalized amplitude can be expanded in terms of $\hat a_s=\hat g_s^2/16 \pi^2$ as $$\label{unm}
 |{\cal M} \rangle = \frac{\hat{\lambda}}{\mu_0^\epsilon} S_\epsilon \unas^{\frac{1}{2}} \left\{  |\unM^{(0)} \rangle + \unas |\unM^{(1)} \rangle + \unas^2 |\unM^{(2)} \rangle + {\cal O}(\hat{a}_s^3) \right\}$$ where, $S_{\ep} = \exp[\frac{\ep}{2} (\gamma_E - \ln 4\pi)]$ with Euler constant $\gamma_E = 0.5772 \ldots$ , results from loop integrals beyond leading order. $|\unM^{(i)} \rangle$ is the unrenormalized color-space vector which represents the $i^{th}$ loop amplitude. In $\overline {MS}$ scheme, the renormalized coupling constant $a_s \equiv a_s (\mu_R^2)$ at the renormalization scale $\mu_R$ is related to unrenormalized coupling constant $\hat{a}_s$ by $$\begin{aligned}
 \label{renas}
 \frac{\hat{a}_s}{\mu_0^{\epsilon}} S_{\epsilon} &= \frac{a_s}{\mu_R^{\epsilon}} {Z}(\mu_R^2)
\nonumber\\[1ex]
 &= \frac{a_s}{\mu_R^{\epsilon}} \left[   1 + a_s \left( \frac{1}{\ep} r_{a_{1;1}} \right) 
  + a_s^2 \left( \frac{1}{\ep^2} r_{a_{2;2}} + \frac{1}{\ep} r_{a_{2;1}} \right) + {\cal O}(a_s^3)  \right] \end{aligned}$$ where, $$r_{a_{1;1}} = 2 \b0 \;, \quad
 r_{a_{2;2}} = 4\b0^2 \;, \quad
 r_{a_{2;1}} = \beta_1 \;, \quad$$

$$\b0 = \left(\frac{11}{3} C_A - \frac{4}{3} T_F n_f\right) , \hspace{0.5cm}
 \beta_1 =  \left(\frac{34}{3} C_A^2 - \frac{20}{3} C_A T_F n_f - 4 C_F T_F n_f\right)$$

with $C_A = N$, $C_F = (N^2 -1)/{2 N} $, $T_F = 1/2$ and $n_f$ is the number of active quark flavors. The bare coupling constant $\hat \lambda$ is renormalized using $$\begin{aligned}
\label{renl}
\frac{\hat{\lambda}}{\mu_0^\epsilon} S_\epsilon &= \frac{\lambda}{\mu_R^\epsilon} {Z}_{\lambda}(\mu_R^2)
\nonumber\\[1ex]
 &= \frac {\lambda}{\mu_R^\epsilon} \left[   1 + a_s \left( \frac{1}{\ep} r_{\lambda_{1;1}} \right)
 + a_s^2 \left( \frac{1}{\ep^2} r_{\lambda_{2;2}}
             + \frac{1}{\ep} r_{\lambda_{2;1}} \right) 
 + {\cal O}(a_s^3)  \right] \; ,\end{aligned}$$ with $\lambda=\lambda(\mu_R^2)$ and $$r_{\lambda_{1;1}} = 6 C_F \,, ~  
r_{\lambda_{2;2}} = \Big(18 C_F^2 + 6 \b0 C_F\Big)\, , ~  
r_{\lambda_{2,1}} =  \left( \frac{3}{2} C_F^2
                    + \frac{97}{6} C_F C_A
                    - \frac{10}{3} C_F T_F n_f \right).$$ Using the eqn.(\[renas\]) and eqn.(\[renl\]), we now can express $|{\cal M } \rangle$ (eqn.(\[unm\])) in powers of renormalized $a_s$ with UV finite matrix elements $|\rnM^{(i)} \rangle$ $$\label{renamp}
 |\rnM \rangle = \frac{\lambda}{\mu_R^\epsilon} \rnas^{\frac{1}{2}} \Bigg( |\rnM^{(0)} \rangle + a_s |\rnM^{(1)} \rangle + a_s^{2} |\rnM^{(2)} \rangle + {\cal O}({a}_s^3) \Bigg)$$ where, $$\begin{aligned}
 \label{rln}
 |\rnM^{(0)} \rangle &= \left(  \frac{1}{\mu^{\ep}_R} \right)^{\frac{1}{2}}    |\unM^{(0)} \rangle \; ,
\nonumber\\[1ex]
|\rnM^{(1)} \rangle &= \left(  \frac{1}{\mu^{\ep}_R} \right)^{\frac{3}{2}} 
\left[ ~ |\unM^{(1)} \rangle + \mu^{\ep}_R \Big( \frac{r_{a_1}}{2} + r_{\lambda_1} \Big) |\unM^{(0)} \rangle ~ \right] \; ,
\nonumber\\[1ex]
|\rnM^{(2)} \rangle &= \left( \frac{1}{\mu^{\ep}_R} \right)^{\frac{5}{2}} 
\Big[ ~ |\unM^{(2)} \rangle  
+ \mu^{\ep}_R \Big( \frac{3r_{a_1}}{2} + r_{\lambda_1} \Big)  |\unM^{(1)} \rangle 
\nonumber\\
& \qquad \qquad \quad
+ \mu^{2\ep}_R \left( \frac{r_{a_2}}{2} - \frac{r_{a_1}^2}{8} + \frac{r_{a_1}}{2} r_{\lambda_1}
+ r_{\lambda_2}  \right)  |\unM^{(0)} \rangle ~ \Big]\end{aligned}$$ with $$\begin{aligned}
 r_{a_1} &=  \left(\frac{1}{\ep} r_{a_{1;1}}\right)  \;, \quad \quad
 r_{a_2} =  \left(\frac{1}{\ep^2} r_{a_{2;2}} + \frac{1}{\ep} r_{a_{2;1}}\right)  \;,
\nonumber\\
 r_{\lambda_1} &=  \left(\frac{1}{\ep} r_{\lambda_{1;1}}\right)  \;, \quad \quad
 r_{\lambda_2} =  \left(\frac{1}{\ep^2} r_{\lambda_{2;2}} + \frac{1}{\ep} r_{\lambda_{2;1}} \right) \,.\end{aligned}$$ We describe the computation of unrenormalized amplitudes $|\unM^{(l)} \rangle, l=0,1,2$ in section \[sec:calc\].

Infrared factorization {#sec:infrared}
----------------------

In addition to UV divergences, the amplitude suffers from soft and collinear divergences beyond leading order due to the presence of soft gluons and collinear massless partons in the loops. According to KLN theorem [@Kinoshita:1962ur; @Lee:1964is], to obtain infrared safe observables, we need to include appropriate contributions coming from real emission processes along with mass factorization counter terms and to perform sum over degenerate configurations. Thanks to factorization properties of QCD amplitudes, the infrared divergence structure of the amplitudes is well understood. The earliest account on two loop QCD amplitudes was by Catani [@catani1], who predicted the infrared poles in $\epsilon$ of multi-parton QCD amplitudes in dimensional regularization excluding two loop single pole. In [@sterman], Sterman and Tejeda-Yeomans demonstrated the connection of single pole in $\epsilon$ to a soft anomalous dimension matrix, later computed in [@Aybat:2006wq; @Aybat:2006mz] using factorization properties of the scattering amplitudes along with infrared evolution equations. The decomposition of single pole term into universal collinear and soft anomalous dimensions at two loop level in QCD was first observed in electromagnetic and Higgs form factors [@Ravindran:2004mb]. Becher and Neubert [@Becher:2009cu], using soft collinear effective theory, derived the exact formula for the infra-red divergences of scattering amplitudes with an arbitrary number of loops and legs in massless QCD including single pole in dimensional regularization. Gardi and Magnea also arrived at, a similar all order result [@Gardi:2009qi] using Wilson lines for hard partons and soft and eikonal jet functions in dimensional regularization. Following Catani, we express the renormalized amplitudes $|{\cal M}^{(i)}\rangle$ in terms of the universal subtraction operators ${\bf I}^{(i)}_b(\epsilon)$ as follows[^1]

$$\begin{aligned}
\label{catani1}
 |\rnM^{(1)}  \rangle &=& 2 \hspace{0.1cm} {\bf{I}}_{b}^{(1)} (\ep) \hspace{0.1cm} |\rnM^{(0)} \rangle
+ |\rnM^{(1)fin}  \rangle \, ,
\nonumber\\[2ex] 
|\rnM^{(2)} \rangle &=& 2 \hspace{0.1cm} {\bf{I}}_{b}^{(1)} (\ep) \hspace{0.1cm} |\rnM^{(1)} \rangle
+ 4 \hspace{0.1cm} {\bf{I}}^{(2)}_{b} (\ep) \hspace{0.1cm} |\rnM^{(0)}  \rangle
+ |\rnM^{(2)fin} \rangle \end{aligned}$$

$$\begin{aligned}
\text{where,} ~~~~
 {\bf{I}}^{(1)}_{b} (\ep) &= \frac{1}{2}   \frac{e^{- \frac{\ep}{2} \gamma_E}}{\Gamma(1+\frac{\ep}{2})}
 \Bigg\{
 \Big( \frac{4}{\ep^2} - \frac{3}{\ep} \Big) (C_A - 2 C_F) \Big[ \Big( - \frac{s}{\mu^2_R} \Big)^{\frac{\ep}{2}} \Big]
\nonumber\\[1ex] 
& +
 \Big( - \frac{4 C_A}{\ep^2} + \frac{3 C_A}{ 2 \ep} + \frac{\b0}{2 \ep} \Big)
\left[ \Big( - \frac{t}{\mu^2_R} \Big)^{\frac{\ep}{2}} + \Big( - \frac{u}{\mu^2_R} \Big)^{\frac{\ep}{2}}   \right]
 \Bigg\} \; ,
\nonumber\\[1ex]
{\bf{I}}^{(2)}_{b} (\ep) &= - \hspace{0.1cm} \frac{1}{2} {\bf{I}}^{(1)}_{b} (\ep) \Big[ {\bf{I}}^{(1)}_{b} (\ep) - \frac{2 \b0}{\ep}  \Big]
 +\hspace{0.1cm} \frac{e^{\frac{\ep}{2} \gamma_E} \hspace{0.1cm} \Gamma(1+\ep)}{\Gamma(1+\frac{\ep}{2})} \Big[ -\frac{\b0}{\ep} + K \Big] \hspace{0.1cm} {\bf{I}}^{(1)}_{b} (2\ep) \hspace{1.4cm}
\nonumber\\[1ex]
& +\hspace{0.1cm} \left( 2 {\bf{H}}^{(2)}_{q} (\ep) + {\bf{H}}^{(2)}_{g} (\ep) \right)\end{aligned}$$

with $$K = \left(\frac{67}{18}-\frac{\pi^2}{6} \right) C_A - \frac{10}{9} T_F n_f \, ,$$

$$\begin{aligned}
 {\bf{H}}^{(2)}_q (\ep) &= \frac{1}{\ep}
 \Bigg\{
        C_A C_F \Big( - \frac{245}{432} + \frac{23}{16} \zeta_2 - \frac{13}{4} \zeta_3 \Big)
      + C_F^2 \Big( \frac{3}{16} - \frac{3}{2} \zeta_2 + 3 \zeta_3 \Big)
\nonumber\\[1ex]
& \quad
      + C_F n_f \Big( \frac{25}{216} - \frac{1}{8} \zeta_2 \Big) 
 \Bigg\} \; ,
\nonumber\\[1ex]
 {\bf{H}}^{(2)}_g (\ep) &= \frac{1}{\ep} 
\Bigg\{ 
        C_A^2 \Big( - \frac{5}{24} - \frac{11}{48} \zeta_2 - \frac{1}{4} \zeta_3 \Big) 
      + C_A n_f \Big(\frac{29}{54} +  \frac{1}{24} \zeta_2 \Big) - \frac{1}{4} C_F n_f -\frac{5}{54} n_f^2  
\Bigg\} \, .\end{aligned}$$

The born amplitude $|\rnM^{(0)}\rangle$ and the finite parts $|\rnM^{(l)fin}\rangle, l=1,2$ are process dependent and hence they are determined by explicit computation.

Calculation of the amplitudes {#sec:calc}
=============================

![Planar topologies of master integrals[]{data-label="fig:planar"}](planar.jpg){width="80.00000%"}

![Non-planar topologies of master integrals[]{data-label="fig:nonplanar"}](nplanar.jpg){width="88.00000%"}

We now describe how we compute the coefficients $A_m$ from the amplitudes $|\hat {\cal M}^{(l)}\rangle$ for the process $H \rightarrow b+\overline b+g$ up to two loop level in QCD perturbation theory. QGRAF [@qgraf] is used to generate the Feynman amplitudes for this process. There are 2 diagrams at tree level, 13 at one loop and 251 at two loops excluding tadpole and self energy corrections to the external legs.

Using FORM [@Vermaseren:2000nd] and Mathematica, output of the QGRAF is converted to a form suitable for further symbolic manipulation. Using the projectors given in eqn.(\[projectors\]), we have projected out unrenormalized $\hat{A}_i$ from these amplitudes. They contain only scalar products among internal and external momenta. For the external on-shell gluon leg the physical polarization sum is done using $$\sum_s \varepsilon^{\mu}(p_3,s) \varepsilon^{\nu *}(p_3,s) = -g^{\mu\nu} + \frac{p_3^{\mu} q^{\nu} + q^{\mu} p_3^{\nu}}{p_3.q}$$ where, $p_3$ is the gluon momentum and $q$ is an arbitrary light-like 4-vector for which we choose $q = p_1$. The Lorentz contractions and Dirac algebra are done in $d = 4 +\epsilon$ dimensions. The next step involves the evaluation of one and two loop tensor and scalar integrals. This is done by first reducing them to an irreducible set of MIs using IBP identities and LI identities and substituting the MIs evaluated to desired accuracy in $\epsilon$. We have used a Mathematica package LiteRed [@litered] to use IBP [@chet] and LI identities [@gr] in an efficient manner. The MIs for the kinematic configuration of the problem at hand are analytically known from the seminal works of Gehrmann and Remiddi [@Gehrmann:2000zt]. We use them to obtain the unrenormalized coefficients in a Laurent series in $\epsilon$. In order to optimize the use of LiteRed, we have reduced all the one and two loop integrals to belong to few integral sets. This is done by shifting the loop momenta suitably using an in-house algorithm which uses FORM. We find that the sets for both one and two loop integrals are exactly same as those given in [@Ahmed:2014gla] for the case of massive spin-2 resonance $\rightarrow$ 3 gluons. The topologies of the appearing planar and non-planar master integrals are shown in fig.(\[fig:planar\]) and fig.(\[fig:nonplanar\]) respectively. For one-loop diagrams, the integral belongs to one of the following sets: $$\begin{aligned}
\label{onebasis}
\{ \cD, \hspace{0.1cm} \cD_{1}, \hspace{0.1cm} \cD_{12}, \hspace{0.1cm} \cD_{123} \} \, ,
\{ \cD, \hspace{0.1cm} \cD_{2}, \hspace{0.1cm} \cD_{23}, \hspace{0.1cm} \cD_{123} \} \, ,\{ \cD, \hspace{0.1cm} \cD_{3}, \hspace{0.1cm} \cD_{31}, \hspace{0.1cm} \cD_{123} \}\end{aligned}$$ where, $$\begin{aligned}
 \cD = k_1^2, \hspace{0.1cm} \cD_{i} = (k_1 - p_i)^2,
\hspace{0.1cm} \cD_{ij} = (k_1 - p_i - p_j)^2,
\hspace{0.1cm} \cD_{ijk} = (k_1 - p_i - p_j - p_k)^2 \; .\end{aligned}$$ At two loops, we have nine independent Lorentz invariants involving loop momenta $k_1$ and $k_2$, namely $\{ ( k_\alpha \cdot k_\beta ),
(k_\alpha \cdot p_i) \}, \alpha,\beta = 1, 2;\hspace{0.1cm} i = 1,...,3 $. Shifting of loop momenta allows us to express each two loop Feynman integral to contain terms belonging to one of the following six sets: $$\begin{aligned}
\label{twobasis}
 &&\{ \cD_0, \hspace{0.1cm} \cD_1, \hspace{0.1cm} \cD_2, \hspace{0.1cm} \cD_{1;1}, \hspace{0.1cm} \cD_{2;1}, \hspace{0.1cm}
 \cD_{1;12}, \hspace{0.1cm} \cD_{2;12}, \hspace{0.1cm} \cD_{1;123}, \hspace{0.1cm} \cD_{2;123} \} \,,
\nonumber\\[1.5ex]
 &&\{ \cD_0, \hspace{0.1cm} \cD_1, \hspace{0.1cm} \cD_2, \hspace{0.1cm} \cD_{1;2}, \hspace{0.1cm} \cD_{2;2}, \hspace{0.1cm}
 \cD_{1;23}, \hspace{0.1cm} \cD_{2;23}, \hspace{0.1cm} \cD_{1;123}, \hspace{0.1cm} \cD_{2;123} \} \,,
\nonumber\\[1.5ex]
 &&\{ \cD_0, \hspace{0.1cm} \cD_1, \hspace{0.1cm} \cD_2, \hspace{0.1cm} \cD_{1;3}, \hspace{0.1cm} \cD_{2;3}, \hspace{0.1cm}
 \cD_{1;31}, \hspace{0.1cm} \cD_{2;31}, \hspace{0.1cm} \cD_{1;123}, \hspace{0.1cm} \cD_{2;123} \} \,,
\nonumber\\[1.5ex]
 &&\{ \cD_0, \hspace{0.1cm} \cD_1, \hspace{0.1cm} \cD_2, \hspace{0.1cm} \cD_{1;1}, \hspace{0.1cm} \cD_{2;1}, \hspace{0.1cm}
 \cD_{0;3}, \hspace{0.1cm} \cD_{1;12}, \hspace{0.1cm} \cD_{2;12}, \hspace{0.1cm} \cD_{1;123} \} \,,
\nonumber\\[1.5ex]
 &&\{ \cD_0, \hspace{0.1cm} \cD_1, \hspace{0.1cm} \cD_2, \hspace{0.1cm} \cD_{1;2}, \hspace{0.1cm} \cD_{2;2}, \hspace{0.1cm}
 \cD_{0;1}, \hspace{0.1cm} \cD_{1;23}, \hspace{0.1cm} \cD_{2;23}, \hspace{0.1cm} \cD_{1;123} \} \,,
\nonumber\\[1.5ex]
 &&\{ \cD_0, \hspace{0.1cm} \cD_1, \hspace{0.1cm} \cD_2, \hspace{0.1cm} \cD_{1;3}, \hspace{0.1cm} \cD_{2;3}, \hspace{0.1cm}
 \cD_{0;2}, \hspace{0.1cm} \cD_{1;31}, \hspace{0.1cm} \cD_{2;31}, \hspace{0.1cm} \cD_{1;123} \}\end{aligned}$$ where, $$\begin{aligned}
 && \cD_0 = (k_1 - k_2)^2, \hspace{0.1cm} \cD_{\alpha} = k_{\alpha}^2, \hspace{0.1cm} \cD_{\alpha;i} = (k_{\alpha} - p_i)^2,
\hspace{0.1cm} \cD_{\alpha; ij} = (k_{\alpha} - p_i - p_j)^2,
\nonumber\\[1ex]
&& \cD_{0;i} = (k_1 - k_2 - p_i)^2, \hspace{0.1cm} \cD_{\alpha;ijk} = (k_{\alpha} - p_i - p_j - p_k)^2 \; .\end{aligned}$$ The UV singularities present in the bare coefficients are systematically removed using eqns.(\[renas\] & \[renl\]). The resulting UV finite coefficients do contain divergences from soft and collinear partons. In the next section we will demonstrate that our results correctly reproduce divergences described in the section \[sec:infrared\] at one and two loop level. We will also present the finite parts of the coefficients $A_m$ up to two loop level.

Results {#sec:result}
=======

In this section we present the results up to two loop level in QCD for the amplitude $H \rightarrow b +\overline b+g$ in the ${\overline {MS}}$ scheme. The results are presented after subtracting the one and two loop universal subtraction operators $I_b^{(i)}(\epsilon),i=1,2$ as described in the section \[sec:infrared\]. Following the eqns.(\[tencoeff\], \[coeffexpnd\] & \[renamp\]), the $l^{th}$ loop amplitude can be written as $$\label{rnMxpnd}
 |\rnM^{(l)} \rangle = 4 \pi ~ T^a_{ij} \Big\{ A_1^{(l)} {\rm T_1} + A_2^{(l)} {\rm T_2} \Big\} \;.$$ The renormalised coefficients $A_m^{(l)}$ are related to their bare counterparts $\hat A_m^{(l)}$ through (see eqn.(\[rln\])): $$\begin{aligned}
 \label{rln2}
 A_m^{(0)} &= \left(  \frac{1}{\mu^{\ep}_R} \right)^{\frac{1}{2}}    \hat{A}_m^{(0)} \; ,
\nonumber\\[1ex]
 A_m^{(1)} &= \left(  \frac{1}{\mu^{\ep}_R} \right)^{\frac{3}{2}} 
\left[ ~ \hat{A}_m^{(1)} + \mu^{\ep}_R \Big( \frac{r_{a_1}}{2} + r_{\lambda_1} \Big) \hat{A}_m^{(0)} ~ \right] \; ,
\nonumber\\[1ex]
 A_m^{(2)} &= \left( \frac{1}{\mu^{\ep}_R} \right)^{\frac{5}{2}} 
\Big[ ~ \hat{A}_m^{(2)}  
+ \mu^{\ep}_R \Big( \frac{3r_{a_1}}{2} + r_{\lambda_1} \Big)  \hat{A}_m^{(1)} 
\nonumber\\
& \qquad \qquad \quad
+ \mu^{2\ep}_R \left( \frac{r_{a_2}}{2} - \frac{r_{a_1}^2}{8} + \frac{r_{a_1}}{2} r_{\lambda_1}
+ r_{\lambda_2}  \right)  \hat{A}_m^{(0)} ~ \Big] \;.\end{aligned}$$ Using the procedure discussed in the previous section, we first compute the bare coefficients $\hat A_m^{(l)}$ and the eqns.(\[rln2\]) give the renormalized coefficients. The finite parts of the coefficients $A_m^{(l)}$ are defined after subtracting terms proportional to universal subtraction terms $I_b^{(l)}$ as follows $$\begin{aligned}
\label{catani2}
 A_m^{(1)} &= 2 \hspace{0.1cm} {\bf{I}}_{b}^{(1)} (\ep) \hspace{0.1cm} A_m^{(0)}
+ A_m^{(1)fin} \,,
\nonumber\\[2ex] 
 A_m^{(2)} &= 2 \hspace{0.1cm} {\bf{I}}_{b}^{(1)} (\ep) \hspace{0.1cm} A_m^{(1)}
+ 4 \hspace{0.1cm} {\bf{I}}^{(2)}_{b} (\ep) \hspace{0.1cm} A_m^{(0)}
+ A_m^{(2)fin}\end{aligned}$$ where, we have used eqns.(\[tencoeff\] & \[catani1\]).

Expanding the right sides of eqns.(\[rln2\] & \[catani2\]) in powers of $\ep$, we find that the infrared poles agree exactly, providing a crucial test on the correctness of our computation. The finite parts of the coefficients have the following expansions: $$\begin{aligned}
 & A_m^{(l)fin} = \sum_{n=0}^l ~ A_m^{(0)} ~ {\cal B}^{(l)}_{m ; n} \ln^n \Big( - \frac{Q^2}{\mu^2} \Big) \end{aligned}$$ where, $$A_1^{(0)} = - \frac{4 i}{t ~ u} \hspace{1cm} \text{and} \hspace{1cm} A_2^{(0)} = i \Big( \frac{1}{t} + \frac{1}{u}  \Big)$$ and the remaining coefficients ${\cal B}^{(l)}_{m;n}$ are given in the appendix. We also performed an independent computation of $\langle {\cal M}^{(0)}|{\cal M}^{(l)}\rangle$ for $l=1,2$ without using any projectors and then compared against one obtained using the projectors, i.e using the coefficients $A_m^{(l)}$. We find both give the same result, providing an independent check on our computation.

Following [@Gehrmann:2002zr] [^2], we have obtained results for the crossed reactions given in eqn. (\[subprocess\]) relevant for Higgs+1 jet production at hadron colliders. The corresponding finite coefficients $A_m^{(l)fin}$ are attached with the arXiv submission.

Conclusions {#sec:conc}
===========

We have presented the amplitudes for the partonic subprocess $H \rightarrow b+\overline b+g$ and other subprocesses related by crossing, up to two loop level in QCD that contribute to exclusive observables involving Higgs boson and a jet. The dominant one is from gluon gluon fusion which is already known to this accuracy. We have used dimensional regularization to perform our computation. Using appropriate projectors, the amplitude is expressed in terms of two scalar coefficients $A_m$. We have found that the infrared structure of the amplitude is according to Catani’s prediction on QCD amplitudes upto two loop level. Also, the coefficient of single pole term is found to be in agreement with predictions based on the observation of the universal behavior of poles in the multi-parton QCD amplitudes.

Acknowledgments {#acknowledgments .unnumbered}
===============

TA, MM and NR thank the Institute of Mathematical Sciences (IMSc) for the hospitality during the course of the work. We thank the staff of IMSc computer center for their help. We sincerely thank T. Gehrmann for providing us the master integrals and analytically continued HPLs and 2d-HPLs required for our computation. We thank R. N. Lee for his help with LiteRed. Finally, we would like to thank K. Hasegawa, M. K. Mandal and L. Tancredi for useful discussions. The work of TA, MM and NR has been partially supported by funding from RECAPP, Department of Atomic Energy, Govt. of India.

Harmonic polylogarithms
=======================

Here, we provide the definition of HPL and 2dHPL. HPL is represented by $H(\vec{m}_w;y)$ with a $w$-dimensional vector $\vec{m}_w$ of parameters and its argument $y$. The elements of $\vec{m}_w$ belong to $\{ 1, 0, -1 \}$ through which we define the following rational functions $$f(1;y) \equiv \frac{1}{1-y}, \qquad f(0;y) \equiv \frac{1}{y},  \qquad f(-1;y) \equiv \frac{1}{1+y} \, .$$ The weight 1 $(w = 1)$ HPLs are $$H(1, y) \equiv - \ln (1 - y), \qquad  H(0, y) \equiv \ln y, \qquad  H(-1, y) \equiv \ln (1 + y) \, .$$ For $w > 1$, the definition of $H(m, \vec{m}_{w};y)$ is $$\label{1dhpl}
 H(m, \vec{m}_w;y) \equiv \int_0^y dx ~ f(m, x) ~ H(\vec{m}_w;x),  \qquad \qquad  m \in 0, \pm 1  \, .$$ The 2dHPLs are defined in the same way as eqn.(\[1dhpl\]) with the new elements $\{ 2, 3 \}$ in $\vec{m}_w$ representing a new class of rational functions $$f(2;y) \equiv f(1-z;y) \equiv \frac{1}{1-y-z}, \qquad f(3;y) \equiv f(z;y) \equiv \frac{1}{y+z}$$ and correspondingly with the weight 1 $(w = 1)$ 2dHPLs $$H(2, y) \equiv - \ln \Big(1 - \frac{y}{1-z} \Big), \qquad  H(3, y) \equiv \ln \Big( \frac{y+z}{z} \Big) \, .$$

One-loop coefficients
=====================

[B]{}\^[(1)]{}\_[1 ; 1]{} = (-11 C\_A - 18 C\_F + 2 n\_f)

[B]{}\^[(1)]{}\_[1 ; 0]{} = ( -6 H(0,y) H(0,z)-6 H(0,y) H(1,z)-6 H(2,y) H(0,z)+12 H(3,y) H(1,z)-10 H(0,y)-9 H(2,y)-6 H(0,2,y) -6 H(2,0,y)+12 H(3,2,y)-10 H(0,z)-9 H(1,z)+6 H(0,1,z)-6 H(1,0,z)- 6 \_2 ) + C\_F ( 2 H(0,y) H(1,z)- 4 H(3,y) H(1,z) + 2 H(2,y) H(0,z) + 3 H(2,y)+12 H(0,2,y)- 2 H(1,0,y)+ 2 H(2,0,y)- 4 H(3,2,y)+ 3 H(1,z)- 2 H(0,1,z)- 2 ) + n\_f ( H(0,y)+H(0,z) )

[B]{}\^[(1)]{}\_[2 ; 1]{} = (-11 C\_A-18 C\_F + 2 n\_f)

[B]{}\^[(1)]{}\_[2 ; 0]{} = ( - 6 H(0,y) H(0,z)-6 H(0,y) H(1,z)-6 H(2,y) H(0,z)+12 H(3,y) H(1,z)-10 H(0,y)-9 H(2,y)-6 H(0,2,y) -6 H(2,0,y) +12 H(3,2,y)-10 H(0,z)-9 H(1,z)+6 H(0,1,z)-6 H(1,0,z) - 6 \_2 + 6 ) + C\_F ( 2 H(0,y) H(1,z)-4 H(3,y) H(1,z) +2 H(2,y) H(0,z) +3 H(2,y)+2 H(0,2,y)-2 H(1,0,y)+2 H(2,0,y)-4 H(3,2,y)+3 H(1,z)-2 H(0,1,z)-3 ) + n\_f ( H(0,y)+H(0,z) )

Two-loop coefficients
=====================

[B]{}\^[(2)]{}\_[1 ; 2]{} = ( 121 C\_A\^2 + 44 C\_A ( 6 C\_F - n\_f) + 4 ( 27 C\_F\^2 - 12 C\_F n\_f + n\_f\^2 ) )

[B]{}\^[(2)]{}\_[1 ; 1]{} = C\_A\^2 ( 594 H(0,y) H(0,z)+594 H(0,y) H(1,z)+594 H(2,y) H(0,z)-1188 H(3,y) H(1,z)+990 H(0,y)+891 H(2,y) +594 H(0,2,y)+594 H(2,0,y)-1188 H(3,2,y)+990 H(0,z)+891 H(1,z)-594 H(0,1,z)+594 H(1,0,z)+495 \_2-108 \_3 -702 ) + C\_A C\_F ( -1188 H(0,y) H(1,z)+54 H(0,y) (6 H(0,z)+6 H(1,z)+10)-1188 H(2,y) H(0,z)+54 (6 H(2,y)+10) H(0,z)+1728 H(3,y) H(1,z)-1296 H(2,y)-864 H(0,2,y)+1188 H(1,0,y)-864 H(2,0,y)+1728 H(3,2,y)-1296 H(1,z)+864 H(0,1,z)+324 H(1,0,z)+1566 \_2-2808 \_3-1048 ) + C\_F\^2 ( -648 H(0,y) H(1,z)-648 H(2,y) H(0,z)

[=]{} +1296 H(3,y) H(1,z)-972 H(2,y)-648 H(0,2,y)+648 H(1,0,y)-648 H(2,0,y)+1296 H(3,2,y)-972 H(1,z)+648 H(0,1,z)-1296 \_2 + 2592 \_3 + 648 ) + C\_A n\_f ( - 2 (-27 (4 H(3,y)-3) H(1,z)+81 H(2,y)+54 H(0,2,y)+54 H(2,0,y)-108 H(3,2,y)-54 H(0,1,z)+54 H(1,0,z) + 45 \_2 - 206 ) -9 H(0,y) (12 H(0,z)+12 H(1,z)+31)-9 (12 H(2,y)+31) H(0,z) ) + C\_F n\_f ( 108 H(0,y) H(1,z)+108 H(2,y) H(0,z)-216 H(3,y) H(1,z)-27 H(0,y)+162 H(2,y)+108 H(0,2,y)-108 H(1,0,y)+108 H(2,0,y)-216 H(3,2,y)-27 H(0,z)+162 H(1,z)-108 H(0,1,z) - 54 \_2 + 32 ) + n\_f\^2 ( 9 H(0,y)+9 H(0,z)-20 )

[B]{}\^[(2)]{}\_[1 ; 0]{} = C\_A\^2 { \_4 ( 39/8 ) + \_3 ( - H(0,y) - 6 H(1,y) + 11 H(2,y) - H(0,z) + 5 H(1,z) + 6 /s + 407/36 ) + \_2 ( 108 s\^2 H(0,y) H(0,z)

[=]{} +72 s\^2 H(1,y) H(0,z)+108 s\^2 H(0,y) H(1,z)+72 s\^2 H(1,y) H(1,z)+36 s\^2 H(2,y) H(0,z)-72 s\^2 H(3,y) H(1,z)+147 s\^2 H(0,y)+453 s\^2 H(2,y)+108 s\^2 H(0,2,y)+72 s\^2 H(1,2,y)+36 s\^2 H(2,0,y)-72 s\^2 H(3,2,y)+147 s\^2 H(0,z)+237 s\^2 H(1,z)+36 s\^2 H(0,1,z)+36 s\^2 H(1,0,z)+72 s\^2 H(1,1,z)-216 s H(1,y)-108 s t H(0,z)-108 s u H(0,y)-216 s u H(1,z)-40 s\^2+96 t u ) / (36 s\^2) + ( - H(2,y) H(1,z)+ H(3,y) H(1,z)+7 H(0,0,y) H(1,z)- H(0,2,y) H(1,z)-7 H(0,3,y) H(1,z)- H(2,0,y) H(1,z)+ H(2,3,y) H(1,z)-7 H(3,0,y) H(1,z)+ H(3,2,y) H(1,z)+ H(3,3,y) H(1,z)+2 H(0,0,2,y) H(1,z)+2 H(0,2,0,y) H(1,z)+4 H(0,3,0,y) H(1,z)+4 H(0,3,3,y) H(1,z)+2 H(1,0,3,y) H(1,z)-2 H(1,2,3,y) H(1,z)+2 H(2,0,0,y) H(1,z)+2 H(2,0,3,y) H(1,z)+2 H(2,1,0,y) H(1,z)+4 H(2,2,3,y) H(1,z)+2 H(2,3,0,y) H(1,z)-8 H(2,3,3,y) H(1,z)+4 H(3,0,3,y) H(1,z)+4 H(3,3,0,y) H(1,z)-16 H(3,3,3,y) H(1,z)

[=]{} - H(1,z)- H(2,y)+ H(0,0,y)+7 H(2,y) H(0,0,z)+2 H(0,0,y) H(0,0,z)+ H(0,0,z)+ H(2,y) H(0,1,z)+ H(3,y) H(0,1,z)+2 H(0,0,y) H(0,1,z)+ H(0,1,z)+2 H(0,0,z) H(0,2,y)-2 H(0,1,z) H(0,2,y)- H(0,2,y)-2 H(0,1,z) H(0,3,y)+++6 H(2,y) H(1,0,z)-7 H(3,y) H(1,0,z)+2 H(0,0,y) H(1,0,z)-6 H(0,3,y) H(1,0,z)- H(1,0,z)+ H(3,y) H(1,1,z)+2 H(0,0,y) H(1,1,z)- H(1,1,z)-2 H(0,1,z) H(1,2,y) [}]{}

[=]{} [{]{} +2 H(1,0,z) H(1,2,y)+2 H(0,0,z) H(2,0,y)+2 H(0,1,z) H(2,0,y)+2 H(1,0,z) H(2,0,y)- H(2,0,y)+2 H(0,0,z) H(2,2,y)+4 H(0,1,z) H(2,2,y)+4 H(1,0,z) H(2,2,y)- H(2,2,y)-6 H(0,1,z) H(2,3,y)+2 H(1,0,z) H(2,3,y)+2 H(0,1,z) H(3,0,y)-2 H(1,0,z) H(3,0,y)+ H(3,2,y)-12 H(0,1,z) H(3,3,y)+4 H(1,0,z) H(3,3,y)-2 H(1,y) H(0,0,1,z)-2 H(2,y) H(0,0,1,z)-8 H(3,y) H(0,0,1,z)+ H(0,0,1,z)+7 H(0,0,2,y)+++2 H(1,y) H(0,1,0,z)+6 H(2,y) H(0,1,0,z)+4 H(3,y) H(0,1,0,z)+ H(0,1,1,z)+7 H(0,2,0,y)- H(0,2,2,y)-7 H(0,3,2,y)+2 H(2,y) H(1,0,0,z)+7 H(1,0,0,z)-2 H(1,y) H(1,0,1,z)+6 H(2,y) H(1,0,1,z)+ H(1,0,1,z)--+2 H(1,y) H(1,1,0,z)+6 H(2,y) H(1,1,0,z)+H(0,y) (-++(-) H(0,z)- H(1,z)+7 H(0,0,z)- H(0,1,z)+ H(1,0,z)- H(1,1,z)-2 H(0,0,1,z)-2 H(0,1,1,z)+2 H(1,0,0,z)+2 H(1,1,0,z)-)+7 H(2,0,0,y)- H(2,0,2,y)+ H(2,1,0,y)- H(2,2,0,y)+ H(2,3,2,y)-7 H(3,0,2,y)-7 H(3,2,0,y)+ H(3,2,2,y)+ H(3,3,2,y)+H(0,z) (- H(2,y)+7 H(0,0,y)+ H(0,2,y)-+ H(2,0,y)- H(2,2,y)-7 H(3,2,y)+2 H(0,0,2,y)+2 H(0,2,0,y)-2 H(0,2,2,y)-6 H(0,3,2,y)-2 H(1,0,2,y)+2 H(2,0,0,y)+2 H(2,2,0,y)+2 H(2,3,2,y)-2 H(3,0,2,y)-2 H(3,2,0,y)+4 H(3,3,2,y)-)+8 H(0,0,1,0,z)+2 H(0,0,1,1,z)+2 H(0,0,2,2,y)+2 H(0,1,0,1,z)+6 H(0,1,1,0,z)+2 H(0,2,0,2,y)

[=]{} +4 H(0,2,1,0,y)+2 H(0,2,2,0,y)+4 H(0,3,0,2,y)+4 H(0,3,2,0,y)+4 H(0,3,3,2,y)+6 H(1,0,1,0,z)+2 H(1,0,3,2,y)+2 H(1,1,0,0,z)+4 H(1,1,0,1,z)+6 H(1,1,1,0,z)-2 H(1,2,3,2,y)+2 H(2,0,0,2,y)+2 H(2,0,1,0,y)+2 H(2,0,2,0,y)+2 H(2,0,3,2,y)+2 H(2,1,0,2,y)+2 H(2,1,2,0,y)+2 H(2,2,0,0,y)+4 H(2,2,1,0,y)+4 H(2,2,3,2,y)+2 H(2,3,0,2,y)+2 H(2,3,2,0,y)-8 H(2,3,3,2,y)+4 H(3,0,3,2,y)+4 H(3,3,0,2,y)+4 H(3,3,2,0,y)-16 H(3,3,3,2,y)- ) }

[=]{} + C\_A C\_F { \_4 ( ) + \_3 ( 22 H(1,y)-18 H(2,y)+4 H(1,z)++-+ ) + \_2 ( ----++++-------2 H(1,y) H(0,z)-6 H(0,y) H(1,z)-2 H(1,y) H(1,z)-4 H(2,y) H(0,z)+ H(1,y)+2 H(0,2,y)+6 H(1,0,y)-2 H(1,2,y)-4 H(2,0,y)-8 H(2,1,y)+8 H(2,2,y)+2 H(0,1,z)+2 H(1,0,z)-2 H(1,1,z)---- ) + ( --++

[=]{} ++-++-+--+----++-++++ H(1,z)+2 H(1,z) H(2,y)+ H(2,y)- H(1,z) H(3,y)-14 H(1,z) H(0,0,y)-14 H(2,y) H(0,0,z)- H(2,y) H(0,1,z)- H(3,y) H(0,1,z)-4 H(0,0,y) H(0,1,z)- H(0,1,z)- H(1,z) H(0,2,y)-4 H(0,0,z) H(0,2,y)+6 H(0,1,z) H(0,2,y)+ H(0,2,y)+20 H(1,z) H(0,3,y)+4 H(0,1,z) H(0,3,y)+-+3 H(1,z) H(1,0,y)

[=]{} -2 H(0,1,z) H(1,0,y)---+20 H(3,y) H(1,0,z)-4 H(0,0,y) H(1,0,z)+8 H(0,2,y) H(1,0,z)+16 H(0,3,y) H(1,0,z)+2 H(1,0,y) H(1,0,z)+ H(3,y) H(1,1,z)-8 H(0,0,y) H(1,1,z)+4 H(0,3,y) H(1,1,z)+2 H(1,1,z)+6 H(0,1,z) H(1,2,y)-2 H(1,0,z) H(1,2,y)- H(1,z) H(2,0,y)-4 H(0,0,z) H(2,0,y)-2 H(0,1,z) H(2,0,y)-6 H(1,0,z) H(2,0,y)+ H(2,0,y)-8 H(0,0,z) H(2,2,y)-16 H(0,1,z) H(2,2,y)-8 H(1,0,z) H(2,2,y)+2 H(2,2,y)- H(1,z) H(2,3,y)+12 H(0,1,z) H(2,3,y)-4 H(1,0,z) H(2,3,y)+20 H(1,z) H(3,0,y)-4 H(0,1,z) H(3,0,y)+4 H(1,1,z) H(3,0,y)+ H(1,z) H(3,2,y)-4 H(0,1,z) H(3,2,y)-4 H(1,0,z) H(3,2,y)- H(3,2,y)- H(1,z) H(3,3,y)+28 H(0,1,z) H(3,3,y) -12 H(1,0,z) H(3,3,y)-8 H(1,1,z) H(3,3,y)+6 H(1,y) H(0,0,1,z)+16 H(3,y) H(0,0,1,z)- H(0,0,1,z) -8 H(1,z) H(0,0,2,y)-14 H(0,0,2,y)+4 H(1,z) H(0,0,3,y)++-+2 H(1,z) H(0,1,0,y)+-2 H(1,y) H(0,1,0,z)-12 H(2,y) H(0,1,0,z)-4 H(3,y) H(0,1,0,z)

[=]{} -4 H(3,y) H(0,1,1,z)+ H(0,1,1,z)-8 H(1,z) H(0,2,0,y)-14 H(0,2,0,y)- H(0,2,2,y)+2 H(1,z) H(0,2,3,y)-8 H(1,z) H(0,3,0,y)+4 H(1,z) H(0,3,2,y)+20 H(0,3,2,y)-12 H(1,z) H(0,3,3,y)+4 H(1,z) H(1,0,0,y)+14 H(1,0,0,y)-4 H(2,y) H(1,0,0,z)+6 H(1,y) H(1,0,1,z)-18 H(2,y) H(1,0,1,z)-4 H(3,y) H(1,0,1,z)- H(1,0,1,z)+3 H(1,0,2,y)-6 H(1,z) H(1,0,3,y)+ H(1,1,0,y)+H(0,y) ((-++2) H(0,z)+(+) H(1,z)+ (---+-+13 H(0,1,z)--+-14 H(1,1,z)+12 H(0,0,1,z)+6 H(0,1,0,z)+24 H(0,1,1,z)+6 H(1,0,1,z)-12 H(1,1,0,z)+))-++-2 H(1,y) H(1,1,0,z)-18 H(2,y) H(1,1,0,z)-4 H(3,y) H(1,1,0,z)+3 H(1,2,0,y)+6 H(1,z) H(1,2,3,y)-8 H(1,z) H(2,0,0,y)-14 H(2,0,0,y)- H(2,0,2,y)-4 H(1,z) H(2,0,3,y)

[=]{} +-2 H(1,z) H(2,1,0,y)- H(2,2,0,y)-16 H(1,z) H(2,2,3,y)-4 H(1,z) H(2,3,0,y)- H(2,3,2,y)+16 H(1,z) H(2,3,3,y)+4 H(1,z) H(3,0,2,y)+20 H(3,0,2,y)-12 H(1,z) H(3,0,3,y)+4 H(1,z) H(3,2,0,y)+20 H(3,2,0,y)+ H(3,2,2,y)-8 H(1,z) H(3,2,3,y)-12 H(1,z) H(3,3,0,y)+H(0,z) ((+) H(2,y)+ (----+-47 H(0,2,y)+(-++28) H(1,0,y)--14 H(2,2,y)+60 H(3,2,y)-24 H(0,0,2,y)+6 H(0,1,0,y)+12 H(0,2,2,y)+48 H(0,3,2,y)+12 H(1,0,0,y)+12 H(1,0,2,y)+6 H(1,2,0,y)-12 H(2,0,0,y)+6 H(2,0,2,y)-6 H(2,1,0,y)-12 H(2,2,0,y)-12 H(2,3,2,y)+24 H(3,0,2,y)+12 H(3,2,2,y)-36 H(3,3,2,y)+))-8 H(1,z) H(3,3,2,y)- H(3,3,2,y)+40 H(1,z) H(3,3,3,y)+4 H(0,0,1,0,y)-8 H(0,0,1,0,z)-8 H(0,0,1,1,z)-8 H(0,0,2,2,y)+4 H(0,0,3,2,y)-8 H(0,1,0,1,z)+2 H(0,1,0,2,y)-12 H(0,1,1,0,z)+2 H(0,1,2,0,y)-8 H(0,2,0,2,y)-2 H(0,2,1,0,y)-8 H(0,2,2,0,y)+2 H(0,2,3,2,y)-8 H(0,3,0,2,y)-8 H(0,3,2,0,y)+4 H(0,3,2,2,y)-12 H(0,3,3,2,y)-2 H(1,0,0,1,z)+4 H(1,0,0,2,y)+4 H(1,0,1,0,y)-4 H(1,0,1,0,z)+4 H(1,0,2,0,y)-6 H(1,0,3,2,y)-12 H(1,1,0,1,z)-14 H(1,1,1,0,z)+4 H(1,2,0,0,y)+4 H(1,2,1,0,y)+6 H(1,2,3,2,y)-8 H(2,0,0,2,y)-8 H(2,0,2,0,y)-4 H(2,0,3,2,y)+4 H(2,1,0,0,y)-2 H(2,1,0,2,y)-8 H(2,1,1,0,y)-2 H(2,1,2,0,y)-8 H(2,2,0,0,y)-8 H(2,2,1,0,y)-16 H(2,2,3,2,y)-4 H(2,3,0,2,y)-4 H(2,3,2,0,y)+16 H(2,3,3,2,y)+8 H(3,0,1,0,y)+4 H(3,0,2,2,y)-12 H(3,0,3,2,y)+4 H(3,2,0,2,y)-8 H(3,2,1,0,y)+4 H(3,2,2,0,y)-8 H(3,2,3,2,y)-12 H(3,3,0,2,y)-12 H(3,3,2,0,y)-8 H(3,3,2,2,y)+40 H(3,3,3,2,y)---+-+- ) }

[=]{} + C\_F\^2 { \_4 (-22) + \_3 ( -16 H(1,y)+8 H(2,y)-8 H(1,z)--+-30 ) + \_2 ( +++----++++++-8 H(2,y) H(1,z)+8 H(0,1,y)-8 H(0,2,y)+8 H(1,1,y)+8 H(2,1,y)-16 H(2,2,y)++++16 ) + ( -+-+--+-++-++-++-+---

[=]{} -18 H(1,z)+9 H(1,z) H(2,y)-18 H(2,y)+8 H(1,z) H(3,y)+12 H(3,y) H(0,1,z)+4 H(0,1,z)+12 H(1,z) H(0,2,y)-4 H(0,1,z) H(0,2,y)-4 H(0,2,y)-12 H(1,z) H(0,3,y)-+-6 H(1,z) H(1,0,y)

[=]{} +4 H(0,1,z) H(1,0,y)++-12 H(3,y) H(1,0,z)-8 H(0,2,y) H(1,0,z)-8 H(0,3,y) H(1,0,z)-24 H(3,y) H(1,1,z)+8 H(0,0,y) H(1,1,z)-8 H(0,3,y) H(1,1,z)+9 H(1,1,z)-4 H(0,1,z) H(1,2,y)+12 H(1,z) H(2,0,y)-4 H(0,1,z) H(2,0,y)-4 H(2,0,y)+8 H(0,0,z) H(2,2,y)+16 H(0,1,z) H(2,2,y)+9 H(2,2,y)-12 H(1,z) H(2,3,y)-12 H(1,z) H(3,0,y)-8 H(1,1,z) H(3,0,y)-24 H(1,z) H(3,2,y)+8 H(0,1,z) H(3,2,y)+8 H(3,2,y)+24 H(1,z) H(3,3,y)-8 H(0,1,z) H(3,3,y)+8 H(1,0,z) H(3,3,y)+16 H(1,1,z) H(3,3,y)-4 H(1,y) H(0,0,1,z)+8 H(2,y) H(0,0,1,z)+8 H(1,z) H(0,0,2,y)-8 H(1,z) H(0,0,3,y)-++-4 H(1,z) H(0,1,0,y)+4 H(2,y) H(0,1,0,z)+8 H(3,y) H(0,1,1,z)-12 H(0,1,1,z)+8 H(1,z) H(0,2,0,y)+12 H(0,2,2,y)-4 H(1,z) H(0,2,3,y)

[=]{} -8 H(1,z) H(0,3,2,y)-12 H(0,3,2,y)+8 H(1,z) H(0,3,3,y)-8 H(1,z) H(1,0,0,y)+2 H(0,y) (--+++(--2) H(1,z)+6 H(1,1,z)-4 H(0,1,1,z)-2 H(1,0,1,z))-4 H(1,y) H(1,0,1,z)+12 H(2,y) H(1,0,1,z)+8 H(3,y) H(1,0,1,z)-6 H(1,0,1,z)-6 H(1,0,2,y)+4 H(1,z) H(1,0,3,y)+--+4 H(2,y) H(1,1,0,z)-6 H(1,2,0,y)-4 H(1,z) H(1,2,3,y)+8 H(1,z) H(2,0,0,y)+12 H(2,0,2,y)--4 H(1,z) H(2,1,0,y)+12 H(2,2,0,y)+16 H(1,z) H(2,2,3,y)-12 H(2,3,2,y)-8 H(1,z) H(3,0,2,y)-12 H(3,0,2,y)+8 H(1,z) H(3,0,3,y)-8 H(1,z) H(3,2,0,y)-12 H(3,2,0,y)-24 H(3,2,2,y)+16 H(1,z) H(3,2,3,y)+8 H(1,z) H(3,3,0,y)+16 H(1,z) H(3,3,2,y)+24 H(3,3,2,y)+2 H(0,z) (++++(--2) H(2,y)+6 H(0,2,y)-

[=]{} ++6 H(2,2,y)-6 H(3,2,y)+4 H(0,0,2,y)-4 H(0,3,2,y)-2 H(2,0,2,y)-4 H(3,0,2,y)-4 H(3,2,2,y)+4 H(3,3,2,y))-16 H(1,z) H(3,3,3,y)+8 H(0,0,1,1,z)+8 H(0,0,2,2,y)-8 H(0,0,3,2,y)+8 H(0,1,0,1,z)-4 H(0,1,0,2,y)+8 H(0,1,1,0,y)+8 H(0,1,1,0,z)-4 H(0,1,2,0,y)+8 H(0,2,0,2,y)-4 H(0,2,1,0,y)+8 H(0,2,2,0,y)-4 H(0,2,3,2,y)-8 H(0,3,2,2,y)+8 H(0,3,3,2,y)+4 H(1,0,0,1,z)-8 H(1,0,0,2,y)+4 H(1,0,1,0,y)-8 H(1,0,2,0,y)+4 H(1,0,3,2,y)+8 H(1,1,0,0,y)+8 H(1,1,0,1,z)+8 H(1,1,1,0,y)+8 H(1,1,1,0,z)-8 H(1,2,0,0,y)-4 H(1,2,1,0,y)-4 H(1,2,3,2,y)+8 H(2,0,0,2,y)-4 H(2,0,1,0,y)+8 H(2,0,2,0,y)-8 H(2,1,0,0,y)-4 H(2,1,0,2,y)+8 H(2,1,1,0,y)-4 H(2,1,2,0,y)+8 H(2,2,0,0,y)+16 H(2,2,3,2,y)-8 H(3,0,1,0,y)-8 H(3,0,2,2,y)+8 H(3,0,3,2,y)-8 H(3,2,0,2,y)+8 H(3,2,1,0,y)-8 H(3,2,2,0,y)+16 H(3,2,3,2,y)+8 H(3,3,0,2,y)+8 H(3,3,2,0,y)+16 H(3,3,2,2,y)-16 H(3,3,3,2,y)+6++-+-+ ) }

[=]{} + C\_A n\_f { \_3 ( - ) + \_2 ( - H(0,y)+ H(2,y)- H(0,z)+ H(1,z)+- ) + ( H(0,y) (H(0,z) (+20)+31 H(1,z)-36 H(0,0,z)+24 H(0,1,z)-36 H(1,0,z)+24 H(1,1,z)++86)+++H(0,z) ( H(2,y)-H(0,0,y)- H(0,2,y)-H(2,0,y)+ H(2,2,y)+H(3,2,y)++)+H(2,y) H(1,z)- H(3,y) H(1,z)-H(0,0,y) H(1,z)+ H(0,2,y) H(1,z)+H(0,3,y) H(1,z)+ H(2,0,y) H(1,z)- H(2,3,y) H(1,z)+H(3,0,y) H(1,z)- H(3,2,y) H(1,z)- H(3,3,y) H(1,z)-H(2,y) H(0,0,z)- H(2,y) H(0,1,z)- H(3,y) H(0,1,z)+H(3,y) H(1,0,z)- H(3,y) H(1,1,z)+ H(2,y)- H(0,0,y)+ H(0,2,y)+ H(2,0,y)+H(2,2,y)- H(3,2,y)-H(0,0,2,y)-H(0,2,0,y)+ H(0,2,2,y)+H(0,3,2,y)-H(2,0,0,y)+ H(2,0,2,y)- H(2,1,0,y)+ H(2,2,0,y)- H(2,3,2,y)+H(3,0,2,y)+H(3,2,0,y)- H(3,2,2,y)- H(3,3,2,y)+ H(1,z)- H(0,0,z)- H(0,1,z)+ H(1,0,z)+H(1,1,z)- H(0,0,1,z)- H(0,1,1,z)-H(1,0,0,z)- H(1,0,1,z)+ ) }

[=]{} + C\_F n\_f { \_3 ( - ) + \_2 ( (4 H(1,y)-5 H(2,y)-H(1,z)+1) ) + ( - H(0,y) H(1,z)-108 H(0,y) H(0,1,z)+27 H(0,y) H(1,0,z)-108 H(0,y) H(1,1,z)- H(2,y) H(0,z)-162 H(2,y) H(1,z)+360 H(3,y) H(1,z)+162 H(0,0,y) H(1,z)+162 H(2,y) H(0,0,z)+108 H(2,y) H(0,1,z)+54 H(3,y) H(0,1,z)+54 H(0,2,y) H(0,z)-108 H(0,2,y) H(1,z)-162 H(0,3,y) H(1,z)-27 H(1,0,y) H(0,z)-108 H(2,y) H(1,0,z)-162 H(3,y) H(1,0,z)+216 H(3,y) H(1,1,z)+54 H(2,0,y) H(0,z)-108 H(2,0,y) H(1,z)-108 H(2,2,y) H(0,z)+216 H(2,3,y) H(1,z)-162 H(3,0,y) H(1,z)-162 H(3,2,y) H(0,z)+216 H(3,2,y) H(1,z)+216 H(3,3,y) H(1,z)-27 H(0,y)-102 H(2,y)- H(0,2,y)+180 H(1,0,y)- H(2,0,y)-162 H(2,2,y)+360 H(3,2,y)+162 H(0,0,2,y)-27 H(0,1,0,y)+162 H(0,2,0,y)-108 H(0,2,2,y)-162 H(0,3,2,y)-162 H(1,0,0,y)+108 H(1,1,0,y)+162 H(2,0,0,y)-108 H(2,0,2,y)-108 H(2,2,0,y)+216 H(2,3,2,y)-162 H(3,0,2,y)-162 H(3,2,0,y)+216 H(3,2,2,y)+216 H(3,3,2,y)-27 H(0,z)-102 H(1,z)+ H(0,1,z)+ H(1,0,z)-162 H(1,1,z)+54 H(0,0,1,z)-27 H(0,1,0,z)+108 H(0,1,1,z)+108 H(1,0,1,z)+200 ) / 81 }

[=]{} + n\_f\^2 { ( H(0,y) (3 H(0,z)-20)+15 H(0,0,y)-20 H(0,z)+15 H(0,0,z)+6 \_2 ) }

[B]{}\^[(2)]{}\_[2 ; 2]{} = (121 [C\_A]{}\^2+44 [C\_A]{} (6 [C\_F]{}-[n\_f]{})+4 (27 [C\_F]{}\^2-12 [C\_F]{} [n\_f]{}+[n\_f]{}\^2))

[B]{}\^[(2)]{}\_[2 ; 1]{} = C\_A\^2 { ( 594 H(0,y) H(0,z)+594 H(0,y) H(1,z)+594 H(2,y) H(0,z)-1188 H(3,y) H(1,z)+990 H(0,y)+891 H(2,y)+594 H(0,2,y)+594 H(2,0,y)-1188 H(3,2,y)+990 H(0,z)+891 H(1,z)-594 H(0,1,z)+594 H(1,0,z)+495 \_2-108 \_3 - 1296 ) } + C\_A C\_F { (-1188 H(0,y) H(1,z)+54 H(0,y) (6 H(0,z)+6 H(1,z)+10)-1188 H(2,y) H(0,z)+54 (6 H(2,y)+10) H(0,z)+1728 H(3,y) H(1,z)-1296 H(2,y)-864 H(0,2,y)+1188 H(1,0,y)-864 H(2,0,y)+1728 H(3,2,y)-1296 H(1,z)+864 H(0,1,z)+324 H(1,0,z)+1566 [\_2]{}-2808 [\_3]{}-778 ) } + C\_F\^2 { ( -648 H(0,y) H(1,z)-648 H(2,y) H(0,z)+1296 H(3,y) H(1,z)-972 H(2,y)-648 H(0,2,y)+648 H(1,0,y)-648 H(2,0,y)+1296 H(3,2,y)-972 H(1,z)+648 H(0,1,z)

[=]{} -1296 [\_2]{}+2592 [\_3]{}+972 ) } + C\_A n\_f { ( -2 ( -27 (4 H(3,y)-3) H(1,z)+81 H(2,y)+54 H(0,2,y)+54 H(2,0,y)-108 H(3,2,y)-54 H(0,1,z)+54 H(1,0,z)+45 [\_2]{}-260 )-9 H(0,y) (12 H(0,z)+12 H(1,z)+31)-9 (12 H(2,y)+31) H(0,z) ) } + C\_F n\_f { ( -2 (54 ( 4 H(3,y)-3 ) H(1,z)-162 H(2,y)-108 H(0,2,y)+108 H(1,0,y)-108 H(2,0,y)+216 H(3,2,y)+108 H(0,1,z)+54 [\_2]{}+22 )-9 H(0,y) (6-24 H(1,z))-9 (6-24 H(2,y)) H(0,z) ) } + n\_f\^2 { H(0,y)+ H(0,z)- }

[B]{}\^[(2)]{}\_[2 ; 0]{} = C\_A\^2 { \_4 ( ) + \_3 ( -H(0,y)-6 H(1,y)+11 H(2,y)-H(0,z)+5 H(1,z)- ( 25 s\^4 +2 s\^3 (79 t\^2+266 t u+79 u\^2)+s\^2 (133 t\^3+1097 t\^2 u+1097 t u\^2+133 u\^3)+2 s t u (241 t\^2+806 t u+241 u\^2)+457 t\^2 u\^2 ) / ( 36 \^2 \^2 ) ) + \_2 ( ( H(0,z) ( s\^2 (109 t+49 u)+2 s t (115 t+61 u)+t\^2 (121 t+85 u))) / ( 12 \^2 ) + ( H(0,y) (s\^2 (49 t+109 u)+2 s u (61 t+115 u)+u\^2 (85 t+121 u)) ) / ( 12 \^2 ) +( H(2,y) (-144 s\^5-257 s\^4 -2 s\^3 (35 t\^2+202 t u+35 u\^2)+s\^2 (43 t\^3-25 t\^2 u-25 t u\^2+43 u\^3)+2 s t u (55 t\^2+74 t u+55 u\^2)+79 t\^2 u\^2) ) / ( 12 \^2 \^2 ) +( H(1,z) (-72 s\^5-s\^4 (113 t+41 u)+2 s\^3 (t\^2+14 t u+73 u\^2)+s\^2 (43 t\^3+191 t\^2 u+407 t u\^2+115 u\^3)+2 s t u (55 t\^2+182 t u+127 u\^2)+t\^2 u\^2 (79 t+151 u)) ) / ( 12 \^2 \^2 ) ++3 H(0,y) H(0,z)+2 H(1,y) H(0,z)+3 H(0,y) H(1,z)+2 H(1,y) H(1,z)+H(2,y) H(0,z)-2 H(3,y) H(1,z)+3 H(0,2,y)+2 H(1,2,y)+H(2,0,y)-2 H(3,2,y)+H(0,1,z)+H(1,0,z)+2 H(1,1,z) -( 28 s\^3 +s\^2 (28 t\^2+71 t u+28 u\^2)+34 s t u-3 t\^2 u\^2 ) / ( 9 s ) )

[=]{} + ( - + (-17) H(0,z) H(0,y) -( (179 u\^2+161 s u+107 t u+107 s t) H(1,z) H(0,y) ) / ( 18 ) +7 H(0,0,z) H(0,y) -( ((13 t+25 u) s\^2+2 u (19 t+28 u) s+31 u\^2) H(0,1,z) H(0,y)) / ( 6 \^2 ) +( (17 t+23 u) H(1,0,z) H(0,y) ) / ( 2 ) - H(1,1,z) H(0,y)-2 H(0,0,1,z) H(0,y)-2 H(0,1,1,z) H(0,y)+2 H(1,0,0,z) H(0,y)+2 H(1,1,0,z) H(0,y) - - H(1,z) -( (179 t\^2+161 s t+107 u t+107 s u) H(0,z) H(2,y) ) / ( 18 ) - H(1,z) H(2,y) - H(2,y) +( (188 s\^2+197 s+206 t u) H(1,z) H(3,y) ) / (9 ) + 7 H(0,z) H(0,0,y)+7 H(1,z) H(0,0,y)

[=]{} + H(0,0,y)+7 H(2,y) H(0,0,z)+2 H(0,0,y) H(0,0,z)+ H(0,0,z) +( ((215 t+269 u) s\^2+(215 t\^2+520 u t+287 u\^2) s+t u (233 t+305 u)) H(0,1,z)) / ( 18 ) +( (18 s\^5+59 s\^4+(61 t\^2+152 u t+61 u\^2) s\^3+4 (5 t\^3+31 u t\^2+31 u\^2 t+5 u\^3) s\^2+34 t \^2 u s+11 t\^2 u\^2) H(2,y) H(0,1,z) ) / ( 3 \^2 \^2 ) + H(3,y) H(0,1,z)+2 H(0,0,y) H(0,1,z) -( (179 u\^2+161 s u+107 t u+107 s t) H(0,2,y)) / ( 18 ) + ( ((17 t+29 u) s\^2+2 t (14 t+23 u) s+11 t\^2 ) H(0,z) H(0,2,y) ) / ( 6 \^2 ) - H(1,z) H(0,2,y)+2 H(0,0,z) H(0,2,y)-2 H(0,1,z) H(0,2,y) -( ((5 t+9 u) s\^4+(9 t\^2+28 u t+19 u\^2) s\^3+2 (2 t\^3+14 u t\^2+21 u\^2 t+5 u\^3) s\^2+2 t u (5 t\^2+14 u t+9 u\^2) s+7 t\^2 u\^2) H(1,z) H(0,3,y) ) / ( \^2 \^2 ) -2 H(0,1,z) H(0,3,y)+( u (6 s\^2-2 t s+9 u s+t u) H(1,0,y)) / ( 3 s ) - - ( 3 (s+t+u) H(1,z) H(1,0,y) ) / ( ) -( ((125 t+107 u) s\^2+t (125 t+119 u) s-6 t\^2 u) H(1,0,z) ) / ( 18 s ) +( (-3 s\^3-(t-3 u) s\^2+8 t s+6 t\^2 ) H(2,y) H(1,0,z) ) / ( \^2 ) -7 H(3,y) H(1,0,z)+2 H(0,0,y) H(1,0,z)-6 H(0,3,y) H(1,0,z)

[=]{} + H(3,y) H(1,1,z)+2 H(0,0,y) H(1,1,z)- H(1,1,z)-2 H(0,1,z) H(1,2,y)+2 H(1,0,z) H(1,2,y) -( (179 u\^2+161 s u+107 t u+107 s t) H(2,0,y) ) / ( 18 ) + H(0,z) H(2,0,y)- H(1,z) H(2,0,y)+2 H(0,0,z) H(2,0,y)+2 H(0,1,z) H(2,0,y)+2 H(1,0,z) H(2,0,y)- H(0,z) H(2,2,y)+2 H(0,0,z) H(2,2,y)+4 H(0,1,z) H(2,2,y)+4 H(1,0,z) H(2,2,y)- H(2,2,y) +( (18 s\^5+61 s\^4+5 (13 t\^2+32 u t+13 u\^2) s\^3+2 (11 t\^3+67 u t\^2+67 u\^2 t+11 u\^3) s\^2+38 t \^2 u s+13 t\^2 u\^2) H(1,z) H(2,3,y) ) / ( 3 \^2 \^2 )-6 H(0,1,z) H(2,3,y)+2 H(1,0,z) H(2,3,y)-7 H(1,z) H(3,0,y)+2 H(0,1,z) H(3,0,y)-2 H(1,0,z) H(3,0,y)+ ( (188 s\^2+197 s+206 t u) H(3,2,y) ) / ( 9 ) -7 H(0,z) H(3,2,y)+ H(1,z) H(3,2,y)+ H(1,z) H(3,3,y)-12 H(0,1,z) H(3,3,y)+4 H(1,0,z) H(3,3,y)+ ( ((7 u-5 t) s\^4+(-13 t\^2+4 u t+17 u\^2) s\^3-2 (4 t\^3+8 u t\^2-13 u\^2 t-5 u\^3) s\^2+2 t u (-5 t\^2+2 u t+7 u\^2) s+t\^2 u\^2) H(0,0,1,z) ) / ( 3 \^2 \^2 ) - 2 H(1,y) H(0,0,1,z)-2 H(2,y) H(0,0,1,z)-8 H(3,y) H(0,0,1,z)+2 H(0,z) H(0,0,2,y)+2 H(1,z) H(0,0,2,y)+7 H(0,0,2,y) - ( u (s\^2+(u-2 t) s-3 t u) H(0,1,0,y) ) / ( \^2 ) - ( t (s\^2+(t-2 u) s-3 t u) H(0,1,0,z)) / ( \^2 ) +2 H(1,y) H(0,1,0,z)+6 H(2,y) H(0,1,0,z)+4 H(3,y) H(0,1,0,z)+ H(0,1,1,z)+2 H(0,z) H(0,2,0,y)+2 H(1,z) H(0,2,0,y)+7 H(0,2,0,y)-2 H(0,z) H(0,2,2,y)- H(0,2,2,y)+4 H(1,z) H(0,3,0,y)- ( ((5 t+9 u) s\^4+(9 t\^2+28 u t+19 u\^2) s\^3+2 (2 t\^3+14 u t\^2+21 u\^2 t+5 u\^3) s\^2+2 t u (5 t\^2+14 u t+9 u\^2) s+7 t\^2 u\^2) H(0,3,2,y) ) / ( \^2 \^2 ) - 6 H(0,z) H(0,3,2,y)+4 H(1,z) H(0,3,3,y)+2 H(2,y) H(1,0,0,z)

[=]{} +7 H(1,0,0,z) +( (9 s\^5+32 s\^4+(34 t\^2+80 u t+34 u\^2) s\^3+(11 t\^3+61 u t\^2+61 u\^2 t+11 u\^3) s\^2+t u (16 t\^2+23 u t+16 u\^2) s+2 t\^2 u\^2) H(1,0,1,z) ) / ( 3 \^2 \^2 ) -2 H(1,y) H(1,0,1,z)+6 H(2,y) H(1,0,1,z) - ( 3 (s+t+u) H(1,0,2,y) ) / ( ) -2 H(0,z) H(1,0,2,y)+2 H(1,z) H(1,0,3,y)+ ( 6 H(1,1,0,y) ) / ( ) + ( (2 (t+3 u) s\^2+t (5 t+14 u) s+3 t\^2 (t+3 u)) H(1,1,0,z) ) / ( \^2 ) +2 H(1,y) H(1,1,0,z)+6 H(2,y) H(1,1,0,z)--2 H(1,z) H(1,2,3,y)+2 H(0,z) H(2,0,0,y)+2 H(1,z) H(2,0,0,y)+7 H(2,0,0,y)- H(2,0,2,y)+2 H(1,z) H(2,0,3,y)+ ( (-9 s\^3+(11 t-u) s\^2+28 u s+20 u\^2) H(2,1,0,y) ) / ( 3 \^2 ) +2 H(1,z) H(2,1,0,y)+2 H(0,z) H(2,2,0,y)- H(2,2,0,y)+4 H(1,z) H(2,2,3,y)+2 H(1,z) H(2,3,0,y)+ ( (18 s\^5+61 s\^4+5 (13 t\^2+32 u t+13 u\^2) s\^3+2 (11 t\^3+67 u t\^2+67 u\^2 t+11 u\^3) s\^2+38 t \^2 u s+13 t\^2 u\^2) H(2,3,2,y) ) / ( 3 \^2 \^2 ) +2 H(0,z) H(2,3,2,y)-8 H(1,z) H(2,3,3,y)-2 H(0,z) H(3,0,2,y)-7 H(3,0,2,y)+4 H(1,z) H(3,0,3,y)-2 H(0,z) H(3,2,0,y)-7 H(3,2,0,y)+ H(3,2,2,y)+4 H(1,z) H(3,3,0,y)+4 H(0,z) H(3,3,2,y)+ H(3,3,2,y)-16 H(1,z) H(3,3,3,y)+8 H(0,0,1,0,z)+2 H(0,0,1,1,z)+2 H(0,0,2,2,y)+2 H(0,1,0,1,z)+6 H(0,1,1,0,z)+2 H(0,2,0,2,y)+4 H(0,2,1,0,y)+2 H(0,2,2,0,y)+4 H(0,3,0,2,y)+4 H(0,3,2,0,y)+4 H(0,3,3,2,y)+6 H(1,0,1,0,z)+2 H(1,0,3,2,y)+2 H(1,1,0,0,z)+4 H(1,1,0,1,z)+6 H(1,1,1,0,z)-2 H(1,2,3,2,y)+2 H(2,0,0,2,y)+2 H(2,0,1,0,y)+2 H(2,0,2,0,y)+2 H(2,0,3,2,y)+2 H(2,1,0,2,y)+2 H(2,1,2,0,y)+2 H(2,2,0,0,y)+4 H(2,2,1,0,y)+4 H(2,2,3,2,y)+2 H(2,3,0,2,y)+2 H(2,3,2,0,y)-8 H(2,3,3,2,y)+4 H(3,0,3,2,y)+4 H(3,3,0,2,y)+4 H(3,3,2,0,y)-16 H(3,3,3,2,y)+ ) }

[=]{} + C\_A C\_F { \_4 ( ) + \_3 ( 22 H(1,y)-18 H(2,y)+4 H(1,z)+ ( 803 s\^5 +4 s\^4 (424 t\^2+965 t u+424 u\^2)+s\^3 (821 t\^3+5329 t\^2 u+5329 t u\^2+821 u\^3)+s\^2 (-72 t\^4+1966 t\^3 u+5336 t\^2 u\^2+1966 t u\^3-72 u\^4)+s t u (-144 t\^3+1307 t\^2 u+1307 t u\^2-144 u\^3)-72 t\^2 u\^2 (t\^2+u\^2) ) / ( 18 s \^2 \^2 ) ) + \_2 ( -( 2 H(1,y) (18 s\^2+23 s t+11 s u-6 t\^2) ) / ( 3 s ) +( t H(0,z) (-10 s\^3-3 s\^2 (7 t+2 u)-9 s t +2 t\^3) ) / ( s \^2 ) +( u H(0,y) (-10 s\^3-3 s\^2 (2 t+7 u)-9 s u+2 u\^3) ) / ( s \^2 ) +( H(2,y) (144 s\^5 +s\^4 (289 t\^2+554 t u+289 u\^2)+4 s\^3 (32 t\^3+147 t\^2 u+147 t u\^2+32 u\^3)+s\^2 (-17 t\^4+90 t\^3 u+166 t\^2 u\^2+90 t u\^3-17 u\^4)-2 s t u (35 t\^3+111 t\^2 u+111 t u\^2+35 u\^3)-t\^2 u\^2 (71 t\^2+166 t u+71 u\^2)) ) / ( 6 \^2 \^2 \^2 ) -( H(1,z) (-72 s\^6 \^2-s\^5 (101 t\^3+231 t\^2 u+183 t u\^2+53 u\^3)+4 s\^4 \^2 (8 t\^2+25 t u+26 u\^2)+s\^3 (61 t\^5+427 t\^4 u+1120 t\^3 u\^2+1216 t\^2 u\^3+523 t u\^4+61 u\^5)+2 s\^2 \^2 u (79 t\^3+236 t\^2 u+79 t u\^2-12 u\^3)+s t u\^2 (115 t\^4+369 t\^3 u+321 t\^2 u\^2+19 t u\^3-48 u\^4)-24 t\^2 \^2 u\^4) ) / ( 6 s \^2 \^2 \^3 ) -2 H(1,y) H(0,z)-6 H(0,y) H(1,z)-2 H(1,y) H(1,z)-4 H(2,y) H(0,z)+2 H(0,2,y)+6 H(1,0,y)-2 H(1,2,y)-4 H(2,0,y)-8 H(2,1,y)+8 H(2,2,y)+2 H(0,1,z)+2 H(1,0,z)-2 H(1,1,z) +( 19 s\^3 +s\^2 (19 t\^2+182 t u+19 u\^2)+145 s t u+108 t\^2 u\^2 ) / ( 6 s ) )

[=]{} + ( ( (26 t\^2+26 s t+33 u t+15 s u) H(0,y) ) / ( 3 ) + (( 18 t u ) / ( s ) +3) H(0,z) H(0,y)+ ( (242 u\^2+215 s u+134 t u+134 s t) H(1,z) H(0,y) ) / ( 9 ) + ( ((13 t+31 u) s\^2+u (44 t+71 u) s+40 u\^2) H(0,1,z) H(0,y) ) / ( 3 \^2 ) - ( (19 s t\^2+62 s u t+31 s u\^2-6 u\^2) H(1,0,z) H(0,y) ) / ( 3 s \^2 ) - H(1,1,z) H(0,y)+4 H(0,0,1,z) H(0,y)+2 H(0,1,0,z) H(0,y)+8 H(0,1,1,z) H(0,y)+2 H(1,0,1,z) H(0,y)-4 H(1,1,0,z) H(0,y)+ ( (26 u\^2+26 s u+33 t u+15 s t) H(0,z) ) / ( 3 ) + H(1,z) + ( (242 t\^2+215 s t+134 u t+134 s u) H(0,z) H(2,y) ) / ( 9 ) +2 H(1,z) H(2,y)+ H(2,y) - ( (466 s\^2+493 s+520 t u) H(1,z) H(3,y) ) / ( 9 ) - 14 H(1,z) H(0,0,y)-14 H(2,y) H(0,0,z)- ( ((251 t+332 u) s\^2+(251 t\^2+637 u t+359 u\^2) s+2 t u (139 t+193 u)) H(0,1,z) ) / ( 9 ) - ( (36 s\^5+112 s\^4+(107 t\^2+268 u t+107 u\^2) s\^3+(31 t\^3+191 u t\^2+191 u\^2 t+31 u\^3) s\^2+2 t u (22 t\^2+35 u t+22 u\^2) s+4 t\^2 u\^2) H(2,y) H(0,1,z) ) / ( 3 \^2 \^2 ) - H(3,y) H(0,1,z)-4 H(0,0,y) H(0,1,z) + ( (242 u\^2+215 s u+134 t u+134 s t) H(0,2,y) ) / ( 9 ) - ( ((29 t+47 u) s\^2+t (49 t+76 u) s+20 t\^2 ) H(0,z) H(0,2,y) ) / ( 3 \^2 ) - H(1,z) H(0,2,y)-4 H(0,0,z) H(0,2,y)+6 H(0,1,z) H(0,2,y)+ ( (2 (7 t+13 u) s\^4+5 (5 t\^2+16 u t+11 u\^2) s\^3+(11 t\^3+79 u t\^2+121 u\^2 t+29 u\^3) s\^2+4 t u (7 t\^2+20 u t+13 u\^2) s+20 t\^2 u\^2) H(1,z) H(0,3,y) ) / ( \^2 \^2 )

[=]{} + 4 H(0,1,z) H(0,3,y)+ ( (-2 (76 t+103 u) s\^2+(37 t-233 u) u s+162 t u\^2) H(1,0,y) ) / ( 9 s ) + ( 2 (3 t\^2+8 s t+14 s u) H(0,z) H(1,0,y) ) / ( 3 s ) + ( + 9) H(1,z) H(1,0,y)-2 H(0,1,z) H(1,0,y)+ ( ((t-2 u) s\^2+t (t+19 u) s+18 t\^2 u) H(1,0,z) ) / ( s ) + ( (18 s\^3+(31 t\^2+32 u t+13 u\^2) s\^2-t (t\^2+17 u t-8 u\^2) s-2 t\^2 (7 t\^2+20 u t+7 u\^2)) H(2,y) H(1,0,z) ) / ( 3 \^2 \^2 ) +20 H(3,y) H(1,0,z)-4 H(0,0,y) H(1,0,z)+8 H(0,2,y) H(1,0,z)+16 H(0,3,y) H(1,0,z)+2 H(1,0,y) H(1,0,z)+ H(3,y) H(1,1,z)-8 H(0,0,y) H(1,1,z)+4 H(0,3,y) H(1,1,z)+2 H(1,1,z)+6 H(0,1,z) H(1,2,y)-2 H(1,0,z) H(1,2,y)+ ( (242 u\^2+215 s u+134 t u+134 s t) H(2,0,y) ) / ( 9 ) - ( (47 t\^3+153 u t\^2+153 u\^2 t+47 u\^3) H(0,z) H(2,0,y) ) / ( 3 \^3 ) - H(1,z) H(2,0,y)-4 H(0,0,z) H(2,0,y)-2 H(0,1,z) H(2,0,y)-6 H(1,0,z) H(2,0,y)- H(0,z) H(2,2,y)-8 H(0,0,z) H(2,2,y)-16 H(0,1,z) H(2,2,y)-8 H(1,0,z) H(2,2,y)+2 H(2,2,y)-( (36 s\^5+98 s\^4+(79 t\^2+212 u t+79 u\^2) s\^3+(17 t\^3+121 u t\^2+121 u\^2 t+17 u\^3) s\^2+2 t u (8 t\^2+7 u t+8 u\^2) s-10 t\^2 u\^2) H(1,z) H(2,3,y) ) / ( 3 \^2 \^2 ) +12 H(0,1,z) H(2,3,y)-4 H(1,0,z) H(2,3,y)+20 H(1,z) H(3,0,y)-4 H(0,1,z) H(3,0,y)+4 H(1,1,z) H(3,0,y)- ( (466 s\^2+493 s+520 t u) H(3,2,y) ) / ( 9 ) +20 H(0,z) H(3,2,y)+ H(1,z) H(3,2,y)-4 H(0,1,z) H(3,2,y)-4 H(1,0,z) H(3,2,y)- H(1,z) H(3,3,y)+28 H(0,1,z) H(3,3,y)-12 H(1,0,z) H(3,3,y)-8 H(1,1,z) H(3,3,y)+ ( (4 (4 t-5 u) s\^4+(41 t\^2-8 u t-49 u\^2) s\^3+(25 t\^3+53 u t\^2-73 u\^2 t-29 u\^3) s\^2+8 t u (4 t\^2-u t-5 u\^2) s-2 t\^2 u\^2) H(0,0,1,z) ) / ( 3 \^2 \^2 )

[=]{} +6 H(1,y) H(0,0,1,z)+16 H(3,y) H(0,0,1,z)-8 H(0,z) H(0,0,2,y)-8 H(1,z) H(0,0,2,y)-14 H(0,0,2,y)+4 H(1,z) H(0,0,3,y)+ ( (-6 u\^4-s (17 t\^2+22 u t+17 u\^2) u\^2+s\^2 (2 t\^2+19 u t-7 u\^2) u+2 s\^3 (5 t\^2+13 u t+2 u\^2)) H(0,1,0,y) ) / ( 3 s \^2 \^2 ) +2 H(0,z) H(0,1,0,y)+2 H(1,z) H(0,1,0,y)+ ( (-6 t\^4+s (t\^2+14 u t+u\^2) t\^2+s\^2 (29 t\^2+91 u t+38 u\^2) t+s\^3 (22 t\^2+62 u t+28 u\^2)) H(0,1,0,z) ) / ( 3 s \^2 \^2 ) -2 H(1,y) H(0,1,0,z)-12 H(2,y) H(0,1,0,z)-4 H(3,y) H(0,1,0,z)-4 H(3,y) H(0,1,1,z)+ H(0,1,1,z)-8 H(1,z) H(0,2,0,y)-14 H(0,2,0,y)+4 H(0,z) H(0,2,2,y)- H(0,2,2,y)+2 H(1,z) H(0,2,3,y)-8 H(1,z) H(0,3,0,y)+ ( (2 (7 t+13 u) s\^4+5 (5 t\^2+16 u t+11 u\^2) s\^3+(11 t\^3+79 u t\^2+121 u\^2 t+29 u\^3) s\^2+4 t u (7 t\^2+20 u t+13 u\^2) s+20 t\^2 u\^2) H(0,3,2,y) ) / ( \^2 \^2 ) +16 H(0,z) H(0,3,2,y)+4 H(1,z) H(0,3,2,y)-12 H(1,z) H(0,3,3,y)+4 H(0,z) H(1,0,0,y)+4 H(1,z) H(1,0,0,y)+14 H(1,0,0,y)-4 H(2,y) H(1,0,0,z)- ( (18 s\^5+49 s\^4+(35 t\^2+88 u t+35 u\^2) s\^3+4 (t\^3+5 u t\^2+5 u\^2 t+u\^3) s\^2-2 t u (5 t\^2+28 u t+5 u\^2) s-23 t\^2 u\^2) H(1,0,1,z) ) / ( 3 \^2 \^2 ) +6 H(1,y) H(1,0,1,z)-18 H(2,y) H(1,0,1,z)-4 H(3,y) H(1,0,1,z) +( + 9) H(1,0,2,y)+4 H(0,z) H(1,0,2,y)-6 H(1,z) H(1,0,3,y)- ( 2 (18 s\^2+23 t s+11 u s-6 t\^2) H(1,1,0,y) ) / ( 3 s ) - ( (2 u (3 t+u) s\^3+(3 t\^3+21 u t\^2+6 u\^2 t-4 u\^3) s\^2+t (3 t\^3+18 u t\^2+3 u\^2 t-8 u\^3) s-4 t\^2 u\^2) H(1,1,0,z) ) / ( s \^2 \^2 ) -2 H(1,y) H(1,1,0,z)-18 H(2,y) H(1,1,0,z)-4 H(3,y) H(1,1,0,z) + (+9) H(1,2,0,y)+2 H(0,z) H(1,2,0,y)

[=]{} +6 H(1,z) H(1,2,3,y)-4 H(0,z) H(2,0,0,y)-8 H(1,z) H(2,0,0,y)-14 H(2,0,0,y)+2 H(0,z) H(2,0,2,y)- H(2,0,2,y)-4 H(1,z) H(2,0,3,y)+ ( (6 s\^3+(9 t\^2+20 u t+15 u\^2) s\^2+u (12 t\^2+13 u t+9 u\^2) s-4 t u\^3) H(2,1,0,y) ) / ( \^2 \^2 ) -2 H(0,z) H(2,1,0,y)-2 H(1,z) H(2,1,0,y)-4 H(0,z) H(2,2,0,y)- H(2,2,0,y)-16 H(1,z) H(2,2,3,y)-4 H(1,z) H(2,3,0,y)- ( (36 s\^5+98 s\^4+(79 t\^2+212 u t+79 u\^2) s\^3+(17 t\^3+121 u t\^2+121 u\^2 t+17 u\^3) s\^2+2 t u (8 t\^2+7 u t+8 u\^2) s-10 t\^2 u\^2) H(2,3,2,y) ) / ( 3 \^2 \^2 ) -4 H(0,z) H(2,3,2,y)+16 H(1,z) H(2,3,3,y)+8 H(0,z) H(3,0,2,y)+4 H(1,z) H(3,0,2,y)+20 H(3,0,2,y)-12 H(1,z) H(3,0,3,y)+4 H(1,z) H(3,2,0,y)+20 H(3,2,0,y)+4 H(0,z) H(3,2,2,y)+ H(3,2,2,y)-8 H(1,z) H(3,2,3,y)-12 H(1,z) H(3,3,0,y)-12 H(0,z) H(3,3,2,y)-8 H(1,z) H(3,3,2,y)- H(3,3,2,y)+40 H(1,z) H(3,3,3,y)+4 H(0,0,1,0,y)-8 H(0,0,1,0,z)-8 H(0,0,1,1,z)-8 H(0,0,2,2,y)+4 H(0,0,3,2,y)-8 H(0,1,0,1,z)+2 H(0,1,0,2,y)-12 H(0,1,1,0,z)+2 H(0,1,2,0,y)-8 H(0,2,0,2,y)-2 H(0,2,1,0,y)-8 H(0,2,2,0,y)+2 H(0,2,3,2,y)-8 H(0,3,0,2,y)-8 H(0,3,2,0,y)+4 H(0,3,2,2,y)-12 H(0,3,3,2,y)-2 H(1,0,0,1,z)+4 H(1,0,0,2,y)+4 H(1,0,1,0,y)-4 H(1,0,1,0,z)+4 H(1,0,2,0,y)-6 H(1,0,3,2,y)-12 H(1,1,0,1,z)-14 H(1,1,1,0,z)+4 H(1,2,0,0,y)+4 H(1,2,1,0,y)+6 H(1,2,3,2,y)-8 H(2,0,0,2,y)-8 H(2,0,2,0,y)-4 H(2,0,3,2,y)+4 H(2,1,0,0,y)-2 H(2,1,0,2,y)-8 H(2,1,1,0,y)-2 H(2,1,2,0,y)-8 H(2,2,0,0,y)-8 H(2,2,1,0,y)-16 H(2,2,3,2,y)-4 H(2,3,0,2,y)-4 H(2,3,2,0,y)+16 H(2,3,3,2,y)+8 H(3,0,1,0,y)+4 H(3,0,2,2,y)-12 H(3,0,3,2,y)+4 H(3,2,0,2,y)-8 H(3,2,1,0,y)+4 H(3,2,2,0,y)-8 H(3,2,3,2,y)-12 H(3,3,0,2,y)-12 H(3,3,2,0,y)-8 H(3,3,2,2,y)+40 H(3,3,3,2,y)- ) }

[=]{} + C\_F\^2 { \_4 ( -22 ) + \_3 ( -16 H(1,y)+8 H(2,y)-8 H(1,z)- ( 2 (25 s\^5 +s\^4 (51 t\^2+112 t u+51 u\^2)+4 s\^3 (6 t\^3+37 t\^2 u+37 t u\^2+6 u\^3)-2 s\^2 (t\^4-27 t\^3 u-69 t\^2 u\^2-27 t u\^3+u\^4)+s t u (-4 t\^3+33 t\^2 u+33 t u\^2-4 u\^3)-2 t\^2 u\^2 (t\^2+u\^2)) ) / ( s \^2 \^2 ) ) + \_2 ( ( 2 t H(0,z) (4 s\^3+2 s\^2 (4 t+u)+3 s t -t\^3) ) / ( s \^2 ) + ( 2 u H(0,y) (4 s\^3+2 s\^2 (t+4 u)+3 s u-u\^3) ) / ( s \^2 ) + ( 2 H(2,y) (s\^4 (5 t\^2+12 t u+5 u\^2)+s\^3 (11 t\^3+39 t\^2 u+39 t u\^2+11 u\^3)+s\^2 (6 t\^4+42 t\^3 u+76 t\^2 u\^2+42 t u\^3+6 u\^4)+2 s t u (7 t\^3+26 t\^2 u+26 t u\^2+7 u\^3)+t\^2 u\^2 (9 t\^2+20 t u+9 u\^2)) ) / ( \^2 \^2 \^2 ) + ( 2 H(1,z) (s\^5 (5 t\^2+16 t u+9 u\^2)+s\^4 (11 t\^3+47 t\^2 u+53 t u\^2+17 u\^3)+s\^3 (6 t\^4+46 t\^3 u+92 t\^2 u\^2+54 t u\^3+6 u\^4)+2 s\^2 u (7 t\^4+29 t\^3 u+29 t\^2 u\^2+6 t u\^3-u\^4)+s t u\^2 (9 t\^3+20 t\^2 u+5 t u\^2-4 u\^3)-2 t\^2 u\^4) ) / ( s \^2 \^2 \^2 ) + ( 4 t (2 s-t) H(1,y) ) / ( s ) -8 H(2,y) H(1,z)+8 H(0,1,y)-8 H(0,2,y)+8 H(1,1,y)+8 H(2,1,y)-16 H(2,2,y)+ ( 2 (6 s\^3 +s\^2 (6 t\^2+4 t u+6 u\^2)-s t u-6 t\^2 u\^2) ) / ( s ) )

[=]{} + ( - ( 12 t u H(0,y) H(0,z) ) / ( s ) - ( 2 (3 s (2 t+u)+t (7 t+3 u)) H(2,y) H(0,z) ) / ( ) + ( 2 ((4 t+6 u) s\^2+t (7 t+10 u) s+3 t\^2 ) H(0,2,y) H(0,z) ) / ( \^2 ) - ( 2 (t\^2+s (t+3 u)) H(1,0,y) H(0,z) ) / ( s ) + ( (6 t\^3+22 u t\^2+22 u\^2 t+6 u\^3) H(2,0,y) H(0,z) ) / ( \^3 ) +12 H(2,2,y) H(0,z)-12 H(3,2,y) H(0,z)+8 H(0,0,2,y) H(0,z)-8 H(0,3,2,y) H(0,z)-4 H(2,0,2,y) H(0,z)-8 H(3,0,2,y) H(0,z)-8 H(3,2,2,y) H(0,z)+8 H(3,3,2,y) H(0,z)- ( 2 (3 s (t+2 u)+u (3 t+7 u)) H(0,y) H(1,z) ) / ( ) -9 H(1,z)+9 H(1,z) H(2,y)-9 H(2,y)+ ( (20 s\^2+22 s+24 t u) H(1,z) H(3,y) ) / ( ) + ( 2 ((4 t+7 u) s\^2+(4 t\^2+13 u t+8 u\^2) s+t u (5 t+9 u)) H(0,1,z) ) / ( ) - ( 2 u (+u) (2 s+3 u) H(0,y) H(0,1,z) ) / ( \^2 ) - ( 2 (2 s\^4+(5 t\^2+12 u t+5 u\^2) s\^3+(3 t\^3+19 u t\^2+19 u\^2 t+3 u\^3) s\^2+2 t u (4 t\^2+11 u t+4 u\^2) s+6 t\^2 u\^2) H(2,y) H(0,1,z) ) / ( \^2 \^2 ) +12 H(3,y) H(0,1,z) - ( 2 (3 s (t+2 u)+u (3 t+7 u)) H(0,2,y) ) / ( ) + 12 H(1,z) H(0,2,y)-4 H(0,1,z) H(0,2,y)- ( 2 (4 (t+2 u) s\^4+(7 t\^2+24 u t+17 u\^2) s\^3+(3 t\^3+23 u t\^2+37 u\^2 t+9 u\^3) s\^2+8 t u (t\^2+3 u t+2 u\^2) s+6 t\^2 u\^2) H(1,z) H(0,3,y) ) / ( \^2 \^2 ) + ( 2 ((3 t+5 u) s\^2+(6 u\^2-4 t u) s-6 t u\^2) H(1,0,y) ) / ( s ) -6 H(1,z) H(1,0,y)+4 H(0,1,z) H(1,0,y)

[=]{} - ( 2 t (s\^2+(t+7 u) s+6 t u) H(1,0,z) ) / ( s ) + ( 2 u (2 s (2 t+u)- u) H(0,y) H(1,0,z) ) / ( s \^2 ) + ( 2 t (2 (t+2 u) s\^2+(5 t\^2+11 u t+2 u\^2) s+t (3 t\^2+8 u t+3 u\^2)) H(2,y) H(1,0,z) ) / ( \^2 \^2 ) -12 H(3,y) H(1,0,z)-8 H(0,2,y) H(1,0,z)-8 H(0,3,y) H(1,0,z)+12 H(0,y) H(1,1,z)-24 H(3,y) H(1,1,z)+8 H(0,0,y) H(1,1,z)-8 H(0,3,y) H(1,1,z)+9 H(1,1,z)-4 H(0,1,z) H(1,2,y)- ( 2 (3 s (t+2 u)+u (3 t+7 u)) H(2,0,y) ) / ( ) +12 H(1,z) H(2,0,y)-4 H(0,1,z) H(2,0,y)+8 H(0,0,z) H(2,2,y)+16 H(0,1,z) H(2,2,y)+9 H(2,2,y)

[=]{} - ( 2 (8 s\^4+(17 t\^2+36 u t+17 u\^2) s\^3+(9 t\^3+49 u t\^2+49 u\^2 t+9 u\^3) s\^2+2 t u (10 t\^2+23 u t+10 u\^2) s+12 t\^2 u\^2) H(1,z) H(2,3,y) ) / ( \^2 \^2 ) -12 H(1,z) H(3,0,y)-8 H(1,1,z) H(3,0,y)+ ( (20 s\^2+22 s+24 t u) H(3,2,y) ) / ( ) -24 H(1,z) H(3,2,y)+8 H(0,1,z) H(3,2,y)+24 H(1,z) H(3,3,y)-8 H(0,1,z) H(3,3,y)+8 H(1,0,z) H(3,3,y)+16 H(1,1,z) H(3,3,y)- ( 2 s (t-u) (2 s\^3+5 s\^2+(3 t\^2+10 u t+3 u\^2) s+4 t u) H(0,0,1,z) ) / ( \^2 \^2 ) - 4 H(1,y) H(0,0,1,z)+8 H(2,y) H(0,0,1,z)+8 H(1,z) H(0,0,2,y)-8 H(1,z) H(0,0,3,y)+ ( 2 ( u\^4+s (9 t\^2+16 u t+9 u\^2) u\^2+2 s\^2 (7 t\^2+12 u t+7 u\^2) u+2 s\^3 (3 t\^2+5 u t+3 u\^2)) H(0,1,0,y) ) / ( s \^2 \^2 ) -4 H(1,z) H(0,1,0,y)+ ( 2 t (-2 u s\^3+2 (t\^2+u\^2) s\^2+t (3 t\^2+4 u t+3 u\^2) s+t\^3 ) H(0,1,0,z) ) / ( s \^2 \^2 ) +4 H(2,y) H(0,1,0,z)-8 H(0,y) H(0,1,1,z)+8 H(3,y) H(0,1,1,z)-12 H(0,1,1,z)+8 H(1,z) H(0,2,0,y)+12 H(0,2,2,y)-4 H(1,z) H(0,2,3,y)- ( 2 (4 (t+2 u) s\^4+(7 t\^2+24 u t+17 u\^2) s\^3+(3 t\^3+23 u t\^2+37 u\^2 t+9 u\^3) s\^2+8 t u (t\^2+3 u t+2 u\^2) s+6 t\^2 u\^2) H(0,3,2,y) ) / ( \^2 \^2 ) -8 H(1,z) H(0,3,2,y)+8 H(1,z) H(0,3,3,y)-8 H(1,z) H(1,0,0,y)- ( 2 (5 s\^4+(11 t\^2+24 u t+11 u\^2) s\^3+(6 t\^3+34 u t\^2+34 u\^2 t+6 u\^3) s\^2+2 t u (7 t\^2+17 u t+7 u\^2) s+9 t\^2 u\^2) H(1,0,1,z) ) / ( \^2 \^2 ) -4 H(0,y) H(1,0,1,z)-4 H(1,y) H(1,0,1,z)

[=]{} +12 H(2,y) H(1,0,1,z)+8 H(3,y) H(1,0,1,z)-6 H(1,0,2,y)+4 H(1,z) H(1,0,3,y)+ ( 4 (2 s-t) t H(1,1,0,y) ) / ( s ) - ( 2 ((t\^2-2 u t-u\^2) s\^3+(t\^3-7 u t\^2-2 u\^2 t+2 u\^3) s\^2+(4 t u\^3-6 t\^3 u) s+2 t\^2 u\^2) H(1,1,0,z) ) / ( s \^2 \^2 ) +4 H(2,y) H(1,1,0,z)-6 H(1,2,0,y)-4 H(1,z) H(1,2,3,y)+8 H(1,z) H(2,0,0,y)+12 H(2,0,2,y)- ( 2 ((6 t\^2+8 u t+4 u\^2) s\^2+u (10 t\^2+13 u t+7 u\^2) s+u\^2 (3 t\^2+4 u t+3 u\^2)) H(2,1,0,y) ) / ( \^2 \^2 ) -4 H(1,z) H(2,1,0,y)+12 H(2,2,0,y)+16 H(1,z) H(2,2,3,y)- ( 2 (8 s\^4+(17 t\^2+36 u t+17 u\^2) s\^3+(9 t\^3+49 u t\^2+49 u\^2 t+9 u\^3) s\^2+2 t u (10 t\^2+23 u t+10 u\^2) s+12 t\^2 u\^2) H(2,3,2,y) ) / ( \^2 \^2 ) -8 H(1,z) H(3,0,2,y)-12 H(3,0,2,y)+8 H(1,z) H(3,0,3,y)-8 H(1,z) H(3,2,0,y)-12 H(3,2,0,y)-24 H(3,2,2,y)+16 H(1,z) H(3,2,3,y)+8 H(1,z) H(3,3,0,y)+16 H(1,z) H(3,3,2,y)+24 H(3,3,2,y)-16 H(1,z) H(3,3,3,y)+8 H(0,0,1,1,z)+8 H(0,0,2,2,y)-8 H(0,0,3,2,y)+8 H(0,1,0,1,z)-4 H(0,1,0,2,y)+8 H(0,1,1,0,y)+8 H(0,1,1,0,z)-4 H(0,1,2,0,y)+8 H(0,2,0,2,y)-4 H(0,2,1,0,y)+8 H(0,2,2,0,y)-4 H(0,2,3,2,y)-8 H(0,3,2,2,y)+8 H(0,3,3,2,y)+4 H(1,0,0,1,z)-8 H(1,0,0,2,y)+4 H(1,0,1,0,y)-8 H(1,0,2,0,y)+4 H(1,0,3,2,y)+8 H(1,1,0,0,y)+8 H(1,1,0,1,z)+8 H(1,1,1,0,y)+8 H(1,1,1,0,z)-8 H(1,2,0,0,y)-4 H(1,2,1,0,y)-4 H(1,2,3,2,y)+8 H(2,0,0,2,y)-4 H(2,0,1,0,y)+8 H(2,0,2,0,y)-8 H(2,1,0,0,y)-4 H(2,1,0,2,y)+8 H(2,1,1,0,y)-4 H(2,1,2,0,y)+8 H(2,2,0,0,y)+16 H(2,2,3,2,y)-8 H(3,0,1,0,y)-8 H(3,0,2,2,y)+8 H(3,0,3,2,y)-8 H(3,2,0,2,y)+8 H(3,2,1,0,y)-8 H(3,2,2,0,y)+16 H(3,2,3,2,y)+8 H(3,3,0,2,y)+8 H(3,3,2,0,y)+16 H(3,3,2,2,y)-16 H(3,3,3,2,y)+ ) }

[=]{} + C\_A n\_f { \_3 ( - ) + \_2 ( - H(0,y)+ H(2,y)- H(0,z)+ H(1,z)-- ) + ( H(0,y) H(0,z) (-)-+H(1,0,z) (-)+++ H(0,y) H(1,z)-H(0,y) H(0,0,z)+ H(0,y) H(0,1,z)-H(0,y) H(1,0,z)+ H(0,y) H(1,1,z)+ H(2,y) H(0,z)+H(2,y) H(1,z)- H(3,y) H(1,z)-H(0,0,y) H(0,z)-H(0,0,y) H(1,z)-H(2,y) H(0,0,z)- H(2,y) H(0,1,z)- H(3,y) H(0,1,z)- H(0,2,y) H(0,z)+ H(0,2,y) H(1,z)+H(0,3,y) H(1,z)+H(3,y) H(1,0,z)- H(3,y) H(1,1,z)-H(2,0,y) H(0,z)+ H(2,0,y) H(1,z)+ H(2,2,y) H(0,z)- H(2,3,y) H(1,z)+H(3,0,y) H(1,z)+H(3,2,y) H(0,z)- H(3,2,y) H(1,z)- H(3,3,y) H(1,z)+ H(2,y)- H(0,0,y)+ H(0,2,y)+ H(2,0,y)+H(2,2,y)- H(3,2,y)-H(0,0,2,y)-H(0,2,0,y)+ H(0,2,2,y)+H(0,3,2,y)-H(2,0,0,y)+ H(2,0,2,y)- H(2,1,0,y)+ H(2,2,0,y)- H(2,3,2,y)+H(3,0,2,y)+H(3,2,0,y)- H(3,2,2,y)- H(3,3,2,y)+ H(1,z)- H(0,0,z)- H(0,1,z)+H(1,1,z)- H(0,0,1,z)- H(0,1,1,z)-H(1,0,0,z)- H(1,0,1,z)- ) }

[=]{} + C\_F n\_f { \_3 ( - ) + \_2 ( H(1,y)- H(2,y)- H(1,z)+ ) + ( --- H(0,y) H(1,z)- H(0,y) H(0,1,z)+ H(0,y) H(1,0,z)- H(0,y) H(1,1,z)- H(2,y) H(0,z)-2 H(2,y) H(1,z)+ H(3,y) H(1,z)+2 H(0,0,y) H(1,z)+2 H(2,y) H(0,0,z)+ H(2,y) H(0,1,z)+ H(3,y) H(0,1,z)+ H(0,2,y) H(0,z)- H(0,2,y) H(1,z)-2 H(0,3,y) H(1,z)- H(1,0,y) H(0,z)- H(2,y) H(1,0,z)-2 H(3,y) H(1,0,z)+ H(3,y) H(1,1,z)+ H(2,0,y) H(0,z)- H(2,0,y) H(1,z)- H(2,2,y) H(0,z)+ H(2,3,y) H(1,z)-2 H(3,0,y) H(1,z)-2 H(3,2,y) H(0,z)+ H(3,2,y) H(1,z)+ H(3,3,y) H(1,z)- H(2,y)- H(0,2,y)+ H(1,0,y)- H(2,0,y)-2 H(2,2,y)+ H(3,2,y)+2 H(0,0,2,y)- H(0,1,0,y)+2 H(0,2,0,y)- H(0,2,2,y)-2 H(0,3,2,y)-2 H(1,0,0,y)+ H(1,1,0,y)+2 H(2,0,0,y)- H(2,0,2,y)- H(2,2,0,y)+ H(2,3,2,y)-2 H(3,0,2,y)-2 H(3,2,0,y)+ H(3,2,2,y)+ H(3,3,2,y)- H(1,z)+ H(0,1,z)+ H(1,0,z)-2 H(1,1,z)+ H(0,0,1,z)- H(0,1,0,z)+ H(0,1,1,z)+ H(1,0,1,z)+ ) }

[=]{} + n\_f\^2 { H(0,y) H(0,z)- H(0,y)+ H(0,0,y)- H(0,z)+ H(0,0,z)+ }

[99]{}

ATLAS Collaboration, ATLAS-CONF-2011-163.

CMS Collaboration, CMS-PAS-HIG-11-032.

R. Barate [*et al.*]{} \[ LEP Working Group for Higgs boson searches and ALEPH and DELPHI and L3 and OPAL Collaborations \], Phys. Lett.  [**B565**]{} (2003) 61.

CDF and D0 Collaborations, Phys.  Rev.  Lett. [**104**]{} (2010) 061802.

ALEPH,CDF,D0, DELPHI,L3,OPAL,SLD Collaborations, CERN PH-EP-2010-095,(2010).

S. Dawson, Nucl. Phys. B [**359**]{} (1991) 283 ; A. Djouadi, M. Spira and P. M. Zerwas, Phys. Lett. B [**264**]{} (1991) 440 ; M. Spira, A. Djouadi, D. Graudenz and P. M. Zerwas, Nucl. Phys. B [**453**]{} (1995) 17 ; R. V. Harlander and W. B. Kilgore, Phys. Rev. D [**64**]{} (2001) 013015 ; S. Catani, D. de Florian and M. Grazzini, JHEP [**0105**]{} (2001) 025 ; S. Catani, D. de Florian, M. Grazzini and P. Nason, JHEP [**0307**]{} (2003) 028 ; R. V. Harlander and W. B. Kilgore, Phys. Rev. Lett.  [**88**]{} (2002) 201801 ; C. Anastasiou and K. Melnikov, Nucl. Phys. B [**646**]{} (2002) 220 ; V. Ravindran, J. Smith and W. L. van Neerven, Nucl. Phys. B [**665**]{} (2003) 325 .

P. Bolzoni, F. Maltoni, S. -O. Moch and M. Zaro, Phys. Rev. Lett.  [**105**]{} (2010) 011801 .

T. Han and S. Willenbrock, Phys. Lett. B [**273**]{} (1991) 167; O. Brein, A. Djouadi and R. Harlander, Phys. Lett. B [**579**]{} (2004) 149 .

D. A. Dicus and S. Willenbrock, Phys. Rev. D [**39**]{} (1989) 751 ; D. Dicus, T. Stelzer, Z. Sullivan and S. Willenbrock, Phys. Rev. D [**59**]{} (1999) 094016 ; F. Maltoni, Z. Sullivan and S. Willenbrock, Phys. Rev. D [**67**]{} (2003) 093005 ; F. I. Olness and W. -K. Tung, Nucl. Phys. B [**308**]{} (1988) 813 ; J. F. Gunion, H. E. Haber, F. E. Paige, W. -K. Tung and S. S. D. Willenbrock, Nucl. Phys. B [**294**]{} (1987) 621 ; R. V. Harlander and W. B. Kilgore, Phys. Rev. D [**68**]{} (2003) 013001 .

L. Reina and S. Dawson, Phys. Rev. Lett.  [**87**]{} (2001) 201804 ; W. Beenakker, S. Dittmaier, M. Kramer, B. Plumper, M. Spira and P. M. Zerwas, Phys. Rev. Lett.  [**87**]{} (2001) 201805 ; S. Dawson, L. H. Orr, L. Reina and D. Wackeroth, Phys. Rev. D [**67**]{} (2003) 071503 ; W. Beenakker, S. Dittmaier, M. Kramer, B. Plumper, M. Spira and P. M. Zerwas, Nucl. Phys. B [**653**]{} (2003) 151 ; R. Raitio and W. W. Wada, Phys. Rev. D [**19**]{} (1979) 941 ; Z. Kunszt, Nucl. Phys. B [**247**]{} (1984) 339.

C. Anastasiou, K. Melnikov and F. Petriello, Nucl. Phys. B [**724**]{} (2005) 197 .

C. Anastasiou, G. Dissertori and F. Stöckli, JHEP [**0709**]{} (2007) 018 . M. Grazzini, JHEP [**0802**]{} (2008) 043 . R. Boughezal, F. Caola, K. Melnikov, F. Petriello and M. Schulze, JHEP [**1306**]{} (2013) 072 .

V. Del Duca, A. Frizzo and F. Maltoni, JHEP [**0405**]{} (2004) 064 ; L. J. Dixon, E. W. N. Glover and V. V. Khoze, JHEP [**0412**]{} (2004) 015 ; S. Badger, E. W. Nigel Glover, P. Mastrolia and C. Williams, JHEP [**1001**]{} (2010) 036 .

T. Gehrmann, M. Jaquier, E. W. N. Glover and A. Koukoutsakis, JHEP [**1202**]{} (2012) 056 .

T. Gehrmann and E. Remiddi, Nucl. Phys. B [**601**]{} (2001) 248 ; Nucl. Phys. B [**601**]{} (2001) 287.

T. Binoth and G. Heinrich, Nucl. Phys.  B [**585**]{} (2000) 741 ; Nucl. Phys. B [**693**]{} (2004) 134 ; C. Anastasiou, K. Melnikov and F. Petriello, Phys. Rev. D [**69**]{} (2004) 076010 ; G. Heinrich, Int. J. Mod. Phys.  A [**23**]{} (2008) 1457 ; J. Carter and G. Heinrich, Comput. Phys. Commun.  [**182**]{} (2011) 1566 ; C. Anastasiou, F. Herzog, A. Lazopoulos, JHEP [**1103**]{} (2011) 038 ; JHEP [**1203**]{} (2012) 035 .

S. Catani, M. Grazzini, Phys. Rev. Lett.  [**98** ]{} (2007) 222002 .

A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, JHEP [**0509**]{} (2005) 056 ; A. Gehrmann-De Ridder, T. Gehrmann, E. W. N. Glover, G. Heinrich, JHEP [**0711**]{} (2007) 058 . A. Daleo, T. Gehrmann, D. Maitre, JHEP [**0704** ]{} (2007) 016 ; A. Daleo, A. Gehrmann-De Ridder, T. Gehrmann, G. Luisoni, JHEP [**1001**]{} (2010) 118 ; R. Boughezal, A. Gehrmann-De Ridder, M. Ritzmann, JHEP [**1102**]{} (2011) 098 ; T. Gehrmann, P. F. Monni, \[arXiv:1107.4037\] ; E. W. N.  Glover, J. Pires, JHEP [**1006**]{} (2010) 096 .

S. Frixione, Z. Kunszt and A. Signer, Nucl. Phys. B [**467**]{} (1996) 399 .

P. Baernreuther, M. Czakon and A. Mitov, Phys. Rev. Lett. [**109**]{}, 132001 (2012); M. Czakon and A. Mitov, JHEP [**1212**]{} (2012) 054 ; JHEP [**1301**]{} (2013) 080 .

R. Boughezal, K. Melnikov and F. Petriello, Phys. Rev. D [**85**]{} (2012) 034025 .

E. Remiddi and J.A.M. Vermaseren, Int. J. Mod. Phys. A [**15**]{} (2000) 725 .

T. Kinoshita, J. Math. Phys.  [**3**]{} (1962) 650.

T. D. Lee and M. Nauenberg, Phys. Rev.  [**133**]{} (1964) B1549.

S. Catani, Phys. Lett. B427 (1998) 161 . G. Sterman and M.E. Tejeda-Yeomans, Phys. Lett. B [**552**]{} (2003) 48 .

S. M. Aybat, L. J. Dixon and G. F. Sterman, Phys. Rev. Lett.  [**97**]{} (2006) 072001 . S. M. Aybat, L. J. Dixon and G. F. Sterman, Phys. Rev. D [**74**]{} (2006) 074004 . V. Ravindran, J. Smith and W. L. van Neerven, Nucl. Phys. B [**704**]{} (2005) 332 .

T. Becher and M. Neubert, Phys. Rev. Lett.  [**102**]{} (2009) 162001 .

E. Gardi and L. Magnea, JHEP [**0903**]{} (2009) 079 .

P. Nogueira, J. Comput. Phys.  [**105**]{} (1993) 279.

J. A. M. Vermaseren, math-ph/0010025 ; Nucl. Phys. Proc. Suppl.  [**183**]{} (2008) 19 . R. N. Lee, arXiv:1212.2685 \[hep-ph\] ; arXiv:1310.1145 \[hep-ph\].

F.V. Tkachov, Phys. Lett. [**100B**]{} (1981) 65;\
K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. [**B192**]{} (1981) 159.

T. Gehrmann and E. Remiddi, Nucl. Phys. B [**580**]{} (2000) 485 . T. Ahmed, M. Mahakhud, P. Mathews, N. Rana and V. Ravindran, JHEP [**1405**]{} (2014) 107.

T. Gehrmann and E. Remiddi, Nucl. Phys. B [**640**]{} (2002) 379 .

[^1]: 

[^2]: We thank T. Gehrmann for providing relevant analytically continued HPLs and 2d HPLs.
