format: text

text:
# 1 Introduction

## 1 Introduction

In our previous paper , referred to as I in the following, we presented the two–loop calculation of the bosonic contributions to the $`W`$– and $`Z`$–self–energies on the mass shell and discussed the main properties and possible applications. In the present paper II we extend this calculation to a full Standard Model (SM) calculation by including the missing fermion–loop and mixed contributions. Throughout the paper, we shall adopt the terminology calling “bosonic corrections” the one’s represented by diagrams without any fermions and “fermionic corrections” the remaining one’s given by diagrams exhibiting at least one fermion loop.

In I, in particular, we have shown by an explicit calculation that up to two-loop order the purely bosonic corrections to the position of the pole of the gauge-boson propagators are real for arbitrary values of the Higgs-boson mass. In contrast, since the $`W`$– and $`Z`$–bosons decay, at leading order, into light fermion pairs, the light (massless) fermion corrections give rise to a non–zero imaginary part at the one loop level already. The problem of gauge (in)dependence of the complex pole has been extensively discussed in literature . Only recently, two important results have been proven to all orders in perturbation theory: i) the position of the complex pole is a gauge independent quantity ; ii) the branching ratios and partial widths associated with the pole residues are gauge independent . Moreover, it has been shown that the pole mass of the W-boson is an infrared finite quantity with respect to massless photonic corrections. An alternative proof of the infrared finiteness of the two-loop bosonic contributions to the pole of the gauge bosons was presented in I by explicite calculation. It is based on the fact that within dimensional regularization , in $`d=42\epsilon `$ space-time dimensions, the singular $`1/\epsilon `$ terms, which regularize both ultraviolet (UV) and infrared (IR) singularities, are absent after UV renormalization of the position of the pole in the propagators.

In the present paper, besides from completing our previous calculation by including the missing fermion contributions, we will discuss in some detail general features and technical problems which are specifically related to these contributions.

The paper is organized as follows. In Section 2 we briefly reconsider the definition of the pole-mass of the massive gauge-bosons within the SM and remind the reader of some notation given in I. The required analytical results for the massless fermion two-loop master-integrals are presented in Section 3. In Section 4 we discuss the UV renormalization of the pole mass and the interrelation between our results and the one’s familiar from the standard renormalization group approach. In particular, we performed several cross-checks of the singular $`1/\epsilon ^2`$\- and $`1/\epsilon `$-terms. General aspects as well as numerical results for the finite parts are discussed in Section 5. Some technical details and a number of our analytical results will be presented in Appendices. In Appendix A we present a set of non–standard binomial sums which are needed for the $`\epsilon `$–expansion of some of the hypergeometric functions entering the master–integrals. The one–loop fermion contribution to the pole masses of the gauge–bosons are reproduced in Appendix B. The reducible two–loop corrections obtained by mass renormalization of the one–loop fermion contributions may be found in Appendix C. In Appendix D the bare two–loop on–shell self–energy contribution of the massless fermions are given in exact analytical form. The corresponding contributions involving the top quark are presented in terms of the first few coefficients of the expansion in $`\mathrm{sin}^2\theta _W`$ and specific mass ratios in Appendix E.

## 2 Pole mass: definition and calculation

The position of the pole $`s_P`$ of the propagator of a massive gauge-boson in a quantum field theory is a solution for $`p^2`$ at which the inverse of the connected full propagator equals zero, i.e.,

$$s_Pm^2\mathrm{\Pi }(s_P,m^2,\mathrm{})=0,$$
(2.1)
where $`\mathrm{\Pi }(p^2,\mathrm{})`$ is the transversal part of the one-particle irreducible self-energy. The latter depends on all SM parameters but, in order to the keep notation simple, we have indicated explicitly only the dependence on the external momentum $`p`$ and in some cases also $`m`$, where $`m`$ is the mass of the particle under consideration. This can be either the bare mass or the renormalized mass defined in some particular renormalization scheme.

Generally, the pole $`s_P`$ is located in the complex plane of $`p^2`$ and has a real and an imaginary part. By writing

$$s_PM^2\mathrm{i}M\mathrm{\Gamma },$$
(2.2)
the real part defines $`M`$ which we call the pole mass while the imaginary part is related to the width $`\mathrm{\Gamma }`$ of the particle. This is the natural generalization of the physical mass of a stable particle, which is defined by the mass of its asymptotic scattering state .

For the remainder of the paper we will adopt the following notation: capital $`M`$ always denotes the pole mass; lower case $`m`$ stands for the renormalized mass in the $`\overline{\mathrm{MS}}`$ scheme, while $`m_0`$ denotes the bare mass. In addition we use $`e`$, $`g`$ and $`g_s`$ to denote the $`U(1)_{\mathrm{em}}`$, $`SU(2)_\mathrm{L}`$ and $`SU(3)_\mathrm{c}`$ couplings of the SM in the $`\overline{\mathrm{MS}}`$ scheme.

In perturbation theory (2.1) is to be solved order by order. To two loops we have the solution666Similarly, it is easy to find the on-shell wave-function renormalization constant $`Z_2`$. Up to two loops it is given by the following equation :
$$Z_2^1=1\mathrm{\Pi }^{(1)}{}_{}{}^{}(m^2,m^2,\mathrm{})\mathrm{\Pi }^{(2)}{}_{}{}^{}(m^2,m^2,\mathrm{})\mathrm{\Pi }^{(1)}(m^2,m^2,\mathrm{})\mathrm{\Pi }^{(1)}{}_{}{}^{\prime \prime }(m^2,m^2,\mathrm{}),$$
where $`m^2`$ is a bare or renormalized mass. For a recent discussion see also .

$`s_P`$ $`=`$ $`m^2+\mathrm{\Pi }^{(1)}(m^2,m^2,\mathrm{})+\mathrm{\Pi }^{(2)}(m^2,m^2,\mathrm{})+\mathrm{\Pi }^{(1)}(m^2,m^2,\mathrm{})\mathrm{\Pi }^{(1)}{}_{}{}^{}(m^2,m^2,\mathrm{}),`$ (2.3)
which yields the pole mass $`M^2`$ and the width $`\mathrm{\Gamma }`$ at this order. $`\mathrm{\Pi }^{(L)}`$ is the bare ($`m=m_0`$) or $`\overline{\mathrm{MS}}`$ -renormalized ($`m`$ the $`\overline{\mathrm{MS}}`$ -mass) $`L`$-loop contribution to $`\mathrm{\Pi }`$, and the prime denotes the derivative with respect to $`p^2`$.

In this paper we show by explicit calculation at the two-loop level that the fermion contribution (including mixed terms) to the propagator pole $`s_P`$ of a gauge boson is a gauge invariant and infrared stable quantity. For more details concerning the tensor decomposition of propagators and abbreviations adopted we refer to Sec. 2 of I.

In order to find the relationship between the poles of the gauge boson-propagators and the $`\overline{\mathrm{MS}}`$ masses $`m_Z^2,m_W^2`$ we have to compute the one- and two-loop self-energies for the $`Z`$\- and $`W`$-bosons at $`p^2=m_Z^2`$ and $`p^2=m_W^2`$, respectively777There exist a number of programs for analytical and/or numerical calculations of two–loop self-energies. A selection may be found in ..

For the calculation of the two-loop fermion contributions we again use the strategy described in our previous paper I. Let us present here some basic features. In order to be able to work with manifestly gauge independent parameter renormalization constants we have to include the Higgs tadpole diagrams. To keep control of gauge invariance we work in the $`R_\xi `$ gauge with three different gauge parameters $`\xi _W,\xi _Z`$ and $`\xi _\gamma `$. However, in order to avoid additional mass parameters like $`\sqrt{\xi _W}m_W`$ and $`\sqrt{\xi _Z}m_Z`$ (ghost masses), we expand the original propagators in an asymptotic series at $`\xi _i=1`$. For the purpose of checking the gauge invariance of our results it turned out to be sufficient to keep the first two terms of the expansion. As an additional cross–check we actually kept one more term.

For calculating diagrams with massless fermion loops we utilize Tarasov’s recurrence relations  (a detailed discussion we postpone to Sec. 3). For diagrams with top-quark and/or Higgs-boson propagators we apply asymptotic expansions with respect to the heavy masses. In these cases we firstly expand the propagator in the weak mixing parameter $`\mathrm{sin}^2\theta _W=1m_W^2/m_Z^20.25`$ and get rid in this way of $`m_W`$ (or $`m_Z`$). All diagrams with top-quarks and/or Higgs-bosons are divided into several prototypes which are presented in Fig. 1. The large mass expansion has been performed with the help of the packages TLAMM  or by the program described in . In the heavy top limit, the leading and next–to–leading contributions to the gauge boson self–energies have been calculated in  and , respectively.

As usual, in order to preserve gauge invariance of the regularized theory, we utilize dimensional regularization  for our calculation. In order to avoid spurious anomalies we adopt an anti-commuting $`\gamma _5`$ in dimensional regularization  and impose vector current conservation by hand in triangle fermion loops which contribute to the two–loop two–point functions . The axial vector current anomalies  cancel for each fermion family by the lepton–quark duality . In contrast to calculations which involve vertex- and/or box-type diagrams , in our calculation of self-energies we do not encounter problems with violation of Ward identities (see also ).

## 3 Master integrals for the massless fermion contributions

This part is devoted to the calculation of the two-loop massless fermion corrections to the pole masses of the gauge bosons. For this class of corrections it is possible to work out the exact analytical results without expansion888In contrast to the bosonic corrections where all diagrams can be expanded from the very beginning in the small parameter $`\mathrm{sin}^2\theta _W`$ the individual diagrams with massless fermion-loops develop threshold singularities which behave like powers of $`\mathrm{ln}\mathrm{sin}^2\theta _W`$. To control these terms we need the exact analytical result..

In contrast to the previous calculations performed in , where results are expressed in terms of the non–minimal set of scalar two–loop integrals, we use here Tarasov’s recurrence relations  which allow us to reduce the number of integrals to a minimal set of master–integrals in our results. The analytical results for these diagrams are presented in . Here we present an independent analytical calculation of the relevant master integrals shown in Fig. 2. Besides the known one’s, we consider here two new master-integrals, shown in Fig. 3, which contribute to the two-loop on-shell top-quark propagator in the limit of massless gauge-bosons. They may also be considered as the “naive” parts of the asymptotic expansion of the original diagrams with massive $`W`$\- and $`Z`$-bosons in the limit, when the external momentum approaches the heavy mass-shell .

For the analytical calculation of these integrals up to the finite parts in $`\epsilon `$ we use the method developed in  (see also ). This approach is based on the possibility to retrieve the analytical results in terms of sums, as predicted by the differential equation method , from several of the first coefficients of the small momentum expansion . Whenever it is possible, we present the exact analytical results in terms of hypergeometric functions. Thereby, we apply Mellin–Barnes techniques developed in .

Let us introduce now the following notation for the finite sums which show up in the master integrals:

$$S_a(n)=\underset{j=1}{\overset{n}{}}\frac{1}{j^a},V_a(n)=\underset{j=1}{\overset{n}{}}\frac{1}{\left(\genfrac{}{}{0pt}{}{2j}{j}\right)}\frac{1}{j^a},W_a(n)=\underset{j=1}{\overset{n}{}}\left(\genfrac{}{}{0pt}{}{2j}{j}\right)\frac{1}{j^a}.$$

Throughout the rest of this paper, the symbol $`S_a`$ will always stand for $`S_a(n1)`$ and $`\overline{S}_a`$ will stand for $`S_a(2n1)`$. The same notation applies for all the other type of sums, e.g., $`V_a=V_a(n1)`$. We denote all master integrals as $`T_{AB\mathrm{}}`$, where the first letter $`T=F,V,J`$ indicates the topology in accordance with notation introduced in ; indices $`A,B,\mathrm{}=0,1,2`$ characterize the relation of the corresponding internal mass to the external momentum: $`0`$ indicates a massless line, $`1`$ corresponds to “internal mass equal to external momentum” and $`2`$ means that mass and momentum are different (see Fig. 2 for details). In our normalization each loop is divided by $`(\pi )^{2\epsilon }\mathrm{\Gamma }(1+\epsilon ).`$

* $`𝐅_{\mathrm{𝟎𝟎𝟎𝟎𝟐}}`$

For this integral we use the representation999Another representation is given by Eq. (104) in .  given in

$`p^2F_{00002}`$ $`=`$ $`2\zeta _2\mathrm{ln}(1+z)+2\mathrm{ln}(z)\mathrm{Li}_2(z)+\mathrm{ln}^2(z)\mathrm{ln}(1+z)+4\mathrm{ln}(1+z)\mathrm{Li}_2(z)`$ (3.4)
$`2\mathrm{L}\mathrm{i}_3(z)2\mathrm{L}\mathrm{i}_3(z)+2S_{1,2}(z^2)4S_{1,2}(z)4S_{1,2}(z),`$
where $`z=\frac{p^2}{M^2}`$ and $`S_{n,p}`$ are Nielsen polylogarithms . For the on–shell gauge-boson propagators two points in the variable $`z`$ are interesting: (i) $`z=\mathrm{cos}^2\theta _W[p^2=m_W^2]`$; (ii)$`z=1/\mathrm{cos}^2\theta _W[p^2=m_Z^2]`$.

* $`𝐅_{\mathrm{𝟐𝟎𝟐𝟎𝟎}}`$

For this integral the following representation101010The general case with two different internal masses is given by Eq. (106) in  . is valid

$`p^2F_{20200}`$ $`=`$ $`\text{Li}_3\left(z\right)6\zeta _3\zeta _2\mathrm{ln}y{\displaystyle \frac{1}{6}}\mathrm{ln}^3y4\mathrm{ln}y\text{Li}_2\left(y\right)`$ (3.5)
$`+4\text{Li}_3\left(y\right)3\text{Li}_3\left(y\right)+{\displaystyle \frac{1}{3}}\text{Li}_3\left(y^3\right),`$
where we have introduced the new “conformal” variable

$$y=\frac{1\sqrt{\frac{z}{z4}}}{1+\sqrt{\frac{z}{z4}}},z=\frac{(1y)^2}{y}\frac{p^2}{M^2}.$$
(3.6)
Here we are interested in the case $`z=1/\mathrm{cos}^2\theta _W`$ only.

* $`𝐕_{\mathrm{𝟐𝟎𝟎𝟐}}`$

The off-shell integral $`V_{X00Y}`$ with two different masses and unit powers of propagators was considered in . The result, Eq. (46) of  was presented as a sum of $`F_2`$ and $`F_4`$ hypergeometric functions of two variables. It turns out, however, that the extraction of the finite part from this result is difficult. The analytical result up to finite parts for this diagram (see Fig. 2) is given by111111The general case with two different internal masses is given by Eq. (95) in .

$`(M^2)^{2\epsilon }V_{2002}={\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{1}{2\epsilon }}{\displaystyle \frac{1}{2}}\zeta _2+{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}z^n\{{\displaystyle \frac{1}{\epsilon }}[{\displaystyle \frac{1}{n}}{\displaystyle \frac{1}{\left(\genfrac{}{}{0pt}{}{2n}{n}\right)}}+4{\displaystyle \frac{1}{n+1}}{\displaystyle \frac{1}{\left(\genfrac{}{}{0pt}{}{2n+2}{n+1}\right)}}]`$
$`{\displaystyle \frac{1}{n}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{1}{n^2}}{\displaystyle \frac{1}{n}}{\displaystyle \frac{1}{\left(\genfrac{}{}{0pt}{}{2n}{n}\right)}}[{\displaystyle \frac{3}{2}}W_1+S_1+4]+{\displaystyle \frac{1}{n+1}}+{\displaystyle \frac{1}{n+1}}{\displaystyle \frac{1}{\left(\genfrac{}{}{0pt}{}{2n+2}{n+1}\right)}}[6W_1(n)+4S_1(n)+16]\},`$
(3.7)
where $`z=p^2/M^2`$ with $`M`$ an internal mass. This result can be written also as

$`\left({\displaystyle \frac{p^2}{z}}\right)^{2\epsilon }V_{2002}={\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{1}{\epsilon }}\left[{\displaystyle \frac{5}{2}}+{\displaystyle \frac{1+y}{1y}}\mathrm{ln}y\right]{\displaystyle \frac{19}{2}}\zeta _2{\displaystyle \frac{1}{2}}\mathrm{Li}_2(z){\displaystyle \frac{1z}{z}}\mathrm{ln}(1z)`$
$`+{\displaystyle \frac{1+y}{1y}}\left[{\displaystyle \frac{7}{2}}\mathrm{Li}_2(y){\displaystyle \frac{1}{2}}\mathrm{Li}_2(y^3){\displaystyle \frac{5}{4}}\mathrm{ln}^2y+2\mathrm{ln}(1+y)\mathrm{ln}y4\mathrm{ln}y+{\displaystyle \frac{3}{2}}\zeta _2\right],`$ (3.8)
where $`y`$ is defined in (3.6). If we compare this with the cases considered previously in , as a novelty, in Eq. (3.7) we encounter series with shifted index of summation ($`n+1`$ instead of $`n`$). The latter give rise to the appearance of terms of different weight in the finite part. They are connected with the presence of UV-divergent subgraphs. For this diagram only the point $`z=1/\mathrm{cos}^2\theta _W`$ $`[p^2=m_Z^2]`$ is of interest here.

* $`𝐉_{\mathrm{𝟎𝟎𝟐}}`$

For the off-shell case, we have two master–integrals: one with all indices equal to one, and a second one with indices 1,1,2. Their $`\epsilon `$-expansion up to order $`O(\epsilon ^2)`$ are given by Eqs. (A1) and (A2) in Appendix A of . For our task only the following finite parts are needed:

$`J_{002}(1,1,1,p^2,m^2)`$ $`=`$ $`(m^2)^{12\epsilon }\left[{\displaystyle \frac{1}{2\epsilon ^2}}+{\displaystyle \frac{z6}{4\epsilon }}\zeta _23+{\displaystyle \frac{13}{8}}z+{\displaystyle \frac{1z^2}{2z}}\mathrm{ln}(1z)\text{Li}_2\left(z\right)\right],`$
$`J_{002}(1,1,2,p^2,m^2)`$ $`=`$ $`(m^2)^{2\epsilon }\left[{\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{1}{2\epsilon }}+{\displaystyle \frac{1}{2}}\zeta _2+{\displaystyle \frac{1z}{z}}\mathrm{ln}(1z)\text{Li}_2\left(z\right)\right].`$
For this diagram again two expansion points are important: $`z=\mathrm{cos}^2\theta _W`$ and $`z=1/\mathrm{cos}^2\theta _W`$.

* $`𝐕_{\mathrm{𝟏𝟎𝟎𝟐}}`$

Let us present some useful details for this integral. Its Mellin-Barnes representation is

$`V_{1002}(\alpha ,\sigma _1,\sigma _2,\beta ,p^2,m^2,M^2)|_{p^2=m^2}=`$
$`{\displaystyle \frac{d^dk_1}{\mathrm{\Gamma }(1+\epsilon )}\frac{d^dk_2}{\mathrm{\Gamma }(1+\epsilon )}P^{(\alpha )}(k_2p,m)P^{(\sigma _1)}(k_1k_2,0)P^{(\sigma _2)}(k_1,0)P^{(\beta )}(k_2,M)}=`$
$`{\displaystyle \frac{(\text{i}^{1d})^2(m^2)^{d\alpha \beta \sigma _1\sigma _2}\mathrm{\Gamma }(\sigma _1+\sigma _2\frac{d}{2})\mathrm{\Gamma }(\frac{d}{2}\sigma _1)\mathrm{\Gamma }(\frac{d}{2}\sigma _2)}{\mathrm{\Gamma }^2(1+\epsilon )\mathrm{\Gamma }(\sigma _1)\mathrm{\Gamma }(\sigma _2)\mathrm{\Gamma }(\alpha )\mathrm{\Gamma }(\beta )\mathrm{\Gamma }(d\sigma _1\sigma _2)}}{\displaystyle \frac{1}{2\pi i}}{\displaystyle _i\mathrm{}^i\mathrm{}}𝑑s\left({\displaystyle \frac{M^2}{m^2}}\right)^s`$
$`{\displaystyle \frac{\mathrm{\Gamma }(s)\mathrm{\Gamma }(\beta +s)\mathrm{\Gamma }(\alpha +\beta +\sigma _1+\sigma _2+sd)\mathrm{\Gamma }(2d2\sigma _12\sigma _22\beta \alpha 2s)}{\mathrm{\Gamma }(\frac{3d}{2}\sigma _1\sigma _2\alpha \beta s)}},`$ (3.9)
where $`P^{(\sigma )}(p,m)=1/(p^2m^2+\mathrm{i}ϵ)^\sigma .`$ Closing the integration contour to the left $`(s0)`$ we find the following result for an arbitrary set of indices:

$`V_{1002}(\alpha ,\sigma _1,\sigma _2,\beta ,p^2,m^2,M^2)|_{p^2=m^2}={\displaystyle \frac{(\text{i}^{1d})^2\mathrm{\Gamma }(\sigma _1+\sigma _2\frac{d}{2})\mathrm{\Gamma }(\frac{d}{2}\sigma _1)\mathrm{\Gamma }(\frac{d}{2}\sigma _2)}{\mathrm{\Gamma }^2(1+\epsilon )\mathrm{\Gamma }(\sigma _1)\mathrm{\Gamma }(\sigma _2)\mathrm{\Gamma }(\alpha )\mathrm{\Gamma }(\beta )\mathrm{\Gamma }(d\sigma _1\sigma _2)}}\times `$
$`(M^2)^{d\alpha \beta \sigma _1\sigma _2}[\left({\displaystyle \frac{m^2}{M^2}}\right)^{d\alpha \sigma _1\sigma _2}{\displaystyle \frac{\mathrm{\Gamma }(\beta )\mathrm{\Gamma }(\alpha +\sigma _1+\sigma _2d)\mathrm{\Gamma }(2d2\sigma _12\sigma _2\alpha )}{\mathrm{\Gamma }(\frac{3d}{2}\alpha \sigma _1\sigma _2)}}\times `$
$`{}_{3}{}^{}F_{2}^{}\left(\begin{array}{c}\beta ,d\sigma _1\sigma _2\frac{\alpha }{2},d\sigma _1\sigma _2+\frac{1\alpha }{2}\\ 1+d\alpha \sigma _1\sigma _2,\frac{3d}{2}\alpha \sigma _1\sigma _2\end{array}{\displaystyle \frac{4m^2}{M^2}}\right)`$ (3.12)
$`+{\displaystyle \frac{\mathrm{\Gamma }(\alpha )\mathrm{\Gamma }(\alpha +\beta +\sigma _1+\sigma _2d)\mathrm{\Gamma }(d\alpha \sigma _1\sigma _2)}{\mathrm{\Gamma }(\frac{d}{2})}}{}_{3}{}^{}F_{2}^{}\left(\begin{array}{c}\frac{\alpha }{2},\frac{1+\alpha }{2},\alpha +\beta +\sigma _1+\sigma _2d\\ 1+\alpha +\sigma _1+\sigma _2d,\frac{d}{2}\end{array}{\displaystyle \frac{4m^2}{M^2}}\right)].`$ (3.15)
(3.16)
We introduce the short notation

$$V_{1002}(1,1,1,1,p^2,m^2,M^2)|_{p^2=m^2}V_{1002}$$

for the master integral $`\alpha =\beta =\sigma _1=\sigma _2=1`$. The result for this particular integral reads

$`V_{1002}={\displaystyle \frac{\left(M^2\right)^{2\epsilon }}{2(1\epsilon )(12\epsilon )\epsilon ^2}}{\displaystyle \frac{\mathrm{\Gamma }(1\epsilon )\mathrm{\Gamma }(1+2\epsilon )}{\mathrm{\Gamma }(1+\epsilon )}}\times `$
$`\left[{}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\frac{1}{2}\\ 2\epsilon \end{array}4x\right)x^{12\epsilon }{\displaystyle \frac{2(1\epsilon )(14\epsilon )}{(23\epsilon )(13\epsilon )}}{\displaystyle \frac{\mathrm{\Gamma }(14\epsilon )\mathrm{\Gamma }(1\epsilon )}{\mathrm{\Gamma }(13\epsilon )\mathrm{\Gamma }(12\epsilon )}}{}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\frac{3}{2}2\epsilon \\ 33\epsilon \end{array}4x\right)\right]`$ (3.21)
where $`x=\frac{m^2}{M^2}`$. Using the quadratic transformation for hypergeometric functions

$${}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\frac{1}{2}\\ 2\epsilon \end{array}4x\right)=\frac{2}{1+\sqrt{14x}}{}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\epsilon \\ 2\epsilon \end{array}\frac{1\sqrt{14x}}{1+\sqrt{14x}}\right)(1+\chi ){}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\epsilon \\ 2\epsilon \end{array}|\chi \right)$$

where $`\chi `$ is defined in Eq. (A.2), the first term can be reduced to a new $`{}_{2}{}^{}F_{1}^{}`$ function, whose all-order $`\epsilon `$-expansion is given by Eq. (2.14) of :

$`{}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\epsilon \\ 2\epsilon \end{array}|z\right)`$ $`=`$ $`{\displaystyle \frac{1\epsilon }{2(12\epsilon )z}}\{1+z(1z)^{12\epsilon }`$ (3.26)
$`2(1z)^{12\epsilon }\epsilon {\displaystyle \underset{j=1}{\overset{\mathrm{}}{}}}\epsilon ^j{\displaystyle \underset{k=1}{\overset{j}{}}}(2)^{jk}S_{k,jk+1}(z)\}.`$
Note, that up to order $`\epsilon ^3`$ the expansion of the given hypergeometric function can be extracted from Eq. (A.3) of . The second term can be reduced to a $`{}_{2}{}^{}F_{1}^{}`$ function of the type considered in Appendix A (see Eq. (A.31)). To this end the relation

$$\frac{\left[1+2(ab)\epsilon \right]z}{2(2b\epsilon )}{}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\frac{3}{2}a\epsilon \\ 3b\epsilon \end{array}z\right)=1(1z){}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\frac{3}{2}a\epsilon \\ 2b\epsilon \end{array}z\right)$$

has to be applied. Combining the relations given above we may represent the integral as a series:

$`(M^2)^{2\epsilon }V_{1002}={\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{3}{2\epsilon }}+{\displaystyle \frac{1}{\epsilon }}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\left(\genfrac{}{}{0pt}{}{2n}{n}\right)}{n+1}}x^n\left[\mathrm{ln}x+2S_1(n)2S_1(2n)+{\displaystyle \frac{1}{n+1}}\right]`$
$`+{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\left(\genfrac{}{}{0pt}{}{2n}{n}\right)}{n+1}}x^n\{\mathrm{ln}^2x+[5S_1(n)+4S_1(2n){\displaystyle \frac{3}{n+1}}2]\mathrm{ln}x+\zeta _2`$
$`+3S_2(n)4S_2(2n)+6S_1^2(n)10S_1(n)S_1(2n)+4S_1^2(2n)`$
$`+7{\displaystyle \frac{S_1(n)}{n+1}}6{\displaystyle \frac{S_1(2n)}{n+1}}+{\displaystyle \frac{4}{(n+1)^2}}+4S_1(n)4S_1(2n)+{\displaystyle \frac{2}{n+1}}\}\zeta _2{\displaystyle \frac{7}{2}}`$ (3.27)
or in analytical form:

$`\left({\displaystyle \frac{p^2}{x}}\right)^{2\epsilon }V_{1002}={\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{1}{\epsilon }}\left[{\displaystyle \frac{5}{2}}+\chi \mathrm{ln}x{\displaystyle \frac{(1\chi ^2)}{\chi }}\mathrm{ln}(1+\chi )\right]`$
$`+\chi \mathrm{ln}^2x\mathrm{ln}x\left(4\chi 1+{\displaystyle \frac{1\chi ^2}{\chi }}\left[\mathrm{ln}(1\chi )+2\mathrm{ln}(1+\chi )\right]\right)\zeta _2{\displaystyle \frac{19}{2}}+\chi \zeta _2`$
$`+{\displaystyle \frac{1\chi ^2}{\chi }}\left[4\mathrm{ln}(1+\chi )2\mathrm{ln}(1+\chi )\mathrm{ln}(1\chi )3\text{Li}_2\left(\chi \right)\text{Li}_2\left(\chi \right)2\mathrm{ln}^2(1+\chi )\right].`$
The all order $`\epsilon `$–expansion for this diagram can be deduced for $`x1/4`$. Closing the contour of the integral (3.9) on the right $`(s0)`$ we obtain an alternative representation. The latter also may be obtained from (3.16) by analytical continuation of the hypergeometric function to the inverse argument:

$`4\epsilon ^2(12\epsilon )\left(m^2\right)^{2\epsilon }V_{1002}=\left({\displaystyle \frac{1}{x}}\right)^{12\epsilon }{\displaystyle \frac{\mathrm{\Gamma }(1\epsilon )\mathrm{\Gamma }(1+2\epsilon )}{\mathrm{\Gamma }(1+\epsilon )}}{}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,\epsilon \\ \frac{3}{2}\end{array}{\displaystyle \frac{1}{4x}}\right)`$ (3.31)
$`+{\displaystyle \frac{2}{(13\epsilon )}}{\displaystyle \frac{\mathrm{\Gamma }^2(1\epsilon )\mathrm{\Gamma }(1+2\epsilon )\mathrm{\Gamma }(14\epsilon )}{\mathrm{\Gamma }(1+\epsilon )\mathrm{\Gamma }(12\epsilon )\mathrm{\Gamma }(13\epsilon )}}{}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}1,1+3\epsilon \\ \frac{1}{2}+2\epsilon \end{array}{\displaystyle \frac{1}{4x}}\right)`$ (3.34)
$`{\displaystyle \frac{4\pi \epsilon }{(12\epsilon )}}{\displaystyle \frac{\mathrm{\Gamma }^3(1\epsilon )\mathrm{\Gamma }(14\epsilon )\mathrm{\Gamma }(1+4\epsilon )}{\mathrm{\Gamma }^3(12\epsilon )\mathrm{\Gamma }(1+\epsilon )\mathrm{\Gamma }(1+2\epsilon )}}\sqrt{{\displaystyle \frac{1}{x}}}\left({\displaystyle \frac{2}{x}}\right)^{2\epsilon }\left(1{\displaystyle \frac{1}{4x}}\right)^{\frac{1}{2}\epsilon }.`$ (3.35)
The $`\epsilon `$–expansion of the first hypergeometric function $`{}_{2}{}^{}F_{1}^{}`$ is given in . It only contains log-sine integrals . The second term at first may be transformed using the Kummer relation to a new function (see Eq. (4.13) in ) whose $`\epsilon `$–expansion is given by Eq. (B.13) in  and includes a Lsc function (A.26). The analytical continuation of the $`\epsilon `$–expanded result to $`x>1/4`$ can be constructed in the manner described in Sec. 2.2. of .

The point of interest here is $`x=\mathrm{cos}^2\theta _W`$ $`[p^2=m_W^2]`$. In the rest of this section the expansion parameter $`x`$ will denote $`x=m^2/M^2`$.

* $`𝐕_{\mathrm{𝟐𝟎𝟎𝟏}}`$

For this integral we have the series

$`(M^2)^{2\epsilon }V_{2001}={\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{3}{2\epsilon }}+{\displaystyle \frac{1}{\epsilon }}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\left(\genfrac{}{}{0pt}{}{2n}{n}\right)}{n+1}}x^n\left[\mathrm{ln}x+2S_1(n)2S_1(2n)+{\displaystyle \frac{1}{n+1}}\right]`$
$`+{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\left(\genfrac{}{}{0pt}{}{2n}{n}\right)}{n+1}}x^n\{\mathrm{ln}^2x[3S_1(n)2S_1(2n)+{\displaystyle \frac{2}{n+1}}+2]\mathrm{ln}x3V_2(n)+S_2(n)`$
$`+2S_1^2(n)2S_1(n)S_1(2n)+3{\displaystyle \frac{S_1(n)}{n+1}}2{\displaystyle \frac{S_1(2n)}{n+1}}+{\displaystyle \frac{2}{(n+1)^2}}+4S_1(n)4S_1(2n)+{\displaystyle \frac{2}{n+1}}\}`$
$`+{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}x^n[{\displaystyle \frac{1}{(n+1)^2}}+{\displaystyle \frac{1}{n+1}}{\displaystyle \frac{1}{n}}]\zeta _2{\displaystyle \frac{7}{2}}`$ (3.36)
and the analytical solution

$`\left({\displaystyle \frac{p^2}{x}}\right)^{2\epsilon }V_{2001}={\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{1}{\epsilon }}\left[{\displaystyle \frac{5}{2}}+\chi \mathrm{ln}x{\displaystyle \frac{(1\chi ^2)}{\chi }}\mathrm{ln}(1+\chi )\right]`$
$`+\chi \mathrm{ln}^2x\mathrm{ln}x\left[4\chi +{\displaystyle \frac{(1\chi ^2)}{\chi }}\mathrm{ln}(1\chi ^2)\right]\zeta _2{\displaystyle \frac{19}{2}}`$
$`{\displaystyle \frac{(1\chi ^2)}{2\chi }}\left[4\mathrm{ln}(1\chi )\mathrm{ln}(1+\chi )8\mathrm{ln}(1+\chi )\right]{\displaystyle \frac{(1x)}{x}}\mathrm{ln}(1x)`$
$`{\displaystyle \frac{(1\chi ^2)}{2\chi }}\left[5\text{Li}_2\left(\chi \right)+8\text{Li}_2\left(\chi \right)\text{Li}_2\left(\chi ^3\right)\right]{\displaystyle \frac{1}{2x}}\text{Li}_2\left(x\right).`$ (3.37)
The physical point is $`x=\mathrm{cos}^2\theta _W`$.

* $`𝐅_{\mathrm{𝟐𝟎𝟏𝟎𝟎}}`$

The series and the analytical solution for this integral are given by

$`M^2F_{20100}`$ $`=`$ $`{\displaystyle \frac{1}{x}}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n}}\{{\displaystyle \frac{1}{4}}\mathrm{ln}^2x+{\displaystyle \frac{1}{2}}\mathrm{ln}(x)\mathrm{ln}x+[2S_12\overline{S}_1+{\displaystyle \frac{1}{n}}]\mathrm{ln}x`$ (3.38)
$`+\left[S_1+\overline{S}_1{\displaystyle \frac{1}{n}}\right]\mathrm{ln}(x)\zeta _2{\displaystyle \frac{3}{2}}V_2{\displaystyle \frac{3}{2}}S_2+3\overline{S}_2`$
$`3S_1^2+6S_1\overline{S}_13\overline{S}_1^24{\displaystyle \frac{S_1}{n}}+4{\displaystyle \frac{\overline{S}_1}{n}}{\displaystyle \frac{3}{2}}{\displaystyle \frac{1}{n^2}}2{\displaystyle \frac{1}{\left(\genfrac{}{}{0pt}{}{2n}{n}\right)}}{\displaystyle \frac{1}{n^2}}\}`$
and

$`p^2F_{20100}={\displaystyle \frac{2}{3}}\mathrm{ln}^3(1+\chi )2\zeta _2\mathrm{ln}(1+\chi )\mathrm{ln}^2(1+\chi )\mathrm{ln}(\chi )+\mathrm{ln}\chi \mathrm{ln}(\chi )\mathrm{ln}(1+\chi )`$
$`{\displaystyle \frac{1}{2}}\mathrm{ln}^2\chi \mathrm{ln}(1+\chi )+2\mathrm{ln}\chi \text{Li}_2\left(\chi \right)2\text{Li}_3\left(\chi \right){\displaystyle \frac{1}{6}}\text{Li}_3\left(\chi ^3\right)+{\displaystyle \frac{3}{2}}\text{Li}_3\left(\chi \right){\displaystyle \frac{1}{2}}\text{Li}_3\left(x\right),`$
(3.39)
respectively. The physical point is again $`x=\mathrm{cos}^2\theta _W`$.

* ON-SHELL2

We are interested also in the integrals $`F_{00200}(1,1,1,1,1,p^2,m^2)`$, $`V_{0002}(1,1,1,1,p^2,m^2)`$ and $`J_{000}(1,1,1,p^2)`$ when the external momentum is on-shell. Up to finite parts these diagrams have been calculated in  for arbitrary values of the external momentum. For our purpose it is important, that the limit $`p^2m^2`$ exists and is smooth, such that we can substitute the on-shell values of these diagrams121212We note that the sign of the imaginary parts of the following diagrams: $`(𝐅_{\mathrm{𝟎𝟎𝟏𝟎𝟎}},𝐕_{\mathrm{𝟎𝟏𝟏𝟎}},𝐕_{\mathrm{𝟎𝟎𝟏𝟎}},𝐕_{\mathrm{𝟎𝟎𝟎𝟎}},𝐉_{\mathrm{𝟎𝟎𝟎}})`$ collected in Appendix C of  must be changed in order to get the correct answer. without problems:

$`F_{00100}(1,1,1,1,1,m^2,m^2)={\displaystyle \frac{1}{m^2}}\left[3\zeta _3+2i\pi \zeta _2\right],`$
$`V_{0001}(1,1,1,1,m^2,m^2)=(m^2)^{2\epsilon }\left[{\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{5}{2\epsilon }}\zeta _2{\displaystyle \frac{19}{2}}i\pi \right],`$
$`J_{001}(1,1,1,m^2)=(m^2)^{12\epsilon }\left[{\displaystyle \frac{1}{2\epsilon ^2}}{\displaystyle \frac{5}{4\epsilon }}{\displaystyle \frac{11}{8}}2\zeta _2\right]`$
$`J_{000}(1,1,1,m^2)=(m^2)^{12\epsilon }\left[{\displaystyle \frac{1}{4\epsilon }}+{\displaystyle \frac{13}{18}}i\pi {\displaystyle \frac{1}{2}}\right].`$
* $`𝐅_{\mathrm{𝟐𝟎𝟏𝟏𝟎}}`$

This master-integral shows up as a Higgs contribution to the pole of the top-quark propagator in the approximation of massless $`W`$\- and $`Z`$-bosons and a massive Higgs-boson $`(m_H^2=M^2,m_t^2=m^2M_W^2,M_Z^2)`$. It may be obtained also as the leading term of a Taylor expansion of the original diagram with massive $`W`$\- and $`Z`$-bosons in the limit, when the external momentum is on the heavy mass-shell . For this integral we find the following series representation

$`M^2F_{20110}`$ $`=`$ $`{\displaystyle \frac{1}{x}}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n}}[{\displaystyle \frac{1}{4}}\mathrm{ln}^2x{\displaystyle \frac{1}{n}}\mathrm{ln}x{\displaystyle \frac{1}{2}}\zeta _2`$ (3.40)
$`{\displaystyle \frac{1}{2}}S_2+\overline{S}_2S_1^2+2S_1\overline{S}_1\overline{S}_1^2+{\displaystyle \frac{3}{2}}{\displaystyle \frac{1}{n^2}}],`$
and the analytical result

$`m^2F_{20110}=\zeta _2\mathrm{ln}(1+\chi )+{\displaystyle \frac{1}{2}}\mathrm{ln}^2\chi \mathrm{ln}(1+\chi )+2\mathrm{ln}\chi \text{Li}_2\left(\chi \right)3\text{Li}_3\left(\chi \right).`$ (3.41)
The on-shell result coincides with .

* $`𝐅_{\mathrm{𝟎𝟎𝟏𝟏𝟐}}`$

This master-integral again is important for the calculation of the Higgs correction to the pole-mass of the top-quark. Its representation is

$`M^2F_{00112}`$ $`=`$ $`{\displaystyle \frac{1}{x}}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n}}\{{\displaystyle \frac{1}{2}}\mathrm{ln}^2x+[2S_1+2\overline{S}_13{\displaystyle \frac{1}{n}}]\mathrm{ln}x+\zeta _2`$ (3.42)
$`+S_22\overline{S}_2+2S_1^24S_1\overline{S}_1+2\overline{S}_1^2+6{\displaystyle \frac{S_1}{n}}6{\displaystyle \frac{\overline{S}_1}{n}}+6{\displaystyle \frac{1}{n^2}}\},`$
or, in analytical form,

$`m^2F_{00112}=2\zeta _2\mathrm{ln}(1+\chi )+\mathrm{ln}^2\chi \mathrm{ln}(1+\chi )+2\mathrm{ln}\chi \text{Li}_2\left(\chi \right).`$ (3.43)
The on-shell result again coincides with the one given in .

## 4 Renormalization

The pole mass is a gauge independent and infrared stable quantity. For the SM this has been shown to all orders of perturbation theory . In order to calculate the pole mass in terms of renormalized quantities at the two-loop level we need to calculate the one-loop renormalization constants for all physical parameters (charge and masses) as well as the two-loop mass renormalization constant. Not needed are the wave-function renormalizations and the unphysical sector renormalizations. In order to obtain a gauge invariant result in the SM, however, we have to add in a proper way the tadpole contributions .

We first perform the UV-renormalization within the $`\overline{\mathrm{MS}}`$ scheme in order to obtain finite results. In a next step we work out the relation between the on–shell and $`\overline{\mathrm{MS}}`$ parameters. We adopt the convention that the $`\overline{\mathrm{MS}}`$ parameters are defined by multiplying each $`L`$-loop integral by the factor $`(\mathrm{exp}(\gamma )/4\pi )^{\epsilon L}`$. As usual we denote by $`\mu `$ the $`\overline{\mathrm{MS}}`$ renormalization scale. For the RG functions we use the following definitions : for all dimensionless coupling constants, like $`g,g^{},g_s,e,\lambda `$, the $`\beta `$-function is given by $`\mu ^2\frac{}{\mu ^2}g=\beta _g`$ and for all mass parameters (a mass or the Higgs v.e.v. $`v`$) the anomalous dimension $`\gamma _{m^2}`$ is given by $`\mu ^2\frac{}{\mu ^2}\mathrm{ln}m^2=\gamma _{m^2}.`$

### 4.1 One-loop charge renormalization

The relationship between the bare charge $`e_0`$ and the $`\overline{\mathrm{MS}}`$ charge131313 All $`\overline{\mathrm{MS}}`$–parameters, like $`e,g,g^{}`$ are $`\mu `$-dependent quantities. $`e`$ reads

$$e_0=\mu ^\epsilon e\left(1+Z_{\overline{\mathrm{MS}}}^{(1)}/\epsilon +O(e^4)\right),$$
(4.44)
with

$$Z_{\overline{\mathrm{MS}}}^{(1)}=\frac{e^2}{16\pi ^2}\left(\frac{7}{2}+\frac{2}{3}n_F\left[\frac{5}{9}N_c+1\right]\right),$$
(4.45)
where $`n_F`$ is the number of fermion families and $`N_c`$ the number of colors. In the SM we have $`n_F=3`$ and $`N_c=3`$. The corresponding $`\beta `$-function reads

$$\mu ^2\frac{}{\mu ^2}e=\beta _e=e^3\left(\frac{\beta _g}{g^3}+\frac{\beta _g^{}}{g^3}\right)=\frac{e^3}{16\pi ^2}\left(\frac{7}{2}+\frac{2}{3}n_F\left[\frac{5}{9}N_c+1\right]\right)+O(e^5).$$
(4.46)
where $`g`$ and $`g^{}`$ denote the $`SU(2)_L`$ and $`U(1)_Y`$ gauge couplings, respectively.

### 4.2 Mass renormalization

The mass renormalization constants $`(Z_m)_V`$ at two loops may be written in the form

$$m_{V,0}^2=m_V^2(\mu )(Z_m)_V=m_V^2(\mu )\left(1+\frac{g^2(\mu )}{16\pi ^2\epsilon }Z_V^{(1,1)}+\frac{g^4(\mu )}{(16\pi ^2)^2\epsilon }Z_V^{(2,1)}+\frac{g^4(\mu )}{(16\pi ^2)^2\epsilon ^2}Z_V^{(2,2)}\right),$$
(4.47)
where $`V`$ stands for any of the bosons $`Z`$, $`W`$ or $`H`$. We shall use the same notation for fermions $`Vf`$, where we are interested in particular in the top quark $`f=t`$. In addition to the masses, we have one coupling constant as a free parameter of the SM which we have chosen above to be the electric charge strength $`e=g\mathrm{sin}\theta _W`$. The one-loop mass counter-terms are well known . We divide all corrections into bosonic and fermionic parts, $`Z_V=Z_{V,\text{boson}}+Z_{V,\text{fermion}}`$. The purely bosonic contributions have been given in I. The fermion-loop corrections yield the following $`\overline{\mathrm{MS}}`$ renormalization constants

$`Z_{H,\text{fermion}}^{(1,1)}`$ $`=`$ $`{\displaystyle \underset{\text{leptons}}{}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_l^2}{m_W^2}}+N_c{\displaystyle \underset{\text{quarks}}{}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_q^2}{m_W^2}},`$ (4.48)
$`Z_{W,\text{fermion}}^{(1,1)}`$ $`=`$ $`{\displaystyle \underset{\text{leptons}}{}}\left(2{\displaystyle \frac{m_l^4}{m_W^2m_H^2}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_l^2}{m_W^2}}\right)+N_c{\displaystyle \underset{\text{quarks}}{}}\left(2{\displaystyle \frac{m_q^4}{m_W^2m_H^2}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_q^2}{m_W^2}}\right)`$ (4.49)
$`+{\displaystyle \frac{1}{3}}n_F\left(N_c+1\right),`$
$`Z_{Z,\text{fermion}}^{(1,1)}`$ $`=`$ $`Z_{W,\text{fermion}}^{(1,1)}+{\displaystyle \frac{1}{3}}n_F\left[{\displaystyle \frac{m_Z^2}{m_W^2}}\left({\displaystyle \frac{11}{9}}N_c+3\right)+{\displaystyle \frac{m_W^2}{m_Z^2}}\left({\displaystyle \frac{20}{9}}N_c+4\right){\displaystyle \frac{31}{9}}N_c7\right]`$
(4.50)
$`Z_{t\mathrm{quark}}^{(1,1)}`$ $`=`$ $`{\displaystyle \frac{2}{3}}{\displaystyle \frac{2}{3}}{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{3}{2}}{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}{\displaystyle \frac{3}{4}}{\displaystyle \frac{m_H^2}{m_W^2}}3{\displaystyle \frac{m_W^2}{m_H^2}}+{\displaystyle \frac{3}{4}}{\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{3}{4}}{\displaystyle \frac{m_b^2}{m_W^2}}`$ (4.51)
$`+2{\displaystyle \underset{\mathrm{leptons}}{}}{\displaystyle \frac{m_l^4}{m_W^2m_H^2}}+2N_c{\displaystyle \underset{\mathrm{quarks}}{}}{\displaystyle \frac{m_q^4}{m_W^2m_H^2}},`$
where sums run over all leptons $`(e,\mu ,\tau )`$ (leptons) or quarks $`(u,d,s,c,b,t)`$ (quarks). All masses here are $`\overline{\mathrm{MS}}`$ masses and depend on the renormalization scale $`\mu `$: $`m_V^2=m_V^2(\mu )`$. In general, we will divide the two-loop fermion contribution to the mass-renormalization constants into two parts: the first one includes the contribution from $`N_m`$ massless fermion-families and the second one the contribution of $`N_h`$ lepton-quark families exhibiting one massive u-quark of mass $`m_t`$ (massive top-quark). The number of fermion families is $`n_F=N_m+N_h`$, where in the SM we have $`N_m=2`$ and $`N_h=1`$. The results are given by the following exact expressions

$`Z_{W,\text{fermion}}^{(2,1)}`$ $`=`$ $`[{\displaystyle \frac{m_Z^4}{m_W^4}}({\displaystyle \frac{55}{432}}N_c+{\displaystyle \frac{5}{16}})+{\displaystyle \frac{m_Z^2}{m_W^2}}({\displaystyle \frac{29}{108}}N_c+{\displaystyle \frac{3}{4}})+{\displaystyle \frac{343}{216}}N_c+{\displaystyle \frac{31}{24}}`$ (4.52)
$`{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}\left({\displaystyle \frac{44}{27}}N_c+4\right)+{\displaystyle \frac{m_Z^2}{m_H^2}}\left({\displaystyle \frac{40}{27}}N_c+{\displaystyle \frac{8}{3}}\right)+{\displaystyle \frac{m_W^2}{m_H^2}}{\displaystyle \frac{4}{3}}(N_c+1)`$
$`+{\displaystyle \frac{m_Z^6}{m_W^4m_H^2}}({\displaystyle \frac{22}{27}}N_c+2)]n_F+Z_{W,\text{top}}^{(2,1)}`$
$`Z_{W,\text{fermion}}^{(2,2)}`$ $`=`$ $`[{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_H^2}{m_W^2}}(N_c+1)+{\displaystyle \frac{m_Z^4}{m_W^4}}({\displaystyle \frac{11}{72}}N_c+{\displaystyle \frac{3}{8}}){\displaystyle \frac{m_Z^2}{m_W^2}}({\displaystyle \frac{1}{18}}N_c+{\displaystyle \frac{1}{2}})`$
$`+{\displaystyle \frac{1}{9}}(N_c+1)^2n_F{\displaystyle \frac{15}{4}}N_c{\displaystyle \frac{127}{36}}+{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}\left({\displaystyle \frac{13}{18}}N_c+{\displaystyle \frac{5}{2}}\right)`$
$`{\displaystyle \frac{m_Z^2}{m_H^2}}({\displaystyle \frac{10}{9}}N_c+2){\displaystyle \frac{m_W^2}{m_H^2}}2(N_c+1){\displaystyle \frac{m_Z^6}{m_H^2m_W^2}}({\displaystyle \frac{11}{18}}N_c+{\displaystyle \frac{3}{2}})]n_F`$
$`+Z_{W,\text{top}}^{(2,2)}`$
$`Z_{Z,\text{fermion}}^{(2,1)}`$ $`=`$ $`[{\displaystyle \frac{1}{8}}N_c+{\displaystyle \frac{17}{8}}{\displaystyle \frac{m_Z^2}{m_W^2}}({\displaystyle \frac{73}{216}}N_c+{\displaystyle \frac{19}{8}})+{\displaystyle \frac{m_W^2}{m_Z^2}}({\displaystyle \frac{301}{162}}N_c+{\displaystyle \frac{7}{6}})`$ (4.54)
$`+{\displaystyle \frac{m_W^2}{m_H^2}}{\displaystyle \frac{4}{3}}(N_c+1)+{\displaystyle \frac{m_Z^4}{m_W^4}}\left({\displaystyle \frac{109}{1296}}N_c+{\displaystyle \frac{13}{16}}\right)+{\displaystyle \frac{m_Z^2}{m_H^2}}\left({\displaystyle \frac{40}{27}}N_c+{\displaystyle \frac{8}{3}}\right)`$
$`{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}({\displaystyle \frac{44}{27}}N_c+4)+{\displaystyle \frac{m_Z^6}{m_W^4m_H^2}}({\displaystyle \frac{22}{27}}N_c+2)]n_F+Z_{Z,\text{top}}^{(2,1)}`$
$`Z_{Z,\text{fermion}}^{(2,2)}`$ $`=`$ $`[{\displaystyle \frac{m_H^2}{m_W^2}}({\displaystyle \frac{11}{18}}N_c+{\displaystyle \frac{3}{2}}){\displaystyle \frac{m_H^2m_Z^2}{m_W^4}}({\displaystyle \frac{11}{36}}N_c+{\displaystyle \frac{3}{4}}){\displaystyle \frac{m_H^2}{m_Z^2}}({\displaystyle \frac{5}{9}}N_c+1)+{\displaystyle \frac{29}{108}}N_c`$
$`+{\displaystyle \frac{1}{243}}(27+11N_c)^2n_F{\displaystyle \frac{m_Z^2}{m_W^2}}\left({\displaystyle \frac{77}{108}}N_c+{\displaystyle \frac{7}{4}}+{\displaystyle \frac{1}{243}}(27+11N_c)^2n_F\right)`$
$`+{\displaystyle \frac{m_Z^4}{m_W^4}}\left({\displaystyle \frac{385}{648}}N_c+{\displaystyle \frac{35}{24}}+{\displaystyle \frac{1}{729}}(27+11N_c)^2n_F\right){\displaystyle \frac{m_Z^2}{m_H^2}}\left({\displaystyle \frac{31}{9}}N_c+7\right){\displaystyle \frac{1}{4}}`$
$`{\displaystyle \frac{m_W^2}{m_Z^2}}\left({\displaystyle \frac{308}{81}}N_c+{\displaystyle \frac{28}{9}}+{\displaystyle \frac{8}{729}}(N_c+9)(5N_c+9)\right)+{\displaystyle \frac{m_W^2}{m_H^2}}\left({\displaystyle \frac{13}{9}}N_c+5\right)`$
$`{\displaystyle \frac{m_W^4}{m_Z^2m_H^2}}({\displaystyle \frac{20}{9}}N_c+4)+{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}({\displaystyle \frac{22}{9}}N_c+6){\displaystyle \frac{m_Z^6}{m_W^4m_H^2}}({\displaystyle \frac{11}{9}}N_c+3)]n_F`$
$`+Z_{Z,\text{top}}^{(2,2)},`$
where

$`Z_{W,\text{top}}^{(2,1)}`$ $`=`$ $`[{\displaystyle \frac{5}{2}}{\displaystyle \frac{m_t^6}{m_W^4m_H^2}}{\displaystyle \frac{7}{32}}{\displaystyle \frac{m_t^4}{m_W^4}}{\displaystyle \frac{4}{9}}{\displaystyle \frac{m_t^4}{m_W^2m_H^2}}+{\displaystyle \frac{4}{9}}{\displaystyle \frac{m_t^4m_Z^2}{m_W^4m_H^2}}{\displaystyle \frac{85}{288}}{\displaystyle \frac{m_t^2m_Z^2}{m_W^4}}`$ (4.56)
$`+{\displaystyle \frac{19}{12}}{\displaystyle \frac{m_t^2m_Z^4}{m_W^4m_H^2}}+{\displaystyle \frac{3}{8}}{\displaystyle \frac{m_t^2m_H^2}{m_W^4}}{\displaystyle \frac{20}{3}}{\displaystyle \frac{m_t^2m_Z^2}{m_W^2m_H^2}}{\displaystyle \frac{43}{144}}{\displaystyle \frac{m_t^2}{m_W^2}}+{\displaystyle \frac{35}{6}}{\displaystyle \frac{m_t^2}{m_H^2}}]N_c`$
$`Z_{W,\text{top}}^{(2,2)}`$ $`=`$ $`[4N_c{\displaystyle \frac{m_t^8}{m_W^4m_H^2}}+({\displaystyle \frac{3}{2}}N_c){\displaystyle \frac{m_t^6}{m_W^4m_H^2}}{\displaystyle \frac{39}{16}}{\displaystyle \frac{m_t^4}{m_W^4}}+{\displaystyle \frac{1}{6}}{\displaystyle \frac{m_t^4m_Z^2}{m_W^4m_H^2}}6{\displaystyle \frac{m_t^4m_Z^4}{m_W^4m_H^4}}`$ (4.57)
$`+\left({\displaystyle \frac{2}{3}}(N_c+1)n_F10\right){\displaystyle \frac{m_t^4}{m_W^2m_H^2}}{\displaystyle \frac{1}{48}}{\displaystyle \frac{m_t^2m_Z^2}{m_W^4}}+{\displaystyle \frac{3}{2}}{\displaystyle \frac{m_t^2m_Z^4}{m_W^4m_H^2}}12{\displaystyle \frac{m_t^4}{m_H^4}}`$
$`+3{\displaystyle \frac{m_t^2}{m_H^2}}({\displaystyle \frac{1}{6}}n_F(N_c+1){\displaystyle \frac{73}{24}}){\displaystyle \frac{m_t^2}{m_W^2}}]N_c`$
$`Z_{Z,\text{top}}^{(2,1)}`$ $`=`$ $`[{\displaystyle \frac{5}{2}}{\displaystyle \frac{m_t^6}{m_W^4m_H^2}}{\displaystyle \frac{7}{32}}{\displaystyle \frac{m_t^4}{m_W^4}}{\displaystyle \frac{4}{9}}{\displaystyle \frac{m_t^4}{m_W^2m_H^2}}+{\displaystyle \frac{4}{9}}{\displaystyle \frac{m_t^4m_Z^2}{m_W^4m_H^2}}{\displaystyle \frac{17}{32}}{\displaystyle \frac{m_t^2m_Z^2}{m_W^4}}{\displaystyle \frac{13}{36}}{\displaystyle \frac{m_t^2}{m_Z^2}}`$ (4.58)
$`+{\displaystyle \frac{19}{12}}{\displaystyle \frac{m_t^2m_Z^4}{m_W^4m_H^2}}+{\displaystyle \frac{3}{8}}{\displaystyle \frac{m_t^2m_H^2}{m_W^4}}{\displaystyle \frac{20}{3}}{\displaystyle \frac{m_t^2m_Z^2}{m_W^2m_H^2}}+{\displaystyle \frac{43}{144}}{\displaystyle \frac{m_t^2}{m_W^2}}+{\displaystyle \frac{35}{6}}{\displaystyle \frac{m_t^2}{m_H^2}}]N_c`$
$`Z_{Z,\text{top}}^{(2,2)}`$ $`=`$ $`[4N_c{\displaystyle \frac{m_t^8}{m_W^4m_H^4}}+({\displaystyle \frac{3}{2}}N_c){\displaystyle \frac{m_t^6}{m_W^4m_H^2}}{\displaystyle \frac{39}{16}}{\displaystyle \frac{m_t^4}{m_W^4}}6{\displaystyle \frac{m_t^4m_Z^4}{m_W^4m_H^4}}12{\displaystyle \frac{m_t^4}{m_H^4}}`$ (4.59)
$`\left(4n_F+{\displaystyle \frac{44}{27}}n_FN_c{\displaystyle \frac{11}{3}}\right){\displaystyle \frac{m_t^4}{m_W^2m_H^2}}+\left(2n_F+{\displaystyle \frac{22}{27}}n_FN_C+{\displaystyle \frac{1}{2}}\right){\displaystyle \frac{m_t^4m_Z^2}{m_W^4m_H^2}}`$
$`+\left({\displaystyle \frac{8}{3}}n_F+{\displaystyle \frac{40}{27}}n_FN_c14\right){\displaystyle \frac{m_t^4}{m_Z^2m_H^2}}\left({\displaystyle \frac{1}{2}}n_F+{\displaystyle \frac{11}{54}}n_FN_c+{\displaystyle \frac{5}{48}}\right){\displaystyle \frac{m_t^2m_Z^2}{m_W^4}}`$
$`+{\displaystyle \frac{3}{2}}{\displaystyle \frac{m_t^2m_Z^4}{m_W^4m_H^2}}+3{\displaystyle \frac{m_t^2}{m_H^2}}+\left(n_F+{\displaystyle \frac{11}{27}}n_FN_c{\displaystyle \frac{3}{8}}\right){\displaystyle \frac{m_t^2}{m_W^2}}`$
$`({\displaystyle \frac{2}{3}}n_F+{\displaystyle \frac{10}{27}}n_FN_c{\displaystyle \frac{7}{2}}){\displaystyle \frac{m_t^2}{m_Z^2}}]N_c.`$
The renormalization group equation may be utilized to verify the higher order pole terms $`1/\epsilon ^2`$. It is given by Eq. (4.16) of I

$$\left(Z_V^{(1,1)}\right)^2+\frac{16\pi ^2}{g^2}2\frac{\beta _g^{(1)}Z_V^{(1,1)}}{g}+\underset{i}{}Z_{m_i}^{(1,1)}m_i^2\frac{}{m_i^2}Z_V^{(1,1)}=2Z_V^{(2,2)}.$$

The value of $`\beta _g^{(1)}=\frac{g^3}{16\pi ^2}\left(\frac{43}{12}+\frac{1}{6}(N_c+1)n_F\right)`$ may be calculated from the relation

$$\beta _g^{(1)}\mathrm{sin}\theta _W=\frac{1}{2}g\frac{\mathrm{cos}^2\theta _W}{\mathrm{sin}\theta _W}\left(\frac{g^2}{16\pi ^2}\right)\left(Z_W^{(1,1)}Z_Z^{(1,1)}\right)+\beta _e^{(1)},$$

where $`\beta _e^{(1)}`$ is given in (4.46), and $`\mathrm{cos}^2\theta _W=m_W^2/m_Z^2`$. An additional relation which holds for the $`1/\epsilon `$–terms is the following (see Eq. (4.19) in I):

$$g^2\frac{\beta _g^{}^{(2)}}{g^3}\mathrm{sin}^4\theta _W\frac{\beta _g^{(2)}}{g}\mathrm{sin}^2\theta _W\mathrm{cos}^2\theta _W=\left(\frac{g^2}{16\pi ^2}\right)^2\left(Z_W^{(2,1)}Z_Z^{(2,1)}\right)\mathrm{cos}^2\theta _W.$$

The two-loop $`\beta `$-functions for $`g`$ and $`g^{}`$ are given in  and read

$`\beta _g^{}|_{N_c=3}=\left({\displaystyle \frac{1}{12}}+{\displaystyle \frac{10}{9}}n_F\right){\displaystyle \frac{g^3}{16\pi ^2}}+\left({\displaystyle \frac{1}{4}}+{\displaystyle \frac{95}{54}}n_F\right){\displaystyle \frac{g^5}{(16\pi ^2)^2}}+\left({\displaystyle \frac{3}{4}}+{\displaystyle \frac{1}{2}}n_F\right){\displaystyle \frac{g^3g^2}{(16\pi ^2)^2}}`$
$`{\displaystyle \frac{17}{24}}{\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{g^3g^2}{(16\pi ^2)^2}},`$
$`\beta _g|_{N_c=3}=\left({\displaystyle \frac{43}{12}}+{\displaystyle \frac{2}{3}}n_F\right){\displaystyle \frac{g^3}{16\pi ^2}}+\left({\displaystyle \frac{259}{12}}+{\displaystyle \frac{49}{6}}n_F\right){\displaystyle \frac{g^5}{(16\pi ^2)^2}}+\left({\displaystyle \frac{1}{4}}+{\displaystyle \frac{n_F}{6}}\right){\displaystyle \frac{g^3g^2}{(16\pi ^2)^2}}`$
$`{\displaystyle \frac{3}{8}}{\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{g^5}{(16\pi ^2)^2}}.`$
In all cases we could verify our results to satisfy the above RG equations. We also have established that our two-loop RG equations (4.52)-(4.59), calculated in the broken phase, are related to the ones found in the unbroken theory a long time ago 141414Only recently, these results have been confirmed in  by an independent calculation.. On the one hand the relationship between the UV counterterms and the RG equations in the broken and the unbroken phases for massless coupling constants is a trivial consequence of the Higgs mechanism (spontaneous symmetry breaking), which by definition must preserve the Ward–Takhhashi– and the Slavonv–Taylor–identities. On the other hand the Higgs field vacuum expectation value $`v`$ appears only in the broken phase and its renormalization is looking very different from the ones one has to perform in the symmetric phase. In the broken phase the $`\overline{\mathrm{MS}}`$ parameter $`v(\mu ^2)`$ may be defined by the ratio $`v^2(\mu ^2)=4m_W^2(\mu ^2)\mathrm{sin}^2\theta _W(\mu ^2)/e^2(\mu ^2)`$ and the $`\overline{\mathrm{MS}}`$ Fermi constant151515Note that our definition of the $`\overline{\mathrm{MS}}`$ Fermi constant is different from the one used in . by $`\sqrt{2}G_F(\mu ^2)=1/v^2(\mu ^2)`$. The relationship with the unbroken phase follows if we define, alternatively, $`v^2(\mu ^2)=m^2(\mu ^2)/\lambda (\mu ^2)`$, where $`m^2`$ and $`\lambda `$ are the $`\overline{\mathrm{MS}}`$ parameters of the scalar potential in the symmetric phase. As a consequence the following relations hold between the RG functions of the broken and the unbroken phase:

$$\gamma _{v^2}=\gamma _W\frac{\mathrm{cos}^2\theta }{\mathrm{sin}^2\theta }\left(\gamma _W\gamma _Z\right)2\frac{\beta _e}{e}=\gamma _{m^2}\frac{\beta _\lambda }{\lambda }.$$
(4.60)
Our results allow us to confirm the $`\beta `$–function for the combination $`\gamma _{m^2}\beta _\lambda /\lambda `$, not however, the two-loop RG equations for $`\lambda `$ and $`m^2`$ independently. Given the RG equation for $`v^2`$, we are able to write the two–loop RG equation for the top–quark and Higgs masses using the two–loop result  for the corresponding coupling, i.e. $`\gamma _{m_H^2}=\gamma _{m^2}`$ (for details we refer to ).

These considerations also shed light on the role of the tadpole contributions in parameter renormalizations. The zero momentum transfer Higgs propagator which multiplies the tadpole loops implies (in the perturbation expansion) non–analytic terms proportional to $`1/m_H^21/\lambda `$. In view of (4.60) and the mass coupling relations, like $`m_W=gv/2`$ etc., it is obvious that such terms must appear in mass counter–terms (see above and Sec. 4.2 of I) as well as in the coefficients of the corresponding $`\beta `$–functions.

### 4.3 QCD corrections

At the two-loop level there are also diagrams exhibiting gluon exchange which are contributing to the gauge-boson propagators. The exact analytical result for the $`O(\alpha \alpha _s)`$ contribution of quark-pairs to the vacuum polarization function of the gauge-bosons have been calculated in  (see also ). The corresponding set of 1PI diagrams yields a gauge independent and infrared stable (with respect to taking the limit of vanishing fermion’s masses) contribution161616It is interesting to note, that in contrast to the two–loop electroweak massless fermion contributions, where it suffices to expand all master-integrals up to the finite parts, the master-integrals $`J_{002}`$ and $`J_{022}`$, showing up as additional basic integrals when considering the QCD corrections to the $`W`$\- and $`Z`$-propagators, in this case must be expanded up to terms proportional to $`\epsilon ^2`$ and $`\epsilon `$, respectively. The situation here is similar the one discussed in .. For the sake of RG invariance, in contrast to the usual conventions, we include tadpole contributions which cancel in physical observables like $`\mathrm{\Delta }r`$, $`\mathrm{\Delta }\rho `$ etc.  (see also ). For each quark species, there is one $`O(\alpha \alpha _s)`$ tadpole diagram which is gauge invariant. For a quark with mass $`m_q`$ its bare contribution to the location of the pole of the gauge boson propagators reads

$$\mathrm{\Delta }_t=\frac{g}{16\pi ^2}\frac{g_s^2}{16\pi ^2}16N_cC_f\frac{1}{(d4)^2}\frac{(d1)}{(d3)}\frac{m_q^4}{m_W}(m_q^2)^{2\epsilon },$$

with $`C_f=4/3`$ in the SM. Taking into account the $`O(\alpha _s)`$ term of the $`\overline{\mathrm{MS}}`$ one–loop top–quark mass renormalization constant

$$Z_{\mathrm{top}}^{\alpha _s}=12\frac{\alpha _s}{4\pi }3C_f\frac{1}{\epsilon }$$

we find

$`Z_W^{\alpha _s}`$ $`=`$ $`1+{\displaystyle \frac{g^2}{(16\pi ^2)}}{\displaystyle \frac{\alpha _s}{4\pi }}N_cC_f\left[{\displaystyle \frac{1}{\epsilon }}\left(4{\displaystyle \frac{m_t^4}{m_H^2m_W^2}}{\displaystyle \frac{5}{4}}{\displaystyle \frac{m_t^2}{m_H^2}}+{\displaystyle \frac{1}{2}}n_F\right)+{\displaystyle \frac{1}{\epsilon ^2}}\left(12{\displaystyle \frac{m_t^4}{m_H^2m_W^2}}+{\displaystyle \frac{3}{2}}{\displaystyle \frac{m_t^2}{m_H^2}}\right)\right]`$
$`Z_Z^{\alpha _s}`$ $`=`$ $`1+{\displaystyle \frac{g^2}{(16\pi ^2)}}{\displaystyle \frac{\alpha _s}{4\pi }}N_cC_f[{\displaystyle \frac{1}{\epsilon }}(4{\displaystyle \frac{m_t^4}{m_H^2m_W^2}}{\displaystyle \frac{5}{4}}{\displaystyle \frac{m_t^2}{m_H^2}}+{\displaystyle \frac{10}{9}}n_F{\displaystyle \frac{m_W^2}{m_Z^2}}+{\displaystyle \frac{11}{18}}n_F{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{11}{9}}n_F)`$ (4.61)
$`+{\displaystyle \frac{1}{\epsilon ^2}}(12{\displaystyle \frac{m_t^4}{m_H^2m_W^2}}+{\displaystyle \frac{3}{2}}{\displaystyle \frac{m_t^2}{m_H^2}})],`$
where $`\alpha _s=g_s^2/4\pi `$ and the first five quarks are treated as massless. The terms proportional to $`m_t^4`$ come from the tadpole contribution and will cancel in observable quantities. Again, the renormalization group equations can be used for cross-checking the $`1/\epsilon ^2`$– and $`1/\epsilon `$–terms. The $`1/\epsilon `$–part satisfies the following relation (see Eq. (4.19) of I):

$$\mathrm{sin}^2\theta _W\left[\frac{\beta _g^{}^{\alpha _s}}{g^3}\mathrm{sin}^2\theta _W\frac{\beta _g^{\alpha _s}}{g^3}\mathrm{cos}^2\theta _W\right]=\frac{1}{g^2}\left(Z_{W,\epsilon }^{\alpha _s}Z_{Z,\epsilon }^{\alpha _s}\right)\mathrm{cos}^2\theta _W,$$

where

$`\beta _g^{}^{\alpha _s}|_{N_c=3}={\displaystyle \frac{22}{9}}n_F{\displaystyle \frac{g^3g_s^2}{(16\pi ^2)^2}},`$
$`\beta _g^{\alpha _s}|_{N_c=3}=2n_F{\displaystyle \frac{g^3g_s^2}{(16\pi ^2)^2}}`$
and the $`Z_{V,\epsilon ^j}^{\alpha _s}`$ denote the $`1/\epsilon ^j`$–parts of the $`Z_V^{\alpha _s}`$ renormalization constants (4.61). The terms proportional to $`1/\epsilon ^2`$ may be calculated from the relation

$$\frac{g^2}{16\pi ^2}\left(Z_{\mathrm{top}}^{\alpha _s}1\right)m_t^2\frac{}{m_t^2}Z_V^{(1,1)}=2Z_{V,\epsilon ^2}^{\alpha _s},$$
(4.62)
where the definition of $`Z_V^{(i,j)}`$ has been given in (4.47). Here, we have to take into account that, in the presence of several dimensionless coupling constants, Eqs. (4.14) and (4.15) in I have to be modified accordingly:

$`\gamma _V`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \underset{j}{}}g_j{\displaystyle \frac{}{g_j}}Z_V^{(1)},`$ (4.63)
$`\left(\gamma _V+{\displaystyle \underset{j}{}}\beta _{g_j}{\displaystyle \frac{}{g_j}}+{\displaystyle \underset{i}{}}\gamma _im_i^2{\displaystyle \frac{}{m_i^2}}\right)Z_V^{(n)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \underset{j}{}}g_j{\displaystyle \frac{}{g_j}}Z_V^{(n+1)}.`$ (4.64)
Note that our Eq. (4.62) explicitly reveals, that the systematic inclusion of the tadpole contributions is important for the self–consistency of the renormalization group equations.

### 4.4 $`\overline{\mathrm{MS}}`$ renormalization of the propagator pole

The calculation of the one-loop $`\overline{\mathrm{MS}}`$ renormalized on–shell amplitude $`\widehat{\mathrm{\Pi }}_V^{(1)}`$ is well known (see e.g. ). We get it by rewriting the bare expression in terms of $`\overline{\mathrm{MS}}`$ parameters

$`\widehat{\mathrm{\Pi }}_V^{(1)}`$ $``$ $`\underset{\epsilon 0}{lim}\left(m_{0,V}^2m_V^2+m_{0,V}^2{\displaystyle \frac{g_0^2}{16\pi ^2}}X_{0,V}^{(1)}\right)`$
$`=`$ $`m_V^2(\mu ){\displaystyle \frac{e^2}{16\pi ^2\mathrm{sin}^2\theta _W}}\underset{\epsilon 0}{lim}\left({\displaystyle \frac{1}{\epsilon }}Z_V^{(1,1)}+X_{0,V}^{(1)}\right)=m_V^2(\mu ){\displaystyle \frac{e^2}{16\pi ^2\mathrm{sin}^2\theta _W}}X_V^{(1)}.`$
As explained in Sec. 5 of I, we may avoid the consideration of wave-function renormalization as well as the renormalization of the ghost sector and of the gauge parameters if we look directly at the full two-loop $`\overline{\mathrm{MS}}`$ renormalized on–shell amplitude171717In contrast, one could attempt to renormalize the individual amplitudes $`\mathrm{\Pi }_{0,V}^{(2)}`$ and $`\mathrm{\Pi }_{0,V}^{(1)}^{}`$ in addition to the known $`\mathrm{\Pi }_{0,V}^{(1)}`$.. The latter can be written in the form

$`\{\mathrm{\Pi }_{0,V}^{(2)}+\mathrm{\Pi }_{0,V}^{(1)}\mathrm{\Pi }_{0,V}^{(1)}{}_{}{}^{}\}_{\overline{MS}}=\underset{\epsilon 0}{lim}(\mathrm{\Pi }_{0,V}^{(2)}+\mathrm{\Pi }_{0,V}^{(1)}\mathrm{\Pi }_{0,V}^{(1)}^{}`$
$`+m_V^2(\mu ){\displaystyle \frac{1}{\epsilon }}\left({\displaystyle \frac{e^2}{16\pi ^2\mathrm{sin}^2\theta _W}}\right)^2\left[Z_V^{(1,1)}+\left[{\displaystyle \frac{\mathrm{\Delta }g^2}{g^2}}\right]+{\displaystyle \underset{j}{}}Z_{m_j^2}^{(1,1)}{\displaystyle \frac{}{m_j^2}}\right]X_{0,V}^{(1)}`$
$`+m_V^2(\mu )\left({\displaystyle \frac{e^2}{16\pi ^2\mathrm{sin}^2\theta _W}}\right)^2[{\displaystyle \frac{1}{\epsilon }}Z_V^{(2,1)}+{\displaystyle \frac{1}{\epsilon ^2}}Z_V^{(2,2)}]),`$ (4.66)
where the sum runs over all species of particles $`j=Z,W,H,t`$ and

$$\left[\frac{\mathrm{\Delta }g^2}{g^2}\right]=\frac{\mathrm{cos}^2\theta _W}{\mathrm{sin}^2\theta _W}\left(Z_W^{(1,1)}Z_Z^{(1,1)}\right)\left(7\frac{4}{3}n_F\left[\frac{5}{9}N_c+1\right]\right)\mathrm{sin}^2\theta _W.$$

The functions $`Z_V^{(i,j)},X_V^{(1)}`$ and $`X_{0,V}^{(1)}`$ are defined in (4.49-4.50), (4.52-4.2), (B.2-B.3) and Appendix C of I, respectively. The derivatives of $`X_{0,V}^{(1)}`$ with respect to the masses may be found in Appendix C.

It is interesting to observe that (LABEL:MS1:subtracted) and (4.66) account for the renormalization of the real and the imaginary parts, i.e., for the mass and the width, simultaneously in a unified form.

## 5 Results and discussions

We have calculated the location of the pole in terms of the $`\overline{\mathrm{MS}}`$ –mass for the massive gauge-bosons $`Z`$ and $`W`$. As in I we will write the result in the form

$$\frac{s_P}{m_V^2}=1+\left(\frac{e^2}{16\pi ^2\mathrm{sin}^2\theta _W}\right)X_V^{(1)}+\left(\frac{e^2}{16\pi ^2\mathrm{sin}^2\theta _W}\right)^2X_V^{(2)},$$
(5.1)
where both $`e`$ and $`\mathrm{sin}\theta _W`$ are to be taken in the $`\overline{\mathrm{MS}}`$ scheme. The one-loop coefficients $`X_V^{(1)}`$ for $`Z,W`$ and $`H`$ are known of course as exact results. We wrote them down for completeness in Appendix B. The light (massless) fermion contributions to the coefficients $`X_V^{(2)}`$ are given in Appendix D in exact analytical from in terms of the master–integrals discussed in Sec. 3. Alternatively, these contributions may be represented in a $`\mathrm{sin}^2\theta _W`$–expanded form181818Details concerning this expansion can be found in . in the same way as the other contributions. In the latter case, after summing all diagrams with massless fermion loops, all singularities of type $`\mathrm{ln}^j\mathrm{sin}^2\theta _W`$ cancel which infers the infrared finiteness of the massless fermion contribution to the pole–mass. When the top quark is involved we perform an expansion and present the coefficients as series with respect to three small parameters. One “small” parameter is the weak mixing parameter $`\mathrm{sin}^2\theta _W`$, the others are the mass ratios $`m_V^2/m_H^2`$ and $`m_V^2/m_t^2`$ where $`V=W\mathrm{or}Z`$. The analytical values of the first three coefficients are presented in Appendix E.

Again we have checked, that the relations between the pole– and the $`\overline{\mathrm{MS}}`$ –masses, when calculated in the so-called modified $`\overline{\mathrm{MS}}`$ scheme ($`\overline{\mathrm{MMS}}`$), coincide with the ones obtained in the ordinary $`\overline{\mathrm{MS}}`$ scheme.

We have developed and implemented into a computer program an efficient algorithm for the calculation of the massive Feynman diagrams with external momentum on the mass shell of one of the internal particles. The necessary analytical continuations which allow us to cover the hole region of parameter space has also been constructed and implemented numerically.

After UV renormalization the pole–mass is a finite expression in the limit $`\epsilon 0`$. Since also the IR singularities have been regularized by dimensional regularization, this result implies the infrared finiteness of the SM contributions to the pole–mass. As a consequence, the pole–masses of the gauge–bosons are infrared finite quantities. Note that, in order to establish the gauge invariance of the location of the pole $`s_p`$, the tadpole contributions have to be taken into account. Also RG-invariance requires the inclusion of the tadpoles.

Our calculation proves that the $`\overline{\mathrm{MS}}`$ renormalization scheme, comprised in (4.66), is self consistent and works properly in case of unstable particles.

In contrast to the mass of a stable particle, the definition of a mass of an unstable particle is not unambiguous. The natural definition by $`\sqrt{\mathrm{Re}s_P}`$ (see 2.1,2.2), which we adopted in this paper, is not always used. In particular LEP/SLC experiments have adopted another definition, which has been used by most of the theoretical papers on the $`Z`$ line–shape before LEP/SLC started to take data (see  for a pre–LEP status report). As a consequence the $`Z`$– and $`W`$–masses determined by LEP/SLC experiments and listed in the particle data tables  correspond to $`M_V^{}_{}{}^{}2=M_V^2+\mathrm{\Gamma }_V^2`$ rather than to the pole–mass $`M_V`$191919Note that $`M_Z^{}M_Z35\text{MeV}`$ and thus in order to obtain the pole–mass from the “experimental” $`Z`$–mass one has to subtract about 35 MeV.. In this case the width $`\mathrm{\Gamma }_Z^{}`$ entering the $`Z`$ line–shape must be defined by $`\mathrm{Im}\mathrm{\Pi }_Z(s)s\mathrm{\Gamma }_Z/M_Z`$ and hence is a function of the c.m. energy square $`s`$. While at one–loop order the width is determined solely by the gauge invariant fermion contributions at higher orders a $`s`$–dependent widths, which is not an on–shell quantity, will cause troubles with gauge invariance. However, one still may understand the LEP measurements to be defined from the quantities defined in terms of the gauge invariant position of the pole $`s_P`$ by the gauge invariant relationship $`M_V^{}_{}{}^{}2=M_V^2+\mathrm{\Gamma }_V^2`$. In the latter case the mass may be considered to be defined via $`M_V^{}/M_V=\left|s_P\right|/\mathrm{Re}s_P`$.

Our results for the two–loop mass renormalization constants in the on-shell and the $`\overline{\mathrm{MS}}`$ scheme can be applied for the calculation of physical quantities in both of these schemes at the two-loop level. Examples of such calculations, where results of our paper I have been used, are the computations of the bosonic two–loop contributions to the muon life–time (usually encoded in the correction $`\mathrm{\Delta }r`$ ) presented in .

The numerical evaluation of our results is of interest in several respects. First of all, even though the calculated quantities are not observables, but only a subset of contributions to those, their numerical magnitude is interesting. It certainly sheds light on the size of possible corrections as well as on the convergence of the perturbation expansion. Last but not least, the numerical evaluation by itself is a challenge: the amplitudes are of remarkable complexity and the numerical stability and the efficiency of the evaluations are far from trivial. Although we are largely working with series expansions already, and not with numerical integration of the master integrals, we have to apply numerical multiple precision techniques to be sure that we get the correct answers. To this end numerical evaluations at the present stage have been performed in MAPLE were we can work at the desired precision only on the expense of CPU time.

Since we are calculating the on–shell counter terms for the gauge bosons, it is natural to consider the on–shell scheme with the fine structure constant $`\alpha `$ and the masses as input parameters. For our numerical calculation we have chosen the following input parameters: $`\alpha =1/137.036`$, $`M_W=80.419`$ GeV, $`M_Z=91.188`$ GeV. Neutrinos, leptons and quarks are taken in the massless approximation and the top–quark mass is taken to be $`m_t=174.3\text{GeV}`$. For the strong coupling constant we are using a fixed value $`\alpha _s=\alpha _s(M_Z)=0.1185`$. Often the two–loop corrections are of comparable size to the mass effects of the “light” fermions at one-loop. We then take the following fermion masses for the numerical evaluation:

Fermion masses in GeV    $`m_e`$ = 0.00051 $`m_u`$ = 0.003 $`m_d`$ = 0.006  $`m_\mu `$ = 0.1056583 $`m_c`$ = 1.2 $`m_s`$ = 0.120  $`m_\tau `$ = 1.77703 $`m_t`$ = 174.3 $`m_b`$ = 4.4

The general features of our expansion have been discussed in I. While the series expansion in $`\mathrm{sin}^2\theta _W`$ for the known fixed value converges very well, the expansions in $`M_V^2/m_H^2`$ for large $`m_H`$ leads to a strong coupling problem and perturbation theory becomes useless above about 800 GeV. In general, as is illustrated by Figs. 4 and 5, the first three coefficients yield a good approximation as long as the asymptotic series behaves well. Of course, the expansion in $`M_V^2/m_H^2`$ for smaller values of $`m_H`$ starts to diverge below about 120 GeV. We have calculated six terms in each of the different expansion parameters, but only write down three of them in Appendix E 202020The full set of coefficients may be found at http://www-zeuthen.desy.de/$`\stackrel{~}{}`$kalmykov/pole/pole.html. Figs. 69 show the corrections $`\mathrm{\Delta }_WM_W^2/m_W^2(M_W)1`$ and $`\mathrm{\Delta }_ZM_Z^2/m_Z^2(M_Z)1`$, respectively, as a function of the Higgs mass $`M_H`$ for intermediate and for heavy Higgs masses. The light fermion contributions are small relative to the other corrections. The same is true for the QCD corrections for Higgs masses above about 300 GeV. In Fig. 10 the imaginary part and the absolute value of the pole position $`s_P`$ are shown.

Very often the inverse of (5.1) is required. To that end we have to solve the real part of (2.3) iteratively for $`m_V^2`$ and to express all $`\overline{\mathrm{MS}}`$ parameters in terms of on-shell ones.

The solution to two loops reads

$`m_V^2`$ $`=`$ $`M_V^2\mathrm{Re}\widehat{\mathrm{\Pi }}_V^{(1)}\mathrm{Re}\{\mathrm{\Pi }_V^{(2)}+\mathrm{\Pi }_V^{(1)}\mathrm{\Pi }_V^{(1)}{}_{}{}^{}\}_{\overline{\mathrm{MS}}}`$ (5.2)
$`{\displaystyle \underset{j}{}}(\mathrm{\Delta }m_j^2)^{(1)}{\displaystyle \frac{}{m_j^2}}\mathrm{Re}\widehat{\mathrm{\Pi }}_V^{(1)}(\mathrm{\Delta }e)^{(1)}{\displaystyle \frac{}{e}}\mathrm{Re}\widehat{\mathrm{\Pi }}_V^{(1)}|_{m_j^2=M_j^2,e=e_{\mathrm{OS}}},`$
where the sum runs over all species of particles $`j=Z,W,H,t`$ and

$$(\mathrm{\Delta }m_j^2)^{(1)}=\mathrm{Re}\widehat{\mathrm{\Pi }}_j^{(1)}|_{m_j^2=M_j^2,e=e_{\mathrm{OS}}}M_V^2\frac{e_{\mathrm{OS}}^2}{16\pi ^2\mathrm{sin}^2\theta _W}X_V^{(1)}|_{m_j^2=M_j^2}$$

stands for the self-energy of the $`j`$th particle at $`p^2=m_j^2`$ in the $`\overline{\mathrm{MS}}`$ scheme and parameters replaced by the on-shell ones. The analytical results for the derivative term $`\underset{j}{}\delta _j\frac{}{m_j^2}\widehat{\mathrm{\Pi }}_V^{(1)}`$ can be extracted from the results of Appendix C. In (5.2) $`m_V^2M_V^2=\left(\delta M_V^2\right)_{\mathrm{finite}}`$ represents the finite part of the on–shell mass counter–term. The corresponding bare one is obtained by adding the $`\overline{\mathrm{MS}}`$ mass counter–term (see section 4.2) $`\delta m_V^2=\left(Z_{\overline{MS}}1\right)m_V^2=m_{0V}^2m_V^2`$.

The relation (5.2) involves a change from the $`\overline{\mathrm{MS}}`$ to the on-shell (OS) scheme also for the electric charge. Let us consider therefore, in details, the relationship between the fine structure constants $`\alpha _{\overline{\mathrm{MS}}}`$ and $`\alpha _{\mathrm{OS}}=\alpha `$.

For the electroweak couplings we have to calculate the $`\overline{\mathrm{MS}}`$ versions from their commonly used on-shell values212121The QCD coupling is parameterized almost always in the $`\overline{\mathrm{MS}}`$ scheme. Therefore $`\alpha _s(\mu ^2)`$ usually is directly determined experimentally by fitting data to a suitable perturbative QCD prediction in terms of $`\alpha _s(\mu ^2)`$.. The $`\overline{\mathrm{MS}}`$ version of the fine structure constant is defined as a solution of the renormalization group equation (see 4.46) and can be calculated from the UV-counterterms of the electrical charge $`e`$ (see (4.45)). In perturbation theory the naive relation between $`\overline{\mathrm{MS}}`$ and on-shell values of $`\alpha `$ is determined by the Thomson limit of Compton scattering and can be written as

$`\alpha (\mu ^2)`$ $`=`$ $`\alpha \left\{1+{\displaystyle \frac{\alpha }{4\pi }}\left[7\mathrm{ln}\left({\displaystyle \frac{m_W^2}{\mu ^2}}\right){\displaystyle \frac{2}{3}}+{\displaystyle \frac{4}{3}}{\displaystyle Q_f^2N_{cf}\mathrm{ln}\frac{\mu ^2}{m_f^2}}\right]\right\},`$ (5.3)
where the sum goes over all fermions $`f`$ with $`N_{cf}=1`$ for leptons and $`N_{cf}=3`$ for quarks and $`Q_f`$ is the fermion charge equal to 1 for leptons, -1/3 and 2/3 for d- and u-quarks, correspondingly.

The terms proportional to logarithms of the fermion masses originate from evaluating the finite part of the derivative of the photon propagator at zero momentum transfer, which we denote as $`\mathrm{\Pi }_{\gamma \gamma }^f(0)`$222222At zero momentum fermions contribute to the charge only via $`\mathrm{\Pi }_{\gamma \gamma }(0)`$.. However, the low energy contribution of the five light quarks $`(u,d,s,c,b)`$ cannot be calculated in perturbation theory. The free quark loops are strongly (non-perturbatively) modified by the strong interaction at low energies. In order to evaluate $`\mathrm{\Pi }_{\gamma \gamma }^f(0)`$ we may write it as

$$\mathrm{\Pi }_{\gamma \gamma }^f(0)=\mathrm{Re}\mathrm{\Pi }_{\gamma \gamma }^f(q^2)\left[\mathrm{Re}\mathrm{\Pi }_{\gamma \gamma }^f(q^2)\mathrm{\Pi }_{\gamma \gamma }^f(0)\right],$$
(5.4)
where $`q^2`$ is chosen sufficiently large (typically $`M_Z^2`$) such that the first term on the r.h.s. can be calculated in perturbative QCD. In the limit of large momentum transfer $`sm_f^2`$ the result is

$$\mathrm{Re}\mathrm{\Pi }_{\gamma \gamma }^f(s)=\frac{4}{3}\underset{f}{}Q_f^2N_{cf}\left(\mathrm{ln}\frac{s}{\mu ^2}\frac{5}{3}\right),$$
(5.5)
which we use in the relation (5.3) for the perturbative light quarks contributions $`(s=M_Z^2)`$.

On the grounds of analyticity and unitarity, the second non–perturbative term in (5.4) can be determined by evaluating a dispersion integral over the known experimental $`e^+e^{}\mathrm{hadrons}`$ cross–section (for details see ). Usually, the cross–sections are represented in terms of the cross–section ratio

$$R(s)=\frac{\sigma _{tot}(e^+e^{}\gamma ^{}\mathrm{hadrons})}{\sigma (e^+e^{}\gamma ^{}\mu ^+\mu ^{})},$$

where $`\sigma (e^+e^{}\gamma ^{}\mu ^+\mu ^{})=\frac{4\pi \alpha ^2}{3s}`$ at tree level. In terms of $`R(s)`$ we obtain

$$\mathrm{\Delta }\mathrm{\Pi }_{\gamma \gamma }(q^2)\left[\mathrm{Re}\mathrm{\Pi }_{\gamma \gamma }(q^2)\mathrm{\Pi }_{\gamma \gamma }(0)\right]=\frac{4q^2}{3}_{4m_\pi ^2}^{\mathrm{}}𝑑s\frac{R(s)}{s(sq^2iϵ)}.$$

The shift in the fine structure constant in then $`\mathrm{\Delta }\alpha (s)=4\pi \alpha \times \frac{1}{16\pi ^2}\mathrm{\Delta }\mathrm{\Pi }_{\gamma \gamma }(s)`$. Using the experimental data for $`R(s)`$ up to $`\sqrt{s}=E_{cut}=5`$ GeV and for the $`\mathrm{{\rm Y}}`$ resonances region between 9.6 and 13 GeV and perturbative QCD from 5.0 to 9.6 GeV and for the high energy tail above 13 GeV one gets

$`\mathrm{\Delta }\alpha _{\mathrm{hadrons}}^{(5)}(M_Z^2)`$ $`=`$ $`0.027572\pm 0.000359;\alpha ^1(M_Z^2)=128.952\pm 0.049`$ (5.6)
at $`M_Z=`$ 91.19 GeV. For numerical estimations we use the results of .

The shift in $`\alpha `$ calculated by the dispersion relation corresponds to the on-shell scheme. Accordingly, since $`\widehat{\mathrm{\Pi }}^{(1)}`$ depends on $`e`$ by an overall factor $`e^2`$ only, we have

$$(\mathrm{\Delta }e)^{(1)}\frac{}{e}\widehat{\mathrm{\Pi }}_V^{(1)}=\left\{\delta \alpha _{\mathrm{bos}}+\delta \alpha _{\mathrm{lep}}+\delta \alpha _{\mathrm{top}}+\mathrm{\Delta }\alpha _{\mathrm{hadrons}}^{(5)}(M_Z^2)\delta \mathrm{\Delta }\alpha _{\mathrm{udscb}}(M_Z^2)\right\}\widehat{\mathrm{\Pi }}_V^{(1)},$$

where

$$\delta \alpha _{\mathrm{bos}}=\frac{\alpha }{4\pi }\left(7\mathrm{ln}\frac{M_W^2}{\mu ^2}\frac{2}{3}\right),\delta \alpha _{\mathrm{lep}}=\frac{\alpha }{3\pi }\underset{\mathrm{}=e,\mu ,\tau }{}\mathrm{ln}\frac{m_{\mathrm{}}^2}{\mu ^2},\delta \alpha _{\mathrm{top}}=\frac{4\alpha }{9\pi }\mathrm{ln}\frac{m_t^2}{\mu ^2}$$

and

$$\delta \mathrm{\Delta }\alpha _{\mathrm{udscb}}(M_Z^2)=\frac{11\alpha }{9\pi }\left(\mathrm{ln}\frac{M_Z^2}{\mu ^2}\frac{5}{3}\right)$$

the perturbative subtraction-term (5.5) form the 5 light quarks. The quark loops contributing to the on–shell gauge boson self–energies are evaluated here at a high energy scale and hence it should be save to calculate them in perturbation theory232323Possibilities to evaluate them by non–perturbative methods via dispersion relations have been discussed in ..

Finally, we analyze the Higgs mass dependence of the $`\overline{\mathrm{MS}}`$ parameter $`\mathrm{sin}^2\theta `$ defined as

$$\mathrm{sin}^2\theta _W=1\frac{m_W^2}{m_Z^2}=\frac{g^2}{g^2+g^2}.$$
(5.7)
The Fig. 11 shows the correction to $`\delta _{\mathrm{sin}^2\mathrm{\Theta }}=\mathrm{sin}^2\theta _W^{\overline{\mathrm{MS}}}/\mathrm{sin}^2\theta _W^{\mathrm{OS}}1`$ as a function of the Higgs mass. Since the unphysical $`m_H^4`$ terms drop out by virtue of Veltman’s screening theorem, the corrections do not blow up so dramatically in the strong coupling regime of large Higgs masses. In the region displayed between 150 and 800 GeV the two–loop correction stays below about 0.6% and has a sign opposite to the one–loop correction. The latter is one order of magnitude larger.

## 6 Conclusion

We have calculated the full SM two–loop radiative corrections to the pole masses of the gauge bosons $`W`$ and $`Z`$. A number of conceptual problems dealing with the renormalization of unstable particles could be studied by explicite calculations: i) The position of the complex pole $`s_P`$ of the gauge boson propagators is manifestly gauge invariant after taking into account the Higgs tadpole contributions. ii) The renormalized on-shell self–energies are infrared finite. iii) Our calculation proves that the $`\overline{\mathrm{MS}}`$ renormalization scheme, comprised in (4.66), is self consistent and works properly in case of unstable particles. iv) Up to two–loops we explicitely confirm that the UV singularities and the related RG equations of the broken phase are completely determined by the unbroken phase . v) The inclusion of the tadpoles is also required from the point of view of the renormalization group invariance. vi) Our results for the 2-loop mass renormalization constants in the on-shell and the $`\overline{\mathrm{MS}}`$ scheme can be applied in calculations of physical quantities in both of these schemes at the 2-loop level. vi) A compact expression for the massless fermion contributions in terms of an irreducible set of master–integrals in given. vii) A new technique for calculating the $`\epsilon `$–expansion of some types of hypergeometric functions has been developed.

Note that all results concerning general properties (like gauge invariance etc.) have been checked analytically. Since exact analytical results for many of the multi–scale master–integrals are not yet available, in this work we had to resort to asymptotic expansion techniques. The exact results in terms of a basis of master–integrals will by given elsewhere. For the latter one dimensional integral representations are available which are suitable for the numerical evaluations . However, for numerical evaluations in the heavy Higgs region above about 200 GeV, the computation in terms of our series expansions is numerically much more stable and much more efficient.

Acknowledgments. We are grateful to D. Bardin, A. Davydychev, O. V. Tarasov and B. Tausk for useful discussions. M. K.’s research was supported in part by INTAS-CERN grant No. 99-0377 and by the Australian Research Council grant No. A00000780. We thank Oleg Tarasov for carefully reading the manuscript.

## Appendix

## Appendix A Sums with binomial coefficients in the numerator

Let us present in the following some useful formulae for the binomial sums which occur in the master-integrals (3.38) and (3.40). The sums which involve only harmonic coefficients, were investigated in , while the sums with binomial coefficients in the denominator242424In  sums of this type were called multiple binomial sums. Here we will use the more adequate name inverse multiple binomial sums. were considered in . In this Appendix we present some results for sums with binomial coefficients in the numerator. A sum $`F_a`$ of weight $`a`$ has form

$`F_a(x){\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right)f(k){\displaystyle \frac{x^k}{k^a}}`$ (A.1)
where $`f(k)`$ is a product of some finite sums (e.g. $`S_a,\overline{S}_b`$, etc.). It is convenient to introduce a new variable $`\chi `$ related to $`x`$ via

$`\chi `$ $`=`$ $`{\displaystyle \frac{1\sqrt{14x}}{1+\sqrt{14x}}},x={\displaystyle \frac{\chi }{(1+\chi )^2}},\sqrt{14x}={\displaystyle \frac{1\chi }{1+\chi }}.`$ (A.2)
A general relation for sums of type (A.1) is then given by

$`F_{a1}(x)=x{\displaystyle \frac{d}{dx}}F_a(x)={\displaystyle \frac{1+\chi }{1\chi }}\chi {\displaystyle \frac{d}{d\chi }}\stackrel{~}{F}_a(\chi ).`$ (A.3)
Relevant analytical results for sums of the type considered are:

$`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right){\displaystyle \frac{x^k}{k^4}}=8\text{S}_{1,2}\left(\chi \right)\mathrm{ln}(1+\chi )+4\text{Li}_3\left(\chi \right)\mathrm{ln}(1+\chi )4\text{Li}_2\left(\chi \right)\mathrm{ln}^2(1+\chi )`$
$`{\displaystyle \frac{2}{3}}\mathrm{ln}^4(1+\chi )\text{Li}_4\left(\chi \right)+4\text{S}_{2,2}\left(\chi \right)8\text{S}_{1,3}\left(\chi \right),`$ (A.4)
$`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right){\displaystyle \frac{x^k}{k^2}}S_1={\displaystyle \frac{4}{3}}\mathrm{ln}^3(1+\chi )+{\displaystyle \frac{1}{2}}\text{Li}_3\left(\chi ^2\right)2\mathrm{ln}(1+\chi )\text{Li}_2\left(\chi ^2\right)`$
$`+4\text{S}_{1,2}\left(\chi \right)2\text{S}_{1,2}\left(\chi ^2\right),`$ (A.5)
$`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right){\displaystyle \frac{x^k}{k^2}}\overline{S}_1={\displaystyle \frac{2}{3}}\mathrm{ln}^3(1+\chi )+2\text{Li}_3\left(\chi \right)4\mathrm{ln}(1+\chi )\text{Li}_2\left(\chi \right)`$
$`+4\text{S}_{1,2}\left(\chi \right)+6\text{S}_{1,2}\left(\chi \right)2\text{S}_{1,2}\left(\chi ^2\right),`$ (A.6)
$`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right){\displaystyle \frac{x^k}{k}}S_2={\displaystyle \frac{4}{3}}\mathrm{ln}^3(1+\chi )4\mathrm{ln}(1+\chi )\text{Li}_2\left(\chi \right)4\text{S}_{1,2}\left(\chi \right),`$ (A.7)
$`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right){\displaystyle \frac{x^k}{k}}S_1^2={\displaystyle \frac{4}{3}}\mathrm{ln}^3(1+\chi )+8\mathrm{ln}(1+\chi )\text{Li}_2\left(\chi \right)+4\mathrm{ln}(1+\chi )\text{Li}_2\left(\chi \right)`$
$`4\text{S}_{1,2}\left(\chi \right)+4\text{S}_{1,2}\left(\chi ^2\right),`$ (A.8)
$`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right){\displaystyle \frac{x^k}{k}}S_1\overline{S}_1={\displaystyle \frac{2}{3}}\mathrm{ln}^3(1+\chi )+6\mathrm{ln}(1+\chi )\text{Li}_2\left(\chi \right)6\text{S}_{1,2}\left(\chi \right)`$
$`+4\text{S}_{1,2}\left(\chi \right)+2\text{S}_{1,2}\left(\chi ^2\right),`$ (A.9)
$`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right){\displaystyle \frac{x^k}{k}}\left[\overline{S}_2\overline{S}_1^2\right]={\displaystyle \frac{2}{3}}\mathrm{ln}^3(1+\chi )4\mathrm{ln}(1+\chi )\text{Li}_2\left(\chi \right)8\text{S}_{1,2}\left(\chi \right),`$ (A.10)
$`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2k}{k}}\right){\displaystyle \frac{x^k}{k}}V_2=\text{Li}_3\left(\chi \right)2\text{Li}_3\left(\chi \right)+{\displaystyle \frac{1}{9}}\text{Li}_3\left(\chi ^3\right)\text{Li}_3\left(x\right).`$ (A.11)
All these low-weights sums may be obtained by applying (A.3).

The following representation is often useful for sums of the type just considered:

$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^a}}f(n)={\displaystyle \frac{1}{(a2)!}}{\displaystyle _0^x}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{z^n}{n}}f(n)\left[\mathrm{ln}x\mathrm{ln}z\right]^{a2}{\displaystyle \frac{dz}{z}},`$ (A.12)
where $`f(n)`$ is an arbitrary combination of finite sums. For particular cases, $`f(n)=1,S_1,\overline{S}_1,`$ $`S_2,S_1\overline{S}_1,\overline{S}_2\overline{S}_1^2`$ or $`V_2`$, the corresponding expressions extracted from (A)-(A.11) allow us to perform the sum under the integral in (A.12) exactly and to obtain one-fold integral representations for a large class of binomial sums. For $`0x1/4`$ , a further simplification is possible by substituting, $`z=\frac{1}{4}\mathrm{sin}^2\varphi `$. It leads to the following representation:

$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^a}}f(n)`$ $`=`$ $`{\displaystyle \frac{2}{(a2)!}}{\displaystyle _0^{\mathrm{arcsin}\sqrt{4x}}}d\varphi {\displaystyle \frac{\mathrm{cos}\varphi }{\mathrm{sin}\varphi }}[\mathrm{ln}(4x)\mathrm{ln}(\mathrm{sin}^2\varphi )]^{a2}\times `$ (A.13)
$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{f(n)}{n}}\left({\displaystyle \frac{\mathrm{sin}^2\varphi }{4}}\right)^n`$
$``$ $`\mathrm{\Sigma }(f,\mathrm{arcsin}\sqrt{4x}),`$
where $`\mathrm{\Sigma }(f,\theta )`$ is the short notation for this sum. Unfortunately, we could not find corresponding analytical results for these sums for arbitrary $`x(x1/4)`$ in terms of known functions, like generalized Nielsen polylogarithms. Not even for the simplest case $`f(n)=1`$, which yields

$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^a}}`$ $`=`$ $`{\displaystyle \frac{2^{a+1}}{(a2)!}}{\displaystyle _0^{\mathrm{arcsin}\sqrt{4x}}}\mathrm{ln}\left(\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)\left[\mathrm{ln}(2\sqrt{x})\mathrm{ln}\left(\mathrm{sin}\theta \right)\right]^{a2}d\left(\mathrm{ln}\mathrm{sin}\theta \right),`$
a solution could be found. The problem here may be compared with the one encountered in the context of the inverse binomial sums, considered in .

For $`x>1/4`$ the representation

$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^a}}f(n)`$ $`=`$ $`\mathrm{\Sigma }(f,{\displaystyle \frac{\pi }{2}})+{\displaystyle \frac{1}{(a2)!}}{\displaystyle _{\frac{1}{4}}^x}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{z^k}{n}}f(n)\left[\mathrm{ln}x\mathrm{ln}z\right]^{a2}{\displaystyle \frac{dz}{z}}`$ (A.14)
$`=`$ $`\mathrm{\Sigma }(f,{\displaystyle \frac{\pi }{2}})+{\displaystyle \frac{2^{a2}}{(a2)!}}{\displaystyle _0^{2\mathrm{arccos}\frac{1}{\sqrt{4x}}}}d\theta {\displaystyle \frac{\mathrm{sin}\frac{\theta }{2}}{\mathrm{cos}\frac{\theta }{2}}}[\mathrm{ln}\sqrt{x}+\mathrm{ln}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)]^{a2}\times `$
$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{f(n)}{n}}\left({\displaystyle \frac{1}{4\mathrm{cos}^2\frac{\theta }{2}}}\right)^n`$
is available, where we have introduced a new variable $`z=1/(4\mathrm{cos}^2\frac{\theta }{2})`$ and $`\mathrm{\Sigma }(f,\frac{\pi }{2})`$ is a constant, related to some (combination of) hypergeometric functions $`{}_{P}{}^{}F_{Q}^{}`$ of argument unity. For each given choice of $`f(n)`$ these constants are different, but in all cases they are expressible in terms of “even basis” elements (see  and Appendix B.1 in ). We have checked this statement by high-precision calculations, which allow for “numerical proofs” in the manner explained in . The second term of (A.14) can be calculated analytically for some particular cases. For example, for $`f(n)=1`$ we have

$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^a}}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{}_{a+2}{}^{}F_{a+1}^{}\left(\begin{array}{c}\frac{3}{2},\{1\}_{a+1}\\ \{2\}_{a+1}\end{array}1\right){\displaystyle \frac{2^a}{(a2)!}}{\displaystyle \underset{j=0}{\overset{a2}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{a2}{j}}\right)\left(\mathrm{ln}\sqrt{x}\right)^{a2j}\times `$
$`\{{\displaystyle \frac{1}{j+2}}[\mathrm{ln}^{j+2}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)\mathrm{ln}^{j+2}2]`$
$`+{\displaystyle \frac{\text{i}\sigma }{2(j+1)}}[\text{Ls}_{j+2}(\pi \theta )\text{Ls}_{j+2}\left(\pi \right)\theta \mathrm{ln}^{j+1}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)]\},`$
where we have used $`\underset{k=1}{\overset{\mathrm{}}{}}\left(\genfrac{}{}{0pt}{}{2k}{k}\right)\frac{x^k}{k}=2\mathrm{ln}(1+\chi )`$, $`a2`$, $`\theta =2\mathrm{arccos}\frac{1}{\sqrt{4x}},`$ $`\sigma =\pm 1`$ and $`\text{Ls}_n\left(\theta \right)`$ is so-called log-sine integral  defined by

$$\text{Ls}_j^{(k)}\left(\theta \right)=\underset{0}{\overset{\theta }{}}d\varphi \varphi ^k\mathrm{ln}^{jk1}\left|2\mathrm{sin}\frac{\varphi }{2}\right|,\text{Ls}_j\left(\theta \right)\text{Ls}_j^{(0)}\left(\theta \right).$$
(A.19)
The values of $`\text{Ls}_j\left(\pi \right)`$ can be expressed in terms of $`\zeta `$–function for any $`j`$ . The sign $`\sigma `$ is defined via the relation $`\chi e^{\mathrm{i}\sigma \theta }`$. The choice of the sign derives from the causal “+i0”–prescription for the propagator or as the sign of the square root of $`\mathrm{cos}^2\frac{\theta }{2}`$. In the physical region of interest here we have $`\sigma =1`$.

For $`x>1/4`$, an non–zero imaginary part develops for the sum (LABEL:f=1). It can be calculated in closed form:

$`\mathrm{𝐼𝑚}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^a}}`$ $`=`$ $`\sigma {\displaystyle \frac{2^{a1}}{(a2)!}}{\displaystyle \underset{j=0}{\overset{a2}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{a2}{j}}\right){\displaystyle \frac{\left(\mathrm{ln}\sqrt{x}\right)^{a2j}}{j+1}}\times `$
$`\left[\text{Ls}_{j+2}\left(\pi \theta \right)\text{Ls}_{j+2}\left(\pi \right)\theta \mathrm{ln}^{j+1}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)\right]_{\theta =2\mathrm{arccos}\frac{1}{\sqrt{4x}}}.`$
On the mass–shell $`(x=1,\theta =\frac{2\pi }{3},\sigma =1)`$ we find the result

$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{1}{n^a}}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{}_{a+2}{}^{}F_{a+1}^{}\left(\begin{array}{c}\frac{3}{2},\{1\}_{a+1}\\ \{2\}_{a+1}\end{array}1\right)+{\displaystyle \frac{2^a}{(a2)!}}{\displaystyle \frac{\mathrm{ln}^a2}{a}}`$ (A.24)
$`+\text{i}{\displaystyle \frac{2^{a1}}{(a1)!}}\left[\text{Ls}_a\left({\displaystyle \frac{\pi }{3}}\right)\text{Ls}_a\left(\pi \right)\right].`$
In general, however, analytical results for these types of sums are not expressible in terms of generalized Nielsen polylogarithms. For example, we have

$`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^a}}\overline{S}_1=\mathrm{\Sigma }(\overline{S}_1,{\displaystyle \frac{\pi }{2}})+{\displaystyle \frac{2^{a2}}{(a2)!}}{\displaystyle _0^{2\mathrm{arccos}\frac{1}{\sqrt{4x}}}}d\theta {\displaystyle \frac{\mathrm{sin}\frac{\theta }{2}}{\mathrm{cos}\frac{\theta }{2}}}[\mathrm{ln}\sqrt{x}+\mathrm{ln}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)]^{a2}\times `$
$`\left\{\left[2\zeta _2\pi \theta +{\displaystyle \frac{\theta ^2}{4}}+\mathrm{ln}^2\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)\right]\text{i}\left[2\text{Cl}_2\left(\theta \right)+\theta \mathrm{ln}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)\right]\right\}`$
$`=\mathrm{\Sigma }(\overline{S}_1,{\displaystyle \frac{\pi }{2}}){\displaystyle \frac{2^{a1}}{(a2)!}}{\displaystyle \underset{j=0}{\overset{a2}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{a2}{j}}\right){\displaystyle \frac{\left(\mathrm{ln}\sqrt{x}\right)^{a2j}}{j+1}}\times `$
$`\{[2\zeta _2\pi \theta +{\displaystyle \frac{\theta ^2}{4}}]\mathrm{ln}^{j+1}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)+{\displaystyle \frac{\pi }{2}}[\text{Ls}_{j+2}(\pi \theta )\text{Ls}_{j+2}\left(\pi \right)]2\zeta _2\mathrm{ln}^{j+1}2`$
$`+{\displaystyle \frac{1}{2}}\left[\text{Ls}_{j+3}^{(1)}\left(\pi \theta \right)\text{Ls}_{j+3}^{(1)}\left(\pi \right)\right]+{\displaystyle \frac{j+1}{j+3}}\left[\mathrm{ln}^{j+3}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)\mathrm{ln}^{j+3}2\right]`$
$`2\text{i}\sigma \left[\text{Cl}_2\right(\theta \left)\mathrm{ln}^{j+1}\right(2\mathrm{cos}{\displaystyle \frac{\theta }{2}})\text{Lsc}_{2,j+2}(\theta \left)\right]`$
$`\text{i}\sigma {\displaystyle \frac{j+1}{j+2}}[\text{Ls}_{j+3}\left(\pi \right)\text{Ls}_{j+3}(\pi \theta )+\theta \mathrm{ln}^{j+2}\left(2\mathrm{cos}{\displaystyle \frac{\theta }{2}}\right)]\}_{\theta =2\mathrm{arccos}\frac{1}{\sqrt{4x}}},`$ (A.25)
where $`\underset{k=1}{\overset{\mathrm{}}{}}\left(\genfrac{}{}{0pt}{}{2k}{k}\right)\frac{x^k}{k}\overline{S}_1=2\text{Li}_2\left(\chi \right)+\mathrm{ln}^2(1+\chi )`$ and $`\text{Lsc}_{i,j}\left(\theta \right)`$ are new type of log-sine integrals introduced in  (where the properties of these functions are given in Appendix A.2) and defined by

$$\text{Lsc}_{i,j}\left(\theta \right)=\underset{0}{\overset{\theta }{}}\text{d}\varphi \mathrm{ln}^{i1}\left|2\mathrm{sin}\frac{\varphi }{2}\right|\mathrm{ln}^{j1}\left|2\mathrm{cos}\frac{\varphi }{2}\right|.$$
(A.26)
For $`i,j>2`$ these functions cannot be written in terms of Nielsen polylogarithms, but are related to harmonic polylogarithms which have been considered in .

Using the results for the binomial series given above, we are able to calculate analytically the required first few coefficients of the $`\epsilon `$–expansion of hypergeometric functions of the following type

$${}_{Q}{}^{}F_{P}^{}\left(\begin{array}{c}\frac{3}{2}+b_1\epsilon ,\mathrm{},\frac{3}{2}+b_I\epsilon ,\mathrm{\hspace{0.33em}1}+a_1\epsilon ,\mathrm{},1+a_K\epsilon ,\mathrm{\hspace{0.33em}2}+d_1\epsilon \mathrm{},2+d_L\epsilon \\ \frac{3}{2}+f_1\epsilon ,\mathrm{},\frac{3}{2}+f_J\epsilon ,\mathrm{\hspace{0.33em}1}+e_1\epsilon ,\mathrm{},1+e_R\epsilon ,\mathrm{\hspace{0.33em}2}+c_1\epsilon ,\mathrm{},2+c_{QJR}\epsilon \end{array}z\right).$$
(A.27)
Rewriting this function as an infinite series and using the well-know representation

$$\frac{\mathrm{\Gamma }(j+a\epsilon )}{\mathrm{\Gamma }(1+a\epsilon )}=(j1)!\mathrm{exp}\left[\underset{k=1}{\overset{\mathrm{}}{}}\frac{(a\epsilon )^k}{k}S_k(j1)\right],$$

we obtain (for details we refer to Appendix B of )

$`{}_{P+1}{}^{}F_{P}^{}\left(\begin{array}{c}\{\frac{3}{2}+b_i\epsilon \}_0^{J+1},\{1+a_i\epsilon \}_0^K,\{2+d_i\epsilon \}_0^L\\ \{\frac{3}{2}+f_i\epsilon \}_0^J,\{1+e_q\epsilon \}_0^R,\{2+c_i\epsilon \}_0^{K+LR}\end{array}4z\right)={\displaystyle \frac{1}{2z}}\times `$ (A.30)
$`\mathrm{\Pi }_{s=1}^{K+LR}(1+c_s\epsilon )\mathrm{\Pi }_{i=1}^L{\displaystyle \frac{1}{(1+d_i\epsilon )}}\mathrm{\Pi }_{r=1}^{J+1}{\displaystyle \frac{1}{(1+2b_r\epsilon )}}\mathrm{\Pi }_{k=1}^J(1+2f_k\epsilon ){\displaystyle \underset{j=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2j}{j}}\right){\displaystyle \frac{z^j}{j^{KR1}}}\times \mathrm{\Delta },`$
(A.31)
where $`P=K+L+J`$ and

$`\mathrm{\Delta }=\mathrm{exp}\left({\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(\epsilon )^k}{k}}\{S_kT_k+2^kU_k\overline{S}_k+{\displaystyle \frac{Y_k}{j^k}}\}\right)=(1\epsilon \{S_1T_1+{\displaystyle \frac{Y_1}{j}}+2U_1\overline{S}_1\}`$
$`+\epsilon ^2\{{\displaystyle \frac{1}{2j^2}}[Y_2+Y_1^2]+{\displaystyle \frac{Y_1}{j}}[S_1T_1+2U_1\overline{S}_1]+2U_1S_1\overline{S}_1T_1+2U_2\overline{S}_2+2U_1^2\overline{S}_1^2`$
$`+{\displaystyle \frac{1}{2}}S_2T_2+{\displaystyle \frac{1}{2}}S_1^2T_1^2\}+\mathrm{O}(\epsilon ^3)).`$ (A.32)
Here, we introduced new constants $`A_k,B_k,C_k,D_k,E_k,F_k,T_k,U_k`$: $`A_k_{i=1}^Ka_i^k`$, $`B_k_{i=1}^{J1}b_i^k`$, $`C_k_{i=1}^{K+LR2}c_i^k`$, $`D_k_{i=1}^Ld_i^k`$, $`E_k_{i=1}^Re_i^k`$, $`F_k_{i=1}^Jf_i^k`$, $`T_kB_k+C_k+E_kA_kD_kF_k`$, $`U_kF_kB_k`$, $`Y_kC_kD_k`$. With the expansion just worked out, analytical results are available for the first three coefficients of the $`\epsilon `$–expansion of the hypergeometric function (A.27). In particular, we find

$`{}_{2}{}^{}F_{1}^{}\left(\begin{array}{c}\frac{3}{2}2\epsilon ,1\\ 23\epsilon \end{array}4z\right)={\displaystyle \frac{1}{2z}}{\displaystyle \frac{13\epsilon }{14\epsilon }}{\displaystyle \frac{1+\chi }{1\chi }}[{\displaystyle \frac{2\chi }{1+\chi }}2\{\mathrm{ln}(1\chi )+2\mathrm{ln}(1+\chi )\}\epsilon `$ (A.35)
$`+2\epsilon ^2\{3(\text{Li}_2\left(\chi \right)+2\text{Li}_2(\chi ))+(\mathrm{ln}(1\chi )+2\mathrm{ln}(1+\chi ))^2\}+\mathrm{O}(\epsilon ^3)].`$ (A.36)
For the analytical representation of the $`\epsilon `$–expansion of V-type integrals we need also sums with shifted arguments (see Eqs. (3.27) and (3.36) ). The following two types of sums have to be considered:

$`I_{a_1,\mathrm{},a_p;b_1,\mathrm{},b_q;c}^{i_1,\mathrm{},i_p;j_1,\mathrm{},j_q}(x)`$ $``$ $`{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^c}}[S_{a_1}(n1)]^{i_1}\mathrm{}[S_{a_p}(n1)]^{i_p}[S_{b_1}(2n1)]^{j_1}\mathrm{}[S_{b_q}(2n1)]^{j_q},`$
$`J_{a_1,\mathrm{},a_p;b_1,\mathrm{},b_q;c}^{i_1,\mathrm{},i_p;j_1,\mathrm{},j_q}(x)`$ $``$ $`{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{(n+1)^c}}[S_{a_1}(n)]^{i_1}\mathrm{}[S_{a_p}(n)]^{i_p}[S_{b_1}(2n)]^{j_1}\mathrm{}[S_{b_q}(2n)]^{j_q}.`$ (A.37)
Cases when there are no sums of the type $`S_a(n1)`$ or $`S_b(2n1)`$ on the r.h.s. of (A.37) we characterize by a “$``$”–sign replacing indices $`(a,i)`$ or $`(b,j)`$ of $`I,J`$, respectively.

With this notation we may write the following relations:

$`{\displaystyle \frac{1}{2x}}I_{a;;c}^{j;}(x)=\left[2J_{a;;c}^{j;}(x)J_{a;;c+1}^{j;}(x)\right],`$
$`{\displaystyle \frac{1}{2x}}I_{a;1;c}^{j;1}(x)=2J_{a;1;c}^{j;1}(x)J_{a;1;c+1}^{j;1}(x)+J_{a;;c+1}^{j;}(x),`$
$`{\displaystyle \frac{1}{2x}}\left[I_{;2;c}^{;1}(x)I_{;1;c}^{;2}(x)\right]=2J_{;2;c}^{;1}(x)2J_{;1;c}^{;2}(x)J_{;2;c+1}^{;1}(x)+J_{;1;c+1}^{;2}(x)2J_{;1;c+1}^{;1},`$
$`{\displaystyle \frac{1}{2x}}{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right){\displaystyle \frac{x^n}{n^c}}V_a(n1)={\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}\left({\displaystyle \genfrac{}{}{0pt}{}{2n}{n}}\right)x^n\left[{\displaystyle \frac{2V_a(n)}{(n+1)^c}}{\displaystyle \frac{V_a(n)}{(n+1)^{c+1}}}\right]`$ (A.38)
## Appendix B Fermion corrections to the pole masses at one-loop

In this Appendix we present, for completeness, the well know  one-loop relations between pole- and $`\overline{\mathrm{MS}}`$-masses of the gauge-bosons. We divide all corrections into bosonic (diagrams without any fermions) and fermionic (diagrams exhibiting a fermion loop) one’s: $`X_V^{(1)}=X_{V,\text{boson}}^{(1)}+X_{V,\text{fermion}}^{(1)}`$ where the purely bosonic contributions are given in Appendix B of I. Using the notation

$$\frac{M_V^2}{m_V^2}=1+\left(\frac{e^2}{16\pi ^2\mathrm{sin}^2\theta _W}\right)X_V^{(1)}$$
(B.1)
we may write the fermion corrections in the following form252525For simplicity we assume a diagonal Cabibbo-Kobayashi-Maskawa matrix.

$`X_{W,\text{fermion}}^{(1)}`$ $`=`$ $`{\displaystyle \frac{1}{3}}n_F[{\displaystyle \frac{1}{3}}+{\displaystyle \frac{1}{3}}N_c]+{\displaystyle \underset{lepton}{}}[{\displaystyle \frac{1}{6}}B_0(0,m_l^2;m_W^2)(2{\displaystyle \frac{m_l^2}{m_W^2}}{\displaystyle \frac{m_l^4}{m_W^4}})`$ (B.2)
$`{\displaystyle \frac{1}{6}}{\displaystyle \frac{m_l^4}{m_W^4}}(\mathrm{ln}{\displaystyle \frac{m_l^2}{\mu ^2}}1)(1+12{\displaystyle \frac{m_W^2}{m_H^2}}){\displaystyle \frac{1}{3}}{\displaystyle \frac{m_l^2}{m_W^2}}\mathrm{ln}{\displaystyle \frac{m_l^2}{\mu ^2}}]`$
$`+N_c{\displaystyle \underset{\{u,d\}}{}}[{\displaystyle \frac{1}{3}}{\displaystyle \frac{m_u^2}{m_W^2}}\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}{\displaystyle \frac{1}{3}}{\displaystyle \frac{m_d^2}{m_W^2}}\mathrm{ln}{\displaystyle \frac{m_d^2}{\mu ^2}}`$
$`+{\displaystyle \frac{1}{6}}\left(2+2{\displaystyle \frac{m_u^2m_d^2}{m_W^4}}{\displaystyle \frac{m_u^4}{m_W^4}}{\displaystyle \frac{m_d^4}{m_W^4}}{\displaystyle \frac{m_u^2}{m_W^2}}{\displaystyle \frac{m_d^2}{m_W^2}}\right)B_0(m_d^2,m_u^2;m_W^2)`$
$`+\left({\displaystyle \frac{1}{6}}{\displaystyle \frac{m_u^2m_d^2}{m_W^4}}{\displaystyle \frac{1}{6}}{\displaystyle \frac{m_u^4}{m_W^4}}2{\displaystyle \frac{m_u^4}{m_W^2m_H^2}}\right)\left(1\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}\right)`$
$`+({\displaystyle \frac{1}{6}}{\displaystyle \frac{m_u^2m_d^2}{m_W^4}}{\displaystyle \frac{1}{6}}{\displaystyle \frac{m_d^4}{m_W^4}}2{\displaystyle \frac{m_d^4}{m_W^2m_H^2}})(1\mathrm{ln}{\displaystyle \frac{m_d^2}{\mu ^2}})]`$
$`X_{Z,\text{fermion}}^{(1)}`$ $`=`$ $`{\displaystyle \frac{1}{3}}n_F[2+{\displaystyle \frac{m_Z^2}{m_W^2}}+{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_Z^2}{m_W^2}}B_0(0,0;m_Z^2)+{\displaystyle \frac{4}{3}}{\displaystyle \frac{m_W^2}{m_Z^2}}`$ (B.3)
$`+N_c({\displaystyle \frac{11}{27}}{\displaystyle \frac{m_Z^2}{m_W^2}}+{\displaystyle \frac{20}{27}}{\displaystyle \frac{m_W^2}{m_Z^2}}{\displaystyle \frac{22}{27}})]`$
$`+{\displaystyle \underset{lepton}{}}[4{\displaystyle \frac{m_l^2}{m_Z^2}}(1{\displaystyle \frac{2}{3}}{\displaystyle \frac{m_W^2}{m_Z^2}})(\mathrm{ln}{\displaystyle \frac{m_l^2}{\mu ^2}}B_0(m_l^2,m_l^2;m_Z^2))`$
$`{\displaystyle \frac{m_l^2}{m_W^2}}\left({\displaystyle \frac{5}{3}}\mathrm{ln}{\displaystyle \frac{m_l^2}{\mu ^2}}{\displaystyle \frac{7}{6}}B_0(m_l^2,m_l^2;m_Z^2)\right)2{\displaystyle \frac{m_l^4}{m_H^2m_W^2}}\left(1\mathrm{ln}{\displaystyle \frac{m_l^2}{\mu ^2}}\right)`$
$`+B_0(m_l^2,m_l^2;m_Z^2)({\displaystyle \frac{5}{6}}{\displaystyle \frac{m_Z^2}{m_W^2}}+{\displaystyle \frac{4}{3}}{\displaystyle \frac{m_W^2}{m_Z^2}}2)]`$
$`+N_c{\displaystyle \underset{u}{}}[2{\displaystyle \frac{m_u^4}{m_H^2m_W^2}}(\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}1)`$
$`+\left({\displaystyle \frac{17}{54}}{\displaystyle \frac{m_Z^2}{m_W^2}}+{\displaystyle \frac{16}{27}}{\displaystyle \frac{m_W^2}{m_Z^2}}{\displaystyle \frac{20}{27}}\right)B_0(m_u^2,m_u^2;m_Z^2)`$
$`{\displaystyle \frac{m_u^2}{m_W^2}}\left({\displaystyle \frac{17}{27}}\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}{\displaystyle \frac{7}{54}}B_0(m_u^2,m_u^2;m_Z^2)\right)`$
$`+({\displaystyle \frac{40}{27}}{\displaystyle \frac{m_u^2}{m_Z^2}}{\displaystyle \frac{32}{27}}{\displaystyle \frac{m_u^2m_W^2}{m_Z^4}})(\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}B_0(m_u^2,m_u^2;m_Z^2))]`$
$`+N_c{\displaystyle \underset{d}{}}[2{\displaystyle \frac{m_d^4}{m_H^2m_W^2}}(\mathrm{ln}{\displaystyle \frac{m_d^2}{\mu ^2}}1)`$
$`+\left({\displaystyle \frac{5}{54}}{\displaystyle \frac{m_Z^2}{m_W^2}}+{\displaystyle \frac{4}{27}}{\displaystyle \frac{m_W^2}{m_Z^2}}{\displaystyle \frac{2}{27}}\right)B_0(m_d^2,m_d^2;m_Z^2)`$
$`{\displaystyle \frac{m_d^2}{m_W^2}}\left({\displaystyle \frac{5}{27}}\mathrm{ln}{\displaystyle \frac{m_d^2}{\mu ^2}}+{\displaystyle \frac{17}{54}}B_0(m_d^2,m_d^2;m_Z^2)\right)`$
$`+\left({\displaystyle \frac{4}{27}}{\displaystyle \frac{m_d^2}{m_Z^2}}{\displaystyle \frac{8}{27}}{\displaystyle \frac{m_d^2m_W^2}{m_Z^4}}\right)\left(\mathrm{ln}{\displaystyle \frac{m_d^2}{\mu ^2}}B_0(m_d^2,m_d^2;m_Z^2)\right)`$
$`X_{H,\text{fermion}}^{(1)}`$ $`=`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_l^2}{m_W^2}}{\displaystyle \underset{lepton}{}}\left[B_0(m_l^2,m_l^2;m_H^2)\left(14{\displaystyle \frac{m_l^2}{m_H^2}}\right)4{\displaystyle \frac{m_l^2}{m_H^2}}\left(1\mathrm{ln}{\displaystyle \frac{m_l^2}{\mu ^2}}\right)\right]`$
$`+`$ $`{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_q^2}{m_W^2}}N_c{\displaystyle \underset{quark}{}}\left[B_0(m_q^2,m_q^2;m_H^2)\left(14{\displaystyle \frac{m_q^2}{m_H^2}}\right)4{\displaystyle \frac{m_q^2}{m_H^2}}\left(1\mathrm{ln}{\displaystyle \frac{m_q^2}{\mu ^2}}\right)\right]`$
$`X_{top}^{(1)}`$ $`=`$ $`\left({\displaystyle \frac{20}{9}}{\displaystyle \frac{7}{36}}{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{17}{36}}{\displaystyle \frac{m_Z^4}{m_W^2m_u^2}}+{\displaystyle \frac{10}{9}}{\displaystyle \frac{m_Z^2}{m_u^2}}{\displaystyle \frac{8}{9}}{\displaystyle \frac{m_W^2}{m_u^2}}{\displaystyle \frac{16}{9}}{\displaystyle \frac{m_W^2}{m_Z^2}}\right)B_0(m_Z^2,m_u^2;m_u^2)`$ (B.5)
$`+\left({\displaystyle \frac{1}{4}}+{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_d^4}{m_W^2m_u^2}}+{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_d^2}{m_u^2}}+{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_u^2}{m_W^2}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_d^2}{m_W^2}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_W^2}{m_u^2}}\right)B_0(m_W^2,m_d^2;m_u^2)`$
$`+\left({\displaystyle \frac{m_u^2}{m_W^2}}{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_H^2}{m_W^2}}\right)B_0(m_H^2,m_u^2;m_u^2)+{\displaystyle \frac{155}{36}}{\displaystyle \frac{1}{4}}\mathrm{ln}{\displaystyle \frac{m_W^2}{\mu ^2}}{\displaystyle \frac{14}{9}}\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}`$
$`+{\displaystyle \frac{m_Z^2}{m_W^2}}\left({\displaystyle \frac{1}{2}}{\displaystyle \frac{17}{36}}\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}\right)+{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}\left(13\mathrm{ln}{\displaystyle \frac{m_Z^2}{\mu ^2}}\right)`$
$`+{\displaystyle \frac{m_Z^2}{m_u^2}}\left({\displaystyle \frac{10}{9}}{\displaystyle \frac{17}{36}}{\displaystyle \frac{m_Z^2}{m_W^2}}\right)\left(1\mathrm{ln}{\displaystyle \frac{m_Z^2}{\mu ^2}}\right)+{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_H^2}{m_W^2}}\left(1\mathrm{ln}{\displaystyle \frac{m_H^2}{\mu ^2}}\right)`$
$`+{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_d^2}{m_W^2}}\left(1+{\displaystyle \frac{m_d^2}{m_u^2}}\right)\left(1\mathrm{ln}{\displaystyle \frac{m_d^2}{\mu ^2}}\right)+{\displaystyle \frac{m_W^2}{m_H^2}}\left(13\mathrm{ln}{\displaystyle \frac{m_W^2}{\mu ^2}}\right)`$
$`+{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_d^2}{m_u^2}}\left(1+\mathrm{ln}{\displaystyle \frac{m_W^2}{\mu ^2}}2\mathrm{ln}{\displaystyle \frac{m_d^2}{\mu ^2}}\right)+{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_u^2}{m_W^2}}\left(1\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}\right)`$
$`+{\displaystyle \frac{m_W^2}{m_u^2}}\left({\displaystyle \frac{25}{18}}+{\displaystyle \frac{1}{2}}\mathrm{ln}{\displaystyle \frac{m_W^2}{\mu ^2}}+{\displaystyle \frac{8}{9}}\mathrm{ln}{\displaystyle \frac{m_Z^2}{\mu ^2}}\right){\displaystyle \frac{16}{9}}{\displaystyle \frac{m_W^2}{m_Z^2}}\left(2\mathrm{ln}{\displaystyle \frac{m_u^2}{\mu ^2}}\right)`$
$`2N_c{\displaystyle \underset{quark}{}}{\displaystyle \frac{m_q^4}{m_W^2m_H^2}}\left(1\mathrm{ln}{\displaystyle \frac{m_q^2}{\mu ^2}}\right)2{\displaystyle \underset{lepton}{}}{\displaystyle \frac{m_l^4}{m_W^2m_H^2}}\left(1\mathrm{ln}{\displaystyle \frac{m_l^2}{\mu ^2}}\right)`$
where

$$B_0(m_1^2,m_2^2;p^2)=\underset{0}{\overset{1}{}}𝑑x\mathrm{ln}\left(\frac{m_1^2}{\mu ^2}x+\frac{m_2^2}{\mu ^2}(1x)\frac{p^2}{\mu ^2}x(1x)\text{i}0\right)$$

denotes a scalar two–point function.

## Appendix C Mass renormalization contributions

For the parameter renormalization of a two–loop amplitude one requires the first order derivatives with respect to all relevant parameters of the one–loop amplitude as may be seen in (4.66). Since, we are interested in the on-shell amplitudes mass derivations implicitly involve a differentiation with respect to the external momentum, as the “on–shell momentum” has been given the value of a mass. We restrict ourselves to consider the effect of the mass renormalization contribution, since for the charge renormalization the derivative with respect to the latter is trivial and it is included in (4.66) as a separate term. We thus consider

$$D_V=\underset{j}{}\delta _j\frac{}{m_j^2}X_{0,V}^{(1)}$$
(C.1)
where the coefficients $`\delta _j`$ are the one–loop mass renormalization counter–terms, which depend on the renormalization scheme. For the $`\overline{\mathrm{MS}}`$ scheme they have to be identified with the $`Z_{m_j^2}^{(1,1)}`$ of (4.66). The explicit expressions presented below are written down for the third fermion family $`(\nu _\tau ,\tau ,t,b)`$ in the approximation of vanishing $`\tau `$– and $`b`$–mass. Accordingly in (C.1), the relevant masses are indexed by $`j=W,Z,H,t`$. For simplicity we give the results assuming the Cabibbo-Kobayashi-Maskawa matrix to be the unit matrix. The contribution of a massless family may be obtained by putting $`m_t=0(A_0(m_t)=0)`$. The results read

$`D_W=\delta _t{\displaystyle \frac{N_c}{d1}}\{A_0(m_u^2){\displaystyle \frac{(d2)m_t^4}{4\mathrm{\Delta }(0,m_u,m_W)}}[{\displaystyle \frac{m_W^2}{m_t^2}}(2d5){\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{m_W^4}{m_t^4}}(d2)(d4)]`$
$`+B_0(0,m_u^2;m_W^2){\displaystyle \frac{(d3)m_t^4}{2\mathrm{\Delta }(0,m_u,m_W)}}\left[{\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{m_W^2}{m_t^2}}(2d5)+{\displaystyle \frac{m_W^4}{m_t^2}}(d2)+(d4)\right]`$
$`+B_0(0,m_u^2;m_W^2)({\displaystyle \frac{3d}{2}}{\displaystyle \frac{m_t^2}{m_W^2}})+A_0(m_u^2)[{\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{d}{4}}+{\displaystyle \frac{m_t^2}{m_H^2}}d(d1){\displaystyle \frac{(d2)^2}{4}}]\}`$
$`+\delta _H\{{\displaystyle \frac{(d2)m_H^2m_W^2}{2\mathrm{\Delta }(m_W,m_W,m_H)}}[A_0(m_H^2)B_0(m_H^2,m_W^2;m_W^2){\displaystyle \frac{d3}{d2}}{\displaystyle \frac{1}{2}}A_0(m_W^2)]`$
$`+B_0(m_H^2,m_W^2;m_W^2)\left[{\displaystyle \frac{m_H^2}{m_W^2}}{\displaystyle \frac{1}{4}}+{\displaystyle \frac{m_W^2}{m_H^2}}{\displaystyle \frac{d3}{2}}{\displaystyle \frac{1}{2}}\right]`$
$`2A_0(m_u^2)N_c{\displaystyle \frac{m_t^4}{m_H^4}}+A_0(m_Z^2){\displaystyle \frac{m_Z^4}{m_H^4}}{\displaystyle \frac{d1}{2}}`$
$`+A_0(m_W^2)[{\displaystyle \frac{m_W^2}{m_H^2}}{\displaystyle \frac{2d}{4}}+{\displaystyle \frac{m_W^4}{m_H^4}}(d1)+{\displaystyle \frac{1}{4}}]+A_0(m_H^2)[{\displaystyle \frac{2d}{4}}{\displaystyle \frac{m_H^2}{m_W^2}}{\displaystyle \frac{1}{4}}]\}`$
$`+\delta _Z{\displaystyle \frac{1}{d1}}\{B_0(m_Z^2,m_W^2;m_W^2)[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{d1}{4}}{\displaystyle \frac{m_W^2}{m_Z^2}}(3d5)(d3)`$
$`{\displaystyle \frac{m_W^4}{m_Z^4}}2(d1)(d5)+{\displaystyle \frac{4d^215d+13}{2}}]`$
$`+A_0(m_W^2)\left[{\displaystyle \frac{m_W^2}{m_Z^2}}(2d3)(d2)+{\displaystyle \frac{m_W^4}{m_Z^4}}{\displaystyle \frac{(d1)(d2)(d5)}{d3}}+{\displaystyle \frac{d1}{4}}\right]`$
$`+A_0(m_Z^2)[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{1d}{4}}+{\displaystyle \frac{m_W^2}{m_Z^2}}{\displaystyle \frac{(d5)(d2)(d1)}{2}}{\displaystyle \frac{m_Z^2}{m_H^2}}{\displaystyle \frac{d(d1)^2}{4}}{\displaystyle \frac{(3d5)(d2)}{2}}]\}`$
$`+\delta _W{\displaystyle \frac{1}{d1}}\{A_0(m_u^2){\displaystyle \frac{(d2)m_t^4}{\mathrm{\Delta }(0,m_u,m_W)}}{\displaystyle \frac{N_c}{4}}[{\displaystyle \frac{m_t^4}{m_W^4}}+{\displaystyle \frac{m_t^2}{m_W^2}}(d4)+{\displaystyle \frac{m_W^2}{m_t^2}}(d2)(2d5)]`$
$`+B_0(0,m_u^2;m_W^2){\displaystyle \frac{(d3)m_t^4}{\mathrm{\Delta }(0,m_u,m_W)}}{\displaystyle \frac{N_c}{4}}[(d1){\displaystyle \frac{m_t^4}{m_W^4}}{\displaystyle \frac{m_W^4}{m_t^4}}(d2)`$
$`+({\displaystyle \frac{m_W^2}{m_t^2}}{\displaystyle \frac{m_t^2}{m_W^2}})(d3)]`$
$`+{\displaystyle \frac{(d1)m_H^2m_W^2}{\mathrm{\Delta }(m_W,m_W,m_H)}}[B_0(m_H^2,m_W^2;m_W^2)2(d3)+(A_0(m_W^2)2A_0(m_H^2))(d2)]`$
$`+B_0(0,m_u^2;m_W^2)N_c\left[{\displaystyle \frac{3}{4}}{\displaystyle \frac{m_t^4}{m_W^4}}+{\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{d3}{4}}+{\displaystyle \frac{d2}{4}}\right]B_0(0,0;m_W^2){\displaystyle \frac{(d2)(d4)}{4}}`$
$`B_0(m_H^2,m_W^2;m_W^2)\left[{\displaystyle \frac{d}{8}}{\displaystyle \frac{m_H^4}{m_W^4}}{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_H^2}{m_W^2}}+{\displaystyle \frac{(d1)(d4)}{2}}\right]{\displaystyle \frac{1}{2}}N_cA_0(m_u^2){\displaystyle \frac{m_t^4}{m_W^4}}`$
$`B_0(m_Z^2,m_W^2;m_W^2)\left[{\displaystyle \frac{d}{8}}{\displaystyle \frac{m_Z^4}{m_W^4}}+{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{(2d5)(d1)}{2}}+6(d1){\displaystyle \frac{m_W^2}{m_Z^2}}+{\displaystyle \frac{d^2+3d8}{2}}\right]`$
$`A_0(m_W^2)[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{d}{8}}+{\displaystyle \frac{m_H^2}{m_W^2}}{\displaystyle \frac{d}{8}}+{\displaystyle \frac{m_W^2}{m_H^2}}{\displaystyle \frac{d(d1)^2}{2}}{\displaystyle \frac{m_W^2}{m_Z^2}}3{\displaystyle \frac{(d1)(d2)}{d3}}`$
$`{\displaystyle \frac{d2}{d3}}{\displaystyle \frac{d^313d^2+39d31}{2}}]+A_0(m_H^2)[{\displaystyle \frac{m_H^4}{m_W^4}}{\displaystyle \frac{d}{8}}+{\displaystyle \frac{(d1)(d2)}{2}}]`$
$`+A_0(m_Z^2)[{\displaystyle \frac{m_Z^4}{m_W^4}}{\displaystyle \frac{d}{8}}+{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{(2d3)(d2)}{2}}+{\displaystyle \frac{3(d1)(d2)}{2}}]\}`$ (C.2)
$`D_Z=\delta _t{\displaystyle \frac{N_c}{d1}}\{A_0(m_u^2){\displaystyle \frac{m_t^4(d2)}{\mathrm{\Delta }(m_u,m_u,m_Z)}}[{\displaystyle \frac{m_Z^4}{m_W^2m_t^2}}{\displaystyle \frac{9d43}{18}}{\displaystyle \frac{m_Z^6}{m_W^2m_t^4}}{\displaystyle \frac{17(d2)}{36}}`$
$`{\displaystyle \frac{m_W^2m_Z^2}{m_t^4}}{\displaystyle \frac{8(d2)}{9}}{\displaystyle \frac{m_W^2}{m_t^2}}{\displaystyle \frac{32}{9}}+{\displaystyle \frac{m_Z^2}{m_t^2}}{\displaystyle \frac{40}{9}}+{\displaystyle \frac{m_Z^4}{m_t^4}}{\displaystyle \frac{10(d2)}{9}}]`$
$`+B_0(m_u^2,m_u^2;m_Z^2){\displaystyle \frac{m_t^4(d3)}{\mathrm{\Delta }(m_u,m_u,m_Z)}}[{\displaystyle \frac{m_Z^6}{m_W^2m_t^4}}{\displaystyle \frac{17(d2)}{18}}{\displaystyle \frac{m_Z^4}{m_W^2m_t^2}}{\displaystyle \frac{(9d43)}{9}}`$
$`+{\displaystyle \frac{m_Z^2m_W^2}{m_t^4}}{\displaystyle \frac{16(d2)}{9}}+{\displaystyle \frac{m_W^2}{m_t^2}}{\displaystyle \frac{64}{9}}{\displaystyle \frac{m_Z^2}{m_t^2}}{\displaystyle \frac{80}{9}}{\displaystyle \frac{m_Z^4}{m_t^4}}{\displaystyle \frac{20(d2)}{9}}]`$
$`+A_0(m_u^2)\left[{\displaystyle \frac{m_Z^2m_t^2}{m_W^2m_H^2}}d(d1){\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{17(d2)^2}{36}}{\displaystyle \frac{m_W^2}{m_Z^2}}{\displaystyle \frac{8(d2)^2}{9}}+{\displaystyle \frac{10(d2)^2}{9}}\right]`$
$`B_0(m_u^2,m_u^2;m_Z^2)[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{9d43}{18}}{\displaystyle \frac{32}{9}}{\displaystyle \frac{m_W^2}{m_Z^2}}+{\displaystyle \frac{40}{9}}]\}`$
$`+\delta _H\{B_0(m_H^2,m_Z^2;m_Z^2){\displaystyle \frac{m_Z^4}{2\mathrm{\Delta }(m_Z,m_Z,m_H)}}{\displaystyle \frac{m_H^2}{m_W^2}}(3d)2A_0(m_u^2)N_c{\displaystyle \frac{m_t^4m_Z^2}{m_H^4m_W^2}}`$
$`B_0(m_H^2,m_Z^2;m_Z^2)\left[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{1}{2}}+{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}{\displaystyle \frac{3d}{2}}{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_H^2}{m_W^2}}\right]+A_0(m_W^2){\displaystyle \frac{m_W^2m_Z^2}{m_H^4}}(d1)`$
$`+A_0(m_H^2){\displaystyle \frac{m_Z^4}{\mathrm{\Delta }(m_Z,m_Z,m_H)}}{\displaystyle \frac{m_H^2}{m_W^2}}{\displaystyle \frac{d2}{2}}A_0(m_Z^2){\displaystyle \frac{m_Z^4}{\mathrm{\Delta }(m_Z,m_Z,m_H)}}{\displaystyle \frac{m_H^2}{m_W^2}}{\displaystyle \frac{d2}{4}}`$
$`+A_0(m_Z^2)[{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_Z^2}{m_W^2}}+{\displaystyle \frac{d1}{2}}{\displaystyle \frac{m_Z^6}{m_H^4m_W^2}}{\displaystyle \frac{d2}{4}}{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}]A_0(m_H^2)[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{d2}{4}}+{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_H^2}{m_W^2}}]\}`$
$`+\delta _W{\displaystyle \frac{1}{d1}}\{B_0(m_H^2,m_Z^2;m_Z^2)[{\displaystyle \frac{m_Z^2m_H^2}{m_W^4}}(d1){\displaystyle \frac{m_Z^4}{m_W^4}}{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_H^4}{m_W^4}}]`$
$`+B_0(0,0;m_Z^2)(d2)\left[{\displaystyle \frac{2}{9}}N_c{\displaystyle \frac{5}{36}}N_c{\displaystyle \frac{m_Z^4}{m_W^4}}{\displaystyle \frac{3}{2}}{\displaystyle \frac{m_Z^4}{m_W^4}}+2\right]`$
$`A_0(m_u^2)N_c\left[2(d1){\displaystyle \frac{m_t^4m_Z^2}{m_W^4m_H^2}}{\displaystyle \frac{17}{18}}(d2){\displaystyle \frac{m_Z^2m_t^2}{m_W^4}}+{\displaystyle \frac{m_t^2}{m_Z^2}}{\displaystyle \frac{16(d2)}{9}}\right]`$
$`+B_0(m_u^2,m_u^2;m_Z^2)N_c\left[{\displaystyle \frac{m_Z^2m_t^2}{m_W^4}}{\displaystyle \frac{9d43}{18}}{\displaystyle \frac{17(d2)}{36}}{\displaystyle \frac{m_Z^4}{m_W^4}}+{\displaystyle \frac{32}{9}}{\displaystyle \frac{m_t^2}{m_Z^2}}+{\displaystyle \frac{8(d2)}{9}}\right]`$
$`+A_0(m_W^2)\left[{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_Z^2}{m_W^2}}+{\displaystyle \frac{m_W^2}{m_Z^2}}(d1)(d2)(d+1){\displaystyle \frac{m_Z^2}{m_H^2}}{\displaystyle \frac{(d1)^2(d2)}{2}}+(d1)(d2)\right]`$
$`+A_0(m_Z^2)[{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_Z^4}{m_W^4}}{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_Z^2m_H^2}{m_W^4}}+{\displaystyle \frac{m_Z^6}{m_H^2m_W^4}}{\displaystyle \frac{(d1)^2}{2}}]+A_0(m_H^2)[{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_H^4}{m_W^4}}+{\displaystyle \frac{(d1)}{2}}{\displaystyle \frac{m_Z^2m_H^2}{m_W^4}}]`$
$`B_0(m_W^2,m_W^2;m_Z^2)[{\displaystyle \frac{1}{4}}{\displaystyle \frac{m_Z^4}{m_W^4}}+{\displaystyle \frac{d3}{2}}{\displaystyle \frac{m_Z^2}{m_W^2}}+2(d^21){\displaystyle \frac{m_W^2}{m_Z^2}}+(4d7)(d1)]\}`$
$`+\delta _Z{\displaystyle \frac{1}{d1}}\{{\displaystyle \frac{m_Z^4(d1)(d3)}{\mathrm{\Delta }(m_Z,m_Z,m_H)}}{\displaystyle \frac{m_H^2}{m_W^2}}[2B_0(m_H^2,m_Z^2;m_Z^2)+A_0(m_Z^2)2A_0(m_H^2)]`$
$`B_0(m_H^2,m_Z^2;m_Z^2)\left[{\displaystyle \frac{m_H^4}{m_W^2m_Z^2}}{\displaystyle \frac{d2}{8}}+{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{(d1)(d6)}{2}}+{\displaystyle \frac{1}{2}}{\displaystyle \frac{m_H^2}{m_W^2}}\right]`$
$`+B_0(0,0;m_Z^2)[{\displaystyle \frac{m_Z^2}{m_W^2}}N_c{\displaystyle \frac{5d(d2)}{72}}+{\displaystyle \frac{m_W^2}{m_Z^2}}N_c{\displaystyle \frac{(d2)(d4)}{9}}N_c{\displaystyle \frac{(d2)^2}{18}}`$
$`+{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{3d(d2)}{4}}+{\displaystyle \frac{m_W^2}{m_Z^2}}(d2)(d4){\displaystyle \frac{3(d2)^2}{2}}]`$
$`A_0(m_u^2){\displaystyle \frac{N_cm_t^4}{\mathrm{\Delta }(m_u,m_u,m_Z)}}[{\displaystyle \frac{m_Z^4}{m_W^2m_t^2}}{\displaystyle \frac{(d2)(9d43)}{72}}{\displaystyle \frac{m_Z^6}{m_W^2m_t^4}}{\displaystyle \frac{17(d2)^2}{144}}`$
$`{\displaystyle \frac{m_W^2m_Z^2}{m_t^4}}{\displaystyle \frac{2(d2)^2}{9}}{\displaystyle \frac{m_W^2}{m_t^2}}{\displaystyle \frac{8(d2)}{9}}+{\displaystyle \frac{m_Z^2}{m_t^2}}{\displaystyle \frac{10(d2)}{9}}+{\displaystyle \frac{m_Z^4}{m_t^4}}{\displaystyle \frac{5(d2)^2}{18}}]`$
$`+B_0(m_u^2,m_u^2;m_Z^2){\displaystyle \frac{N_cm_t^4}{\mathrm{\Delta }(m_u,m_u,m_Z)}}[{\displaystyle \frac{m_Z^4}{m_W^2m_t^2}}{\displaystyle \frac{(d3)(9d43)}{36}}{\displaystyle \frac{m_Z^6}{m_W^2m_t^4}}{\displaystyle \frac{17(d2)(d3)}{72}}`$
$`{\displaystyle \frac{m_W^2m_Z^2}{m_t^4}}{\displaystyle \frac{4(d2)(d3)}{9}}{\displaystyle \frac{m_W^2}{m_t^2}}{\displaystyle \frac{16(d3)}{9}}+{\displaystyle \frac{m_Z^2}{m_t^2}}{\displaystyle \frac{20(d3)}{9}}+{\displaystyle \frac{m_Z^4}{m_t^4}}{\displaystyle \frac{5(d2)(d3)}{9}}]`$
$`+A_0(m_u^2)N_c[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{17(d2)^2}{144}}+{\displaystyle \frac{m_t^4}{m_W^2m_H^2}}2(d1){\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{(d2)(9d+25)}{72}}`$
$`+{\displaystyle \frac{m_W^2m_t^2}{m_Z^4}}{\displaystyle \frac{8(d2)}{3}}+{\displaystyle \frac{m_W^2}{m_Z^2}}{\displaystyle \frac{2(d2)^2}{9}}{\displaystyle \frac{m_t^2}{m_Z^2}}{\displaystyle \frac{10(d2)}{9}}{\displaystyle \frac{5(d2)^2}{18}}]`$
$`+B_0(m_u^2,m_u^2;m_Z^2)N_c[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{17(d2)}{24}}{\displaystyle \frac{m_t^2}{m_W^2}}{\displaystyle \frac{9d43}{36}}{\displaystyle \frac{16}{3}}{\displaystyle \frac{m_W^2m_t^2}{m_Z^4}}{\displaystyle \frac{m_W^2}{m_Z^2}}{\displaystyle \frac{4(d2)}{9}}`$
$`+{\displaystyle \frac{20}{9}}{\displaystyle \frac{m_t^2}{m_Z^2}}{\displaystyle \frac{5(d2)}{9}}]`$
$`A_0(m_W^2)\left[{\displaystyle \frac{d}{4}}+{\displaystyle \frac{m_W^2}{m_Z^2}}(2d3)(d2)+{\displaystyle \frac{m_W^2}{m_H^2}}(d1)^2+{\displaystyle \frac{m_W^4}{m_Z^4}}3(d1)(d2)\right]`$
$`A_0(m_Z^2)\left[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{d}{4}}+{\displaystyle \frac{m_Z^4}{m_W^2m_H^2}}{\displaystyle \frac{(d1)^2(d+2)}{4}}+{\displaystyle \frac{m_H^2}{m_W^2}}{\displaystyle \frac{d2}{8}}\right]`$
$`+A_0(m_H^2)\left[{\displaystyle \frac{m_H^4}{m_W^2m_Z^2}}{\displaystyle \frac{d2}{8}}+{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{(d1)(d2)}{2}}{\displaystyle \frac{m_H^2}{m_W^2}}{\displaystyle \frac{d1}{2}}\right]`$
$`+B_0(m_W^2,m_W^2;m_Z^2)[{\displaystyle \frac{m_Z^2}{m_W^2}}{\displaystyle \frac{d}{8}}+{\displaystyle \frac{m_W^2}{m_Z^2}}{\displaystyle \frac{d^2+3d8}{2}}+6(d1){\displaystyle \frac{m_W^4}{m_Z^4}}+{\displaystyle \frac{(2d5)(d1)}{2}}]\}`$ (C.3)
where

$$\mathrm{\Delta }(m_1,m_2,m_3)=2m_1^2m_2^2+2m_1^2m_3^2+2m_2^2m_3^3m_1^4m_2^4m_3^4.$$

and

$`A_0(m_1^2)`$ $`=`$ $`{\displaystyle \frac{1}{m_1^2}}{\displaystyle \frac{\text{d}^dq}{\pi ^{d/2}\mathrm{\Gamma }(1+\epsilon )}\frac{1}{\left(q^2+m_1^2\right)}},`$
$`B_0(m_1^2,m_2^2;p^2)`$ $`=`$ $`{\displaystyle \frac{\text{d}^dq}{\pi ^{d/2}\mathrm{\Gamma }(1+\epsilon )}\frac{1}{\left(q^2+m_1^2\right)\left((pq)^2+m_2^2\right)}}.`$ (C.4)
## Appendix D Bare two–loop contribution of the massless fermions

As one of our results we present the exact analytic two–loop contribution to the on–shell self–energies of the gauge bosons. After reduction of the set of basic integrals to a minimal set of master–integrals by means of Tarasov’s recurrence relations  we obtain the following expressions:

$`(\mathrm{\Pi }_W^{(2)}+\mathrm{\Pi }_W^{(1)}\mathrm{\Pi }_W^{(1)}{}_{}{}^{})=`$
$`V_{00W}\{{\displaystyle \frac{\left(d2\right)}{4c^2d\left(d4\right)\left(d3\right)\left(d1\right)^2}}[8c^2d^6+12c^4d^5120c^2d^5+596c^2d^4`$
$`+d^496c^4d^41136c^2d^3+268c^4d^311d^3+508c^2d^2312c^4d^2+40d^248d`$
$`+528c^2d+128c^4d384c^2]+{\displaystyle \frac{\left(d2\right)}{4m_H^2m_W^2\left(m_H^24m_W^4\right)\left(d1\right)^2}}[32d^3m_W^6`$
$`8d^3m_W^4m_H^2+28d^2m_W^4m_H^2128d^2m_W^6+m_H^6d4m_H^4dm_W^228dm_W^4m_H^2`$
$`+160dm_W^6+16m_H^4m_W^24m_H^664m_W^6+8m_W^4m_H^2]\}`$
$`+V_{00Z}\{{\displaystyle \frac{\left(d2\right)}{12\left(d1\right)^2c^4d}}[80c^6d^4+64c^8d^4+40c^4d^488c^2d^3+4d^3+520c^6d^3`$
$`480c^8d^3124c^4d^3+434c^2d^2236c^4d^2872c^6d^251d^2+992c^8d^2+78d`$
$`508c^2d576c^8d+632c^4d+272c^6d+32c^21664c^4]`$
$`{\displaystyle \frac{\left(d2\right)^2\left(5m_Z^2+8c^2m_W^210c^2m_Z^2\right)}{3c^2m_H^2}}\}`$
$`F_{0000Z}{\displaystyle \frac{m_W^4\left(d2\right)\left(1+2c^2+2c^4\right)\left(2c^4c^2d+8c^2+2\right)}{6\left(d1\right)c^6}}`$
$`F_{Z0W00}8{\displaystyle \frac{m_W^4\left(c^2+2\right)\left(d2\right)}{d1}}`$
$`+V_{W00Z}{\displaystyle \frac{m_W^2\left(d2\right)}{12\left(d1\right)^2c^6}}[120c^2d^2+500c^4d^2552c^6d^2+96c^8d^2+256c^{10}d^2`$
$`1408c^{10}d15d288c^8d+2856c^6d2444c^4d+570c^2d640c^2+30+2712c^4`$
$`3456c^6+832c^8+1152c^{10}]`$
$`+V_{Z00W}{\displaystyle \frac{m_W^2\left(d2\right)}{4\left(d1\right)^2c^4}}[24c^6d^28c^2d^2+44c^4d^2d80c^6d+42c^2d204c^4d`$
$`+458c^2+208c^4+56c^6]`$
$`+V_{H00W}\{m_W^2{\displaystyle \frac{\left(d2\right)^2}{d1}}{\displaystyle \frac{\left(d2\right)}{4m_W^2\left(m_H^24m_W^2\right)\left(d1\right)^2}}[8d^2m_W^6+32dm_W^6+m_H^6d`$
$`10m_H^4dm_W^2+24dm_W^4m_H^24m_H^624m_W^672m_W^4m_H^2+34m_H^4m_W^2]\}`$
$`J_{00Z}{\displaystyle \frac{1}{12\left(d1\right)^2c^4d}}[256c^2d302c^2d^2432d64c^2+432d^2+128c^420c^4d`$
$`206c^4d^2+1288c^2d^4+103c^2d^356c^4d^4+220c^4d^3+12d^4140d^3+52c^6d^4`$
$`612c^6d^32240c^6d+2008c^6d^2+768c^6]`$
$`J_{00Z}^{}{\displaystyle \frac{m_W^2\left(c1\right)\left(c+1\right)}{3d\left(d1\right)^2c^6}}[72c^6d+4c^6d^3+36c^6d^2+32c^6+24c^4d+20c^4d^3+16c^4`$
$`48c^4d^215c^2d^28c^2+2c^2d^3+18c^2d+18d^224d2d^3+8]`$
$`+J_{00H}\left\{{\displaystyle \frac{d2}{2\left(d1\right)^2}}+{\displaystyle \frac{\left(d2\right)\left(3d8\right)\left(4dm_W^44m_W^2m_H^24m_W^4+m_H^4\right)}{4m_W^2\left(m_H^24m_W^2\right)\left(d1\right)^2}}\right\}`$
$`+J_{00H}^{}\{4m_W^2{\displaystyle \frac{d2}{d1}}+{\displaystyle \frac{\left(d2\right)}{m_W^2\left(m_H^24m_W^2\right)\left(d1\right)^2}}[12dm_W^612m_W^6+m_H^6`$
$`5m_H^4m_W^2+4m_W^4m_H^2]\}`$
$`+J_{00W}{\displaystyle \frac{1}{12c^4d\left(2d7\right)\left(d4\right)\left(d1\right)^2}}[30240c^2d+61000c^2d^2+1680d2580d^2`$
$`53760c^4+249600c^4d361120c^4d^23340c^2d^5+18110c^2d^447720c^2d^3`$
$`744c^4d^6+11988c^4d^575452c^4d^4+233208c^4d^3+240c^2d^6345d^4+1440d^3`$
$`+30d^518432c^8d+10840c^6d^429472c^6d^31992c^6d^536928c^8d^3+43264c^8d^2`$
$`+14624c^8d^42720c^8d^522272c^6d+40352c^6d^2+144c^6d^6+192c^8d^6]`$
$`V_{000}{\displaystyle \frac{2\left(d2\right)\left(d^58d^4+27d^338d^224d+96\right)\left(c1\right)\left(c+1\right)}{3\left(d4\right)^2\left(d3\right)\left(d1\right)^2}}`$
$`+J_WJ_{00}\{{\displaystyle \frac{\left(d2\right)}{2d\left(d4\right)\left(d1\right)^2c^2}}[2c^2d^524c^2d^4+2c^4d^46c^4d^3+60c^2d^3+4c^4d^2`$
$`+34c^2d^2d^2136c^2d+4d+64c^2]{\displaystyle \frac{\left(d2\right)}{2m_H^2m_W^2\left(m_H^24m_W^2\right)\left(d1\right)^2}}[2d^3m_W^4m_H^2`$
$`8d^3m_W^6+32d^2m_W^66d^2m_W^4m_H^240dm_W^6+4dm_W^4m_H^2+m_H^6+16m_W^64m_H^4m_W^2]\}`$
$`J_ZJ_{00}\{{\displaystyle \frac{\left(d2\right)}{6d\left(d1\right)^2c^2}}[6c^4d^4+4c^6d^3+52c^4d^38c^2d^3+36c^6d^2122c^4d^2+19c^2d^2`$
$`+4d^27c^2d+92c^4d7d72c^6d+8c^2+32c^616c^4]{\displaystyle \frac{m_Z^2\left(d2\right)^2}{2m_H^2}}\}`$
$`+J_HJ_{00}\{(1{\displaystyle \frac{d}{2}})+{\displaystyle \frac{\left(d2\right)}{2m_H^2m_W^2\left(m_H^24m_W^2\right)\left(d1\right)^2}}[4d^2m_W^6+2d^2m_W^4m_H^2+12dm_{W}^{}{}_{}{}^{6}`$
$`6dm_W^4m_H^2+m_H^64m_H^4m_W^28m_W^6+4m_W^4m_H^2]\}`$
$`+J_{00}J_{00}{\displaystyle \frac{m_W^2\left(d2\right)}{6c^4\left(d4\right)\left(d1\right)^2}}[16c^4d^3+c^2d^38c^2d^22d^2+8c^6d^2+140c^4d^2+19c^2d`$
$`384c^4d+10d32c^6d12c^28+296c^4+24c^6]`$
$`+J_{00}J_{ZW}{\displaystyle \frac{m_W^2\left(d2\right)}{2c^4\left(d1\right)^2}}\left[4c^6d^2+8c^4d^252c^4d+5c^2d16c^6d+60c^4+113c^2+12c^6\right]`$
$`+J_{00}J_{HW}{\displaystyle \frac{\left(d2\right)}{2\left(m_H^24m_W^2\right)\left(d1\right)^2m_W^2}}[4d^2m_W^616dm_W^64dm_W^4m_H^2+m_H^4dm_W^2+m_H^6`$
$`9m_H^4m_W^2+12m_W^6+20m_W^4m_H^2],`$ (D.1)
$`(\mathrm{\Pi }_Z^{(2)}+\mathrm{\Pi }_Z^{(1)}\mathrm{\Pi }_Z^{(1)}{}_{}{}^{})=`$
$`V_{00W}\{{\displaystyle \frac{\left(d2\right)}{6c^2d\left(d1\right)^2}}[24c^4d^4+164c^4d^34c^2d^3+8c^6d^3340c^4d^210c^2d^2+72c^6d^2`$
$`+11d^2+232c^4d20d+34c^2d144c^6d+16c^2+64c^632c^4]+2{\displaystyle \frac{m_{Z}^{}{}_{}{}^{2}\left(d2\right)^2}{m_{H}^{}{}_{}{}^{2}}}\}`$
$`+V_{00Z}\{{\displaystyle \frac{1}{54}}{\displaystyle \frac{\left(d2\right)}{c^4d\left(d1\right)^2}}[976c^6d^3760c^2d^3+190d^3+352c^8d^3+1248c^4d^33872c^8d^2`$
$`+10736c^6d^213872c^4d^22180d^2+8540c^2d^211090c^2d+4928c^8d+17832c^4d`$
$`13664c^6d+2885d+3904c^64992c^4760+3040c^21408c^8]`$
$`+{\displaystyle \frac{\left(d2\right)\left(5+8c^410c^2\right)}{12m_{H}^{}{}_{}{}^{2}c^4m_{Z}^{}{}_{}{}^{2}\left(m_H^24m_Z^2\right)\left(d1\right)^2}}[16d^3m_{Z}^{}{}_{}{}^{6}4d^3m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}+12d^2m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}`$
$`64d^2m_{Z}^{}{}_{}{}^{6}+80dm_{Z}^{}{}_{}{}^{6}+m_{H}^{}{}_{}{}^{6}d4m_{H}^{}{}_{}{}^{4}dm_{Z}^{}{}_{}{}^{2}8dm_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}+16m_{H}^{}{}_{}{}^{4}m_{Z}^{}{}_{}{}^{2}`$
$`4m_{H}^{}{}_{}{}^{6}32m_{Z}^{}{}_{}{}^{6}]\}`$
$`F_{0000W}{\displaystyle \frac{m_Z^4\left(d2\right)\left(2c^2+2c^41\right)\left(2c^4c^2d+8c^2+2\right)}{3\left(d1\right)c^2}}`$
$`8F_{W0W00}{\displaystyle \frac{m_Z^4c^4\left(2+c^2\right)\left(d2\right)}{d1}}V_{H00Z}\{{\displaystyle \frac{m_Z^2\left(d2\right)^2\left(10c^2+5+8c^4\right)}{3c^4\left(d1\right)}}`$
$`+{\displaystyle \frac{\left(d2\right)\left(5+8c^410c^2\right)}{12m_{Z}^{}{}_{}{}^{2}c^4\left(m_H^24m_Z^2\right)\left(d1\right)^2}}[8d^2m_{Z}^{}{}_{}{}^{6}+m_{H}^{}{}_{}{}^{6}d+24dm_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}10m_{H}^{}{}_{}{}^{4}dm_{Z}^{}{}_{}{}^{2}`$
$`+32dm_{Z}^{}{}_{}{}^{6}4m_{H}^{}{}_{}{}^{6}+34m_{H}^{}{}_{}{}^{4}m_{Z}^{}{}_{}{}^{2}24m_{Z}^{}{}_{}{}^{6}72m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}]\}`$
$`+V_{W00W}{\displaystyle \frac{m_Z^2\left(d2\right)}{2\left(d1\right)^2c^2}}[44c^4d^2+24c^6d^28c^2d^280c^6d204c^4dd+42c^2d`$
$`+458c^2+56c^6+208c^4]`$
$`+J_{00W}{\displaystyle \frac{1}{6d\left(d1\right)^2c^2}}[1921008c^6d^228c^2d358c^2d^2+608c^4d+512d+256c^2`$
$`512c^4512c^688c^4d^2478d^2+1344c^6d+159d^3148c^4d^3+152c^6d^3+266c^2d^3`$
$`+24c^6d^452c^2d^4+44c^4d^416d^4]`$
$`+J_{00W}^{}{\displaystyle \frac{m_Z^22\left(c1\right)\left(c+1\right)}{3d\left(d1\right)^2c^2}}[32c^6+36c^6d^272c^6d+4c^6d^3+16c^448c^4d^2`$
$`+20c^4d^3+24c^4d+18c^2d15c^2d^28c^2+2c^2d^324d2d^3+8+18d^2]`$
$`+J_{00H}{\displaystyle \frac{\left(8c^4+510c^2\right)\left(d2\right)}{6c^4\left(d1\right)^2}}+J_{00H}^{}m_Z^24{\displaystyle \frac{\left(8c^4+510c^2\right)\left(d2\right)}{3\left(d1\right)c^4}}`$
$`+J_{00H}{\displaystyle \frac{\left(d2\right)\left(5+8c^410c^2\right)\left(3d8\right)\left(4dm_{Z}^{}{}_{}{}^{4}4m_{H}^{}{}_{}{}^{2}m_{Z}^{}{}_{}{}^{2}+m_{H}^{}{}_{}{}^{4}4m_{Z}^{}{}_{}{}^{4}\right)}{12c^4m_{Z}^{}{}_{}{}^{2}\left(m_H^24m_Z^2\right)\left(d1\right)^2}}`$
$`+J_{00H}^{}{\displaystyle \frac{\left(d2\right)\left(12dm_{Z}^{}{}_{}{}^{6}+4m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}+m_{H}^{}{}_{}{}^{6}5m_{H}^{}{}_{}{}^{4}m_{Z}^{}{}_{}{}^{2}12m_{Z}^{}{}_{}{}^{6}\right)\left(5+8c^410c^2\right)}{3c^4m_{Z}^{}{}_{}{}^{2}\left(m_H^24m_Z^2\right)\left(d1\right)^2}}`$
$`J_{00Z}{\displaystyle \frac{\left(7d^3+168d70d^280\right)\left(380c^2+95+624c^4488c^6+176c^8\right)}{27c^4d\left(d1\right)}}`$
$`V_{000}{\displaystyle \frac{8}{27}}{\displaystyle \frac{\left(d2\right)\left(d3\right)\left(d^24d+8\right)\left(c1\right)\left(c+1\right)\left(44c^461c^2+26\right)}{c^2\left(d4\right)^2\left(d1\right)}}`$
$`F_{0000Z}{\displaystyle \frac{m_Z^4}{54}}{\displaystyle \frac{\left(d2\right)\left(d12\right)\left(380c^2+95+624c^4488c^6+176c^8\right)}{\left(d1\right)c^4}}`$
$`+J_WJ_{00}\{{\displaystyle \frac{\left(d2\right)}{3\left(d1\right)^2c^4d}}[16c^8d^420c^6d^4+10c^4d^434c^2d^3+176c^6d^332c^4d^3`$
$`+2d^3160c^8d^3+368c^8d^223d^2114c^4d^2+172c^2d^2340c^6d^2+120c^6d`$
$`+268c^4d+34d204c^2d224c^8d832c^4+16c^2]{\displaystyle \frac{m_{Z}^{}{}_{}{}^{2}\left(d2\right)^2\left(5+8c^410c^2\right)}{3c^2m_{H}^{}{}_{}{}^{2}}}\}`$
$`+J_ZJ_{00}\{{\displaystyle \frac{\left(d2\right)}{54\left(d1\right)^2c^4d}}[760c^2d^3976c^6d^3+352c^8d^3+1248c^4d^3+190d^3`$
$`13800c^4d^22135d^23872c^8d^2+10736c^6d^2+8450c^2d^213664c^6d+2795d`$
$`+17688c^4d10910c^2d+4928c^8d7601408c^8+3040c^24992c^4+3904c^6]`$
$`{\displaystyle \frac{\left(d2\right)\left(5+8c^410c^2\right)}{6c^4m_{H}^{}{}_{}{}^{2}\left(m_H^24m_Z^2\right)\left(d1\right)^2m_{Z}^{}{}_{}{}^{2}}}[d^3m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}4d^3m_{Z}^{}{}_{}{}^{6}2d^2m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}+16d^2m_{Z}^{}{}_{}{}^{6}`$
$`20dm_{Z}^{}{}_{}{}^{6}dm_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}+2m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}4m_{H}^{}{}_{}{}^{4}m_{Z}^{}{}_{}{}^{2}+m_{H}^{}{}_{}{}^{6}+8m_{Z}^{}{}_{}{}^{6}]\}`$
$`J_HJ_{00}\{{\displaystyle \frac{\left(8c^4+510c^2\right)\left(d2\right)}{6c^4}}{\displaystyle \frac{\left(d2\right)\left(5+8c^410c^2\right)}{6c^4m_{H}^{}{}_{}{}^{2}\left(m_H^24m_Z^2\right)\left(d1\right)^2m_{Z}^{}{}_{}{}^{2}}}[2d^2m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}`$
$`4d^2m_{Z}^{}{}_{}{}^{6}+12dm_{Z}^{}{}_{}{}^{6}6dm_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}+m_{H}^{}{}_{}{}^{6}4m_{H}^{}{}_{}{}^{4}m_{Z}^{}{}_{}{}^{2}+4m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}8m_{Z}^{}{}_{}{}^{6}]\}`$
$`J_{00}J_{00}{\displaystyle \frac{m_Z^2}{54}}{\displaystyle \frac{\left(d2\right)}{c^4\left(d4\right)\left(d1\right)^2}}[1920c^4d^31816c^6d^3+768c^8d^3+205d^3942c^2d^3`$
$`14244c^4d^21070d^26696c^8d^2+14840c^6d^2+5928c^2d^243120c^6d+1345d`$
$`+36972c^4d12074c^2d+20360c^8d138021344c^8+12488c^2+46224c^638688c^4]`$
$`+J_{00}J_{WW}{\displaystyle \frac{m_Z^2\left(d2\right)}{6\left(d1\right)^2c^4}}[40c^2d^2144c^6d^2+16c^8d^2+140c^4d^2+64c^{10}d^2+920c^6d`$
$`808c^4d5d+200c^2d+16c^8d512c^{10}d+10230c^2+448c^{10}+224c^8`$
$`1240c^6+968c^4]`$
$`+J_{00}J_{HZ}{\displaystyle \frac{\left(d2\right)\left(5+8c^410c^2\right)}{6c^4\left(m_H^24m_Z^2\right)\left(d1\right)^2m_{Z}^{}{}_{}{}^{2}}}[4d^2m_{Z}^{}{}_{}{}^{6}+m_{H}^{}{}_{}{}^{4}dm_{Z}^{}{}_{}{}^{2}16dm_{Z}^{}{}_{}{}^{6}4dm_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}`$
$`+m_{H}^{}{}_{}{}^{6}+12m_{Z}^{}{}_{}{}^{6}+20m_{Z}^{}{}_{}{}^{4}m_{H}^{}{}_{}{}^{2}9m_{H}^{}{}_{}{}^{4}m_{Z}^{}{}_{}{}^{2}],`$ (D.2)
where

$`F_{ABCDE}`$ $`=`$ $`{\displaystyle \frac{d^dk_1}{\mathrm{\Gamma }(1+\epsilon )}\frac{d^dk_2}{\mathrm{\Gamma }(1+\epsilon )}P^{(1)}(k_1,m_A)P^{(1)}(k_2,m_B)}`$
$`P^{(1)}(k_1p,m_C)P^{(1)}(k_2p,m_D)P^{(1)}(k_1k_2,m_E),`$
$`V_{ABCD}`$ $`=`$ $`{\displaystyle \frac{d^dk_1}{\mathrm{\Gamma }(1+\epsilon )}\frac{d^dk_2}{\mathrm{\Gamma }(1+\epsilon )}P^{(1)}(k_1p,m_A)P^{(1)}(k_1k_2,m_B)}`$
$`P^{(1)}(k_2,m_C)P^{(1)}(k_1,m_D),`$
$`J_{00A}`$ $`=`$ $`{\displaystyle \frac{d^dk_1}{\mathrm{\Gamma }(1+\epsilon )}\frac{d^dk_2}{\mathrm{\Gamma }(1+\epsilon )}P^{(1)}(k_1p,0)P^{(1)}(k_1k_2,0)P^{(1)}(k_2,m_A)},`$
$`J_{00A}^{}`$ $`=`$ $`{\displaystyle \frac{d^nk_1}{\mathrm{\Gamma }(1+\epsilon )}\frac{d^nk_2}{\mathrm{\Gamma }(1+\epsilon )}P^{(1)}(k_1p,0)P^{(1)}(k_1k_2,0)P^{(2)}(k_2,m_A)},`$
$`V_{ABC}`$ $`=`$ $`{\displaystyle \frac{d^dk_1}{\mathrm{\Gamma }(1+\epsilon )}\frac{d^dk_2}{\mathrm{\Gamma }(1+\epsilon )}P^{(1)}(k_1,m_A)P^{(1)}(k_1k_2,m_B)P^{(1)}(k_2,m_C)},`$
$`J_{AB}`$ $`=`$ $`{\displaystyle \frac{d^dk_1}{\mathrm{\Gamma }(1+\epsilon )}P^{(1)}(k_1p,m_A)P^{(1)}(k_1,m_B)},`$
$`J_A`$ $`=`$ $`{\displaystyle \frac{d^dk_1}{\mathrm{\Gamma }(1+\epsilon )}P^{(1)}(k_1,m_A)},`$ (D.3)
with $`c=m_W/m_Z=\mathrm{cos}\mathrm{\Theta }_W`$ and $`P^{(\sigma )}(p,m)=1/(p^2+m^2)^\sigma `$ and set $`\{A,B,C,D,E\}`$ denotes the particles. The external momentum belongs to the mass shell.

## Appendix E $`\overline{\mathrm{MS}}`$ vs. pole masses at two-loop order

After expansion of the diagrams with respect to $`\mathrm{sin}^2\theta _W`$ we get rid of one of the boson masses and write the functions $`X_V^{(2)}`$ introduced in (5.1) in the form

$$X_V^{(2)}=\frac{m_H^4}{m_W^4}\underset{k=0}{\overset{4}{}}\mathrm{sin}^{2k}\theta _WA_k^V.$$
(E.1)
In particular for the $`Z`$ boson propagator we eliminate $`m_W`$ and vice versa. Consequently, the coefficients $`A_i^V`$ in the above formula are functions of the Higgs and top masses and one of the boson masses. We expand this function with respect to $`m_V^2/m_{H,t}^2`$

$$A_i^V=\underset{j=0}{\overset{5}{}}B_{i,j}^V(m_H^2,m_t^2)(m_V^2)^j$$

and calculate analytically the first sixth coefficients 262626 In parameterization (E.1) $`B_{j,0}=0`$ for $`j>0`$.. This is not a naive Taylor expansion. The general rules for asymptotic expansions  allow us to extract also logarithmic dependences, or in other words, to preserve all analytical properties of the original diagrams. In the result of the asymptotic expansion all propagator diagrams are reduced to single scale massive diagrams (including the two-loop bubbles). In the finite part, we meet the following constants and functions:

$`S_0`$ $`=`$ $`{\displaystyle \frac{\pi }{\sqrt{3}}}1.813799365\mathrm{},S_1={\displaystyle \frac{\pi }{\sqrt{3}}}\mathrm{ln}31.992662272\mathrm{},`$
$`S_2`$ $`=`$ $`{\displaystyle \frac{4}{9}}{\displaystyle \frac{\text{Cl}_2\left(\frac{\pi }{3}\right)}{\sqrt{3}}}0.260434137632162\mathrm{},S_3=\pi \text{Cl}_2\left(\frac{\pi }{3}\right)3.188533097\mathrm{}`$ (E.2)
Furthermore, $`\mathrm{ln}(m_V^2)`$ denotes $`\mathrm{ln}\left(m_V^2/\mu ^2\right)`$ where $`\mu `$ is the ’t Hooft scale. We also introduce the notation

$`z_H={\displaystyle \frac{m_Z^2}{m_H^2}},z_t={\displaystyle \frac{m_Z^2}{m_t^2}},\omega _H={\displaystyle \frac{m_W^2}{m_H^2}},\omega _t={\displaystyle \frac{m_W^2}{m_t^2}},`$
$`t_H={\displaystyle \frac{m_t^2}{m_H^2}},y={\displaystyle \frac{m_t^2}{m_H^2m_t^2}},r={\displaystyle \frac{m_t^2}{m_H^24m_t^2}},`$
$`\text{L}_2=\text{Li}_2\left({\displaystyle \frac{m_H^2}{m_t^2}}\right),\text{L}_1=\mathrm{ln}\left(1{\displaystyle \frac{m_H^2}{m_t^2}}\right),F=\mathrm{\Phi }\left({\displaystyle \frac{m_H^2}{4m_t^2}}\right),`$ (E.3)
where $`N_m`$ is a number of massless fermion families. and $`\mathrm{\Phi }(z)`$ is the finite part of two-loop bubble master integrals defined in . We rewrite is as follow (see -):

$`\mathrm{\Phi }(z)=2{\displaystyle \frac{1\eta }{1+\eta }}\left[\text{Li}_2\left(\eta \right)\text{Li}_2\left({\displaystyle \frac{1}{\eta }}\right)\right],\eta ={\displaystyle \frac{1\sqrt{\frac{z}{z1}}}{1+\sqrt{\frac{z}{z1}}}}.`$ (E.4)
Below, we present the coefficients for $`\mathrm{Re}A_{i,fermion}^V`$ including fermionic contribution only, which are defined as

$$A_{i,full}^V=A_{i,boson}^V+A_{i,fermion}^V$$

and $`A_{i,boson}^V`$ given in Appendix D of . In the following we present the coefficients in the form $`A_{i,j}^V=B_{i,j}^V(m_H^2,m_t^2)(m_V^2)^j`$.

### E.1 The expansion coefficients for the $`W`$

$`A_{0,0}^W={\displaystyle \frac{9}{8}}\mathrm{ln}(\omega _H)\mathrm{ln}(m_W^2)t_H{\displaystyle \frac{9}{8}}\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_H+3/8\pi ^2t_H{\displaystyle \frac{11}{16}}\pi ^2t_{H}^{}{}_{}{}^{2}+1/2\pi ^2t_{H}^{}{}_{}{}^{3}`$
$`{\displaystyle \frac{9}{16}}\text{L}_1\mathrm{ln}(\omega _H)+{\displaystyle \frac{9}{16}}\text{L}_1\mathrm{ln}(\omega _t){\displaystyle \frac{57}{16}}\mathrm{ln}(\omega _H)t_H+{\displaystyle \frac{57}{8}}\mathrm{ln}(\omega _H)t_{H}^{}{}_{}{}^{2}+{\displaystyle \frac{9}{16}}\mathrm{ln}(\omega _t)t_H`$
$`{\displaystyle \frac{57}{4}}\mathrm{ln}(\omega _t)t_{H}^{}{}_{}{}^{2}+15\mathrm{ln}(\omega _t)t_{H}^{}{}_{}{}^{3}+36\mathrm{ln}(\omega _t)t_{H}^{}{}_{}{}^{4}+{\displaystyle \frac{9}{8}}\left(\mathrm{ln}(\omega _t)\right)^2t_H{\displaystyle \frac{117}{16}}\left(\mathrm{ln}(\omega _t)\right)^2t_{H}^{}{}_{}{}^{2}`$
$`9/2\left(\mathrm{ln}(\omega _t)\right)^2t_{H}^{}{}_{}{}^{3}+36\left(\mathrm{ln}(\omega _t)\right)^2t_{H}^{}{}_{}{}^{4}+3\mathrm{ln}(m_W^2)t_H+{\displaystyle \frac{57}{8}}\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{2}15\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{3}`$
$`36\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{4}{\displaystyle \frac{117}{16}}\left(\mathrm{ln}(m_W^2)\right)^2t_{H}^{}{}_{}{}^{2}9/2\left(\mathrm{ln}(m_W^2)\right)^2t_{H}^{}{}_{}{}^{3}+36\left(\mathrm{ln}(m_W^2)\right)^2t_{H}^{}{}_{}{}^{4}`$
$`+{\displaystyle \frac{9}{16}}Ft_H+{\displaystyle \frac{45}{8}}Ft_{H}^{}{}_{}{}^{2}{\displaystyle \frac{27}{2}}Ft_{H}^{}{}_{}{}^{3}9/4\text{L}_2t_H+{\displaystyle \frac{45}{16}}\text{L}_2t_{H}^{}{}_{}{}^{2}{\displaystyle \frac{9}{8}}\text{L}_2t_{H}^{}{}_{}{}^{3}9/4\text{L}_1\mathrm{ln}(\omega _t)t_H`$
$`{\displaystyle \frac{81}{16}}t_H+{\displaystyle \frac{9}{16}}\text{L}_2+9/4\text{L}_1\mathrm{ln}(\omega _H)t_H{\displaystyle \frac{45}{16}}\text{L}_1\mathrm{ln}(\omega _H)t_{H}^{}{}_{}{}^{2}+{\displaystyle \frac{9}{8}}\text{L}_1\mathrm{ln}(\omega _H)t_{H}^{}{}_{}{}^{3}`$
$`+{\displaystyle \frac{45}{16}}\text{L}_1\mathrm{ln}(\omega _t)t_{H}^{}{}_{}{}^{2}{\displaystyle \frac{9}{8}}\text{L}_1\mathrm{ln}(\omega _t)t_{H}^{}{}_{}{}^{3}{\displaystyle \frac{9}{8}}\mathrm{ln}(\omega _h)\mathrm{ln}(\omega _t)t_H+{\displaystyle \frac{117}{8}}\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{2}`$
$`+9\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{3}72\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{4}+27t_{H}^{}{}_{}{}^{3}{\displaystyle \frac{3}{32}}\pi ^2+{\displaystyle \frac{213}{64}}t_{H}^{}{}_{}{}^{2}{\displaystyle \frac{9}{32}}F`$ (E.1)
$`A_{0,1}^W=1/4\text{L}_2\omega _Ht_{H}^{}{}_{}{}^{2}+{\displaystyle \frac{177}{2}}\mathrm{ln}(m_W^2)\omega _Ht_{H}^{}{}_{}{}^{2}{\displaystyle \frac{339}{16}}\mathrm{ln}(m_W^2)\omega _ht_H5F\omega _Ht_{H}^{}{}_{}{}^{3}`$
$`+{\displaystyle \frac{925}{48}}\mathrm{ln}(\omega _t)\omega _Ht_H3/4\mathrm{ln}(\omega _H)\mathrm{ln}(\omega _t)\omega _H2/3\mathrm{ln}(\omega _H)\omega _HN_m+3/4\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)\omega _H`$
$`+{\displaystyle \frac{113}{16}}\left(\mathrm{ln}(m_W^2)\right)^2\omega _ht_H+16/3\omega _Ht_{H}^{}{}_{}{}^{2}N_m1/2\text{L}_1\mathrm{ln}(\omega _H)\omega _Ht_H1/4\text{L}_1\mathrm{ln}(\omega _H)\omega _Ht_{H}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{37}{24}}\text{L}_1\mathrm{ln}(\omega _H)\omega _Ht_{H}^{}{}_{}{}^{3}+1/2\text{L}_1\mathrm{ln}(\omega _t)\omega _Ht_H+1/4\text{L}_1\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{2}{\displaystyle \frac{37}{24}}\text{L}_1\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{3}`$
$`+3/8\text{L}_1\mathrm{ln}(\omega _h)\omega _t{\displaystyle \frac{79}{48}}\mathrm{ln}(\omega _H)\omega _Ht_H+{\displaystyle \frac{11}{12}}\text{L}_1\mathrm{ln}(\omega _t)\omega _H{\displaystyle \frac{11}{72}}\pi ^2\omega _H+1/16\pi ^2\omega _t+{\displaystyle \frac{977}{192}}\omega _ht_H`$
$`{\displaystyle \frac{1451}{24}}\omega _Ht_{H}^{}{}_{}{}^{2}{\displaystyle \frac{7}{9}}\omega _HN_m{\displaystyle \frac{1621}{24}}\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{2}+{\displaystyle \frac{11}{6}}\mathrm{ln}(m_W^2)\omega _HN_m+6\left(\mathrm{ln}(\omega _t)\right)^2\omega _Ht_{H}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{33}{16}}\omega _Ht_HS_0{\displaystyle \frac{11}{12}}\text{L}_1\mathrm{ln}(\omega _h)\omega _H+{\displaystyle \frac{41}{12}}F\omega _Ht_{H}^{}{}_{}{}^{2}+{\displaystyle \frac{23}{36}}\omega _H{\displaystyle \frac{43}{2}}\left(\mathrm{ln}(m_W^2)\right)^2\omega _Ht_{H}^{}{}_{}{}^{2}`$
$`+\mathrm{ln}(\omega _H)\mathrm{ln}(m_W^2)\omega _H\left(\mathrm{ln}(m_W^2)\right)^2\omega _H{\displaystyle \frac{7}{48}}\pi ^2\omega _ht_H2/3\omega _Ht_HN_m+{\displaystyle \frac{99}{2}}\omega _ht_{H}^{}{}_{}{}^{2}S_0`$
$`3/8\text{L}_1\mathrm{ln}(\omega _t)\omega _t+{\displaystyle \frac{277}{24}}\mathrm{ln}(\omega _H)\omega _Ht_{H}^{}{}_{}{}^{2}+3/16\left(\mathrm{ln}(\omega _t)\right)^2\omega _Ht_H{\displaystyle \frac{5}{24}}F\omega _ht_H+1/2\text{L}_2\omega _Ht_H`$
$`{\displaystyle \frac{37}{24}}\text{L}_2\omega _Ht_{H}^{}{}_{}{}^{3}+\mathrm{ln}(\omega _H)\mathrm{ln}(m_W^2)\omega _HN_m\left(\mathrm{ln}(m_W^2)\right)^2\omega _HN_m1/24\mathrm{ln}(\omega _H)\omega _H`$
$`1/2\mathrm{ln}(\omega _t)\omega _H+5/6\mathrm{ln}(m_W^2)\omega _H5/6F\omega _H+3/16F\omega _t+{\displaystyle \frac{11}{12}}\text{L}_2\omega _H3/8\text{L}_2\omega _t`$
$`5/4\mathrm{ln}(\omega _H)\mathrm{ln}(\omega _t)\omega _ht_H+5\mathrm{ln}(\omega _H)\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{2}+5/4\mathrm{ln}(\omega _H)\mathrm{ln}(m_W^2)\omega _Ht_H`$
$`5\mathrm{ln}(\omega _h)\mathrm{ln}(m_W^2)\omega _Ht_{H}^{}{}_{}{}^{2}{\displaystyle \frac{99}{8}}\mathrm{ln}(\omega _t)\omega _Ht_HS_04/3\mathrm{ln}(\omega _t)\omega _Ht_HN_m+{\displaystyle \frac{99}{2}}\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{2}S_0`$
$`+16/3\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{2}N_m{\displaystyle \frac{29}{4}}\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)\omega _Ht_H+{\displaystyle \frac{31}{2}}\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)\omega _Ht_{H}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{99}{8}}\mathrm{ln}(m_W^2)\omega _Ht_HS_0+7/3\mathrm{ln}(m_W^2)\omega _Ht_HN_m{\displaystyle \frac{99}{2}}\mathrm{ln}(m_W^2)\omega _Ht_{H}^{}{}_{}{}^{2}S_0`$
$`{\displaystyle \frac{40}{3}}\mathrm{ln}(m_W^2)\omega _ht_{H}^{}{}_{}{}^{2}N_m2\left(\mathrm{ln}(m_W^2)\right)^2\omega _ht_HN_m+8\left(\mathrm{ln}(m_W^2)\right)^2\omega _ht_{H}^{}{}_{}{}^{2}N_m`$
$`+2\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)\omega _Ht_HN_m8\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)\omega _Ht_{H}^{}{}_{}{}^{2}N_m`$ (E.2)
$`A_{0,2}^W={\displaystyle \frac{1}{10368}}\omega _H(12960\text{L}_2\omega _Ht_{H}^{}{}_{}{}^{2}95904\mathrm{ln}(m_W^2)\omega _Ht_{H}^{}{}_{}{}^{2}98172\mathrm{ln}(m_W^2)\omega _Ht_H`$
$`+139968F\omega _Ht_{H}^{}{}_{}{}^{3}+162\omega _Hy5292\mathrm{ln}(\omega _t)\omega _Ht_H17280\mathrm{ln}(\omega _h)\mathrm{ln}(\omega _t)\omega _H`$
$`26880\mathrm{ln}(\omega _H)\omega _HN_m97632\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)\omega _H139968\left(\mathrm{ln}(m_W^2)\right)^2\omega _Ht_H`$
$`6480\text{L}_1\mathrm{ln}(\omega _H)\omega _Ht_H+12960\text{L}_1\mathrm{ln}(\omega _H)\omega _Ht_{H}^{}{}_{}{}^{2}28080\text{L}_1\mathrm{ln}(\omega _H)\omega _Ht_{H}^{}{}_{}{}^{3}`$
$`+6480\text{L}_1\mathrm{ln}(\omega _t)\omega _Ht_H12960\text{L}_1\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{2}+28080\text{L}_1\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{3}`$
$`+1296\text{L}_1\mathrm{ln}(\omega _H)\omega _t+39312\mathrm{ln}(\omega _h)\omega _Ht_H+4320\text{L}_1\mathrm{ln}(\omega _t)\omega _H12024\pi ^2\omega _H`$
$`+216\pi ^2\omega _t+276426\omega _Ht_H516672\omega _Ht_{H}^{}{}_{}{}^{2}+1144080\omega _HN_m866052\omega _HS_2`$
$`15984\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{2}666816\mathrm{ln}(m_W^2)\omega _HN_m4320\text{L}_1\mathrm{ln}(\omega _h)\omega _H69552F\omega _Ht_{H}^{}{}_{}{}^{2}`$
$`+886115\omega _H+3888\omega _t+559872\left(\mathrm{ln}(m_W^2)\right)^2\omega _Ht_{H}^{}{}_{}{}^{2}116424\omega _HS_0`$
$`+23040\mathrm{ln}(\omega _H)\mathrm{ln}(m_W^2)\omega _H+133056\left(\mathrm{ln}(m_W^2)\right)^2\omega _H+15552\pi ^2\omega _Ht_H1296\text{L}_1\mathrm{ln}(\omega _t)\omega _t`$
$`214704\mathrm{ln}(\omega _H)\omega _Ht_{H}^{}{}_{}{}^{2}+69984\left(\mathrm{ln}(\omega _t)\right)^2\omega _Ht_H+27648F\omega _Ht_H+6480\text{L}_2\omega _Ht_H`$
$`+28080\text{L}_2\omega _Ht_{H}^{}{}_{}{}^{3}+23040\mathrm{ln}(\omega _H)\mathrm{ln}(m_W^2)\omega _HN_m+114624\left(\mathrm{ln}(m_W^2)\right)^2\omega _HN_m`$
$`20742\mathrm{ln}(\omega _H)\omega _H+287154\mathrm{ln}(\omega _t)\omega _H608340\mathrm{ln}(m_W^2)\omega _H7344F\omega _H+648F\omega _t`$
$`+4320\text{L}_2\omega _h1296\text{L}_2\omega _t+34992\mathrm{ln}(\omega _H)\mathrm{ln}(\omega _t)\omega _ht_H93312\mathrm{ln}(\omega _H)\mathrm{ln}(\omega _t)\omega _Ht_{H}^{}{}_{}{}^{2}`$
$`34992\mathrm{ln}(\omega _H)\mathrm{ln}(m_W^2)\omega _Ht_H+93312\mathrm{ln}(\omega _H)\mathrm{ln}(m_W^2)\omega _Ht_{H}^{}{}_{}{}^{2}+69984\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)\omega _Ht_H`$
$`559872\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)\omega _Ht_{H}^{}{}_{}{}^{2}55296\omega _HS_3+43008\omega _hN_{m}^{}{}_{}{}^{2}\mathrm{ln}(m_W^2)`$
$`18432\omega _HN_{m}^{}{}_{}{}^{2}\left(\mathrm{ln}(m_W^2)\right)^26156\omega _H\mathrm{ln}(\omega _t)y+103680\omega _H\zeta (3)1296\omega _H\left(\mathrm{ln}(\omega _H)\right)^2`$
$`+23760\omega _H\left(\mathrm{ln}(\omega _t)\right)^220480\omega _HN_{m}^{}{}_{}{}^{2}+2916\mathrm{ln}(\omega _H)\omega _t2916\mathrm{ln}(m_W^2)\omega _t`$
$`+228096\omega _HN_mS_0\mathrm{ln}(m_W^2)+27648\omega _H\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)N_m99792\omega _h\mathrm{ln}(\omega _t)S_0`$
$`32256\omega _H\mathrm{ln}(\omega _t)N_m42624\omega _HN_m\pi ^25184\omega _H\left(\mathrm{ln}(\omega _h)\right)^2N_m+6156\omega _H\mathrm{ln}(\omega _H)y`$
$`+162\omega _h\mathrm{ln}(\omega _H)y^2+18432\omega _HN_{m}^{}{}_{}{}^{2}\pi ^2+414720\omega _HN_m\zeta (3)162\omega _H\mathrm{ln}(\omega _t)y^2`$
$`+228096\omega _HS_0\mathrm{ln}(m_W^2)38016\omega _HN_mS_0221184\omega _HN_mS_31539648\omega _hN_mS_2)`$ (E.3)
$`A_{1,1}^W={\displaystyle \frac{1}{576}}t_H\omega _H(476\pi ^2972S_0+82599504t_HS_03552t_H+2376\mathrm{ln}(\omega _t)S_0`$
$`+336\mathrm{ln}(\omega _t)9504\mathrm{ln}(\omega _t)t_HS_03552\mathrm{ln}(\omega _t)t_H+612\left(\mathrm{ln}(\omega _t)\right)^22304\left(\mathrm{ln}(\omega _t)\right)^2t_H`$
$`576\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)+2016\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_H2376\mathrm{ln}(m_W^2)S_0660\mathrm{ln}(m_W^2)`$
$`+9504\mathrm{ln}(m_W^2)t_HS_0+960\mathrm{ln}(m_W^2)t_H36\left(\mathrm{ln}(m_W^2)\right)^2+288\left(\mathrm{ln}(m_W^2)\right)^2t_H)`$ (E.4)
$`A_{1,2}^W={\displaystyle \frac{1}{1296}}\omega _{H}^{}{}_{}{}^{2}(1128\pi ^2N_m+576\mathrm{ln}(m_W^2)N_m1296\left(\mathrm{ln}(m_W^2)\right)^2N_m+10368\zeta (3)N_m`$
$`+5832\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_H24300\mathrm{ln}(\omega _t)t_H31104\mathrm{ln}(\omega _t)t_{H}^{}{}_{}{}^{2}+5832\left(\mathrm{ln}(\omega _t)\right)^2t_H`$
$`+15552\mathrm{ln}(m_W^2)t_H+7776\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{2}+46656\left(\mathrm{ln}(m_W^2)\right)^2t_{H}^{}{}_{}{}^{2}+32238t_H+5364S_0`$
$`+1452N_m46656\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{2}+3969\mathrm{ln}(\omega _t)+1116\mathrm{ln}(m_W^2)899\pi ^227216t_{H}^{}{}_{}{}^{2}`$
$`+1458S_2+2304S_3+1728\zeta (3)342\left(\mathrm{ln}(\omega _t)\right)^21296\left(\mathrm{ln}(m_W^2)\right)^211664\left(\mathrm{ln}(m_W^2)\right)^2t_H`$
$`+972\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)6336S_0N_m17496S_2N_m+9216S_3N_m8029`$
$`4752\mathrm{ln}(m_W^2)S_04752\mathrm{ln}(m_W^2)S_0N_m+2700\mathrm{ln}(\omega _t)S_0)`$ (E.5)
$`A_{2,1}^W={\displaystyle \frac{1}{576}}t_H\omega _H(476\pi ^2+1872S_08619+8640t_HS_06528t_H2160\mathrm{ln}(\omega _t)S_0`$
$`+2184\mathrm{ln}(\omega _t)+8640\mathrm{ln}(\omega _t)t_HS_06528\mathrm{ln}(\omega _t)t_H612\left(\mathrm{ln}(\omega _t)\right)^2+2304\left(\mathrm{ln}(\omega _t)\right)^2t_H`$
$`+576\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)2016\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_H+2160\mathrm{ln}(m_W^2)S_01860\mathrm{ln}(m_W^2)`$
$`8640\mathrm{ln}(m_W^2)t_HS_0xsy+9120\mathrm{ln}(m_W^2)t_H+36\left(\mathrm{ln}(m_W^2)\right)^2288\left(\mathrm{ln}(m_W^2)\right)^2t_H)`$ (E.6)
$`A_{2,2}^W={\displaystyle \frac{1}{2592}}\omega _{H}^{}{}_{}{}^{2}(17496\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_H+3564\mathrm{ln}(\omega _t)t_H+139968\mathrm{ln}(\omega _t)t_{H}^{}{}_{}{}^{2}`$
$`17496\left(\mathrm{ln}(\omega _t)\right)^2t_H+22680\mathrm{ln}(m_W^2)t_H69984\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{2}139968\left(\mathrm{ln}(m_W^2)\right)^2t_{H}^{}{}_{}{}^{2}`$
$`70686t_H63114S_013708N_m+139968\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)t_{H}^{}{}_{}{}^{2}9561\mathrm{ln}(\omega _t)`$
$`10632\mathrm{ln}(m_W^2)+4991\pi ^2+104976t_{H}^{}{}_{}{}^{2}113238S_2+3888\zeta (3)+90\left(\mathrm{ln}(\omega _t)\right)^2`$
$`+4752\left(\mathrm{ln}(m_W^2)\right)^2+34992\left(\mathrm{ln}(m_W^2)\right)^2t_H7488\mathrm{ln}(\omega _t)S_01944\mathrm{ln}(\omega _t)\mathrm{ln}(m_W^2)`$
$`+25344\mathrm{ln}(m_W^2)S_0103152S_0N_m141912S_2N_m+4752\left(\mathrm{ln}(m_W^2)\right)^2N_m+5768\pi ^2N_m`$
$`+25344\mathrm{ln}(m_W^2)S_0N_m+1728\zeta (3)N_m9552\mathrm{ln}(m_W^2)N_m228)`$ (E.7)
### E.2 The expansion coefficients for the $`Z`$

$`A_{0,0}^Z={\displaystyle \frac{9}{8}}\mathrm{ln}(z_t)\mathrm{ln}(z_h)t_h+{\displaystyle \frac{9}{8}}\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)t_h{\displaystyle \frac{9}{8}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_h{\displaystyle \frac{5}{16}}\pi ^2t_{h}^{}{}_{}{}^{2}`$
$`+1/2\pi ^2t_{h}^{}{}_{}{}^{3}{\displaystyle \frac{75}{16}}\mathrm{ln}(z_h)t_h+6\mathrm{ln}(z_h)t_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{57}{16}}\mathrm{ln}(m_Z^2)t_h+{\displaystyle \frac{105}{16}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{39}{2}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{3}36\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{4}{\displaystyle \frac{117}{16}}\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}9/2\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{3}`$
$`+36\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{4}+{\displaystyle \frac{9}{8}}\mathrm{ln}(z_t)t_h{\displaystyle \frac{201}{16}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{39}{2}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{3}+36\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{4}`$
$`+{\displaystyle \frac{9}{8}}\left(\mathrm{ln}(z_t)\right)^2t_h{\displaystyle \frac{117}{16}}\left(\mathrm{ln}(z_t)\right)^2t_{h}^{}{}_{}{}^{2}9/2\left(\mathrm{ln}(z_t)\right)^2t_{h}^{}{}_{}{}^{3}+36\left(\mathrm{ln}(z_t)\right)^2t_{h}^{}{}_{}{}^{4}{\displaystyle \frac{9}{8}}Ft_h`$
$`+{\displaystyle \frac{135}{16}}Ft_{h}^{}{}_{}{}^{2}{\displaystyle \frac{27}{2}}Ft_{h}^{}{}_{}{}^{3}{\displaystyle \frac{81}{16}}t_h+{\displaystyle \frac{117}{8}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}`$
$`+9\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{3}72\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{4}+{\displaystyle \frac{9}{64}}t_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{63}{2}}t_{h}^{}{}_{}{}^{3}`$ (E.8)
$`A_{0,1}^Z=3/8\mathrm{ln}(z_t)\mathrm{ln}(z_h)z_h+{\displaystyle \frac{5}{48}}t_hz_hr1/2Ft_hz_h{\displaystyle \frac{49}{48}}\mathrm{ln}(z_h)t_hz_h+10\mathrm{ln}(z_h)t_{h}^{}{}_{}{}^{2}z_h`$
$`+5/8\mathrm{ln}(z_h)t_hz_hr^2+5/4\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)t_hz_h5\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_h+{\displaystyle \frac{99}{8}}\mathrm{ln}(m_Z^2)t_hz_hS_0`$
$`+4/3\mathrm{ln}(m_Z^2)t_hz_hN_m+{\displaystyle \frac{579}{64}}t_hz_h{\displaystyle \frac{121}{2}}t_{h}^{}{}_{}{}^{2}z_h{\displaystyle \frac{7}{9}}z_hN_m+{\displaystyle \frac{7}{48}}\mathrm{ln}(z_h)z_h+{\displaystyle \frac{49}{48}}\mathrm{ln}(m_Z^2)z_h`$
$`{\displaystyle \frac{7}{16}}\mathrm{ln}(z_t)z_h+{\displaystyle \frac{409}{128}}Ft_{h}^{}{}_{}{}^{2}z_h+{\displaystyle \frac{31}{288}}z_h+\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_h\left(\mathrm{ln}(m_Z^2)\right)^2z_h+3/8\pi ^2t_hz_h`$
$`{\displaystyle \frac{33}{4}}t_hz_hS_0+{\displaystyle \frac{99}{2}}t_{h}^{}{}_{}{}^{2}z_hS_0+\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_hN_m\left(\mathrm{ln}(m_Z^2)\right)^2z_hN_m+16/3t_{h}^{}{}_{}{}^{2}z_hN_m`$
$`{\displaystyle \frac{55}{48}}\mathrm{ln}(z_h)t_hz_hr2/3\mathrm{ln}(z_h)z_hN_m{\displaystyle \frac{281}{16}}\mathrm{ln}(m_Z^2)t_hz_h+87\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_h`$
$`+{\displaystyle \frac{11}{6}}\mathrm{ln}(m_Z^2)z_hN_m+{\displaystyle \frac{113}{16}}\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_h{\displaystyle \frac{43}{2}}\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_h+{\displaystyle \frac{205}{12}}\mathrm{ln}(z_t)t_hz_h`$
$`{\displaystyle \frac{135}{2}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_h+3/8\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_h+{\displaystyle \frac{15}{16}}\left(\mathrm{ln}(z_t)\right)^2t_hz_h+3\left(\mathrm{ln}(z_t)\right)^2t_{h}^{}{}_{}{}^{2}z_h`$
$`5Ft_{h}^{}{}_{}{}^{3}z_h{\displaystyle \frac{99}{2}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hS_0{\displaystyle \frac{40}{3}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hN_m2\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_hN_m`$
$`+8\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_hN_m+{\displaystyle \frac{55}{48}}\mathrm{ln}(z_t)t_hz_hr5/8\mathrm{ln}(z_t)t_hz_hr^2{\displaystyle \frac{99}{8}}\mathrm{ln}(z_t)t_hz_hS_0`$
$`4/3\mathrm{ln}(z_t)t_hz_hN_m+{\displaystyle \frac{99}{2}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_hS_0+16/3\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_hN_m5/4\mathrm{ln}(z_t)\mathrm{ln}(z_h)t_hz_h`$
$`+5\mathrm{ln}(z_t)\mathrm{ln}(z_h)t_{h}^{}{}_{}{}^{2}z_h8\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_h+{\displaystyle \frac{37}{2}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_h{\displaystyle \frac{115}{384}}Ft_hz_hr`$
$`+{\displaystyle \frac{5}{32}}Ft_hz_hr^2+2\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_hN_m8\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hN_m`$ (E.9)
$`A_{0,2}^Z={\displaystyle \frac{16}{9}}\pi ^2z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}+11/3z_{h}^{}{}_{}{}^{2}S_0N_m+{\displaystyle \frac{49}{128}}Fz_{h}^{}{}_{}{}^{2}r^4+{\displaystyle \frac{18911}{23040}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r^2`$
$`{\displaystyle \frac{112}{27}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}3/2\pi ^2t_hz_{h}^{}{}_{}{}^{2}+1/2\left(\mathrm{ln}(z_h)\right)^2z_{h}^{}{}_{}{}^{2}N_m+{\displaystyle \frac{5673}{10240}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r`$
$`{\displaystyle \frac{27}{4}}\left(\mathrm{ln}(z_t)\right)^2t_hz_{h}^{}{}_{}{}^{2}{\displaystyle \frac{1865}{768}}Ft_hz_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{25979}{30720}}\mathrm{ln}(z_t)t_hz_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{297}{2}}z_{h}^{}{}_{}{}^{2}S_2N_m`$
$`+18\mathrm{ln}(z_h)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{121}{144}}\pi ^2z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{487163}{15360}}t_hz_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{377}{8}}t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{7729}{46080}}z_{h}^{}{}_{}{}^{2}r{\displaystyle \frac{2557}{11520}}z_{h}^{}{}_{}{}^{2}r^2`$
$`22\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}S_0+{\displaystyle \frac{64}{3}}S_3z_{h}^{}{}_{}{}^{2}N_m+{\displaystyle \frac{25}{9}}\pi ^2z_{h}^{}{}_{}{}^{2}N_m40\zeta (3)z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{26047}{6144}}\mathrm{ln}(z_h)t_hz_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{433}{384}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r^3+{\displaystyle \frac{49}{32}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r^4+{\displaystyle \frac{70}{27}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{20}{9}}\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{9}{160}}\mathrm{ln}(z_h)z_tz_h+{\displaystyle \frac{2343}{160}}\mathrm{ln}(m_Z^2)t_hz_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{37}{4}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{3473}{54}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m`$
$`+{\displaystyle \frac{27}{2}}\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_{h}^{}{}_{}{}^{2}54\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{199}{18}}\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}N_m`$
$`+{\displaystyle \frac{16}{9}}\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}+{\displaystyle \frac{9}{160}}\mathrm{ln}(m_Z^2)z_tz_h+{\displaystyle \frac{27}{8}}\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)t_hz_{h}^{}{}_{}{}^{2}`$
$`9\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{20}{9}}\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m22\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}S_0N_m`$
$`{\displaystyle \frac{27}{8}}\mathrm{ln}(z_t)\mathrm{ln}(z_h)t_hz_{h}^{}{}_{}{}^{2}+9\mathrm{ln}(z_t)\mathrm{ln}(z_h)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{27}{4}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_{h}^{}{}_{}{}^{2}`$
$`+54\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}4/3\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m+16/3S_3z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{2746861}{51840}}z_{h}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{49}{192}}z_{h}^{}{}_{}{}^{2}r^3+{\displaystyle \frac{297}{8}}z_{h}^{}{}_{}{}^{2}S_2+{\displaystyle \frac{55}{6}}z_{h}^{}{}_{}{}^{2}S_0{\displaystyle \frac{7969}{72}}z_{h}^{}{}_{}{}^{2}N_m+{\displaystyle \frac{160}{81}}z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}7\zeta (3)z_{h}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{8413}{4320}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{5}{16}}\left(\mathrm{ln}(z_h)\right)^2z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{50209}{864}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{77}{6}}\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{38411}{1440}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{35}{48}}\left(\mathrm{ln}(z_t)\right)^2z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{3}{40}}z_tz_h+3/16Fz_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{17}{4}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{5673}{10240}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r{\displaystyle \frac{18911}{23040}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r^2+{\displaystyle \frac{433}{384}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r^3{\displaystyle \frac{49}{32}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r^4`$
$`+{\displaystyle \frac{55}{4}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}S_0+{\displaystyle \frac{14}{9}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}N_m+5/6\mathrm{ln}(z_t)\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{113}{24}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{53695}{6144}}Ft_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{27}{2}}Ft_{h}^{}{}_{}{}^{3}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{89}{768}}Fz_{h}^{}{}_{}{}^{2}r+{\displaystyle \frac{1441}{6144}}Fz_{h}^{}{}_{}{}^{2}r^2{\displaystyle \frac{241}{768}}Fz_{h}^{}{}_{}{}^{2}r^3`$ (E.10)
$`A_{1,1}^Z=\mathrm{ln}(z_t)\mathrm{ln}(z_h)z_h+1/2\mathrm{ln}(z_h)t_hz_h{\displaystyle \frac{165}{4}}\mathrm{ln}(m_Z^2)t_hz_hS_08/3\mathrm{ln}(m_Z^2)t_hz_hN_m`$
$`{\displaystyle \frac{9919}{288}}t_hz_h+{\displaystyle \frac{607}{3}}t_{h}^{}{}_{}{}^{2}z_h+{\displaystyle \frac{14}{9}}z_hN_m1/3\mathrm{ln}(z_h)z_h2\mathrm{ln}(m_Z^2)z_h+7/6\mathrm{ln}(z_t)z_h`$
$`+1/4Ft_{h}^{}{}_{}{}^{2}z_h2/9z_h2\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_h+2\left(\mathrm{ln}(m_Z^2)\right)^2z_h3/8\pi ^2t_hz_h+{\displaystyle \frac{201}{8}}t_hz_hS_0`$
$`165t_{h}^{}{}_{}{}^{2}z_hS_02\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_hN_m+2\left(\mathrm{ln}(m_Z^2)\right)^2z_hN_m{\displaystyle \frac{32}{3}}t_{h}^{}{}_{}{}^{2}z_hN_m`$
$`+3\mathrm{ln}(z_h)t_hz_hr+4/3\mathrm{ln}(z_h)z_hN_m+{\displaystyle \frac{1267}{24}}\mathrm{ln}(m_Z^2)t_hz_h{\displaystyle \frac{790}{3}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_h`$
$`11/3\mathrm{ln}(m_Z^2)z_hN_m{\displaystyle \frac{127}{8}}\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_h+57\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_h{\displaystyle \frac{625}{12}}\mathrm{ln}(z_t)t_hz_h`$
$`+{\displaystyle \frac{583}{3}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_h\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_h+{\displaystyle \frac{11}{8}}\left(\mathrm{ln}(z_t)\right)^2t_hz_h12\left(\mathrm{ln}(z_t)\right)^2t_{h}^{}{}_{}{}^{2}z_h`$
$`+165\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hS_0+{\displaystyle \frac{80}{3}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hN_m+4\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_hN_m`$
$`16\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_hN_m3\mathrm{ln}(z_t)t_hz_hr+{\displaystyle \frac{165}{4}}\mathrm{ln}(z_t)t_hz_hS_0+8/3\mathrm{ln}(z_t)t_hz_hN_m`$
$`165\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_hS_0{\displaystyle \frac{32}{3}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_hN_m+{\displaystyle \frac{29}{2}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_h45\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_h`$
$`+3/4Ft_hz_hr4\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_hN_m+16\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hN_m`$ (E.11)
$`A_{1,2}^Z=16/3\pi ^2z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}+{\displaystyle \frac{374}{9}}z_{h}^{}{}_{}{}^{2}S_0N_m{\displaystyle \frac{103}{54}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r^2+{\displaystyle \frac{32}{3}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}`$
$`+3\pi ^2t_hz_{h}^{}{}_{}{}^{2}\left(\mathrm{ln}(z_h)\right)^2z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{1321}{1080}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r+9\left(\mathrm{ln}(z_t)\right)^2t_hz_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{391}{288}}Ft_hz_{h}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{17969}{360}}\mathrm{ln}(z_t)t_hz_{h}^{}{}_{}{}^{2}720z_{h}^{}{}_{}{}^{2}S_2N_m{\displaystyle \frac{167}{54}}\pi ^2z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{1913}{180}}t_hz_{h}^{}{}_{}{}^{2}42t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{233}{540}}z_{h}^{}{}_{}{}^{2}r`$
$`+{\displaystyle \frac{23}{27}}z_{h}^{}{}_{}{}^{2}r^2+{\displaystyle \frac{242}{3}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}S_0{\displaystyle \frac{832}{9}}S_3z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{107}{9}}\pi ^2z_{h}^{}{}_{}{}^{2}N_m+176\zeta (3)z_{h}^{}{}_{}{}^{2}N_m`$
$`{\displaystyle \frac{577}{72}}\mathrm{ln}(z_h)t_hz_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{46}{9}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r^3{\displaystyle \frac{140}{27}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}N_m+{\displaystyle \frac{40}{9}}\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}`$
$`+1/5\mathrm{ln}(z_h)z_tz_h{\displaystyle \frac{509}{10}}\mathrm{ln}(m_Z^2)t_hz_{h}^{}{}_{}{}^{2}+12\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{595}{3}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m`$
$`18\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_{h}^{}{}_{}{}^{2}+72\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{107}{3}}\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}N_m`$
$`16/3\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}1/5\mathrm{ln}(m_Z^2)z_tz_h+{\displaystyle \frac{40}{9}}\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m`$
$`+{\displaystyle \frac{242}{3}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}S_0N_m+9\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_{h}^{}{}_{}{}^{2}72\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{40}{9}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{208}{9}}S_3z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{475837}{3240}}z_{h}^{}{}_{}{}^{2}180z_{h}^{}{}_{}{}^{2}S_2{\displaystyle \frac{289}{18}}z_{h}^{}{}_{}{}^{2}S_0`$
$`+{\displaystyle \frac{4091}{12}}z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{80}{27}}z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}+{\displaystyle \frac{112}{3}}\zeta (3)z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{548}{135}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}1/2\left(\mathrm{ln}(z_h)\right)^2z_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{39517}{216}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}+41\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{51127}{540}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}7/6\left(\mathrm{ln}(z_t)\right)^2z_{h}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{4}{15}}z_tz_h1/2Fz_{h}^{}{}_{}{}^{2}48\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{1321}{1080}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r+{\displaystyle \frac{103}{54}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r^2`$
$`{\displaystyle \frac{46}{9}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r^3{\displaystyle \frac{619}{12}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}S_0{\displaystyle \frac{116}{27}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{20}{9}}\mathrm{ln}(z_t)\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{565}{36}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{577}{72}}Ft_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{71}{288}}Fz_{h}^{}{}_{}{}^{2}r{\displaystyle \frac{7}{12}}Fz_{h}^{}{}_{}{}^{2}r^2+{\displaystyle \frac{23}{18}}Fz_{h}^{}{}_{}{}^{2}r^3`$ (E.12)
$`A_{2,1}^Z=4/3\mathrm{ln}(z_t)\mathrm{ln}(z_h)z_h2/3\mathrm{ln}(z_h)t_hz_h+{\displaystyle \frac{183}{4}}\mathrm{ln}(m_Z^2)t_hz_hS_0+{\displaystyle \frac{32}{9}}\mathrm{ln}(m_Z^2)t_hz_hN_m`$
$`+{\displaystyle \frac{1567}{36}}t_hz_h{\displaystyle \frac{2138}{9}}t_{h}^{}{}_{}{}^{2}z_h{\displaystyle \frac{56}{27}}z_hN_m+4/9\mathrm{ln}(z_h)z_h+8/3\mathrm{ln}(m_Z^2)z_h{\displaystyle \frac{14}{9}}\mathrm{ln}(z_t)z_h`$
$`1/3Ft_{h}^{}{}_{}{}^{2}z_h+{\displaystyle \frac{8}{27}}z_h+8/3\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_h8/3\left(\mathrm{ln}(m_Z^2)\right)^2z_h{\displaystyle \frac{109}{4}}t_hz_hS_0`$
$`+183t_{h}^{}{}_{}{}^{2}z_hS_0+8/3\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_hN_m8/3\left(\mathrm{ln}(m_Z^2)\right)^2z_hN_m+{\displaystyle \frac{128}{9}}t_{h}^{}{}_{}{}^{2}z_hN_m`$
$`4\mathrm{ln}(z_h)t_hz_hr{\displaystyle \frac{16}{9}}\mathrm{ln}(z_h)z_hN_m{\displaystyle \frac{518}{9}}\mathrm{ln}(m_Z^2)t_hz_h+{\displaystyle \frac{2324}{9}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_h`$
$`+{\displaystyle \frac{44}{9}}\mathrm{ln}(m_Z^2)z_hN_m+{\displaystyle \frac{31}{6}}\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_h{\displaystyle \frac{62}{3}}\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_h+{\displaystyle \frac{524}{9}}\mathrm{ln}(z_t)t_hz_h`$
$`{\displaystyle \frac{2042}{9}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_h+4/3\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_h8/3\left(\mathrm{ln}(z_t)\right)^2t_hz_h+{\displaystyle \frac{32}{3}}\left(\mathrm{ln}(z_t)\right)^2t_{h}^{}{}_{}{}^{2}z_h`$
$`183\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hS_0{\displaystyle \frac{320}{9}}\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hN_m16/3\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_hN_m`$
$`+{\displaystyle \frac{64}{3}}\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_hN_m+4\mathrm{ln}(z_t)t_hz_hr{\displaystyle \frac{183}{4}}\mathrm{ln}(z_t)t_hz_hS_0{\displaystyle \frac{32}{9}}\mathrm{ln}(z_t)t_hz_hN_m`$
$`+183\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_hS_0+{\displaystyle \frac{128}{9}}\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_hN_m5/2\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_h+10\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_h`$
$`Ft_hz_hr+16/3\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_hN_m{\displaystyle \frac{64}{3}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_hN_m`$ (E.13)
$`A_{2,2}^Z=16/3\pi ^2z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}{\displaystyle \frac{7295}{27}}z_{h}^{}{}_{}{}^{2}S_0N_m+{\displaystyle \frac{206}{81}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r^2{\displaystyle \frac{32}{27}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}`$
$`3/2\pi ^2t_hz_{h}^{}{}_{}{}^{2}+4/3\left(\mathrm{ln}(z_h)\right)^2z_{h}^{}{}_{}{}^{2}N_m+{\displaystyle \frac{1321}{810}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r9/2\left(\mathrm{ln}(z_t)\right)^2t_hz_{h}^{}{}_{}{}^{2}{\displaystyle \frac{391}{216}}Ft_hz_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{34453}{540}}\mathrm{ln}(z_t)t_hz_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{2655}{2}}z_{h}^{}{}_{}{}^{2}S_2N_m+{\displaystyle \frac{995}{216}}\pi ^2z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{4208}{135}}t_hz_{h}^{}{}_{}{}^{2}+39t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{233}{405}}z_{h}^{}{}_{}{}^{2}r`$
$`{\displaystyle \frac{92}{81}}z_{h}^{}{}_{}{}^{2}r^2{\displaystyle \frac{638}{9}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}S_0+{\displaystyle \frac{1408}{9}}S_3z_{h}^{}{}_{}{}^{2}N_m+{\displaystyle \frac{1025}{54}}\pi ^2z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{860}{3}}\zeta (3)z_{h}^{}{}_{}{}^{2}N_m`$
$`+{\displaystyle \frac{577}{54}}\mathrm{ln}(z_h)t_hz_{h}^{}{}_{}{}^{2}{\displaystyle \frac{184}{27}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}r^3+{\displaystyle \frac{560}{81}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{160}{27}}\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{4}{15}}\mathrm{ln}(z_h)z_tz_h+{\displaystyle \frac{3457}{60}}\mathrm{ln}(m_Z^2)t_hz_{h}^{}{}_{}{}^{2}42\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{3583}{27}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m`$
$`+9\left(\mathrm{ln}(m_Z^2)\right)^2t_hz_{h}^{}{}_{}{}^{2}36\left(\mathrm{ln}(m_Z^2)\right)^2t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{97}{3}}\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}N_m`$
$`+16/3\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}+{\displaystyle \frac{4}{15}}\mathrm{ln}(m_Z^2)z_tz_h{\displaystyle \frac{160}{27}}\mathrm{ln}(z_h)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m`$
$`{\displaystyle \frac{638}{9}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}S_0N_m9/2\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_hz_{h}^{}{}_{}{}^{2}+36\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{112}{27}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}N_m+{\displaystyle \frac{352}{9}}S_3z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{1126103}{19440}}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{2655}{8}}z_{h}^{}{}_{}{}^{2}S_2{\displaystyle \frac{1744}{27}}z_{h}^{}{}_{}{}^{2}S_0`$
$`{\displaystyle \frac{68035}{324}}z_{h}^{}{}_{}{}^{2}N_m{\displaystyle \frac{1040}{81}}z_{h}^{}{}_{}{}^{2}N_{m}^{}{}_{}{}^{2}{\displaystyle \frac{533}{9}}\zeta (3)z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{2192}{405}}\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}+2/3\left(\mathrm{ln}(z_h)\right)^2z_{h}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{3458}{27}}\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{113}{3}}\left(\mathrm{ln}(m_Z^2)\right)^2z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{206869}{2160}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{23}{108}}\left(\mathrm{ln}(z_t)\right)^2z_{h}^{}{}_{}{}^{2}`$
$`{\displaystyle \frac{16}{45}}z_tz_h+2/3Fz_{h}^{}{}_{}{}^{2}+60\mathrm{ln}(z_t)t_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{1321}{810}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r{\displaystyle \frac{206}{81}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r^2`$
$`+{\displaystyle \frac{184}{27}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}r^3+{\displaystyle \frac{983}{18}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}S_0{\displaystyle \frac{64}{81}}\mathrm{ln}(z_t)z_{h}^{}{}_{}{}^{2}N_m+{\displaystyle \frac{80}{27}}\mathrm{ln}(z_t)\mathrm{ln}(z_h)z_{h}^{}{}_{}{}^{2}`$
$`+{\displaystyle \frac{437}{27}}\mathrm{ln}(z_t)\mathrm{ln}(m_Z^2)z_{h}^{}{}_{}{}^{2}{\displaystyle \frac{577}{54}}Ft_{h}^{}{}_{}{}^{2}z_{h}^{}{}_{}{}^{2}+{\displaystyle \frac{71}{216}}Fz_{h}^{}{}_{}{}^{2}r+{\displaystyle \frac{7}{9}}Fz_{h}^{}{}_{}{}^{2}r^2{\displaystyle \frac{46}{27}}Fz_{h}^{}{}_{}{}^{2}r^3`$ (E.14)