format: text

text:
# Cavendish-HEP-02/04DESY 02-114hep-ph/0208130 NLO Corrections to the Photon Impact Factor: Combining Real and Virtual Corrections

## I Introduction

This paper represents the third part of our investigation of the photon impact factor in next-to-leading order QCD. First results of the calculation of NLO corrections to the photon impact factor have been published in  and . In  we have computed the QCD NLO corrections to the process $`\gamma ^{}+q(q\overline{q})+q`$ which lead to the virtual corrections of the photon impact factor. The second paper, , contains the process $`\gamma ^{}+q(q\overline{q}g)+q`$, leading to the real corrections of the photon impact factor. Results are given for the helicity-summed squared matrix elements, with the longitudinally polarized photon in the initial state. In the present paper we will complete this part of our calculation by listing also the results for the transverse polarization of the photon.

The main purpose of the present paper is the combination of singular pieces of virtual and real corrections. To begin with the former ones, in  we have listed the results for the virtual corrections which contain both the finite and the singular pieces. As the first step we therefore have to separate the sum of all singular $`1/ϵ`$ terms from the finite part. In the next step we turn to the real corrections. Infrared singularities will come from those kinematic regions of the $`q\overline{q}g`$ system where the gluon is either soft or collinear with the quark or with the antiquark. We compute the soft and the collinear approximations of our matrix elements, and, by subtracting these divergent contributions from the full matrix elements, we form finite combinations. Finally, by combining the singular pieces of the real corrections with those coming from the virtual corrections, we obtain further finite terms which have the form of the Born approximation and can be computed analytically. As the main result of this paper, we present finite expressions for the $`q\overline{q}g`$ contributions to the photon impact factor. As the fourth and remaining part of our program we will be left with the task of evaluating numerically the finite integrals obtained in this paper, as well as the finite pieces from the virtual corrections.

This paper will be organized as follows. After the definition of the impact factor (section II) we first (section III) list the infrared singular pieces of the virtual corrections. We then turn to the real corrections and present the complete results for the real corrections, both for the longitudinal and for the transverse photon (section IV). In the following section V we compute the collinear and the soft limits of the real corrections. We discuss in detail the definition of the impact factor correction in section VI, where the subtraction of the central region, i.e. the leading $`\mathrm{log}s`$ term, is discussed. In section VII we complete our program, by combining the singular pieces of the real corrections with those of the virtual corrections and by demonstrating their infrared finiteness. In a final section we summarize our results and give a brief outline of the remaining part of our program.

## II Definition of the Impact Factor

We consider the elastic scattering of two particles $`A`$ and $`B`$ with momenta $`p_A`$ and $`p_B`$ in the Regge limit $`s\mathrm{}`$, where the momentum transfer $`t=q^2`$ is kept fixed. In leading order in $`\alpha _s`$ we write the scattering amplitude $`T_{AB}`$ in the form:

$$T_{AB}^{(0)}(s,t=q^2)=is\frac{d^{D2}𝒓}{(2\pi )^{D2}}\stackrel{~}{\mathrm{\Phi }}_A^{(0)}\frac{1}{𝒓^2}\frac{1}{(𝒒𝒓)^2}\stackrel{~}{\mathrm{\Phi }}_B^{(0)}.$$
(1)
The impact factors $`\mathrm{\Phi }_A^{(0)}`$ and $`\mathrm{\Phi }_B^{(0)}`$ for the external particles (quarks, gluons, or photons) are of order $`\alpha _s`$. The exchanged gluons carry transverse momenta $`𝒓`$ and $`(𝒒𝒓)`$ respectively. Taking the cut amplitude we define the total cross section for the scattering of $`A`$ and $`B`$ via the optical theorem as

$$\sigma _{AB}^{(0)}=\frac{1}{s}\mathrm{Im}T_{AB}^{(0)}(s,t=0)=\frac{d^{D2}𝒓}{(2\pi )^{D2}}\mathrm{\Phi }_A^{(0)}\frac{1}{𝒓^4}\mathrm{\Phi }_B^{(0)}.$$
(2)
Here we have dropped the tilde symbol on the impact factors to denote the additional $`s`$-channel cut. This equation can be used as the definition of the leading order impact factors $`\mathrm{\Phi }_A^{(0)}`$ and $`\mathrm{\Phi }_B^{(0)}`$. Note that we have summed over the colour indices of the $`t`$-channel gluons implicitly. The total cross section can also be written as

$$\sigma _{AB}^{(0)}=\frac{1}{2s}|_{AB}^{(0)}|^2𝑑\varphi ,$$
(3)
where $`1/2s`$ denotes the flux factor (in the high energy limit), $`_{AB}^{(0)}`$ the matrix element for the scattering process $`A+BA^{}+B^{}`$ with the exchange of a gluon, and $`d\varphi `$ stands for the phase space of the final state $`A^{}+B^{}`$. For our case of interest, $`\gamma ^{}q`$ (or $`\gamma ^{}\gamma ^{}`$) scattering, $`A^{}`$ denotes a $`q\overline{q}`$ pair, $`B^{}`$ a quark (or a $`q\overline{q}`$ pair).

In the high energy limit, the scattering amplitude $`_{AB}`$ for the scattering $`A+BA^{}+B^{}`$ (with the invariant masses of the particles or systems $`A^{}`$ and $`B^{}`$ being finite) (Fig. 1 without the extra produced gluon) is described by the exchange of a reggeized gluon and can be written in the following form:

$$_{AB}=\frac{s}{t}\mathrm{\Gamma }_{AA^{}}^a\left[\left(\frac{s}{t}\right)^{\omega (t)}+\left(\frac{s}{t}\right)^{\omega (t)}\right]\mathrm{\Gamma }_{BB^{}}^a.$$
(4)
Here, $`\omega (t)`$ stands for the Regge trajectory of the gluon, and $`\mathrm{\Gamma }_{AA^{}}^a`$ and $`\mathrm{\Gamma }_{BB^{}}^a`$ are the particle-particle-reggeon vertices. The index $`a`$ denotes the colour of the exchanged reggeized gluon. Eq.(4) exhibits the factorization property of Regge theory: when written in complex angular plane, the residue of the gluon Regge pole can be written as a product of two vertex functions, one for the incoming particle $`AA^{}`$, the other one for $`BB^{}`$. In the total cross section formula (3) this factorization property leads to the factorization of impact factors. Since (4) is valid not only in leading order, also the validity of the impact factor representation is rather general. However, when going beyond the leading order, we have to replace the $`t`$-channel propagators of the exchanged gluons by reggeon propagators. In leading order, we put $`\omega (t)`$ in (4) equal to zero and obtain

$$_{AB}^{(0)}=\mathrm{\Gamma }_{AA^{}}^{(0),a}\frac{2s}{t}\mathrm{\Gamma }_{BB^{}}^{(0),a}.$$
(5)
Inserting this expression into (3), we arrive at (2).

In the next order, $`\alpha _s^3`$, eq.(2) receives several corrections. In order to illustrate the general pattern it is instructive to start from (3). There are three new contributions on the right-hand-side of (3). The first one comes from the same intermediate state $`A^{}`$, $`B^{}`$ as in the leading order case, but there are higher order corrections inside the vertex functions $`\mathrm{\Gamma }_{AA^{}}`$ or $`\mathrm{\Gamma }_{BB^{}}`$ of the scattering amplitude $`A+BA^{}+B^{}`$ (or its complex conjugate). A second correction is due to the reggeization of the exchanged gluon: from (4) it follows that in next-to-leading order reggeization provides an extra term in $`_{AB}`$ proportional to $`\mathrm{\Gamma }_{AA^{}}^{(0),a}\omega (t)\mathrm{ln}(s/t)\mathrm{\Gamma }_{BB^{}}^{(0),a}`$. These two corrections will be referred to as ‘virtual’ corrections to the impact factor. The third correction comes from the production of an extra gluon in the intermediate state (‘real’ corrections): $`A+BA^{}+g+B^{}`$ (Fig. 1). In the high energy limit, we divide this intermediate state into configurations with one or two large rapidity gaps: in the former case, the extra gluon belongs to the fragmentation region of $`A`$ or $`B`$, in the latter case to the central (or: multiperipheral) region.

Correspondingly, we separate the integral over the phase space $`d\varphi `$. Begin with the configuration with one single rapidity gap. If $`N_A^{}`$ ($`N_B^{}`$) particles in the final state belong to the system $`A^{}`$ ($`B^{}`$) and carry momenta $`k_A^{}^i`$ ($`k_B^{}^j`$), we may write the phase space in $`D`$ dimensions explicitly as

$$d\varphi =d\stackrel{~}{\varphi }_A^{}d\stackrel{~}{\varphi }_B^{}(2\pi )^D\delta ^{(D)}(p_A+p_B_ik_A^{}^i_jk_B^{}^j)$$
(6)
with

$$d\stackrel{~}{\varphi }_A^{}=\underset{i=1}{\overset{N_A}{}}\frac{d^Dk_A^{}^i}{(2\pi )^{D1}}\delta ^+(k_A^{}^i{}_{}{}^{2}m_A^i{}_{}{}^{2})$$
(7)
and

$$d\stackrel{~}{\varphi }_B^{}=\underset{j=1}{\overset{N_B}{}}\frac{d^Dk_B^{}^j}{(2\pi )^{D1}}\delta ^+(k_B^{}^j{}_{}{}^{2}m_B^j{}_{}{}^{2}).$$
(8)
Introducing the explicit integration over the reggeon momentum $`r`$ with a Sudakov decomposition $`r=\alpha _rq^{}+\beta _rp+r_{}`$ via one extra delta function and

$$\frac{d^Dr}{(2\pi )^D}=\frac{s}{2}\frac{d\alpha _r}{2\pi }\frac{d\beta _r}{2\pi }\frac{d^{D2}𝒓}{(2\pi )^{D2}}$$
(9)
we may write

$$d\varphi =d\varphi _A^{}\frac{sd\beta _r}{2\pi }d\varphi _B^{}\frac{sd\alpha _r}{2\pi }\frac{1}{2s}\frac{d^{D2}𝒓}{(2\pi )^{D2}},$$
(10)
where $`d\varphi _A^{}=d\stackrel{~}{\varphi }_A^{}(2\pi )^D\delta ^{(D)}(p_A+r_ik_A^{}^i)`$ and $`d\varphi _B^{}=(2\pi )^Dd\stackrel{~}{\varphi }_B^{}\delta ^{(D)}(p_Br_jk_B^{}^j)`$ are the usual phase space measures for the final state particles belonging to the systems $`A^{}`$ and $`B^{}`$ respectively. For an intermediate state with two rapidity gaps (gluon in the central region) we have two $`t`$-channel reggeon momenta, $`r`$ and $`r^{}`$. Denoting the momentum of the produced gluon by $`\mathrm{}rr^{}=\alpha _{\mathrm{}}q^{}+\beta _{\mathrm{}}p+\mathrm{}_{}`$ eq. (10) will be generalized to:

$$d\varphi =(\frac{1}{2s})^2d\varphi _A^{}\frac{sd\beta _r}{2\pi }\frac{sd\alpha _r}{2\pi }d\varphi _g\frac{sd\beta _r^{}}{2\pi }\frac{sd\alpha _r^{}}{2\pi }d\varphi _B^{}\frac{d^{D2}𝒓}{(2\pi )^{D2}}\frac{d^{D2}𝒓^{}}{(2\pi )^{D2}}$$
(11)
where

$$d\varphi _g=\frac{d^D\mathrm{}}{(2\pi )^{D1}}\delta ^+(\mathrm{}^2)(2\pi )^D\delta ^{(D)}(rr^{}+\mathrm{})$$
(12)
is the phase space measure for the gluon in the central region.

Returning to the impact factor representation (2), the NLO corrections due to the extra gluon in the $`s`$-channel can be grouped into two parts. The gluon in the fragmentation region of $`A`$ or $`B`$ leads to corrections of the impact factor $`\mathrm{\Phi }_A`$ or $`\mathrm{\Phi }_B`$, respectively, whereas the gluon in the central region represents the first rung of the BFKL ladder (in the leading-logarithmic approximation) . Together with the other NLO corrections mentioned before, the entire sum of NLO corrections to the total cross section can be cast into the following form:

$`\sigma _{AB}^{(1)}`$ $`={\displaystyle \frac{1}{s}}\mathrm{Im}T_{AB}^{(1)}(s,t=0)`$ (13)
$`={\displaystyle \frac{d^{D2}𝒓}{(2\pi )^{D2}}\mathrm{\Phi }_A^{(1)}\frac{1}{𝒓^4}\mathrm{\Phi }_B^{(0)}}+{\displaystyle \frac{d^{D2}𝒓}{(2\pi )^{D2}}\mathrm{\Phi }_A^{(0)}\frac{1}{𝒓^4}\mathrm{\Phi }_B^{(1)}}`$ (14)
$`+{\displaystyle \frac{d^{D2}𝒓}{(2\pi )^{D2}}\mathrm{\Phi }_A^{(0)}\mathrm{ln}(s/𝒓^2)2\omega ^{(1)}(𝒓^2)\frac{1}{𝒓^4}\mathrm{\Phi }_B^{(0)}}`$ (15)
$`+\mathrm{ln}(s/s_0){\displaystyle \frac{d^{D2}𝒓}{(2\pi )^{D2}}\frac{d^{D2}𝒓^{}}{(2\pi )^{D2}}\mathrm{\Phi }_A^{(0)}\frac{1}{𝒓^4}𝒦_{\mathrm{r}𝒶𝓁}(𝒓,𝒓^{})\frac{1}{𝒓_{}^{}{}_{}{}^{4}}\mathrm{\Phi }_B^{(0)}}.`$ (16)
Here $`𝒦_{\mathrm{r}𝒶𝓁}(𝒓,𝒓^{})`$ denotes the square of the gluon production vertex, the BFKL kernel (sometimes referred to as the ‘real’ part of the BFKL kernel, i.e. it does not contain the gluon reggeization). Later on it will be convenient to rewrite the reggeization of the gluon (third term on the right-hand-side): putting $`\mathrm{ln}(s/𝒓^2)=\mathrm{ln}(s/s_0)+\mathrm{ln}(s_0/𝒓^2)`$ we combine the $`\mathrm{ln}(s/s_0)`$ piece with the real part of the BFKL kernel: in the sum of both terms the infrared singularity in the limit $`𝒓𝒓^{}0`$ drops out. Note the factorization of eq. (16): as indicated after (4), this feature is a general consequence of Regge theory and (with suitable generalizations of the $`t`$-channel gluon propagators) is expected to hold in arbitrary order perturbation theory. For the impact factor we state the general definition:

$$\mathrm{\Phi }_A=\frac{\delta ^{ab}}{\sqrt{N_c^21}}\underset{A^{}}{}\mathrm{\Gamma }_{AA^{}}^a\mathrm{\Gamma }_{AA^{}}^b𝑑\varphi _A^{}\frac{sd\beta _r}{2\pi }$$
(17)
The average symbol $``$ stands for the sum over color and helicity in the intermediate state $`A^{}`$. In the following we will use the short hand notation

$$\mathrm{\Phi }_A=\underset{A^{}}{}|\mathrm{\Gamma }_{AA^{}}|^2𝑑\varphi _A^{}\frac{sd\beta _r}{2\pi },$$
(18)
In particular, in the square of the matrix elements we have included the colour projector $`\delta ^{ab}/\sqrt{N_c^21}`$.

In this paper we are interested in the NLO corrections to the photon impact factor. The strategy follows from the discussion outlined in this section. The NLO corrections can be divided into the $`q\overline{q}`$ and $`q\overline{q}g`$ intermediate states (named virtual and real corrections, respectively). In the former case we need the NLO corrections to the $`\gamma ^{}q\overline{q}`$ vertex, $`\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(1)}`$. They contain infrared divergences (in dimensional regularization: poles in $`ϵ`$) and are listed in . An additional (infrared divergent) contribution comes from the gluon trajectory function in the third term on the right-hand-side of (16). For the $`q\overline{q}g`$ contribution we need the vertex $`\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}`$, in the Born approximation: since, by definition, the gluon lies in the fragmentation region of the virtual photon, we have to divide the phase space of the produced gluon. In section VI we present a detailed discussion of this rather subtle issue: we will show that the separation of the fragmentation region introduces the scale $`s_0`$, in the last term of (16). The fragmentation region contains soft and collinear divergences. When combining them with the infrared divergences of the $`q\overline{q}`$ state (including those from the gluon trajectory function), the infrared divergences cancel . It is the main purpose of this paper to find a separation of the fragmentation from the central region which allows to define infrared finite combinations of virtual and real corrections to the impact factor.

## III Singular Terms of the Virtual Corrections

For future reference, and in order to define our notation we briefly summarize the Born level impact factor and list the divergent terms of the virtual corrections. The amplitude $`T_{\gamma ^{}q}`$ for the scattering $`\gamma ^{}+q^{}q\overline{q}+q^{}`$ (Fig. 2) has been calculated in  to next-to-leading order in $`\alpha _s`$. We write $`T_{\gamma ^{}q}`$ as an expansion in the strong coupling $`g`$

$$T_{\gamma ^{}q}=g^2T_{\gamma ^{}q}^{(0)}+g^4T_{\gamma ^{}q}^{(1)}.$$
(19)
The momenta are labelled as in Fig. 2. The kinematical variables we use are the centre-of-mass energy $`s`$ of the incoming virtual photon and the incoming quark, the virtuality of the photon, $`Q^2`$, the Bjorken variable $`x=Q^2/2pq`$, the momentum transfers $`t=r^2`$, $`t_a=k^2`$, $`t_b=(qkr)^2`$, and the invariant mass $`M^2`$ of the outgoing $`q\overline{q}`$-system. In addition we use a Sudakov decomposition of the momentum $`k`$ with respect to the light cone momenta $`q^{}=qxp`$ and $`p`$ with $`2pq^{}=s`$ , $`k=\alpha q^{}+\beta p+k_{}`$. Two dimensional transverse momenta are denoted as $`𝒌^2=k_{}^2`$.

The Born impact factor is obtained from $`T_{\gamma ^{}q}^{(0)}`$. To leading order we may write

$$T_{\gamma ^{}q}^{(0)}=\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0),a}\frac{2s}{t}\mathrm{\Gamma }_{qq}^{(0),a}$$
(20)
with the Born level vertices

$`\mathrm{\Gamma }_{qq}^{(0),a}`$ $`={\displaystyle \frac{1}{s}}\overline{u}(pr,\lambda _q^{})\overline{)}q^{}\lambda ^au(p,\lambda _q),`$ (21)
$`\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0),a}`$ $`=iee_f\left({\displaystyle \frac{H_T^a}{st_a}}{\displaystyle \frac{\overline{H}_T^a}{st_b}}\right),`$ (22)
where

$`H_T^a`$ $`=`$ $`\overline{u}(k+r,\lambda )\overline{)}p\overline{)}k\overline{)}\epsilon \lambda ^av(qk,\lambda ^{}),`$ (23)
$`\overline{H}_T^a`$ $`=`$ $`\overline{u}(k+r,\lambda )\overline{)}\epsilon (\overline{)}q\overline{)}k\overline{)}r)\overline{)}p\lambda ^av(qk,\lambda ^{})`$ (24)
and $`\lambda ^a`$ are the generators of the colour group $`SU(N_c)`$. Squaring the Born level vertices (21) and averaging (summing) over incoming (outgoing) colours and helicities we obtain for the coupling to the reggeized gluons with colours $`(a,b)`$ , contracted with the colour singlet projector (cf. eqs. (17) and (18)):

$$|\mathrm{\Gamma }_{qq}^{(0)}|^2=\frac{\sqrt{N_c^21}}{2N_c}.$$
(25)
The $`\gamma ^{}q\overline{q}`$-vertex $`|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}|^2=|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}(𝒌,\alpha )|^2`$ for a longitudinally polarized $`\gamma ^{}`$ coupling to a quark flavour $`f`$ reads

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}|_L^2=4\sqrt{N_c^21}e^2e_f^2\alpha ^3(1\alpha )^3Q^2\left(\frac{1}{D(𝒌)}\frac{1}{D(𝒌+𝒓)}\right)^2,$$
(26)
where $`D(𝒌)=𝒌^2+\alpha (1\alpha )Q^2`$. Here the exchanged reggeized gluons must be in the colour singlet. For the sum over the two transverse photon polarizations we have

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}|_T^2=2\sqrt{N_c^21}e^2e_f^2\alpha (1\alpha )\left[\alpha ^2+(1\alpha )^2\right]\left(\frac{𝒌}{D(𝒌)}\frac{𝒌+𝒓}{D(𝒌+𝒓)}\right)^2.$$
(27)
The LO impact factor follows from the definition in eq. (18), from the two-particle phase space measure

$$d\varphi _{q\overline{q}}\frac{sd\beta _r}{2\pi }=\frac{d\alpha }{2\alpha (1\alpha )}\frac{d^{D2}𝒌}{(2\pi )^{D1}}$$
(28)
and from the Born level vertices (26, 27). Expanding $`\mathrm{\Phi }_\gamma ^{}`$ in the strong coupling,

$$\mathrm{\Phi }_\gamma ^{}=g^2\mathrm{\Phi }_\gamma ^{}^{(0)}+g^4\mathrm{\Phi }_\gamma ^{}^{(1)},$$
(29)
we write the Born level impact factor as

$$\mathrm{\Phi }_{\gamma ^{};T,L}^{(0)}=\frac{d^{D2}𝒌}{(2\pi )^{D1}}\frac{d\alpha }{2\alpha (1\alpha )}|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}|_{T,L}^2d^{D2}𝒌𝑑\alpha _{2;T,L}(\alpha ,𝒌;𝒓,Q).$$
(30)
The last equation defines the integrand of the LO impact factor:

$$_{2;L}(\alpha ,𝒌;𝒓,Q)=\frac{2e^2e_f^2\sqrt{N_c^21}}{(2\pi )^{D1}}\alpha ^2(1\alpha )^2Q^2\left(\frac{1}{D(𝒌)}\frac{1}{D(𝒌+𝒓)}\right)^2$$
(31)
$$_{2;T}(\alpha ,𝒌;𝒓,Q)=\frac{e^2e_f^2\sqrt{N_c^21}}{(2\pi )^{D1}}\left[\alpha ^2+(1\alpha )^2\right]\left(\frac{𝒌}{D(𝒌)}\frac{𝒌+𝒓}{D(𝒌+𝒓)}\right)^2$$
(32)
Because of the symmetry of the squared matrix element under the exchange $`q\overline{q}`$, namely $`\alpha 1\alpha `$, $`𝒌𝒓𝒌`$, we have the relation $`_2(\alpha ,𝒌;𝒓,Q)=_2(1\alpha ,𝒓𝒌;𝒓,Q)`$.

Turning to the NLO part of (19), we split $`T_{\gamma ^{}q}`$ into a finite part $`F`$ and a divergent part $`D`$ as

$$T_{\gamma ^{}q}^{(1)}=F_{\gamma ^{}q}+\frac{c_\mathrm{\Gamma }}{(4\pi )^{2ϵ}}T_{\gamma ^{}q}^{(0)}D,c_\mathrm{\Gamma }\frac{\mathrm{\Gamma }(1+ϵ)\mathrm{\Gamma }^2(1ϵ)}{\mathrm{\Gamma }(12ϵ)}.$$
(33)
The divergences are given in dimensional regularization and appear as poles when the space-time dimension $`D=42ϵ`$ approaches $`D=4`$. Here we explicitly include the divergent contribution from self-energy insertions on all four external quark lines

$$D_{\mathrm{SE}}=2C_F\left(\frac{1}{ϵ_{\mathrm{UV}}}\frac{1}{ϵ}\right),$$
(34)
labelling poles from ultraviolet divergences explicitly by $`ϵ_{\mathrm{UV}}`$. As usual, the $`SU(N_c)`$ invariants are denoted by $`C_A=N_c`$ and $`C_F=(N_c^21)/(2N_c)`$.

The divergent part can then be written as the sum of three terms

$$D=D_{\gamma ^{}q\overline{q}}+D_\omega +D_{qq},$$
(35)
following the expansion of the Regge ansatz (4) in powers of $`\alpha _s`$. The three terms denote the divergent parts of the vertex $`\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(1)}`$, of the Regge trajectory of the gluon, $`\omega ^{(1)}(t)`$, which is accompanied by the large leading logarithms in $`s`$, and of the vertex $`\mathrm{\Gamma }_{qq}^{(1)}`$, respectively. They have the form:

$$D_{\gamma ^{}q\overline{q}}=\frac{2C_F}{ϵ^2}\frac{3C_F}{ϵ}+\frac{\beta _0}{ϵ_{\mathrm{UV}}}+\frac{2C_F}{ϵ}\mathrm{ln}\left(M^2\right)+\frac{C_A}{ϵ}\left(\mathrm{ln}\frac{\alpha (1\alpha )t}{M^2}+\mathrm{ln}\frac{s_0}{t}\right),$$
(36)
$$D_\omega =\frac{C_A}{ϵ}\left[\mathrm{ln}\frac{s}{s_0}+\mathrm{ln}\frac{s}{s_0}\right]$$
(37)
and

$$D_{qq}=\frac{2C_F}{ϵ^2}\frac{3C_F}{ϵ}+\frac{\beta _0}{ϵ_{\mathrm{UV}}}+\frac{2C_F}{ϵ}\mathrm{ln}\left(t\right)+\frac{C_A}{ϵ}\mathrm{ln}\frac{s_0}{t}.$$
(38)
Here, we have separated infrared and ultraviolet divergences, making use of the formulae given in . Note that, compared to the Regge ansatz (4), we have changed the energy scale. In (37), we have introduced the energy scale $`s_0`$: instead of $`\mathrm{ln}s/(t)`$ we write $`\mathrm{ln}s/(t)=\mathrm{ln}s/s_0+\mathrm{ln}s_0/(t)`$ and absorb the second term into $`D_{\gamma ^{}q\overline{q}}`$ and $`D_{qq}`$ in (36, 38). Let us first consider the infrared singular terms. For the virtual NLO corrections to the photon impact factor we start from

$$\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{v}irtual)}=\frac{\delta _{ab}}{\sqrt{N_c^21}}\frac{d^{D2}𝒌}{(2\pi )^{D1}}\frac{d\alpha }{2\alpha (1\alpha )}\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(1),a}\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0),b}𝑑\varphi _{q\overline{q}}\frac{sd\beta _r}{2\pi }+\text{c.c.}.$$
(39)
Making use of eqs. (33) and (36) — the latter one without the UV pole which we will discuss later — we can immediately deduce the IR divergent part of the virtual correction to the impact factor:

$`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{v}irtual)}|^{\mathrm{I}Rdivergent}={\displaystyle \frac{c_\mathrm{\Gamma }}{(4\pi )^{2ϵ}}}\mathrm{\Phi }_\gamma ^{}^{(0)}[D_{\gamma ^{}q\overline{q}}^{\mathrm{I}R}+\text{c.c.}]`$ (40)
$`={\displaystyle 𝑑𝒌_0^1𝑑\alpha _2(\alpha ,𝒌;𝒓,Q)\frac{c_\mathrm{\Gamma }}{(4\pi )^{2ϵ}}}`$ (41)
$`\times \left\{{\displaystyle \frac{C_A}{ϵ}}\left[2\mathrm{log}s_0+2\mathrm{log}\alpha (1\alpha )2\mathrm{log}M^2\right]+{\displaystyle \frac{C_F}{ϵ}}\left[{\displaystyle \frac{4}{ϵ}}+4\mathrm{log}M^26\right]\right\}.`$ (42)
Below we will show that these divergent terms cancel against the real corrections.

Turning to the ultraviolet divergent pieces in $`D_{\gamma ^{}q\overline{q}}`$ and $`D_{qq}`$, they are proportional to $`\beta _0=(11N_c2n_f)/6`$, and in  it has been shown that they are removed by the renormalization of the strong coupling. It is instructive to repeat the argument in more detail. The ultraviolet part of (36), including constants, is obtained by adding, in , eqs. (56), (57), (61), (62), (67), plus one half of eq. (60):

$$\frac{g^4c_\mathrm{\Gamma }}{(4\pi )^{2ϵ}}\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0),a}\left[\frac{(𝒓^2)^ϵ}{ϵ_{UV}}\left(\frac{14N_c}{6}\frac{1}{3}n_f\frac{1}{2N_c}\right)\frac{C_F}{ϵ_{UV}}\right]$$
(43)
Except for the last term which is due to the wave function renormalization of the massless outgoing quarks, all ultraviolet poles are proportional to $`(𝒓^2)^ϵ`$. It will be convenient to slightly modify the last term: instead of using, for the renormalization of the outgoing massless quarks, the decomposition (34), we write

$$D_{\mathrm{SE}}^{}=2C_F\left(\frac{(𝒓^2)^ϵ}{ϵ_{\mathrm{UV}}}\frac{(𝒓^2)^ϵ}{ϵ}\right).$$
(44)
This allows to simplify (43):

$$\frac{g^4c_\mathrm{\Gamma }}{(4\pi )^{2ϵ}}\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0),a}\frac{(𝒓^2)^ϵ}{ϵ_{UV}}\beta _0$$
(45)
When addressing, in the next step of our program, the calculation of the finite terms of the virtual corrections, we will have to remember that the change (44) in the wave function renormalization introduces a new finite term

$$\frac{g^4}{(4\pi )^2}\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0),a}C_F\mathrm{ln}(𝒓^2).$$
(46)
In order to remove the UV divergences we have to renormalize, i.e. we add the Born approximation with the replacement:

$$g\frac{Z_1g_\mu \mu ^ϵ}{Z_2\sqrt{Z_3}}=g_\mu \mu ^ϵ\left[1\frac{\alpha _s}{4\pi }\beta _0\left(\frac{1}{ϵ_{UV}}\gamma _E+\mathrm{ln}(4\pi )\right)+\mathrm{}\right]$$
(47)
with the renormalization constants taken in the $`\overline{\mathrm{MS}}`$ scheme. When both contributions are added, the ultraviolet poles cancel, and we are left with

$$g^2\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(1;\mathrm{UV})}=\frac{\alpha _s}{4\pi }\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}\beta _0\mathrm{ln}\frac{𝒓^2}{\mu ^2}.$$
(48)
Note that all dependence upon the renormalization scale $`\mu `$ is contained in the logarithm which introduces the scale dependence of the strong coupling constant, $`\alpha _s(𝒓^2)`$. Correspondingly, in the photon impact factor the $`\mu `$-dependence is contained in the contribution:

$$g^2\mathrm{\Phi }_\gamma ^{}^{(1;\mathrm{UV})}=\frac{\alpha _s}{2\pi }\mathrm{\Phi }_\gamma ^{}^{(0)}\beta _0\mathrm{ln}\frac{𝒓^2}{\mu ^2}$$
(49)
as it is required by the renormalization group equation. Our derivation might suggest that $`𝒓^2`$ is the ‘natural’ momentum scale of $`\alpha _s`$ in the photon impact factor: however, it is clear that we can write

$$g^2\mathrm{\Phi }_\gamma ^{}^{(1;\mathrm{UV})}=\frac{\alpha _s}{2\pi }\mathrm{\Phi }_\gamma ^{}^{(0)}\left(\beta _0\mathrm{ln}\frac{Q^2}{\mu ^2}+\beta _0\mathrm{ln}\frac{𝒓^2}{Q^2}\right)$$
(50)
i.e. by simply redefining the $`\mu `$-independent part of the impact factor we can switch from one momentum scale of $`\alpha _s`$, $`\alpha _s(𝒓^2)`$, to another scale, $`\alpha _s(Q^2)`$.

Stating that (49) is the only $`\mu `$-dependent term in the impact factor implies that the $`\mu `$-dependence drops out of all the remaining contributions to the impact factor. In fact, we silently dropped a $`\mu `$ already in eq. (46), anticipating the scale invariance of the finite terms which are left after the cancellation of the infrared divergencies of real and virtual corrections. This can be understood as follows. Without having it written explicitly in , the infrared divergences in dimensional regularization appear as poles in $`ϵ`$, accompanied by the power of some dimensionful scale, $`s_i`$. In a one-loop calculation the sum of infrared divergent terms from real (R) and virtual (V) corrections generally reads

$$T_{\mathrm{IR}}=\underset{i}{}A_i\frac{(s_i)^ϵ}{ϵ^2}+B_i\frac{(s_i)^ϵ}{ϵ}.$$
(51)
Adding real and virtual corrections, the coefficients $`A_i`$ and $`B_i`$ are sums as well,

$$A_i=A_i^V+A_i^R,B_i=B_i^V+B_i^R.$$
(52)
Expanding (51) in $`ϵ`$ we obtain

$$T_{\mathrm{IR}}=\underset{i}{}\frac{A_i}{ϵ^2}\frac{A_i}{ϵ}\mathrm{ln}s_i+\frac{A_i}{2}\mathrm{ln}^2s_i+\frac{B_i}{ϵ}B_i\mathrm{ln}s_i.$$
(53)
Introducing the scale $`\mu `$ with the renormalized coupling as a factor $`\mu ^{2ϵ}`$ and reexpanding results in

$`T_{\mathrm{IR}}=`$ $`{\displaystyle \underset{i}{}}\left({\displaystyle \frac{1}{ϵ^2}}+{\displaystyle \frac{\mathrm{ln}\mu ^2}{ϵ}}+{\displaystyle \frac{\mathrm{ln}^2\mu ^2}{2}}\right)A_i\left({\displaystyle \frac{1}{ϵ}}+\mathrm{ln}\mu ^2\right)A_i\mathrm{ln}s_i+{\displaystyle \frac{A_i}{2}}\mathrm{ln}^2s_i`$ (54)
$`+\left({\displaystyle \frac{1}{ϵ}}+\mathrm{ln}\mu ^2\right)B_iB_i\mathrm{ln}s_i.`$ (55)
Further below (section VII) we will show that in the sum of virtual and real corrections all $`ϵ`$-poles cancel: this means that $`A_i=0`$ and $`_iB_i=0`$. Hence, we see that the cancellation of infrared divergences implies that also the $`\mathrm{ln}\mu ^2`$ terms cancel. This allows us to write e.g. the logarithm in (46) with a dimensionful argument. We did the same with the results of diagram 14 in  after expanding divergent terms of the form $`s_i^ϵ`$.

## IV Real Corrections

We now turn to the real corrections. We study the process $`\gamma ^{}+q(q\overline{q}g)+q`$ with transverse photon polarization (Fig.3); the case of the longitudinal photon has been studied in , and we follow the procedure outlined in this earlier paper. The diagrams are shown in Fig. 4. We use the Sudakov decomposition of the gluon momentum $`\mathrm{}=\alpha _{\mathrm{}}q^{}+\beta _{\mathrm{}}p+\mathrm{}_{}`$ and introduce the following abbreviations :

$`\alpha _1`$ $`\alpha ,`$ (56)
$`\overline{\alpha }_1`$ $`(1\alpha \alpha _{\mathrm{}}),`$ (57)
$`\alpha _2`$ $`(1\alpha ),`$ (58)
$`\overline{\alpha }_2`$ $`(\alpha +\alpha _{\mathrm{}}).`$ (59)
The propagator denominators, that occur in the calculation, have the following form:

$`D_1`$ $`=(k+\mathrm{}+rq)^2`$ (60)
$`=\overline{\alpha }_1Q^2(𝒌+\mathbf{}+𝒓)^2{\displaystyle \frac{\overline{\alpha }_1}{\alpha _1}}(𝒌+𝒓)^2{\displaystyle \frac{\overline{\alpha }_1}{\alpha _{\mathrm{}}}}\mathbf{}^2`$ (61)
$`D_2`$ $`=(k+rq)^2`$ (62)
$`=\alpha _2Q^2{\displaystyle \frac{1}{\alpha _1}}(𝒌+𝒓)^2`$ (63)
$`D_3`$ $`=(k+\mathrm{}+r)^2`$ (64)
$`=(𝒌+\mathbf{}+𝒓)^2+{\displaystyle \frac{\overline{\alpha }_2}{\alpha _1}}(𝒌+𝒓)^2+{\displaystyle \frac{\overline{\alpha }_2}{\alpha _{\mathrm{}}}}\mathbf{}^2`$ (65)
$`={\displaystyle \frac{1}{\alpha _{\mathrm{}}\alpha _1}}(\alpha _{\mathrm{}}(𝒌+𝒓)+\alpha _1\mathbf{})^2`$ (66)
$`D_4`$ $`=(kq)^2`$ (67)
$`=𝒌^2+{\displaystyle \frac{\alpha _2}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2+{\displaystyle \frac{\alpha _2}{\alpha _{\mathrm{}}}}\mathbf{}^2`$ (68)
$`={\displaystyle \frac{1}{\alpha _{\mathrm{}}\overline{\alpha }_1}}(\alpha _{\mathrm{}}(𝒌+\mathbf{})+\overline{\alpha }_1\mathbf{})^2`$ (69)
$`D_5`$ $`=(\mathrm{}r)^2`$ (70)
$`=\alpha _{\mathrm{}}Q^2{\displaystyle \frac{\alpha _{\mathrm{}}}{\alpha _1}}(𝒌+𝒓)^2{\displaystyle \frac{\alpha _{\mathrm{}}}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2(\mathbf{}𝒓)^2`$ (71)
$`D_6`$ $`=k^2`$ (72)
$`=\alpha _1Q^2𝒌^2{\displaystyle \frac{\alpha _1}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2{\displaystyle \frac{\alpha _1}{\alpha _{\mathrm{}}}}\mathbf{}^2`$ (73)
$`D_7`$ $`=(k+\mathrm{})^2`$ (74)
$`=\overline{\alpha }_2Q^2{\displaystyle \frac{1}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2.`$ (75)
We present the results for $`|_{q\overline{q}g}|^2`$, the squared matrix element, averaged (summed) over the helicities and colors of the incoming (outgoing) quarks and summed over equal photon polarizations; thereby we made use of $`_{i=j=1,2}(ϵ_T^ik_{})(ϵ_T^j\mathrm{}_{})=𝒌\mathbf{}`$. Our results can be simplified by expressing them in terms of the matrix elements for longitudinal photons; for this reason we also recall the results for longitudinal photon polarization . The matrix elements for the process $`\gamma ^{}q(q\overline{q}g)q`$ can be written in the following factorized form:

$$|_{q\overline{q}g}|^2=\frac{1}{N_c^21}g^2\delta ^{ab}|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2\left(\frac{2s}{t}\right)^2\delta ^{ab}|\mathrm{\Gamma }_{qq}^{(0)}|^2,$$
(76)
and by using eq.(25) we extract the squared vertex $`|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2`$ (summed over helicities and color in the $`q\overline{q}g`$ state). Following the notation used in  we list our results in correspondence with the Feynman diagrams:

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|_L^2=4e^2e_f^2g^4\sqrt{N_c^21}Q^2\underset{i,j=1}{\overset{5}{}}\left(𝒜𝒜_{ij}^L+𝒜_{ij}^L+𝒜_{ij}^L+_{ij}^L\right),$$
(77)
$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|_T^2=e^2e_f^2g^4\sqrt{N_c^21}\underset{i,j=1}{\overset{5}{}}\left(𝒜𝒜_{ij}^T+𝒜_{ij}^T+𝒜_{ij}^T+_{ij}^T\right);$$
(78)
$`𝒜𝒜_{11}^T=`$ $`C_F{\displaystyle \frac{4}{D_{1}^{}{}_{}{}^{2}D_{2}^{}{}_{}{}^{2}}}\alpha _1\overline{\alpha }_1(1ϵ)((𝒌+\mathbf{}+𝒓)^2(𝒌+𝒓)^2`$ (79)
$`+{\displaystyle \frac{1}{\alpha _1^2}}\overline{\alpha }_1((𝒌+𝒓)^2+\alpha _1Q^2)(2\alpha _1(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+\overline{\alpha }_1((𝒌+𝒓)^2+\alpha _1Q^2)))`$ (80)
$`+{\displaystyle \frac{2(𝒌+𝒓)^2}{\alpha _{1}^{}{}_{}{}^{2}}}𝒜𝒜_{11}^L`$ (81)
$`𝒜𝒜_{12}^T=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{2}{D_1D_{2}^{}{}_{}{}^{2}D_4}}[\overline{\alpha }_1(𝒌+𝒓)^2(𝒌𝒓+\alpha _2𝒓^2)+\alpha _1\alpha _2\overline{\alpha }_1𝒓(𝒌+𝒓)Q^2`$ (82)
$`+(1ϵ)((\alpha _1\overline{\alpha }_1𝒌(𝒌+\mathbf{}+𝒓)\overline{\alpha }_{1}^{}{}_{}{}^{2}𝒌(𝒌+𝒓)+\alpha _1\alpha _2(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)`$ (83)
$`+\alpha _2\overline{\alpha }_1(𝒌+\mathbf{})(𝒌+𝒓))(𝒌+𝒓)^2`$ (84)
$`\alpha _1\overline{\alpha }_1(\overline{\alpha }_1𝒌(𝒌+𝒓)\alpha _2(𝒌+\mathbf{})(𝒌+𝒓))Q^2)]+{\displaystyle \frac{2(𝒌+𝒓)^2}{\alpha _{1}^{}{}_{}{}^{2}}}𝒜𝒜_{12}^L`$ (85)
$`𝒜𝒜_{13}^T=`$ $`C_F{\displaystyle \frac{4}{\alpha _{1}^{}{}_{}{}^{2}\alpha _{\mathrm{}}D_{1}^{}{}_{}{}^{2}D_2D_3}}\overline{\alpha }_1[\alpha _{1}^{}{}_{}{}^{3}(\overline{\alpha }_2\alpha _{\mathrm{}}(𝒌+\mathbf{}+𝒓)^2\overline{\alpha }_{1}^{}{}_{}{}^{2}\mathbf{}^2)Q^2`$ (86)
$`+\alpha _{1}^{}{}_{}{}^{2}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\left(\alpha _{\mathrm{}}(𝒌+\mathbf{}+𝒓)^2+\overline{\alpha }_1\mathbf{}^2\overline{\alpha }_1\alpha _{\mathrm{}}Q^2\right)+\alpha _{\mathrm{}}(𝒌+𝒓)^2`$ (87)
$`(\alpha _1\overline{\alpha }_2(𝒌+\mathbf{}+𝒓)^2+\overline{\alpha }_{1}^{}{}_{}{}^{2}(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_1(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)`$ (88)
$`\alpha _1\overline{\alpha }_{1}^{}{}_{}{}^{2}({\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}}\alpha _1Q^2))`$ (89)
$`+\alpha _1(1ϵ)(\alpha _{1}^{}{}_{}{}^{2}\overline{\alpha }_{1}^{}{}_{}{}^{2}\mathbf{}^2Q^2\alpha _1\alpha _{\mathrm{}}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)(\alpha _2(𝒌+\mathbf{}+𝒓)^2`$ (90)
$`4\overline{\alpha }_1(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+{\displaystyle \frac{\alpha _1\overline{\alpha }_1\mathbf{}^2}{\alpha _{\mathrm{}}}}\overline{\alpha }_{1}^{}{}_{}{}^{2}Q^2)`$ (91)
$`+\alpha _{\mathrm{}}(𝒌+𝒓)^2((2\alpha _1\overline{\alpha }_1+\alpha _2\overline{\alpha }_2)(𝒌+\mathbf{}+𝒓)^2\overline{\alpha }_{1}^{}{}_{}{}^{2}(𝒌+𝒓)^2`$ (92)
$`+\overline{\alpha }_1(\alpha _1+\overline{\alpha }_1\overline{\alpha }_2)(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+\alpha _1\overline{\alpha }_{1}^{}{}_{}{}^{2}({\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}}Q^2)))]`$ (93)
$`𝒜𝒜_{14}^T=`$ $`N_c{\displaystyle \frac{1}{D_1D_{2}^{}{}_{}{}^{2}D_5}}[\overline{\alpha }_{1}^{}{}_{}{}^{2}(𝒌+𝒓)^2((𝒌+𝒓)^2+{\displaystyle \frac{\alpha _1\mathbf{}^2}{\alpha _{\mathrm{}}}}\alpha _1Q^2)`$ (94)
$`+{\displaystyle \frac{2\alpha _{\mathrm{}}}{\alpha _1}}(1ϵ)(\alpha _1\overline{\alpha }_1((𝒌+𝒓)^2+\alpha _1Q^2)(\mathrm{\hspace{0.17em}2}(𝒌+\mathbf{})(𝒌+𝒓)+𝒓(𝒌+𝒓)+\overline{\alpha }_1Q^2)`$ (95)
$`+(𝒌+𝒓)^2(\alpha _{1}^{}{}_{}{}^{2}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)+\overline{\alpha }_{1}^{}{}_{}{}^{2}(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_{1}^{}{}_{}{}^{2}Q^2))`$ (96)
$`+\alpha _1\overline{\alpha }_1(𝒌+𝒓)^2\left(𝒌(𝒌+\mathbf{})+(\mathbf{}𝒓)(𝒌+\mathbf{}+𝒓)𝒓^2\right)`$ (97)
$`+\overline{\alpha }_1((𝒌+𝒓)^2\alpha _1Q^2)(\alpha _{\mathrm{}}(\mathrm{\hspace{0.17em}2}(𝒌+\mathbf{})(𝒌+𝒓)+𝒓(𝒌+𝒓))`$ (98)
$`+\overline{\alpha }_1(2\mathbf{}𝒓)(𝒌+𝒓))]`$ (99)
$`+{\displaystyle \frac{2(𝒌+𝒓)^2}{\alpha _{1}^{}{}_{}{}^{2}}}𝒜𝒜_{14}^L`$ (100)
$`𝒜𝒜_{15}^T=`$ $`N_c{\displaystyle \frac{2}{\alpha _1\alpha _{\mathrm{}}D_1D_{2}^{}{}_{}{}^{2}D_5}}(\alpha _{1}^{}{}_{}{}^{2}+\alpha _{2}^{}{}_{}{}^{2})\overline{\alpha }_{1}^{}{}_{}{}^{2}(𝒌+𝒓)^2𝒓^2`$ (101)
$`𝒜𝒜_{22}^T=`$ $`{\displaystyle \frac{2(\alpha _{1}^{}{}_{}{}^{2}+\alpha _{2}^{}{}_{}{}^{2})(𝒌+𝒓)^2}{\alpha _{1}^{}{}_{}{}^{2}\alpha _{2}^{}{}_{}{}^{2}}}𝒜𝒜_{22}^L`$ (102)
$`𝒜𝒜_{23}^T=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{2}{D_1D_2D_3D_4}}[(\alpha _1+\overline{\alpha }_1)(𝒌(𝒌+𝒓)(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)`$ (103)
$`+𝒌(𝒌+\mathbf{})(𝒌+𝒓)(𝒌+\mathbf{}+𝒓))`$ (104)
$`+(\alpha _1(𝒌+\mathbf{})(𝒌+𝒓)+\overline{\alpha }_1(𝒌+𝒓)^2)({\displaystyle \frac{(\alpha _1+\overline{\alpha }_1)}{\alpha _1}}𝒌(𝒌+\mathbf{}+𝒓)`$ (105)
$`\alpha _2(𝒌+\mathbf{}+𝒓)^2+{\displaystyle \frac{\alpha _{2}^{}{}_{}{}^{2}\overline{\alpha }_1}{\alpha _1}}({\displaystyle \frac{(𝒌+𝒓)^2}{\alpha _1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}}))`$ (106)
$`+ϵ((\overline{\alpha }_1𝒌(𝒌+\mathbf{}+𝒓)+\alpha _2(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓))({\displaystyle \frac{(\alpha _1\overline{\alpha }_2)}{\alpha _1}}(𝒌+𝒓)^2`$ (107)
$`+2(𝒌+𝒓)(𝒌+\mathbf{}+𝒓))`$ (108)
$`+(\overline{\alpha }_1𝒌(𝒌+𝒓)+\alpha _2(𝒌+\mathbf{})(𝒌+𝒓))(\alpha _2(𝒌+\mathbf{}+𝒓)^2`$ (109)
$`2(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+\alpha _1\overline{\alpha }_1({\displaystyle \frac{(𝒌+𝒓)^2}{\alpha _1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}})))]`$ (110)
$`+{\displaystyle \frac{4(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)}{\alpha _2\overline{\alpha }_1}}𝒜𝒜_{23}^L`$ (111)
$`𝒜𝒜_{24}^T=`$ $`N_c{\displaystyle \frac{1}{D_{2}^{}{}_{}{}^{2}D_4D_5}}[{\displaystyle \frac{2\alpha _{\mathrm{}}(1ϵ)}{\overline{\alpha }_1}}(\alpha _1(\overline{\alpha }_1𝒌(𝒌+\mathbf{})+\alpha _2(𝒌+\mathbf{})^2)(𝒌+𝒓)^2`$ (112)
$`+\overline{\alpha }_1(\overline{\alpha }_1𝒌(𝒌+𝒓)+\alpha _2(𝒌+\mathbf{})(𝒌+𝒓))((𝒌+𝒓)^2+\alpha _1Q^2))`$ (113)
$`+\alpha _1(𝒌+𝒓)^2(\alpha _{\mathrm{}}𝒌(𝒌+\mathbf{})\alpha _2(𝒌+\mathbf{})^2+\overline{\alpha }_1𝒌(\mathbf{}2𝒓)+{\displaystyle \frac{\alpha _2\overline{\alpha }_1\mathbf{}^2}{\alpha _{\mathrm{}}}}`$ (114)
$`+\alpha _2(2\mathbf{}𝒓)(𝒌+\mathbf{}))`$ (115)
$`\alpha _2((𝒌+𝒓)^2+\alpha _1Q^2)(\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+𝒓)+\overline{\alpha }_1(𝒌+𝒓)^2+\overline{\alpha }_1(\mathbf{}+𝒓)(𝒌+𝒓))]`$ (116)
$`+{\displaystyle \frac{2(𝒌+𝒓)^2}{\alpha _{1}^{}{}_{}{}^{2}}}𝒜𝒜_{24}^L`$ (117)
$`𝒜𝒜_{25}^T=`$ $`N_c{\displaystyle \frac{2}{\alpha _1\alpha _{\mathrm{}}D_{2}^{}{}_{}{}^{2}D_4D_5}}\alpha _2(\alpha _{1}^{}{}_{}{}^{2}+\alpha _{2}^{}{}_{}{}^{2})\overline{\alpha }_1(𝒌+𝒓)^2𝒓^2`$ (118)
$`𝒜𝒜_{33}^T=`$ $`\mathrm{\hspace{0.33em}2}\left[{\displaystyle \frac{(𝒌+\mathbf{}+𝒓)^2}{\overline{\alpha }_{1}^{}{}_{}{}^{2}}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}^{}{}_{}{}^{2}}}\right]𝒜𝒜_{33}^L`$ (119)
$`𝒜𝒜_{34}^T=`$ $`N_c{\displaystyle \frac{1}{D_1D_2D_3D_5}}[2ϵ({\displaystyle \frac{2}{\alpha _2}}\overline{\alpha }_1\alpha _{\mathrm{}}((𝒌+𝒓)(𝒌+\mathbf{}+𝒓))^2`$ (120)
$`+(𝒌+\mathbf{})(𝒌+𝒓)\left(\alpha _2\alpha _{\mathrm{}}(𝒌+\mathbf{}+𝒓)^2+\overline{\alpha }_1\alpha _{\mathrm{}}(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_1\mathbf{}^2\right)`$ (121)
$`+\alpha _1\overline{\alpha }_{1}^{}{}_{}{}^{2}\mathbf{}^2Q^2+{\displaystyle \frac{\alpha _{\mathrm{}}}{\alpha _1\alpha _2}}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)(\overline{\alpha }_1(2\alpha _1\overline{\alpha }_2+\alpha _2(\overline{\alpha }_1+\overline{\alpha }_2))(𝒌+𝒓)^2`$ (122)
$`+2\alpha _1\alpha _2\mathbf{}(𝒌+\mathbf{})\alpha _1\alpha _2\overline{\alpha }_{1}^{}{}_{}{}^{2}Q^2)`$ (123)
$`+\overline{\alpha }_1(𝒌+𝒓)^2({\displaystyle \frac{(\alpha _1\overline{\alpha }_2)\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)}{\alpha _1\overline{\alpha }_1}}+\alpha _{\mathrm{}}(𝒌+\mathbf{}+𝒓)^2`$ (124)
$`+{\displaystyle \frac{\overline{\alpha }_1\alpha _{\mathrm{}}(𝒌+𝒓)^2}{\alpha _1}}+\overline{\alpha }_1\mathbf{}^2+\overline{\alpha }_1\alpha _{\mathrm{}}Q^2))`$ (125)
$`+{\displaystyle \frac{3(\alpha _1\overline{\alpha }_1)\overline{\alpha }_1(𝒌+𝒓)^2𝒓(𝒌+\mathbf{}+𝒓)}{\alpha _1}}`$ (126)
$`+((\alpha _1\overline{\alpha }_1)(𝒌+\mathbf{})(𝒌+𝒓)+{\displaystyle \frac{\overline{\alpha }_{1}^{}{}_{}{}^{2}(𝒌+𝒓)^2}{\alpha _1}}+\alpha _1\overline{\alpha }_1\overline{\alpha }_2Q^2)(\mathbf{}+𝒓)(𝒌+\mathbf{}+𝒓)`$ (127)
$`{\displaystyle \frac{1}{\alpha _{1}^{}{}_{}{}^{2}\alpha _{\mathrm{}}}}\left(\alpha _{1}^{}{}_{}{}^{2}\alpha _{\mathrm{}}(𝒌+\mathbf{}+𝒓)^2\alpha _2\overline{\alpha }_1\alpha _{\mathrm{}}(𝒌+𝒓)^2\alpha _1\alpha _2\overline{\alpha }_1\mathbf{}^2\right)`$ (128)
$`\left(2\alpha _1\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+𝒓)+\overline{\alpha }_1\alpha _{\mathrm{}}(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_1(\mathbf{}2𝒓)(𝒌+𝒓)\right)`$ (129)
$`(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)\left(\alpha _{\mathrm{}}(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_2\alpha _{\mathrm{}}Q^2+(\alpha _1\overline{\alpha }_1)(\mathbf{}+𝒓)(𝒌+𝒓)\right)`$ (130)
$`+(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)({\displaystyle \frac{\overline{\alpha }_1}{\alpha _1\alpha _2}}(\alpha _{2}^{}{}_{}{}^{2}\overline{\alpha }_1\alpha _1\alpha _2\overline{\alpha }_2+2\alpha _2(\overline{\alpha }_1\overline{\alpha }_2)\alpha _{\mathrm{}}`$ (131)
$`4\alpha _1\overline{\alpha }_2\alpha _{\mathrm{}})(𝒌+𝒓)^22\overline{\alpha }_1\overline{\alpha }_2\alpha _{\mathrm{}}Q^2`$ (132)
$`{\displaystyle \frac{(\alpha _1+\overline{\alpha }_1)}{\alpha _{\mathrm{}}}}(\alpha _{\mathrm{}}(𝒌+\mathbf{})^2+\overline{\alpha }_1\mathbf{}^2\overline{\alpha }_1\alpha _{\mathrm{}}Q^2)`$ (133)
$`+{\displaystyle \frac{(\alpha _1+\overline{\alpha }_1)}{\alpha _2}}(2\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+𝒓)\alpha _2(\mathbf{}+𝒓)(𝒌+\mathbf{})`$ (134)
$`2\overline{\alpha }_1(\mathbf{}+𝒓)(𝒌+𝒓)))]+{\displaystyle \frac{4(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)}{\alpha _2\overline{\alpha }_1}}𝒜𝒜_{34}^L`$ (135)
$`𝒜𝒜_{35}^T=`$ $`N_c{\displaystyle \frac{2}{\alpha _{\mathrm{}}D_1D_2D_3D_5}}\overline{\alpha }_1(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2)(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)𝒓^2`$ (136)
$`𝒜𝒜_{44}^T=`$ $`N_c{\displaystyle \frac{2}{D_{2}^{}{}_{}{}^{2}D_{5}^{}{}_{}{}^{2}}}[{\displaystyle \frac{2\alpha _{\mathrm{}}^{}{}_{}{}^{2}}{\alpha _1\overline{\alpha }_1}}(1ϵ)(\alpha _{1}^{}{}_{}{}^{2}(𝒌+\mathbf{})^2(𝒌+𝒓)^2+\overline{\alpha }_{1}^{}{}_{}{}^{2}(𝒌+𝒓)^4`$ (137)
$`+\alpha _{1}^{}{}_{}{}^{2}\overline{\alpha }_1(𝒌+\mathbf{})(𝒌+𝒓)Q^2`$ (138)
$`+\alpha _1\overline{\alpha }_1(2(𝒌+𝒓)^2+\alpha _1Q^2)((𝒌+\mathbf{})(𝒌+𝒓)+\overline{\alpha }_1Q^2))`$ (139)
$`+\alpha _1(\alpha _{\mathrm{}}(𝒌+\mathbf{})^2(𝒌+𝒓)^2+2\overline{\alpha }_1(𝒌+𝒓)^2((\mathbf{}𝒓)^2+𝒓^2)`$ (140)
$`+\alpha _{\mathrm{}}(𝒌+𝒓)^2({\displaystyle \frac{\overline{\alpha }_1(𝒌+𝒓)^2}{\alpha _1}}\overline{\alpha }_1Q^2+3(\mathbf{}𝒓)(𝒌+\mathbf{})))`$ (141)
$`3\alpha _{\mathrm{}}((𝒌+𝒓)^2+\alpha _1Q^2)(\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+𝒓)+\overline{\alpha }_1(\mathbf{}𝒓)(𝒌+𝒓))]`$ (142)
$`+{\displaystyle \frac{2(𝒌+𝒓)^2}{\alpha _{1}^{}{}_{}{}^{2}}}𝒜𝒜_{44}^L`$ (143)
$`𝒜𝒜_{45}^T=`$ $`N_c{\displaystyle \frac{2}{\alpha _1D_{2}^{}{}_{}{}^{2}D_{5}^{}{}_{}{}^{2}}}(\alpha _{1}^{}{}_{}{}^{2}+\alpha _{2}^{}{}_{}{}^{2})\overline{\alpha }_1(𝒌+𝒓)^2𝒓^2`$ (144)
$`𝒜𝒜_{55}^T=`$ $`\mathrm{\hspace{0.33em}0}`$ (145)
$`𝒜_{11}^T=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{2}{D_1D_2D_6D_7}}[(\overline{\alpha }_1𝒌(𝒌+𝒓)+\alpha _1(𝒌+𝒓)(𝒌+\mathbf{}+𝒓))`$ (146)
$`\left((𝒌+\mathbf{})^22\overline{\alpha }_1(𝒌+\mathbf{})(𝒌+𝒓)+\alpha _2\overline{\alpha }_1Q^2\right)`$ (147)
$`+(\overline{\alpha }_1𝒌(𝒌+\mathbf{})+\alpha _1(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓))(2\alpha _1(𝒌+\mathbf{})(𝒌+𝒓)`$ (148)
$`+(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_2Q^2)`$ (149)
$`+ϵ((𝒌+\mathbf{})(𝒌+𝒓)((\alpha _1\overline{\alpha }_1+\alpha _2\overline{\alpha }_2)𝒌(𝒌+\mathbf{}+𝒓)2\overline{\alpha }_1\overline{\alpha }_2𝒌(𝒌+𝒓)`$ (150)
$`2\alpha _1\alpha _2(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)+2\alpha _1\overline{\alpha }_1(𝒌+\mathbf{})(𝒌+𝒓))`$ (151)
$`+(\alpha _1\overline{\alpha }_1+\alpha _2\overline{\alpha }_2)(𝒌(𝒌+𝒓)(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)`$ (152)
$`𝒌(𝒌+\mathbf{})(𝒌+𝒓)(𝒌+\mathbf{}+𝒓))`$ (153)
$`+\overline{\alpha }_1(𝒌+𝒓)^2\left(𝒌(𝒌+\mathbf{})+\alpha _1\overline{\alpha }_1Q^2\right)`$ (154)
$`+\alpha _1(𝒌+\mathbf{})^2\left((𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+\alpha _1\overline{\alpha }_1Q^2\right)`$ (155)
$`+\alpha _1\overline{\alpha }_1Q^2(\overline{\alpha }_1𝒌(𝒌+\mathbf{})+\alpha _1(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+\alpha _1\overline{\alpha }_1Q^2))]`$ (156)
$`𝒜_{12}^T=`$ $`C_F{\displaystyle \frac{4}{D_1D_2D_3D_7}}[{\displaystyle \frac{1}{\alpha _1}}(\alpha _2\overline{\alpha }_{1}^{}{}_{}{}^{2}2\alpha _1\overline{\alpha }_1\overline{\alpha }_2\alpha _1\overline{\alpha }_{2}^{}{}_{}{}^{2})(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)(𝒌+𝒓)^2`$ (157)
$`\alpha _1\overline{\alpha }_{2}^{}{}_{}{}^{2}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)Q^2`$ (158)
$`+{\displaystyle \frac{1}{\alpha _1}}(𝒌+\mathbf{})(𝒌+𝒓)(\alpha _1(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2)(𝒌+\mathbf{}+𝒓)^2`$ (159)
$`+\overline{\alpha }_1(\overline{\alpha }_1\overline{\alpha }_2)\overline{\alpha }_2(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_1\overline{\alpha }_2(4(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\overline{\alpha }_2Q^2))`$ (160)
$`+ϵ({\displaystyle \frac{(\alpha _1\overline{\alpha }_{1}^{}{}_{}{}^{2}+2\alpha _1\overline{\alpha }_1\overline{\alpha }_2\alpha _2\overline{\alpha }_{2}^{}{}_{}{}^{2})}{\alpha _1}}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)(𝒌+𝒓)^2`$ (161)
$`+\alpha _1\overline{\alpha }_{1}^{}{}_{}{}^{2}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)Q^2`$ (162)
$`+{\displaystyle \frac{1}{\alpha _1}}(𝒌+\mathbf{})(𝒌+𝒓)(\alpha _1(\alpha _1\overline{\alpha }_1+\alpha _2\overline{\alpha }_2)(𝒌+\mathbf{}+𝒓)^2`$ (163)
$`\overline{\alpha }_1(\overline{\alpha }_1\overline{\alpha }_2)\overline{\alpha }_2(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_1\overline{\alpha }_2(4(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\overline{\alpha }_1Q^2)))]`$ (164)
$`𝒜_{13}^T=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{2}{D_1D_2D_4D_6}}[(\overline{\alpha }_1𝒌^2\alpha _2𝒌(𝒌+\mathbf{}+𝒓))(\alpha _2(𝒌+\mathbf{})(𝒌+𝒓)\overline{\alpha }_1(𝒌+𝒓)^2)`$ (165)
$`+(\alpha _{1}^{}{}_{}{}^{2}\alpha _{2}^{}{}_{}{}^{2})\left(𝒌(𝒌+𝒓)(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)𝒌(𝒌+\mathbf{})(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\right)`$ (166)
$`\alpha _1\left(𝒌(𝒌+\mathbf{}+𝒓)+(𝒌+\mathbf{})^2+{\displaystyle \frac{\overline{\alpha }_1}{\alpha _{\mathrm{}}}}\mathbf{}^2\right)`$ (167)
$`\left(\alpha _1(𝒌+\mathbf{})(𝒌+𝒓)+\overline{\alpha }_1(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_1Q^2\right)`$ (168)
$`+ϵ(\overline{\alpha }_1(\overline{\alpha }_1𝒌^2+\alpha _2𝒌(𝒌+\mathbf{}))(𝒌+𝒓)^2`$ (169)
$`+(\alpha _{1}^{}{}_{}{}^{2}\alpha _{2}^{}{}_{}{}^{2})(𝒌(𝒌+𝒓)(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)`$ (170)
$`+𝒌(𝒌+\mathbf{}+𝒓)(𝒌+\mathbf{})(𝒌+𝒓)𝒌(𝒌+\mathbf{})(𝒌+𝒓)(𝒌+\mathbf{}+𝒓))`$ (171)
$`+(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\left(\alpha _2\overline{\alpha }_1𝒌^2+\alpha _{1}^{}{}_{}{}^{2}(𝒌+\mathbf{})^2+{\displaystyle \frac{\alpha _{1}^{}{}_{}{}^{2}\overline{\alpha }_1}{\alpha _{\mathrm{}}}}\mathbf{}^2\right)`$ (172)
$`+{\displaystyle \frac{\alpha _1\overline{\alpha }_1}{\alpha _{\mathrm{}}}}(\overline{\alpha }_1\mathbf{}^2+\alpha _{\mathrm{}}\mathbf{}(𝒌+\mathbf{}))((𝒌+𝒓)^2+\alpha _1Q^2))]`$ (173)
$`+{\displaystyle \frac{4𝒌(𝒌+𝒓)}{\alpha _{1}^{}{}_{}{}^{2}}}𝒜_{13}^L`$ (174)
$`𝒜_{14}^T=`$ $`N_c{\displaystyle \frac{1}{D_1D_2D_5D_7}}[2ϵ({\displaystyle \frac{2}{\alpha _2}}\overline{\alpha }_1\alpha _{\mathrm{}}[(𝒌+\mathbf{})(𝒌+𝒓)]^2`$ (175)
$`+{\displaystyle \frac{(\alpha _1\overline{\alpha }_1\alpha _2\overline{\alpha }_2)\alpha _{\mathrm{}}}{\alpha _1}}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)(𝒌+𝒓)^2`$ (176)
$`\alpha _{\mathrm{}}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\left((𝒌+\mathbf{})^2+\alpha _1\overline{\alpha }_1Q^2\right)`$ (177)
$`\overline{\alpha }_1\alpha _{\mathrm{}}Q^2\left(\alpha _1(𝒌+\mathbf{})^2+\overline{\alpha }_1(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_1Q^2\right)`$ (178)
$`+\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+𝒓)(2(\mathbf{}𝒓)(𝒌+\mathbf{}+𝒓)`$ (179)
$`+({\displaystyle \frac{2\overline{\alpha }_1\overline{\alpha }_2}{\alpha _2}}+{\displaystyle \frac{\overline{\alpha }_1(\overline{\alpha }_1+\overline{\alpha }_2)}{\alpha _1}})(𝒌+𝒓)^2\overline{\alpha }_{1}^{}{}_{}{}^{2}Q^2))`$ (180)
$`((𝒌+𝒓)^2+\alpha _1\overline{\alpha }_2Q^2)\left(\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)+\overline{\alpha }_1(\mathbf{}2𝒓)(𝒌+\mathbf{})\right)`$ (181)
$`{\displaystyle \frac{1}{\alpha _1}}(𝒌+𝒓)^2\left(\overline{\alpha }_1\alpha _{\mathrm{}}(𝒌+\mathbf{})^2\alpha _2\overline{\alpha }_{1}^{}{}_{}{}^{2}\alpha _{\mathrm{}}Q^2+\overline{\alpha }_1(\alpha _1+\overline{\alpha }_1)(\mathbf{}+𝒓)(𝒌+\mathbf{})\right)`$ (182)
$`+(𝒌+\mathbf{})(𝒌+𝒓)[{\displaystyle \frac{\overline{\alpha }_1(\alpha _1+\overline{\alpha }_1)(\alpha _2+2\alpha _{\mathrm{}})}{\alpha _1\alpha _2}}(𝒌+𝒓)^2`$ (183)
$`+{\displaystyle \frac{2\overline{\alpha }_1\overline{\alpha }_2\alpha _{\mathrm{}}}{\alpha _2}}({\displaystyle \frac{1}{\alpha _1}}(𝒌+𝒓)^2+\alpha _2Q^2)`$ (184)
$`+{\displaystyle \frac{(\alpha _1+\overline{\alpha }_1)}{\alpha _2\alpha _{\mathrm{}}}}(\alpha _2(\alpha _{\mathrm{}}(𝒌+\mathbf{})^2+\overline{\alpha }_1\mathbf{}^2\overline{\alpha }_1\alpha _{\mathrm{}}Q^2+\alpha _{\mathrm{}}(\mathbf{}2𝒓)(𝒌+\mathbf{}+𝒓))`$ (185)
$`2\alpha _{\mathrm{}}(\alpha _{\mathrm{}}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\overline{\alpha }_1(\mathbf{}2𝒓)(𝒌+𝒓)))]`$ (186)
$`+(\alpha _1+\overline{\alpha }_1)((𝒌+𝒓)(𝒌+\mathbf{}+𝒓)(\mathbf{}2𝒓)(𝒌+\mathbf{})`$ (187)
$`(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)(\mathbf{}2𝒓)(𝒌+𝒓))`$ (188)
$`+((𝒌+\mathbf{})^2+\alpha _2\overline{\alpha }_1Q^2)(2\alpha _{\mathrm{}}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+\overline{\alpha }_1(\mathbf{}+𝒓)(𝒌+𝒓))]`$ (189)
$`+{\displaystyle \frac{4(𝒌+\mathbf{})(𝒌+𝒓)}{\alpha _2\overline{\alpha }_1}}𝒜_{14}^L`$ (190)
$`𝒜_{15}^T=`$ $`N_c{\displaystyle \frac{2}{\alpha _{\mathrm{}}D_1D_2D_5D_7}}\overline{\alpha }_1(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2)(𝒌+\mathbf{})(𝒌+𝒓)𝒓^2`$ (191)
$`𝒜_{22}^T=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{2}{D_2D_3D_4D_7}}[{\displaystyle \frac{(\alpha _1\overline{\alpha }_1)}{\alpha _1}}(𝒌+𝒓)^2(\overline{\alpha }_2𝒌(𝒌+\mathbf{})+\alpha _2(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓))`$ (192)
$`{\displaystyle \frac{(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2)}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2\left(\overline{\alpha }_2𝒌(𝒌+𝒓)+\alpha _2(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\right)`$ (193)
$`2(𝒌+\mathbf{})(𝒌+𝒓)(\overline{\alpha }_{2}^{}{}_{}{}^{2}𝒌(𝒌+𝒓)+\alpha _{2}^{}{}_{}{}^{2}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)`$ (194)
$`+\alpha _2\overline{\alpha }_2(𝒌(𝒌+\mathbf{})+(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)))`$ (195)
$`+ϵ((𝒌+\mathbf{})(𝒌+𝒓)((\alpha _1\overline{\alpha }_1+\alpha _2\overline{\alpha }_2)𝒌(𝒌+\mathbf{}+𝒓)2\overline{\alpha }_1\overline{\alpha }_2𝒌(𝒌+𝒓)`$ (196)
$`2\alpha _1\alpha _2(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓))+2\alpha _2\overline{\alpha }_2[(𝒌+\mathbf{})(𝒌+𝒓)]^2`$ (197)
$`+(\alpha _1\overline{\alpha }_1\alpha _2\overline{\alpha }_2)\left(𝒌(𝒌+𝒓)(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)𝒌(𝒌+\mathbf{})(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\right)`$ (198)
$`+(\alpha _1\overline{\alpha }_1\alpha _2\overline{\alpha }_2)({\displaystyle \frac{\overline{\alpha }_2}{\alpha _1}}𝒌(𝒌+\mathbf{})(𝒌+𝒓)^2{\displaystyle \frac{\alpha _2\overline{\alpha }_2}{\alpha _1\overline{\alpha }_1}}(𝒌+\mathbf{})^2(𝒌+𝒓)^2`$ (199)
$`+{\displaystyle \frac{\alpha _2}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)))]`$ (200)
$`𝒜_{23}^T=`$ $`C_F{\displaystyle \frac{4}{\overline{\alpha }_1D_2D_{4}^{}{}_{}{}^{2}D_6}}\alpha _{1}^{}{}_{}{}^{2}(1ϵ)\left(\overline{\alpha }_1𝒌^2\alpha _2(𝒌+\mathbf{})^2{\displaystyle \frac{\alpha _2\overline{\alpha }_1\mathbf{}^2}{\alpha _{\mathrm{}}}}\right)\mathbf{}(𝒌+𝒓)`$ (201)
$`+{\displaystyle \frac{2𝒌(𝒌+𝒓)}{\alpha _{1}^{}{}_{}{}^{2}}}𝒜_{23}^L`$ (202)
$`𝒜_{24}^T=`$ $`N_c{\displaystyle \frac{1}{D_2D_4D_5D_7}}[2ϵ[{\displaystyle \frac{\alpha _2(\alpha _1\overline{\alpha }_1+\alpha _2\overline{\alpha }_2)\alpha _{\mathrm{}}}{\alpha _1\overline{\alpha }_1}}(𝒌+\mathbf{})^2(𝒌+𝒓)^2`$ (203)
$`+\alpha _{\mathrm{}}𝒌(𝒌+\mathbf{})\left(2(𝒌+\mathbf{})(𝒌+𝒓)+{\displaystyle \frac{(\alpha _1\overline{\alpha }_1\alpha _2\overline{\alpha }_2)}{\alpha _1}}(𝒌+𝒓)^2\right)`$ (204)
$`\alpha _{\mathrm{}}𝒌(𝒌+𝒓)\left((𝒌+\mathbf{})^2+2(𝒌+\mathbf{})(𝒌+𝒓)+\alpha _1\overline{\alpha }_1Q^2\right)`$ (205)
$`{\displaystyle \frac{\alpha _2\alpha _{\mathrm{}}}{\overline{\alpha }_1}}(𝒌+\mathbf{})(𝒌+𝒓)((𝒌+\mathbf{})^22(𝒌+\mathbf{})(𝒌+𝒓)\alpha _1\overline{\alpha }_1Q^2)]`$ (206)
$`+(𝒌+\mathbf{})(𝒌+𝒓)({\displaystyle \frac{2\alpha _2\alpha _{\mathrm{}}((𝒌+\mathbf{})^2+\alpha _2\overline{\alpha }_1Q^2)}{\overline{\alpha }_1}}`$ (207)
$`+(\overline{\alpha }_1\alpha _1)(\alpha _2({\displaystyle \frac{(𝒌+\mathbf{})^2}{\overline{\alpha }_1}}{\displaystyle \frac{(𝒌+𝒓)^2}{\alpha _1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}}Q^2)`$ (208)
$`+{\displaystyle \frac{1}{\overline{\alpha }_1}}(2\alpha _{\mathrm{}}𝒌(𝒌+\mathbf{})+\overline{\alpha }_1𝒌(\mathbf{}2𝒓)+2\alpha _2(\mathbf{}2𝒓)(𝒌+\mathbf{}))))`$ (209)
$`+(\overline{\alpha }_1\alpha _1)(𝒌(𝒌+𝒓)(\mathbf{}2𝒓)(𝒌+\mathbf{})+𝒌(𝒌+\mathbf{})(\mathbf{}2𝒓)(𝒌+𝒓)`$ (210)
$`{\displaystyle \frac{1}{\alpha _1}}(𝒌+𝒓)^2\left(2\alpha _{\mathrm{}}𝒌(𝒌+\mathbf{})+\alpha _2(\mathbf{}+𝒓)(𝒌+\mathbf{})\right)`$ (211)
$`+{\displaystyle \frac{1}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2(\alpha _{\mathrm{}}𝒌(𝒌+𝒓)\alpha _2(\mathbf{}2𝒓)(𝒌+𝒓)))`$ (212)
$`+{\displaystyle \frac{\alpha _2}{\alpha _1}}((𝒌+\mathbf{})^2+\alpha _2\overline{\alpha }_1Q^2)(\alpha _{\mathrm{}}(𝒌+𝒓)^2+\alpha _1(\mathbf{}+𝒓)(𝒌+𝒓))]`$ (213)
$`+{\displaystyle \frac{4(𝒌+\mathbf{})(𝒌+𝒓)}{\alpha _2\overline{\alpha }_1}}𝒜_{24}^L`$ (214)
$`𝒜_{25}^T=`$ $`N_c{\displaystyle \frac{2}{\alpha _{\mathrm{}}D_2D_4D_5D_7}}\alpha _2(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2)(𝒌+\mathbf{})(𝒌+𝒓)𝒓^2`$ (215)
$`𝒜_{33}^T=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{2}{D_1D_3D_4D_6}}[(\alpha _1(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)+\overline{\alpha }_1(𝒌+𝒓)(𝒌+\mathbf{}+𝒓))`$ (216)
$`\left(\overline{\alpha }_1𝒌^2+2\alpha _1𝒌(𝒌+\mathbf{}+𝒓)+\alpha _1\overline{\alpha }_2({\displaystyle \frac{(𝒌+\mathbf{})^2}{\overline{\alpha }_1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}})\right)`$ (217)
$`(\alpha _1𝒌(𝒌+\mathbf{})+\overline{\alpha }_1𝒌(𝒌+𝒓))(2\overline{\alpha }_1𝒌(𝒌+\mathbf{}+𝒓)\alpha _1(𝒌+\mathbf{}+𝒓)^2`$ (218)
$`+\alpha _2\overline{\alpha }_1({\displaystyle \frac{(𝒌+𝒓)^2}{\alpha _1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}}))`$ (219)
$`+ϵ((\alpha _1\overline{\alpha }_1+\alpha _2\overline{\alpha }_2)𝒌(𝒌+\mathbf{}+𝒓)(𝒌+\mathbf{})(𝒌+𝒓)`$ (220)
$`+\overline{\alpha }_2\overline{\alpha }_1\left(2𝒌(𝒌+\mathbf{}+𝒓)𝒌(𝒌+𝒓)+𝒌^2(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\right)`$ (221)
$`+\overline{\alpha }_2\alpha _2\left(𝒌(𝒌+𝒓)(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)𝒌(𝒌+\mathbf{})(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\right)`$ (222)
$`+\alpha _1\alpha _2\left(2𝒌(𝒌+\mathbf{}+𝒓)(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)+𝒌(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)^2\right)`$ (223)
$`\alpha _1\overline{\alpha }_1\left(𝒌(𝒌+𝒓)(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)𝒌(𝒌+\mathbf{})(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\right)`$ (224)
$`+{\displaystyle \frac{\alpha _1\overline{\alpha }_{1}^{}{}_{}{}^{2}}{\alpha _{\mathrm{}}}}\left(\alpha _{\mathrm{}}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+\alpha _1\mathbf{}^2\right)\left({\displaystyle \frac{(𝒌+\mathbf{})^2}{\overline{\alpha }_1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}}\right)`$ (225)
$`+{\displaystyle \frac{\alpha _{1}^{}{}_{}{}^{2}\overline{\alpha }_1}{\alpha _{\mathrm{}}}}\left(\alpha _{\mathrm{}}𝒌(𝒌+\mathbf{})+\overline{\alpha }_1\mathbf{}^2\right)\left({\displaystyle \frac{(𝒌+𝒓)^2}{\alpha _1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}}\right)`$ (226)
$`\alpha _1\overline{\alpha }_1(2(𝒌(𝒌+\mathbf{}+𝒓))^2+𝒌^2(𝒌+\mathbf{}+𝒓)^2`$ (227)
$`(𝒌+\mathbf{})^2(𝒌+𝒓)^2+{\displaystyle \frac{\alpha _1\overline{\alpha }_1\mathbf{}^4}{\alpha _{\mathrm{}}^{}{}_{}{}^{2}}}))]`$ (228)
$`𝒜_{34}^T=`$ $`N_c{\displaystyle \frac{1}{D_1D_3D_5D_7}}[2ϵ(\alpha _{\mathrm{}}(\alpha _1(𝒌+\mathbf{})^2\overline{\alpha }_2(𝒌+\mathbf{})(𝒌+𝒓))(𝒌+\mathbf{}+𝒓)^2`$ (229)
$`{\displaystyle \frac{(\overline{\alpha }_1\overline{\alpha }_2)\alpha _{\mathrm{}}}{\alpha _1}}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)(𝒌+𝒓)^2`$ (230)
$`\alpha _{\mathrm{}}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)\left((𝒌+\mathbf{})^2+\overline{\alpha }_{1}^{}{}_{}{}^{2}Q^2\right)`$ (231)
$`+\overline{\alpha }_1\alpha _{\mathrm{}}({\displaystyle \frac{(𝒌+𝒓)^2}{\alpha _1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}})(\overline{\alpha }_1(𝒌+\mathbf{})(𝒌+𝒓)+\alpha _1((𝒌+\mathbf{})^2+\overline{\alpha }_1Q^2)))`$ (232)
$`\overline{\alpha }_1({\displaystyle \frac{(𝒌+𝒓)^2}{\alpha _1}}+{\displaystyle \frac{\mathbf{}^2}{\alpha _{\mathrm{}}}})(\overline{\alpha }_1\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+𝒓)`$ (233)
$`+2\alpha _1\alpha _{\mathrm{}}((𝒌+\mathbf{})^2+\overline{\alpha }_1Q^2)+\alpha _1\overline{\alpha }_1(\mathbf{}2𝒓)(𝒌+\mathbf{}))`$ (234)
$`+(𝒌+\mathbf{}+𝒓)^2(2\alpha _1\alpha _{\mathrm{}}(𝒌+\mathbf{})^2+\overline{\alpha }_2\alpha _{\mathrm{}}(𝒌+\mathbf{})(𝒌+𝒓)`$ (235)
$`+\alpha _1\overline{\alpha }_2((\mathbf{}2𝒓)(𝒌+\mathbf{})))`$ (236)
$`+3(\overline{\alpha }_1\overline{\alpha }_2)\left((𝒌+𝒓)(𝒌+\mathbf{}+𝒓)𝒓(𝒌+\mathbf{})+(𝒌+\mathbf{})(𝒌+𝒓)𝒓(𝒌+\mathbf{}+𝒓)\right)`$ (237)
$`+((𝒌+\mathbf{})^2+\overline{\alpha }_{1}^{}{}_{}{}^{2}Q^2)\left(\alpha _{\mathrm{}}(𝒌+𝒓)(𝒌+\mathbf{}+𝒓)+\alpha _1(\mathbf{}+𝒓)(𝒌+\mathbf{}+𝒓)\right)`$ (238)
$`+{\displaystyle \frac{(\overline{\alpha }_1\overline{\alpha }_2)}{\alpha _1\overline{\alpha }_1\alpha _{\mathrm{}}}}(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)(2\overline{\alpha }_1\alpha _{\mathrm{}}^{}{}_{}{}^{2}(𝒌+𝒓)^2+\alpha _{1}^{}{}_{}{}^{2}\overline{\alpha }_1\mathbf{}^2`$ (239)
$`\alpha _{1}^{}{}_{}{}^{2}\alpha _{\mathrm{}}((𝒌+\mathbf{})^2+\overline{\alpha }_1Q^2)\alpha _1\overline{\alpha }_1\alpha _{\mathrm{}}((𝒌+\mathbf{})(𝒌+𝒓)+\mathbf{}(𝒌+𝒓)))]`$ (240)
$`+{\displaystyle \frac{4(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)}{\overline{\alpha }_{1}^{}{}_{}{}^{2}}}𝒜_{34}^L`$ (241)
$`𝒜_{35}^T=`$ $`N_c{\displaystyle \frac{2}{\alpha _{\mathrm{}}D_1D_3D_5D_7}}\alpha _1(\overline{\alpha }_{1}^{}{}_{}{}^{2}+\overline{\alpha }_{2}^{}{}_{}{}^{2})(𝒌+\mathbf{})(𝒌+\mathbf{}+𝒓)𝒓^2`$ (242)
$`𝒜_{44}^T=`$ $`N_c{\displaystyle \frac{1}{D_2D_{5}^{}{}_{}{}^{2}D_7}}[4\alpha _{\mathrm{}}^{}{}_{}{}^{2}((𝒌+\mathbf{})(𝒌+𝒓))^2+{\displaystyle \frac{(\alpha _1+\overline{\alpha }_1)\alpha _{\mathrm{}}^{}{}_{}{}^{2}}{\alpha _1\overline{\alpha }_1}}(𝒌+\mathbf{})^2(𝒌+𝒓)^2`$ (243)
$`+4ϵ(4\alpha _{\mathrm{}}^{}{}_{}{}^{2}((𝒌+\mathbf{})(𝒌+𝒓))^2{\displaystyle \frac{(\alpha _1\overline{\alpha }_2)\alpha _{\mathrm{}}^{}{}_{}{}^{2}}{\alpha _1\overline{\alpha }_1}}(𝒌+\mathbf{})^2(𝒌+𝒓)^2`$ (244)
$`+\alpha _{\mathrm{}}^{}{}_{}{}^{2}Q^2\left(\alpha _1(𝒌+\mathbf{})^2+\overline{\alpha }_1(𝒌+𝒓)^2+\alpha _1\overline{\alpha }_1Q^2\right)`$ (245)
$`{\displaystyle \frac{\alpha _{\mathrm{}}^{}{}_{}{}^{2}}{\alpha _1\overline{\alpha }_1}}(𝒌+\mathbf{})(𝒌+𝒓)(\alpha _1(\alpha _1+\alpha _2)(𝒌+\mathbf{})^2+\overline{\alpha }_1(\overline{\alpha }_1+\overline{\alpha }_2)(𝒌+𝒓)^2`$ (246)
$`\alpha _1\overline{\alpha }_1(\alpha _1+\overline{\alpha }_1)Q^2))`$ (247)
$`+{\displaystyle \frac{\overline{\alpha }_1\alpha _{\mathrm{}}}{\alpha _1}}(𝒌+𝒓)^2\left(\alpha _2\alpha _{\mathrm{}}Q^23(\mathbf{}𝒓)(𝒌+\mathbf{})\right)`$ (248)
$`+{\displaystyle \frac{\alpha _1\alpha _{\mathrm{}}}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2\left(\overline{\alpha }_2\alpha _{\mathrm{}}Q^2+3(\mathbf{}𝒓)(𝒌+𝒓)\right)`$ (249)
$`+(𝒌+\mathbf{})(𝒌+𝒓)({\displaystyle \frac{2\alpha _{\mathrm{}}(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2+2(\alpha _1\alpha _2)\alpha _{\mathrm{}})}{\overline{\alpha }_1}}(𝒌+\mathbf{})^2`$ (250)
$`{\displaystyle \frac{2\alpha _{\mathrm{}}(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2+2(\overline{\alpha }_1\overline{\alpha }_2)\alpha _{\mathrm{}})}{\alpha _1}}(𝒌+𝒓)^2`$ (251)
$`+4(\alpha _2+\overline{\alpha }_2)\alpha _{\mathrm{}}^{}{}_{}{}^{2}Q^2+2(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2)\left(\alpha _{\mathrm{}}Q^2+2𝒓^2+2(\mathbf{}𝒓)^2\right)`$ (252)
$`6(\alpha _1\alpha _2)\alpha _{\mathrm{}}(\mathbf{}𝒓)(𝒌+\mathbf{})+6(\overline{\alpha }_1\overline{\alpha }_2)\alpha _{\mathrm{}}(\mathbf{}𝒓)(𝒌+𝒓))`$ (253)
$`+3\alpha _2\overline{\alpha }_1\alpha _{\mathrm{}}Q^2(\mathbf{}𝒓)(𝒌+𝒓)3\alpha _1\overline{\alpha }_2\alpha _{\mathrm{}}Q^2(\mathbf{}𝒓)(𝒌+\mathbf{})]`$ (254)
$`𝒜_{45}^T=`$ $`N_c{\displaystyle \frac{2}{D_2D_{5}^{}{}_{}{}^{2}D_7}}(\alpha _2\overline{\alpha }_1+\alpha _1\overline{\alpha }_2)(𝒌+\mathbf{})(𝒌+𝒓)𝒓^2`$ (255)
$`𝒜_{55}^T=`$ $`\mathrm{\hspace{0.33em}0}`$ (256)
At several places we have expressed the matrix elements $`𝒜𝒜_{(ij)}^T`$, $`𝒜_{(ij)}^T`$ in terms of the corresponding matrix elements for the longitudinally polarized photon; we therefore list here also the expressions for longitudinal polarizations :

$`𝒜𝒜_{11}^L=`$ $`C_F{\displaystyle \frac{2}{D_1^2D_2^2}}(1ϵ)\alpha _1\overline{\alpha }_1𝑨_1^2`$ (257)
$`𝒜𝒜_{12}^L=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{1}{D_1D_2^2D_4}}\left((1ϵ)\alpha _1\alpha _2𝑨_1𝑨_2\alpha _1\alpha _2^2\overline{\alpha }_1𝒓^2\right)`$ (258)
$`𝒜𝒜_{13}^L=`$ $`C_F{\displaystyle \frac{2}{D_1^2D_2D_3}}\left((1ϵ)\overline{\alpha }_1^2𝑨_1𝑨_3+\overline{\alpha }_1^2\left(𝒌+𝒓\right)\left(𝒌+\mathbf{}+𝒓\right)\right)`$ (259)
$`𝒜𝒜_{14}^L=`$ $`N_c{\displaystyle \frac{1}{D_1D_2^2D_5}}({\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\alpha _1(\alpha _2\overline{\alpha }_1+2(1ϵ)\alpha _{\mathrm{}}^2)𝑨_1𝑨_4`$ (260)
$`+{\displaystyle \frac{1}{2}}\alpha _1\alpha _2\overline{\alpha }_1^2D_2+{\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\alpha _1\alpha _2^2\overline{\alpha }_1𝑨_1𝒓\alpha _1\alpha _2^2\overline{\alpha }_1𝒓^2)`$ (261)
$`𝒜𝒜_{15}^L=`$ $`N_c{\displaystyle \frac{1}{\alpha _{\mathrm{}}D_1D_2^2D_5}}\alpha _1\alpha _2^2\overline{\alpha }_1^2𝒓^2`$ (262)
$`𝒜𝒜_{22}^L=`$ $`C_F{\displaystyle \frac{2}{\overline{\alpha }_1D_2^2D_4^2}}(1ϵ)\alpha _1\alpha _2^2𝑨_2^2`$ (263)
$`𝒜𝒜_{23}^L=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{1}{D_1D_2D_3D_4}}\left((1ϵ)\alpha _2\overline{\alpha }_1𝑨_2𝑨_3+\alpha _2\overline{\alpha }_1\left(𝒌+𝒓\alpha _1𝒓\right)\left(𝒌+\mathbf{}+\overline{\alpha }_1𝒓\right)\right)`$ (264)
$`𝒜𝒜_{24}^L=`$ $`N_c{\displaystyle \frac{1}{D_2^2D_4D_5}}({\displaystyle \frac{1}{2\alpha _{\mathrm{}}\overline{\alpha }_1}}\alpha _1\alpha _2(\alpha _2\overline{\alpha }_1+2(1ϵ)\alpha _{\mathrm{}}^2)𝑨_2𝑨_4`$ (265)
$`+{\displaystyle \frac{1}{2}}\alpha _1\alpha _2^2\overline{\alpha }_1D_2+{\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\alpha _1\alpha _2^2\overline{\alpha }_1𝑨_2𝒓+\alpha _1\alpha _2^2\overline{\alpha }_1𝒓^2)`$ (266)
$`𝒜𝒜_{25}^L=`$ $`N_c{\displaystyle \frac{1}{\alpha _{\mathrm{}}D_2^2D_4D_5}}\alpha _1\alpha _2^3\overline{\alpha }_1𝒓^2`$ (267)
$`𝒜𝒜_{33}^L=`$ $`C_F{\displaystyle \frac{2}{\alpha _1D_1^2D_3^2}}(1ϵ)\overline{\alpha }_1^3𝑨_3^2`$ (268)
$`𝒜𝒜_{34}^L=`$ $`N_c{\displaystyle \frac{1}{D_1D_2D_3D_5}}({\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\overline{\alpha }_1(\alpha _1\overline{\alpha }_1\alpha _{\mathrm{}}2(1ϵ)\alpha _{\mathrm{}}^2)𝑨_3𝑨_4`$ (269)
$`{\displaystyle \frac{1}{2}}\alpha _1\overline{\alpha }_1^2\overline{\alpha }_2D_2+\alpha _1\alpha _2\overline{\alpha }_1^2D_7+{\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\alpha _1\alpha _2\overline{\alpha }_1^2𝑨_3𝒓+\alpha _1\alpha _2\overline{\alpha }_1^2𝒓^2`$ (270)
$`+{\displaystyle \frac{3}{2}}\alpha _2\overline{\alpha }_1𝑨_3(𝒌+\mathbf{}))`$ (271)
$`𝒜𝒜_{35}^L=`$ $`N_c{\displaystyle \frac{1}{\alpha _{\mathrm{}}D_1D_2D_3D_5}}\alpha _1\alpha _2\overline{\alpha }_1^2\overline{\alpha }_2𝒓^2`$ (272)
$`𝒜𝒜_{44}^L=`$ $`N_c{\displaystyle \frac{1}{D_2^2D_5^2}}\left({\displaystyle \frac{2}{\overline{\alpha }_1}}𝑨_4^2\alpha _1\left((1ϵ)\alpha _{\mathrm{}}^2+\alpha _2\overline{\alpha }_1\right)+\alpha _1\alpha _2\alpha _{\mathrm{}}\overline{\alpha }_1D_2+2\alpha _1\alpha _2^2\overline{\alpha }_1𝒓^2\right)`$ (273)
$`𝒜𝒜_{45}^L=`$ $`N_c{\displaystyle \frac{1}{D_2^2D_5^2}}\alpha _1\alpha _2^2\overline{\alpha }_1𝒓^2`$ (274)
$`𝒜𝒜_{55}^L=`$ $`\mathrm{\hspace{0.33em}0}`$ (275)
$`𝒜_{11}^L=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{1}{D_1D_2D_6D_7}}\left((1ϵ)\alpha _1\overline{\alpha }_1𝑨_1𝑩_1+\alpha _1\overline{\alpha }_1\left(𝒌+\alpha _1𝒓\right)\left(𝒌+\mathbf{}+\overline{\alpha }_2𝒓\right)\right)`$ (276)
$`𝒜_{12}^L=`$ $`C_F{\displaystyle \frac{2}{D_1D_2D_3D_7}}\left((1ϵ)\overline{\alpha }_1\overline{\alpha }_2𝑨_1𝑩_2+\overline{\alpha }_1\overline{\alpha }_2\left(𝒌+𝒓\right)\left(𝒌+\mathbf{}+𝒓\right)\right)`$ (277)
$`𝒜_{13}^L=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{1}{D_1D_2D_4D_6}}\left((1ϵ)\alpha _1^2𝑨_1𝑩_3\alpha _1^2\alpha _2\overline{\alpha }_1𝒓^2\right)`$ (278)
$`𝒜_{14}^L=`$ $`N_c{\displaystyle \frac{1}{D_1D_2D_5D_7}}({\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\overline{\alpha }_1(\alpha _1\overline{\alpha }_1\alpha _{\mathrm{}}(1ϵ)\mathrm{\hspace{0.17em}2}\alpha _{\mathrm{}}^2)𝑨_1𝑩_4`$ (279)
$`+\alpha _1\overline{\alpha }_1^2\overline{\alpha }_2D_2{\displaystyle \frac{1}{2}}\alpha _1\alpha _2\overline{\alpha }_1^2D_7+{\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\alpha _1\alpha _2\overline{\alpha }_1\overline{\alpha }_2𝑨_1𝒓\alpha _1\alpha _2\overline{\alpha }_1\overline{\alpha }_2𝒓^2`$ (280)
$`{\displaystyle \frac{3}{2}}\overline{\alpha }_1\overline{\alpha }_2𝑨_1(𝒌+𝒓))`$ (281)
$`𝒜_{15}^L=`$ $`N_c{\displaystyle \frac{1}{\alpha _{\mathrm{}}D_1D_2D_5D_7}}\alpha _1\alpha _2\overline{\alpha }_1^2\overline{\alpha }_2𝒓^2`$ (282)
$`𝒜_{22}^L=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{1}{D_2D_3D_4D_7}}\left((1ϵ)\alpha _2\overline{\alpha }_2𝑨_2𝑩_2+\alpha _2\overline{\alpha }_2\left(𝒌+𝒓\alpha _1𝒓\right)\left(𝒌+\mathbf{}+\overline{\alpha }_1𝒓\right)\right)`$ (283)
$`𝒜_{23}^L=`$ $`C_F{\displaystyle \frac{2}{\overline{\alpha }_1D_2D_4^2D_6}}(1ϵ)\alpha _1^2\alpha _2𝑨_2𝑩_3`$ (284)
$`𝒜_{24}^L=`$ $`N_c{\displaystyle \frac{1}{D_2D_4D_5D_7}}({\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\alpha _2(\alpha _1\overline{\alpha }_1\alpha _{\mathrm{}}2(1ϵ)\alpha _{\mathrm{}}^2)𝑨_2𝑩_4`$ (285)
$`+\alpha _1\alpha _2\overline{\alpha }_1\overline{\alpha }_2D_2{\displaystyle \frac{1}{2}}\alpha _1\alpha _2^2\overline{\alpha }_1D_7+{\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\alpha _1\alpha _2\overline{\alpha }_1\overline{\alpha }_2𝑨_2𝒓+\alpha _1\alpha _2\overline{\alpha }_1\overline{\alpha }_2𝒓^2`$ (286)
$`{\displaystyle \frac{3}{2}}\alpha _2\overline{\alpha }_2𝑨_2(𝒌+𝒓))`$ (287)
$`𝒜_{25}^L=`$ $`N_c{\displaystyle \frac{1}{\alpha _{\mathrm{}}D_2D_4D_5D_7}}\alpha _1\alpha _2^2\overline{\alpha }_1\overline{\alpha }_2𝒓^2`$ (288)
$`𝒜_{33}^L=`$ $`{\displaystyle \frac{1}{N_c}}{\displaystyle \frac{1}{D_1D_3D_4D_6}}\left((1ϵ)\alpha _1\overline{\alpha }_1𝑨_3𝑩_3+\alpha _1\overline{\alpha }_1\left(𝒌+𝒓\alpha _1𝒓\right)\left(𝒌+\mathbf{}+\overline{\alpha }_1𝒓\right)\right)`$ (289)
$`𝒜_{34}^L=`$ $`N_c{\displaystyle \frac{1}{D_1D_3D_5D_7}}({\displaystyle \frac{1}{2\alpha _1\alpha _{\mathrm{}}}}\overline{\alpha }_1^2(\alpha _1\overline{\alpha }_2+2(1ϵ)\alpha _{\mathrm{}}^2)𝑨_3𝑩_4`$ (290)
$`+{\displaystyle \frac{1}{2}}\alpha _1\overline{\alpha }_1^2\overline{\alpha }_2D_7+{\displaystyle \frac{1}{2\alpha _{\mathrm{}}}}\alpha _1\overline{\alpha }_1^2\overline{\alpha }_2𝑨_3𝒓+\alpha _1\overline{\alpha }_1^2\overline{\alpha }_2𝒓^2)`$ (291)
$`𝒜_{35}^L=`$ $`N_c{\displaystyle \frac{1}{\alpha _{\mathrm{}}D_1D_3D_5D_7}}\alpha _1\overline{\alpha }_1^2\overline{\alpha }_2^2𝒓^2`$ (292)
$`𝒜_{44}^L=`$ $`N_c{\displaystyle \frac{1}{D_2D_5^2D_7}}((2\alpha _1\overline{\alpha }_1+\alpha _{\mathrm{}}2(1ϵ)\alpha _{\mathrm{}}^2)𝑨_4𝑩_4`$ (293)
$`+{\displaystyle \frac{1}{2}}\alpha _1\alpha _{\mathrm{}}\overline{\alpha }_1\overline{\alpha }_2D_2+{\displaystyle \frac{1}{2}}\alpha _1\alpha _2\alpha _{\mathrm{}}\overline{\alpha }_1D_7+2\alpha _1\alpha _2\overline{\alpha }_1\overline{\alpha }_2𝒓^2\alpha _{\mathrm{}}^2(𝒌+\mathbf{})(𝒌+𝒓))`$ (294)
$`𝒜_{45}^L=`$ $`N_c{\displaystyle \frac{1}{D_2D_5^2D_7}}\alpha _1\alpha _2\overline{\alpha }_1\overline{\alpha }_2𝒓^2`$ (295)
$`𝒜_{55}^L=`$ $`\mathrm{\hspace{0.33em}0}.`$ (296)
Here we have used the following abbreviations:

$`𝑨_1=`$ $`\alpha _2\mathbf{}+\alpha _{\mathrm{}}(𝒌+𝒓)`$ (297)
$`𝑨_2=`$ $`\overline{\alpha }_1\mathbf{}+\alpha _{\mathrm{}}(𝒌+\mathbf{})`$ (298)
$`𝑨_3=`$ $`\alpha _1\mathbf{}\alpha _{\mathrm{}}(𝒌+𝒓)`$ (299)
$`𝑨_4=`$ $`\alpha _2(\mathbf{}𝒓)+\alpha _{\mathrm{}}(𝒌+𝒓)`$ (300)
$`𝑩_1=`$ $`\overline{\alpha }_2\mathbf{}\alpha _{\mathrm{}}(𝒌+\mathbf{})`$ (301)
$`𝑩_2=`$ $`\alpha _1\mathbf{}\alpha _{\mathrm{}}(𝒌+𝒓)`$ (302)
$`𝑩_3=`$ $`\overline{\alpha }_1\mathbf{}+\alpha _{\mathrm{}}(𝒌+\mathbf{})`$ (303)
$`𝑩_4=`$ $`\overline{\alpha }_2(\mathbf{}𝒓)\alpha _{\mathrm{}}(𝒌+\mathbf{}).`$ (304)
We conclude our listing of the squared matrix elements with a few remarks. First, similar to the longitudinal case, the sum of the squared matrix elements for the transverse photon is ultraviolet finite: individual terms diverge as $`\mathbf{}`$ becomes large (for fixed $`𝒌`$), but in the sum these divergences cancel. Secondly, it is easy to show that the factorization properties in configuration space which for the longitudinal photon have been discussed in section V of  also hold for the transverse photon: this supports the validity of the photon wave function picture beyond the LO approximation. Finally, the limit of small $`\alpha _{\mathrm{}}`$ agrees with the LO BFKL calculation; we shall return to this question further below.

## V Singular Limits of the Real Corrections

In order to formulate finite expressions for the $`q\overline{q}g`$ contribution to the photon impact factor we examine the divergent limits of the squared matrix element $`|_{q\overline{q}g}|^2`$ for the process $`\gamma ^{}+qq\overline{q}g+q`$. For the phase space of the $`q\overline{q}g`$-final state we have from (11), (12):

$$d\varphi _{q\overline{q}g}\frac{sd\beta _r}{2\pi }=\frac{d\alpha }{2\alpha (1\alpha \alpha _{\mathrm{}})}\frac{d^{D2}𝒌}{(2\pi )^{D1}}\frac{d\alpha _{\mathrm{}}}{2\alpha _{\mathrm{}}}\frac{d^{D2}\mathbf{}}{(2\pi )^{D1}}$$
(305)
In the following we suppress the superscript $`(D2)`$ in the integration measure of the transverse momenta. The contribution of the $`q\overline{q}g`$ state to the impact factor reads:

$$\frac{d\alpha }{2\alpha (1\alpha \alpha _{\mathrm{}})}\frac{d𝒌}{(2\pi )^{D1}}\frac{d\alpha _{\mathrm{}}}{2\alpha _{\mathrm{}}}\frac{d\mathbf{}}{(2\pi )^{D1}}|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2,$$
(306)
where the square of the vertex function, $`|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2`$, follows from (77) and (78).

We begin by examining the limit in which the gluon is either collinear to the quark or to the antiquark of the $`q\overline{q}g`$-system. In the latter case, the gluon momentum $`\mathrm{}`$ is collinear with the momentum $`qk\mathrm{}`$ of the antiquark, resulting in $`D_4=0`$ or

$$\mathbf{}^{}=\mathbf{}+\frac{\alpha _{\mathrm{}}}{1\alpha }𝒌=0.$$
(307)
We approximate our results for $`|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0),a}|^2`$ in the region around the singularity, i.e. we calculate the residue of the pole at $`D_4=0`$. The collinear limit of QCD matrix elements is well known and can be taken from the literature . As an additional check, we have calculated the divergent limits directly from our matrix elements. For both polarizations of the photon, the $`g\overline{q}`$ collinear limit of (306) can be written as

$`|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2d\varphi _{q\overline{q}g}{\displaystyle \frac{sd\beta _r}{2\pi }}|^{\mathrm{c}oll,\overline{q}}=|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}|^2{\displaystyle \frac{d\alpha }{2\alpha (1\alpha )}}{\displaystyle \frac{d𝒌}{(2\pi )^{32ϵ}}}`$ (308)
$`\times C_FP_{gq}({\displaystyle \frac{\alpha _{\mathrm{}}}{1\alpha }},ϵ){\displaystyle \frac{d\alpha _{\mathrm{}}}{(1\alpha )}}{\displaystyle \frac{d\mathbf{}}{(2\pi )^{32ϵ}\mathbf{}_{}^{\mathbf{}}{}_{}{}^{2}}}`$ (309)
$`=_2(\alpha ,𝒌)d\alpha d𝒌C_FP_{gq}({\displaystyle \frac{\alpha _{\mathrm{}}}{1\alpha }},ϵ){\displaystyle \frac{d\alpha _{\mathrm{}}}{(1\alpha )}}{\displaystyle \frac{d\mathbf{}}{(2\pi )^{32ϵ}\mathbf{}_{}^{\mathbf{}}{}_{}{}^{2}}}`$ (310)
where $`_2(\alpha ,𝒌)_2(\alpha ,𝒌;𝒓,Q)`$ was defined in (30), (31), (32), and $`C_FP_{gq}(z,ϵ)`$ denotes the $`qg`$ Altarelli-Parisi splitting function in $`D=42ϵ`$ dimensions:

$$P_{gq}(z,ϵ)=\frac{1+(1z)^2}{z}ϵz.$$
(311)
The color coefficient $`C_F`$ indicates that, for the examination of singular limits it will be convenient to use the color basis $`C_F`$ and $`C_A`$ rather than the color factors $`N_c`$, $`C_F`$, and $`1/N_c`$ which appear in the matrix elements of the previous section. In this basis, the $`C_A`$ part is free from collinear singularities.

In the case where the gluon is collinear to the quark, the infrared limit is obtained from (310) by substituting<sup>*</sup><sup>*</sup>*This is nothing but the $`q\overline{q}`$ variables exchange at fixed gluon momentum. $`\alpha `$ by $`\overline{\alpha }=1\alpha \alpha _{\mathrm{}}`$ and $`𝒌`$ by $`\overline{𝒌}=𝒌𝒓\mathbf{}`$. The “soft-collinear” limit of (310) that is needed later is obtained by retaining from the splitting function only the pole $`2/z`$:

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2d\varphi _{q\overline{q}g}\frac{sd\beta _r}{2\pi }|^{\mathrm{s}oft,coll,\overline{q}}=_2(\alpha ,𝒌)d\alpha d𝒌\mathrm{\hspace{0.17em}2}C_F\frac{d\alpha _{\mathrm{}}}{\alpha _{\mathrm{}}}\frac{d\mathbf{}}{(2\pi )^{32ϵ}\mathbf{}_{}^{\mathbf{}}{}_{}{}^{2}}.$$
(312)
Next, we examine the ‘soft limit’ where the gluon momentum $`\alpha _{\mathrm{}}0`$ and $`\mathbf{}0`$. Formally we substitute $`\mathrm{}\rho \mathrm{}`$ or $`\mathbf{}\rho \mathbf{},\alpha _{\mathrm{}}\rho \alpha _{\mathrm{}}`$, expand around $`\rho =0`$ and keep only the leading term. As in the collinear case, the result is proportional to the Born approximation, for both polarizations of the photon. We separate the result according to different colour structures:

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2|^{\mathrm{soft}}=C_A|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2|_{C_A}^{\mathrm{soft}}+C_F|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2|_{C_F}^{\mathrm{soft}}$$
(313)
with

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2|_{C_A}^{\mathrm{soft}}=|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}|^2\frac{4\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)}{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)^2\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)^2}$$
(314)
and

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2|_{C_F}^{\mathrm{soft}}=|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}|^2\frac{4\alpha _{\mathrm{}}^2(𝒂+𝒃)^2}{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)^2\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)^2},$$
(315)
where, for convenience, we have defined

$$𝒂=\frac{𝒌+𝒓}{\alpha },𝒃=\frac{𝒌}{1\alpha }.$$
(316)
In the soft limit, the phase space measure (305) factorizes in the product of the $`q\overline{q}`$ measure (28) and differentials of the gluon variables

$$d\varphi _{q\overline{q}g}\frac{sd\beta _r}{2\pi }\left(d\varphi _{q\overline{q}}\frac{sd\beta _r}{2\pi }\right)\left(\frac{d\alpha _{\mathrm{}}}{2\alpha _{\mathrm{}}}\frac{d\mathbf{}}{(2\pi )^{D1}}\right).$$
(317)
Our result for the soft limit agrees with . We also note that, in the limit $`\mathbf{}^{}=\mathbf{}+\alpha _{\mathrm{}}𝒃0`$, we reproduce from (315) and (317) the “soft-collinear” limit (312).

Finally, we have to consider the limit in which the gluon is emitted in the central rapidity region. Formally, this corresponds to the singular limit $`\alpha _{\mathrm{}}0`$: the matrix elements listed in the previous section are finite for $`\alpha _{\mathrm{}}0`$, but the phase space element (12) introduces a logarithmic singularity. Consequently, when doing the subsequent $`\alpha _{\mathrm{}}`$ integration this divergence $`1/\alpha _{\mathrm{}}`$ will produce the leading-$`\mathrm{ln}s`$ logarithm. As discussed in section II, this part of the gluon phase space belongs to final states with two large rapidity gaps and has to be subtracted from our result. This will be the topic of section VI. For the moment we only state that in the limit $`\alpha _{\mathrm{}}0`$, for both polarizations of the photon, the sums of our matrix elements split into the product of the Born approximation and the BFKL kernel:

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2|_{\mathrm{LL}}=|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}|^2|_{\begin{array}{c}𝒓𝒓\mathbf{}\hfill \\ 𝒌𝒌+\mathbf{}\hfill \end{array}}\frac{1}{\left(\mathbf{}𝒓\right)^4}\frac{4C_A𝒓^2\left(\mathbf{}𝒓\right)^2}{\mathbf{}^2}.$$
(318)
Here the last factor is the BFKL kernel (without the gluon trajectory function)

$$𝒦_{\mathrm{r}𝒶𝓁}(𝒓\mathbf{},𝒓)=\frac{C_A}{\pi }\frac{𝒓^2\left(\mathbf{}𝒓\right)^2}{\mathbf{}^2},$$
(319)
and $`1/\left(\mathbf{}𝒓\right)^4`$ are the propagators of the (reggeized) gluons. Note that in the Born level vertex functions the transverse momenta are shifted from $`𝒓`$ to $`𝒓\mathbf{}`$ and from $`𝒌`$ to $`𝒌+\mathbf{}`$. The color factor $`C_A`$ signals that the $`C_F`$ part of our matrix elements does not contribute to our subtraction of the central region. In the limit where the centrally emitted gluon is soft ($`\mathbf{}𝒓`$), we obtain

$$|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2|_{\mathrm{LL}}^{\mathrm{s}oft}=|\mathrm{\Gamma }_{\gamma ^{}q\overline{q}}^{(0)}(\alpha ,𝒌)|^2\frac{4C_A}{\mathbf{}^2}.$$
(320)
The same result is obtained from the soft limit (314), if we go into the central region limit $`\alpha _{\mathrm{}}1`$.

## VI Subtraction of the Central Region

After having calculated the various divergent limits of our matrix elements, we now turn to the subtraction of the central region. As we have discussed in section II, contributions to the total cross section are divided according to the number of large rapidity gaps. What we have presented in the last two sections are the results for the process $`\gamma ^{}+q(q\overline{q}g)+q`$, imposing the constraint of a large rapidity gap between the $`(q\overline{q}g)`$ system and the quark. These results therefore contain, as a special case, still the configuration where the $`q\overline{q}g`$ system contains a second rapidity gap between the $`q\overline{q}`$ pair and the gluon: this piece has to be removed since it counts as a configuration with two gaps and, in order $`\alpha _s^3`$ of the total cross section, it entirely belongs to the leading-$`\mathrm{ln}s`$ approximation (16).

In order to make the separation more explicit, we first divide the phase space of the produced gluon into ’upper’ and ’lower’ halves. Introducing

$$\alpha _{\mathrm{}}^{cut}=\frac{|\mathbf{}|}{\sqrt{s}},$$
(321)
the ’upper’ region $`\alpha _{\mathrm{}}^{cut}<\alpha _{\mathrm{}}<1\alpha `$ consists of the ’upper’ half of the central region plus the fragmentation region of the incoming photon, whereas the ’lower’ region $`\alpha _{\mathrm{}}<\alpha _{\mathrm{}}^{cut}`$ contains the lower half of the central region plus the fragmentation region of the quark. In the upper region we write:

$`{\displaystyle |\mathrm{\Gamma }_{\gamma ^{}q\overline{q}g}^{(0)}|^2𝑑\varphi _{q\overline{q}g}\frac{sd\beta _r}{2\pi }\mathrm{\Theta }(\alpha _{\mathrm{}}\alpha _{\mathrm{}}^{\mathrm{c}ut})}`$ (322)
$`={\displaystyle 𝑑𝒌𝑑\mathbf{}_0^1𝑑\alpha _0^{1\alpha }𝑑\alpha _{\mathrm{}}_3(\alpha _{\mathrm{}},\mathbf{};\alpha ,𝒌;𝒓,Q)\mathrm{\Theta }(\alpha _{\mathrm{}}\alpha _{\mathrm{}}^{\mathrm{c}ut})}`$ (323)
Here $`_3`$ denotes the full integrand involving the sum of the matrix elements of the three-particle final state.

The division into central region and fragmentation region is closely related to the choice of the energy scale $`s_0`$. In (16) we have introduced the energy scale $`s_0`$: since a change in $`s_0`$ can be absorbed into the NLO impact factors $`\mathrm{\Phi }_A^{(1)}`$ and $`\mathrm{\Phi }_B^{(1)}`$, one might expect that there exists some freedom in choosing this scale. We adopt the simplest choice, namely a constant value which is independent of kinematic variables inside the impact factors. Examples for $`s_0`$ include the (negative) virtuality of the $`t`$-channel gluon, $`𝒓^2`$, or the (negative) mass of the photon, $`Q^2`$.

In order to define the $`q\overline{q}g`$ contribution to the photon impact factor, $`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{real})}`$, we have to subtract, starting from (322), the leading-$`\mathrm{ln}s`$ piece with the energy scale $`s_0`$. This can be achieved by considering the central region limit (318) of the matrix element which implies

$$_3(\alpha _{\mathrm{}},\mathbf{};\alpha ,𝒌;𝒓,Q)|_{\mathrm{L}L}=_2(\alpha ,𝒌+\mathbf{};𝒓𝒍,Q)\frac{1}{(2\pi )^{22ϵ}(\mathbf{}𝒓)^4}𝒦_{\mathrm{r}𝒶𝓁}(𝒓\mathbf{},𝒓)\frac{1}{\alpha _{\mathrm{}}}.$$
(324)
We recognize the integrand $`_2`$ occurring in the LO impact factor (30) with the value of the reggeon momentum being shifted to $`𝒓𝒓\mathbf{}`$ (since the momentum $`\mathbf{}`$ is carried by the outgoing gluon), the propagators associated to the $`t`$-channel (reggeized) gluons, the measure factor $`1/(2\pi )^{22ϵ}`$ required by the definition (2), the real part of the BFKL kernel, and, finally, a factor $`1/\alpha _{\mathrm{}}`$. It is this factor $`1/\alpha _{\mathrm{}}`$ which, when the integration over $`\alpha _{\mathrm{}}`$ is done, provides the logarithm of the energy. Our aim of finding the LL piece with energy scale $`s_0`$ leads to the following $`\alpha _{\mathrm{}}`$ integral:

$$_0^{\alpha _{\mathrm{}}^{\mathrm{m}ax}}\frac{d\alpha _{\mathrm{}}}{\alpha _{\mathrm{}}}\mathrm{\Theta }(\alpha _{\mathrm{}}\alpha _{\mathrm{}}^{\mathrm{c}ut})=\frac{1}{2}\mathrm{log}\frac{(\alpha _{\mathrm{}}^{\mathrm{m}ax})^2s}{\mathbf{}^2}.$$
(325)
with $`\alpha _{\mathrm{}}^{\mathrm{m}ax}|\mathbf{}|/\sqrt{s_0}`$. Therefore, we define as genuine real emission contribution to the photon impact factor the expression

$`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{r}eal)}`$ $`={\displaystyle }d𝒌d\mathbf{}{\displaystyle _0^1}d\alpha {\displaystyle _{\alpha _{\mathrm{}}^{\mathrm{c}ut}}^{\mathrm{}}}d\alpha _{\mathrm{}}[_3(\alpha _{\mathrm{}},\mathbf{};\alpha ,𝒌;𝒓,Q)\mathrm{\Theta }(1\alpha \alpha _{\mathrm{}})`$ (326)
$`_2(\alpha ,𝒌+\mathbf{};𝒓𝒍,Q){\displaystyle \frac{1}{(2\pi )^{22ϵ}(\mathbf{}𝒓)^4}}𝒦_{\mathrm{r}𝒶𝓁}(𝒓\mathbf{},𝒓){\displaystyle \frac{1}{\alpha _{\mathrm{}}}}\mathrm{\Theta }({\displaystyle \frac{|\mathbf{}|}{\sqrt{s_0}}}\alpha _{\mathrm{}})].`$ (327)
In this definition of our impact factor, the following remark is in place. With our choice of a constant energy scale, which has led to the subtraction, we are adding (in the LL term) and subtracting (in $`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{r}eal)}`$) a region of the phase space which is kinematically forbidden. In fact, the LL subtraction (last line of Eq. (327)) includes the region of phase space where $`\alpha _{\mathrm{}}>1\alpha `$. It is important to make sure that this “forbidden” region $`|\mathbf{}|>(1\alpha )\sqrt{s_0},\mathrm{\hspace{0.33em}1}\alpha <\alpha _{\mathrm{}}<|\mathbf{}|/\sqrt{s_0}`$ contributes a finite term and does not introduce new divergences. The finiteness results from the “good” UV behaviour of the LO impact factor, namely $`\mathrm{\Phi }_\gamma ^{}^{(0)}(𝒓\mathbf{})\text{const}`$ for $`\mathbf{}\mathrm{}`$ (up to logarithms). After multiplication with the propagator $`1/\mathbf{}^4`$ and with the factor $`\mathrm{log}[|\mathbf{}|/(1\alpha )\sqrt{s_0}]`$ (which is due to the $`\alpha _{\mathrm{}}`$ integral), the resulting $`\mathbf{}`$ integral is finite. This situation is quite different from the case of partonic impact factors , where $`\mathrm{\Phi }_{q,g}^{(0)}(𝒓\mathbf{})\mathbf{}^2`$, and the $`\mathbf{}`$ integral, after including the gluon propagators and the logarithmic factor, leads to a double pole in $`1/ϵ`$.

## VII Finite Combinations

In the previous section we have suggested a suitable form of the LL subtraction and thus specified the real emission contribution to the impact factor. What remains are the infrared divergences which have to cancel once we combine real and virtual corrections. In this final section of our paper we show explicitly the cancellation of all the IR singularities, and we present expressions consisting of well defined finite terms.

Given that all the singular terms of the virtual corrections are proportional to the LO impact factor (30), as one can see from Eq. (33), it is convenient to study the singular regions of the phase space integrals of the real corrections by integrating first only over the gluon variables. In fact, the gluon integration generates singular terms which are in turn proportional to the LO impact factor. They will be combined with the singular virtual corrections after a proper identification of the remaining integration variables.

It is also convenient to separate the different colour structure contributions into a $`C_A`$ term and a $`C_F`$ term

$$_3C_A_3^{C_A}+C_F_3^{C_F},$$
(328)
because the LL subtraction is contained only in the former term, whilst the collinear singularities are present only in the latter. Soft singularities are found in both terms.

### A $`C_A`$ term

Let us consider the $`C_A`$ contribution to the real emission impact factor corrections, as specified in Eq. (327). We want to compute first the $`\mathbf{}`$ and $`\alpha _{\mathrm{}}`$ integrals at fixed $`𝒌`$ and $`\alpha `$:

$`{\displaystyle \frac{d\mathrm{\Phi }_\gamma ^{}^{(1)}}{d\alpha d𝒌}}|_{C_A}={\displaystyle 𝑑\mathbf{}_{\alpha _{\mathrm{}}^{\mathrm{c}ut}}^{\mathrm{}}𝑑\alpha _{\mathrm{}}}`$ $`[_3^{C_A}(\alpha _{\mathrm{}},\mathbf{};\alpha ,𝒌;𝒓,Q)\mathrm{\Theta }(1\alpha \alpha _{\mathrm{}})`$ (329)
$`_2(\alpha ,𝒌+\mathbf{};𝒓𝒍,Q){\displaystyle \frac{2}{(2\pi )^{32ϵ}}}{\displaystyle \frac{𝒓^2}{(𝒓\mathbf{})^2\mathbf{}^2}}{\displaystyle \frac{1}{\alpha _{\mathrm{}}}}\mathrm{\Theta }({\displaystyle \frac{|\mathbf{}|}{\sqrt{s_0}}}\alpha _{\mathrm{}})].`$ (330)
This integral is divergent in the soft region $`\alpha _{\mathrm{}}0,\mathbf{}0`$ when $`ϵ=0`$. We evaluate it in dimensional regularization by means of the subtraction method: we subtract from the integrand its soft approximation, in such a way that the resulting integral will be $`ϵ`$-finite. Subsequently we re-add what we have subtracted, and we integrate this term analytically. Of course, the soft limit is not uniquely defined outside the soft region. We choose the soft approximation according to Eqs. (314, 320, 317); we extend the soft subtraction up to the kinematic limit of the gluon longitudinal momentum, i.e., $`\alpha _{\mathrm{}}<1\alpha `$:

$`{\displaystyle \frac{d\mathrm{\Phi }_\gamma ^{}^{(1)}}{d\alpha d𝒌}}|_{C_A}^{\mathrm{s}oft}`$ $`=_2(\alpha ,𝒌;𝒓,Q){\displaystyle \frac{2}{(2\pi )^{32ϵ}}}`$ (331)
$`\times {\displaystyle }d\mathbf{}{\displaystyle _{\alpha _{\mathrm{}}^{\mathrm{c}ut}}^{\mathrm{}}}{\displaystyle \frac{d\alpha _{\mathrm{}}}{\alpha _{\mathrm{}}}}[{\displaystyle \frac{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)}{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)^2\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)^2}}{\displaystyle \frac{1}{\mathbf{}^2}}\mathrm{\Theta }({\displaystyle \frac{|\mathbf{}|}{\sqrt{s_0}}}\alpha _{\mathrm{}})]\mathrm{\Theta }(1\alpha \alpha _{\mathrm{}}).`$ (332)
Accordingly, we define the finite $`C_A`$ part of the real corrections to the impact factor as

$`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{r}eal)}|_{C_A}^{\mathrm{f}inite}{\displaystyle }d𝒌d\mathbf{}{\displaystyle _0^1}d\alpha {\displaystyle _0^{\mathrm{}}}d\alpha _{\mathrm{}}\{[_3^{C_A}(\alpha _{\mathrm{}},\mathbf{};\alpha ,𝒌;𝒓,Q)\mathrm{\Theta }(1\alpha \alpha _{\mathrm{}})`$ (333)
$`_2(\alpha ,𝒌+\mathbf{};𝒓𝒍,Q){\displaystyle \frac{2}{(2\pi )^{32ϵ}}}{\displaystyle \frac{𝒓^2}{(𝒓\mathbf{})^2\mathbf{}^2}}{\displaystyle \frac{1}{\alpha _{\mathrm{}}}}\mathrm{\Theta }({\displaystyle \frac{|\mathbf{}|}{\sqrt{s_0}}}\alpha _{\mathrm{}})]`$ (334)
$`_2(\alpha ,𝒌;𝒓,Q){\displaystyle \frac{2}{(2\pi )^{32ϵ}}}{\displaystyle \frac{1}{\alpha _{\mathrm{}}}}[{\displaystyle \frac{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)}{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)^2\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)^2}}{\displaystyle \frac{1}{\mathbf{}^2}}\mathrm{\Theta }({\displaystyle \frac{|\mathbf{}|}{\sqrt{s_0}}}\alpha _{\mathrm{}})]\mathrm{\Theta }(1\alpha \alpha _{\mathrm{}})\}.`$ (335)
Note that we have shifted to zero the lower limit of integration of $`\alpha _{\mathrm{}}`$ because, thanks to the LL subtraction, there is no contribution to the integral coming from the central region. Note also that the singularities of the integrand due to the denominators $`\mathbf{}\alpha _{\mathrm{}}𝒂`$ and $`\mathbf{}+\alpha _{\mathrm{}}𝒃`$ are integrable, because of the scalar product in the numerator. This reflects the fact that the $`C_A`$ part is free of collinear singularities.

The expression written in Eq. (332) can be computed analytically. Since the lower limit of integration in $`\alpha _{\mathrm{}}`$ depends on the energy $`s`$, one gets an $`s`$ dependent result. However, because of the LL-soft subtraction to the soft term, the central region contribution is suppressed, and we expect a finite expression in the $`s\mathrm{}`$ limit. We can compute the integrals of the two terms in Eq. (332) in this limit ($`\alpha _{\mathrm{}}^{\mathrm{c}ut}0`$) separately because dimensional regularization promotes the divergent $`\mathrm{log}s`$ contributions into UV $`ϵ`$-poles, which will cancel in the difference. We have:

$$J_1𝑑\mathbf{}_0^{1\alpha }\frac{d\alpha _{\mathrm{}}}{\alpha _{\mathrm{}}}\frac{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)}{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)^2\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)^2}=c_\mathrm{\Gamma }\pi ^{1ϵ}\frac{(1\alpha )^{2ϵ}}{2ϵ^2}\left((𝒂+𝒃)^2\right)^ϵ$$
(336)
and

$$J_2𝑑\mathbf{}_0^{1\alpha }\frac{d\alpha _{\mathrm{}}}{\alpha _{\mathrm{}}}\frac{1}{\mathbf{}^2}\mathrm{\Theta }(\frac{|\mathbf{}|}{\sqrt{s_0}}\alpha _{\mathrm{}})=\frac{\pi ^{1ϵ}}{\mathrm{\Gamma }(1ϵ)}\frac{(1\alpha )^{2ϵ}}{2ϵ^2}s_0^ϵ.$$
(337)
These two integrals have to be multiplied by $`_2`$ and integrated over $`\alpha `$ and $`𝒌`$. After Eq. (30) we have noticed that $`_2`$ is symmetric under the change of variables $`\alpha 1\alpha `$, $`𝒌\mathbf{}𝒌`$, corresponding to the $`q\overline{q}`$ exchange; the term $`(𝒂+𝒃)^2=M^2/\alpha (1\alpha )`$ is also invariant under this exchange. Therefore, on the right-hand-side of Eqs. (336) and (337), we can replace $`(1\alpha )^{2ϵ}\frac{1}{2}\alpha ^{2ϵ}+\frac{1}{2}(1\alpha )^{2ϵ}`$, which casts the expressions into explicitly symmetric forms w.r.t. $`\alpha 1\alpha `$. This manipulation will be useful when combining real and virtual corrections.

In conclusion, the divergent $`C_A`$ part of the real corrections to the impact factor is

$`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{r}eal)}|_{C_A}^{\mathrm{d}ivergent}`$ $`{\displaystyle 𝑑𝒌_0^1𝑑\alpha _2(\alpha ,𝒌;𝒓,Q)\frac{2}{(2\pi )^{32ϵ}}(J_1J_2)}`$ (338)
$`={\displaystyle }d𝒌{\displaystyle _0^1}d\alpha _2(\alpha ,𝒌;𝒓,Q){\displaystyle \frac{c_\mathrm{\Gamma }}{(4\pi )^{2ϵ}}}\{{\displaystyle \frac{2}{ϵ}}[\mathrm{log}M^2\mathrm{log}s_0\mathrm{log}\alpha (1\alpha )]`$ (339)
$`+[(\mathrm{log}s_0+\mathrm{log}\alpha (1\alpha ))^2\mathrm{log}^2M^2]+𝒪(ϵ)\}.`$ (340)
For consistency, the variables $`\alpha `$ and $`𝒌`$ of the above equation have to be identified with the analogous variables of Eq. (42).

### B $`C_F`$ term

Let us now consider the $`C_F`$ contribution to the real emission impact factor corrections. The first thing to note is that, since the BFKL kernel is proportional to $`C_A`$, the LL subtraction of Eq. (327) does not affect the $`C_F`$ part. Therefore, the whole $`C_F`$ part of $`_3`$ contributes to the impact factor. In addition, we can shift the lower limit $`\alpha _{\mathrm{}}^{\mathrm{c}ut}`$ for the $`\alpha _{\mathrm{}}`$ integral to zero, just because the central region does not contribute in this case.

The $`C_F`$ term contains both soft and collinear divergences. We separate the singular contributions by means of the subtraction method. In what follows, the last two parameters of $`_2`$ are always $`𝒓,Q`$. Therefore we simplify the notation by indicating $`_2(\alpha ,𝒌)_2(\alpha ,𝒌;𝒓,Q)`$.

We start by identifying the soft singular term. This is done by subtracting from $`_3^{C_F}`$ its soft limit and then re-adding it again. We remark again that, beyond the soft region, the limit of integration is not uniquely defined. In practice, by using the expression (315), we adopt an asymmetric function with respect to the $`q\overline{q}`$ exchange at fixed non-vanishing gluon momentum:

$`{\displaystyle \frac{d\mathrm{\Phi }_\gamma ^{}^{(1)}}{d\alpha d𝒌}}|_{C_F}^{\mathrm{s}oft}`$ $`=_2(\alpha ,𝒌){\displaystyle \frac{2(𝒂+𝒃)^2}{(2\pi )^{32ϵ}}}{\displaystyle 𝑑\mathbf{}_0^{1\alpha }𝑑\alpha _{\mathrm{}}\frac{\alpha _{\mathrm{}}}{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)^2\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)^2}}`$ (341)
$`=_2(\alpha ,𝒌){\displaystyle \frac{c_\mathrm{\Gamma }}{(4\pi )^{2ϵ}}}{\displaystyle \frac{2}{ϵ^2}}\left[\alpha ^{2ϵ}+(1\alpha )^{2ϵ}\right]\left[{\displaystyle \frac{M^2}{\alpha (1\alpha )}}\right]^ϵ,`$ (342)
where in the last equality a symmetrization in $`\alpha `$, as indicated after Eq. (337), has been performed.

Both the full $`C_F`$ term $`_3^{C_F}`$ and its soft limit (342) have a divergent behaviour when the gluon is either collinear to the outgoing quark or antiquark. Also for the evaluation of the collinear divergences we adopt the subtraction method. For the case of the gluon being collinear to the antiquark ($`g\overline{q}`$), we derive the limiting expression of the original integrand (i.e., without soft subtraction yet) from Eq. (310), which yields ($`\mathbf{}^{}=\mathbf{}+\alpha _{\mathrm{}}𝒃`$)

$`{\displaystyle \frac{d\mathrm{\Phi }_\gamma ^{}^{(1)}}{d\alpha d𝒌}}|_{C_F}^{\mathrm{c}oll,\overline{q}}`$ $`=_2(\alpha ,𝒌){\displaystyle \frac{d\mathbf{}}{(2\pi )^{32ϵ}}_0^{1\alpha }\frac{d\alpha _{\mathrm{}}}{(1\alpha )}P_{gq}(\frac{\alpha _{\mathrm{}}}{1\alpha },ϵ)\frac{\mathrm{\Theta }(\alpha _{\mathrm{}}\mathrm{\Lambda }|\mathbf{}^{\mathbf{}}|)}{\mathbf{}_{}^{\mathbf{}}{}_{}{}^{2}}}`$ (343)
$`=_2(\alpha ,𝒌){\displaystyle \frac{\alpha ^{2ϵ}+(1\alpha )^{2ϵ}}{(4\pi )^{2ϵ}\mathrm{\Gamma }(1ϵ)}}{\displaystyle \frac{\mathrm{\Lambda }^{2ϵ}}{ϵ}}\left[{\displaystyle \frac{1}{ϵ}}+{\displaystyle \frac{2}{12ϵ}}{\displaystyle \frac{1}{2}}\right],`$ (344)
where we have decided to perform the collinear subtraction only in a cone — specified by the $`\mathrm{\Theta }`$ function — whose axis coincides with the gluon momentum being parallel to the antiquark momentum, and the vertex sits on the soft point. The “opening angle” of the cone is parametrized by a cutoff $`\mathrm{\Lambda }`$: the bigger $`\mathrm{\Lambda }`$, the wider the domain of the collinear subtraction.

In order to end up with a finite integral, we have to subtract also the collinear singularities of the soft subtraction term (342). For the $`g\overline{q}`$ we have (see Eq. (312):

$`{\displaystyle \frac{d\mathrm{\Phi }_\gamma ^{}^{(1)}}{d\alpha d𝒌}}|_{C_F}^{\mathrm{s}oft,coll,\overline{q}}`$ $`=_2(\alpha ,𝒌){\displaystyle \frac{2}{(2\pi )^{32ϵ}}}{\displaystyle _0^{1\alpha }}{\displaystyle \frac{d\alpha _{\mathrm{}}}{\alpha _{\mathrm{}}}}{\displaystyle 𝑑\mathbf{}\frac{\mathrm{\Theta }(\alpha _{\mathrm{}}\mathrm{\Lambda }|\mathbf{}^{\mathbf{}}|)}{\mathbf{}_{}^{\mathbf{}}{}_{}{}^{2}}}.`$ (345)
$`=_2(\alpha ,𝒌){\displaystyle \frac{1}{(4\pi )^{2ϵ}\mathrm{\Gamma }(1ϵ)}}\left[\alpha ^{2ϵ}+(1\alpha )^{2ϵ}\right]{\displaystyle \frac{\mathrm{\Lambda }^{2ϵ}}{ϵ^2}}.`$ (346)
The soft-collinear subtraction uses the same cone as defined for the collinear one. It is important that Eq. (346) can be obtained both as a collinear limit of Eq. (342) and as a soft limit of Eq. (344). This guarantees not only that the collinear singularities introduced by the soft subtraction are subtracted, but also that the soft singularities introduced by the collinear subtraction are properly removed.

Note that the collinear term (344) is proportional to $`_2(\alpha ,𝒌)`$. This is a consequence of the convenient choice of variable employed to parametrize the outgoing particles. In fact, when the gluon is collinear to the antiquark, the $`g\overline{q}`$ system has to be treated as a single entity carrying the sum of quantum numbers and momenta of the constituting particles. In this case, the $`g\overline{q}`$ longitudinal momentum fraction and transverse momentum are $`1\alpha `$ and $`𝒌`$ respectively. The integration over the gluon variables can therefore be independently performed leaving the LO integrand $`_2`$ factored out.

The treatment of the $`gq`$ collinear singularity can be performed in the same way as for the $`g\overline{q}`$ collinear case. There are only two small differences: first of all, the collinear limit of the integrand $`_3`$ is proportional to $`_2(\alpha +\alpha _{\mathrm{}},𝒌+\mathbf{})`$, which involves explicitly the gluon variables $`\alpha _{\mathrm{}}`$ and $`\mathbf{}`$. This is due to our set of variables that does not describe the $`qg`$ system and the single gluon independently. A more suitable set of variable in this case is simply obtained by changing the $`q\overline{q}`$ labels:

$`\overline{\alpha }`$ $`1\alpha \alpha _{\mathrm{}}`$ (347a)
$`\overline{𝒌}`$ $`𝒓𝒌\mathbf{}`$ (347b)
$`\mathbf{}^{\prime \prime }`$ $`\mathbf{}\alpha _{\mathrm{}}𝒂.`$ (347c)
From Eq. (344) one gets

$`{\displaystyle \frac{d\mathrm{\Phi }_\gamma ^{}^{(1)}}{d\overline{\alpha }d\overline{𝒌}}}|_{C_F}^{\mathrm{c}oll,q}`$ $`=_2(\overline{\alpha },\overline{𝒌}){\displaystyle \frac{d\mathbf{}}{(2\pi )^{32ϵ}}_0^{1\overline{\alpha }}\frac{d\alpha _{\mathrm{}}}{(1\overline{\alpha })}P_{gq}(\frac{\alpha _{\mathrm{}}}{1\overline{\alpha }},ϵ)}`$ (348)
$`\times \left({\displaystyle \frac{1\overline{\alpha }}{1\overline{\alpha }\alpha _{\mathrm{}}}}\right)^2{\displaystyle \frac{\mathrm{\Theta }(\alpha _{\mathrm{}}\mathrm{\Lambda }\frac{1\overline{\alpha }\alpha _{\mathrm{}}}{1\overline{\alpha }}|\mathbf{}^{\prime \prime }|)}{\mathbf{}_{}^{\prime \prime }{}_{}{}^{2}}}`$ (349)
$`=_2(1\overline{\alpha },𝒓\overline{𝒌}){\displaystyle \frac{\overline{\alpha }^{2ϵ}+(1\overline{\alpha })^{2ϵ}}{(4\pi )^{2ϵ}\mathrm{\Gamma }(1ϵ)}}{\displaystyle \frac{\mathrm{\Lambda }^{2ϵ}}{ϵ}}\left[{\displaystyle \frac{1}{ϵ}}+{\displaystyle \frac{2}{12ϵ}}{\displaystyle \frac{1}{2}}\right].`$ (350)
The additional factors in the second line on the rhs of Eq. (350) (when compared to the corresponding line of Eq. (344)) are due to the fact that

$$\frac{}{}\mathbf{}+\alpha _{\mathrm{}}𝒃|_{\begin{array}{c}\alpha \overline{\alpha }\hfill \\ 𝒌\overline{𝒌}\hfill \end{array}}=\frac{1\overline{\alpha }\alpha _{\mathrm{}}}{1\overline{\alpha }}(\mathbf{}\alpha _{\mathrm{}}𝒂).$$
(351)
The second difference is that, considering the soft-collinear subtraction, we cannot use the soft limit of Eq. (350) at fixed $`\overline{\alpha }`$ and $`\overline{𝒌}`$, because this does not correspond to the collinear limit of the soft approximation (342). The right way of subtracting the soft singularity of (350) and at the same time the $`gq`$ collinear singularity of (342) is to take the collinear limit of the latter:

$`{\displaystyle \frac{d\mathrm{\Phi }_\gamma ^{}^{(1)}}{d\alpha d𝒌}}|_{C_F}^{\mathrm{s}oft,coll,q}`$ $`=_2(\alpha ,𝒌){\displaystyle \frac{2}{(2\pi )^{32ϵ}}}{\displaystyle _0^{1\alpha }}{\displaystyle \frac{d\alpha _{\mathrm{}}}{\alpha _{\mathrm{}}}}{\displaystyle 𝑑\mathbf{}\frac{\mathrm{\Theta }(\alpha _{\mathrm{}}\mathrm{\Lambda }|\mathbf{}^{\prime \prime }|)}{\mathbf{}_{}^{\prime \prime }{}_{}{}^{2}}}.`$ (352)
$`=_2(\alpha ,𝒌){\displaystyle \frac{1}{(4\pi )^{2ϵ}\mathrm{\Gamma }(1ϵ)}}\left[\alpha ^{2ϵ}+(1\alpha )^{2ϵ}\right]{\displaystyle \frac{\mathrm{\Lambda }^{2ϵ}}{ϵ^2}}.`$ (353)
In conclusion, we write the finite part of the $`C_F`$ term of the impact factor by subtracting, from $`_3^{C_F}`$, the unintegrated expressions of Eqs. (342, 344, 346, 350, 353):

$`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{r}eal)}|_{C_F}^{\mathrm{f}inite}{\displaystyle }d𝒌d\mathbf{}{\displaystyle _0^1}d\alpha {\displaystyle _0^{1\alpha }}d\alpha _{\mathrm{}}\{_3^{C_F}(\alpha _{\mathrm{}},\mathbf{};\alpha ,𝒌;𝒓,Q)`$ (354)
$`_2(\alpha ,𝒌){\displaystyle \frac{2(𝒂+𝒃)^2}{(2\pi )^{32ϵ}}}{\displaystyle \frac{\alpha _{\mathrm{}}}{\left(\mathbf{}\alpha _{\mathrm{}}𝒂\right)^2\left(\mathbf{}+\alpha _{\mathrm{}}𝒃\right)^2}}`$ (355)
$`_2(\alpha ,𝒌){\displaystyle \frac{1}{(2\pi )^{32ϵ}}}\left[{\displaystyle \frac{1}{1\alpha }}P_{gq}({\displaystyle \frac{\alpha _{\mathrm{}}}{1\alpha }},ϵ){\displaystyle \frac{2}{\alpha _{\mathrm{}}}}\right]{\displaystyle \frac{\mathrm{\Theta }(\alpha _{\mathrm{}}\mathrm{\Lambda }|\mathbf{}+\alpha _{\mathrm{}}𝒃|)}{(\mathbf{}+\alpha _{\mathrm{}}𝒃)^2}}`$ (356)
$`[_2(\alpha +\alpha _{\mathrm{}},𝒌+\mathbf{}){\displaystyle \frac{\alpha +\alpha _{\mathrm{}}}{\alpha ^2}}P_{gq}({\displaystyle \frac{\alpha _{\mathrm{}}}{\alpha +\alpha _{\mathrm{}}}},ϵ)\mathrm{\Theta }(\alpha _{\mathrm{}}\mathrm{\Lambda }{\displaystyle \frac{\alpha }{\alpha +\alpha _{\mathrm{}}}}|\mathbf{}\alpha _{\mathrm{}}𝒂|)`$ (357)
$`_2(\alpha ,𝒌){\displaystyle \frac{2}{\alpha _{\mathrm{}}}}\mathrm{\Theta }(\alpha _{\mathrm{}}\mathrm{\Lambda }|\mathbf{}\alpha _{\mathrm{}}𝒂|)]{\displaystyle \frac{1}{(2\pi )^{32ϵ}(\mathbf{}\alpha _{\mathrm{}}𝒂)^2}}.`$ (358)
The divergent part of the $`C_F`$ term is obtained by re-adding the subtracted pieces in their integrated form. In doing this, we identify the $`1\overline{\alpha }`$ and $`𝒓\overline{𝒌}`$ variables of Eq. (350) with the variables $`\alpha `$ and $`𝒌`$ of Eqs. (344, 346, 353) (and later with $`\alpha `$ and $`𝒌`$ of Eq. (42)), so that

$`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{r}eal)}|_{C_F}^{\mathrm{d}ivergent}`$ $`={\displaystyle }d𝒌{\displaystyle _0^1}d\alpha _2(\alpha ,𝒌){\displaystyle \frac{c_\mathrm{\Gamma }}{(4\pi )^{2ϵ}}}\{{\displaystyle \frac{4}{ϵ^2}}+{\displaystyle \frac{1}{ϵ}}[64\mathrm{log}M^2]`$ (359)
$`+2[8+\mathrm{log}^2M^23\mathrm{log}\mathrm{\Lambda }^23\mathrm{log}\alpha (1\alpha )+\mathrm{log}^2{\displaystyle \frac{\alpha }{1\alpha }}]+𝒪(ϵ)\}.`$ (360)
The final expression for the one-loop correction to the photon impact factor can be obtained by summing the real contributions given in Eqs. (334, 340, 358, 360) to the virtual corrections $`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{v}irtual)}`$ of . One can check the cancellation of soft and collinear singularities by combining the singular piece of the virtual corrections (42) with the $`ϵ`$-poles of Eqs. (340) and (360). The cancellation of double and single pole is straightforward.

In conclusion, the final expression of the impact factor at NLO can be obtained by summing the finite part of the virtual corrections, the scale dependent term (49) and the extra term (46) as discussed in section III, the finite terms of Eqs. (340) and (360), and the integrals in Eqs. (334) and (358) (evaluated at $`ϵ=0`$):

$`\mathrm{\Phi }_\gamma ^{}^{(1)}=`$ $`\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{v}irtual)}|^{\mathrm{f}inite}{\displaystyle \frac{2\mathrm{\Phi }_\gamma ^{}^{(0)}}{(4\pi )^2}}\{\beta _0\mathrm{ln}{\displaystyle \frac{𝒓^2}{\mu ^2}}+C_F\mathrm{ln}(𝒓^2)\}`$ (361)
$`+{\displaystyle \frac{1}{(4\pi )^2}}{\displaystyle }d𝒌{\displaystyle _0^1}d\alpha _2(\alpha ,𝒌)\{C_A[\mathrm{ln}^2\alpha (1\alpha )s_0\mathrm{ln}^2M^2]`$ (362)
$`+2C_F[83\mathrm{ln}\alpha (1\alpha )\mathrm{\Lambda }^2+\mathrm{ln}^2M^2+\mathrm{ln}^2{\displaystyle \frac{\alpha }{1\alpha }}]\}`$ (363)
$`+C_A\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{r}eal)}|_{C_A}^{\mathrm{f}inite}+C_F\mathrm{\Phi }_\gamma ^{}^{(1,\mathrm{r}eal)}|_{C_F}^{\mathrm{f}inite}.`$ (364)
## VIII Conclusions

In this third part of our program of calculating the NLO corrections to the photon impact factor we first have presented a list of the real corrections, the $`q\overline{q}g`$ intermediate state inside the impact factor. This completes the results of the paper  which contained the longitudinal photon only. We then have identified and computed the divergent parts of the phase space integral of the real corrections: in addition to the infrared singularities which are due to collinear and soft configurations of the produced gluon there is a logarithmic divergence related to the $`\mathrm{ln}s`$-piece of the central region. Removal of the latter part of the real corrections introduces the energy scale $`s_0`$ which represents a central element of the NLO corrections. Finally, we have combined the real corrections with the divergent parts of the virtual corrections: these infrared finite combinations are the main results of this paper.

With these calculations we have essentially completed the analytic part of our program. What remains are numerical steps: the phase space integrals both for the $`q\overline{q}`$ state (which includes the finite pieces of the virtual corrections) and the $`q\overline{q}g`$ state. Certain parts of these integrals can be done analytically, but, as a general strategy, we define standard integrals and express our integrands in terms of these expressions. The final evaluations have to be done with the computer.

Once this is solved technically, the photon impact factor allows us to compute the total cross section for $`\gamma ^{}\gamma ^{}`$-scattering to next-to-leading logarithmic accuracy in $`s`$. This will be the first consistent calculation of a cross section to this order since the NLO corrections to the BFKL kernel were first published . Taking advantage of the recently obtained impact factor for forward jets  in the same limit we may apply the photon impact factor also to the calculation of the forward jet cross section at HERA.

We have emphasized already before, that interest in the NLO calculation of the photon impact factor comes from various directions. A prominent example is the photon wave function picture. In our analysis  of the real corrections we have extracted the new $`q\overline{q}g`$ Fock component of the photon wave function which, most conveniently, is expressed in configuration space. In the present paper we have studied the cancellation of divergences, but this was done on momentum space. It will therefore be necessary to translate these calculations into configuration space. Furthermore, in order to obtain a clearer geometrical picture of how color charge is distributed inside the virtual photon, we need to take a closer look into the various pieces found in  and in this paper. Another question of interest of the photon wave function picture is the form of the NLO correction to the $`q\overline{q}`$ Fock component and of the $`q\overline{q}`$ dipole cross section: do they destroy some of the main features of the LO calculations, for example the conservation of the transverse dipole size during the interaction with the target? Answers to these questions are to be found in our virtual corrections, but they require further theoretical efforts.

###### Acknowledgements.

Helpful conversations with C.-F. Qiao are gratefully acknowledged. S.G. wishes to thank the II. Institut für Theoretische Physik of the Universität Hamburg where large parts of this work were done for kind hospitality.