text:
# 1 Introduction

## 1 Introduction

The Higgs boson, which is the corner stone of the standard model, is the only particle which has not been discovered yet. Its discovery or its absence will shed light on the mechanism how particles acquire mass as well as answer questions about supersymmetric extensions of the standard model or about compositeness of the existing particles and the Higgs boson. The LEP experiments  give a lower mass limit of about $`m_H114\mathrm{GeV}/\mathrm{c}^2`$ and fits to the data using precision calculations in the electro-weak sector of the standard model indicate an upper limit $`m_H<200\mathrm{GeV}/\mathrm{c}^2`$ with $`95\%`$ confidence level. After the end of the LEP program the search for the Higgs will be continued at hadron colliders in particular at the TEVATRON and the LHC.

In this paper we concentrate on Higgs production where the lowest order reaction proceeds via the gluon-gluon fusion mechanism. In the standard model the Higgs boson couples to the gluons via heavy quark loops. Since the coupling of the scalar Higgs boson $`\mathrm{H}`$ to a fermion loop is proportional to the mass of the fermion (for a review see ), the top-quark loop is the most important. The latter contribution is also dominant for the pseudo-scalar Higgs boson $`\mathrm{A}`$ provided $`\mathrm{tan}\beta `$ is small where $`\beta `$ denotes the mixing angle in the Two-Higgs-Doublet model (2HDM). In lowest order (LO) the gluon-gluon fusion process $`g+g\mathrm{B}`$ with $`\mathrm{B}=\mathrm{H},\mathrm{A}`$, represented by the top-quark triangle graph, was computed in . The two-to-two body tree graphs, given by gluon bremsstrahlung $`g+gg+\mathrm{B}`$, $`g+q(\overline{q})q(\overline{q})+\mathrm{B}`$ and $`q+\overline{q}g+\mathrm{B}`$ were computed for $`\mathrm{B}=\mathrm{H}`$ in  and for $`\mathrm{B}=\mathrm{A}`$ in . From these reactions one can derive the transverse momentum ($`p_T`$) and rapidity ($`y`$) distributions of the (pseudo-) scalar Higgs boson. The total integrated cross section, which also involves the computation of the QCD corrections to the top-quark loop, has been calculated in . This calculation is rather cumbersome since it involves the computation of two-loop triangular graphs with massive quarks. Furthermore also the two-to-three parton reactions have been computed in  using helicity methods. From the experience gained from the next-to-leading (NLO) corrections presented in  it is clear that it will be very difficult to obtain the exact NLO corrections to one-particle inclusive distributions as well as the NNLO corrections to the total cross section.

Fortunately one can simplify the calculations if one takes the large top-quark mass limit $`m_t\mathrm{}`$. In this case the Feynman graphs are obtained from an effective Lagrangian describing the direct coupling of the (pseudo-) scalar Higgs boson to the gluons. The LO and NLO contributions to the total cross section in this approximation were computed in . A thorough analysis ,  reveals that the error introduced by taking the $`m_t\mathrm{}`$ limit is less than about $`5\%`$ provided $`m_H2m_t`$. The two-to-three body processes were computed with the effective Lagrangian approach for the scalar and pseudo-scalar Higgs bosons in  and  respectively using helicity methods. The one-loop corrections to the two-to-two body reactions above were computed for the scalar Higgs boson in  and the pseudo-scalar Higgs boson in . These matrix elements were used to compute the transverse momentum and rapidity distributions of the scalar Higgs boson up to NLO in ,  and the pseudo-scalar Higgs boson in . The effective Lagrangian method was also applied to obtain the NNLO total cross section for scalar Higgs production by the calculation of the two-loop corrections to the Higgs-gluon-gluon vertex in , the soft-plus-virtual gluon corrections in ,  and the computation of the two-to-three body processes in , . These calculations were repeated for pseudo-scalar Higgs production in , .

In this paper we recalculate the NNLO corrections to the total cross section for (pseudo-) scalar Higgs production which have been computed recently in - and we found complete agreement with their analytic results when the number of colours is set $`N=3`$. However our method differs from the one presented in ,  and the approach followed in , . The authors in the latter references compute the total cross section using the Cutkosky  technique where one- and two-loop Feynman integrals are cut in certain ways. These Feynman integrals can be computed using various techniques (for more details see the references in ). Furthermore this method been also applied to compute rapidity distributions in NLO which is presented in . The authors in  choose a more conventional method which was already used in , ,  to compute the coefficient functions for the Drell-Yan process. However instead of an exact computation of the $`22`$ and $`23`$ body phase space integrals they expand them around $`x=1`$. Here $`x=m^2/s`$ where $`m`$ and $`\sqrt{s}`$ denote the (pseudo-) scalar Higgs mass and the partonic centre-of-mass energy respectively. Since the coefficient functions can be expressed into a finite number of Polylogarithms of the types $`\mathrm{Li}_n(x)`$, $`\mathrm{S}_{n,p}(x)`$ (see e.g. ) and logarithms of the types $`\mathrm{ln}^ix\mathrm{ln}^j(1x)`$, which are all multiplied by polynomials in $`x`$, one can expand them in the limit $`x1`$ and match them with the expressions coming from the phase space integrals. In this way the coefficients of the above functions are determined. However as is shown in  this procedure is not necessary since one can obtain the Polylogarithms in a direct way by expanding the hypergeometric functions  in $`\epsilon =n4`$ where $`n`$ is the number of dimensions characteristic of $`n`$-dimensional regularization. The point is that there is no essential difference between Higgs production in the effective Lagrangian approach and the Drell-Yan production mechanisms except that the matrix elements in the former case contain higher powers of the invariants in the numerators. In  a program was made using the algebraic manipulation program FORM  to evaluate the phase space integrals analytically. For this calculation we have extended the program so that it can accommodate integrals having higher powers of invariants in the numerator. Our approach can be also used for differential distributions in particular for jet production (see e.g. ). Here one cannot avoid multi-particle phase space integrals. Of course two-to-four body processes are even more cumbersome to deal with but one can at least try to compute the pole terms $`(1/\epsilon )^k`$, which arise from infrared and collinear divergences, analytically or numerically (for numerical methods see ). Another difference with the previous results is that we present the radiative parts of the coefficient functions for general colour factors of the local gauge group SU(N). The latter is important if one wants to resum the dominant contributions. (see , ).

Our paper will be organized as follows. In Section 2 we give definitions and discuss the regularization method to compute the partonic cross sections in the effective Lagrangian approach. In Section 3 we present the calculation of the coefficient functions. In particular we discuss the various frames in which the three-particle phase space integrals are computed. In Section 4 we present the total cross sections for scalar Higgs production at the LHC and the TEVATRON. The long expressions for the coefficient functions expressed in the various colour factors are presented in Appendix A (scalar Higgs) and Appendix B (pseudo-scalar Higgs).

## 2 Application of the effective Lagrangian approach to Higgs production

In the large top-quark mass limit the Feynman rules (see e.g. ) for scalar Higgs production ($`\mathrm{H}`$) can be derived from the following effective Lagrangian density

$`_{eff}^\mathrm{H}=G_\mathrm{H}\mathrm{\Phi }^\mathrm{H}(x)O(x)\text{with}O(x)={\displaystyle \frac{1}{4}}G_{\mu \nu }^a(x)G^{a,\mu \nu }(x),`$ (2.1)
whereas pseudo-scalar Higgs ($`\mathrm{A}`$) production is obtained from

$`_{eff}^\mathrm{A}=\mathrm{\Phi }^\mathrm{A}(x)\left[G_\mathrm{A}O_1(x)+\stackrel{~}{G}_\mathrm{A}O_2(x)\right]\text{with}`$
$`O_1(x)={\displaystyle \frac{1}{8}}ϵ_{\mu \nu \lambda \sigma }G_a^{\mu \nu }G_a^{\lambda \sigma }(x),`$
$`O_2(x)={\displaystyle \frac{1}{2}}^\mu {\displaystyle \underset{i=1}{\overset{n_f}{}}}\overline{q}_i(x)\gamma _\mu \gamma _5q_i(x),`$ (2.2)
where $`\mathrm{\Phi }^\mathrm{H}(x)`$ and $`\mathrm{\Phi }^\mathrm{A}(x)`$ represent the scalar and pseudo-scalar fields respectively and $`n_f`$ denotes the number of light flavours. Furthermore the gluon field strength is given by $`G_a^{\mu \nu }`$ and the quark field is denoted by $`q_i`$. The factors multiplying the operators are chosen in such a way that the vertices are normalised to the effective coupling constants $`G_\mathrm{H}`$, $`G_\mathrm{A}`$ and $`\stackrel{~}{G}_\mathrm{A}`$. The latter are determined by the top-quark triangular loop graph, including all QCD corrections, taken in the limit $`m_t\mathrm{}`$ which describes the decay process $`\mathrm{B}g+g`$ with $`\mathrm{B}=\mathrm{H},\mathrm{A}`$ namely

$`G_\mathrm{B}`$ $`=`$ $`2^{5/4}a_s(\mu _r^2)G_F^{1/2}\tau _\mathrm{B}F_\mathrm{B}(\tau _\mathrm{B})𝒞_\mathrm{B}(a_s(\mu _r^2),{\displaystyle \frac{\mu _r^2}{m_t^2}}),`$
$`\stackrel{~}{G}_\mathrm{A}`$ $`=`$ $`\left[a_s(\mu _r^2)C_F\left({\displaystyle \frac{3}{2}}3\mathrm{ln}{\displaystyle \frac{\mu _r^2}{m_t^2}}\right)+\mathrm{}\right]G_\mathrm{A},`$ (2.3)
and $`a_s(\mu _r^2)`$ is defined by

$`a_s(\mu _r^2)={\displaystyle \frac{\alpha _s(\mu _r^2)}{4\pi }},`$ (2.4)
where $`\alpha _s(\mu _r^2)`$ is the running coupling constant and $`\mu _r`$ denotes the renormalization scale. Further $`G_F`$ represents the Fermi constant and the functions $`F_\mathrm{B}`$ are given by

$`F_\mathrm{H}(\tau )=1+(1\tau )f(\tau ),F_\mathrm{A}(\tau )=f(\tau )\mathrm{cot}\beta ,`$
$`\tau ={\displaystyle \frac{4m_t^2}{m^2}},`$
$`f(\tau )=\mathrm{arcsin}^2{\displaystyle \frac{1}{\sqrt{\tau }}},\text{for}\tau 1,`$
$`f(\tau )={\displaystyle \frac{1}{4}}\left(\mathrm{ln}{\displaystyle \frac{1\sqrt{1\tau }}{1+\sqrt{1\tau }}}+\pi i\right)^2\text{for}\tau <1,`$ (2.5)
where $`\mathrm{cot}\beta `$ denotes the mixing angle in the Two-Higgs-Doublet Model. Further $`m`$ and $`m_t`$ denote the masses of the (pseudo-) scalar Higgs boson and the top quark respectively. In the large $`m_t`$-limit we have

$`\underset{\tau \mathrm{}}{\text{lim}}F_\mathrm{H}(\tau )={\displaystyle \frac{2}{3\tau }},\underset{\tau \mathrm{}}{\text{lim}}F_\mathrm{A}(\tau )={\displaystyle \frac{1}{\tau }}\mathrm{cot}\beta .`$ (2.6)
The coefficient functions $`𝒞_\mathrm{B}`$ originate from the corrections to the top-quark triangular graph provided one takes the limit $`m_t\mathrm{}`$. We have presented the Born level couplings $`G_\mathrm{B}`$ in Eq.(2) for general $`m_t`$ for on-shell gluons only in order to keep some part of the top-quark mass dependence. This is an approximation because the gluons which couple to the H and A bosons via the top-quark loop in partonic subprocesses are very often virtual. The virtual-gluon momentum dependence is neither described by $`F_\mathrm{B}(\tau )`$ nor by $`𝒞_\mathrm{B}`$. However for total cross sections the main contribution comes from the region where the gluons are almost on-shell so that this approximation is better than it is for differential cross sections with large transverse momentum. The coefficient functions are computed up to order $`\alpha _s^2`$ in ,  for the H and in  for the A. They read as follows

$`𝒞_\mathrm{H}(a_s(\mu _r^2),{\displaystyle \frac{\mu _r^2}{m_t^2}})=1+a_s^{(5)}(\mu _r^2)[5C_A3C_F]+\left(a_s^{(5)}(\mu _r^2)\right)^2[{\displaystyle \frac{27}{2}}C_F^2`$
$`{\displaystyle \frac{100}{3}}C_AC_F+{\displaystyle \frac{1063}{36}}C_A^2{\displaystyle \frac{4}{3}}C_FT_f{\displaystyle \frac{5}{6}}C_AT_f+(7C_A^2`$
$`11C_AC_F)\mathrm{ln}{\displaystyle \frac{\mu _r^2}{m_t^2}}+n_fT_f(5C_F{\displaystyle \frac{47}{9}}C_A+8C_F\mathrm{ln}{\displaystyle \frac{\mu _r^2}{m_t^2}})],`$ (2.7)
$`𝒞_\mathrm{A}(a_s(\mu _r^2),{\displaystyle \frac{\mu _r^2}{m_t^2}})=1,`$ (2.8)
where $`a_s^{(5)}`$ is presented in the five-flavour number scheme. Notice that the coefficient function in Eq. (2.7) is derived for general colour factors of the group SU(N) from Eq. (6) in . These factors are given by

$`C_A=N,C_F={\displaystyle \frac{N^21}{2N}},T_f={\displaystyle \frac{1}{2}}.`$ (2.9)
Notice that $`T_f`$ is also incorporated into $`G_\mathrm{B}`$ in Eq. (2) where it is set to the value $`T_f=1/2`$.

Using the effective Lagrangian approach we will calculate the total cross section of the reaction

$`H_1(P_1)+H_2(P_2)\mathrm{B}(p_5)+^{}X^{},`$ (2.10)
where $`H_1`$ and $`H_2`$ denote the incoming hadrons and $`X`$ represents an inclusive hadronic state. The total cross section is given by

$`\sigma _{\mathrm{tot}}={\displaystyle \frac{\pi G_\mathrm{B}^2}{8(N^21)}}{\displaystyle \underset{a,b=q,\overline{q},g}{}}{\displaystyle _x^1}𝑑x_1{\displaystyle _{x/x_1}^1}𝑑x_2f_a(x_1,\mu ^2)f_b(x_2,\mu ^2)`$
$`\times \mathrm{\Delta }_{ab,\mathrm{B}}({\displaystyle \frac{x}{x_1x_2}},{\displaystyle \frac{m^2}{\mu ^2}}),\mathrm{B}=\mathrm{H},\mathrm{A},`$
$`\text{with}x={\displaystyle \frac{m^2}{S}},S=(P_1+P_2)^2,p_5^2=m^2,`$ (2.11)
where the factor $`1/(N^21)`$ originates from the colour average in the case of the local gauge group SU(N). Further we have assumed that the (pseudo-) scalar Higgs boson is mainly produced on-shell i.e. $`p_5^2m^2`$. The parton densities denoted by $`f_a(y,\mu ^2)`$ ($`a,b=q,\overline{q},g`$) depend on the mass factorization/renormalization scale $`\mu `$. The same scales also enter the coefficient functions $`\mathrm{\Delta }_{ab,\mathrm{B}}`$ which are derived from the partonic cross sections

$`\sigma _{ab,\mathrm{B}}(z,{\displaystyle \frac{m^2}{\mu ^2}})={\displaystyle \frac{\pi }{s}}{\displaystyle \frac{d^np_5}{(2\pi )^n}\delta ^+(p_5^2m^2)T_{ab,\mathrm{B}}(p_5,p_1,p_2)},`$ (2.12)
where the incoming parton momenta, denoted by $`p_1`$ and $`p_2`$, are related to the hadron momenta by

$`p_1=x_1P_1,p_2=x_2P_2,`$
$`s=(p_1+p_2)^2,s=x_1x_2S,z={\displaystyle \frac{m^2}{s}}.`$ (2.13)
The amplitude $`T_{ab,\mathrm{B}}`$ can be written as

$`T_{ab,\mathrm{H}}(p_5,p_1,p_2)`$ $`=`$ $`G_\mathrm{H}^2{\displaystyle d^4ye^{ip_5y}a,b|O(y)O(0)|a,b},`$ (2.14)
$`T_{ab,\mathrm{A}}(p_5,p_1,p_2)`$ $`=`$ $`{\displaystyle d^4ye^{ip_5y}a,b|\left(G_\mathrm{A}O_1(y)+\stackrel{~}{G}_\mathrm{A}O_2(y)\right)}`$ (2.15)
$`\times \left(G_\mathrm{A}O_1(0)+\stackrel{~}{G}_\mathrm{A}O_2(0)\right)|a,b.`$
The expressions above for the amplitude $`T_{ab}`$ are similar to those given for the Drell-Yan process except that the conserved electroweak currents are replaced by the operators $`O`$ and $`O_1,O_2`$. The latter are not conserved so that they acquire additional renormalization constants defined by

$`O(y)=Z_O\widehat{O}(y),O_i(y)=Z_{ij}\widehat{O}_j(y).`$ (2.16)
where the hat indicates that the quantities under consideration are unrenormalized. Insertion of the above equations into Eqs. (2.14),(2.15) leads to the renormalized expressions

$`T_{ab,\mathrm{H}}(p_5,p_1,p_2)=G_\mathrm{H}^2Z_O^2{\displaystyle d^4ye^{ip_5y}a,b|\widehat{O}(y)\widehat{O}(0)|a,b},`$ (2.17)
$`T_{ab,\mathrm{A}}(p_5,p_1,p_2)={\displaystyle }d^4ye^{ip_5y}[\{G_\mathrm{A}^2Z_{11}^2+\stackrel{~}{G}_\mathrm{A}^2Z_{21}^2+2G_\mathrm{A}\stackrel{~}{G}_\mathrm{A}Z_{11}Z_{21}\}`$
$`\times a,b|\widehat{O}_1(y)\widehat{O}_1(0)|a,b+\left\{G_\mathrm{A}^2Z_{12}^2+\stackrel{~}{G}_\mathrm{A}^2Z_{22}^2+2G_\mathrm{A}\stackrel{~}{G}_\mathrm{A}Z_{12}Z_{22}\right\}`$
$`\times a,b|\widehat{O}_2(y)\widehat{O}_2(0)|a,b+\{G_\mathrm{A}^2Z_{11}Z_{12}+\stackrel{~}{G}_\mathrm{A}^2Z_{22}Z_{21}+G_\mathrm{A}\stackrel{~}{G}_\mathrm{A}`$
$`\times (Z_{11}Z_{22}+Z_{12}Z_{21})\}a,b|\widehat{O}_1(y)\widehat{O}_2(0)+\widehat{O}_2(y)\widehat{O}_1(0)|a,b].`$ (2.18)
The operator renormalization constants depend on the regularization scheme in particular on the prescription for the $`\gamma _5`$-matrix and the Levi-Civita tensor in Eq. (2). The computation of $`T_{ab,\mathrm{H}}`$ will be carried out by choosing n-dimensional regularization where in the case of $`T_{ab,\mathrm{A}}`$ we adopt the HVBM prescription ,  for the $`\gamma _5`$-matrix. For this choice the contraction of the Levi-Civita tensors proceeds as

$`ϵ_{\mu _1\nu _1\lambda _1\sigma _1}ϵ^{\mu _2\nu _2\lambda _2\sigma _2}=\left|\begin{array}{cccc}\delta _{\mu _1}^{\mu _2}& \delta _{\mu _1}^{\nu _2}& \delta _{\mu _1}^{\lambda _2}& \delta _{\mu _1}^{\sigma _2}\\ \delta _{\nu _1}^{\mu _2}& \delta _{\nu _1}^{\nu _2}& \delta _{\nu _1}^{\lambda _2}& \delta _{\nu _1}^{\sigma _2}\\ \delta _{\lambda _1}^{\mu _2}& \delta _{\lambda _1}^{\nu _2}& \delta _{\lambda _1}^{\lambda _2}& \delta _{\lambda _1}^{\sigma _2}\\ \delta _{\sigma _1}^{\mu _2}& \delta _{\sigma _1}^{\nu _2}& \delta _{\sigma _1}^{\lambda _2}& \delta _{\sigma _1}^{\sigma _2}\end{array}\right|,`$ (2.23)
where all Lorentz indices are taken to be n-dimensional. To facilitate the calculation one can replace $`\gamma _\mu \gamma _5`$ in Eq. (2) by (see )

$`\gamma _\mu \gamma _5={\displaystyle \frac{i}{6}}ϵ_{\mu \rho \sigma \tau }\gamma ^\rho \gamma ^\sigma \gamma ^\tau ,`$ (2.24)
which is equivalent to the HVBM scheme. Choosing the $`\overline{\mathrm{MS}}`$ subtraction scheme the renormalization constant corresponding to the operator $`O`$ becomes

$`Z_O`$ $`=`$ $`1+a_s(\mu _r^2)S_\epsilon {\displaystyle \frac{2}{\epsilon }}\beta _0+a_s^2(\mu _r^2)S_\epsilon ^2\left[{\displaystyle \frac{4}{\epsilon ^2}}\beta _0^2+{\displaystyle \frac{2}{\epsilon }}\beta _1\right]+\mathrm{}`$ (2.25)
where $`S_\epsilon `$ denotes the spherical factor characteristic of n-dimensional regularization. It is defined by

$`S_\epsilon =\mathrm{exp}\left\{{\displaystyle \frac{\epsilon }{2}}\left[\gamma _E\mathrm{ln}4\pi \right]\right\}.`$ (2.26)
The lowest order coefficients $`\beta _0`$ and $`\beta _1`$ originate from the beta-function given by

$`\beta (\alpha _s)=a_s^2(\mu _r^2)\beta _0a_s^3(\mu _r^2)\beta _1+\mathrm{},`$
$`\beta _0={\displaystyle \frac{11}{3}}C_A{\displaystyle \frac{4}{3}}n_fT_f,\beta _1={\displaystyle \frac{34}{3}}C_A^24n_fT_fC_F{\displaystyle \frac{20}{3}}n_fT_fC_A.`$
The operator renormalization constants corresponding to $`O_1`$ and $`O_2`$ are computed in  and they read

$`Z_{11}=Z_{\alpha _s}=1+a_s(\mu _r^2)S_\epsilon {\displaystyle \frac{2}{\epsilon }}\beta _0+a_s^2(\mu _r^2)S_\epsilon ^2\left[{\displaystyle \frac{4}{\epsilon ^2}}\beta _0^2+{\displaystyle \frac{1}{\epsilon }}\beta _1\right]+\mathrm{}`$ (2.28)
where $`Z_{\alpha _s}`$ denotes the coupling constant renormalization factor defined by

$`\widehat{a}_s=Z_{\alpha _s}a_s(\mu _r^2).`$ (2.29)
The remaining constants are

$`Z_{21}`$ $`=`$ $`0,`$ (2.30)
$`Z_{12}`$ $`=`$ $`a_s(\mu _r^2)S_\epsilon {\displaystyle \frac{1}{\epsilon }}\left[6C_F\right],`$ (2.31)
$`Z_{22}`$ $`=`$ $`Z_{\mathrm{MS}}^sZ_5^s,`$ (2.32)
where $`Z_{\mathrm{MS}}^s`$ and $`Z_5^s`$ are the constants characteristic of the HVBM scheme. They are given by

$`Z_{\mathrm{MS}}^s`$ $`=`$ $`1+a_s^2(\mu _r^2)S_\epsilon {\displaystyle \frac{1}{\epsilon }}\left[{\displaystyle \frac{44}{3}}C_AC_Fn_fT_fC_F{\displaystyle \frac{20}{3}}\right],`$
$`Z_5^s`$ $`=`$ $`1+a_s(\mu _r^2)\left[4C_F\right]+a_s^2(\mu _r^2)\left[22C_F^2{\displaystyle \frac{107}{9}}C_AC_F+{\displaystyle \frac{31}{9}}n_fT_fC_F\right].`$
The latter renormalization constant is determined in such a way that the Adler-Bell-Jackiw anomaly

$`O_2(y)=4a_s(\mu _r^2)n_fT_fO_1(y),`$ (2.34)
is preserved in all orders in perturbation theory according to the Adler-Bardeen theorem .

## 3 Computation of the partonic cross section up to NNLO

In this section we will give a short outline of the computation of the partonic cross sections $`\sigma _{ab,\mathrm{B}}`$ in Eq. (2.12) up to next-to-next-leading order (NNLO). Insertion of a complete set of intermediate states and using translational invariance in Eq. (2) we obtain for scalar Higgs production

$`\sigma _{ab,\mathrm{H}}=K_{ab}{\displaystyle \frac{\pi G_\mathrm{H}^2}{s}}Z_O^2{\displaystyle d^np_5\delta ^+(p_5^2m^2)\underset{m=3}{\overset{\mathrm{}}{}}\underset{i=3,i5}{\overset{m}{}}\frac{d^np_i}{(2\pi )^{n1}}\delta ^+(p_i^2)}`$
$`\times \delta ^{(n)}({\displaystyle \underset{j=1}{\overset{m}{}}}p_j)|M_{abX\mathrm{H}}|^2,`$
$`M_{abX\mathrm{H}}=p_1,p_2|\widehat{O}(0)|X,p_5\text{with}|X,p_5=|p_3,p_4,p_6\mathrm{}p_m,p_5,`$
where $`K_{ab}`$ represents the spin and colour average over the initial states. Further one can write down a similar expression for the pseudo-scalar Higgs (see Eq. (2.15)). Up to NNLO we have to compute the following subprocesses. On the Born level we have the reaction

$`g+g\mathrm{B}.`$ (3.2)
In NLO we have in addition to the one-loop virtual corrections to the above reaction the following two-to-two body processes

$`g+g\mathrm{B}+g,g+q(\overline{q})\mathrm{B}+q(\overline{q}),q+\overline{q}\mathrm{B}+g.`$ (3.3)
In NNLO we receive contributions from the two-loop virtual corrections to the Born process in Eq. (3.2) and the one-loop corrections to the reactions in Eq. (3.3). To these contribution one has to add the results obtained from the following two-to-three body reactions

$`g+g\mathrm{B}+g+g,g+g\mathrm{B}+q_i+\overline{q}_i,`$ (3.4)
$`g+q(\overline{q})\mathrm{B}+q(\overline{q})+g,`$ (3.5)
$`q+\overline{q}\mathrm{B}+g+g,q+\overline{q}\mathrm{B}+q_i+\overline{q}_i,`$ (3.6)
$`q_1+q_2\mathrm{B}+q_1+q_2,q_1+\overline{q}_2\mathrm{B}+q_1+\overline{q}_2,`$ (3.7)
$`q+q\mathrm{B}+q+q.`$ (3.8)
In the case of pseudo-scalar Higgs production one also has to add the contributions due to interference terms coming from the operators $`\widehat{O}_1`$ and $`\widehat{O}_2`$ in Eq. (2.18). Up to order $`\alpha _s^2`$ we have to compute the following expression for the reactions in (3.4)

$`\sigma _{gg,\mathrm{A}}=K_{ab}{\displaystyle \frac{\pi }{s}}{\displaystyle }d^np_5\delta ^+(p_5^2m^2)[\{G_\mathrm{A}^2Z_{11}^2{\displaystyle \underset{m=3}{\overset{4}{}}}{\displaystyle \underset{i=3}{\overset{m}{}}}{\displaystyle }{\displaystyle \frac{d^np_i}{(2\pi )^{n1}}}\delta ^+(p_i^2)`$ (3.9)
$`\times \delta ^{(n)}({\displaystyle \underset{j=1}{\overset{m}{}}}p_j)|p_1,p_2|\widehat{O}_1|X,p_5|^2\}+\{G_\mathrm{A}\stackrel{~}{G}_\mathrm{A}\delta ^{(n)}(p_1+p_2+p_5)`$
$`\times (p_1,p_2|\widehat{O}_1|p_5p_5|\widehat{O}_2|p_1,p_2+p_1,p_2|\widehat{O}_2|p_5p_5|\widehat{O}_1|p_1,p_2)\}].`$
The Feynman graphs for $`p_1,p_2|\widehat{O}_1|p_5`$ and $`p_1,p_2|\widehat{O}_2|p_5`$ can be found in Fig. 1a and Fig. 2a of  respectively. From the above equation one infers that the interference term contributes as a delta-function $`\delta (1z)`$ to $`\sigma _{gg,\mathrm{A}}`$. Furthermore since $`\stackrel{~}{G}_\mathrm{A}`$ is proportional to $`G_\mathrm{A}`$ (see Eq. (2)) one can extract the latter constant as an overall factor from the above equation. For the reactions in Eq. (3.5) we have to compute

$`\sigma _{gq,\mathrm{A}}=K_{ab}{\displaystyle \frac{\pi G_\mathrm{A}^2}{s}}{\displaystyle }d^np_5\delta ^+(p_5^2m^2)[\{Z_{11}^2{\displaystyle \underset{m=3}{\overset{4}{}}}{\displaystyle \underset{i=3}{\overset{m}{}}}{\displaystyle }{\displaystyle \frac{d^np_i}{(2\pi )^{n1}}}\delta ^+(p_i^2)`$ (3.10)
$`\times \delta ^{(n)}({\displaystyle \underset{j=1}{\overset{m}{}}}p_i)|p_1,p_2|\widehat{O}_1|X,p_5|^2\}+\{Z_{12}{\displaystyle }{\displaystyle \frac{d^np_3}{(2\pi )^{n1}}}\delta ^+(p_3^2)`$
$`\times \delta ^{(n)}(p_1+p_2+p_3+p_5)(p_1,p_2|\widehat{O}_1|p_5,p_3p_5,p_3|\widehat{O}_2|p_1,p_2`$
$`+p_1,p_2|\widehat{O}_2|p_5,p_3p_5,p_3|\widehat{O}_1|p_1,p_2)\}].`$
The Feynman graphs for $`p_1,p_2|\widehat{O}_1|p_3,p_5`$ and $`p_1,p_2|\widehat{O}_2|p_3,p_5`$ can be found in Fig. 1b and Fig. 2b of  respectively. In the HVBM-scheme the latter matrix element is proportional to $`\epsilon =n4`$. However this contribution does not vanish in the limit $`\epsilon 0`$ because of the single pole term present in $`Z_{12}`$ in Eq. (2.31). The same expression as in Eq. (3.10) also exists for $`\sigma _{q\overline{q},\mathrm{A}}`$ originating from the first reaction in Eq. (3.6). In this case the Feynman graphs corresponding to $`p_1,p_2|\widehat{O}_1|p_3,p_5`$ and $`p_1,p_2|\widehat{O}_2|p_3,p_5`$ are shown in Fig. 1d and Fig. 2d of  respectively. Notice that these contributions also survive in the calculation of differential distributions, see .

The computation of the phase space integrals proceeds as follows. First we take over the expressions from the one- and two-loop corrections to the Born reaction in Eq.(3.2) which are computed in  (one-loop) and  (two-loop). Unfortunately the two-loop result is not presented for arbitrary colour factors so that we cannot distinguish between $`n_fC_AT_f`$ and $`n_fC_FT_f`$. Hence we have shuffled all contributions proportional to $`n_f`$ into the term proportional to $`n_fC_AT_f`$. To compute the two-to-two body processes including the virtual corrections we have chosen the centre-of-mass frame of the incoming partons given by

$`p_1`$ $`=`$ $`{\displaystyle \frac{1}{2}}\sqrt{P_{12}}(1,0,\mathrm{},0,1),`$
$`p_2`$ $`=`$ $`{\displaystyle \frac{1}{2}}\sqrt{P_{12}}(1,0,\mathrm{},0,1),`$
$`p_3`$ $`=`$ $`{\displaystyle \frac{P_{12}m^2}{2\sqrt{P_{12}}}}(1,0,\mathrm{},\mathrm{sin}\theta ,\mathrm{cos}\theta ),`$
$`p_5`$ $`=`$ $`{\displaystyle \frac{1}{2\sqrt{P_{12}}}}(P_{12}+m^2,0,\mathrm{},(P_{12}m^2)\mathrm{sin}\theta ,(P_{12}m^2)\mathrm{cos}\theta ),`$
where we have introduced the invariants

$`P_{ij}=(p_i+p_j)^2,P_{12}=s.`$ (3.12)
In this frame the two-to-two body phase space integral becomes

$`\sigma _{abc\mathrm{H}}`$ $`=`$ $`K_{ab}Z_O^2{\displaystyle \frac{\pi G_\mathrm{H}^2}{P_{12}}}{\displaystyle \frac{2^{3n}}{(4\pi )^{n/2}}}{\displaystyle \frac{1}{\mathrm{\Gamma }(n/21)}}\left({\displaystyle \frac{P_{12}}{\mu ^2}}\right)^{1n/2}\left({\displaystyle \frac{P_{12}m^2}{\mu ^2}}\right)^{n3}`$ (3.13)
$`\times {\displaystyle _0^\pi }d\theta \mathrm{sin}^{n3}\theta |M_{abc\mathrm{B}}|^2.`$
Here $`\mu `$ indicates the dimension of the strong coupling constant $`gg\mu ^{(4n)/2}`$ which is characteristic of n-dimensional regularization. In order to make all ratios dimensionless we have already included the factor $`\mu ^{(4n)/2}`$ in the phase space integrals. Notice that in principle the scale $`\mu `$ has nothing to do with the factorization or renormalization scale. However for convenience one puts the latter scales equal to $`\mu `$. The evaluation of the integral in Eq. (3.13) is rather easy even when $`|M_{abc\mathrm{B}}|`$ contains virtual contributions. Since the integrals can be evaluated analytically we have made a routine using the algebraic manipulation program FORM which provides us with the results. A similar routine is made for the two-to-three body phase space integrals but here the computation is not so easy unless one chooses a suitable frame. Since we integrate over the total phase space the integrals are Lorentz invariant and therefore frame independent. The matrix elements of the partonic reactions can be partial fractioned in such a way that maximally two factors $`P_{ij}`$ in Eq. (3.12) depend on the polar angle $`\theta _1`$ and the azimuthal angle $`\theta _2`$. Furthermore one factor only depends on $`\theta _1`$ whereas the other one depends both on $`\theta _1`$ and $`\theta _2`$. Therefore the following combinations show up in the matrix elements

$`P_{ij}^kP_{mn}^l,P_{ij}^kP_{m5}^l,P_{i5}^kP_{m5}^l,4k2,4l2,`$
$`p_i^2=p_j^2=p_m^2=p_n^2=0,p_5^2=m^2.`$ (3.14)
For the first combination it is easy to perform the angular integration since all momenta represent massless particles and one obtains a hypergeometric function (see e.g. Eq. (4.19) in ). The angular integral of the second combination is more difficult to compute because one particle is massive and the result is an one dimensional integral over a hypergeometric function which however can be expanded around $`\epsilon `$. Examples of these types of integrals can be found in Appendix C of . The last combination is very difficult to compute in n dimensions because both factors contain the massive particle indicated by $`p_5`$. Therefore one has to avoid this combination at any cost. This is possible if one chooses the following three frames. The first one is the centre-of-mass frame of the incoming partons. The momenta are

$`p_1`$ $`=`$ $`{\displaystyle \frac{1}{2}}\sqrt{P_{12}}(1,0,\mathrm{},0,0,1),`$
$`p_2`$ $`=`$ $`{\displaystyle \frac{1}{2}}\sqrt{P_{12}}(1,0,\mathrm{},0,0,1),`$
$`p_3`$ $`=`$ $`{\displaystyle \frac{P_{12}P_{45}}{2\sqrt{P_{12}}}}(1,0,\mathrm{},0,\mathrm{sin}\theta _1,\mathrm{cos}\theta _1),`$
$`p_4`$ $`=`$ $`{\displaystyle \frac{P_{12}P_{35}}{2\sqrt{P_{12}}}}(1,0,\mathrm{},\mathrm{sin}\psi \mathrm{sin}\theta _2,\mathrm{cos}\psi \mathrm{sin}\theta _1+\mathrm{sin}\psi \mathrm{cos}\theta _2\mathrm{cos}\theta _1,`$
$`\mathrm{cos}\psi \mathrm{cos}\theta _1\mathrm{sin}\psi \mathrm{cos}\theta _2\mathrm{sin}\theta _1),`$
$`\mathrm{sin}^2{\displaystyle \frac{\psi }{2}}`$ $`=`$ $`{\displaystyle \frac{P_{12}P_{34}}{(P_{12}P_{35})(P_{12}P_{45})}},`$ (3.15)
$`\sigma _{abcd\mathrm{H}}`$ $`=`$ $`K_{ab}Z_O^2{\displaystyle \frac{G_\mathrm{H}^2}{2P_{12}}}{\displaystyle \frac{1}{(4\pi )^n}}{\displaystyle \frac{1}{\mathrm{\Gamma }(n3)}}\left({\displaystyle \frac{P_{12}}{\mu ^2}}\right)^{1n/2}{\displaystyle _{m^2}^{P_{12}}}𝑑P_{35}`$ (3.16)
$`\times {\displaystyle _{P_{12}m^2/P_{35}}^{P_{12}+m^2P_{35}}}dP_{45}\left({\displaystyle \frac{P_{35}P_{45}P_{12}m^2}{\mu ^4}}\right)^{n/22}`$
$`\times \left({\displaystyle \frac{P_{12}+m^2P_{35}P_{45}}{\mu ^2}}\right)^{n/22}{\displaystyle _0^\pi }𝑑\theta _1\mathrm{sin}^{n3}\theta _1`$
$`\times {\displaystyle _0^\pi }d\theta _2\mathrm{sin}^{n4}\theta _2|M_{abcd\mathrm{B}}|^2.`$
For the next frame we choose the centre-of-mass frame of the two outgoing particles indicated by the momenta $`p_3`$ and $`p_4`$.

$`p_1`$ $`=`$ $`\omega _1(1,0,\mathrm{},0,0,1),`$
$`p_2`$ $`=`$ $`(\omega _2,0,\mathrm{},0,|\stackrel{}{p}_5|\mathrm{sin}\psi ,|\stackrel{}{p}_5|\mathrm{cos}\psi \omega _1),`$
$`p_3`$ $`=`$ $`{\displaystyle \frac{1}{2}}\sqrt{P_{34}}(1,0,\mathrm{},\mathrm{sin}\theta _1\mathrm{sin}\theta _2,\mathrm{sin}\theta _1\mathrm{cos}\theta _2,\mathrm{cos}\theta _1),`$
$`p_4`$ $`=`$ $`{\displaystyle \frac{1}{2}}\sqrt{P_{34}}(1,0,\mathrm{},\mathrm{sin}\theta _1\mathrm{sin}\theta _2,\mathrm{sin}\theta _1\mathrm{cos}\theta _2,\mathrm{cos}\theta _1),`$
$`p_5`$ $`=`$ $`(\omega _5,0,\mathrm{},0,|\stackrel{}{p}_5|\mathrm{sin}\psi ,|\stackrel{}{p}_5|\mathrm{cos}\psi ),`$
$`\omega _1`$ $`=`$ $`{\displaystyle \frac{P_{12}+P_{15}m^2}{2\sqrt{P_{34}}}},\omega _2={\displaystyle \frac{P_{12}+P_{25}m^2}{2\sqrt{P_{34}}}}`$
$`\omega _5={\displaystyle \frac{P_{15}+P_{25}}{2\sqrt{P_{34}}}},`$
$`\mathrm{cos}\psi `$ $`=`$ $`{\displaystyle \frac{(P_{34}m^2)(P_{15}m^2)P_{12}(P_{25}+m^2)}{(P_{12}+P_{15}m^2)\sqrt{(P_{15}+P_{25})^24m^2P_{34}}}},`$ (3.17)
$`\sigma _{abcd\mathrm{H}}=K_{ab}Z_O^2{\displaystyle \frac{G_\mathrm{H}^2}{2P_{12}}}{\displaystyle \frac{1}{(4\pi )^n}}{\displaystyle \frac{1}{\mathrm{\Gamma }(n3)}}\left({\displaystyle \frac{P_{12}}{\mu ^2}}\right)^{1n/2}{\displaystyle _{m^2P_{12}}^0}𝑑P_{25}`$
$`\times {\displaystyle _{m^2P_{12}P_{25}}^{m^2+P_{12}m^2/(P_{25}m^2)}}dP_{15}\left({\displaystyle \frac{(P_{15}m^2)(P_{25}m^2)P_{12}m^2}{\mu ^4}}\right)^{n/22}`$
$`\times \left({\displaystyle \frac{P_{12}+P_{15}+P_{25}m^2}{\mu ^2}}\right)^{n/22}{\displaystyle _0^\pi }𝑑\theta _1\mathrm{sin}^{n3}\theta _1`$
$`\times {\displaystyle _0^\pi }d\theta _2\mathrm{sin}^{n4}\theta _2|M_{abcd\mathrm{B}}|^2.`$ (3.18)
For the last frame we choose the centre-of-mass frame of one of the outgoing partons and the (pseudo-)scalar Higgs boson indicated by the momenta $`p_4`$ and $`p_5`$ respectively

$`p_1`$ $`=`$ $`\omega _1(1,0,\mathrm{},0,0,1),`$
$`p_2`$ $`=`$ $`(\omega _2,0,\mathrm{},0,\omega _3\mathrm{sin}\psi ,\omega _3\mathrm{cos}\psi \omega _1),`$
$`p_3`$ $`=`$ $`\omega _3(1,0,\mathrm{},0,\mathrm{sin}\psi ,\mathrm{cos}\psi ),`$
$`p_4`$ $`=`$ $`\omega _4(1,0,\mathrm{},\mathrm{sin}\theta _1\mathrm{sin}\theta _2,\mathrm{sin}\theta _1\mathrm{cos}\theta _2,\mathrm{cos}\theta _1)`$
$`p_5`$ $`=`$ $`(\omega _5,0,\mathrm{},\omega _4\mathrm{sin}\theta _1\mathrm{sin}\theta _2,\omega _4\mathrm{sin}\theta _1\mathrm{cos}\theta _2,\omega _4\mathrm{cos}\theta _1),`$
$`\omega _1`$ $`=`$ $`{\displaystyle \frac{P_{12}+P_{13}}{2\sqrt{P_{45}}}},\omega _2={\displaystyle \frac{P_{12}+P_{23}}{2\sqrt{P_{45}}}}`$
$`\omega _3`$ $`=`$ $`{\displaystyle \frac{P_{13}+P_{23}}{2\sqrt{P_{45}}}},\omega _4={\displaystyle \frac{P_{45}m^2}{2\sqrt{P_{45}}}}`$
$`\omega _5`$ $`=`$ $`{\displaystyle \frac{P_{45}+m^2}{2\sqrt{P_{45}}}},\mathrm{cos}\psi ={\displaystyle \frac{P_{12}P_{23}P_{45}P_{13}}{(P_{12}+P_{13})(P_{13}+P_{23})}},\mathrm{`}`$ (3.19)
$`\sigma _{abcd\mathrm{H}}=K_{ab}Z_O^2{\displaystyle \frac{G_\mathrm{H}^2}{2P_{12}}}{\displaystyle \frac{1}{(4\pi )^n}}{\displaystyle \frac{1}{\mathrm{\Gamma }(n3)}}\left({\displaystyle \frac{P_{12}}{\mu ^2}}\right)^{1n/2}{\displaystyle _{m^2P_{12}}^0}𝑑P_{23}`$
$`\times {\displaystyle _{m^2P_{12}P_{23}}^0}dP_{13}\left({\displaystyle \frac{P_{13}P_{23}}{\mu ^4}}\right)^{n/22}\left({\displaystyle \frac{P_{12}+P_{13}+P_{23}m^2}{\mu ^2}}\right)^{n3}`$
$`\times \left({\displaystyle \frac{P_{12}+P_{13}+P_{23}}{\mu ^2}}\right)^{1n/2}{\displaystyle _0^\pi }𝑑\theta _1\mathrm{sin}^{n3}\theta _1{\displaystyle _0^\pi }𝑑\theta _2\mathrm{sin}^{n4}\theta _2`$
$`\times |M_{abcd\mathrm{B}}|^2.`$ (3.20)
The integration over the angular independent invariants $`P_{ij}`$ can be performed by rescaling them with respect to $`P_{12}`$ (see Appendix E in ). The integrals are such that they can be performed in an algebraic way by a program based on FORM  which has been also used to compute the NNLO coefficient functions of the Drell-Yan process in , ,  (for more details see ). The partonic cross sections Eq. (2.12) can be written as follows

$`\sigma _{ab,\mathrm{H}}(z,{\displaystyle \frac{m^2}{\mu ^2}})`$ $`=`$ $`{\displaystyle \frac{\pi G_\mathrm{H}^2}{8(N^21)}}{\displaystyle \frac{1}{1+\epsilon /2}}Z_O^2\widehat{W}_{ab,\mathrm{H}}(z,{\displaystyle \frac{m^2}{\mu ^2}})`$
$`\sigma _{ab,\mathrm{A}}(z,{\displaystyle \frac{m^2}{\mu ^2}})`$ $`=`$ $`{\displaystyle \frac{\pi G_\mathrm{A}^2}{8(N^21)}}{\displaystyle \frac{1+\epsilon }{1+\epsilon /2}}Z_{11}^2\widehat{W}_{ab,\mathrm{A}}(z,{\displaystyle \frac{m^2}{\mu ^2}})`$ (3.21)
$`z^1\widehat{W}_{gg,\mathrm{B}}`$ $`=`$ $`\delta (1z)+\widehat{a}_sS_\epsilon \left({\displaystyle \frac{m^2}{\mu ^2}}\right)^{\epsilon /2}\left[{\displaystyle \frac{2}{\epsilon }}(P_{gg}^{(0)}2\beta _0)+w_{gg,\mathrm{B}}^{(1)}+\epsilon \overline{w}_{gg,\mathrm{B}}^{(1)}\right]`$ (3.22)
$`+\widehat{a}_s^2S_\epsilon ^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)^\epsilon [{\displaystyle \frac{1}{\epsilon ^2}}(2P_{gg}^{(0)}P_{gg}^{(0)}+P_{gq}^{(0)}P_{qg}^{(0)}10\beta _0P_{gg}^{(0)}`$
$`+12\beta _0^2)+{\displaystyle \frac{1}{\epsilon }}(P_{gg}^{(1)}6\beta _0w_{gg,\mathrm{B}}^{(1)}+2P_{gg}^{(0)}w_{gg,\mathrm{B}}^{(1)}`$
$`+2w_{qg,\mathrm{B}}^{(1)}P_{qg}^{(0)}+z_\mathrm{B})6\beta _0\overline{w}_{gg,\mathrm{B}}^{(1)}+2P_{gg}^{(0)}\overline{w}_{gg,\mathrm{B}}^{(1)}`$
$`+2\overline{w}_{qg,\mathrm{B}}^{(1)}P_{qg}^{(0)}+w_{gg,\mathrm{B}}^{(2)}]`$
$`\text{with}z_\mathrm{H}=4\beta _1z_\mathrm{A}=2\beta _1,`$
The coefficients $`z_\mathrm{B}`$, multiplying the single pole terms in the equations above, originate from the operators $`O`$ in Eq. (2.1) and $`O_1`$ in Eq. (2) and they are removed via the operator renormalization constants $`Z_O`$ in Eq. (2.25) and $`Z_{11}`$ in Eq. (2.28).

$`z^1\widehat{W}_{gq,\mathrm{B}}`$ $`=`$ $`\widehat{a}_sS_\epsilon \left({\displaystyle \frac{m^2}{\mu ^2}}\right)^{\epsilon /2}\left[{\displaystyle \frac{2}{\epsilon }}P_{gq}^{(0)}+w_{gq,\mathrm{B}}^{(1)}+\epsilon \overline{w}_{gq,\mathrm{B}}^{(1)}\right]+\widehat{a}_s^2S_\epsilon ^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)^\epsilon `$ (3.23)
$`\times [{\displaystyle \frac{1}{\epsilon ^2}}(+{\displaystyle \frac{3}{2}}P_{gg}^{(0)}P_{gq}^{(0)}+{\displaystyle \frac{1}{2}}P_{qq}^{(0)}P_{gq}^{(0)}5\beta _0P_{gq}^{(0)})`$
$`+{\displaystyle \frac{1}{\epsilon }}({\displaystyle \frac{1}{2}}P_{gq}^{(1)}6\beta _0w_{gq,\mathrm{B}}^{(1)}+{\displaystyle \frac{1}{2}}P_{qg}^{(0)}w_{q\overline{q},\mathrm{B}}^{(1)}+P_{gq}^{(0)}w_{gg,\mathrm{B}}^{(1)}`$
$`+(P_{gg}^{(0)}+P_{qq}^{(0)})w_{gq,\mathrm{B}}^{(1)})6\beta _0\overline{w}_{gq,\mathrm{B}}^{(1)}+{\displaystyle \frac{1}{2}}P_{qg}^{(0)}\overline{w}_{q\overline{q},\mathrm{B}}^{(1)}`$
$`+P_{gq}^{(0)}\overline{w}_{gg,\mathrm{B}}^{(1)}+(P_{gg}^{(0)}+P_{qq}^{(0)})\overline{w}_{gq,\mathrm{B}}^{(1)}+w_{gq,\mathrm{B}}^{(2)}],`$
$`z^1\widehat{W}_{q\overline{q},\mathrm{B}}`$ $`=`$ $`\widehat{a}_sS_\epsilon \left({\displaystyle \frac{m^2}{\mu ^2}}\right)^{\epsilon /2}[w_{q\overline{q},\mathrm{B}}^{(1)}+\epsilon \overline{w}_{q\overline{q},\mathrm{B}}^{(1)}]+\widehat{a}_s^2S_\epsilon ^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)^\epsilon [`$ (3.24)
$`{\displaystyle \frac{1}{\epsilon ^2}}P_{gq}^{(0)}P_{gq}^{(0)}+{\displaystyle \frac{1}{\epsilon }}(6\beta _0w_{q\overline{q}}^{(1)}+2P_{gq}^{(0)}w_{gq,\mathrm{B}}^{(1)}`$
$`+2P_{qq}^{(0)}w_{q\overline{q},\mathrm{B}}^{(1)})6\beta _0\overline{w}_{q\overline{q},\mathrm{B}}^{(1)}+2P_{gq}^{(0)}\overline{w}_{gq,\mathrm{B}}^{(1)}`$
$`+2P_{qq}^{(0)}\overline{w}_{q\overline{q},\mathrm{B}}^{(1)}+w_{q\overline{q},\mathrm{B}}^{(2)}],`$
$`z^1\widehat{W}_{q_1q_2,\mathrm{B}}`$ $`=`$ $`z^1\widehat{W}_{q_1\overline{q}_2,\mathrm{B}},`$ (3.25)
$`=`$ $`\widehat{a}_s^2S_\epsilon ^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)^\epsilon [{\displaystyle \frac{1}{\epsilon ^2}}P_{gq}^{(0)}P_{gq}^{(0)}+{\displaystyle \frac{2}{\epsilon }}P_{gq}^{(0)}w_{gq,\mathrm{B}}^{(1)}`$
$`+2P_{gq}^{(0)}\overline{w}_{gq,\mathrm{B}}^{(1)}+w_{q_1q_2,\mathrm{B}}^{(2)}],`$
where $``$ denotes the convolution symbol defined by

$`fg(z)={\displaystyle _0^1}𝑑z_1{\displaystyle _0^1}𝑑z_2\delta (zz_1z_2)f(z_1)g(z_2).`$ (3.26)
For identical quark-quark scattering we have the same formula as in Eq. (3.25) except that $`w_{qq,\mathrm{B}}^{(2)}w_{q_1q_2,\mathrm{B}}^{(2)}`$. The expressions above follow from the renormalization group equations. They are constructed in such a way that become finite after coupling constant renormalization, operator renormalization and mass factorization are carried out. The splitting functions $`P_{ab}(z)`$ and the coefficients $`w_{ab}(z)`$ with $`z=m^2/s`$ also occur in the coefficient functions given below except for the NLO terms $`\overline{w}_{ab}^{(1)}(z)`$, which are proportional to $`\epsilon `$. They are given by

$`\overline{w}_{gg,\mathrm{H}}^{(1)}`$ $`=`$ $`C_A[8𝒟_2(z)6\zeta (2)𝒟_0(z)+({\displaystyle \frac{8}{z}}16+8z8z^2)(\mathrm{ln}^2(1z)`$
$`{\displaystyle \frac{3}{4}}\zeta (2))+({\displaystyle \frac{8}{1z}}+{\displaystyle \frac{8}{z}}16+8z8z^2\left)\right({\displaystyle \frac{1}{4}}\mathrm{ln}^2z`$
$`\mathrm{ln}z\mathrm{ln}(1z)){\displaystyle \frac{22}{3}}{\displaystyle \frac{(1z)^3}{z}}(\mathrm{ln}(1z){\displaystyle \frac{1}{2}}\mathrm{ln}z){\displaystyle \frac{55}{3}}+{\displaystyle \frac{67}{9z}}`$
$`+{\displaystyle \frac{55}{3}}z{\displaystyle \frac{67}{9}}z^2+2\delta (1z)],`$
$`\overline{w}_{gq,\mathrm{H}}^{(1)}`$ $`=`$ $`C_F[({\displaystyle \frac{4}{z}}4+2z)(\mathrm{ln}^2(1z)\mathrm{ln}z\mathrm{ln}(1z)+{\displaystyle \frac{1}{4}}\mathrm{ln}^2z`$
$`{\displaystyle \frac{3}{4}}\zeta (2))+({\displaystyle \frac{3}{z}}+6z)(\mathrm{ln}(1z){\displaystyle \frac{1}{2}}\mathrm{ln}z)+{\displaystyle \frac{7}{2z}}5+{\displaystyle \frac{3}{2}}z],`$
$`\overline{w}_{q\overline{q},\mathrm{H}}^{(1)}`$ $`=`$ $`C_F^2\left[{\displaystyle \frac{8}{3}}{\displaystyle \frac{(1z)^3}{z}}\left(\mathrm{ln}(1z){\displaystyle \frac{1}{2}}\mathrm{ln}z+{\displaystyle \frac{1}{6}}\right)\right],`$ (3.27)
where $`𝒟_i`$ denotes the distribution

$`𝒟_i=\left({\displaystyle \frac{\mathrm{ln}^i(1z)}{1z}}\right)_+.`$ (3.28)
The difference between scalar and pseudo-scalar Higgs production is given by

$`\overline{w}_{gg,\mathrm{A}\mathrm{H}}^{(1)}`$ $`=`$ $`C_A\left[8(1z)14\delta (1z)\right],`$
$`\overline{w}_{gq,\mathrm{A}\mathrm{H}}^{(1)}`$ $`=`$ $`0,`$
$`\overline{w}_{q\overline{q},\mathrm{A}\mathrm{H}}^{(1)}`$ $`=`$ $`0.`$ (3.29)
To render the partonic cross sections finite one has first to perform coupling constant renormalization. This is done by replacing the bare coupling constant by the renormalized one see Eq. (2.29). Then one has to carry out operator renormalization which is achieved by multiplying the expressions in Eqs. (3.22)-(3.25) with the operator renormalization constants in Eqs. (2.25), (2.28). The remaining divergences, which are of collinear origin, are removed by mass factorization

$`z^1Z_O^2\widehat{W}_{ab,\mathrm{B}}({\displaystyle \frac{1}{\epsilon }},{\displaystyle \frac{m^2}{\mu ^2}})={\displaystyle \underset{c,d=q,\overline{q},g}{}}\mathrm{\Gamma }_{ca}\left({\displaystyle \frac{1}{\epsilon }}\right)\mathrm{\Gamma }_{db}\left({\displaystyle \frac{1}{\epsilon }}\right)z^1\mathrm{\Delta }_{cd,\mathrm{B}}\left({\displaystyle \frac{m^2}{\mu ^2}}\right),`$
where $`\mathrm{\Gamma }_{ca}(z)`$ denote the kernels containing the splitting functions which multiply the collinear divergences represented by the pole terms $`1/\epsilon ^k`$. Note that in the case of pseudo-scalar Higgs production the $`Z_O`$ (see Eq. (2.25)) in the above equation has to be replaced by $`Z_{11}`$ in Eq. (2.28). For the four different subprocesses the mass factorization relations in Eq. (3) become equal to

$`z^1Z_O^2\widehat{W}_{gg,\mathrm{B}}`$ $`=`$ $`\mathrm{\Gamma }_{gg}\mathrm{\Gamma }_{gg}z^1\mathrm{\Delta }_{gg,\mathrm{B}}+4\mathrm{\Gamma }_{gg}\mathrm{\Gamma }_{qg}z^1\mathrm{\Delta }_{gq,\mathrm{B}},`$
$`z^1Z_O^2\widehat{W}_{gq,\mathrm{B}}`$ $`=`$ $`\mathrm{\Gamma }_{qq}\mathrm{\Gamma }_{gg}z^1\mathrm{\Delta }_{gq,\mathrm{B}}+\mathrm{\Gamma }_{gq}\mathrm{\Gamma }_{gg}z^1\mathrm{\Delta }_{gg,\mathrm{B}}`$ (3.32)
$`+\mathrm{\Gamma }_{qq}\mathrm{\Gamma }_{qg}z^1\mathrm{\Delta }_{q\overline{q},\mathrm{B}},`$
$`z^1Z_O^2\widehat{W}_{q\overline{q},\mathrm{B}}`$ $`=`$ $`\mathrm{\Gamma }_{gq}\mathrm{\Gamma }_{gq}z^1\mathrm{\Delta }_{gg,\mathrm{B}}+\mathrm{\Gamma }_{qq}\mathrm{\Gamma }_{qq}z^1\mathrm{\Delta }_{q\overline{q},\mathrm{B}}`$ (3.33)
$`+2\mathrm{\Gamma }_{gq}\mathrm{\Gamma }_{qq}z^1\mathrm{\Delta }_{gq,\mathrm{B}},`$
$`z^1Z_O^2\widehat{W}_{q_1q_2,\mathrm{B}}`$ $`=`$ $`\mathrm{\Gamma }_{gq}\mathrm{\Gamma }_{gq}z^1\mathrm{\Delta }_{gg,\mathrm{B}}+2\mathrm{\Gamma }_{gq}\mathrm{\Gamma }_{qq}z^1\mathrm{\Delta }_{gq,\mathrm{B}}`$ (3.34)
$`+\mathrm{\Gamma }_{qq}\mathrm{\Gamma }_{qq}z^1\mathrm{\Delta }_{q_1q_2,\mathrm{B}},`$
where we have identified

$`\mathrm{\Gamma }_{qg}=\mathrm{\Gamma }_{\overline{q}g},\mathrm{\Gamma }_{gq}=\mathrm{\Gamma }_{g\overline{q}},`$
$`\widehat{W}_{gq,\mathrm{B}}=\widehat{W}_{g\overline{q},\mathrm{B}}=\widehat{W}_{qg,\mathrm{B}}=\widehat{W}_{\overline{q}g,\mathrm{B}},`$
$`\mathrm{\Delta }_{gq,\mathrm{B}}=\mathrm{\Delta }_{g\overline{q},\mathrm{B}}=\mathrm{\Delta }_{qg,\mathrm{B}}=\mathrm{\Delta }_{\overline{q}g,\mathrm{B}}.`$ (3.35)
Since we need the finite expressions up to order $`\alpha _s^2`$ it is sufficient to expand the kernels $`\mathrm{\Gamma }_{ca}`$ up to the following order in the renormalized coupling constant

$`\mathrm{\Gamma }_{qq}`$ $`=`$ $`\delta (1z)+a_s(\mu ^2)S_\epsilon \left[{\displaystyle \frac{1}{\epsilon }}P_{qq}^{(0)}\right],`$ (3.36)
$`\mathrm{\Gamma }_{qg}`$ $`=`$ $`a_s(\mu ^2)S_\epsilon \left[{\displaystyle \frac{1}{2\epsilon }}P_{qg}^{(0)}\right],`$ (3.37)
$`\mathrm{\Gamma }_{qq}^{\mathrm{PS}}`$ $`=`$ $`a_s^2(\mu ^2)S_\epsilon ^2\left[{\displaystyle \frac{1}{4\epsilon ^2}}P_{qg}^{(0)}P_{gq}^{(0)}+{\displaystyle \frac{1}{4\epsilon ^2}}P_{qq}^{(1),\mathrm{PS}}\right],`$ (3.38)
$`\mathrm{\Gamma }_{gq}`$ $`=`$ $`a_s(\mu ^2)S_\epsilon \left[{\displaystyle \frac{1}{\epsilon }}P_{gq}^{(0)}\right]+a_s^2(\mu ^2)S_\epsilon ^2[{\displaystyle \frac{1}{\epsilon ^2}}({\displaystyle \frac{1}{2}}P_{gq}^{(0)}P_{gg}^{(0)}`$ (3.39)
$`+{\displaystyle \frac{1}{2}}P_{qq}^{(0)}P_{gq}^{(0)}+\beta _0P_{gq}^{(0)})+{\displaystyle \frac{1}{2\epsilon }}P_{gq}^{(1)}],`$
$`\mathrm{\Gamma }_{gg}`$ $`=`$ $`\delta (1z)+a_s(\mu ^2)S_\epsilon \left[{\displaystyle \frac{1}{\epsilon }}P_{gg}^{(0)}\right]+a_s^2(\mu ^2)S_\epsilon ^2[{\displaystyle \frac{1}{\epsilon ^2}}({\displaystyle \frac{1}{2}}P_{gg}^{(0)}P_{gg}^{(0)}`$ (3.40)
$`+{\displaystyle \frac{1}{2}}P_{gq}^{(0)}P_{qg}^{(0)}+\beta _0P_{gg}^{(0)})+{\displaystyle \frac{1}{2\epsilon }}P_{gg}^{(1)}],`$
(PS denotes pure-singlet). After renormalization and mass factorization the coefficient functions have the following algebraic form

$`z^1\mathrm{\Delta }_{gg,\mathrm{B}}`$ $`=`$ $`\delta (1z)+a_s(\mu ^2)[P_{gg}^{(0)}\mathrm{ln}{\displaystyle \frac{m^2}{\mu ^2}}+w_{gg,\mathrm{B}}^{(1)}]+a_s^2(\mu ^2)[`$ (3.41)
$`\left\{{\displaystyle \frac{1}{2}}P_{gg}^{(0)}P_{gg}^{(0)}+{\displaystyle \frac{1}{4}}P_{gq}^{(0)}P_{qg}^{(0)}{\displaystyle \frac{5}{2}}\beta _0P_{gg}^{(0)}+3\beta _0^2\right\}\mathrm{ln}^2{\displaystyle \frac{m^2}{\mu ^2}}`$
$`+\{P_{gg}^{(1)}3\beta _0w_{gg,\mathrm{B}}^{(1)}+P_{gg}^{(0)}w_{gg,\mathrm{B}}^{(1)}+w_{gq,\mathrm{B}}^{(1)}P_{qg}^{(0)}`$
$`+z_\mathrm{B}\}\mathrm{ln}{\displaystyle \frac{m^2}{\mu ^2}}+w_{gg,\mathrm{B}}^{(2)}],`$
$`z^1\mathrm{\Delta }_{gq,\mathrm{B}}`$ $`=`$ $`a_s(\mu ^2)[{\displaystyle \frac{1}{2}}P_{gq}^{(0)}\mathrm{ln}{\displaystyle \frac{m^2}{\mu ^2}}+w_{gq,\mathrm{B}}^{(1)}]+a_s^2(\mu ^2)[\{{\displaystyle \frac{3}{8}}P_{gg}^{(0)}P_{gq}^{(0)}`$ (3.42)
$`+{\displaystyle \frac{1}{8}}P_{qq}^{(0)}P_{gq}^{(0)}{\displaystyle \frac{5}{4}}\beta _0P_{gq}^{(0)}\}\mathrm{ln}^2{\displaystyle \frac{m^2}{\mu ^2}}+\{{\displaystyle \frac{1}{2}}P_{gq}^{(1)}3\beta _0w_{gq,\mathrm{B}}^{(1)}`$
$`+{\displaystyle \frac{1}{4}}P_{qg}^{(0)}w_{q\overline{q},\mathrm{B}}^{(1)}+{\displaystyle \frac{1}{2}}P_{gq}^{(0)}w_{gg,\mathrm{B}}^{(1)}`$
$`+{\displaystyle \frac{1}{2}}(P_{gg}^{(0)}+P_{qq}^{(0)})w_{gq,\mathrm{B}}^{(1)}\}\mathrm{ln}{\displaystyle \frac{m^2}{\mu ^2}}+w_{gq,\mathrm{B}}^{(2)}],`$
$`z^1\mathrm{\Delta }_{q\overline{q},\mathrm{B}}`$ $`=`$ $`a_s(\mu ^2)\left[w_{q\overline{q},\mathrm{B}}^{(1)}\right]+a_s^2(\mu ^2)[{\displaystyle \frac{1}{4}}P_{gq}^{(0)}P_{gq}^{(0)}\mathrm{ln}^2{\displaystyle \frac{m^2}{\mu ^2}}`$ (3.43)
$`+\left\{3\beta _0w_{q\overline{q},\mathrm{B}}^{(1)}+P_{gq}^{(0)}w_{gq,\mathrm{B}}^{(1)}+P_{qq}^{(0)}w_{q\overline{q},\mathrm{B}}^{(1)}\right\}\mathrm{ln}{\displaystyle \frac{m^2}{\mu ^2}}`$
$`+w_{q\overline{q},\mathrm{B}}^{(2)}],`$
$`z^1\mathrm{\Delta }_{q_1q_2,\mathrm{B}}`$ $`=`$ $`z^1\mathrm{\Delta }_{q_1\overline{q}_2,\mathrm{B}},`$ (3.44)
$`=`$ $`a_s^2(\mu ^2)[{\displaystyle \frac{1}{4}}P_{gq}^{(0)}P_{gq}^{(0)}\mathrm{ln}^2{\displaystyle \frac{m^2}{\mu ^2}}+P_{gq}^{(0)}w_{gq,\mathrm{B}}^{(1)}\mathrm{ln}{\displaystyle \frac{m^2}{\mu ^2}}`$
$`+w_{q_1q_2,\mathrm{B}}^{(2)}].`$
In Appendix A and B we give the explicit expressions for the coefficient functions so that one can determine the coefficients $`P_{ab}^{(k)}`$ and $`w_{ab,\mathrm{B}}^{(k)}`$. Our results which are expressed in the colour factors $`C_A`$, $`C_F`$ and $`T_f`$ agree with those published in , ,  for $`N=3`$. In the representation of the coefficient functions above we have put the renormalization scale $`\mu _r`$ equal to the mass factorization scale $`\mu `$. If one wants to distinguish between both scales one can make the simple substitution

$`\alpha _s(\mu ^2)=\alpha _s(\mu _r^2)\left[1+{\displaystyle \frac{\alpha _s(\mu _r^2)}{4\pi }}\beta _0\mathrm{ln}{\displaystyle \frac{\mu _r^2}{\mu ^2}}\right].`$ (3.45)
## 4 Total cross sections for the process <br>$`p+pH+^{}X^{}`$

In this section we will present total cross sections (see Eq. (2.1)) for Higgs-boson production in proton-proton collisions at the LHC and in proton-anti-proton collisions at the TEVATRON. The cross section can be written in two ways (for the definitions see , )

$`\sigma _{\mathrm{tot}}={\displaystyle \frac{\pi G_\mathrm{B}^2}{8(N^21)}}{\displaystyle \underset{a,b=q,\overline{q},g}{}}{\displaystyle _x^1}𝑑y\stackrel{~}{\mathrm{\Phi }}_{ab}(y,\mu ^2)\mathrm{\Delta }_{ab}({\displaystyle \frac{x}{y}},{\displaystyle \frac{m^2}{\mu ^2}}),`$ (4.1)
where $`x=m^2/S`$ and $`\stackrel{~}{\mathrm{\Phi }}_{ab}`$ is the momentum fraction luminosity defined by

$`\stackrel{~}{\mathrm{\Phi }}_{ab}(y,\mu ^2)={\displaystyle _y^1}{\displaystyle \frac{du}{u}}f_a(u,\mu ^2)f_b({\displaystyle \frac{y}{u}},\mu ^2).`$ (4.2)
However Eq. (4.1) can be also cast in the form

$`\sigma _{\mathrm{tot}}={\displaystyle \frac{\pi G_\mathrm{B}^2}{8(N^21)}}{\displaystyle \underset{a,b=q,\overline{q},g}{}}x{\displaystyle _x^1}𝑑y\mathrm{\Phi }_{ab}(y,\mu ^2)\left({\displaystyle \frac{y}{x}}\mathrm{\Delta }_{ab}({\displaystyle \frac{x}{y}},{\displaystyle \frac{m^2}{\mu ^2}})\right),`$
$`\text{with}\mathrm{\Phi }_{ab}(y,\mu ^2)=y^1\stackrel{~}{\mathrm{\Phi }}_{ab}(y,\mu ^2),`$ (4.3)
where $`\mathrm{\Phi }_{ab}`$ denotes the parton luminosity. For the exact cross section the expressions in Eqs. (4.1) and (4) lead to the same results but if one wants to study the various approximations in the literature like the soft-plus-virtual (SV) gluon approximation then it makes a difference if $`\mathrm{\Delta }_{ab}(z)`$ or $`z^1\mathrm{\Delta }_{ab}(z)`$ will be expanded around $`z=1`$. This difference will show up in the less singular terms as we will show at the end of this section. In the subsequent part of this section we study the dependence of the cross section on the input parameters like the QCD scale $`\mathrm{\Lambda }`$, the renormalization/factorization scale $`\mu `$, the mass $`m`$ of the Higgs boson and the input parton densities. We are also interested in the theoretical uncertainty of the NNLO cross section. One uncertainty is due to the corrections which show up beyond NNLO. They can be guessed from the rate of convergence represented by the $`K`$-factor and the variation with respect to the scale $`\mu `$. We will show that the coefficient functions $`\mathrm{\Delta }_{ab}(z)`$ in Eq. (4.1) are dominated by the logarithmic terms of the type $`\mathrm{ln}^i(1z)`$. A determination of these logarithms is of importance to estimate the higher order corrections. Another uncertainty is caused by the input parton densities. Parton density sets differ from each other in several aspects. They are based on fitting data from different experiments. Furthermore the shapes of the densities at the specific input scale chosen for $`\mu `$ differ from each other. Finally these sets are based on different theoretical assumptions, for instance in the treatment of heavy flavours. In some sets the heavy flavour is treated as a massless particle which is described by a parton density similar to the treatment of the light quarks. In other sets the mass of the heavy flavour is considered to be on the same order of magnitude as the other large scales. Then the heavy flavour production is described by exact perturbation theory at a certain order so that threshold effects are fully taken into account. Although the coefficient functions in the effective Lagrangian approach are computed exactly in NNLO this is not the case for the parton densities because the exact three-loop splitting functions (anomalous dimensions) are not known yet. Until now only a finite number of moments are available (see ) which are used in  together with other constraints to approximate the splitting functions. These approximations are very reliable as long as $`y>10^4`$ in Eq. (4.1). The approximated splitting functions were used in  to obtain NNLO parton density sets. One of them called MRST02 (see Table 1) will be used in our paper. For the approximations in LO and NLO we use the sets in  and  respectively. For the LO, NLO and NNLO plots we employ the one-, two-, and three-loop asymptotic forms of the running coupling constant as given in Eq. (3) of . In order to make a comparison with other parametrizations for the parton densities we also choose the sets made by CTEQ  and GRV . However these sets do not contain NNLO versions so that we can present the NLO results only. All sets are listed in Table 1 <sup>2</sup><sup>2</sup>2Notice that according to the prescription of the CTEQ group  also the LO corrected quantities have to be computed with the NLO running coupling constant. together with the corresponding QCD scale $`\mathrm{\Lambda }_{n_f}`$ determined for $`n_f=5`$. The same number of flavours is also chosen for the coefficient functions.

Notice that the GRV sets do not contain densities for charm and bottom quarks. For simplicity the factorization scale is set equal to the renormalization scale $`\mu _r`$. For our plots we take $`\mu =m`$ and use the MRST sets for the LO, NLO and NNLO computations (see Table 1) unless mentioned otherwise. Here we want to emphasize that the magnitude of the $`\sigma _{\mathrm{tot}}`$ is extremely sensitive to the renormalization scale because the effective coupling constant $`G_\mathrm{B}\alpha _s(\mu _r)`$, which implies that $`\sigma ^{\mathrm{LO}}\alpha _s^2`$. The sensitivity to the factorization scale is much smaller.

For the computation of the effective coupling constant $`G_\mathrm{B}`$ in Eq. (2) we choose the top quark mass $`m_t=173.4\mathrm{GeV}/\mathrm{c}^2`$ and the Fermi constant $`G_F=1.16639\mathrm{GeV}^2=4541.68\mathrm{pb}`$. In this paper we will only study scalar Higgs production and omit pseudo-scalar Higgs production. The cross sections of the latter are about $`9/4\mathrm{cot}^2\beta `$ larger than those for the standard model Higgs so that all our conclusions also apply to pseudo-scalar Higgs production. We give results for both proton-proton collisions at the centre-of-mass energy $`\sqrt{S}=14\mathrm{TeV}`$ and proton-anti-proton collisions at $`\sqrt{S}=2\mathrm{TeV}`$.

In Fig. 1 we have plotted the contributions coming from the various subprocesses up to NLO. From this figure we infer that $`\sigma _{\mathrm{tot}}`$ is completely dominated by the $`gg`$ reaction whereas the other contributions are down by two ($`gq+g\overline{q}`$-subprocess) or even by three ($`q\overline{q}`$-subprocess) orders of magnitude. This picture is different from the behaviour of the transverse momentum $`p_T`$\- and rapidity $`y`$-distributions in  where the $`gq+g\overline{q}`$-reaction amounted to about $`1/3`$ of that for the $`gg`$-reaction when $`p_T>30\mathrm{GeV}/\mathrm{c}`$. This is because $`\sigma _{\mathrm{tot}}`$ receives all its contribution from the small $`p_T`$-region ($`p_T<30\mathrm{GeV}/\mathrm{c}`$) where the $`gg`$ subprocess overwhelms all other reactions. This picture is not changed when we study the cross section in NNLO (see Fig. 2) where new subprocesses contribute given by the $`qq`$, $`q\overline{q}`$ and $`\overline{q}\overline{q}`$ reactions. Notice that the incoming (anti-) quarks can be identical and non-identical. The relative order of magnitude is the same as in NLO and the new reactions are down by three orders of magnitude. In Fig. 3 we have plotted $`\sigma _{\mathrm{tot}}`$ in LO, NLO and NNLO as a function of $`m`$. The cross section decreases when $`m`$ increases due to a reduction in the available phase space. In NNLO the curve shows a steeper decrease than in NLO and LO. Furthermore the growth of the cross section slows down while going to higher orders. To show more clearly how the cross sections change with respect to $`m`$ we have put them in Table 2. We have also shown them for the TEVATRON (proton-anti-proton collisions at $`\sqrt{S}=2\mathrm{TeV}`$) in Table 3 where they are appreciably smaller than in the case of the LHC.

The NNLO cross sections in Tables 2, 3 have some theoretical errors which can be estimated in various ways. First one can study their variation with respect to the scale $`\mu `$. This can be achieved by computing the ratio

$`N\left({\displaystyle \frac{\mu }{\mu _0}}\right)={\displaystyle \frac{\sigma _{\mathrm{tot}}(\mu )}{\sigma _{\mathrm{tot}}(\mu _0)}},`$ (4.4)
where $`\mu _0=m`$ and $`\mu `$ is varied in the range $`0.1<\mu /\mu _0<10`$. In Fig. 4 a plot of $`N`$ is shown for $`m=100\mathrm{GeV}/\mathrm{c}^2`$. Here one observes a clear improvement while going from LO to NNLO. In particular the curve for NNLO is flatter than that for NLO. The improvement becomes even better when we choose a larger Higgs mass see e.g. Fig. 5 ($`m=200\mathrm{GeV}/\mathrm{c}^2`$). However for still larger masses there is hardly any improvement anymore (see Fig. 6 where $`m=300\mathrm{GeV}/\mathrm{c}^2`$). This is in contrast to the NNLO corrected cross section for vector boson production in  which is insensitive to the choice of scale chosen in the range above.

A second way to study the reliability of our prediction is to study the rate of convergence of the perturbation series, which is represented by the $`K`$-factor. We choose the following definitions

$`K^{\mathrm{NLO}}={\displaystyle \frac{\sigma _{\mathrm{tot}}^{\mathrm{NLO}}}{\sigma _{\mathrm{tot}}^{\mathrm{LO}}}},K^{\mathrm{NNLO}}={\displaystyle \frac{\sigma _{\mathrm{tot}}^{\mathrm{NNLO}}}{\sigma _{\mathrm{tot}}^{\mathrm{NLO}}}}.`$ (4.5)
Notice that in the above expression the definition of $`K^{\mathrm{NNLO}}`$ differs from the usual one which is given by $`K^{\mathrm{NNLO}}=\sigma _{\mathrm{tot}}^{\mathrm{NNLO}}/\sigma _{\mathrm{tot}}^{\mathrm{LO}}`$. We have chosen this definition because the rate of convergence can be shown in a better way. In Fig. 7 one observes that both $`K`$-factors vary slowly as $`m`$ increases. Moreover there is considerable improvement in the rate of convergence if one goes from NLO to NNLO. At $`m=100\mathrm{GeV}/\mathrm{c}^2`$ we have $`K^{\mathrm{NLO}}1.74`$ whereas $`K^{\mathrm{NNLO}}1.15`$. Still the corrections for Higgs boson production are larger than those obtained from $`Z`$-boson production at the LHC where one gets $`K^{\mathrm{NLO}}1.22`$ and $`K^{\mathrm{NNLO}}0.95`$ (see ). In Fig. 8 we compare the $`K`$-factor in NLO obtained from the MRST with the results computed by using the parton density sets given by the CTEQ  and GRV  collaborations (see Table 1). From this figure we infer that the $`K`$-factor computed using GRV98 and CTEQ6 is larger than the one obtained from MRST01(NLO) at $`m=100\mathrm{GeV}/\mathrm{c}^2`$. However at large $`m`$ there is a cross over and we observe $`K^{\mathrm{CTEQ}}>K^{\mathrm{MRST}}>K^{\mathrm{GRV98}}`$ . The large result shown by the CTEQ6-set is due to the behaviour of their gluon density in LO. In NLO this density behaves in the same way as those in the MRST and GRV sets. We also investigated older sets provided by the CTEQ-collaboration, namely CTEQ5 and CTEQ6. It turns out that all these LO CTEQ sets lead to the same large value for $`K^{\mathrm{NLO}}`$ as shown by Fig. 8. A study of $`K^{\mathrm{NLO}}`$, using the most recent parton densities, was also made in  but now for the TEVATRON.

After the K-factor we now make a comparison between $`\sigma _{\mathrm{tot}}`$ computed from the various sets which leads to the third uncertainty in the theoretical prediction. For that purpose we plot the ratios

$`R^{\mathrm{CTEQ}}={\displaystyle \frac{\sigma _{\mathrm{tot}}^{\mathrm{CTEQ}}}{\sigma _{\mathrm{tot}}^{\mathrm{MRST}}}},R^{\mathrm{GRV}}={\displaystyle \frac{\sigma _{\mathrm{tot}}^{\mathrm{GRV}}}{\sigma _{\mathrm{tot}}^{\mathrm{MRST}}}},`$ (4.6)
which are shown in NLO in Fig. 9. In this figure we see that the CTEQ6 set leads to the same cross section as the one obtained by MRST01. The cross section given by the GRV98 is $`20\%`$ above the MRST01 result and the ratio decreases with increasing $`m`$ to $`5\%`$ above at $`m=300\mathrm{GeV}/\mathrm{c}^2`$. The reason that the CTEQ6 result is closer to the MRST sets can be attributed to the fact that their parton density sets are more recent and are constructed to fit the same experiments. Notice that a study of the ratios in Eq. (4.6) was also made for the TEVATRON in  using the most recent parton densities.

The factorization scale dependence in Figs. 4-6 and the $`K`$-factor in Fig. 7 of the NNLO cross section can be used to give an error estimate of the latter for the LHC. Since $`(\sigma ^{\mathrm{NNLO}}\sigma ^{\mathrm{NLO}})/\sigma ^{\mathrm{NLO}}=K^{\mathrm{NNLO}}1`$ the error on $`\sigma ^{\mathrm{NNLO}}`$ is equal to $`15\%`$, $`19\%`$ and $`21\%`$ for $`m=100\mathrm{GeV}/\mathrm{c}^2`$, $`m=200\mathrm{GeV}/\mathrm{c}^2`$ and $`m=300\mathrm{GeV}/\mathrm{c}^2`$ respectively (see Fig. 7). This is corroborated by the scale variation if we choose the range $`0.25<\mu /\mu _0<4.0`$. For $`\mu /\mu _0=0.25`$ the errors above become $`18\%`$, $`15\%`$ and $`14\%`$ respectively. In the case $`\mu /\mu _0=0.25`$ they are given by $`16\%`$ for all $`m`$ given above. Notice that the error due the $`K`$-factor increases when $`m`$ gets larger whereas the opposite happens if the error estimate is derived from the mass factorization/renormalisation scale. On top of this we have to add the error due to the chosen parton density. For this we choose the NLO plots in Fig. 9 where the cross section predicted by CTEQ6 hardly differs from the one given by MRST01 but the GRV98 plot deviates from the latter by $`20\%`$ for small $`m`$ to $`5\%`$ for large $`m`$.

An important feature of the Higgs cross section, which is very often emphasised in the literature, is that the integral in Eq. (4.1) is dominated by $`yx=m^2/S`$. Since $`x`$ is small the $`gg`$ part of $`\sigma _{\mathrm{tot}}`$ dominates the contributions coming from the other partonic subprocesses because of the steep rise of the gluon-gluon luminosity in Eq. (4.2) at small $`y`$. The importance of the small $`y`$-region is revealed when we impose an upper cut on the integral in Eq. (4.1). Hence we compute

$`\sigma _{\mathrm{tot}}(x_{\mathrm{max}})`$ $`=`$ $`{\displaystyle \frac{\pi G_\mathrm{B}^2}{8(N^21)}}[\stackrel{~}{\mathrm{\Phi }}_{gg}(x,\mu ^2)+{\displaystyle \underset{a,b=q,g}{}}{\displaystyle \underset{i=1}{\overset{\mathrm{}}{}}}a_s^i(\mu ^2){\displaystyle _x^{x_{\mathrm{max}}}}dy\stackrel{~}{\mathrm{\Phi }}_{ab}(y,\mu ^2)`$ (4.7)
$`\times \mathrm{\Delta }_{ab}^{(i)}({\displaystyle \frac{x}{y}},{\displaystyle \frac{m^2}{\mu ^2}})],`$
and plot the ratio

$`R(x_{\mathrm{max}})={\displaystyle \frac{\sigma _{\mathrm{tot}}(x_{\mathrm{max}})}{\sigma _{\mathrm{tot}}(1)}},\sigma _{\mathrm{tot}}(1)=\sigma _{\mathrm{tot}}^{\mathrm{EXACT}},`$ (4.8)
which is shown for NLO and NNLO in Fig. 10a for the choice $`x_{\mathrm{max}}=5x`$. The figure reveals that more than 95 $`\%`$ of the cross section comes from the integration region $`xy5x`$ where $`x<5.1\times 10^5`$. When $`m`$ gets larger $`R(x_{\mathrm{max}})`$ will increase because the available phase space for Higgs boson production decreases. The latter also happens when the centre-of-mass energy decreases like in the case for the TEVATRON where $`\sqrt{S}=2\mathrm{TeV}`$ (Here $`x<2.5\times 10^3`$) see Fig. 10b. (Note Figs. 10a and 10b have different scales.) In the latter case more than 99 $`\%`$ of the cross section receives its contribution from $`xy5x`$ which becomes 100 $`\%`$ when $`m=300\mathrm{GeV}/\mathrm{c}^2`$.

Finally we study the validity of the soft-plus-virtual gluon approximation including subleading terms. The result of this approximation depends on the definitions for the total cross section given in Eq. (4.1) and Eq. (4) as has been pointed out in . If we follow the definition in , , where one choses the expression in Eq. (4.1), the sof-plus-virtual coefficient function becomes

$`\mathrm{\Delta }_{ab}^{\mathrm{SV}}(z,{\displaystyle \frac{m^2}{\mu ^2}})`$ $`=`$ $`\delta (1z)+{\displaystyle \underset{i=1}{\overset{\mathrm{}}{}}}a_s^i(\mu ^2)\left[{\displaystyle \underset{j=0}{\overset{2i1}{}}}a_{i,j}𝒟_j+b_i\delta (1z)\right].`$ (4.9)
where the distributions $`𝒟_i(z)`$ are defined in Eq. (3.28). Furthermore the logarithmic terms $`\mathrm{ln}^j(1z)`$ including the constants ($`j=0`$) are important. The latter are denoted by

$`\mathrm{\Delta }_{ab}^\mathrm{L}(z,{\displaystyle \frac{m^2}{\mu ^2}})={\displaystyle \underset{i=1}{\overset{\mathrm{}}{}}}a_s^i(\mu ^2){\displaystyle \underset{j=0}{\overset{2i1}{}}}c_{i,j}\mathrm{ln}^j(1z).`$ (4.10)
Notice that the coefficients $`a_{i,j},b_i,c_{i,j}`$ also contain terms in $`\mathrm{ln}m^2/\mu ^2`$. The soft-plus-virtual gluon contributions only emerge from the $`gg`$-subprocess and they can be found in Eqs. (A.2) and (A.6). This subprocess also leads to the logarithms of the type $`c_{i,j}\mathrm{ln}^j(1z)`$. The latter also show up in the $`gq+g\overline{q}`$-channel but they are absent in the $`q\overline{q}`$ and $`qq`$ subprocesses provided one neglects the terms $`(1z)^k`$ for $`k>0`$. The coefficients $`c_{i,j}`$ can be found in Eqs. (A.26)-(A.32) (see also ) (scalar Higgs) and Eqs. (B.21), (B.22) (pseudo-scalar Higgs). If we use the definition of the cross section in Eq. (4) advocated in ,  the expansion of $`z^1\mathrm{\Delta }_{ab}^{\mathrm{SV}}`$ will not differ from the one shown in Eq. (4.9) but $`\mathrm{\Delta }_{ab}^\mathrm{L}`$ will change due to contribution of subleading logarithms and it becomes equal to

$`\mathrm{\Delta }_{ab}^\mathrm{L}(z,{\displaystyle \frac{m^2}{\mu ^2}})={\displaystyle \underset{i=1}{\overset{\mathrm{}}{}}}a_s^i(\mu ^2){\displaystyle \underset{j=0}{\overset{2i1}{}}}\left(a_{i,j}+c_{i,j}\right)\mathrm{ln}^j(1z).`$ (4.11)
Because of the different luminosities $`\stackrel{~}{\mathrm{\Phi }}(z)`$ and $`\mathrm{\Phi }(z)`$ in Eqs. (4.1) and (4) respectively the value of the soft-plus-virtual cross section $`\sigma _{\mathrm{tot}}^{\mathrm{SV}}`$ will change. This will be only partially compensated by $`\sigma _{\mathrm{tot}}^\mathrm{L}`$ coming from Eqs. (4.10) and (4.11). In fact the two results for $`\sigma _{\mathrm{tot}}^{\mathrm{SVL}}`$ coming from definitions (4.1) and (4) will differ by terms in $`(1z)^k\mathrm{ln}^j(1z)`$ for $`k1`$ convoluted by one of the luminosities. In order to study the effect of the two definitions on the approximation we will plot

$`R^{\mathrm{SV}}={\displaystyle \frac{\sigma _{\mathrm{tot}}^{\mathrm{SV}}}{\sigma _{\mathrm{tot}}}},R^{\mathrm{SVL}}={\displaystyle \frac{\sigma _{\mathrm{tot}}^{\mathrm{SVL}}}{\sigma _{\mathrm{tot}}}}.`$ (4.12)
The same study is made in , ,  and . However in these papers the complete result for the NNLO coefficient function was not known yet. Therefore one could only study the leading logarithmic terms in Eqs. (4.10), (4.11) which are of collinear origin. In Fig. 11a and Fig. 11b we have shown the ratios above for the LHC in NLO and NNLO respectively. The figures reveal that $`\sigma _{\mathrm{tot}}^{\mathrm{SV}}`$ using Eq. (4) gives a much better approximation than Eq. (4.1). This also holds if we include the logarithms in Eqs. (4.10) and (4.11) denoted by $`\sigma _{\mathrm{tot}}^{\mathrm{SVL}}`$. The latter overestimates the exact cross section a little bit. Further we observe that both SV approximations become worse in NNLO with respect to NLO. In Figs. 12a,b we have performed the same analysis but now for the TEVATRON. The features are the same as in Figs. 11a,b for LHC except that for the TEVATRON the SV approximations derived from Eq. (4.1) and Eq. (4) become better. This is no surprise because for the TEVATRON the Higgs mass is much closer to the boundary of phase space. This also explains why at larger masses the SV approximation improves. From this study we infer that the soft-plus-virtual approximation works much better for Eq. (4) than for Eq. (4.1). The reason is that the cross section is dominated by the small $`x`$-region (see Eq. (4.7) and Fig. 10). By choosing the gluon luminosity in Eq. (4) this region will become even more important than the gluon luminosity chosen in Eq. (4.2) which is used in Eq. (4.7). For NLO this was already shown in  but it is now also confirmed in NNLO. Further this approximation becomes worse in NNLO which holds for both definitions. In spite of this drawback one should still choose Eq. (4) to perform the soft-plus-virtual gluon resummation (see ) in order to get a better estimate of the cross section. We also noticed that the contribution of logarithmic terms in Eqs. (4.10), Eqs. (4.11) is substantial. The determination of the coefficient of the leading term, which is of collinear origin, is done in , . A systematic approach to determine $`c_{i,j}`$ can be found in .

Summarising our results we have recalculated the NNLO corrections to the total cross section for (pseudo-) scalar Higgs production using a different method than those presented in  and . In particular our approach differs from the one given in  since we directly evaluated the multi-particle phase space integrals without fitting the coefficient functions to an expansion in the terms $`(1z)^k\mathrm{ln}^l(1z)`$. Further we presented the radiative parts of the coefficient functions for general colour factors and checked that for $`N=3`$ our answers are in agreement with the literature. It turns out that the total cross section is almost completely determined by the $`gg`$-subprocess which is in contrast to the differential cross section where also the $`gq+g\overline{q}`$-channel gives a considerable contribution. We studied the $`K`$-factors and the dependence of the cross sections on the chosen mass factorization/renormalization scale and the adopted parton density set from which one can infer an error estimate of the cross section. If we exclude the GRV98 parton density the error can be mainly attributed to the missing higher order contributions to the coefficient functions which we estimate to lie between $`14\%`$ and $`21\%`$. Finally we studied the validity of the soft-plus-virtual gluon approximation. Depending on the definition of the total cross section (Eq. (4.1) versus Eq. (4)) this approximation is excellent in NLO but it becomes less good in NNLO but Eq. (4) is the best way to resum the soft-plus-virtual corrections.
Acknowledgement: V. Ravindran would like to thank C. Anastasiou and K. Melnikov for providing him with the unpublished results. In particular he could check that we had the same partonic cross section for $`g+gg+g+\mathrm{H}(\mathrm{A})`$ prior to renormalization and mass factorization were carried out.

## Appendix A

In this Appendix we present the coefficient functions for scalar Higgs boson production. In the case of the $`gg`$ subprocess one can split the corresponding coefficient function into

$`\mathrm{\Delta }_{gg,\mathrm{B}}(z,{\displaystyle \frac{m^2}{\mu ^2}})=\mathrm{\Delta }_{gg,\mathrm{B}}^{\mathrm{S}+\mathrm{V}}(z,{\displaystyle \frac{m^2}{\mu ^2}})+\mathrm{\Delta }_{gg,\mathrm{B}}^\mathrm{H}(z,{\displaystyle \frac{m^2}{\mu ^2}}),`$ (A.1)
where the superscripts $`S+V`$ and $`H`$ denote the soft-plus-virtual and hard gluon parts respectively. The former contains the distributions $`\delta (1z)`$ and $`𝒟_i(z)`$ where the latter is defined in Eq.(3.28). We find

$`\mathrm{\Delta }_{gg,\mathrm{H}}^{(1),\mathrm{S}+\mathrm{V}}=C_A\left[8𝒟_0(z)\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+16𝒟_1(z)+8\zeta (2)\delta (1z)\right],`$ (A.2)
$`\mathrm{\Delta }_{gg,\mathrm{H}}^{(1)}`$ $`=`$ $`C_A[\{16z+8z^28z^3\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+(32z+16z^216z^3)`$ (A.3)
$`\times (\mathrm{ln}(1z){\displaystyle \frac{1}{2}}\mathrm{ln}z){\displaystyle \frac{8}{1z}}\mathrm{ln}z{\displaystyle \frac{22}{3}}(1z)^3],`$
$`\mathrm{\Delta }_{gq,\mathrm{H}}^{(1)}`$ $`=`$ $`C_F[\{44z+2z^2\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+(88z+4z^2)(\mathrm{ln}(1z)`$ (A.4)
$`{\displaystyle \frac{1}{2}}\mathrm{ln}z)3+6zz^2],`$
$`\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(1)}`$ $`=`$ $`C_F^2\left[{\displaystyle \frac{8}{3}}(1z)^3\right],`$ (A.5)
$`\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{S}+\mathrm{V}}`$ $`=`$ $`C_A^2[\{64𝒟_1(z){\displaystyle \frac{44}{3}}𝒟_0(z)32\zeta (2)\delta (1z)\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (A.6)
$`+\{192𝒟_2(z){\displaystyle \frac{176}{3}}𝒟_1(z)+({\displaystyle \frac{536}{9}}80\zeta (2))𝒟_0(z)`$
$`+\delta (1z)(24{\displaystyle \frac{88}{3}}\zeta (2)+152\zeta (3))\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+128𝒟_3(z){\displaystyle \frac{176}{3}}𝒟_2(z)+\left({\displaystyle \frac{1072}{9}}160\zeta (2)\right)𝒟_1(z)`$
$`+\left({\displaystyle \frac{1616}{27}}+{\displaystyle \frac{176}{3}}\zeta (2)+312\zeta (3)\right)𝒟_0(z)`$
$`+\delta (1z)(93+{\displaystyle \frac{536}{9}}\zeta (2){\displaystyle \frac{220}{3}}\zeta (3){\displaystyle \frac{4}{5}}\zeta ^2(2))]`$
$`+n_fT_fC_A[\left\{{\displaystyle \frac{16}{3}}𝒟_0(z)\right\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+\{{\displaystyle \frac{64}{3}}𝒟_1(z){\displaystyle \frac{160}{9}}𝒟_0(z)`$
$`+\delta (1z)(16+{\displaystyle \frac{32}{3}}\zeta (2))\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})+{\displaystyle \frac{64}{3}}𝒟_2(z){\displaystyle \frac{320}{9}}𝒟_1(z)`$
$`+({\displaystyle \frac{448}{27}}{\displaystyle \frac{64}{3}}\zeta (2))𝒟_0(z)+\delta (1z)({\displaystyle \frac{2024}{27}}{\displaystyle \frac{160}{9}}\zeta (2)`$
$`+{\displaystyle \frac{80}{9}}\zeta (3))]`$
$`+n_fT_fC_F\left[8\delta (1z)\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+4\delta (1z)\right].`$
Notice that the colour decomposition of the $`\delta (1z)`$ term into $`n_fT_fC_A`$ and $`n_fT_fC_F`$ is arbitrary due to our ignorance of the same parts appearing in the two-loop virtual corrections computed in . They are only correct when $`N=3`$.

$`\mathrm{\Delta }_{gg,\mathrm{H}}^{(2)}=C_A^2\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}^2}+n_fT_fC_A\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{T}_\mathrm{f}}+n_fT_fC_F\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}\mathrm{T}_\mathrm{f}},`$ (A.7)
$`\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}^2}`$ $`=`$ $`[\{(128z+64z^264z^3)\mathrm{ln}(1z)+(96z^2+32z^3`$ (A.8)
$`{\displaystyle \frac{32}{1z}})\mathrm{ln}z{\displaystyle \frac{352}{3}}+{\displaystyle \frac{376}{3}}z{\displaystyle \frac{332}{3}}z^2+132z^3\left\}\mathrm{ln}^2\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+\{(64z+32z^2+32z^3+{\displaystyle \frac{32}{1+z}})(\mathrm{Li}_2(z)+\mathrm{ln}z\mathrm{ln}(1+z))`$
$`\left(256z+256z^2\right)\mathrm{Li}_2(1z)+\left(384z+192z^2192z^3\right)`$
$`\times \mathrm{ln}^2(1z)+\left(192z480z^2+224z^3{\displaystyle \frac{224}{1z}}\right)`$
$`\times \mathrm{ln}z\mathrm{ln}(1z)+\left(128z^248z^3+{\displaystyle \frac{40}{1z}}{\displaystyle \frac{8}{1+z}}\right)\mathrm{ln}^2z`$
$`+(528+{\displaystyle \frac{2032}{3}}z{\displaystyle \frac{1856}{3}}z^2+{\displaystyle \frac{1760}{3}}z^3)\mathrm{ln}(1z)+(176`$
$`{\displaystyle \frac{904}{3}}z+{\displaystyle \frac{1328}{3}}z^2{\displaystyle \frac{1496}{3}}z^3+{\displaystyle \frac{88}{3}}{\displaystyle \frac{1}{1z}})\mathrm{ln}z+(192z64z^2`$
$`+96z^3+{\displaystyle \frac{16}{1+z}})\zeta (2)+{\displaystyle \frac{1262}{3}}{\displaystyle \frac{3562}{9}}z+{\displaystyle \frac{3026}{9}}z^2{\displaystyle \frac{4322}{9}}z^3\}`$
$`\times \mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$
$`+(128z+64z^2+64z^3+{\displaystyle \frac{64}{1+z}})(\mathrm{ln}z\mathrm{ln}(1z)\mathrm{ln}(1+z)`$
$`+\mathrm{ln}(1z)\mathrm{Li}_2(z)+\mathrm{Li}_3({\displaystyle \frac{1z}{1+z}})\mathrm{Li}_3\left({\displaystyle \frac{1z}{1+z}}\right))`$
$`+(24+8z32z^232z^3+{\displaystyle \frac{16}{1+z}})(2\mathrm{S}_{1,2}(z)`$
$`+\zeta (2)\mathrm{ln}(1+z)+2\mathrm{ln}(1+z)\mathrm{Li}_2(z)+\mathrm{ln}z\mathrm{ln}^2(1+z))`$
$`+(8+24z64z^3+{\displaystyle \frac{16}{1+z}})\mathrm{Li}_3(z)+(16208z`$
$`624z^2+112z^3{\displaystyle \frac{144}{1z}}{\displaystyle \frac{64}{1+z}})\mathrm{S}_{1,2}(1z)+(16`$
$`+672z+536z^2+88z^3{\displaystyle \frac{8}{1z}}+{\displaystyle \frac{64}{1+z}})\mathrm{Li}_3(1z)+(24`$
$`104z16z^2+16z^3{\displaystyle \frac{64}{1+z}})\mathrm{ln}z\mathrm{Li}_2(z)+(8+8z`$
$`88z^2+8z^3{\displaystyle \frac{56}{1z}}{\displaystyle \frac{48}{1+z}})\mathrm{ln}z\mathrm{Li}_2(1z)`$
$`512z(1+z)\mathrm{ln}(1z)\mathrm{Li}_2(1z)+(240z504z^2`$
$`+248z^3{\displaystyle \frac{248}{1z}})\mathrm{ln}z\mathrm{ln}^2(1z)+(96z+304z^2144z^3`$
$`+{\displaystyle \frac{128}{1z}}{\displaystyle \frac{16}{1+z}})\mathrm{ln}^2z\mathrm{ln}(1z)+(2092z16z^216z^3`$
$`{\displaystyle \frac{56}{1+z}})\mathrm{ln}^2z\mathrm{ln}(1+z)+(256z+128z^2128z^3)`$
$`\times \mathrm{ln}^3(1z)+({\displaystyle \frac{8}{3}}+{\displaystyle \frac{16}{3}}z{\displaystyle \frac{136}{3}}z^2+24z^3{\displaystyle \frac{56}{3}}{\displaystyle \frac{1}{1z}}`$
$`+{\displaystyle \frac{16}{3}}{\displaystyle \frac{1}{1+z}})\mathrm{ln}^3z+(384z128z^2+192z^3+{\displaystyle \frac{32}{1+z}})`$
$`\times \zeta (2)\mathrm{ln}(1z)+(872z+392z^2168z^3+{\displaystyle \frac{144}{1z}}`$
$`{\displaystyle \frac{24}{1+z}})\zeta (2)\mathrm{ln}z+({\displaystyle \frac{112}{3}}32z+56z^2+{\displaystyle \frac{88}{3}}z^3\left)\right(\mathrm{Li}_2(z)`$
$`+\mathrm{ln}z\mathrm{ln}(1+z))+({\displaystyle \frac{1112}{3}}+{\displaystyle \frac{464}{3}}z+{\displaystyle \frac{548}{3}}z^2{\displaystyle \frac{572}{3}}z^3`$
$`{\displaystyle \frac{44}{3}}{\displaystyle \frac{1}{1z}})\mathrm{Li}_2(1z)+({\displaystyle \frac{1232}{3}}{\displaystyle \frac{2144}{3}}z+{\displaystyle \frac{2848}{3}}z^2{\displaystyle \frac{2992}{3}}z^3`$
$`+{\displaystyle \frac{176}{3}}{\displaystyle \frac{1}{1z}})\mathrm{ln}z\mathrm{ln}(1z)+(528+{\displaystyle \frac{2032}{3}}z{\displaystyle \frac{1856}{3}}z^2`$
$`+{\displaystyle \frac{1760}{3}}z^3)\mathrm{ln}^2(1z)+({\displaystyle \frac{286}{3}}+164z302z^2+{\displaystyle \frac{748}{3}}z^3`$
$`{\displaystyle \frac{22}{3}}{\displaystyle \frac{1}{1z}})\mathrm{ln}^2z+({\displaystyle \frac{8108}{9}}{\displaystyle \frac{8636}{9}}z+{\displaystyle \frac{7576}{9}}z^21020z^3)`$
$`\times \mathrm{ln}(1z)+({\displaystyle \frac{3602}{9}}+{\displaystyle \frac{2188}{3}}z{\displaystyle \frac{1582}{3}}z^2+{\displaystyle \frac{8108}{9}}z^3`$
$`{\displaystyle \frac{536}{9}}{\displaystyle \frac{1}{1z}})\mathrm{ln}z+(608z+320z^2352z^3+{\displaystyle \frac{8}{1+z}})\zeta (3)`$
$`+({\displaystyle \frac{1220}{3}}{\displaystyle \frac{1156}{3}}z+{\displaystyle \frac{1016}{3}}z^2{\displaystyle \frac{1408}{3}}z^3)\zeta (2){\displaystyle \frac{16705}{27}}`$
$`+{\displaystyle \frac{18523}{27}}z{\displaystyle \frac{19751}{27}}z^2+{\displaystyle \frac{7055}{9}}z^3],`$
$`\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{T}_\mathrm{f}}`$ $`=`$ $`[\{{\displaystyle \frac{32}{3}}z+{\displaystyle \frac{16}{3}}z^2{\displaystyle \frac{16}{3}}z^3\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (A.9)
$`+\{({\displaystyle \frac{128}{3}}z+{\displaystyle \frac{64}{3}}z^2{\displaystyle \frac{64}{3}}z^3)\mathrm{ln}(1z)+({\displaystyle \frac{32}{3}}z{\displaystyle \frac{64}{3}}z^2`$
$`+{\displaystyle \frac{32}{3}}z^3{\displaystyle \frac{32}{3}}{\displaystyle \frac{1}{1z}})\mathrm{ln}z{\displaystyle \frac{296}{9}}+{\displaystyle \frac{728}{9}}z{\displaystyle \frac{568}{9}}z^2`$
$`+{\displaystyle \frac{152}{3}}z^3\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+(16+16z8z^2)(\mathrm{S}_{1,2}(1z)+\mathrm{Li}_3(1z))+({\displaystyle \frac{16}{3}}`$
$`{\displaystyle \frac{160}{3}}z+{\displaystyle \frac{8}{3}}z^2{\displaystyle \frac{32}{3}}z^3+{\displaystyle \frac{16}{3}}{\displaystyle \frac{1}{1z}})\mathrm{Li}_2(1z)+({\displaystyle \frac{128}{3}}z`$
$`+{\displaystyle \frac{64}{3}}z^2{\displaystyle \frac{64}{3}}z^3\left)\right(\mathrm{ln}^2(1z)\zeta (2))+({\displaystyle \frac{64}{3}}z{\displaystyle \frac{128}{3}}z^2`$
$`+{\displaystyle \frac{64}{3}}z^3{\displaystyle \frac{64}{3}}{\displaystyle \frac{1}{1z}})\mathrm{ln}z\mathrm{ln}(1z)+({\displaystyle \frac{8}{3}}8z+20z^2`$
$`{\displaystyle \frac{16}{3}}z^3+{\displaystyle \frac{8}{3}}{\displaystyle \frac{1}{1z}})\mathrm{ln}^2z+({\displaystyle \frac{592}{9}}+{\displaystyle \frac{1456}{9}}z{\displaystyle \frac{1160}{9}}z^2`$
$`+{\displaystyle \frac{304}{3}}z^3)\mathrm{ln}(1z)+({\displaystyle \frac{296}{9}}{\displaystyle \frac{208}{3}}z+{\displaystyle \frac{268}{3}}z^264z^3`$
$`+{\displaystyle \frac{160}{9}}{\displaystyle \frac{1}{1z}})\mathrm{ln}z+{\displaystyle \frac{3142}{27}}{\displaystyle \frac{6188}{27}}z+{\displaystyle \frac{5218}{27}}z^2`$
$`{\displaystyle \frac{3068}{27}}z^3],`$
$`\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}\mathrm{T}_\mathrm{f}}`$ $`=`$ $`[\{16z(1+z)\mathrm{ln}z+{\displaystyle \frac{32}{3}}+8z8z^2{\displaystyle \frac{32}{3}}z^3\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (A.10)
$`+\{32z(1+z)(2\mathrm{Li}_2(1z)+2\mathrm{ln}z\mathrm{ln}(1z)\mathrm{ln}^2z)`$
$`+({\displaystyle \frac{128}{3}}+32z32z^2{\displaystyle \frac{128}{3}}z^3)\mathrm{ln}(1z)+({\displaystyle \frac{64}{3}}`$
$`64z48z^2+{\displaystyle \frac{128}{3}}z^3)\mathrm{ln}z{\displaystyle \frac{256}{9}}{\displaystyle \frac{448}{3}}z+{\displaystyle \frac{352}{3}}z^2`$
$`+{\displaystyle \frac{544}{9}}z^3\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+\left(32+32z+80z^2\right)\mathrm{S}_{1,2}(1z)+\left(32160z112z^2\right)`$
$`\times \mathrm{Li}_3(1z)+({\displaystyle \frac{64}{3}}48z160z^2+{\displaystyle \frac{64}{3}}z^3)\mathrm{Li}_2(1z)`$
$`+z(1+z)(128\mathrm{ln}(1z)\mathrm{Li}_2(1z)32\mathrm{ln}z\mathrm{Li}_2(1z)`$
$`64\zeta (2)\mathrm{ln}z+64\mathrm{ln}z\mathrm{ln}^2(1z)64\mathrm{ln}^2z\mathrm{ln}(1z)`$
$`+{\displaystyle \frac{40}{3}}\mathrm{ln}^3z)+({\displaystyle \frac{128}{3}}+32z32z^2{\displaystyle \frac{128}{3}}z^3\left)\right(\mathrm{ln}^2(1z)`$
$`\zeta (2))+({\displaystyle \frac{128}{3}}128z96z^2+{\displaystyle \frac{256}{3}}z^3)\mathrm{ln}z\mathrm{ln}(1z)`$
$`+({\displaystyle \frac{32}{3}}+52z+28z^232z^3)\mathrm{ln}^2z+({\displaystyle \frac{512}{9}}{\displaystyle \frac{896}{3}}z`$
$`+{\displaystyle \frac{704}{3}}z^2+{\displaystyle \frac{1088}{9}}z^3)\mathrm{ln}(1z)+({\displaystyle \frac{256}{9}}+192z128z^2`$
$`{\displaystyle \frac{1088}{9}}z^3)\mathrm{ln}z{\displaystyle \frac{608}{27}}+{\displaystyle \frac{4144}{9}}z{\displaystyle \frac{3280}{9}}z^2{\displaystyle \frac{1984}{27}}z^3],`$
$`\mathrm{\Delta }_{gq,\mathrm{H}}^{(2)}=C_F^2\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}+C_AC_F\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}}+n_fT_fC_F\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}\mathrm{T}_\mathrm{f}},`$ (A.11)
$`\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`[\{(88z+4z^2)\mathrm{ln}(1z)+(4z2z^2)\mathrm{ln}z+4zz^2\}`$ (A.12)
$`\times \mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$
$`+\{16\mathrm{Li}_2(1z)+(88z+4z^2)(3\mathrm{ln}^2(1z)4\zeta (2))`$
$`+\left(16+32z16z^2\right)\mathrm{ln}z\mathrm{ln}(1z)+\left(8z+4z^2\right)\mathrm{ln}^2z`$
$`+\left(36+64z28z^2\right)\mathrm{ln}(1z)+\left(8z+38z^2+{\displaystyle \frac{16}{3}}z^3\right)`$
$`\times \mathrm{ln}z+{\displaystyle \frac{106}{9}}12z4z^2{\displaystyle \frac{124}{9}}z^3\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+\left(32+48z24z^2\right)\left(\mathrm{S}_{1,2}(1z)\zeta (2)\mathrm{ln}z\right)`$
$`+\left(1616z+8z^2\right)\left(\mathrm{Li}_3(1z)\mathrm{ln}(1z)\mathrm{Li}_2(1z)\right)`$
$`+\left(44z+2z^2\right)\left({\displaystyle \frac{13}{3}}\mathrm{ln}^3(1z)16\zeta (2)\mathrm{ln}(1z)+8\zeta (3)\right)`$
$`+\left(16+8z4z^2\right)\mathrm{ln}z\mathrm{Li}_2(1z)+\left(824z+12z^2\right)`$
$`\times \mathrm{ln}^2z\mathrm{ln}(1z)+\left(24+40z20z^2\right)\mathrm{ln}z\mathrm{ln}^2(1z)`$
$`+\left({\displaystyle \frac{10}{3}}z{\displaystyle \frac{5}{3}}z^2\right)\mathrm{ln}^3z+\left({\displaystyle \frac{32}{3}}+16z+16z^2+{\displaystyle \frac{32}{3}}z^3\right)`$
$`\times (\mathrm{Li}_2(z)+\mathrm{ln}z\mathrm{ln}(1+z))+({\displaystyle \frac{92}{3}}+56z+42z^2`$
$`+{\displaystyle \frac{16}{3}}z^3)\mathrm{Li}_2(1z)+(60+94z45z^2)\mathrm{ln}^2(1z)`$
$`+(3656z+88z^2+{\displaystyle \frac{32}{3}}z^3)\mathrm{ln}z\mathrm{ln}(1z)+(14z`$
$`{\displaystyle \frac{63}{2}}z^2{\displaystyle \frac{40}{3}}z^3)\mathrm{ln}^2z+({\displaystyle \frac{232}{3}}104z+60z^2+{\displaystyle \frac{16}{3}}z^3)\zeta (2)`$
$`+({\displaystyle \frac{878}{9}}144z+48z^2{\displaystyle \frac{248}{9}}z^3)\mathrm{ln}(1z)+({\displaystyle \frac{106}{9}}`$
$`+69z{\displaystyle \frac{214}{3}}z^2+{\displaystyle \frac{112}{9}}z^3)\mathrm{ln}z{\displaystyle \frac{1393}{54}}{\displaystyle \frac{130}{9}}z+{\displaystyle \frac{17}{18}}z^2`$
$`+{\displaystyle \frac{1304}{27}}z^3],`$
$`\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}}`$ $`=`$ $`[\{(2424z+12z^2)\mathrm{ln}(1z)+(2424z24z^2)\mathrm{ln}z`$ (A.13)
$`{\displaystyle \frac{230}{3}}+{\displaystyle \frac{188}{3}}z{\displaystyle \frac{4}{3}}z^2+8z^3\left\}\mathrm{ln}^2\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+\{(48144z72z^2)\mathrm{Li}_2(1z)+(16+16z+8z^2)`$
$`\times \left(\mathrm{Li}_2(z)+\mathrm{ln}z\mathrm{ln}(1+z)\right)+\left(7272z+36z^2\right)`$
$`\times \mathrm{ln}^2(1z)+\left(12864z112z^2\right)\mathrm{ln}z\mathrm{ln}(1z)`$
$`+(24+32z+28z^2)\mathrm{ln}^2z+16z\zeta (2)+({\displaystyle \frac{868}{3}}+{\displaystyle \frac{736}{3}}z`$
$`+{\displaystyle \frac{40}{3}}z^2+32z^3)\mathrm{ln}(1z)+(124{\displaystyle \frac{448}{3}}z+{\displaystyle \frac{8}{3}}z^232z^3)`$
$`\times \mathrm{ln}z+{\displaystyle \frac{2422}{9}}{\displaystyle \frac{1724}{9}}z{\displaystyle \frac{362}{9}}z^2{\displaystyle \frac{32}{9}}z^3\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+(136248z148z^2)\mathrm{S}_{1,2}(1z)+(104+344z`$
$`+148z^2)\mathrm{Li}_3(1z)+(8+8z+4z^2\left)\right(4\mathrm{Li}_3({\displaystyle \frac{1z}{1+z}})`$
$`4\mathrm{Li}_3\left({\displaystyle \frac{1z}{1+z}}\right)+2\mathrm{Li}_3(z)4\mathrm{ln}z\mathrm{Li}_2(z)`$
$`+4\mathrm{ln}(1z)\mathrm{Li}_2(z)+4\mathrm{ln}z\mathrm{ln}(1z)\mathrm{ln}(1+z)`$
$`3\mathrm{ln}^2z\mathrm{ln}(1+z))+(80304z136z^2)\mathrm{ln}(1z)`$
$`\times \mathrm{Li}_2(1z)+(4832z32z^2)\mathrm{ln}z\mathrm{Li}_2(1z)+({\displaystyle \frac{140}{3}}`$
$`{\displaystyle \frac{140}{3}}z+{\displaystyle \frac{70}{3}}z^2)\mathrm{ln}^3(1z)+(812z10z^2)\mathrm{ln}^3z`$
$`+(64+48z+64z^2)\mathrm{ln}^2z\mathrm{ln}(1z)+(13260z`$
$`114z^2)\mathrm{ln}z\mathrm{ln}^2(1z)+32z\zeta (2)\mathrm{ln}(1z)+(64+112z`$
$`+80z^2)\zeta (2)\mathrm{ln}z+(136112z+68z^2)\zeta (3)+({\displaystyle \frac{538}{3}}`$
$`+40z+34z^2{\displaystyle \frac{40}{3}}z^3)\mathrm{Li}_2(1z)+({\displaystyle \frac{88}{3}}40z12z^2`$
$`{\displaystyle \frac{16}{3}}z^3\left)\right(\mathrm{Li}_2(z)+\mathrm{ln}z\mathrm{ln}(1+z))+({\displaystyle \frac{784}{3}}+{\displaystyle \frac{634}{3}}z`$
$`+{\displaystyle \frac{97}{3}}z^2+32z^3)\mathrm{ln}^2(1z)+({\displaystyle \frac{692}{3}}{\displaystyle \frac{832}{3}}z{\displaystyle \frac{40}{3}}z^2`$
$`64z^3)\mathrm{ln}z\mathrm{ln}(1z)+(62+{\displaystyle \frac{344}{3}}z{\displaystyle \frac{10}{3}}z^2+{\displaystyle \frac{68}{3}}z^3)`$
$`\times \mathrm{ln}^2z+({\displaystyle \frac{574}{3}}128z58z^2{\displaystyle \frac{104}{3}}z^3)\zeta (2)+({\displaystyle \frac{4342}{9}}`$
$`{\displaystyle \frac{904}{3}}z116z^2{\displaystyle \frac{64}{9}}z^3)\mathrm{ln}(1z)+({\displaystyle \frac{2402}{9}}+{\displaystyle \frac{2660}{9}}z`$
$`+{\displaystyle \frac{1079}{9}}z^2+{\displaystyle \frac{224}{9}}z^3)\mathrm{ln}z{\displaystyle \frac{21539}{54}}+{\displaystyle \frac{9962}{27}}z+{\displaystyle \frac{1171}{54}}z^2`$
$`{\displaystyle \frac{238}{27}}z^3],`$
$`\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}\mathrm{T}_\mathrm{f}}`$ $`=`$ $`[\{{\displaystyle \frac{16}{3}}{\displaystyle \frac{16}{3}}z+{\displaystyle \frac{8}{3}}z^2\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (A.14)
$`+\left\{({\displaystyle \frac{32}{3}}{\displaystyle \frac{32}{3}}z+{\displaystyle \frac{16}{3}}z^2)(\mathrm{ln}(1z)\mathrm{ln}z)\right){\displaystyle \frac{232}{9}}+{\displaystyle \frac{304}{9}}z`$
$`{\displaystyle \frac{152}{9}}z^2\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+\left({\displaystyle \frac{8}{3}}{\displaystyle \frac{8}{3}}z+{\displaystyle \frac{4}{3}}z^2\right)\left(\mathrm{ln}^2(1z)+2\mathrm{ln}^2z4\mathrm{ln}z\mathrm{ln}(1z)\right)`$
$`+({\displaystyle \frac{104}{3}}+{\displaystyle \frac{128}{3}}z24z^2)\mathrm{ln}(1z)+({\displaystyle \frac{232}{9}}{\displaystyle \frac{304}{9}}z`$
$`+{\displaystyle \frac{152}{9}}z^2)\mathrm{ln}z+{\displaystyle \frac{1060}{27}}{\displaystyle \frac{1672}{27}}z+{\displaystyle \frac{716}{27}}z^2],`$
$`\mathrm{\Delta }_{q_1q_2,\mathrm{H}}^{(2)}=\mathrm{\Delta }_{q_1\overline{q}_2,\mathrm{H}}^{(2)}=C_F^2\mathrm{\Delta }_{q_1q_2,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2},`$ (A.15)
$`\mathrm{\Delta }_{q_1q_2,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`[\{(1616z4z^2)\mathrm{ln}z24+16z+8z^2\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (A.16)
$`+\{(16+16z+4z^2)(4\mathrm{Li}_2(1z)4\mathrm{ln}z\mathrm{ln}(1z)`$
$`+\mathrm{ln}^2z)+(96+64z+32z^2)\mathrm{ln}(1z)+(4832z`$
$`20z^2)\mathrm{ln}z+10272z30z^2\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+(16+16z+4z^2)(8\mathrm{S}_{1,2}(1z)+8\mathrm{Li}_3(1z){\displaystyle \frac{1}{3}}\mathrm{ln}^3z`$
$`8\mathrm{ln}(1z)\mathrm{Li}_2(1z)2\mathrm{ln}z\mathrm{Li}_2(1z)4\mathrm{ln}z\mathrm{ln}^2(1z)`$
$`+2\mathrm{ln}^2z\mathrm{ln}(1z)+4\zeta (2)\mathrm{ln}z)+(48+32z+8z^2)`$
$`\times \mathrm{Li}_2(1z)+\left(96+64z+32z^2\right)\left(\mathrm{ln}^2(1z)\zeta (2)\right)`$
$`+(9664z40z^2)\mathrm{ln}z\mathrm{ln}(1z)+(24+32z`$
$`+8z^2)\mathrm{ln}^2z+(204144z60z^2)\mathrm{ln}(1z)+(118`$
$`+88z+58z^2)\mathrm{ln}z210+188z+22z^2],`$
$`\mathrm{\Delta }_{qq,\mathrm{H}}^{(2)}=C_AC_F^2\mathrm{\Delta }_{qq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}^2}+C_F^3\mathrm{\Delta }_{qq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^3}+C_F^2\mathrm{\Delta }_{qq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2},`$ (A.17)
$`\mathrm{\Delta }_{qq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`[(16+16z8z^2)\mathrm{S}_{1,2}(1z)+(8z+12z^2)\mathrm{ln}^2z`$ (A.18)
$`+(16z36z^2)\mathrm{ln}z1032z+42z^2],`$
$`\mathrm{\Delta }_{qq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^3}`$ $`=`$ $`[(3232z+16z^2)\mathrm{S}_{1,2}(1z)+(16z24z^2)\mathrm{ln}^2z`$ (A.19)
$`+(32z+72z^2)\mathrm{ln}z+20+64z84z^2],`$
$`\mathrm{\Delta }_{qq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`\mathrm{\Delta }_{q_1q_2,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2},`$ (A.20)
$`\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2)}`$ $`=`$ $`C_AC_F^2\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}^2}+C_F^3\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^3}+C_F^2\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ (A.21)
$`+n_fT_fC_F^2\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2\mathrm{T}_\mathrm{f}},`$
$`\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`[\{{\displaystyle \frac{88}{3}}+88z88z^2+{\displaystyle \frac{88}{3}}z^3\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (A.22)
$`+(8+8z+4z^2)(4\mathrm{S}_{1,2}(z)6\mathrm{Li}_3(z)`$
$`4\mathrm{ln}(1+z)\mathrm{Li}_2(z)+6\mathrm{ln}z\mathrm{Li}_2(z)2\mathrm{ln}z\mathrm{ln}^2(1+z)`$
$`+3\mathrm{ln}^2z\mathrm{ln}(1+z)2\zeta (2)\mathrm{ln}(1+z)4\zeta (3))+({\displaystyle \frac{16}{3}}16z`$
$`+16z^2{\displaystyle \frac{16}{3}}z^3\left)\right(5\mathrm{Li}_2(1z)+2\mathrm{ln}^2(1z)8\mathrm{ln}(1z))`$
$`+\left(32z+16z^2+{\displaystyle \frac{32}{3}}z^3\right)\left(\mathrm{Li}_2(z)+\mathrm{ln}z\mathrm{ln}(1+z)\right)`$
$`+\left({\displaystyle \frac{176}{3}}+192z168z^2+64z^3\right)\zeta (2)+\left(40z{\displaystyle \frac{32}{3}}z^3\right)`$
$`\times \mathrm{ln}^2z+\left({\displaystyle \frac{88}{3}}{\displaystyle \frac{296}{3}}z+{\displaystyle \frac{500}{3}}z^2{\displaystyle \frac{128}{3}}z^3\right)\mathrm{ln}z+130316z`$
$`+{\displaystyle \frac{890}{3}}z^2{\displaystyle \frac{332}{3}}z^3],`$
$`\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^3}`$ $`=`$ $`[\{({\displaystyle \frac{64}{3}}64z+64z^2{\displaystyle \frac{64}{3}}z^3)\mathrm{ln}(1z)+(32z32z^2`$ (A.23)
$`+{\displaystyle \frac{64}{3}}z^3)\mathrm{ln}z+{\displaystyle \frac{64}{3}}z{\displaystyle \frac{64}{3}}z^2\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`+(16+16z+8z^2)(4\mathrm{S}_{1,2}(z)+6\mathrm{Li}_3(z)+4\zeta (3)`$
$`+4\mathrm{ln}(1+z)\mathrm{Li}_2(z)6\mathrm{ln}z\mathrm{Li}_2(z)+2\mathrm{ln}z\mathrm{ln}^2(1+z)`$
$`3\mathrm{ln}^2z\mathrm{ln}(1+z)+2\zeta (2)\mathrm{ln}(1+z))+(32+128z`$
$`128z^2+{\displaystyle \frac{160}{3}}z^3)\mathrm{Li}_2(1z)+(64z32z^2\left)\right(\mathrm{Li}_2(z)`$
$`+\mathrm{ln}z\mathrm{ln}(1+z))+({\displaystyle \frac{32}{3}}32z+32z^2{\displaystyle \frac{32}{3}}z^3)`$
$`\times \mathrm{ln}^2(1z)+\left({\displaystyle \frac{64}{3}}+96z96z^2+{\displaystyle \frac{128}{3}}z^3\right)`$
$`\times \mathrm{ln}z\mathrm{ln}(1z)+(40z+48z^2{\displaystyle \frac{16}{3}}z^3)\mathrm{ln}^2z+({\displaystyle \frac{224}{3}}256z`$
$`+208z^2{\displaystyle \frac{224}{3}}z^3)\zeta (2)+(32+{\displaystyle \frac{352}{3}}z{\displaystyle \frac{352}{3}}z^2+32z^3)`$
$`\times \mathrm{ln}(1z)+\left({\displaystyle \frac{256}{3}}z{\displaystyle \frac{200}{3}}z^232z^3\right)\mathrm{ln}z{\displaystyle \frac{164}{3}}64z`$
$`+{\displaystyle \frac{260}{3}}z^2+32z^3],`$
$`\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`\mathrm{\Delta }_{q_1q_2,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2},`$ (A.24)
$`\mathrm{\Delta }_{q\overline{q},\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2\mathrm{T}_\mathrm{f}}`$ $`=`$ $`[\{{\displaystyle \frac{32}{3}}32z+32z^2{\displaystyle \frac{32}{3}}z^3\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (A.25)
$`+({\displaystyle \frac{64}{9}}{\displaystyle \frac{64}{3}}z+{\displaystyle \frac{64}{3}}z^2{\displaystyle \frac{64}{9}}z^3)\mathrm{ln}(1z)+({\displaystyle \frac{32}{3}}+{\displaystyle \frac{160}{3}}z`$
$`{\displaystyle \frac{128}{3}}z^2+{\displaystyle \frac{128}{9}}z^3)\mathrm{ln}z{\displaystyle \frac{368}{27}}+{\displaystyle \frac{592}{9}}z{\displaystyle \frac{688}{9}}z^2`$
$`+{\displaystyle \frac{656}{27}}z^3].`$
As discussed below Eq. (4.10) the corrections are dominated by the terms $`c_{i,j}\mathrm{ln}^j(1z)`$ in the coefficient functions. They are obtained by taking the limit $`z1`$ in the expressions above. Most of them vanish. Those which survive are given by

$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gg,\mathrm{H}}^{(1),\mathrm{C}_\mathrm{A}}`$ $`=`$ $`\left[16\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)32\mathrm{ln}(1z)+8\right],`$ (A.26)
$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}^2}`$ $`=`$ $`[\{128\mathrm{ln}(1z)+{\displaystyle \frac{184}{3}}\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+\{384\mathrm{ln}^2(1z)`$ (A.27)
$`+{\displaystyle \frac{1024}{3}}\mathrm{ln}(1z)+160\zeta (2){\displaystyle \frac{1336}{9}}\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})`$
$`256\mathrm{ln}^3(1z)+{\displaystyle \frac{1096}{3}}\mathrm{ln}^2(1z)+(320\zeta (2)`$
$`{\displaystyle \frac{2660}{9}})\mathrm{ln}(1z)624\zeta (3){\displaystyle \frac{784}{3}}\zeta (2)+{\displaystyle \frac{4228}{27}}],`$
$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gg,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{T}_\mathrm{f}}`$ $`=`$ $`[{\displaystyle \frac{32}{3}}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+\{{\displaystyle \frac{128}{3}}\mathrm{ln}(1z)+{\displaystyle \frac{416}{9}}\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$
$`{\displaystyle \frac{128}{3}}\mathrm{ln}^2(1z)+{\displaystyle \frac{808}{9}}\mathrm{ln}(1z)+{\displaystyle \frac{128}{3}}\zeta (2){\displaystyle \frac{1232}{27}}],`$
$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gq,\mathrm{H}}^{(1),\mathrm{C}_\mathrm{F}}`$ $`=`$ $`\left[2\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+4\mathrm{ln}(1z)+2\right],`$ (A.29)
$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`[\{4\mathrm{ln}(1z)+3\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+\{12\mathrm{ln}^2(1z)`$
$`16\zeta (2)18\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})+{\displaystyle \frac{26}{3}}\mathrm{ln}^3(1z)11\mathrm{ln}^2(1z)`$
$`+(32\zeta (2)26)\mathrm{ln}(1z)+16\zeta (3)+12\zeta (2)+9],`$
$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}}`$ $`=`$ $`[\{12\mathrm{ln}(1z){\displaystyle \frac{22}{3}}\}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+\{36\mathrm{ln}^2(1z)`$ (A.31)
$`+{\displaystyle \frac{4}{3}}\mathrm{ln}(1z)4\zeta (2)+{\displaystyle \frac{304}{9}}\left\}\mathrm{ln}\right({\displaystyle \frac{m^2}{\mu ^2}})+{\displaystyle \frac{70}{3}}\mathrm{ln}^3(1z)`$
$`+{\displaystyle \frac{43}{3}}\mathrm{ln}^2(1z)+\left(8\zeta (2)+58\right)\mathrm{ln}(1z)+62\zeta (3)`$
$`+14\zeta (2){\displaystyle \frac{460}{27}}],`$
$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gq,\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}\mathrm{T}_\mathrm{f}}`$ $`=`$ $`[{\displaystyle \frac{8}{3}}\mathrm{ln}^2\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+\{{\displaystyle \frac{16}{3}}\mathrm{ln}(1z){\displaystyle \frac{80}{9}}\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (A.32)
$`+{\displaystyle \frac{4}{3}}\mathrm{ln}^2(1z)16\mathrm{ln}(1z)+{\displaystyle \frac{104}{27}}].`$
The nonvanishing coefficients $`c_{i,j}`$ in Eq. (4.9) can be read off from the equations above.

## Appendix B

In this Appendix we present the coefficient functions for pseudo-scalar Higgs boson production. To shorten the expressions we give the difference between the coefficient functions originating from pseudo-scalar and scalar production only, namely,

$`\mathrm{\Delta }_{gg,\mathrm{A}\mathrm{H}}^{(1),\mathrm{S}+\mathrm{V}}`$ $`=`$ $`C_A\left[8\delta (1z)\right],`$ (B.1)
$`\mathrm{\Delta }_{gg,\mathrm{A}\mathrm{H}}^{(1)}`$ $`=`$ $`0,`$ (B.2)
$`\mathrm{\Delta }_{gq,\mathrm{A}\mathrm{H}}^{(1)}`$ $`=`$ $`0,`$ (B.3)
$`\mathrm{\Delta }_{q\overline{q},\mathrm{A}\mathrm{H}}^{(1)}`$ $`=`$ $`0.`$ (B.4)
In the NNLO corrections below we indicate by $`O_{12}`$ the contribution coming from the interferences between the operators $`O_1`$ and $`O_2`$ in Eqs. (3.9), Eqs. (3.10). In the coefficient functions one has to put $`O_{12}=1`$,

$`\mathrm{\Delta }_{gg,\mathrm{A}\mathrm{H}}^{(2),\mathrm{S}+\mathrm{V}}`$ $`=`$ $`C_A^2[\{64𝒟_0(z){\displaystyle \frac{20}{3}}\delta (1z)\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (B.5)
$`+128𝒟_1(z)+\delta (1z)({\displaystyle \frac{215}{3}}+64\zeta (2))]`$
$`+C_AT_fn_f[\{{\displaystyle \frac{8}{3}}\delta (1z)\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$
$`+\delta (1z)({\displaystyle \frac{196}{9}}+O_{12}({\displaystyle \frac{32}{3}}\mathrm{ln}\left({\displaystyle \frac{\mu _r^2}{m_t^2}}\right){\displaystyle \frac{16}{3}})]`$
$`+C_FT_fn_f\left[\left\{8\delta (1z)\right\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)4\delta (1z)\right].`$
For the colour decomposition of the $`n_f`$-part of the expression above see the remark below Eq. (A.6). Also

$`\mathrm{\Delta }_{gg,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}^2}`$ $`=`$ $`[\{128z+64z^264z^3\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (B.6)
$`16z\mathrm{ln}^2z+(256z+128z^2128z^3)\mathrm{ln}(1z)+(32`$
$`+{\displaystyle \frac{440}{3}}z64z^2+64z^3{\displaystyle \frac{64}{1z}})\mathrm{ln}z+{\displaystyle \frac{44}{3}}+{\displaystyle \frac{232}{3}}z{\displaystyle \frac{452}{3}}z^2`$
$`+{\displaystyle \frac{176}{3}}z^3],`$
$`\mathrm{\Delta }_{gg,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{T}_\mathrm{f}}`$ $`=`$ $`\left[{\displaystyle \frac{32}{3}}z\mathrm{ln}z+{\displaystyle \frac{16}{3}}+{\displaystyle \frac{16}{3}}z{\displaystyle \frac{32}{3}}z^2\right],`$ (B.7)
$`\mathrm{\Delta }_{gg,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}\mathrm{T}_\mathrm{f}}`$ $`=`$ $`\left[16z\mathrm{ln}^2z16+32z16z^2\right],`$ (B.8)
$`\mathrm{\Delta }_{gq,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`[8z\mathrm{ln}^2z24z\mathrm{ln}z4+16z12z^2+O_{12}(12`$ (B.9)
$`+12z^2)],`$
$`\mathrm{\Delta }_{gq,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}}`$ $`=`$ $`[\{3232z+16z^2\}\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)`$ (B.10)
$`16z\mathrm{ln}^2z+(6464z+32z^2)\mathrm{ln}(1z)+(80z`$
$`16z^2)\mathrm{ln}z+6048z+4z^2],`$
$`\mathrm{\Delta }_{gq,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}\mathrm{T}_\mathrm{f}}`$ $`=`$ $`0,`$ (B.11)
$`\mathrm{\Delta }_{q_1q_2,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`\left[16z\mathrm{ln}^2z+\left(32+48z\right)\mathrm{ln}z+8896z+8z^2\right],`$ (B.12)
$`\mathrm{\Delta }_{qq,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`\left[16z\mathrm{ln}^2z+32z\mathrm{ln}z+3232z\right],`$ (B.13)
$`\mathrm{\Delta }_{qq,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^3}`$ $`=`$ $`\left[32z\mathrm{ln}^2z64z\mathrm{ln}z64+64z\right],`$ (B.14)
$`\mathrm{\Delta }_{qq,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`\mathrm{\Delta }_{q_1q_2,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2},`$ (B.15)
$`\mathrm{\Delta }_{q_1\overline{q}_2,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`\mathrm{\Delta }_{q_1q_2,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2},`$ (B.16)
$`\mathrm{\Delta }_{q\overline{q},\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`\left[16z\mathrm{ln}^2z+{\displaystyle \frac{80}{3}}z\mathrm{ln}z+{\displaystyle \frac{8}{3}}+{\displaystyle \frac{56}{3}}z^2{\displaystyle \frac{64}{3}}z^3\right],`$ (B.17)
$`\mathrm{\Delta }_{q\overline{q},\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^3}`$ $`=`$ $`[16z\mathrm{ln}^2z32z\mathrm{ln}z+64160z+96z^2`$ (B.18)
$`+O_{12}(48+96z48z^2)],`$
$`\mathrm{\Delta }_{q\overline{q},\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2}`$ $`=`$ $`\mathrm{\Delta }_{q_1q_2,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2},`$ (B.19)
$`\mathrm{\Delta }_{q\overline{q},\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{F}^2\mathrm{T}_\mathrm{f}}`$ $`=`$ $`\left[{\displaystyle \frac{64}{3}}z\mathrm{ln}z{\displaystyle \frac{32}{3}}+{\displaystyle \frac{32}{3}}z^2\right].`$ (B.20)
The nonvanishing coefficients $`c_{i,j}`$ in Eq. (4.9) appearing in the difference between the pseudo-scalar and scalar coefficient functions are determined from the expressions above by taking the limit $`z1`$. Those which survive in this limit are given by

$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gg,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}^2}`$ $`=`$ $`\left[128\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)256\mathrm{ln}(1z)+64\right],`$ (B.21)
$`\underset{z1}{\text{lim}}\mathrm{\Delta }_{gg,\mathrm{A}\mathrm{H}}^{(2),\mathrm{C}_\mathrm{A}\mathrm{C}_\mathrm{F}}`$ $`=`$ $`\left[16\mathrm{ln}\left({\displaystyle \frac{m^2}{\mu ^2}}\right)+32\mathrm{ln}(1z)+16\right].`$ (B.22)
Figure Captions

The total cross section $`\sigma _{\mathrm{tot}}`$ plotted as a function of the Higgs mass at $`\sqrt{S}=14\mathrm{TeV}`$ with $`\mu =m`$. The NLO plots are presented for the total (solid line) and the subprocesses $`gg`$ (long-dashed line), $`10\times \mathrm{abs}(gq+g\overline{q})`$ (dot-dashed line) and $`100\times (q\overline{q})`$ (short-dashed line). Also shown is the LO result (dotted line).

The total cross section $`\sigma _{\mathrm{tot}}`$ plotted as a function of the Higgs mass at $`\sqrt{S}=14\mathrm{TeV}`$ with $`\mu =m`$. The NNLO plots are presented for the total (solid line) and the subprocesses $`gg`$ (long-dashed line), $`\mathrm{abs}(gq+g\overline{q})`$ (dot-dashed line), $`100\times (q\overline{q})`$ (dotted line) and $`100\times (qq+\overline{q}\overline{q})`$ (short-dashed line).

The total cross section $`\sigma _{\mathrm{tot}}`$ with all channels included plotted as a function of the Higgs mass at $`\sqrt{S}=14\mathrm{TeV}`$ with $`\mu =m`$; NNLO (solid line), NLO (dashed line) and LO (dotted line).

The quantity $`N(\mu /\mu _0)`$ (see Eq. (4) at $`\sqrt{S}=14\mathrm{TeV}`$ plotted in the range $`0.1<\mu /\mu _0<10`$ with $`\mu _0=m`$ and $`m=100\mathrm{GeV}/\mathrm{c}^2`$. The results are shown for LO (dotted line), NLO (dashed line) and NNLO (solid line).

Same as Fig. 4 but now for $`m=200\mathrm{GeV}/\mathrm{c}^2`$.

Same as Fig. 4 but now for $`m=300\mathrm{GeV}/\mathrm{c}^2`$.

The $`K`$-factors in NLO and NNLO at $`\sqrt{S}=14\mathrm{TeV}`$ as a function of the Higgs mass using the MRST-sets; $`K^{\mathrm{NNLO}}`$ (solid line), $`K^{\mathrm{NLO}}`$ (dot-dashed line).

The $`K`$-factors in NLO at $`\sqrt{S}=14\mathrm{TeV}`$ as a function of the Higgs mass using the following parton density sets MRST01 (solid line), GRV98 (dashed line) and CTEQ6 (dotted line).

The ratios $`R=\sigma _{\mathrm{tot}}/\sigma _{\mathrm{tot}}^{\mathrm{MRST}}`$ (Eq.(4.5)) in NLO at $`\sqrt{S}=14\mathrm{TeV}`$ and $`\mu =m`$ as a function of the Higgs mass; $`R^{\mathrm{GRV}}`$ (solid line), $`R^{\mathrm{CTEQ}}`$ (dotted line).

The ratio $`R(x_{\mathrm{max}})`$ (see Eq. (4.7)) for proton-proton collisions (LHC) where $`x_{\mathrm{max}}=5x`$ with $`x=m^2/S`$. The CM energy and scale are given by $`\sqrt{S}=14\mathrm{TeV}`$ and $`\mu =m`$ respectively; NNLO (solid line), NLO (dashed line).

Same as in Fig. 10a but for proton-anti-proton collisions (TEVATRON) at $`\sqrt{S}=2\mathrm{TeV}`$ and $`\mu =m`$.

The ratios in Eq. (4.12) in NLO for proton-proton collisions (LHC) as a function of the Higgs mass at $`\sqrt{S}=14\mathrm{TeV}`$ and $`\mu =m`$; $`R^{\mathrm{SV}}`$(Eq. (4.1)) (dot-dashed line), $`R^{\mathrm{SVL}}`$(Eq. (4.1)) (dashed line), $`R^{\mathrm{SV}}`$(Eq. (4)) (dotted line), $`R^{\mathrm{SVL}}`$(Eq. (4)) (solid line).

The ratios in Eq. (4.12) in NNLO for proton-proton collisions (LHC) as a function of the Higgs mass at $`\sqrt{S}=14\mathrm{TeV}`$ and $`\mu =m`$; $`R^{\mathrm{SV}}`$(Eq. (4.1)) (dot-dashed line), $`R^{\mathrm{SVL}}`$(Eq. (4.1)) (dashed line), $`R^{\mathrm{SV}}`$(Eq. (4)) (dotted line), $`R^{\mathrm{SVL}}`$(Eq. (4)) (solid line).

Same as in Fig. 11a but for proton-anti-proton collisions (TEVATRON) at $`\sqrt{S}=2\mathrm{TeV}`$ and $`\mu =m`$.

Same as in Fig. 11b but now for proton-anti-proton collisions (TEVATRON) at $`\sqrt{S}=2\mathrm{TeV}`$ and $`\mu =m`$.