text:
# 1 Introduction

## 1 Introduction

Two years ago the kernel of the BFKL equation  was obtained in the next-to-leading order (NLO)  for the case of forward scattering, i.e. for the momentum transfer $`t=0`$ and the singlet colour representation in the $`t`$-channel. This long awaited event had led to the appearance of a number of papers (see, for instance,  and references therein) devoted to the problem of possible applications of the NLO result in the physics of semi-hard processes.

An important way of development of the BFKL approach became the generalization of the obtained results to the non-forward scattering . For singlet colour representation in the $`t`$-channel, the generalized approach can be used directly for the description of a wide circle of physical processes. The generalization for non-singlet colour states is also of great importance, especially for the antisymmetric colour octet state of two Reggeized gluons in the $`t`$-channel. This case is especially important since the BFKL approach is based on the gluon Reggeization, although the Reggeization is not proved in the NLO. Therefore, at this order the Reggeization is a hypothesis which must be carefully checked. This can be done using the “bootstrap” equations  appearing from the requirement of the compatibility of the gluon Reggeization with the $`s`$-channel unitarity. In the BFKL approach the scattering amplitude for the high energy process $`A+BA^{}+B^{}`$ is presented as the convolution of the impact factors $`\mathrm{\Phi }_{A^{}A}`$ and $`\mathrm{\Phi }_{B^{}B}`$, which describe the $`AA^{}`$ and $`BB^{}`$ transitions in the particle-Reggeon scattering, and the Green’s function $`G`$ for the Reggeon-Reggeon scattering. The bootstrap conditions for the impact factors are already checked both for the gluon and quark scattering, for helicity conserving and helicity non-conserving amplitudes . We remind that the impact factors are infrared finite  for colourless particles only. Considering scattering amplitudes at the parton level, one needs to use an infrared regularization. We use the dimensional regularization, which is commonly adopted. Remarkably, the bootstrap conditions for the impact factors are satisfied at arbitrary space-time dimension $`D=4+2ϵ`$.

The kernel of the BFKL equation in the octet channel must also satisfy the bootstrap condition (which is called “first bootstrap condition”). This condition was checked (also for arbitrary $`D`$) in the part concerning the quark contribution to the kernel . Of course, the most important (and most complicated for verification) is the gluon part of the first bootstrap condition. The gluon part of the kernel in the NLO is expressed in terms of the two-loop gluon trajectory, the one-loop effective vertex for the gluon production in the Reggeon-Reggeon collisions (RRG-vertex) and the two-gluon production contribution. The last contribution was obtained recently . Since the two-loop gluon trajectory  and the one-loop RRG-vertex  are already known, it seems that the bootstrap condition could be checked. But in the bootstrap equation

$$\frac{g^2Nt}{2(2\pi )^{D1}}\frac{d^{D2}q_1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_1\stackrel{}{q})^2}\frac{d^{D2}q_2}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_2\stackrel{}{q})^2}𝒦^{(8)(1)}(\stackrel{}{q}_1,\stackrel{}{q}_2;\stackrel{}{q})=\omega ^{(1)}(t)\omega ^{(2)}(t),$$
(1.1)
where $`g`$ is the coupling constant, $`N`$ is the number of colour, $`\stackrel{}{q}`$ is the transverse momentum transfer, $`t=\stackrel{}{q}^{\mathrm{\hspace{0.25em}2}}`$, $`𝒦^{(8)(1)}`$ is the one-loop contribution to the non-forward octet kernel, $`\omega ^{(1)}`$ and $`\omega ^{(2)}`$ are the one- and two-loop contributions to the gluon trajectory, the integration over the transverse momenta $`\stackrel{}{q}_{1,2}`$ is singular. Therefore, in order to check at least the terms finite for $`ϵ0`$ in $`\omega ^{(2)}`$ in the R.H.S. of Eq. (1.1), one needs to know the kernel in regions of singularities at arbitrary $`ϵ`$ (because, for example, the region of arbitrary small $`\stackrel{}{q}_1`$ does contribute, so that $`(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ`$ cannot be expanded in powers of $`ϵ`$). Unfortunately, the kernel with such accuracy is not known yet, although for the gluon trajectory  and for the two-gluon contribution to the octet kernel  the integral representations for arbitrary $`D`$ are known. The problem is in the one-gluon contribution or, in other words, in the RRG-vertex. Let us note that, from now on, we will talk only about the gluon part, i.e. consider pure gluodynamics, since the quark part of the vertex is known at arbitrary $`D`$ . Instead, in the gluon part firstly only the terms finite for $`ϵ0`$ were kept , but in the process of calculation of the forward BFKL kernel, it was understood that in this kernel the RRG-vertex at small transverse momenta $`\stackrel{}{k}`$ of the produced gluon must be known at arbitrary $`D`$. After this, the vertex at small $`\stackrel{}{k}`$ was calculated  at arbitrary $`D`$. Later the results of  were obtained by another method in . But for the verification of the bootstrap equation (1.1) it must be known at arbitrary $`D`$ in a wider kinematical region. Therefore it became clear the necessity of the knowledge of the RRG-vertex at arbitrary $`D`$, especially taking into account that it can be used not only for the check of the bootstrap, but, for example, in the Odderon problem in the NLO and so on.

In this paper we calculate the RRG-vertex in the NLO and its contribution to the non-forward BFKL kernel at arbitrary $`D`$. In the next Section we calculate the amplitude of the gluon production in the multi-Regge kinematics for gluon-gluon collisions. The RRG-vertices are defined and calculated in Section 3 and the one-gluon contribution to the BFKL kernel in Section 4.

## 2 The gluon production amplitude

The RRG-vertex can be obtained from the amplitudes of the gluon production in the multi-Regge kinematics (in other words, in the central kinematical region) at collisions of any pair of particles. We will consider the gluon-gluon collisions and use intermediate results obtained in  for this process, which are valid at arbitrary $`D`$. We will use the denotations $`p_A`$ and $`p_B`$ ($`p_A^{}`$ and $`p_B^{}`$) for the momenta of the incoming (outgoing) gluons, $`k`$ and $`e(k)`$ for momentum and polarization vector of the produced gluon; $`q_1=p_Ap_A^{}`$ and $`q_2=p_B^{}p_B`$ are the momentum transfers, so that $`k=q_1q_2`$. All the denotations are the same as used in , except that for the momentum of the produced gluon, which was called $`p_D`$ there. The kinematics is defined by the relations

$$ss_1s_2|t_1||t_2|,$$
(2.1)
where

$$s=(p_A+p_B)^2,s_1=(p_A^{}+k)^2,s_2=(p_B^{}+k)^2,t_{1,2}=q_{1,2}^2.$$
(2.2)
In terms of the parameters of the Sudakov decomposition

$$k=\beta p_A+\alpha p_B+k_{},q_i=\beta _ip_A+\alpha _ip_B+q_{i}^{}{}_{}{}^{},$$
(2.3)
the relations (2.1) give

$$1\beta \beta _1\alpha _1\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{s},\mathrm{\hspace{0.33em}\hspace{0.33em}\hspace{0.33em}\hspace{0.33em}\hspace{0.33em}\hspace{0.33em}\hspace{0.33em}\hspace{0.33em}\hspace{0.33em}\hspace{0.33em}1}\alpha \alpha _2\beta _2\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{s},$$

$$s_1s\alpha ,s_2s\beta ,\stackrel{}{k}^2=k_{}^2\frac{s_1s_2}{s}.$$
(2.4)
Here and below the vector sign is used for the components of the momenta transverse to the plane of the momenta of the initial particles $`p_A`$ and $`p_B`$.

Since we are interested in the RRG-vertex, we will consider the amplitudes with conservation of the helicities of the scattered gluons (helicity non-conservation in the NLO is related with the gluon-gluon-Reggeon vertices only, which are already known at arbitrary $`D`$ ). This amplitude can be presented in the form (cf. )

$$A_{23}=2sg^3T_{A^{}A}^{c_1}\frac{1}{t_1}T_{c_2c_1}^d\frac{1}{t_2}T_{B^{}B}^{c_2}e_\mu ^{}(k)𝒜^\mu (q_2,q_1),$$
(2.5)
where $`T_{bc}^a`$ are matrix elements of the colour group generator in the adjoint representation and the amplitude $`𝒜^\mu `$ in the Born approximation is equal to $`C^\mu (q_2,q_1)`$:

$$𝒜_{\mathrm{B}orn}^\mu =C^\mu (q_2,q_1)=q_{1}^{}{}_{}{}^{\mu }q_{2}^{}{}_{}{}^{\mu }+\frac{p_A^\mu }{s_1}(\stackrel{}{k}^22\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})\frac{p_B^\mu }{s_2}(\stackrel{}{k}^22\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}).$$
(2.6)
It was shown in  that at one-loop order the amplitude can be presented as the sum of three contributions, which come from integrations over the momenta of the virtual gluons in three different kinematical regions: the region of fragmentation of the particle $`A`$, the central region and the region of fragmentation of the particle $`B`$. Correspondingly, we represent $`𝒜^\mu `$ as

$$𝒜^\mu =𝒜_{\mathrm{B}orn}^\mu +𝒜_A^\mu +𝒜_{\mathrm{a}s}^\mu +𝒜_B^\mu ,$$
(2.7)
where for the contribution of the fragmentation of the particle $`A`$, we obtain from Eq. (63) of

$$𝒜_A^\mu =C^\mu (q_2,q_1)\overline{g}^2(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ\frac{\mathrm{\Gamma }^2(1+ϵ)}{\mathrm{\Gamma }(4+2ϵ)}\frac{(2)}{ϵ^2}\left(3(1+ϵ)^2+ϵ^2\right),$$
(2.8)
where $`\overline{g}^2=g^2N\mathrm{\Gamma }(1ϵ)/(4\pi )^{2+ϵ}`$. In the derivation of (2.8) from Eq. (63) of , we have used the relations

$$e_{A\lambda }^\alpha e_{A^{}\lambda }^\alpha ^{}\left(\delta _{\alpha \alpha ^{}}(D2)\frac{q_{1}^{}{}_{\alpha }{}^{}q_{1}^{}{}_{\alpha ^{}}{}^{}}{q_{1}^{}{}_{}{}^{2}}\right)=0,e_{A\lambda }^\alpha e_{A^{}\lambda }^\alpha ^{}\delta _{\alpha \alpha ^{}}=1,$$
(2.9)
for the polarization vectors $`e_{A\lambda }`$ and $`e_{A^{}\lambda }`$ of the particles $`A`$ and $`A^{}`$ with the definite helicity $`\lambda `$ in the $`D`$-dimensional space-time. The analogous relations exist for the particle $`B`$. Evidently, the contribution from the fragmentation region of the particle $`B`$ can be obtained from (2.8) by the substitution $`\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}`$.

The contribution of the central region is given by Eqs. (77), (80) and (81) of . At this step in  was not yet made an expansion in $`ϵ`$. Let us pay attention that in these formulae the anti-symmetrization with respect to the substitution $`p_Bp_B^{}`$ must be done. The formulae (80) and (81) of  can be greatly simplified and the contribution $`𝒜_{\mathrm{a}s}^\mu `$, with accuracy up to terms proportional to $`k^\mu `$ and therefore not contributing to (2.7), can be presented as

$$𝒜_{\mathrm{a}s}^\mu =\frac{2g^2N}{(2\pi )^D}\widehat{𝒮}\{(C^\mu (q_2,q_1)2\stackrel{}{k}^2𝒫^\mu )[s_1t_1I_{4A}\frac{4(D2)}{t_1}\frac{p_A^\rho p_B^\sigma }{s}I_2^{\rho \sigma }(q_1)2t_1\frac{p_B^\rho }{s}I_{3A}^\rho $$

$$+4\stackrel{}{k}^2\frac{p_A^\rho }{s_1}I_3^\rho 2t_1I_{3A}+\stackrel{}{k}^2I_32I_2(q_1)]+4𝒫^\mu [\frac{s_1s_2}{4}t_1t_2I_5+\frac{s_1t_1}{2}(\stackrel{}{k}^2t_2)I_{4A}2\stackrel{}{k}^2t_1\frac{p_A^\rho }{s_1}I_3^\rho $$

$$\stackrel{}{k}^2t_2I_3\stackrel{}{k}^2(2I_2(q_1)+t_1I_{3A})t_2I_2(q_1)]2(g^{\mu \rho }2\frac{p_B^\mu k^\rho }{s_2}\left)\right[\frac{t_1s_2}{s}I_{3A}^\rho $$

$$+2(D2)\frac{\stackrel{}{k}^2}{t_1}\frac{p_A^\sigma }{s_1}I_2^{\rho \sigma }(q_1)]4(D2)\frac{p_A^\rho p_B^\sigma }{s}[I_3^{\mu \rho \sigma }+I_3^{\mu \rho }k^\sigma 2\frac{p_A^\mu }{s_1}I_2^{\rho \sigma }(q_1)]$$

$$st_1t_2[2\frac{p_B^\mu }{s_2}I_4^{(1)}+I_5^\mu ]+2(\stackrel{}{k}^2+t_1+t_2)[2\frac{p_A^\mu }{s_1}I_2(q_1)I_3^\mu ]\},$$
(2.10)
where

$$𝒫^\mu =\frac{p_A^\mu }{s_1}\frac{p_B^\mu }{s_2}$$
(2.11)
and $`\widehat{𝒮}`$ is the operator of symmetrization with respect to each of the substitutions $`p_Ap_A`$ and $`p_Bp_B`$ and anti-symmetrization with respect to the substitution

$$p_Ap_B,q_1q_2.$$
(2.12)
The $`I`$’s in Eq. (2.10) are the integrals introduced in :

$$I_2^{\mu _1\mathrm{}\mu _n}(q)=\frac{1}{i}d^Dp\frac{p^{\mu _1}\mathrm{}p^{\mu _n}}{(p^2+i\epsilon )[(p+q)^2+i\epsilon ]},$$

$$I_3^{\mu _1\mathrm{}\mu _n}=\frac{1}{i}d^Dp\frac{p^{\mu _1}\mathrm{}p^{\mu _n}}{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+k)^2+i\epsilon ]},$$

$$I_{3A}^{\mu _1\mathrm{}\mu _n}=\frac{1}{i}d^Dp\frac{p^{\mu _1}\mathrm{}p^{\mu _n}}{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+p_A)^2+i\epsilon ]},$$

$$I_4^{(1)}=\frac{1}{i}d^Dp\frac{1}{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+p_A)^2+i\epsilon ][(pp_B)^2+i\epsilon ]},$$
(2.13)
$$I_{4A}=\frac{1}{i}d^Dp\frac{1}{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+q_2)^2+i\epsilon ][(p+p_A)^2+i\epsilon ]},$$

$$I_5=\frac{1}{i}d^Dp\frac{1}{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+q_2)^2+i\epsilon ][(p+p_A)^2+i\epsilon ][(pp_B)^2+i\epsilon ]},$$

$$I_5^\mu =\frac{1}{i}d^Dp\frac{(p+q_1)^\mu }{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+q_2)^2+i\epsilon ][(p+p_A)^2+i\epsilon ][(pp_B)^2+i\epsilon ]}.$$

In the representation (2.10) the transversality to the momenta $`k^\mu `$ (which guarantees the gauge invariance of the amplitude) of the first three terms with square brackets is evident, while the transversality of the last three terms can be easily checked with account of the anti-symmetrization (2.12).

The integrals with two or three denominators ($`I_2^{\mu _1\mathrm{}\mu _n}`$, $`I_3^{\mu _1\mathrm{}\mu _n}`$, $`I_{3A}^{\mu _1\mathrm{}\mu _n}`$) were calculated in  at arbitrary $`D`$, whereas for the other integrals the expansion in $`ϵ`$ was used. It occurs that, in the multi-Regge kinematics (2.1), the integral $`I_4^{(1)}`$ can also be calculated at arbitrary $`D`$, as it is shown in Appendix A. Unfortunately, the integrals $`I_{4A}`$, $`I_5`$ and $`I_5^\mu `$ cannot be expressed in terms of elementary functions at arbitrary $`D`$. It is possible, however, to perform the integration over the longitudinal variables in the Sudakov decomposition for $`p`$ and to express them in terms of the same $`(D2)`$-dimensional integrals over $`\stackrel{}{p}`$, which appear in the two-gluon contribution to the kernel . This is done in Appendix A, where, for completeness, the integrals with two and three denominators are also given.

The integrals which cannot be expressed through elementary functions, and will be used below, are the following:

$$_{4A}=_0^1\frac{dx}{1x}\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}[\frac{x}{[(1x)\stackrel{}{k}_1^2+x(\stackrel{}{k}_1+\stackrel{}{q}_2)^2](\stackrel{}{k}_1x\stackrel{}{k})^2}\frac{1}{(\stackrel{}{k}_1+\stackrel{}{q}_2)^2(\stackrel{}{k}_1\stackrel{}{k})^2}],$$

$$_5=_0^1\frac{dx}{1x}\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{1}{\stackrel{}{k}_1^2[(1x)\stackrel{}{k}_1^2+x(\stackrel{}{k}_1\stackrel{}{q}_1)^2]}\left(\frac{x^2}{(\stackrel{}{k}_1x\stackrel{}{k})^2}\frac{1}{(\stackrel{}{k}_1\stackrel{}{k})^2}\right),$$

$$_3=\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{1}{\stackrel{}{k}_1^2(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2}\mathrm{ln}\left(\frac{(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2}{\stackrel{}{k}^2\stackrel{}{k}_1^2}\right),$$

$$_3=\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{1}{\stackrel{}{k}_1^2(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2},$$
(2.14)
as well as $`_{4B}`$, which is obtained from $`_{4A}`$ by the substitution $`q_1q_2`$, and $`_5^\mu `$, $`_3^\mu `$ and $`_3^\mu `$, the first of which differs from the integral $`_5`$ by the factor $`(k_1q_1)_{}^\mu `$ and the other two from $`_3`$ and $`_3`$, respectively, by the factor $`k_{1}^{}{}_{}{}^{\mu }`$, in the corresponding integrands.

Using the results of the Appendix A, one gets

$$𝒜_{\mathrm{a}s}^\mu =\overline{g}^2\left\{2t_1t_2_5^\mu +r_{\mathrm{a}s}C^\mu (q_2,q_1)+2t_1r_A\frac{p_A^\mu }{s_1}2t_2r_B\frac{p_B^\mu }{s_2}\right\},$$
(2.15)
where

$$_5^\mu =_5^\mu +_3^\mu +\frac{1}{2}\mathrm{ln}\left(\frac{s(s)(\stackrel{}{k}^2)^2}{s_1(s_1)s_2(s_2)}\right)_3^\mu ,$$

$$r_{\mathrm{a}s}=\{\frac{t_1t_2}{2}_5+t_2_{4B}+\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ[\frac{1}{2}\mathrm{ln}\left(\frac{s_1(s_1)}{t_1^2}\right)+2\psi (ϵ)\psi (2ϵ)\psi (1ϵ)$$

$$+\frac{1}{2ϵ(1+2ϵ)(3+2ϵ)}(\frac{t_1(3+14ϵ+8ϵ^2)t_2(3+3ϵ+ϵ^2)}{t_1t_2}$$

$$+\frac{\stackrel{}{k}^2t_1ϵ}{(t_1t_2)^3}((2+ϵ)t_2ϵt_1))]\}+\{AB\},$$

$$_5=_5_3\frac{1}{2}\mathrm{ln}\left(\frac{s(s)(\stackrel{}{k}^2)^2}{s_1(s_1)s_2(s_2)}\right)_3,$$

$$r_A=t_2(t_1_5+_{4B})+\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ[\frac{1}{2}\mathrm{ln}\left(\frac{s_1(s_1)s_2(s_2)}{s(s)t_2^2}\right)\psi (1)\psi (ϵ)$$

$$.+\psi (1ϵ)+\psi (2ϵ)]+\frac{\mathrm{\Gamma }^2(ϵ)}{2t_1(1+2ϵ)\mathrm{\Gamma }(2ϵ)(3+2ϵ)}\{((\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ+1}[\frac{t_2}{t_1t_2}(11+7ϵ)$$

$$+\frac{\stackrel{}{k}^2}{(t_1t_2)^3}(t_2(t_1+t_2)ϵt_1(t_1t_2))+\frac{(\stackrel{}{k}^2)^2}{(t_1t_2)^3}((2+ϵ)t_2ϵt_1)])+(AB)\},$$

$$r_B=r_A(AB).$$
(2.16)
Here and below $`AB`$ means the substitution (2.12). The amplitude $`𝒜^\mu `$ must be anti-symmetric under this substitution. The amplitude (2.15) has this property since $`C^\mu (q_2,q_1)`$ (Eq. (2.6)) and $`_5^\mu `$ (Eqs. (2.16)) change their sign under this substitution.

Note that in all physical regions we have $`s_1s_2/s=\stackrel{}{k}^2`$, so that

$$\mathrm{ln}\left(\frac{s(s)(\stackrel{}{k}^2)^2}{(s_1(s_1)s_2(s_2))}\right)=i\pi .$$
(2.17)
Below we will work in the physical region of the $`s`$-channel and use this relation.

The total amplitude $`𝒜^\mu `$ is defined by (2.7), (2.8) and (2.15) and can be represented as

$$𝒜^\mu =C^\mu (q_2,q_1)\left(1+\overline{g}^2r\right)+\overline{g}^2\left\{2t_1t_2_5^\mu +2t_1r_A\frac{p_A^\mu }{s_1}2t_2r_B\frac{p_B^\mu }{s_2}\right\},$$
(2.18)
where

$$r=r_{\mathrm{a}s}\frac{2}{ϵ^2}\frac{\mathrm{\Gamma }^2(1+ϵ)}{\mathrm{\Gamma }(4+2ϵ)}\left(3(1+ϵ)^2+ϵ^2\right)\left((\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ+(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ\right)$$
(2.19)
and $`r_{\mathrm{a}s}`$, $`r_A`$ and $`r_B`$ are defined in (2.16). The gauge invariance of the amplitude (2.5) follows from the properties of the amplitudes (2.8) and (2.15):

$$k_\mu 𝒜_A^\mu =k_\mu 𝒜_B^\mu =k_\mu 𝒜_{\mathrm{a}s}^\mu =0.$$
(2.20)
The first two of these relations follow evidently from the transversality of $`C^\mu (q_2,q_1)`$ (Eq. (2.6)), whereas the last of them is not so trivial and is fulfilled due to the equality:

$$2k_\mu _5^\mu f_{_{}}=[t_1_5_{4B}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^{ϵ1}(\frac{1}{2}\mathrm{ln}\left(\frac{s_1(s_1)s_2(s_2)}{s(s)t_2^2}\right)$$

$$..+\psi (1ϵ)\psi (ϵ)+\psi (2ϵ)\psi (1))][AB],$$
(2.21)
which is derived in the Appendix B. The amplitude can be written in an explicit gauge invariant form if we use an analogous relation, also obtained in the Appendix B:

$$2(q_1+q_2)_\mu _5^\mu f__+=\{t_1_5_{4B}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}[(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^{ϵ1}(\frac{1}{2}\mathrm{ln}\left(\frac{s_1(s_1)s_2(s_2)}{s(s)t_2^2}\right)$$

$$.+\psi (1ϵ)\psi (ϵ)+\psi (2ϵ)\psi (1))$$

$$(\stackrel{}{k}^2)^{ϵ1}(\frac{1}{2}\mathrm{ln}\left(\frac{s_1(s_1)s_2(s_2)}{s(s)(\stackrel{}{k}^2)^2}\right)+\psi (1ϵ)\psi (ϵ))]\}+\{AB\}.$$
(2.22)
Since the general form for $`_5^\mu `$ is

$$_5^\mu =r_{_{}}k_{}^\mu +r__+(q_1+q_2)_{}^\mu ,$$
(2.23)
we can express $`r_{_{}}`$ and $`r__+`$ in terms of $`f_{_{}}`$ and $`f__+`$:

$$r__+=\frac{f_{_{}}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})f__+\stackrel{}{k}^2}{8(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_1\stackrel{}{q}_2)^2)},r_{_{}}=\frac{f__+(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})f_{_{}}(\stackrel{}{q}_1+\stackrel{}{q}_2)^2}{8(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_1\stackrel{}{q}_2)^2)}.$$
(2.24)
Rewriting then (2.23) as

$$_5^\mu =r_{_{}}k^\mu r__+\left[C^\mu (q_2,q_1)+2\stackrel{}{q}_1\stackrel{}{q}_2\left(\frac{p_A^\mu }{s_1}\frac{p_B^\mu }{s_2}\right)\right]+\frac{f_{_{}}}{2}\left(\frac{p_A^\mu }{s_1}+\frac{p_B^\mu }{s_2}\right)$$
(2.25)
and using this equality in the expression (2.15) for $`𝒜_{\mathrm{a}s}^\mu `$, we obtain the explicitly gauge invariant form for the amplitude $`𝒜^\mu `$:

$$𝒜^\mu =C^\mu (q_2,q_1)(1+\overline{g}^2r_C)+𝒫^\mu \overline{g}^2\mathrm{\hspace{0.17em}2}t_1t_2r_𝒫,$$
(2.26)
where the terms proportional to $`k^\mu `$ were omitted and

$$r_C=\{t_1t_2(r__++\frac{_5}{2})+t_2_{4B}+\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ[\frac{1}{2}\mathrm{ln}\left(\frac{s_1(s_1)}{t_1^2}\right)+2\psi (ϵ)\psi (2ϵ)$$

$$\psi (1ϵ)+\frac{1}{2ϵ(1+2ϵ)(3+2ϵ)}(3(1+ϵ)\frac{ϵ^2}{1+ϵ}+\frac{t_1(3+14ϵ+8ϵ^2)t_2(3+3ϵ+ϵ^2)}{t_1t_2}$$

$$+\frac{\stackrel{}{k}^2t_1ϵ}{(t_1t_2)^3}((2+ϵ)t_2ϵt_1)]\}+\{AB\},$$

$$r_𝒫=\{(\stackrel{}{q}_1\stackrel{}{q}_2)r__+\frac{t_1+t_2}{4}_5\frac{_{4B}}{2}+\frac{1}{2}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ1}(\frac{1}{2}\mathrm{ln}\left(\frac{s_1(s_1)s_2(s_2)}{s(s)t_1^2}\right)+\psi (1)$$

$$+\psi (ϵ)\psi (1ϵ)\psi (2ϵ)+\frac{t_1}{t_2(1+2ϵ)(3+2ϵ)}[\frac{t_2}{t_1t_2}(11+7ϵ)$$

$$+\frac{\stackrel{}{k}^2}{(t_1t_2)^3}(t_2(t_1+t_2)ϵt_1(t_1t_2))+\frac{(\stackrel{}{k}^2)^2}{(t_1t_2)^3}((2+ϵ)t_2ϵt_1)])\}+\{AB\}.$$
(2.27)
The functions $`r__+`$, $`_5`$ and $`_{4B}`$ entering Eq. (2.27) are given by (2.21)-(2.24), (2.16) and (2.14), respectively.

## 3 The Reggeon-Reggeon-gluon vertices

Since the production amplitude must not have simultaneously discontinuities in the overlapping channels $`s_1`$ and $`s_2`$, it cannot be a simple generalization of the Regge form for the elastic amplitude. Instead, it has the form

$$8g^3𝒜^\mu =\mathrm{\Gamma }(t_1;\stackrel{}{k}^2)\mathrm{\Gamma }(t_2;\stackrel{}{k}^2)\{[\left(\frac{s_1}{\stackrel{}{k}^2}\right)^{\omega _1\omega _2}+\left(\frac{s_1}{\stackrel{}{k}^2}\right)^{\omega _1\omega _2}][\left(\frac{s}{\stackrel{}{k}^2}\right)^{\omega _2}+\left(\frac{s}{\stackrel{}{k}^2}\right)^{\omega _2}]R^\mu $$

$$+[\left(\frac{s_2}{\stackrel{}{k}^2}\right)^{\omega _2\omega _1}+\left(\frac{s_2}{\stackrel{}{k}^2}\right)^{\omega _2\omega _1}][\left(\frac{s}{\stackrel{}{k}^2}\right)^{\omega _1}+\left(\frac{s}{\stackrel{}{k}^2}\right)^{\omega _1}]L^\mu \},$$
(3.1)
where $`\omega _i=\omega (t_i)`$ and we have chosen $`\stackrel{}{k}^2`$ as the scale of energy, since with this choice the one-gluon contribution to the BFKL kernel is expressed through the RRG-vertex in the simplest way . $`\mathrm{\Gamma }(t_i;\stackrel{}{k}^2)`$ are the helicity conserving gluon-gluon-Reggeon vertices; they depend on $`\stackrel{}{k}^2`$ as on the energy scale . $`R^\mu `$ and $`L^\mu `$ are the right and left RRG-vertices, depending on $`\stackrel{}{q}_1`$ and $`\stackrel{}{q}_2`$. They are real in all physical channels, as well as $`\mathrm{\Gamma }(t_i;\stackrel{}{k}^2)`$.

In the one-loop approximation we have

$$\omega (t)=\omega ^{(1)}(t)=\overline{g}^2\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}^{\mathrm{\hspace{0.25em}2}})^ϵ,$$
(3.2)
$$\mathrm{\Gamma }(t;\stackrel{}{k}^2)=g\left(1+\mathrm{\Gamma }^{(1)}(t;\stackrel{}{k}^2)\right),$$
(3.3)
where

$$\mathrm{\Gamma }^{(1)}(t;\stackrel{}{k}^2)=\overline{g}^2\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}^{\mathrm{\hspace{0.25em}2}})^ϵ[\psi (ϵ)\frac{1}{2}\psi (1)\frac{1}{2}\psi (1ϵ)$$

$$+\frac{9(1+ϵ)^2+2}{4(1+ϵ)(1+2ϵ)(3+2ϵ)}\frac{1}{2}\mathrm{ln}\left(\frac{\stackrel{}{k}^2}{\stackrel{}{q}^{\mathrm{\hspace{0.25em}2}}}\right)].$$
(3.4)
In the same approximation we obtain from (3.1)

$$2g\{𝒜^\mu C^\mu (q_2,q_1)[\mathrm{\Gamma }^{(1)}(t_1;\stackrel{}{k}^2)+\mathrm{\Gamma }^{(1)}(t_2;\stackrel{}{k}^2)+\frac{\omega _1}{2}\mathrm{ln}\left(\frac{s_1(s_1)}{(\stackrel{}{k}^2)^2}\right)+\frac{\omega _2}{2}\mathrm{ln}\left(\frac{s_2(s_2)}{(\stackrel{}{k}^2)^2}\right)$$

$$+\frac{\omega _1+\omega _2}{4}\mathrm{ln}\left(\frac{s(s)(\stackrel{}{k}^2)^2}{s_1(s_1)s_2(s_2)}\right)]\}=R^\mu +L^\mu +(R^\mu L^\mu )\frac{\omega _1\omega _2}{4}\mathrm{ln}(\frac{s_1(s_1)s_2(s_2)}{s(s)(\stackrel{}{k}^2)^2}),$$
(3.5)
where now $`\omega _i\omega ^{(1)}(t_i)`$. Since in the physical region of the $`s`$-channel the relation (2.17) holds, the combinations $`R^\mu +L^\mu `$ and $`R^\mu L^\mu `$ are determined by the real and imaginary parts of the L.H.S. of Eq. (3.5) in this region. Using (2.18) and comparing the imaginary parts in (3.5), we obtain:

$$R^\mu L^\mu =\frac{2g\overline{g}^2}{\omega _1\omega _2}\{[2t_1t_2_3^\mu +(C^\mu (q_2,q_1)4t_1\frac{p_A^\mu }{s_1})$$

$$\times (t_1t_2_3\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ)][AB]\}.$$
(3.6)
We can express $`_3^\mu `$ in term of $`_3`$ by comparing the imaginary parts of (2.25). In this way we obtain from (3.6) the explicitly gauge invariant expression:

$$R^\mu L^\mu =\frac{2g\overline{g}^2}{\omega _1\omega _2}\{C^\mu (q_2,q_1)\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{k}^2)^ϵ+\frac{(\stackrel{}{q}_1\stackrel{}{q}_2)C^\mu (q_2,q_1)+2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}𝒫^\mu }{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_1\stackrel{}{q}_2)^2}$$

$$\times [\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2_3+\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}((\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ(\stackrel{}{q}_2\stackrel{}{k})(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ(\stackrel{}{q}_1\stackrel{}{k})(\stackrel{}{k}^2)^ϵ(\stackrel{}{q}_1\stackrel{}{q}_2))]\}.$$
(3.7)
Evidently, the same result is obtained from the imaginary part of (3.5) if the representation (2.26) for $`𝒜^\mu `$ is used.

The integral $`_3`$ (see (2.14)) cannot be expressed through elementary functions at arbitrary $`D`$. For $`ϵ0`$ we have (see  or Appendix C)

$$_3\frac{1}{ϵ}\left[\frac{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^{ϵ1}}{(\stackrel{}{k}^2)^ϵ}+\frac{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2)^{ϵ1}}{(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ}+\frac{(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2)^{ϵ1}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ}\right].$$
(3.8)
Using this expression we obtain from (3.7)

$$R^\mu L^\mu \frac{4g\overline{g}^2}{\omega _1\omega _2}C^\mu (q_2,q_1)\left(\frac{1}{ϵ}+(\stackrel{}{k}^2)^ϵ\right)$$
(3.9)
for $`ϵ0`$, in agreement with  (see also ).

It’s easy to obtain from (3.7) also the limit $`\stackrel{}{k}0`$ at arbitrary $`D`$. In this limit the main contribution to the integral $`_3`$ (see (2.14)) comes from the region $`\stackrel{}{k}_1\stackrel{}{q}_1\stackrel{}{q}_2`$, so that in the integrand the replacement

$$\frac{1}{\stackrel{}{k}_1^{\mathrm{\hspace{0.25em}2}}}\frac{1}{\stackrel{}{Q}^{\mathrm{\hspace{0.25em}2}}},\stackrel{}{Q}=\frac{\stackrel{}{q}_1+\stackrel{}{q}_2}{2}$$
(3.10)
can be made. After this we have gently

$$_3\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}\frac{(\stackrel{}{k}^2)^{ϵ1}}{\stackrel{}{Q}^{\mathrm{\hspace{0.25em}2}}},$$
(3.11)
and using this result we see that only the first term in curly brackets in (3.7) does contribute, so that in the limit $`\stackrel{}{k}0`$ at arbitrary $`ϵ`$ we obtain

$$R^\mu L^\mu =\frac{2g\overline{g}^2}{\omega _1\omega _2}C^\mu (q_2,q_1)\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{k}^2)^ϵ,$$
(3.12)
in agreement with .

The real parts of (3.5), with account of (2.18), give

$$R^\mu +L^\mu =2g[C^\mu (q_2,q_1)(1+\overline{g}^2\overline{r})$$

$$+\overline{g}^2\{2t_1t_2(_5^\mu +_3^\mu )+2t_1\overline{r}_A\frac{p_A^\mu }{s_1}2t_2\overline{r}_B\frac{p_B^\mu }{s_2}\}],$$
(3.13)
where $`_5^\mu `$ and $`_3^\mu `$ differ from the integrals $`_5`$ and $`_3`$ defined in (2.14) by the factors $`(k_1q_1)_{}^\mu `$ and $`k_{1}^{}{}_{}{}^{\mu }`$, respectively, in the corresponding integrands,

$$\overline{r}=\{\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{2}(_5_3)\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}_{4B}+\frac{\mathrm{\Gamma }^2(ϵ)}{2\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ[\mathrm{ln}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}\right)+\psi (1)+2\psi (ϵ)\psi (1ϵ)2\psi (2ϵ)$$

$$+\frac{1}{2(1+2ϵ)(3+2ϵ)}(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}(11+7ϵ)+2ϵ\frac{\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^2}4\frac{\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^3})]\}+\{AB\},$$
(3.14)
$$\overline{r}_A=\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}(_5_3)+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}_{4B}+\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ\left(\mathrm{ln}\left(\frac{\stackrel{}{k}^2}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)\psi (1)\psi (ϵ)+\psi (1ϵ)+\psi (2ϵ)\right)$$

$$\frac{\mathrm{\Gamma }^2(ϵ)}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\mathrm{\Gamma }(2ϵ)}\frac{1}{2(1+2ϵ)(3+2ϵ)}\{[(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}(11+7ϵ)+ϵ\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^2}$$

$$\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}2\stackrel{}{k}^2}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^3})]+[AB]\},$$
(3.15)
and

$$\overline{r}_B=\overline{r}_A(AB).$$
(3.16)
As well as $`R^\mu L^\mu `$, we can present $`R^\mu +L^\mu `$ in the explicitly gauge invariant form. It can be done using the real part of (2.25) to express $`_5^\mu `$ and $`_3^\mu `$ in terms of $`_5`$, $`_3`$ and $`_{4A,B}`$. We obtain:

$$R^\mu +L^\mu =2g\{C^\mu (q_2,q_1)+\overline{g}^2(\frac{(\stackrel{}{q}_1\stackrel{}{q}_2)C^\mu (q_2,q_1)+2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}𝒫^\mu }{2(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_1\stackrel{}{q}_2)^2)}[\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{2}(_5_3)$$

$$\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_1\stackrel{}{k})_{4B}+\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\frac{(\stackrel{}{k}^2)^ϵ}{2}(\stackrel{}{q}_1\stackrel{}{q}_2)(\psi (ϵ)\psi (1ϵ))+(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ(\stackrel{}{q}_2\stackrel{}{k})(\mathrm{ln}\left(\frac{\stackrel{}{k}^2}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)\psi (1)$$

$$\psi (ϵ)+\psi (1ϵ)+\psi (2ϵ)))]+C^\mu (q_2,q_1)[\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{2}_{4B}+\frac{\mathrm{\Gamma }^2(ϵ)}{2\mathrm{\Gamma }(2ϵ)}(\frac{(\stackrel{}{k}^2)^ϵ}{2}(\psi (ϵ)\psi (1ϵ))$$

$$+(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ[\psi (ϵ)\psi (2ϵ)+\frac{1}{2(1+2ϵ)(3+2ϵ)}(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}(11+7ϵ)+2ϵ\frac{\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^2}$$

$$4\frac{\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^3})])]+𝒫^\mu \frac{\mathrm{\Gamma }^2(ϵ)(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ}{\mathrm{\Gamma }(2ϵ)(1+2ϵ)(3+2ϵ)}[\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}(11+7ϵ)+ϵ\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^2}$$

$$\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}2\stackrel{}{k}^2}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^3}])\overline{g}^2(AB)\}.$$
(3.17)
The same result is obtained from the real part of (3.5) if for $`𝒜^\mu `$ the representation (2.26) is used.

The integrals $`_5`$, $`_3`$ and $`_{4A,B}`$ defined in (2.14) cannot be written in terms of elementary functions at arbitrary $`ϵ`$. For $`ϵ0`$ we have, with accuracy up to terms vanishing in this limit (see Appendix C):

$$_5_3\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\left[\frac{1}{ϵ^2}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}\right)^ϵ\frac{\pi ^2}{6}\right]$$

$$+\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}\left[\frac{1}{ϵ^2}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)^ϵ\frac{\pi ^2}{2}2L\left(1\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)\right]+\frac{1}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}\left[\frac{1}{ϵ^2}\left(\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)^ϵ\frac{\pi ^2}{2}2L\left(1\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)\right],$$
(3.18)
$$_{4B}\frac{1}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\left[\frac{(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ2(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ}{ϵ^2}\frac{\pi ^2}{6}+2L\left(1\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)\right],$$
(3.19)
where

$$L(x)=_0^x\frac{dt}{t}\mathrm{ln}(1t).$$
(3.20)
Using these integrals, we obtain from (3.17)

$$R^\mu +L^\mu =2g\{C^\mu (q_2,q_1)(1+\overline{g}^2[\frac{1}{ϵ^2}\frac{1}{ϵ}\mathrm{ln}(\stackrel{}{k}^2)\frac{1}{2}\mathrm{ln}^2(\stackrel{}{k}^2)\frac{1}{2}\mathrm{ln}^2\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)+\frac{\pi ^2}{2}$$

$$+\frac{\stackrel{}{k}^2}{3}\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^2}+\frac{1}{6}(11\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\frac{4\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^3})\mathrm{ln}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)])$$

$$+𝒫^\mu \overline{g}^2\frac{2}{3}[\frac{\stackrel{}{k}^2}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^2}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2)+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2))+(11\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}$$

$$\frac{\stackrel{}{k}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^3}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}2\stackrel{}{k}^2))\mathrm{ln}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)]\},$$
(3.21)
in agreement with , taking into account the difference in the energy scales (see ) and the charge renormalization (remind that in all formulae above the bare charge $`g`$ is used).

The limit $`\stackrel{}{k}0`$ at arbitrary $`D`$ can be also considered without difficulties. The value of the integral $`I_5`$ in this limit is known . In the $`s`$-channel physical region

$$I_5\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}\frac{(\stackrel{}{k}^2)^{ϵ1}}{s\stackrel{}{Q}^{\mathrm{\hspace{0.25em}2}}}\left[\psi (ϵ)\psi (1ϵ)+i\pi \right].$$
(3.22)
Therefore from (C.2) and (3.11) we obtain

$$_5_3\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}\frac{(\stackrel{}{k}^2)^{ϵ1}}{\stackrel{}{Q}^{\mathrm{\hspace{0.25em}2}}}\left[\psi (1ϵ)\psi (ϵ)\right].$$
(3.23)
The integral $`_{4A}`$ (see (2.14)) is finite at $`\stackrel{}{k}0`$ and can be calculated in this limit by standard methods. We have

$$_{4A}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{Q}^2)^{ϵ1}\left[\psi (ϵ)\psi (2ϵ)\right].$$
(3.24)
Using (3.23) and (3.24), we obtain from (3.17)

$$R^\mu +L^\mu =2gC^\mu (q_2,q_1)\left(1+\overline{g}^2\frac{\mathrm{\Gamma }^2(ϵ)}{2\mathrm{\Gamma }(2ϵ)}(\stackrel{}{k}^2)^ϵ\left[\psi (ϵ)\psi (1ϵ)\right]\right),$$
(3.25)
that agrees with , taking into account the difference in the energy scales (see ).

## 4 The one-gluon contribution to the BFKL kernel

The complicated analytical structure (3.1) of the RRG-vertex is irrelevant for the calculation of the contribution of the one-gluon production to the BFKL kernel at the NLO, where only the real parts of the production amplitudes do contribute, because only these parts interfere with the leading order amplitudes, which are real. Therefore, with the required accuracy, the production amplitude can be presented in the same form as in the leading order:

$$A_{23}=2s\mathrm{\Gamma }_{A^{}A}^{c_1}\left(\frac{s_1}{\stackrel{}{k}^2}\right)^{\omega _1}\frac{1}{t_1}T_{A^{}A}^{c_1}\gamma _{c_1c_2}^d(q_1,q_2)\left(\frac{s_2}{\stackrel{}{k}^2}\right)^{\omega _2}\frac{1}{t_2}T_{B^{}B}^{c_2},$$
(4.1)
where $`\omega _i\omega (t_i)`$, $`\stackrel{}{k}^{\mathrm{\hspace{0.25em}2}}`$ is used as the energy scale and

$$\gamma _{c_1c_2}^d(q_1,q_2)=\frac{1}{2}T_{c_2c_1}^de_\mu ^{}(k)(R^\mu +L^\mu ).$$
(4.2)
The contribution of the one-gluon production to the BFKL kernel for non-forward scattering with momentum transfer $`q`$ and irreducible representation $``$ of the colour group in the $`t`$-channel is

$$𝒦_{RRG}^{G()}(\stackrel{}{q}_1,\stackrel{}{q}_2;\stackrel{}{q})=\frac{1}{(2\pi )^{D1}}\frac{c_1c_1^{}|\widehat{𝒫}_{}|c_2c_2^{}}{2n_{}}\underset{d,\lambda }{}\gamma _{c_1c_2}^d(q_1,q_2)\left(\gamma _{c_1^{}c_2^{}}^d(q_1^{},q_2^{})\right)^{},$$
(4.3)
where $`n_{}`$ is the number of independent states in the representation $``$, the sum is performed over colours and polarizations of the produced gluons and

$$q_1^{}=q_1q,q_2^{}=q_2q.$$
(4.4)
The most interesting representations $``$ are the colour singlet (Pomeron channel) and the antisymmetric colour octet (gluon channel). We have for the singlet case

$$c_1c_1^{}|\widehat{𝒫}_1|c_2c_2^{}=\frac{\delta _{c_1c_1^{}}\delta _{c_2c_2^{}}}{N^21},n_1=1,$$
(4.5)
and for the octet case

$$c_1c_1^{}|\widehat{𝒫}_8|c_2c_2^{}=\frac{f_{c_1c_1^{}c}f_{c_2c_2^{}c}}{N},n_8=N^21,$$
(4.6)
where $`f_{abc}`$ are the structure constants of the colour group. The vertex (4.2) is explicitly invariant under the gauge transformation

$$e^\mu (k)e^\mu (k)+k^\mu \chi ,$$
(4.7)
so that we can use the relation

$$\underset{\lambda }{}e_\mu ^{(\lambda )}(k)e_\nu ^{(\lambda )}(k)=g_{\mu \nu }.$$
(4.8)
Substituting (4.2) in (4.3), using (4.8) for the sum over polarizations and

$$\delta _{c_1c_1^{}}\delta _{c_2c_2^{}}T_{c_1c_2}^d\left(T_{c_1^{}c_2^{}}^d\right)^{}=N(N^21),f_{c_1c_1^{}c}f_{c_2c_2^{}c}T_{c_1c_2}^d\left(T_{c_1^{}c_2^{}}^d\right)^{}=\frac{N^2(N^21)}{2}$$
(4.9)
for the sum over colour indices, we obtain with the NLO accuracy:

$$𝒦_{RRG}^{G()}(\stackrel{}{q}_1,\stackrel{}{q}_2;\stackrel{}{q})=\frac{c_{}}{8(2\pi )^{D1}}2g[C_\mu (q_2^{},q_1^{})(R+L)^\mu +C_\mu (q_2,q_1)(R^{}+L^{})^\mu $$

$$2gC_\mu (q_2,q_1)C^\mu (q_2^{},q_1^{})],$$
(4.10)
where for the singlet ($`=1`$) and octet ($`=8`$) cases

$$c_1=N,c_8=\frac{N}{2},$$
(4.11)
$`C_\mu (q_2,q_1)`$ is defined in (2.6), the sum $`(R+L)^\mu `$ is presented in two different forms in (3.13)-(3.16) and (3.17), $`(R^{}+L^{})^\mu `$ is obtained from $`(R+L)^\mu `$ by the substitution

$$\stackrel{}{q}_1\stackrel{}{q}_1^{}=\stackrel{}{q}_1\stackrel{}{q},\stackrel{}{q}_2\stackrel{}{q}_2^{}=\stackrel{}{q}_2\stackrel{}{q}.$$
(4.12)
The convolution over the vector indices in (4.10) is easily performed with the help of the relations

$$C_\mu (q_2,q_1)C^\mu (q_2^{},q_1^{})=2\left(\stackrel{}{q}^{\mathrm{\hspace{0.25em}2}}\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}\right),$$

$$C_\mu (q_2,q_1)𝒫^\mu =\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}1.$$
(4.13)
In general case the convolution does not lead to noticeable simplifications, so that it has no sense to rewrite (4.10) in unfolded form. But at $`\stackrel{}{q}=0`$ there is a huge simplification due to the equality

$$C_\mu (q_2,q_1)\left[(\stackrel{}{q}_1\stackrel{}{q}_2)C^\mu (q_2,q_1)+2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}𝒫^\mu \right]=0.$$
(4.14)
For this case we obtain

$$𝒦_{RRG}^{G()}(\stackrel{}{q}_1,\stackrel{}{q}_2;\stackrel{}{0})=\frac{g^2c_{}}{(2\pi )^{D1}}\{\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}[1+2\overline{g}^2(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}_{4B}+\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\frac{(\stackrel{}{k}^2)^ϵ}{2}(\psi (ϵ)\psi (1ϵ))$$

$$+(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ(\psi (ϵ)\psi (2ϵ))))]\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}\frac{\overline{g}^2(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ+1}}{(1+2ϵ)(3+2ϵ)}[12\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}$$

$$+\frac{\stackrel{}{k}^2\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^3}(3(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})2\stackrel{}{k}^2)+ϵ(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}7\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})}\frac{\stackrel{}{k}^2}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^2}(2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2))]\}$$

$$+\frac{g^2c_{}}{(2\pi )^{D1}}\{AB\}.$$
(4.15)
Acknowledgment: V.S.F. thanks the Dipartimento di Fisica della Università della Calabria and the Istituto Nazionale di Fisica Nucleare - Gruppo collegato di Cosenza for the warm hospitality while part of this work was done.

## A Appendix

In this Appendix we give the expressions of the integrals relevant for this paper, among those having the general form of Eqs. (2.13). The integrals with two or three denominators were already calculated in  at arbitrary $`D`$. Those with four and five denominators were calculated in Ref.  in the $`ϵ`$-expansion. Unfortunately, it is not possible in general to express them in terms of elementary functions at arbitrary $`D`$. It is possible, however, to perform the integration over the longitudinal variables in the Sudakov decomposition for the integration variable $`p`$ and to express the result in terms of $`(D2)`$-dimensional integrals over $`\stackrel{}{p}`$. As an illustration of the method, we will show in some details the steps involved in the calculation of $`I_4^{(1)}`$ for which it turns out that, in the multi-Regge kinematics (2.1), the complete integration can be performed. For the remaining integrals with four of five denominators, we will merely list the results.

Let us start from the integrals with two or three denominators:

$$I_2(q)=\frac{1}{i}d^Dp\frac{1}{(p^2+i\epsilon )[(p+q)^2+i\epsilon ]}=\frac{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}{2(1+2ϵ)}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}^{\mathrm{\hspace{0.25em}2}})^ϵ,$$
(A.1)
$$I_2^{\mu \nu }(q)=\frac{1}{i}d^Dp\frac{p^\mu p^\nu }{(p^2+i\epsilon )[(p+q)^2+i\epsilon ]}=\frac{1}{4(3+2ϵ)}\left[2(2+ϵ)q^\mu q^\nu q^2g^{\mu \nu }\right]I_2(q),$$
(A.2)
$$I_3=\frac{1}{i}d^Dp\frac{1}{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+k)^2+i\epsilon ]}=\frac{1+2ϵ}{ϵ}\frac{I_2(q_2)I_2(q_1)}{t_1t_2},$$
(A.3)
$$I_3^\mu =\frac{1}{i}d^Dp\frac{p^\mu }{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+k)^2+i\epsilon ]}=q_1^\mu \frac{I_2(q_1)I_2(q_2)}{t_1t_2}$$

$$+\frac{k^\mu }{(t_1t_2)^2}\left[\frac{t_1}{ϵ}I_2(q_1)+\left(t_2\frac{1+ϵ}{ϵ}t_1\right)I_2(q_2)\right],$$
(A.4)
$$I_3^{\mu \nu }=\frac{1}{i}d^Dp\frac{p^\mu p^\nu }{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+k)^2+i\epsilon ]}=\frac{g^{\mu \nu }}{4(1+ϵ)}\frac{t_1I_2(q_1)t_2I_2(q_2)}{t_1t_2}$$

$$\frac{q_1^\mu q_1^\nu }{2}\frac{I_2(q_1)I_2(q_2)}{t_1t_2}+\frac{k^\mu q_1^\nu +k^\nu q_1^\mu }{2(1+ϵ)}\frac{t_1I_2(q_1)+[(1+ϵ)t_1ϵt_2)]I_2(q_2)}{(t_1t_2)^2}$$

$$+\frac{k^\mu k^\nu }{(t_1t_2)^3}\left[\frac{t_1^2I_2(q_1)}{ϵ(1+ϵ)}+\left(\frac{t_1^2}{ϵ}\frac{t_1t_2}{1+ϵ}+\frac{(t_1t_2)^2}{2}\right)I_2(q_2)\right],$$
(A.5)
$$I_3^{\mu \nu \rho }=\frac{1}{i}d^Dp\frac{p^\mu p^\nu p^\rho }{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+k)^2+i\epsilon ]}$$

$$=\left(g^{\mu \nu }q_1^\rho +g^{\mu \rho }q_1^\nu +g^{\nu \rho }q_1^\mu \right)\frac{t_2I_2(q_2)t_1I_2(q_1)}{4(3+2ϵ)(t_1t_2)}$$

$$\left(g^{\mu \nu }k^\rho +g^{\mu \rho }k^\nu +g^{\nu \rho }k^\mu \right)\frac{t_1^2I_2(q_1)[(2+ϵ)t_2(t_1t_2)+t_2^2]I_2(q_2)}{4(1+ϵ)(3+2ϵ)(t_1t_2)^2}$$

$$+q_1^\mu q_1^\nu q_1^\rho \frac{(2+ϵ)[I_2(q_1)I_2(q_2)]}{2(3+2ϵ)(t_1t_2)}$$

$$+\left(q_1^\mu q_1^\nu k^\rho +q_1^\mu k^\nu q_1^\rho +k^\mu q_1^\nu q_1^\rho \right)\frac{t_1I_2(q_1)[t_2+(1+ϵ)(t_1t_2)]I_2(q_2)}{2(3+2ϵ)(t_1t_2)^2}$$

$$+\left(q_1^\mu k^\nu k^\rho +k^\mu q_1^\nu k^\rho +k^\mu k^\nu q_1^\rho \right)\frac{2t_1^2I_2(q_1)[2t_2^2+(2+ϵ)(t_1t_2)(2t_2+(1+ϵ)(t_1t_2))]I_2(q_2)}{2(1+ϵ)(3+2ϵ)(t_1t_2)^3}$$

$$+k^\mu k^\nu k^\rho \{\frac{3[t_1^3I_2(q_1)t_2^3I_2(q_2)]}{ϵ(1+ϵ)(3+2ϵ)(t_1t_2)^4}$$

$$[\frac{3(3+ϵ)t_2^2}{ϵ(1+ϵ)(3+2ϵ)(t_1t_2)^3}+\frac{3t_2}{(t_1t_2)^2}(\frac{1}{2(3+2ϵ)}+\frac{1}{ϵ(1+ϵ)})$$

$$+\frac{1}{t_1t_2}(\frac{1}{4(3+2ϵ)}+\frac{1}{ϵ}+\frac{1}{4})]I_2(q_2)\},$$
(A.6)
$$I_{3A}=\frac{1}{i}d^Dp\frac{1}{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+p_A)^2+i\epsilon ]}=\frac{1+2ϵ}{ϵ}\frac{I_2(q_1)}{t_1},$$
(A.7)
$$I_{3A}^\mu =\frac{1}{i}d^Dp\frac{p^\mu }{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+p_A)^2+i\epsilon ]}=\left(q_1^\mu +\frac{p_A^\mu }{ϵ}\right)\frac{I_2(q_1)}{t_1}.$$
(A.8)
Let us consider now integrals with four or five denominators. As anticipated above, we will show in some detail how to calculate $`I_4^{(1)}`$ and simply give the final results for the remaining integrals.

The integral under consideration is

$$I_4^{(1)}=\frac{1}{i}d^Dp\frac{1}{(p^2+i\epsilon )[(p+q_1)^2+i\epsilon ][(p+p_A)^2+i\epsilon ][(pp_B)^2+i\epsilon ]}.$$
(A.9)
The Sudakov decompositions for the integration momentum $`p`$ and for $`q_1`$ are

$$p=\beta p_A+\alpha p_B+p_{},q_1p_Ap_A^{}=\frac{s_2}{s}p_A\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{s}p_B+q_{1}^{}{}_{}{}^{}.$$
(A.10)
After changing the integration variables to $`\alpha `$, $`\beta `$ and $`p_{}`$ and using

$$d^Dp=\frac{s}{2}d\alpha d\beta d^{D2}p,$$
(A.11)
we perform the integration over the variable $`\alpha `$ by the method of residues. This leads in the multi-Regge kinematics to the following result

$$I_4^{(1)}=\pi _0^1𝑑\beta \beta ^2\frac{d^{D2}p}{\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{p}+\beta \stackrel{}{q}_1)^2[\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}+\beta (1\beta )(si\epsilon )]}$$

$$\pi \frac{s_2}{s}_0^1𝑑\beta \beta ^2\frac{d^{D2}p}{\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}[(\stackrel{}{p}+\beta \stackrel{}{q}_1)^2+\beta (1\beta )\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}][(\stackrel{}{p}+\beta \stackrel{}{q}_1)^2+\beta (1\beta )(s_2i\epsilon )]}.$$
(A.12)
The first integral in the R.H.S. of the above equation can be manipulated as follows:

$$\left[I_4^{(1)}\right]_1\pi _0^1𝑑\beta \beta ^2\frac{d^{D2}p}{\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{p}+\beta \stackrel{}{q}_1)^2[\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}+\beta (1\beta )(si\epsilon )]}$$

$$=\pi \frac{d^{D2}p}{\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{p}+\stackrel{}{q}_1)^2}_0^1𝑑\beta \frac{(1\beta )^{D5}}{(1\beta )\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}+\beta (si\epsilon )}$$

$$\pi \frac{d^{D2}p}{\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{p}+\stackrel{}{q}_1)^2}_0^1𝑑\beta \frac{(1\beta )^{D5}}{\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}+\beta (si\epsilon )},$$
(A.13)
where the first equality follows from the change of variables $`\stackrel{}{p}\beta \stackrel{}{p}`$ and $`\beta 1\beta `$, while the last approximated equality holds since in the $`s\mathrm{}`$ limit $`s+\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}=s`$. Then, the integral over $`\beta `$ can be performed:

$$_0^1𝑑\beta \frac{(1\beta )^{D5}}{\delta +\beta }=_0^1𝑑\beta \frac{(1\beta )^{D5}1}{\delta +\beta }+_0^1\frac{d\beta }{\delta +\beta }$$

$$_0^1𝑑\beta \frac{(1\beta )^{D5}1}{\beta }+_0^1\frac{d\beta }{\delta +\beta }\psi (1)\psi (D4)+\mathrm{ln}\left(\frac{1}{\delta }\right),$$
(A.14)
where $`\delta \stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}/(si\epsilon )`$ is a quantity tending to zero and $`\psi (x)`$ is the logarithmic derivative of $`\mathrm{\Gamma }(x)`$. After performing the integration over $`p_{}`$, we obtain for $`\left[I_4^{(1)}\right]_1`$ the following result:

$$\left[I_4^{(1)}\right]_1=\frac{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}{st_1}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ\left[\mathrm{ln}\left(\frac{si\epsilon }{t_1}\right)+\psi (1ϵ)\psi (ϵ)\right].$$
(A.15)
In the last expression, the terms $`i\epsilon `$ which appears in the argument of the logarithm fixes the prescription for the analytic continuation of this function to the region of positive $`s`$. In the following, this term will always be omitted in the final results, but it should be understood to accompany $`s`$, $`s_1`$ or $`s_2`$ every time these quantities appear in the argument of the logarithm. The second integral in the R.H.S. of Eq. (A.12) can be decomposed as follows:

$$\left[I_4^{(1)}\right]_2=\pi \frac{s_2}{s}_0^1𝑑\beta \beta ^2\frac{d^{D2}p}{(\stackrel{}{p}+\beta \stackrel{}{q}_1)^2[\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}+\beta (1\beta )\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}][\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}+\beta (1\beta )(s_2i\epsilon )]}$$

$$=\pi \frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{s}_0^1𝑑\beta \beta ^2\frac{d^{D2}p}{\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{p}+\beta \stackrel{}{q}_1)^2[\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}+\beta (1\beta )\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}]}$$

$$\pi \frac{s_2}{s}_0^1𝑑\beta \beta ^2\frac{d^{D2}p}{\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}(\stackrel{}{p}+\beta \stackrel{}{q}_1)^2[\stackrel{}{p}^{\mathrm{\hspace{0.25em}2}}+\beta (1\beta )(s_2i\epsilon )]},$$
(A.16)
where it has been used that in the multi-Regge kinematics $`s_2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}`$. The second integral in the R.H.S. of the above equation has the same structure as $`\left[I_4^{(1)}\right]_1`$ (see Eq. (A.13)) and does not need to be calculated. The first one can be easily evaluated by the usual Feynman parametrization technique. The final result for $`\left[I_4^{(1)}\right]_2`$ is then

$$\left[I_4^{(1)}\right]_2=\frac{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}{st_1}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ\left[\mathrm{ln}\left(\frac{t_1}{s_2}\right)+\psi (1)\psi (1ϵ)\right].$$
(A.17)
Combining Eqs. (A.15) and (A.17), we obtain the final result for $`I_4^{(1)}`$:

$$I_4^{(1)}=\frac{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}{st_1}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ\left[\mathrm{ln}\left(\frac{s}{s_2}\right)+\psi (1)\psi (ϵ)\right].$$
(A.18)
The same procedure described for the calculation of $`I_4^{(1)}`$ can be applied for the integrals $`I_{4A}`$, $`I_5`$ and $`I_5^\mu `$, except that for them the integrations in the transverse space cannot be performed in complete way. We list here the final results:

$$I_{4A}=\frac{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}{s_1}\left[\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ1}\left(\mathrm{ln}\left(\frac{s_1}{t_1}\right)+\psi (1ϵ)2\psi (ϵ)+\psi (2ϵ)\right)+_{4A}\right],$$
(A.19)
$$I_5=\frac{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}{s}\left[\mathrm{ln}\left(\frac{(s)\stackrel{}{k}^2}{(s_1)(s_2)}\right)_3+_3_5\right],$$
(A.20)
$$I_5^\mu =q_1^\mu I_5p_A^\mu \frac{I_{4A}}{s}+p_B^\mu \frac{I_{4B}}{s}+\frac{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}{s}\left[\mathrm{ln}\left(\frac{(s)\stackrel{}{k}^2}{(s_1)(s_2)}\right)_3^\mu _3^\mu _5^\mu \right],$$
(A.21)
where

$$_{4A}=_0^1\frac{dx}{1x}\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}[\frac{x}{[(1x)\stackrel{}{k}_1^2+x(\stackrel{}{k}_1+\stackrel{}{q}_2)^2](\stackrel{}{k}_1x\stackrel{}{k})^2}\frac{1}{(\stackrel{}{k}_1+\stackrel{}{q}_2)^2(\stackrel{}{k}_1\stackrel{}{k})^2}],$$

$$_5=_0^1\frac{dx}{1x}\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{1}{\stackrel{}{k}_1^2[(1x)\stackrel{}{k}_1^2+x(\stackrel{}{k}_1\stackrel{}{q}_1)^2]}\left(\frac{x^2}{(\stackrel{}{k}_1x\stackrel{}{k})^2}\frac{1}{(\stackrel{}{k}_1\stackrel{}{k})^2}\right),$$

$$_5^\mu =_0^1\frac{dx}{1x}\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{(k_1q_1)_{}^\mu }{\stackrel{}{k}_1^2[(1x)\stackrel{}{k}_1^2+x(\stackrel{}{k}_1\stackrel{}{q}_1)^2]}\left(\frac{x^2}{(\stackrel{}{k}_1x\stackrel{}{k})^2}\frac{1}{(\stackrel{}{k}_1\stackrel{}{k})^2}\right),$$

$$_3=\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{1}{\stackrel{}{k}_1^2(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2}\mathrm{ln}\left(\frac{(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2}{\stackrel{}{k}^2\stackrel{}{k}_1^2}\right),$$

$$_3^\mu =\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{k_{1}^{}{}_{}{}^{\mu }}{\stackrel{}{k}_1^2(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2}\mathrm{ln}\left(\frac{(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2}{\stackrel{}{k}^2\stackrel{}{k}_1^2}\right),$$

$$_3=\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{1}{\stackrel{}{k}_1^2(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2},$$

$$_3^\mu =\frac{d^{2+2ϵ}k_1}{\pi ^{1+ϵ}\mathrm{\Gamma }(1ϵ)}\frac{k_{1}^{}{}_{}{}^{\mu }}{\stackrel{}{k}_1^2(\stackrel{}{k}_1\stackrel{}{q}_1)^2(\stackrel{}{k}_1\stackrel{}{q}_2)^2},$$
(A.22)
$`I_{4B}`$ in the expression for $`I_5^\mu `$ (Eq. (A.21)) is obtained from $`I_{4A}`$ by the replacements $`p_Ap_B`$ and $`q_1q_2`$. In the multi-Regge kinematics these replacements imply that $`I_{4B}=I_{4A}(s_1s_2,\stackrel{}{q}_1\stackrel{}{q}_2)`$.

## B Appendix

In this Appendix we derive the relations (2.21) and (2.22). Let us consider first Eq. (2.21). From the definition of $`_5^\mu `$ given in the first of Eqs. (2.16) and from the expression for $`I_5^\mu `$ given in Eq. (A.21), we obtain

$$_5^\mu =\frac{s}{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}\left[q_1^\mu I_5p_A^\mu \frac{I_{4A}}{s}+p_B^\mu \frac{I_{4B}}{s}I_5^\mu \right]+\frac{1}{2}\mathrm{ln}\left(\frac{s(s_1)(s_2)}{(s)s_1s_2}\right)_3^\mu .$$
(B.1)
If we contract both sides of the above equation with $`2k_\mu `$, we get

$$2k_\mu _5^\mu =\frac{s}{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}\left[(t_1t_2)I_5\frac{s_1}{s}I_{4A}+\frac{s_2}{s}I_{4B}+I_4^{(1)}I_4^{(2)}\right]$$

$$+\frac{1}{2}\mathrm{ln}\left(\frac{s(s_1)(s_2)}{(s)s_1s_2}\right)\left[\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}\left((\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^{ϵ1}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ1}\right)+(t_1t_2)_3\right],$$
(B.2)
where we have used the relations

$$k_\mu I_5^\mu =\frac{I_4^{(2)}I_4^{(1)}}{2},k_\mu _3^\mu =\frac{1}{2}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}\left((\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^{ϵ1}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ1}\right)+\frac{t_1t_2}{2}_3.$$
(B.3)
$`I_4^{(2)}`$ in Eq. (B.2) is obtained from $`I_4^{(1)}`$ by the replacements $`p_Ap_B`$ and $`q_1q_2`$. In the multi-Regge kinematics these replacements imply $`I_4^{(2)}=I_4^{(1)}(t_1t_2,s_2s_1)`$. Using Eqs. (A.18)-(A.20) and the third of Eqs. (2.16) in the R.H.S. of Eq. (B.2), it is easy to obtain the relation (2.21).

Let us consider now the relation (2.22). Starting again from Eq. (B.1) and contracting both sides with $`2(q_1+q_2)_\mu `$, we get

$$2(q_1+q_2)_\mu _5^\mu =\frac{s}{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}\left[(t_1+t_2)I_5+\frac{s_1}{s}I_{4A}+\frac{s_2}{s}I_{4B}I_4^{(1)}I_4^{(2)}+2I_4\right]$$

$$+\frac{1}{2}\mathrm{ln}\left(\frac{s(s_1)(s_2)}{(s)s_1s_2}\right)\left[\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}\left((\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ1}+(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^{ϵ1}2(\stackrel{}{k}^2)^{ϵ1}\right)+(t_1+t_2)_3\right],$$
(B.4)
where we have used the relations

$$(q_1+q_2)_\mu I_5^\mu =\frac{I_4^{(1)}+I_4^{(2)}}{2}I_4+t_1I_5,$$

$$(q_1+q_2)_\mu _3^\mu =\frac{1}{2}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}\left((\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ1}+(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^{ϵ1}2(\stackrel{}{k}^2)^{ϵ1}\right)+\frac{t_1+t_2}{2}_3.$$
(B.5)
In the previous expressions a new integral appeared

$$I_4=\frac{1}{i}d^Dp\frac{1}{[(p+q_1)^2+i\epsilon ][(p+q_2)^2+i\epsilon ][(p+p_A)^2+i\epsilon ][(pp_B)^2+i\epsilon ]}.$$
(B.6)
It can be evaluated using the procedure illustrated in Appendix A for the integral $`I_4^{(1)}`$. The result is

$$I_4=\frac{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}{s}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{k}^2)^{ϵ1}\left[\mathrm{ln}\left(\frac{(s)\stackrel{}{k}^2}{(s_1)(s_2)}\right)+\psi (ϵ)\psi (1ϵ)\right].$$
(B.7)
Using Eqs. (A.18)-(A.20), (B.7) and the third of Eqs. (2.16) in the R.H.S. of Eq. (B.4), it is easy to obtain the relation (2.22).

## C Appendix

In this Appendix we consider the $`ϵ0`$ limit for the integrals $`_3`$, $`_{4A}`$ and for the combination $`_5_3`$. The definitions of these integrals are given in (2.14). The integral $`_{4B}`$, whose expression for $`ϵ0`$ is given in Eq. (3.19), can be obtained from $`_{4A}`$ by the replacement $`\stackrel{}{q}_1\stackrel{}{q}_2`$.

Let us start from $`_3`$. This integral was calculated in the Appendix II of Ref.  to the order $`ϵ^0`$. We report here its expression

$$_3=\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\left[\frac{1}{ϵ}+\mathrm{ln}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}\right)\right]+\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}\left[\frac{1}{ϵ}+\mathrm{ln}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)\right]+\frac{1}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}\left[\frac{1}{ϵ}+\mathrm{ln}\left(\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)\right]$$

$$\frac{1}{ϵ}\left[\frac{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^{ϵ1}}{(\stackrel{}{k}^2)^ϵ}+\frac{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2)^{ϵ1}}{(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ}+\frac{(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2)^{ϵ1}}{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ}\right],$$
(C.1)
where the last approximated equality holds with accuracy $`O(ϵ)`$.

The combination of integrals $`_5_3`$ can be expressed in terms of the $`I_5`$ and $`_3`$, according to Eq. (A.20):

$$_5_3=\frac{s}{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}I_5+\mathrm{ln}\left(\frac{(s)\stackrel{}{k}^2}{(s_1)(s_2)}\right)_3.$$
(C.2)
The integral $`I_5`$ was calculated in the Appendix III of Ref.  to the order $`ϵ^0`$:

$``$ $`{\displaystyle \frac{s}{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}}I_5=`$
$``$ $`{\displaystyle \frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}}\left[{\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{1}{ϵ}}\mathrm{ln}\left({\displaystyle \frac{(s)\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(s_1)(s_2)}}\right)+{\displaystyle \frac{1}{2}}\mathrm{ln}^2\left({\displaystyle \frac{(s)\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(s_1)(s_2)}}\right)+{\displaystyle \frac{\pi ^2}{3}}\right]`$
$`+`$ $`{\displaystyle \frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}}\left[{\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{1}{ϵ}}\mathrm{ln}\left({\displaystyle \frac{(s_1)(s_2)\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{(s)\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}}\right)+{\displaystyle \frac{1}{2}}\mathrm{ln}^2\left({\displaystyle \frac{(s_1)(s_2)\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{(s)\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}}\right)2L\left(1{\displaystyle \frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}}\right)\right]`$
$`+`$ $`{\displaystyle \frac{1}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}}\left[{\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{1}{ϵ}}\mathrm{ln}\left({\displaystyle \frac{(s_1)(s_2)\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(s)\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}}\right)+{\displaystyle \frac{1}{2}}\mathrm{ln}^2\left({\displaystyle \frac{(s_1)(s_2)\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{(s)\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}}\right)2L\left(1{\displaystyle \frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}}\right)\right],`$
where

$$L(x)=_0^x\frac{dt}{t}\mathrm{ln}(1t).$$
(C.4)
Using in (C.2) the expression for $`_3`$ given in the first line of Eq. (C.1), we get

$`_5_3=`$ $``$ $`{\displaystyle \frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}}\left[{\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{1}{ϵ}}\mathrm{ln}\left({\displaystyle \frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}}\right)+{\displaystyle \frac{1}{2}}\mathrm{ln}^2\left({\displaystyle \frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}}\right){\displaystyle \frac{\pi ^2}{6}}\right]`$
$`+`$ $`{\displaystyle \frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}}\left[{\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{1}{ϵ}}\mathrm{ln}\left({\displaystyle \frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}}\right)+{\displaystyle \frac{1}{2}}\mathrm{ln}^2\left({\displaystyle \frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}}\right){\displaystyle \frac{\pi ^2}{2}}2L\left(1{\displaystyle \frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}}\right)\right]`$
$`+`$ $`{\displaystyle \frac{1}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}}\left[{\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{1}{ϵ}}\mathrm{ln}\left({\displaystyle \frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}}\right)+{\displaystyle \frac{1}{2}}\mathrm{ln}^2\left({\displaystyle \frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}}\right){\displaystyle \frac{\pi ^2}{2}}2L\left(1{\displaystyle \frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}}\right)\right]`$
$$\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\left[\frac{1}{ϵ^2}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{k}^2}\right)^ϵ\frac{\pi ^2}{6}\right]+\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}\left[\frac{1}{ϵ^2}\left(\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)^ϵ\frac{\pi ^2}{2}2L\left(1\frac{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)\right]$$

$$+\frac{1}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}\left[\frac{1}{ϵ^2}\left(\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\stackrel{}{k}^2}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)^ϵ\frac{\pi ^2}{2}2L\left(1\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)\right],$$
(C.5)
where the last approximated equality holds with accuracy $`O(ϵ)`$ and we have used $`\mathrm{ln}(s)=\mathrm{ln}si\pi `$ and the analogous relations for $`s_1`$ and $`s_2`$, valid in the $`s`$-channel physical region.

Finally, let us consider $`_{4A}`$. According to (A.19), it can be written in the following way:

$$_{4A}=\frac{s_1}{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}I_{4A}\frac{\mathrm{\Gamma }^2(ϵ)}{\mathrm{\Gamma }(2ϵ)}(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^{ϵ1}\left[\mathrm{ln}\left(\frac{(s_1)}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)+\psi (1ϵ)2\psi (ϵ)+\psi (2ϵ)\right].$$
(C.6)
The integral $`I_{4A}`$ was calculated in the Appendix III of Ref.  to the order $`ϵ^0`$:

$$\frac{s_1}{\pi ^{2+ϵ}\mathrm{\Gamma }(1ϵ)}I_{4A}=\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}[\frac{2}{ϵ^2}+\frac{2}{ϵ}\mathrm{ln}\left(\frac{(s_1)\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}\right)\mathrm{ln}^2\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}+2\mathrm{ln}\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\mathrm{ln}(s_1)$$

$$\pi ^2+2L(1\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}})].$$
(C.7)
Using this expression in the R.H.S. of Eq. (C.6) together with the expansion to the order $`ϵ^0`$ of the second term in the R.H.S. of the same equation, we get

$$_{4A}=\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\left[\frac{1}{ϵ^2}+\frac{1}{ϵ}\mathrm{ln}\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}+\frac{1}{2}\mathrm{ln}^2\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}\frac{2}{ϵ}\mathrm{ln}\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\mathrm{ln}^2\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}\frac{\pi ^2}{6}+2L\left(1\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)\right]$$

$$\frac{1}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\left[\frac{(\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}})^ϵ2(\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}})^ϵ}{ϵ^2}\frac{\pi ^2}{6}+2L\left(1\frac{\stackrel{}{q}_2^{\mathrm{\hspace{0.25em}2}}}{\stackrel{}{q}_1^{\mathrm{\hspace{0.25em}2}}}\right)\right],$$
(C.8)
where the last approximated equality holds with accuracy $`O(ϵ)`$.