text:
# 1 Introduction

## 1 Introduction

The Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation  and the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation  are used now for a theoretical description of structure functions of the deep-inelastic $`ep`$ scattering (DIS) at small values of the Bjorken variable $`x`$. The structure functions are measured by the H1 and ZEUS collaborations . The higher order QCD corrections to the splitting kernels of the DGLAP equation are well known . On the other hand, the calculation of the next-to-leading order (NLO) corrections to the BFKL kernel was completed comparatively recently .

In supersymmetric gauge theories the structure of the BFKL and DGLAP equations is significantly simplified. In the case of an extended $`N=4`$ SUSY the NLO corrections to the BFKL equation were calculated in ref.  for an arbitrary value of the conformal spin $`|n|`$ in a framework of the dimensional regularization (DREG) scheme. Below in the section 3 these results are presented in the dimensional reduction (DRED) scheme  which does not violate the supersymmetry. The analyticity of the eigenvalue of the BFKL kernel as a function of the conformal spin $`|n|`$ gives a possibility to relate the DGLAP and BFKL equations in this model, as we show below (see also ). Further, in both schemes (DREG and DRED) this eigenvalue has the property of the hermitian separability which is similar to the holomorphic separability (see Section 3).

Let us introduce the unintegrated parton distributions (UnPD) $`\phi _a(x,k_{}^2)`$ (hereafter $`a=q,g,\phi `$ for the spinor, vector and scalar particles, respectively) and the (integrated) parton distributions (PD) $`f_a(x,Q^2)`$

$$f_a(x,Q^2)=_{k_{}^2<Q^2}𝑑k_{}^2\phi _a(x,k_{}^2).$$

In DIS $`Q^2=q^2`$ and $`x=Q^2/(2pq)`$ are the Bjorken variables, $`k_{}`$ is the transverse component of the parton momentum and $`q`$ and $`p`$ are the virtual photon and hadron momenta, respectively.

The DGLAP equation has a clear probabilistic interpretation and relates the parton densities with different values of $`Q^2`$:

$`{\displaystyle \frac{d}{d\mathrm{ln}Q^2}}f_a(x,Q^2)`$ $`=`$ $`\stackrel{~}{W}_af_a(x,Q^2)+{\displaystyle _x^1}{\displaystyle \frac{dy}{y}}{\displaystyle \underset{b}{}}\stackrel{~}{W}_{ba}(x/y)f_b(y,Q^2),`$
$`\stackrel{~}{W}_a`$ $`=`$ $`{\displaystyle \underset{b}{}}{\displaystyle _0^1}𝑑xx\stackrel{~}{W}_{ab}(x),`$ (1)
where the second term in the r.h.s. of the equation is the Mellin convolution of the transition probabilities $`\stackrel{~}{W}_{ba}(x)`$ and PD $`f_b(x,Q^2)`$. Usually the first and second terms are unified with the use of the splitting kernel $`W_{ba}`$ :

$`{\displaystyle \frac{d}{d\mathrm{ln}Q^2}}f_a(x,Q^2)={\displaystyle _x^1}{\displaystyle \frac{dy}{y}}{\displaystyle \underset{b}{}}W_{ba}(x/y)f_b(y,Q^2).`$ (2)
It is well known, that equation (2) is simplified after its Mellin transformation to the Lorentz spin $`j`$ representation:

$$\frac{d}{d\mathrm{ln}Q^2}f_a(j,Q^2)=\underset{b}{}\gamma _{ab}(j)f_b(j,Q^2),$$

where

$$f_a(j,Q^2)=_0^1𝑑xx^{j1}f_a(x,Q^2)$$

are the Mellin momenta of parton distributions. The Mellin moment of the splitting kernel

$$\gamma _{ab}(j)=_0^1𝑑xx^{j1}W_{ba}(x)$$

coincides with the anomalous dimension matrix for the twist-2 operators 111As in Ref. , the anomalous dimensions differ from those used usually in the description of DIS by a factor $`(1/2)`$, i.e. $`\gamma _{ab}(j)=(1/2)\gamma _{ab}^{DIS}(j)`$.. These operators are constructed as bilinear combinations of the fields which describe corresponding partons (see Eq. (3) below).

On the other hand, the BFKL equation relates the unintegrated gluon distributions with small values of the Bjorken variable $`x`$:

$$\frac{d}{d\mathrm{ln}(1/x)}\phi _g(x,k_{}^2)=2\omega (k_{}^2)\phi _g(x,k_{}^2)+d^2k_{}^{}K(k_{},k_{}^{})\phi _g(x,k_{}^2),$$

where $`\omega (k_{}^2)<0`$ is the gluon Regge trajectory .

Let us introduce the local gauge-invariant twist-two operators:

$`O_{\mu _1,\mathrm{},\mu _j}^g`$ $`=`$ $`\widehat{S}G_{\rho \mu _1}D_{\mu _2}D_{\mu _3}\mathrm{}D_{\mu _{j1}}G_{\rho \mu _j},`$
$`\stackrel{~}{O}_{\mu _1,\mathrm{},\mu _j}^g`$ $`=`$ $`\widehat{S}G_{\rho \mu _1}D_{\mu _2}D_{\mu _3}\mathrm{}D_{\mu _{j1}}\stackrel{~}{G}_{\rho \mu _j},`$
$`O_{\mu _1,\mathrm{},\mu _j}^q`$ $`=`$ $`\widehat{S}\overline{\mathrm{\Psi }}\gamma _{\mu _1}D_{\mu _2}\mathrm{}D_{\mu _j}\mathrm{\Psi },`$
$`\stackrel{~}{O}_{\mu _1,\mathrm{},\mu _j}^q`$ $`=`$ $`\widehat{S}\overline{\mathrm{\Psi }}\gamma _5\gamma _{\mu _1}D_{\mu _2}\mathrm{}D_{\mu _j}\mathrm{\Psi },`$
$`O_{\mu _1,\mathrm{},\mu _j}^\phi `$ $`=`$ $`\widehat{S}\mathrm{\Phi }D_{\mu _1}D_{\mu _2}\mathrm{}D_{\mu _j}\mathrm{\Phi },`$ (3)
where the spinor $`\mathrm{\Psi }`$ and field tensor $`G_{\rho \mu }`$ describe gluinos and gluons, respectively. The last expression is constructed from the covariant derivatives $`D_\mu `$ of the scalar field $`\mathrm{\Phi }`$ appearing in extended supersymmetric models. The symbol $`\widehat{S}`$ implies a symmetrization of the tensor in the Lorentz indices $`\mu _1,\mathrm{},\mu _j`$ and a subtraction of its traces.

The matrix elements of $`O_{\mu _1,\mathrm{},\mu _j}^a`$ and $`\stackrel{~}{O}_{\mu _1,\mathrm{},\mu _j}^a`$ are related to the moments of the distributions $`f_a(x,Q^2)`$ and $`\mathrm{\Delta }f_a(x,Q^2)`$ for unpolarized and polarized partons in a hadron $`h`$ in the following way

$`{\displaystyle _0^1}𝑑xx^{j1}f_a(x,Q^2)`$ $`=`$ $`<h|\stackrel{~}{n}^{\mu _1}\mathrm{}\stackrel{~}{n}^{\mu _j}O_{\mu _1,\mathrm{},\mu _j}^a|h>,a=(q,g,\phi ),`$
$`{\displaystyle _0^1}𝑑xx^{j1}\mathrm{\Delta }f_a(x,Q^2)`$ $`=`$ $`<h|\stackrel{~}{n}^{\mu _1}\mathrm{}\stackrel{~}{n}^{\mu _j}\stackrel{~}{O}_{\mu _1,\mathrm{},\mu _j}^a|h>,a=(q,g).`$ (4)
Here the vector $`\stackrel{~}{n}^\mu `$ is light-like: $`\stackrel{~}{n}^2=0`$. In the deep-inelastic $`ep`$ scattering we have $`\stackrel{~}{n}^\mu q^\mu +xp^\mu `$.

The conformal spin $`|n|`$ and the quantity $`1+\omega `$ ($`\omega `$ is an eigenvalue of the BFKL kernel) coincide respectively with total numbers of the transversal and longitudinal Lorentz indices of the tensor $`O_{\mu _1,\mathrm{},\mu _J}^a`$ with the rank

$$J=1+\omega +|n|.$$
(5)
Namely, we can introduce the following projectors of this tensor

$`\stackrel{~}{n}^{\mu _1}\mathrm{}\stackrel{~}{n}^{\mu _{1+\omega }}O_{\mu _1,\mathrm{},\mu _{1+\omega },\sigma _1,\mathrm{},\sigma _{|n|}}^al_{}^{\sigma _1}\mathrm{}l_{}^{\sigma _{|n|}},`$ (6)
where the complex transverse vector $`l_{}`$ is given below

$$l_{}^\sigma =\frac{1}{\sqrt{2}}\left(\delta _1^\sigma +i\delta _2^\sigma \right),l_{}^2=l_{}p=l_{}q=0.$$

It is important, that the anomalous dimension matrices $`\gamma _{ab}(j)`$ and $`\stackrel{~}{\gamma }_{ab}(j)`$ for the twist-2 operators $`O_{\mu _1,\mathrm{},\mu _j}^a`$ and $`\stackrel{~}{O}_{\mu _1,\mathrm{},\mu _j}^a`$ do not depend on various projections of indices due to the Lorentz invariance. But generally for the mixed projections the Lorentz spin $`j`$ of the tensor is less than the total number of its Lorentz indices $`J`$.

The matrix elements of the light-cone projections are expressed through solutions $`f_a`$ of the DGLAP equation (see Eqs.(4)). On the other hand the mixed projections

$`\stackrel{~}{n}^{\mu _1}\mathrm{}\stackrel{~}{n}^{\mu _{1+\omega }}l_{}^{\mu _{2+\omega }}\mathrm{}l_{}^{\mu _j}
{\displaystyle _0^1}𝑑xx^\omega {\displaystyle d^2k_{}\left(\frac{k_{}}{|k_{}|}\right)^n\phi _g(x,k_{}^2)},`$ (7)
can be found at small $`\omega `$ in terms of solutions $`\phi _g`$ of the BFKL equation.

Thus, it looks reasonable to extract some additional information concerning the parton $`x`$-distributions satisfying the DGLAP equation from the analogous $`k_{}`$-distributions satisfying the BFKL equation. Moreover, due to the fact, that in an extended $`N=4`$ SUSY the $`\beta `$-function vanishes, the $`4`$-dimensional conformal invariance may allow us to relate the Regge and Bjorken asymptotics of scattering amplitudes.

Our paper is organized as follows 222Some results were presented on the XXXV Winter School (see ). In Section 2 we discuss the relation between DGLAP and BFKL equations in the $`N=4`$ model obtained by an analytic continuation of the eigenvalue $`\omega ^0(n,\nu )`$ in $`|n|`$. In Sections 3 and 4 we review shortly our previous results for the BFKL kernel , rewrite them in a framework of the DRED scheme and investigate their properties. Section 5 is devoted to an independent calculation of anomalous dimensions in the leading and next-to-leading approximations for N=4 SUSY with the use of the renormalization group. In Section 6 we discuss the relation between the DGLAP and BFKL equations in the NLO approximation for this model. Appendices A, B and C are devoted to the derivation of some results used below. A summary is given in Conclusion.

## 2 Anomalous dimensions of twist-2 operators and their singularities

In the leading logarithmic approximation fermions and scalars do not give any contribution to the BFKL equation and therefore its integral kernel is the same at all supersymmetric gauge theories. Due to the Möbius invariance in the impact parameter space $`\stackrel{}{\rho }`$ the solution of the homogeneous BFKL equation has the form (see )

$$E_{\nu ,n}(\stackrel{}{\rho _{10}},\stackrel{}{\rho _{20}})<\varphi (\stackrel{}{\rho _1})O_{m,\stackrel{~}{m}}(\stackrel{}{\rho _0})\varphi (\stackrel{}{\rho _2})>=\left(\frac{\rho _{12}}{\rho _{10}\rho _{20}}\right)^m\left(\frac{\rho _{12}^{}}{\rho _{10}^{}\rho _{20}^{}}\right)^{\stackrel{~}{m}},$$

where

$$m=1/2+i\nu +n/2,\stackrel{~}{m}=1/2+i\nu n/2$$

are the conformal weights simply related to eigenvalues of the Casimir operators of the Möbius group. We use in $`E_{\nu ,n}(\stackrel{}{\rho _{10}},\stackrel{}{\rho _{20}})`$ the complex variables $`\rho _k=x_k+iy_k`$ in the transverse subspace and the notation $`\rho _{kl}=\rho _k\rho _l`$.

For a principal series of the unitary representations the quantities $`\nu `$ and $`n`$ are respectively real and integer numbers. The projection $`n`$ of the conformal spin $`|n|`$ can be positive or negative, but the eigenvalue of the BFKL equation in LLA

$`\omega =\omega ^0(n,\nu )={\displaystyle \frac{g^2N_c}{2\pi ^2}}\left(\mathrm{\Psi }(1)Re\mathrm{\Psi }\left({\displaystyle \frac{1}{2}}+i\nu +{\displaystyle \frac{|n|}{2}}\right)\right)`$ (8)
depends only on $`|n|`$. The Möbius invariance takes place also for the Schrödinger equation describing composite states of several reggeized gluons .

It is important, that the above expression has the property of the holomorphic separability

$$\omega ^0(n,\nu )=\omega ^0(m)+\omega ^0(\stackrel{~}{m}),\omega ^0(m)=\frac{g^2N_c}{8\pi ^2}\left(2\mathrm{\Psi }(1)\mathrm{\Psi }(m)\mathrm{\Psi }(1m)\right)$$
(9)
due to the relation

$$\mathrm{\Psi }(m)+\mathrm{\Psi }(1\stackrel{~}{m})=\mathrm{\Psi }(\stackrel{~}{m})+\mathrm{\Psi }(1m).$$

Moreover, the more general property $`H=h+h^{},[h,h^{}]=0`$ is valid also for the Hamiltonian of an arbitrary number of reggeized gluons in the multi-colour QCD $`N_c\mathrm{}`$ . For the case of the Odderon constructed from three gluons the holomorphic hamiltonian $`h`$ is integrable . Furthermore, in the multi-color QCD the BFKL dynamics for an arbitrary number of reggeized gluons is also completely integrable . It is related to the fact, that in this case the holomorphic Hamiltonian $`h`$ coincides with the local Hamiltonian for an integrable Heisenberg spin model . The theory turns out to be invariant under the duality transformation . Presumably the remarkable mathematical properties of the reggeon dynamics in LLA are consequences of the extended $`N=4`$ supersymmetry , which is one of the reasons for our investigation of the BFKL and DGLAP equations in this model. It was argued in Ref.  that generalized DGLAP equations for matrix elements of quasi-partonic operators  for $`N=4`$ SUSY are also integrable at large $`N_c`$. For the case of the Odderon a non-trivial integral of motion discovered in Ref.  was used to construct the energy spectrum and wave functions of the three-gluon state . A new Odderon solution with the intercept $`j_0`$ equal to 1 was found in . Recently the energy spectrum for reggeon composite states in the multi-color QCD was obtained through the solution of the Baxter equation (see refs.  and ). The effective action for reggeon interactions is known . It gives a possibility to calculate next-to-leading corrections for the gluon trajectory and reggeon couplings -.

The solution of the inhomogeneous BFKL equation in the LLA approximation can be written as the four-point Green function of a two-dimensional field theory

$$<\varphi (\stackrel{}{\rho _1})\varphi (\stackrel{}{\rho _2})\varphi (\stackrel{}{\rho _1^{}})\varphi (\stackrel{}{\rho _2^{}})>=\underset{n}{}_{\mathrm{}}^{\mathrm{}}𝑑\nu C(\nu ,|n|)d^2\rho _0\frac{E_{\nu ,|n|}(\stackrel{}{\rho _{10}},\stackrel{}{\rho _{20}})E_{\nu ,|n|}^{}(\stackrel{}{\rho _{1^{}0}},\stackrel{}{\rho _{2^{}0}})}{\omega \omega ^0(|n|,\nu )},$$

where $`C(\nu ,|n|)`$ is expressed through an inhomogeneous term of the equation with the use of a completeness condition for $`E_{\nu ,|n|}`$ (see ).

For $`\stackrel{}{\rho _1^{}}\stackrel{}{\rho _2^{}}`$ we obtain $`\stackrel{}{\rho _{01^{}}}\stackrel{}{\rho _{1^{}2^{}}}`$ and therefore

$`<\varphi (\stackrel{}{\rho _1})\varphi (\stackrel{}{\rho _2})\varphi (\stackrel{}{\rho _1^{}})\varphi (\stackrel{}{\rho _2^{}})>`$ $``$ $`{\displaystyle \underset{n}{}}{\displaystyle _{\mathrm{}}^{\mathrm{}}}𝑑\nu C(\nu ,|n|){\displaystyle \frac{E_{\nu ,|n|}(\stackrel{}{\rho _{11^{}}},\stackrel{}{\rho _{21^{}}})}{\omega \omega ^0(|n|,\nu )}}\rho _{1^{}2^{}}^m\rho _{1^{}2^{}}^{\stackrel{~}{m}}`$
$``$ $`{\displaystyle \underset{n}{}}C(\nu _\omega ,|n|){\displaystyle \frac{E_{\nu _\omega ,|n|}(\stackrel{}{\rho _{11^{}}},\stackrel{}{\rho _{21^{}}})}{\omega ^{}(|n|,\nu _\omega )}}|\rho _{1^{}2^{}}|^{1+2i\nu _\omega }\left({\displaystyle \frac{\rho _{1^{}2^{}}}{\rho _{1^{}2^{}}^{}}}\right)^{|n|/2},`$
where $`\nu _\omega `$ is a solution of the algebraic equation

$$\omega =\omega ^0(|n|,\nu )$$

with $`Im\nu _\omega <0`$.

The above asymptotics has a simple interpretation in terms of the Wilson operator-product expansion

$$\underset{\rho _1^{}\rho _2^{}}{lim}\varphi (\stackrel{}{\rho _1^{}})\varphi (\stackrel{}{\rho _2^{}})=\underset{n}{}\frac{C(\nu _\omega ,|n|)}{\omega ^{}(|n|,\nu _\omega )}|\rho _{1^{}2^{}}|^{2\mathrm{\Gamma }_\omega }\left(\frac{\rho _{1^{}2^{}}}{\rho _{1^{}2^{}}^{}}\right)^{|n|/2}O_{\nu _\omega ,|n|}(\stackrel{}{\rho _1^{}}),$$

where

$$\mathrm{\Gamma }_\omega =\frac{1}{2}+i\nu _\omega $$

is the transverse dimension of the operator $`O_{\nu _\omega ,|n|}(\stackrel{}{\rho _1^{}})`$ calculated in units of the squared mass. This operator is a mixed projection

$$O_{\nu _\omega ,|n|}(\stackrel{}{\rho _1^{}})=\stackrel{~}{n}^{\mu _1}\mathrm{}\stackrel{~}{n}^{\mu _{1+\omega }}l_{}^{\sigma _1}\mathrm{}l_{}^{\sigma _{|n|}}O_{\mu _1,\mathrm{},\mu _{1+\omega },\sigma _1,\mathrm{},\sigma _{|n|}}$$

of a gauge-invariant tensor $`O_{\mu _1,\mathrm{}\mu _J}`$ with $`J=1+\omega +|n|`$ indices. Note, that because $`\mathrm{\Gamma }_\omega `$ is real in the deep-inelastic regime $`\rho _{12}0`$ the operator $`O_{\nu _\omega ,|n|}(\stackrel{}{\rho _1^{}})`$ belongs to an exceptional series of unitary representations of the Möbius group (see ).

The anomalous dimension $`\gamma (j)`$ obtained from the BFKL equation in LLA (see (8)) has the poles

$$\mathrm{\Gamma }_\omega =1+\frac{|n|}{2}\gamma (j),\gamma (j)|_{\omega 0}=\frac{g^2N_c}{4\pi ^2\omega }.$$
(11)
Note, that there are also solutions of the BFKL equation corresponding to operators with the shifted anomalous dimensions

$$\mathrm{\Gamma }_\omega (k)=\mathrm{\Gamma }_\omega +k(k=1,2,3,\mathrm{}).$$

They contain the covariant laplacians $`D_\mu ^2`$ and have higher twists.

The canonical contribution $`1+|n|/2`$ to the transverse dimension $`\mathrm{\Gamma }_\omega `$ corresponds to the local operator

$$G_{\rho \mu _1}^{||}D_{\mu _2}^{||}\mathrm{}D_{\mu _\omega }^{||}D_{\sigma _1}^{}\mathrm{}D_{\sigma _{|n|}}^{}G_{\rho \mu _{1+\omega }}^{||}$$

in accordance with the fact, that in the light-cone gauge $`A_\mu \stackrel{~}{n}_\mu =0`$ the tensor

$$G_{\rho \mu _1}G_{\rho \mu _2}\stackrel{~}{n}^{\mu _1}\stackrel{~}{n}^{\mu _2}_{\mu _1}A_\rho ^{}_{\mu _2}A_\rho ^{}\stackrel{~}{n}^{\mu _1}\stackrel{~}{n}^{\mu _2}$$

has the transverse dimension equal to 1.

The local operator $`O_{\nu _\omega ,|n|}`$ for $`|n|=1,\mathrm{\hspace{0.17em}2},\mathrm{}`$ has the twist higher than 2 because its anomalous dimension $`\gamma `$ is singular at $`\omega 0`$. Indeed, for integer $`|n|>0`$ such singularity is impossible for twist-2 operators, because in this case for $`\omega 0`$ the total number of their indices $`1+\omega +|n|`$ would coincide with the physical value of the Lorentz spin $`j=1+|n|2`$. Moreover, the twist of $`O_{\nu _\omega ,|n|}`$ should grow at large $`|n|`$. For the case $`\omega >|n|>0`$ this operator is a mixed projection of the tensor

$$G_{\rho \mu _1}G_{\mu _2\sigma _1}\mathrm{}G_{\mu _{|n|+1}\sigma _{|n|}}D_{\mu _{|n|+2}}\mathrm{}D_{\mu _\omega }G_{\rho \mu _{1+\omega }}.$$

In particular, for $`n=0`$ its twist equals 2. For $`\omega <|n|`$ the corresponding operator is

$$G_{\rho \mu _1}G_{\mu _2\sigma _1}\mathrm{}G_{\mu _\omega \sigma _{\omega 1}}D_{\sigma _\omega }D_{\sigma _{\omega +1}}\mathrm{}D_{\sigma _n}G_{\rho \mu _{1+\omega }}$$

and its anomalous dimension has a singularity at $`\omega 0`$.

In the previous paper  we suggested to use an analytic continuation of the BFKL anomalous dimension $`\gamma (|n|,\omega )`$ to the points $`|n|=r1`$ ($`r=0,1,2,\mathrm{}`$) to calculate the singularities of the anomalous dimension of the twist-2 operators at negative integer $`j`$

$$j=1+\omega +|n|r.$$
(12)
Because for fixed $`|n|r1`$ and $`\omega 0`$ the quantity $`\gamma (|n|,\omega )`$ describes higher twist operators, to obtain the anomalous dimension $`\gamma `$ for the twist-2 operators one should push $`|n|`$ to $`r1`$ sufficiently rapidly at $`\omega 0`$

$$\mathrm{\Delta }(|n|)|n|+r+1=\stackrel{~}{C}_1(r)\omega ^2+\stackrel{~}{C}_2(r)\omega ^3+\mathrm{}.$$
(13)
Note, that for $`N=4`$ SUSY there are several multiplicatively renormalizable twist-2 operators (3) constructed from bosonic and fermionic fields. Their anomalous dimensions can be obtained from the universal function $`\gamma ^{uni}(j)`$ (see Eq. (71)) by a shift of its argument $`jj+k`$ with $`k=0,\mathrm{\hspace{0.17em}1},\mathrm{\hspace{0.17em}2},\mathrm{\hspace{0.17em}3},\mathrm{\hspace{0.17em}4}`$. It is related to the fact, that these operators enter in the same supermultiplet. We assume in accordance with Ref. , that from the BFKL equation one is able to calculate the singularities of $`\gamma ^{uni}(j)`$.

In LLA one can easily derive from the BFKL eigenvalue $`\omega ^0(n,\nu )`$ in the limit (12) the following result for $`\gamma (j)`$

$$\gamma (j)|_{jr}=\frac{g^2N_c}{4\pi ^2}\frac{1}{j+r}$$
(14)
for all $`r=1,0,1,\mathrm{}`$. Further, for the BFKL equation in LLA the fermions are not important, but generally they give non-vanishing contributions to the residues of the poles for $`\gamma (j)`$ in the DGLAP equation even in LLA. Therefore the result (14) for the anomalous dimension $`\gamma ^{uni}(j)`$ can be derived only for a certain generalization of QCD.

In our previous paper  it was shown, that only in an extended N=4 supersymmetric Yang-Mills theory the anomalous dimension, calculated in the next-to-leading approximation from the BFKL equation, can be analytically continued to negative $`|n|`$. Therefore the above result for $`\gamma (j)`$ in LLA should be valid only for $`N=4`$ SUSY. Indeed, using the conservation of the energy-momentum stress tensor $`T_{\mu \nu }`$ (corresponding to the condition $`\gamma (2)=0`$) to fix a subtraction constant in the expansion of $`\gamma (j)`$ over the poles at $`j=r`$, we obtain :

$$\gamma (j)=\frac{g^2N_c}{16\pi ^2}\gamma ^{LLA}(j),\gamma ^{LLA}(j)=4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j1)\right)$$

in an agreement with the direct calculation of $`\gamma ^{uni}(j)`$ in this theory (see  and Section 5 below). It is important that this anomalous dimension is the same for all twist-2 operators entering in the $`N=4`$ supermultiplet up to a shift of its argument by an integer number, because this property leads to an integrability of the evolution equation for matrix elements of quasi-partonic operators  in the multi-colour limit $`N_c\mathrm{}`$ (see ).

In the NLO approximation there is a more complicated situation. Namely, the anomalous dimension $`\gamma ^{uni}(j)`$ has the multiple poles $`\mathrm{\Delta }\gamma \alpha ^2(j+r)^3`$ at even $`r`$, which is related to an appearance of the double-logarithmic corrections $`(\alpha \mathrm{ln}^2s)^n`$ in the Regge limit $`s\mathrm{}`$ at upper orders $`n`$ of the perturbation theory.

The origin of the double-logarithmic terms can be understood in a simple way using as an example the process of the forward annihilation of the $`e^+e^{}`$ pair in the $`\mu ^+\mu ^{}`$ pair in QED . Let us write the scattering amplitude for this process in the following form

$$f_{e^+e^{}\mu ^+\mu ^{}}=4\pi \alpha _{em}\frac{\gamma _\sigma ^{}\gamma _\sigma ^{}}{s}A(s,k_{}^2),\alpha _{em}=\frac{e_{em}^2}{4\pi },s=2p_ek,$$

where $`p_e`$ and $`k`$ are momenta of the incoming electron and positron respectively and $`k_{}`$ is the transverse part of $`k`$. $`A(s,k_{})`$ is the Bethe-Salpeter amplitude for the electron interaction with the positron having the large virtuality $`k_{}^2m^2`$. The integral equation for this amplitude in the double-logarithmic approximation ($`\alpha _{em}\mathrm{ln}^2s1,\alpha _{em}1`$) can be presented as follows

$`A(s,k_{}^2)=1+{\displaystyle \frac{\alpha _{em}}{2\pi }}{\displaystyle _{m^2}^s}{\displaystyle \frac{dk_{}^{}{}_{}{}^{2}}{k_{}^{}{}_{}{}^{2}}}{\displaystyle _{k_{}^{}{}_{}{}^{2}}^{\stackrel{~}{s}}}{\displaystyle \frac{ds^{}}{s^{}}}A(s^{},k_{}^{}{}_{}{}^{2}),\left(\stackrel{~}{s}=\mathrm{min}\{s,sk_{}^{}{}_{}{}^{2}/k_{}^2\}\right),`$ (15)
where we took into account, that the double-logarithmic contribution appears from the integration region

$$s^{}s,\frac{k_{}^{}{}_{}{}^{2}}{s^{}}\frac{k_{}^2}{s},s^{}k_{}^2m^2.$$

One can search its solution in a form of the double Mellin transformation

$$A(s,k_{}^2)=_{\sigma i\mathrm{}}^{\sigma +i\mathrm{}}\frac{d\omega }{2\pi i}\left(\frac{s}{k_{}^2}\right)^\omega _i\mathrm{}^i\mathrm{}\frac{d\gamma }{2\pi i}\left(\frac{k_{}^2}{m^2}\right)^\gamma f_\omega (\gamma ).$$

Putting this anzatz in the equation we obtain a simple expression for $`f_\omega (\gamma )`$ 333 The expression for $`f_\omega (\gamma )`$ in the published version of this paper (see Nucl.Phys. B661 (2003) 19) was not right because it did not satisfy all boundary conditions contained in the integral equation (15). The correct expression (16) is derived in the Appendix D of this hep-ph version and is sent by us as an Erratum to Nucl.Phys.B. We thank J. Bartels and M. Lyublinsky who drew our attention to this problem.

$`f_\omega (\gamma )={\displaystyle \frac{(\omega 2\gamma )\gamma \omega \lambda ^2}{\omega \frac{\alpha _{em}}{2\pi }\left[\frac{1}{\omega \gamma }+\frac{1}{\gamma }\right]}}.`$ (16)
The pole of $`f_\omega (\gamma )`$ in $`\gamma `$ is situated at

$`\gamma ={\displaystyle \frac{\omega }{2}}\left(1\sqrt{1{\displaystyle \frac{2\alpha _{em}}{\pi \omega ^2}}}\right).`$ (17)
Expanding this expression for the anomalous dimension $`\gamma `$ in the series over $`\alpha _{em}`$ we shall generate the triple pole $`\mathrm{\Delta }\gamma \alpha _{em}^2/\omega ^3`$.

It is important to note, that the pole contributions

$$\frac{1}{\omega \gamma }+\frac{1}{\gamma }$$

in the denominator of $`f_\omega (\gamma )`$ are similar to singularities of the expression

$`\omega 2\mathrm{\Psi }(1)\mathrm{\Psi }(1+|n|+\omega \gamma )\mathrm{\Psi }(\gamma )`$ (18)
(at $`|n|=1,3,\mathrm{}`$) appearing for $`\omega 0`$ in the eigenvalue of the BFKL kernel of $`N=4`$ SUSY in a modified Born approximation (see Eqs. (44) and (45) below). For odd negative $`j`$ the double-logarithmic contributions are absent because the kernel of the Bethe-Salpeter equation has an odd number of factors $`k_{}`$ leading to a cancellation of one logarithm after the integration of $`d^2k_{}/k_{}^2`$ over the azimuthal angle $`\phi `$.

It is obvious from the above QED example, that the Bethe-Salpeter equation within the double-logarithmic accuracy can be derived both from the BFKL-like and DGLAP-like equations by imposing more accurate constraints on the region of integration. In the first case one should restrict $`k_{}^2`$ by $`s^{}`$ from above, which results in the substitution $`\gamma \omega \gamma `$ in the argument of the first $`\mathrm{\Psi }`$-function (see the above expression for $`\omega _0`$). In the second case one should restrict $`s^{}`$ by $`\stackrel{~}{s}`$ from above, which results in the substitution $`\gamma ^1\gamma ^1+(\omega \gamma )^1`$ and leads to a presence of the triple poles $`\mathrm{\Delta }\gamma \alpha ^2/\omega ^3`$ in the two-loop anomalous dimension. This procedure is similar to the rapidity veto approach invented by Bo Andersson and other authors to explain the appearance of the large NLO corrections to the BFKL kernel .

## 3 NLO corrections to the BFKL kernel in $`N=4`$ SUSY

To begin with, we review shortly the results of Ref. , where NLO corrections to the BFKL integral kernel at $`t=0`$ were calculated in the case of QCD and supersymmetric gauge theories. Only $`N=4`$ SUSY is discussed in details.

### 3.1 The set of eigenvalues

The eigenvalues of the homogeneous BFKL equation for the $`N=4`$ supersymmetric gauge theory can be written as follows

$`\omega =4\overline{a}\left[\chi (n,\gamma )+\left({\displaystyle \frac{1}{3}}\chi (n,\gamma )+\delta (n,\gamma )\right)\overline{a}\right],`$ (19)
where $`\gamma =\frac{1}{2}i\nu `$ and

$`\chi (n,\gamma )`$ $`=`$ $`2\mathrm{\Psi }(1)\mathrm{\Psi }(M)\mathrm{\Psi }(1\stackrel{~}{M}),`$ (20)
$`\delta (n,\gamma )`$ $`=`$ $`\mathrm{\Psi }^{\prime \prime }(M)+\mathrm{\Psi }^{\prime \prime }(1\stackrel{~}{M})+6\zeta (3)2\zeta (2)\chi (n,\gamma )`$ (21)
$``$ $`2\mathrm{\Phi }(|n|,\gamma )2\mathrm{\Phi }(|n|,1\gamma )`$
with

$$M=\gamma +\frac{|n|}{2},\stackrel{~}{M}=\gamma \frac{|n|}{2}$$

Here $`\mathrm{\Psi }(z)`$, $`\mathrm{\Psi }^{}(z)`$ and $`\mathrm{\Psi }^{\prime \prime }(z)`$ are the Euler $`\mathrm{\Psi }`$ -function and its derivatives, respectively; $`\overline{a}=g^2N_c/(16\pi ^2)`$ is expressed in terms of the coupling constant $`g`$ in the DREG scheme. The function $`\mathrm{\Phi }(|n|,\gamma )`$ is given below 444Notice that the representation of $`\mathrm{\Phi }(|n|,\gamma )`$ in terms of the sum over $`k`$ in ref.  contained a misprint: the factors $`(1)^{k+1}`$ in the last sum of (22) were substituted by $`(1)^k`$.

$`\mathrm{\Phi }(|n|,\gamma )`$ $`=`$ $`{\displaystyle _0^1}dx{\displaystyle \frac{x^{M1}}{1+x}}\left[{\displaystyle \frac{1}{2}}\right(\mathrm{\Psi }^{}\left({\displaystyle \frac{|n|+1}{2}}\right)\zeta (2))+\mathrm{Li}_2(x)+\mathrm{Li}_2(x)`$ (22)
$`+`$ $`\mathrm{ln}(x)\left(\mathrm{\Psi }(|n|+1)\mathrm{\Psi }(1)+\mathrm{ln}(1+x)+{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(x)^k}{k+|n|}}\right)`$
$`+`$ $`{\displaystyle \underset{k=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{x^k}{(k+|n|)^2}}(1(1)^k)]`$
$`={\displaystyle \underset{k=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(1)^{k+1}}{k+M}}[\mathrm{\Psi }^{}(k+|n|+1)\mathrm{\Psi }^{}(k+1)+(1)^{k+1}(\beta ^{}(k+|n|+1)+\beta ^{}(k+1)))`$
$``$ $`{\displaystyle \frac{1}{k+M}}(\mathrm{\Psi }(k+|n|+1)\mathrm{\Psi }(k+1))],`$
and

$`\beta ^{}(z)={\displaystyle \frac{1}{4}}\left[\mathrm{\Psi }^{}\left({\displaystyle \frac{z+1}{2}}\right)\mathrm{\Psi }^{}\left({\displaystyle \frac{z}{2}}\right)\right]={\displaystyle \underset{r=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(1)^r}{(z+r)^2}}.`$ (23)
Note, that the term

$$\frac{1}{3}\chi (n,\gamma )$$

appears as a result of the use of the DREG scheme (see eq. (29) in ). It is well known, that DREG scheme violates SUSY. The DREG procedure could be used provided that the numbers of boson and fermion degrees of freedom would coincide (see ). But in its standard version the number of degrees of freedom for vector fields is $`4+2\widehat{\epsilon }`$ 555Here $`\widehat{\epsilon }=(D4)/2`$ and $`D`$ is the space-time dimension. and does not depend on $`\widehat{\epsilon }`$ for other particles. Usually (see discussions in -) it is possible to restore the equality between boson and fermion degrees of freedom in SUSY by adding $`2\widehat{\epsilon }`$ additional scalar fields. Formally it is equivalent to the dimensional reduction (for $`N=4`$ SUSY) from the space-time dimension $`10`$ to the dimension $`4+2\widehat{\epsilon }`$ and we call it below the DRED scheme. The change of the number of scalars from $`6`$ to $`62\widehat{\epsilon }`$ does not violate the cancellation of the singular contribution $`\widehat{\epsilon }^1`$ to the coupling constant (the $`\beta `$-function is zero for $`N=4`$ SUSY in the DRED scheme 666In the DRED scheme the $`\beta `$-function vanishing was demonstrated at first three orders of perturbation theory (see  and references therein).) but leads to an additional contribution to $`\omega `$ of the order of $`\overline{a}^2`$. Indeed, the correct result in the framework of the DRED scheme can be obtained from eq.(19) by a finite renormalization of the coupling constant (see also G. Altarelli et al. and G.A. Shuler et al. in )

$`\overline{a}\widehat{a}=\overline{a}+{\displaystyle \frac{1}{3}}\overline{a}^2,`$ (24)
which eliminates the term $`\widehat{a}^2\chi `$ in (19). For the new coupling constant $`\widehat{a}`$ in the DRED scheme the above expression for $`\omega `$ can be written in the following simple form

$$\omega =4\widehat{a}\left[\chi (n,\gamma )+\delta (n,\gamma )\widehat{a}\right].$$

### 3.2 Hermitian separability of the Bethe-Salpeter kernel

Following the method of Ref.  and using the results of the previous section we can write the above NLO correction to the BFKL equation in the form having the property similar to the holomorphic separability (cf. ). This property can be called the hermitian separability, because it guarantees the reality of the eigenvalue $`\omega `$ for real $`\nu `$.

Indeed, using the results of Appendix A we can obtain for the most complicated contribution to $`\delta (n,\gamma )`$

$`\mathrm{\Phi }(|n|,\gamma )+\mathrm{\Phi }(|n|,1\gamma )=\chi (n,\gamma )\left(\beta ^{}(M)+\beta ^{}(1\stackrel{~}{M})\right)`$
$`+\mathrm{\Phi }_2(M)\beta ^{}(M)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(M)\right]+\mathrm{\Phi }_2(1\stackrel{~}{M})\beta ^{}(1\stackrel{~}{M})\left[\mathrm{\Psi }(1)\mathrm{\Psi }(1\stackrel{~}{M})\right],`$
where $`\chi (n,\gamma )`$ is given by Eq.(20).

Thus, we can rewrite the NLO corrections $`\delta (n,\gamma )`$ in a (generalized) hermitian separable form (providing that $`\omega _0`$ is substituted by $`\omega `$, which is valid with our accuracy):

$`\delta (n,\gamma )`$ $`=`$ $`\varphi (M)+\varphi (1\stackrel{~}{M}){\displaystyle \frac{\omega _0}{2\widehat{a}}}\left(\rho (M)+\rho (1\stackrel{~}{M})\right),`$ (25)
$`\omega _0`$ $`=`$ $`4\widehat{a}\left(2\mathrm{\Psi }(1)\mathrm{\Psi }(M)\mathrm{\Psi }(1\stackrel{~}{M})\right)`$ (26)
and

$`\rho (M)`$ $`=`$ $`\beta ^{}(M)+{\displaystyle \frac{1}{2}}\zeta (2),`$ (27)
$`\varphi (M)`$ $`=`$ $`3\zeta (3)+\mathrm{\Psi }^{^{\prime \prime }}(M)2\mathrm{\Phi }_2(M)+2\beta ^{^{}}(M)\left(\mathrm{\Psi }(1)\mathrm{\Psi }(M)\right).`$ (28)
In a simpler way the hermitian separability takes place in the eigenvalue equation for the corresponding Bethe-Salpeter kernel

$`1`$ $`=`$ $`{\displaystyle \frac{4\widehat{a}}{\omega }}\left(2\mathrm{\Psi }(1)\mathrm{\Psi }(M)\mathrm{\Psi }(1\stackrel{~}{M})+\widehat{a}\left(\varphi (M)+\varphi (1\stackrel{~}{M})\right)\right)`$ (29)
$`2\widehat{a}\left(\rho (M)+\rho (1\stackrel{~}{M})\right).`$
In the right hand side the first contribution corresponds to its singularity at $`l=1`$ in the Gribov-Froissart representation generated by a pole of the Legendre function $`Q_l(z)`$ at $`\omega =1+l0`$ and the last term appears from the regular part of the Born contribution. Note, that $`M`$ and $`1\stackrel{~}{M}`$ coincide with the anomalous dimensions appearing in the asymptotic expressions for the BFKL kernel in the limits when the gluon virtualities are large: $`q_1^2\mathrm{}`$ and $`q_2^2\mathrm{}`$, respectively. Because $`1\stackrel{~}{M}=M^{}`$, the hermitian separability guarantees the symmetry of $`\omega `$ for the principal series of the unitary representations of the Möbius group (with $`m=1/2+i\nu +n/2`$) to the substitution $`\nu \nu `$ and the hermicity of the BFKL hamiltonian.

### 3.3 The violation of a generalized holomorphical separability

Eqs.(25)-(28) show the hermitian separability with two contributions depending on $`M=\gamma +|n|/2`$ and $`M^{}=1\stackrel{~}{M}=1\gamma +|n|/2`$. Although the above expressions are symmetric to $`MM^{}`$, it is not clear, if we can symmetrize $`\delta (n,\gamma )`$ to the substitution $`|n||n|`$ as it was done in the Born approximation to obtain a true holomorphic separability (see Eq.(9)). Using properties of polygamma-functions and Eqs. (A4) and (A5) one can attempt to present Eqs.(25)-(28) in a separable form symmetric to the substitution $`m\stackrel{~}{m}`$ (i.e., respectively, $`M\stackrel{~}{M}`$ for both positive and negative $`n`$). From the results of Appendix B the NLO correction $`\delta (n,\gamma )`$ can be written in the form

$`\delta (n,\gamma )`$ $`=`$ $`\overline{\varphi }(m)+\overline{\varphi }(\stackrel{~}{m}){\displaystyle \frac{\omega _0}{2\widehat{a}}}\left(\overline{\rho }(m)+\overline{\rho }(\stackrel{~}{m})\right)`$ (30)
$``$ $`\pi ^3\left({\displaystyle \frac{\mathrm{cos}^2(m\pi )}{\mathrm{sin}^3(m\pi )}}+{\displaystyle \frac{\mathrm{cos}^2(\stackrel{~}{m}\pi )}{\mathrm{sin}^3(\stackrel{~}{m}\pi )}}\right)+\delta _2^{(2)}(n,\gamma ),`$
where

$`\delta _2^{(2)}(n,\gamma )`$ $`=`$ $`(1+(1)^n){\displaystyle \frac{\pi ^2\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )}}\left[\mathrm{\Psi }(M)\mathrm{\Psi }(\stackrel{~}{M})\right]`$ (31)
$`=`$ $`(1+(1)^n)\epsilon (n)\pi ^2\left({\displaystyle \frac{\mathrm{\Psi }(m)\mathrm{cos}(m\pi )}{\mathrm{sin}^2(m\pi )}}{\displaystyle \frac{\mathrm{\Psi }(\stackrel{~}{m})\mathrm{cos}(\stackrel{~}{m}\pi )}{\mathrm{sin}^2(\stackrel{~}{m}\pi )}}\right)`$
with

$$\epsilon (n)=\{\begin{array}{cc}+1,& n0\\ 1,& n<0\end{array}.$$

The term $`\delta _2^{(2)}(n,\gamma )`$ violates the holomorphic separability. Because the contribution proportional to $`\omega _0`$ in (30) depends not only on $`M`$ and $`1\stackrel{~}{M}`$ we have the situation similar to the case of quantum anomalies. Namely, the property of the hermitian separability of $`\delta (n,\gamma )`$ is responsible for the violation of the holomorphic separability in $`\delta _2^{(2)}(n,\gamma )`$. The anomalous term $`\delta _2^{(2)}(n,\gamma )`$ is zero for odd $`n`$, where $`\delta (n,\gamma )`$ coincides with its analytic continuation to corresponding negative $`|n|`$. Note, that for the Pomeron only even values of $`n`$ are physical due to the Bose symmetry of its two-gluon wave function. The colorless composite state of three reggeized gluons with the $`f`$-coupling in the color space and odd $`n`$ has an anti-symmetric wave function and can give a large contribution to the small-$`x`$ behaviour of the structure function $`g_2(x)`$ (see ).

### 3.4 Asymptotics of “cross-sections” at $`s\mathrm{}`$

As an application of obtained results we consider the cross-section for the inclusive production of two pairs of particles with their mass $`\mu `$ in the polarized $`\gamma \gamma `$ collisions (see ):

$$\sigma (s)=\alpha _{em}^2\widehat{a}^2\frac{1}{\mu ^2}\frac{32}{81}\left(\sigma _0(s)+\left(\mathrm{cos}^2\vartheta \frac{1}{2}\right)\sigma _2(s)\right),$$

where $`\alpha _{em}`$ is the electromagnetic fine structure constant, the coefficient $`\sigma _0(s)`$ is proportional to the cross-section for the interaction of unpolarized photons and $`\sigma _2(s)`$ describes the spin correlation depending on the azimuthal angle $`\vartheta `$ between the polarization vectors of colliding photons in their c.m system.

The asymptotic behavior of the cross-sections $`\sigma _k(s)(k=0,2)`$ at $`s\mathrm{}`$ corresponds to an unmoving singularity of the $`t`$-channel partial wave $`f_\omega (t)(\omega \omega _k)^{1/2}`$ situated at

$`\omega _k`$ $`=`$ $`4\widehat{a}\left[\chi (k,{\displaystyle \frac{1}{2}})+\widehat{a}\delta (k,{\displaystyle \frac{1}{2}})\right]4\widehat{a}\chi (k,{\displaystyle \frac{1}{2}})\left[1\widehat{a}c(k,{\displaystyle \frac{1}{2}})\right]\left(c={\displaystyle \frac{\delta }{\chi }}\right),`$
$`\sigma _0(s)`$ $`=`$ $`{\displaystyle \frac{9\pi ^{5/2}}{32\sqrt{7\zeta (3)}}}{\displaystyle \frac{s^{\omega _0}}{\left(\mathrm{ln}(s/s_0)\right)^{1/2}}}\left(1+O(\widehat{a})\right),`$ (32)
$`\sigma _2(s)`$ $`=`$ $`{\displaystyle \frac{\pi ^{5/2}}{932\sqrt{7\zeta (3)8}}}{\displaystyle \frac{s^{\omega _2}}{\left(\mathrm{ln}(s/s_0)\right)^{1/2}}}\left(1+O(\widehat{a})\right).`$ (33)
Here the symbol $`O(\widehat{a})`$ denotes unknown next-to-leading corrections to the impact factors. Note, that strictly speaking for the $`N=4`$ supersymmetric gauge theory we can not consider electro-magnetic interactions without the SUSY violation and it could be more natural to investigate the interaction of gravitons or their superpartners. But here we want only to illustrate the relative magnitudes of the radiative corrections to the BFKL equation in QCD and in N=4 SUSY.

Using our results (25) and (26), we obtain in the $`N=4`$ case the following numerical values for $`\chi (k,1/2)`$ and $`c(k,1/2)`$ $`(k=0,2)`$:

$`\chi (0,{\displaystyle \frac{1}{2}})`$ $`=`$ $`4\mathrm{ln}2,\chi (2,{\displaystyle \frac{1}{2}})=4(\mathrm{ln}21),`$ (34)
$`c(0,{\displaystyle \frac{1}{2}})`$ $`=`$ $`2\zeta (2)+{\displaystyle \frac{1}{2\mathrm{ln}2}}\left[11\zeta (3)32\mathrm{L}\mathrm{s}_3\left({\displaystyle \frac{\pi }{2}}\right)14\pi \zeta (2)\right]=7.5812,`$ (35)
$`c(2,{\displaystyle \frac{1}{2}})`$ $`=`$ $`2\zeta (2)+{\displaystyle \frac{1}{2(\mathrm{ln}21)}}\left[11\zeta (3)+32\mathrm{L}\mathrm{s}_3\left({\displaystyle \frac{\pi }{2}}\right)+14\pi \zeta (2)32\mathrm{ln}2\right]`$ (36)
$`=`$ $`6.0348,`$
where (see )

$$\mathrm{Ls}_3(x)=_0^x\mathrm{ln}^2\left|2\mathrm{sin}\left(\frac{y}{2}\right)\right|dy.$$

Note, that the function $`\mathrm{Ls}_3(x)`$ appears also in contributions of some massive diagrams (see, for example, the recent papers  and references therein).

The LLA results (34) coincide with ones obtained in ref. . As it was shown in , in the framework of QCD the NLO correction $`c^{QCD}(0,1/2)`$ is large and leads to a strong reduction of the value of the Pomeron intercept (see recent analyses  of various resummations of the large NLO terms). Contrary to $`c^{QCD}(0,1/2)`$, the correction $`c(0,1/2)`$ is not large ($`c^{QCD}(0,1/2)/c(0,1/2)3.5`$). Because in N=4 SUSY the $`\beta `$-function is zero, it is natural to interpret the large correction to the intercept of the BFKL Pomeron in QCD as an effect related to the coupling constant renormalization. It seems to support the results of refs.  (see also the recent review  and references therein) concerning a large value for the argument of the running QCD coupling constant at high energy processes. The corrections $`c^{QCD}(2,1/2)`$ (see ) and $`c(2,1/2)`$ for $`N=4`$ are small and do not change significantly the LO value  of the angle-dependent contribution.

## 4 Non-symmetric choice of the energy normalization

Analogously to refs.  one can calculate eigenvalues of the BFKL kernel in the case of a non-symmetric choice for the energy normalization $`s_0`$ in eq.(15) related to the interpretation of the NLO corrections in the framework of the renormalization group (cf. ). For the scale $`s_0=q^2`$ natural for the deep-inelastic scattering we have respectively

in DREG-scheme

$$\omega =4\overline{a}[\chi (n,\gamma \frac{\omega }{2})+(\frac{1}{3}\chi (n,\gamma )+\delta (n,\gamma )\left)\overline{a}\right]$$

$$=4\overline{a}[\chi (n,\gamma )+(\frac{1}{3}\chi (n,\gamma )+\stackrel{~}{\delta }(n,\gamma )\left)\overline{a}\right],$$

in DRED-scheme

$$\omega =4\widehat{a}[\chi (n,\gamma )+\stackrel{~}{\delta }(n,\gamma ),\widehat{a}],$$

where

$`\stackrel{~}{\delta }(n,\gamma )=\delta (n,\gamma )2\chi (n,\gamma )\chi ^{}(n,\gamma )`$ (37)
and

$$\chi ^{}(n,\gamma )\frac{d}{d\gamma }\chi (n,\gamma )=\mathrm{\Psi }^{}\left(M\right)+\mathrm{\Psi }^{}\left(1\stackrel{~}{M}\right).$$

Note, that here $`\gamma `$ does not coincide with the anomalous dimension of the higher-twist operators with $`|n|1`$ (see below).

### 4.1 Limit $`\gamma 0`$ for $`n=0`$

By considering the limit $`\gamma 0`$ of the BFKL eigenvalue one can obtain for $`n=0`$ (see also the analysis in )

$`\chi (0,\gamma )`$ $`=`$ $`{\displaystyle \frac{1}{\gamma }}+O(\gamma ^2),`$
$`{\displaystyle \frac{1}{3}}\chi (n,\gamma )`$ $`+`$ $`\stackrel{~}{\delta }(0,\gamma )={\displaystyle \frac{B^{DREG}}{\gamma }}+C+O(\gamma ^2),`$
$`\stackrel{~}{\delta }(0,\gamma )`$ $`=`$ $`{\displaystyle \frac{B^{DRED}}{\gamma }}+C+O(\gamma ^2),`$ (38)
where

$`B^{DREG}={\displaystyle \frac{1}{3}},B^{DRED}=0,C=2\zeta (3).`$ (39)
According to refs. with the use of eqs.(38) and (39 ) one can calculate the anomalous dimensions $`\gamma `$ of the twist-2 operators at $`\omega 0`$ (i.e. near $`j=1`$) respectively

in DREG-scheme

$`\gamma =4\overline{a}[({\displaystyle \frac{1}{\omega }}+O(\omega ))+\overline{a}({\displaystyle \frac{B^{DREG}}{\omega }}+O(1))+\overline{a}^2({\displaystyle \frac{C}{\omega ^2}}+O\left(\omega ^1\right))],`$ (40)
in DRED-scheme

$`\gamma =4\widehat{a}[({\displaystyle \frac{1}{\omega }}+O(\omega ))+\widehat{a}({\displaystyle \frac{B^{DRED}}{\omega }}+O(1))+\widehat{a}^2({\displaystyle \frac{C}{\omega ^2}}+O\left(\omega ^1\right))].`$ (41)
Thus, in the framework of the DRED scheme the singular contribution $`\widehat{a}^2/\omega `$ to the anomalous dimension is zero. These results will be used below for the calculation of NLO corrections to $`\gamma `$ from the DGLAP equation (see Eqs. (68) and (71)).

### 4.2 Hermitian separability for a non-symmetric normalization and the symmetry $`\gamma J\gamma `$

For the scale $`s_0=q^2`$ the expression for $`\omega `$ as a function of the anomalous dimension $`\gamma =1/2+i\nu +|n|/2`$ (correctly defined at general $`n`$ (see (11)) can be written in the following form

$$\omega =4\widehat{a}\left(\chi (n,\gamma \frac{|n|}{2})+\widehat{a}\stackrel{~}{\delta }(n,\gamma \frac{|n|}{2})\right),$$

where $`\chi (n,\gamma )`$ and $`\stackrel{~}{\delta }(n,\gamma )`$ (for the non-symmetric choice of $`s_0`$) are given by Eqs. (20), (21) and (37), respectively.

By summing double-logarithmic terms in all orders of the perturbation theory (see discussion in Section 2) we can write $`\omega `$ in the ”Lorentz invariant” form

$`\omega =4\widehat{a}\left(2\mathrm{\Psi }(1)\mathrm{\Psi }(\gamma )\mathrm{\Psi }(1+|n|+\omega \gamma )+\stackrel{~}{\mathrm{\Delta }}(n,\gamma )\widehat{a}\right),`$ (42)
where the Lorentz spin $`j`$ of the corresponding operators (for negative $`|n|`$) is given by (12). For $`n=0`$ the Lorentz spin of the twist-2 operator is $`j=1+\omega `$.

In the above expression

$`\stackrel{~}{\mathrm{\Delta }}(n,\gamma )=\delta (n,\gamma {\displaystyle \frac{|n|}{2}})+2\left[\mathrm{\Psi }^{}(\gamma )+\mathrm{\Psi }^{}(|n|+1\gamma )\right]\chi (n,\gamma {\displaystyle \frac{|n|}{2}}).`$ (43)
Using the analysis of the subsection 3.2 one can present $`\omega `$ as follows

$`\omega =4\widehat{a}\left(2\mathrm{\Psi }(1)\mathrm{\Psi }(\gamma )\mathrm{\Psi }(|n|+1+\omega \gamma )+\epsilon \right),`$ (44)
where $`\epsilon `$ is written for $`\omega 0`$ in the form

$`\epsilon `$ $`=`$ $`{\displaystyle \frac{\omega }{2}}\left(p(\gamma )+p(1+|n|\gamma )\right)+\widehat{a}\left(\varphi (\gamma )+\varphi (1+|n|\gamma )\right),`$
$`\omega `$ $`=`$ $`4\widehat{a}\left(2\mathrm{\Psi }(1)\mathrm{\Psi }(\gamma )\mathrm{\Psi }(|n|+1\gamma )\right)+O(\widehat{a}^2).`$ (45)
Here

$$p(\gamma )=\mathrm{\Psi }^{}(\gamma )\rho (\gamma )=\mathrm{\Psi }^{}(\gamma )\beta ^{}(\gamma )\frac{1}{2}\zeta (2)=2\underset{k=0}{\overset{\mathrm{}}{}}\frac{1}{(\gamma +2k)^2}\frac{1}{2}\zeta (2)$$

and $`\varphi (\gamma )`$ is given by Eq.(28).

In accordance with the hermitian separability of the BFKL kernel established in the subsection 3.2 the eigenvalue equation for the corresponding Bethe-Salpeter equation can be written as follows

$`1`$ $`=`$ $`{\displaystyle \frac{4\widehat{a}}{\omega }}\left(2\mathrm{\Psi }(1)\mathrm{\Psi }(\gamma )\mathrm{\Psi }(J\gamma )+\widehat{a}\left(\varphi (\gamma )+\varphi (J\gamma )\right)\right)`$ (46)
$`+2\widehat{a}\left(p(\gamma )+p(J\gamma )\right),J=1+|n|+\omega ,`$
where $`J`$ is the total number of tensor indices of the local operator which does not coincide generally with its Lorentz spin $`j`$. Note that without the shift $`1+|n|\gamma J\gamma `$ of the argument of $`\mathrm{\Psi }`$-function at the LLA level, the symmetry of the eigenvalue to the substitution $`\gamma J\gamma `$ is violated by the term $`\chi ^{}(n,\gamma |n|/2)`$.

Let us calculate $`\stackrel{~}{\mathrm{\Delta }}(n,\gamma )`$ near its singularities at small $`\gamma `$. To begin with, we consider $`\gamma 0`$ for the physical integer values $`|n|0`$:

$$\stackrel{~}{\mathrm{\Delta }}(n,\gamma )\frac{4}{\gamma ^2}\left(\mathrm{\Psi }(1)\mathrm{\Psi }(|n|+1)\right)+\frac{2}{\gamma }\left(c(n)+2\mathrm{\Psi }^{}(|n|+1)\right),$$
(47)
where

$$c(n)=\mathrm{\Psi }^{}(|n|+1)\mathrm{\Psi }^{}(1)\beta ^{}(|n|+1)+\beta ^{}(1).$$
(48)
Therefore by solving the equation $`\omega =\omega (n,\gamma )`$ one can obtain

$$\gamma =\frac{4\widehat{a}}{\omega }\left(1+\omega \left(\mathrm{\Psi }(1)\mathrm{\Psi }(|n|+1)\right)\right)+\frac{(4\widehat{a})^2}{\omega ^2}\left(\mathrm{\Psi }(1)\mathrm{\Psi }(|n|+1)+\frac{\omega }{2}c(n)\right).$$
(49)
At $`n=0`$ the correction $`\widehat{a}`$ to $`\gamma `$ is absent, but for other $`n`$ we have the large correction

$$\mathrm{\Delta }\gamma =4\widehat{a}\left(\mathrm{\Psi }(1)\mathrm{\Psi }(|n|+1)\right),$$

having the poles at $`1+\omega +|n|r`$, which leads to a contribution changing even the singularities of the Born term. The explanation of this effect is related to the presence of the double-logarithmic terms $`\mathrm{\Delta }\gamma \widehat{a}^2/\omega ^3`$ near the points $`j=0,2,\mathrm{}`$ (see Sections 2 and 6). We remind, that for positive integer $`|n|`$ we calculate the anomalous dimensions $`\gamma `$ of the higher twist operators (with an anti-symmetrization between $`n`$ transversal and $`1+\omega `$ longitudinal indices). The singularities of the anomalous dimensions of the twist-2 operators can be obtained only in the limit when $`\mathrm{\Delta }(|n|)`$ tends to zero more rapidly than $`\omega `$ (see Eq. (13).

## 5 Anomalous dimension matrix in the $`N=4`$ SUSY

The DGLAP evolution equation for the moments of parton distributions for $`N=4`$ SUSY has the form

$`{\displaystyle \frac{d}{d\mathrm{ln}Q^2}}f_a(j,Q^2)`$ $`=`$ $`{\displaystyle \underset{k}{}}\gamma _{ab}(j)f_b(j,Q^2)(a,b=q,g,\phi ),`$ (50)
$`{\displaystyle \frac{d}{d\mathrm{ln}Q^2}}\mathrm{\Delta }f_a(j,Q^2)`$ $`=`$ $`{\displaystyle \underset{k}{}}\stackrel{~}{\gamma }_{ab}(j)\mathrm{\Delta }f_b(j,Q^2)(a,b=q,g),`$ (51)
where the anomalous dimension matrices $`\gamma _{ab}(j)`$ and $`\stackrel{~}{\gamma }_{ab}(j)`$ can be obtained in the form of expansions over the coupling constant $`\widehat{a}`$

$`\gamma _{ab}(j)`$ $`=`$ $`\widehat{a}\gamma _{ab}^{(0)}(j)+\widehat{a}^2\gamma _{ab}^{(1)}(j)+O(\widehat{a}^3),`$
$`\stackrel{~}{\gamma }_{ab}(j)`$ $`=`$ $`\widehat{a}\stackrel{~}{\gamma }_{ab}^{(0)}(j)+\widehat{a}^2\stackrel{~}{\gamma }_{ab}^{(1)}(j)+O(\widehat{a}^3).`$ (52)
In the following subsections we will present the results of exact calculations for the leading order coefficients $`\gamma _{ab}^{(0)}(j)`$ and $`\stackrel{~}{\gamma }_{ab}^{(0)}(j)`$ and construct the anomalous dimensions of the multiplicatively renormalizable twist-2 operators. In the NLO approximation the corresponding coefficients $`\gamma _{ab}^{(1)}(j)`$ and $`\stackrel{~}{\gamma }_{ab}^{(1)}(j)`$ were unknown (see, however, ). It is important, that the form of the LLA anomalous dimension matrix of the multiplicatively renormalizable operators in $`N=4`$ SUSY is rather simple because the result is expressed in terms of one function. Taking into account this universality related to the superconformal invariance, the existing information about the anomalous dimensions of twist-2 operators in the QCD case, the known NLO corrections to the BFKL kernel and an experience in integrating certain types of the Feynman diagrams (see, for example, ), we derive below the expressions for the NLO anomalous dimensions of the multiplicatively renormalizable operators in $`N=4`$ SUSY. This result is checked by direct calculations of the matrix elements $`\gamma _{ab}^{(1)}(j)`$ and $`\stackrel{~}{\gamma }_{ab}^{(1)}(j)`$ .

### 5.1 LLA results for the anomalous dimension matrix in $`N=4`$ SUSY

The elements of the LLA anomalous dimension matrix in the $`N=4`$ SUSY have the following form (see ):

for tensor twist-2 operators

$`\gamma _{gg}^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j1){\displaystyle \frac{2}{j}}+{\displaystyle \frac{1}{j+1}}{\displaystyle \frac{1}{j+2}}\right),`$
$`\gamma _{qg}^{(0)}(j)`$ $`=`$ $`8\left({\displaystyle \frac{1}{j}}{\displaystyle \frac{2}{j+1}}+{\displaystyle \frac{2}{j+2}}\right),\gamma _{\phi g}^{(0)}(j)=12\left({\displaystyle \frac{1}{j+1}}{\displaystyle \frac{1}{j+2}}\right),`$
$`\gamma _{gq}^{(0)}(j)`$ $`=`$ $`2\left({\displaystyle \frac{2}{j1}}{\displaystyle \frac{2}{j}}+{\displaystyle \frac{1}{j+1}}\right),\gamma _{q\phi }^{(0)}(j)={\displaystyle \frac{8}{j}},`$
$`\gamma _{qq}^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j)+{\displaystyle \frac{1}{j}}{\displaystyle \frac{2}{j+1}}\right),\gamma _{\phi q}^{(0)}(j)={\displaystyle \frac{6}{j+1}},`$
$`\gamma _{\phi \phi }^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j+1)\right),\gamma _{g\phi }^{(0)}(j)=4\left({\displaystyle \frac{1}{j1}}{\displaystyle \frac{1}{j}}\right),`$ (53)
for the pseudo-tensor operators:

$`\stackrel{~}{\gamma }_{gg}^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j+1){\displaystyle \frac{2}{j+1}}+{\displaystyle \frac{2}{j}}\right),`$
$`\stackrel{~}{\gamma }_{qg}^{a,(0)}(j)`$ $`=`$ $`8\left({\displaystyle \frac{1}{j}}+{\displaystyle \frac{2}{j+1}}\right),\stackrel{~}{\gamma }_{gq}^{(0)}(j)=2\left({\displaystyle \frac{2}{j}}{\displaystyle \frac{1}{j+1}}\right),`$
$`\stackrel{~}{\gamma }_{qq}^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j+1)+{\displaystyle \frac{1}{j+1}}{\displaystyle \frac{1}{j}}\right).`$ (54)
Note, that in the $`N=4`$ SUSY multiplet there are also twist-2 operators with fermion quantum numbers but their anomalous dimensions coincide up to an integer shift of the argument with the above expressions for the bosonic components (cf. ref. ).

### 5.2 Anomalous dimensions and twist-2 operators with a multiplicative renormalization

It is possible to construct five independent twist-two operators with a multiplicative renormalization. The corresponding parton distributions and their LLA anomalous dimensions have the form (see ):

$`f_I(j)`$ $`=`$ $`f_g(j)+f_q(j)+f_\phi (j)f_{q,g,\phi }^+(j),`$
$`\gamma _I^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j1)\right)4S_1(j2)\gamma _+^{(0)}(j),`$ (55)
$`f_{II}(j)`$ $`=`$ $`2(j1)f_g(j)+f_q(j)+{\displaystyle \frac{2}{3}}(j+1)f_\phi (j)f_{q,g,\phi }^0(j),`$
$`\gamma _{II}^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j+1)\right)4S_1(j)\gamma _0^{(0)}(j),`$ (56)
$`f_{III}(j)`$ $`=`$ $`{\displaystyle \frac{j1}{j+2}}f_g(j)+f_q(j){\displaystyle \frac{j+1}{j}}f_\phi (j)f_{q,g,\phi }^{}(j),`$
$`\gamma _{III}^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j+3)\right)4S_1(j+2)\gamma _{}^{(0)}(j),`$ (57)
$`f_{IV}(j)`$ $`=`$ $`2\mathrm{\Delta }f_g(j)+\mathrm{\Delta }f_q(j)\mathrm{\Delta }f_{q,g}^+(j),`$
$`\gamma _{IV}^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j)\right)4S_1(j1)\stackrel{~}{\gamma }_+^{(0)}(j),`$ (58)
$`f_V(j)`$ $`=`$ $`(j1)\mathrm{\Delta }f_g(j)+{\displaystyle \frac{j+2}{2}}\mathrm{\Delta }f_q(j)\mathrm{\Delta }f_{q,g}^{}(j),`$
$`\gamma _V^{(0)}(j)`$ $`=`$ $`4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j+2)\right)4S_1(j+1)\stackrel{~}{\gamma }_{}^{(0)}(j),`$ (59)
Thus, we have one supermultiplet of operators with the same anomalous dimension $`\gamma ^{LLA}(j)`$ proportional to $`\mathrm{\Psi }(1)\mathrm{\Psi }(j1)`$ 777A similar result for another type of operators was obtained in ref. .. The momenta of the corresponding linear combinations of parton distributions can be obtained from the above expressions $`f_k(j)`$ by an appropriate shift of their argument $`j`$ in accordance with the corresponding shift of the argument of $`\gamma _k(j)`$. Moreover, the coefficients in these linear combinations for $`N=4`$ SUSY can be found from the super-conformal invariance (cf. Ref ) and should be the same for all orders of the perturbation theory in an appropriate renormalization scheme.

The momenta of three multiplicatively renormalizable twist-2 operators for the unpolarized case are

$$f_N(j)=a_gf_g(j)+a_qf_q(j)+a_\phi f_\phi (j)$$

where the coefficients $`a_i`$ can be extracted from above expressions (55)-(59). If we insert this anzatz in the DGLAP equations (50) the following representations for the corresponding anomalous dimensions

$`\gamma _N^{(0)}(j)`$ $`=`$ $`\gamma _{gg}^{(0)}(j)+{\displaystyle \frac{a_q}{a_g}}\gamma _{qg}^{(0)}(j)+{\displaystyle \frac{a_\phi }{a_g}}\gamma _{\phi g}^{(0)}(j)`$ (60)
$`=`$ $`\gamma _{qq}^{(0)}(j)+{\displaystyle \frac{a_g}{a_q}}\gamma _{gq}^{(0)}(j)+{\displaystyle \frac{a_\phi }{a_q}}\gamma _{\phi q}^{(0)}(j)`$
$`=`$ $`\gamma _{\phi \phi }^{(0)}(j)+{\displaystyle \frac{a_g}{a_\phi }}\gamma _{g\phi }^{(0)}(j)+{\displaystyle \frac{a_q}{a_\phi }}\gamma _{q\phi }^{(0)}(j)`$
can be obtained. Eqs.(60) lead to relations among the anomalous dimension matrix $`\gamma _{ab}^{(0)}(j)`$ $`(a,b=g,q,\phi )`$ which should be valid also in the NLO approximation up to effects of breaking the superconformal invariance (see  and references therein).

Analogously for the set of two multiplicatively renormalizable operators in the polarized case

$`\mathrm{\Delta }f_N(j)`$ $`=`$ $`\stackrel{~}{a}_g\mathrm{\Delta }f_g(j)+\stackrel{~}{a}_q\mathrm{\Delta }f_q(j)`$
we can derive the following relations

$`\stackrel{~}{\gamma }_N^{(0)}(j)=\stackrel{~}{\gamma }_{gg}^{(0)}(j)+{\displaystyle \frac{\stackrel{~}{a}_q}{\stackrel{~}{a}_g}}\stackrel{~}{\gamma }_{qg}^{(0)}(j)=\stackrel{~}{\gamma }_{qq}^{(0)}(j)+{\displaystyle \frac{\stackrel{~}{a}_g}{\stackrel{~}{a}_q}}\stackrel{~}{\gamma }_{gq}^{(0)}(j).`$ (61)
So, we have nine equations for the matrix elements in the case of the usual partonic distributions and four equations for the polarized distributions, which determines completely the anomalous dimension matrices $`\gamma _{ab}^{(0)}(j)`$ $`(a,b=g,q,\phi )`$ and $`\stackrel{~}{\gamma }_{ab}^{(0)}(j)`$ $`(a,b=g,q)`$ in terms of their eigenvalues in LLA

$$\gamma _\pm ^{(0)}(j)=4S_1(j2),\gamma _0^{(0)}(j)=4S_1(j),\stackrel{~}{\gamma }_\pm ^{(0)}(j)=4S_1(j1).$$

This procedure is considered in details in Appendix C.

### 5.3 NLO anomalous dimensions and twist-two operators with the multiplicative renormalization

We have the following initial information to predict the NLO anomalous dimensions of twist-two operators with the multiplicative renormalization in N=4 SUSY.

1. As it was shown in the previous subsections, the LLA anomalous dimensions are meromorphic functions having the poles at $`j=r,r=1,0,1,\mathrm{}`$. Moreover, there is only one basic anomalous dimension $`\gamma ^{LLA}(j)`$ and all others can be obtained as $`\gamma ^{LLA}(j+m)`$, where $`m`$ is an integer number. It is useful to choose (see Eq. (14)):

$$\gamma ^{LLA}(j)=4\left(\mathrm{\Psi }(1)\mathrm{\Psi }(j1)\right)4S_1(j2).$$

Then, $`\gamma ^{LLA}(j)`$ has a pole at $`j1`$ and vanishes at $`j=2`$. One should keep the above universality and linear relations among the matrix elements also for the NLO anomalous dimensions $`\gamma _{ab}^{(1)}(j)`$ and $`\stackrel{~}{\gamma }_{ab}^{(1)}(j)`$ because it is a consequence of the super-conformal invariance. It means, that we should construct only the basic NLO anomalous dimension $`\gamma ^{NLO}(j)`$.

2. There are known results for the NLO corrections to the QCD anomalous dimensions.

3. In the $`\overline{MS}`$-scheme with the coupling constant $`\overline{a}`$ and in the $`\overline{MS}`$-like scheme with the coupling constant $`\widehat{a}`$ (i.e. in the scheme based on DRED procedure), the terms $`\zeta (2)`$ should disappear in the final result for the forward Compton scattering (see -). Therefore the terms $`\zeta (2)`$ are cancelled at even $`j`$ in anomalous dimensions for the structure functions $`F_2`$ and $`F_L`$ (related to the unpolarized parton distributions) and at odd $`j`$ in anomalous dimensions for structure functions $`g_1`$ and $`F_3`$ (related to the polarized parton distributions).

4. From the BFKL equation in the framework of DRED scheme (see (20), (21) and (22)) we know, that there is no mixing among the functions of different transcendentality levels $`i`$ 888Note that similar arguments were used also in  to obtain analytic results for contributions of some complicated massive Feynman diagrams. The method was based on a direct calculation of several terms in the series over an inverse mass with taking into account its basic structure found earlier in  by considering few special diagrams with the use of the differential equation method ., i.e. all special functions at the NLO correction contain the sums of the terms $`1/n^i(i=3)`$. More precisely, if one will introduce the transcendentality level of the functions in accordance with the complexity of the terms in the corresponding sums

$$\mathrm{\Psi }1/n,\mathrm{\Psi }^{}\beta ^{}\zeta (2)1/n^2,\mathrm{\Psi }^{\prime \prime }\beta ^{\prime \prime }\zeta (3)1/n^3,$$

then for the BFKL equation in LLA and in NLO the corresponding levels are $`i=1`$ and $`i=3`$, respectively.

Because in N=4 SUSY there is a relation between the BFKL and DGLAP equations, the similar properties are assumed to be valid for anomalous dimensions themselves, i.e. the basic functions $`\gamma ^{LLA}(j)`$ and $`\gamma ^{NLO}(j)`$ should be of the types $`1/j^i`$ with the levels $`i=1`$ and $`i=3`$, respectively. The only exception could be for the terms appearing in the Born approximation, because such contribution can be removed by an approximate finite renormalization of the coupling constant. The LLA basic anomalous dimension is given above in terms of $`S_1(j2)`$. Then, the NLO basic anomalous dimension $`\gamma ^{NLO}(j)`$ can be expressed through the functions:

$$S_{\pm i}(j2),S_{\pm k,l}(j2),\zeta (k)S_{\pm l}(j2),\zeta (i)$$

(here $`i=3`$ and $`k+l=i`$), where

$`S_i(j)`$ $`=`$ $`{\displaystyle \underset{m=1}{\overset{j}{}}}{\displaystyle \frac{1}{m^i}}`$
$`S_{2,1}(j)`$ $`=`$ $`{\displaystyle \underset{m=1}{\overset{j}{}}}{\displaystyle \frac{1}{m^2}}S_1(m),`$ (62)
$`S_2(j)`$ $`=`$ $`(1)^j{\displaystyle \underset{m=1}{\overset{j}{}}}{\displaystyle \frac{(1)^m}{m^2}}\left(1(1)^j\right){\displaystyle \frac{1}{2}}\zeta (2),`$
$`S_3(j)`$ $`=`$ $`(1)^j{\displaystyle \underset{m=1}{\overset{j}{}}}{\displaystyle \frac{(1)^m}{m^3}}\left(1(1)^j\right){\displaystyle \frac{3}{4}}\zeta (3),`$
$`S_{2,1}(j)`$ $`=`$ $`(1)^j{\displaystyle \underset{m=1}{\overset{j}{}}}{\displaystyle \frac{(1)^m}{m^2}}S_1(m)\left(1(1)^j\right){\displaystyle \frac{5}{8}}\zeta (3)`$ (63)
The functions $`S_2(j)`$, $`S_3(j)`$ and $`S_{2,1}(j)`$, which were introduced recently in , coincide (up to an opposite sign) with the functions $`K_2(j)`$, $`K_3(j)`$ and $`K_{2,1}(j)`$ considered in  and used in our previous paper , i.e. $`S_i(j)=K_i(j)`$, $`S_{2,1}(j)=K_{2,1}(j)`$. Note, however, that our definition (63) of $`S_i(j)`$ and $`S_{k,l}(j)`$ coincides with the functions of  only at even values of $`j`$. Expressions (63) after the shift of the summation index $`mmj`$ can be analytically continuated from even to complex values of $`j`$, which reproduces $`\beta ^{}(j+1)`$ (23) and associated functions.

Note that the terms $`\zeta (2)`$ should be absent for even values of $`j`$ and for odd values of $`j`$ in the unpolarized and in polarized cases, respectively, in accordance with item 3.

Moreover, there is an important observation: the function

$$S_1(j)=(1)^j\underset{m=1}{\overset{j}{}}\frac{(1)^m}{m}\left(1(1)^j\right)\mathrm{ln}2$$

does not contribute to the QCD anomalous dimensions and the Wilson coefficient functions (see  and , respectively). Further, as it was shown in , the terms $`S_1(j)`$ cancel in the final results for the diagrams describing the longitudinal Wilson coefficient function.

5. The terms

$`S_{\pm l}(j2)/(j\pm m)^k(j+l=i)`$ (64)
are absent in the anomalous dimension $`\gamma ^{NLO}(j)`$. There are two reasons for this conclusion.

Firstly, these terms have additional poles in the points $`j=m`$. But such poles should cancel, if we start with the BFKL equation and obtain $`\gamma ^{NLO}(j)`$ by the analytic continuation to $`|n|=1r`$. Indeed, using such procedure one can not obtain the doublets of poles.

The second reason comes from the consideration of the multiplicatively renormalizable linear combinations (58) and (59). If, for example, in the polarized case functions (64) contribute, then we shall have the terms $`(j\pm m)^1`$ in one combination and the terms $`(j\pm m\pm 2)^1`$ in another combination. However, from the direct calculation of the NLO anomalous dimensions in the polarized case (see ) we know that only the terms $`j^1`$ and $`(j+1)^1`$ are present in these combinations.

So, terms (64) should be absent in the universal NLO anomalous dimension in the N=4 SUSY case.

Note, that the sums $`S_2(j)`$, $`S_3(j)`$ and $`S_{2.1}(j)`$ (see Eqs.(63)) appear explicitly in calculations only at even values of $`j`$ in the unpolarized case (and/or at odd values of $`j`$, after the replacement $`(1)^j(1)^{j+1}`$, in the polarized case). The analytic continuation to complex values of $`j`$ for the functions $`S_{\pm i}(j)`$ ($`i=1,2,3`$) and $`S_{2.1}(j)`$ in Eqs. (62) and (63) can be done easily (see ) by a replacement of the sums $`_{m=1}^j`$ in the r.h.s with the difference $`_{m=1}^{\mathrm{}}_{m=1+j}^{\mathrm{}}`$. They are expressed in terms of the polygamma- and associated functions:

$`S_1(j)`$ $`=`$ $`\mathrm{\Psi }(j+1)\mathrm{\Psi }(1),`$
$`S_i(j)`$ $`=`$ $`{\displaystyle \frac{(1)^{i1}}{(i1)!}}\left[\mathrm{\Psi }^{i1}(j+1)\mathrm{\Psi }^{i1}(1)\right]={\displaystyle \frac{(1)^{i1}}{(i1)!}}\mathrm{\Psi }^{i1}(j+1)+\zeta (i)(i>1),`$
$`S_{2,1}(j)`$ $`=`$ $`2\zeta (3){\displaystyle \underset{m=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{S_1(m+j)}{(m+j)^2}}=2\zeta (3){\displaystyle \underset{m=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{1}{(m+j)^2}}\left[\mathrm{\Psi }(j+1)\mathrm{\Psi }(1)\right],`$ (65)
$`S_1(j)`$ $`=`$ $`\beta (j+1)\beta (1)=\beta (j+1)+\mathrm{ln}2,`$
$`S_i(j)`$ $`=`$ $`{\displaystyle \frac{(1)^{i1}}{(i1)!}}\left[\beta ^{i1}(j+1)\beta ^{i1}(1)\right]={\displaystyle \frac{(1)^{i1}}{(i1)!}}\beta ^{i1}(j+1)+\zeta (i)(i>1),`$
$`S_{2,1}(j)`$ $`=`$ $`{\displaystyle \underset{m=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(1)^mS_1(m+j+1)}{(m+j+1)^2}}{\displaystyle \frac{5}{8}}\zeta (3)`$ (66)
$`={\displaystyle \underset{m=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(1)^m}{(m+j+1)^2}}\left[\mathrm{\Psi }(m+j+2)\mathrm{\Psi }(1)\right]{\displaystyle \frac{5}{8}}\zeta (3)\stackrel{~}{\beta }^{^{\prime \prime }}(j+1){\displaystyle \frac{5}{8}}\zeta (3),`$
where $`\zeta (r)=\left(2^{1r}1\right)\zeta (r)(r>1),\zeta (1)=\mathrm{ln}2`$.

6. Further, the NLO anomalous dimension $`\gamma ^{NLO}(j)`$ is equal to a combination of the most complicated contributions (i.e. the functions with a maximal value of the transcendentality level $`i=3`$) for the QCD anomalous dimensions (with the SUSY relation for the QCD color factors $`C_F=C_A=N_c`$).

These most complicated contributions (with $`i=3`$) are the same for all QCD anomalous dimensions (coinciding in N=4 SUSY)  (only the NLO scalar-scalar anomalous dimension is not known yet):

$`\gamma _{qq}^{(1)QCD}(j)`$ $``$ $`\gamma _{gg}^{(1)QCD}(j)\stackrel{~}{\gamma }_{qq}^{(1)QCD}(j)\stackrel{~}{\gamma }_{gg}^{(1)QCD}(j)`$ (67)
$`=`$ $`16\left({\displaystyle \frac{\alpha _s}{4\pi }}\right)^2N_c^2Q(j2)+\mathrm{}(C_F=C_A=N_c),`$
where we omit less complicated contributions and

$`Q(j)={\displaystyle \frac{1}{2}}\left(S_3(j)+S_3(j)\right)+S_1(j)\left(S_2(j)+S_2(j)\right)S_{2,1}(j).`$ (68)
Using Eqs. (65) and (66), the function $`Q(j)`$ can be rewritten in terms of the polygamma- and associated functions:

$`Q(j)`$ $`=`$ $`{\displaystyle \frac{3}{4}}\zeta (3)\stackrel{~}{\beta }^{^{\prime \prime }}(j+1)+{\displaystyle \frac{1}{4}}\left(\mathrm{\Psi }^{^{\prime \prime }}(j+1)+\beta ^{^{\prime \prime }}(j+1)\right)`$ (69)
$`+`$ $`\left(\mathrm{\Psi }(j+1)\mathrm{\Psi }(1)\right)\left({\displaystyle \frac{1}{2}}\zeta (2)\mathrm{\Psi }^{^{}}(j+1)\beta ^{^{}}(j+1)\right)`$
$`=`$ $`{\displaystyle \frac{3}{4}}\zeta (3)\stackrel{~}{\beta }^{^{\prime \prime }}(j+1)+{\displaystyle \frac{1}{16}}\mathrm{\Psi }^{^{\prime \prime }}\left({\displaystyle \frac{j+1}{2}}\right)`$
$`+`$ $`{\displaystyle \frac{1}{2}}\left(\mathrm{\Psi }(j+1)\mathrm{\Psi }(1)\right)\left(\zeta (2)\mathrm{\Psi }^{^{}}\left({\displaystyle \frac{j+1}{2}}\right)\right).`$
Thus, for $`N=4`$ SUSY the NLO universal anomalous dimension $`\gamma ^{NLO}(j)`$ has the form

$`\gamma ^{NLO}(j)=16Q(j2)`$ (70)
7. We could add the term $`\zeta (3)`$ to the r.h.s. of (68), but due to the condition
$`\gamma ^{NLO}(j=2)=0`$ it cancels.

So, the universal anomalous dimension $`\gamma (j)`$ in two first orders of the perturbation theory for N=4 SUSY is

$`\gamma (j)\gamma ^{uni}(j)=\widehat{a}\gamma ^{LLA}(j)+\widehat{a}^2\gamma ^{NLO}(j),`$ (71)
where $`\gamma ^{LLA}(j)`$ and $`\gamma ^{NLO}(j)`$ are given by Eqs. (14) and (70), respectively. All other anomalous dimensions can be obtained as $`\gamma ^{uni}(j+m)`$, where $`m`$ is an integer number.

Thus, the above arguments allow us to construct the NLO corrections to anomalous dimensions in N=4 SUSY, which were unknown earlier. We check these results by direct calculations  and reproduce the anomalous dimension $`\gamma ^{NLO}(j)`$. Note, however, that in  the coupling constant $`\overline{a}`$ (see Eq.(24)) is used, which is responsable for an appearence of the additional contribution $`1/3\gamma ^{LO}(j)`$ to the above obtained NLO anomalous dimension (70). This contribution can be absorbed in $`\widehat{a}`$ by a finite renormalization of the coupling constant (see the subsection 3.1). Note also that in the dimensional reqularization scheme the eigenvalues of this matrix are not expressed only in terms od the function $`\gamma ^{uni}(j)`$ with the shift $`jj+m`$ of its argument.

### 5.4 DGLAP evolution

Using our knowledge of the anomalous dimensions we can construct the solution of the DGLAP equation in the Mellin moment space in the framework of N=4 SUSY.

A. Polarized case

The polarized parton distributions are splitted in two contributions:

$`\mathrm{\Delta }f_{q,g}(j,Q^2)=\mathrm{\Delta }f_{q,g}^+(j,Q^2)+\mathrm{\Delta }f_{q,g}^{}(j,Q^2),`$ (72)
where

at LO

$`\mathrm{\Delta }f_{q,g}^\pm (j,Q^2)=\mathrm{\Delta }f_{q,g}^{\pm ,LO}(j,Q_0^2)\left({\displaystyle \frac{Q^2}{Q_0^2}}\right)^{\stackrel{~}{\gamma }_\pm ^{(0)}a}\left(\stackrel{~}{\gamma }_\pm ^{(0)}=4S_1(j1)\right),`$ (73)
at NLO

$`\mathrm{\Delta }f_{q,g}^\pm (j,Q^2)=\mathrm{\Delta }f_{q,g}^{\pm ,NLO}(j,Q_0^2)\left({\displaystyle \frac{Q^2}{Q_0^2}}\right)^{\left(\stackrel{~}{\gamma }_\pm ^{(0)}a+\stackrel{~}{\gamma }_{\pm \pm }^{(1)}a^2\right)}\left(\stackrel{~}{\gamma }_{\pm \pm }^{(1)}=16Q(j1)\right),`$ (74)
where

$`\mathrm{\Delta }f_{q,g}^{\pm ,NLO}(j,Q_0^2)`$ $`=`$ $`\mathrm{\Delta }f_{q,g}^{\pm ,NLO}(j,Q_0^2)\left(1{\displaystyle \frac{\stackrel{~}{\gamma }_\pm ^{(1)}\widehat{a}}{\stackrel{~}{\gamma }_\pm ^{(0)}\stackrel{~}{\gamma }_{}^{(0)}}}\right)`$ (75)
$`+`$ $`{\displaystyle \frac{\stackrel{~}{\gamma }_\pm ^{(1)}\widehat{a}}{\stackrel{~}{\gamma }_{}^{(0)}\stackrel{~}{\gamma }_\pm ^{(0)}}}\mathrm{\Delta }f_{q,g}^{,NLO}(j,Q_0^2).`$
Here the anomalous dimensions $`\stackrel{~}{\gamma }_\pm ^{(1)}`$ and $`\stackrel{~}{\gamma }_\pm ^{(1)}`$ are related to $`\stackrel{~}{\gamma }_{a,b}^{(1)}`$ $`(a,b=q,g)`$ as follows:

$`\left(\begin{array}{cc}\stackrel{~}{\gamma }_{++}^{(1)}(j)& \stackrel{~}{\gamma }_+^{(1)}(j)\\ \stackrel{~}{\gamma }_+^{(1)}(j)& \stackrel{~}{\gamma }_{}^{(1)}(j)\end{array}\right)=\widehat{\stackrel{~}{V}}^1\left(\begin{array}{cc}\stackrel{~}{\gamma }_{gg}^{(1)}(j)& \stackrel{~}{\gamma }_{qg}^{(1)}(j)\\ \stackrel{~}{\gamma }_{gq}^{(1)}(j)& \stackrel{~}{\gamma }_{qq}^{(1)}(j)\end{array}\right)\widehat{\stackrel{~}{V}}`$ (80)
with $`\widehat{\stackrel{~}{V}}`$ and $`\widehat{\stackrel{~}{V}}^1`$ given in eq.(C12).

Notice that only anomalous dimensions $`\stackrel{~}{\gamma }_{\pm \pm }^{(1)}`$ are important for $`N=4`$ SUSY at the order $`O(\widehat{a}^2)`$ because they contribute to the $`Q^2`$-evolution of parton distributions.

The anomalous dimensions $`\stackrel{~}{\gamma }_\pm ^{(1)}`$, which were not calculated in the previous section, can be found from eq.(80) provided that the gluino-gluon polarized anomalous dimensions $`\stackrel{~}{\gamma }_{ab}^{(1)}(j)`$ $`(a,b=g,q)`$ are known. The anomalous dimensions $`\stackrel{~}{\gamma }_\pm ^{(1)}`$ give contributions at the order $`O(\widehat{a}^2)`$ only to the normalization factors $`\mathrm{\Delta }f_{q,g}^{\pm ,NLO}(j,Q_0^2)`$. They appear in the $`Q^2`$-dependent part of $`\mathrm{\Delta }f_{q,g}^\pm (j,Q^2)`$ at the level $`O(\widehat{a}^3)`$ in the following form:

$$\widehat{a}^3\frac{\stackrel{~}{\gamma }_\pm ^{(1)}\stackrel{~}{\gamma }_\pm ^{(1)}}{\stackrel{~}{\gamma }_\pm ^{(0)}\stackrel{~}{\gamma }_{}^{(0)}}$$

B. Unpolarized case

The unpolarized parton distributions are splitted in the three parts:

$`f_{q,g,\phi }(j,Q^2)={\displaystyle \underset{i=+,,0}{}}f_{q,g,\phi }^i(j,Q^2),`$ (81)
where

at LO

$`f_{q,g,\phi }^i(j,Q^2)=f_{q,g,\phi }^{i,LO}(j,Q_0^2)\left({\displaystyle \frac{Q^2}{Q_0^2}}\right)^{\gamma _i^{(0)}a}\left(\gamma _\pm ^{(0)}=4S_1(j2),\gamma _0^{(0)}=4S_1(j)\right),`$ (82)
at NLO

$`f_{q,g,\phi }^i(j,Q^2)=f_{q,g,\phi }^{i,NLO}(j,Q_0^2)\left({\displaystyle \frac{Q^2}{Q_0^2}}\right)^{\left(\gamma _i^{(0)}a+\gamma _{ii}^{(1)}a^2\right)}\left(\gamma _{\pm \pm }^{(1)}=16Q(j2),\gamma _{00}^{(1)}=16Q(j)\right)`$ (83)
As in the previous case A, only the anomalous dimensions $`\gamma _{ii}^{(1)}`$ are important at the level $`O(\widehat{a}^2)`$ in N=4 SUSY. The anomalous dimensions $`\gamma _{ik}^{(1)}`$ $`(ik)`$ contribute at $`O(\widehat{a}^2)`$ level only to the normalization factors $`f_{q,g,\phi }^{i,NLO}(j,Q_0^2)`$. They are related to the qluino-gluon-scalar anomalous dimensions $`\gamma _{ab}^{(1)}(j)`$ $`(a,b=g,q,\phi )`$ as follows:

$$\left(\begin{array}{ccc}\gamma _{++}^{(1)}(j)& \gamma _{0+}^{(1)}(j)& \gamma _+^{(1)}(j)\\ \gamma _{+0}^{(1)}(j)& \gamma _{00}^{(1)}(j)& \gamma _0^{(1)}(j)\\ \gamma _+^{(1)}(j)& \gamma _0^{(1)}(j)& \gamma _{}^{(1)}(j)\end{array}\right)=\widehat{V}^1\left(\begin{array}{ccc}\gamma _{gg}^{(1)}(j)& \gamma _{qg}^{(1)}(j)& \gamma _{\phi g}^{(1)}(j)\\ \gamma _{gq}^{(1)}(j)& \gamma _{qq}^{(1)}(j)& \gamma _{\phi q}^{(1)}(j)\\ \gamma _{g\phi }^{(1)}(j)& \gamma _{q\phi }^{(1)}(j)& \gamma _{\phi \phi }^{(1)}(j)\end{array}\right)\widehat{V},$$

where $`\widehat{V}`$ and $`\widehat{V}^1`$ are given by eq.(C36). The anomalous dimensions $`\gamma _{ik}^{(1)}`$ $`(ik)`$ can be obtained from the above equation provided that the gluino-gluon-scalar anomalous dimension matrix $`\gamma _{ab}^{(1)}(j)`$ $`(a,b=g,q,\phi )`$ is known.

## 6 Relation between the DGLAP and BFKL equations

As we pointed out in subsection 3.6, in the case of $`N=4`$ SUSY the BFKL results (20) and (21) are analytic in $`|n|`$ and one can continue the eigenvalues to the negative values of $`|n|`$. It gives a possibility to find the singular contributions to anomalous dimensions of the twist-2 operators not only at $`j=1`$ but also at other integer non-physical points $`j=0,1,2\mathrm{}`$. In the Born approximation for the anomalous dimension of the supermultiplet of the twist-2 operators we obtain $`\gamma ^{uni}=4\widehat{a}(\mathrm{\Psi }(1)\mathrm{\Psi }(j1))`$ which coincides with the result of the direct calculations (see  and the discussions in subsection 4.2 ). Thus, in the case of $`N=4`$ the BFKL equation presumably contains the information sufficient for restoring the kernel of the DGLAP equation. Below we investigate the relation between these equations in the NLO approximation.

### 6.1 DGLAP approach

Let us start with an investigation of singularities of the anomalous dimensions of twist-2 operators which were obtained in a direct way at the previous Section.

By presenting the Lorentz spin $`j`$ as $`\omega r`$, where $`r=1,0,1,\mathrm{}`$ and pushing $`\omega 0`$ one can calculate the singular behavior of the universal anomalous dimension $`\gamma ^{uni}(j)`$ (71). Note, that in our discussion of the BFKL equation we used the more general definition $`j=1+\omega +|n|`$ for the rank of the Lorentz tensor $`O_{\mu _1\mathrm{}\mu _j}`$ (see (12)), where $`1+\omega `$ is the number of its longitudinal indices and $`|n|`$ is the conformal spin. Therefore strictly speaking in all expressions discussed in this subsection $`\omega `$ should be substituted by $`j+r`$.

For the special functions contributed to $`\gamma ^{LLA}(j)`$ and $`\gamma ^{NLO}(j)`$ at $`j=\omega rr,r0`$ we have the following expansion in $`\omega `$:

$`S_1(j2)`$ $`=`$ $`{\displaystyle \frac{1}{\omega }}+S_1(r+1){\displaystyle \underset{l=1}{\overset{\mathrm{}}{}}}\omega ^l\left[(1)^l\zeta (l+1)S_{l+1}(r+1)\right]`$
$`S_2(j2)`$ $`=`$ $`{\displaystyle \frac{1}{\omega ^2}}S_2(r+1){\displaystyle \underset{l=1}{\overset{\mathrm{}}{}}}(l+1)\omega ^l\left[(1)^l\zeta (l+2)+S_{l+2}(r+1)\right]`$
$`S_3(j2)`$ $`=`$ $`{\displaystyle \frac{1}{\omega ^3}}+S_3(r+1)`$ (84)
$``$ $`{\displaystyle \underset{l=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(l+1)(l+2)}{2}}\omega ^l\left[(1)^l\zeta (l+3)S_{l+3}(r+1)\right]`$
$`S_2(j2)`$ $`=`$ $`(1)^{r+1}\{{\displaystyle \frac{1}{\omega ^2}}+\zeta (2)(1(1)^r)+S_2(r+1)`$
$`+`$ $`{\displaystyle \underset{l=1}{\overset{\mathrm{}}{}}}(l+1)\omega ^l[(1)^l\zeta ((l+2))+S_{(l+2)}(r+1)]\}`$
$`S_3(j2)`$ $`=`$ $`(1)^{r+1}\{{\displaystyle \frac{1}{\omega ^3}}+\zeta (3)(1(1)^r)S_3(r+1)`$
$`+`$ $`{\displaystyle \underset{l=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(l+1)(l+2)}{2}}\omega ^l[(1)^l\zeta ((l+3))S_{(l+3)}(r+1)]\}`$
$`S_{2,1}(j2)`$ $`=`$ $`(1)^{r+1}\{{\displaystyle \frac{1}{\omega }}[\zeta (2)S_2(r+1)]`$ (85)
$`\zeta (3)+{\displaystyle \frac{5}{8}}\zeta (3)(1(1)^r)S_{2,1}(r+1)S_3(r+1)\}.`$
Presenting our results for $`\gamma ^{LLA}(j)`$ and $`\gamma ^{NLO}(j)`$ in terms of these functions, we obtain

$`\gamma ^{LLA}(j)`$ $`=`$ $`4\left[{\displaystyle \frac{1}{\omega }}S_1(r+1)\widehat{S}_2(r+1)\omega +O(\omega ^2)\right]`$ (86)
$`\gamma ^{NLO}(j)`$ $`=`$ $`8[{\displaystyle \frac{(1+(1)^r)}{\omega ^3}}2S_1(r+1){\displaystyle \frac{(1+(1)^r)}{\omega ^2}}`$ (89)
$``$ $`((1+(1)^r)\zeta (2)+2(1)^rS_2(r+1)){\displaystyle \frac{1}{\omega }}+O(\omega ^0)]`$
$`=\{\begin{array}{cc}16\left[\frac{1}{\omega ^3}2S_1(r+1)\frac{1}{\omega ^2}\widehat{S}_2(r+1)\frac{1}{\omega }+O(\omega ^0)\right]& \text{ if }r=2k\hfill \\ 16\left[S_2(r+1)\frac{1}{\omega }+O(\omega ^0)\right]& \text{ if }r=2k+1,\hfill \end{array}`$
where $`\widehat{S}_2(r)=\zeta (2)+S_2(r)`$. So, the double-logarithmic poles $`\omega ^3`$ appear in the case of even values of $`r`$ (see discussions in Section 2).

Note that the functions $`S_1(r+1)`$ and $`S_2(r+1)`$ do not contribute to $`\gamma ^{LLA}(j)`$ in the limit $`jr,r0`$. The absence of $`S_1(r+1)`$ is explained by the fact, that the quantity $`S_1(j2)`$ does not appear in the basic NLO anomalous dimension (see discussion in the subsection 5.3, item 4).

On the other hand, the absence of the term $`S_2(r+1)`$ is related to the following important property of $`\gamma `$. As it follows from eq.(85), the functions $`S_2(j2)`$, $`S_3(j2)`$ and $`S_{2,1}(j2)`$ giving contributions to $`\gamma ^{NLO}(j)`$ (see eqs. (70) and (68)) are responsible for the different asymptotics of it at $`jr,r>0`$ for even and odd $`r`$ (see the r.h.s. of eq.(89)). But the combination of these functions $`S_2(j2)`$, $`S_3(j2)`$ and $`S_{2,1}(j2)`$ contributes to $`\gamma ^{NLO}(j)`$ in such way, that the function $`S_2(r+1)`$ is absent in the r.h.s. of eq.(89).

Thus, for $`\gamma ^{uni}(j)`$ we obtain at $`j=\omega rr`$ ($`r0`$) 999Because really the expansion parameter is $`4\widehat{a}`$ (see, for example, Eqs. (86) and (89), we shall consider below $`G=4\widehat{a}`$ as a coupling constant.

$`\gamma ^{uni}(j)`$ $`=`$ $`G\left[{\displaystyle \frac{1}{\omega }}S_1(r+1)\omega \widehat{S}_2(r+1)+O(\omega ^2)\right]`$ (92)
$`+`$ $`G^2\{\begin{array}{cc}\frac{1}{\omega ^3}2S_1(r+1)\frac{1}{\omega ^2}\widehat{S}_2(r+1)\frac{1}{\omega }+O(\omega ^0)& \text{ if }r=2k\hfill \\ S_2(r+1)\frac{1}{\omega }+O(\omega ^0)& \text{ if }r=2k+1.\hfill \end{array}`$
Note that $`\gamma ^{uni}(j)`$ at $`j=\omega rr`$ can be presented as a solution of the following Bethe-Salpeter equations

$`1`$ $`=`$ $`G\left[{\displaystyle \frac{1\omega S_1(r+1)\omega ^2\widehat{S}_2(r+1)}{\gamma (\omega \gamma )}}S_1^2(r+1)\right]+\mathrm{}`$ (93)
for $`r=2k0`$ and

$`1`$ $`=`$ $`G[{\displaystyle \frac{1}{\gamma \omega }}{\displaystyle \frac{S_1(r+1)+\omega \widehat{S}_2(r+1)}{\gamma }}+S_2(r+1)]+\mathrm{}`$ (94)
for $`r=2k+1>0`$.

As it was shown in the subsection 5.4 the anomalous dimensions, responsible for the $`Q^2`$ evolution of parton distributions in the framework of $`N=4`$ SUSY model, have the following form

in unpolarized case

$`\gamma _\pm (j)`$ $`=`$ $`G\gamma _\pm ^{(0)}(j)+G^2\gamma _{\pm \pm }^{(1)}(j)+\mathrm{}=\gamma ^{uni}(j+22)`$
$`\gamma _0(j)`$ $`=`$ $`G\gamma _0^{(0)}(j)+G^2\gamma _{00}^{(1)}(j)+\mathrm{}=\gamma ^{uni}(j+2)`$ (95)
and in polarized case

$`\stackrel{~}{\gamma }_\pm (j)`$ $`=`$ $`G\stackrel{~}{\gamma }_\pm ^{(0)}(j)+G^2\stackrel{~}{\gamma }_{\pm \pm }^{(1)}(j)+\mathrm{}=\gamma ^{uni}(j+12).`$ (96)
Thus, for unpolarized and polarized cases, the corresponding anomalous dimensions have the double-logarithmic poles $`\omega ^3`$ at even and odd $`r`$, respectively. The arguments of regular terms and of the functions in front of poles $`\omega ^i`$ are shifted by an integer number. By chosing certain values of this number we obtain various $`\gamma _s(j)`$.

### 6.2 BFKL approach

Let us investigate a possibility to obtain the residues of the anomalous dimension $`\gamma ^{uni}(j)`$ in the poles at $`j=r`$ ($`r=0,1,2,\mathrm{}`$ ) from the BFKL equation. Its eigenvalue $`\omega `$ can be analytically continued to the negative points $`|n|=r1`$. One can present $`\omega `$ as a solution of the equation (see (44))

$$\omega =G\left(2\mathrm{\Psi }(1)\mathrm{\Psi }(\gamma )\mathrm{\Psi }(|n|+1+\omega \gamma )+\epsilon \right),$$

having the property of the hermitian separability corresponding to the symmetry between $`\gamma `$ and $`|n|+1+\omega \gamma `$ (see (46)). The expression for $`\epsilon `$ is given in Eq. (45).

To begin with, we note that providing that $`j=|n|+1+\omega `$ (see Eq. (12)) the r.h.s. of the BFKL equation contains the $`\mathrm{\Psi }`$-function with its argument equal to $`j`$ at $`\gamma 0`$, whereas $`\gamma ^{LLA}(j)`$ obtained from the DGLAP equation contains $`\mathrm{\Psi }(j1)`$. The functions contributing to the NLO correction have an analogous shift of their argument in comparison with the functions appearing in the NLO corrections to $`\gamma ^{uni}(j)`$ (92).

In the Born approximation $`\gamma ^{LLA}G`$ the reason for the difference of arguments of $`\mathrm{\Psi }`$-functions appearing in the DGLAP and BFKL approaches can be easily understood. Starting from the BFKL equation one reproduces the universal anomalous dimension $`\gamma ^{uni}`$ of twist-two operators by an analytic continuation of its eigenvalue $`\omega `$ to the points $`j=|n|+1+\omega =\omega r`$ ($`r=0,1,\mathrm{}`$). The residue of the pole of $`\gamma ^{uni}`$ at $`j=1`$ is calculated also in accordance with the DGLAP equation but without any continuation (see (40), (41)). In the BFKL approach the term regular at $`jr`$ does not contain the contribution from the pole at $`j=1`$ and equals $`GS_1(r)`$ instead of $`GS_1(r+1)`$. It is important, that the anomalous dimension obtained from the BFKL equation depends on two parameters $`\omega `$ and $`\left|n\right|`$ and therefore we can impose a certain relation (see (13)) between two asymptotics $`\omega 0`$ and $`\left|n\right|+r+10`$. Another ambiguity is related to the existence of double-logarithmic terms near $`j=r`$ (see below). As a result, in particular in LLA we can obtain an agreement between the regular contribution to $`\gamma ^{uni}`$ and its pole singularities.

In the next-to-leading approximation the situation is similar but not so transparent. As in the case of LLA to obtain the residues of the poles for $`\gamma ^{uni}(j)`$ in the BFKL approach, we continue the eigenvalue $`\omega `$ to $`j=r`$ ($`r=0,1,\mathrm{}`$). But this continuation does not contain the regular contributions from the singularity at $`j=1`$. Another problem is that to calculate the residues of the non-leading poles $`G^2\omega ^2`$ and $`G^2\omega ^1`$ at $`\omega j+r0`$ for $`\gamma ^{uni}`$ one should continue the BFKL equation from the region of its applicapability $`\omega \gamma `$ to the double-logarithmic region $`\gamma \omega `$. The property of the hermitian separability (46) of the BFKL kernel is helpful for this purpose, because it relates the small parameters in one combination $`\omega +|n|r1\gamma `$. It turns out, that in the intermediate region $`\gamma \omega `$ we should rearrange the BFKL equation in such way that the Born contribution will contain the double-logarithmic terms. This resummation procedure leads to a renormalization of the next-to-leading correction $`G^2`$. Namely, to avoid the double-counting we should subtruct from this correction the sum of poles containing the square of the most singular Born contribution. It turns out, that the subtracted terms are small in the BFKL region $`\omega \gamma `$, but they are essential for the non-leading poles $`G^2\omega ^2`$ and $`G^2\omega ^1`$ of the anomalous dimension. These subtructions can be found by calculating the next-to-leading contribution to the BFKL kernel with a greater accuracy in comparison with ref. , but we leave this problem for future publications. Below we discuss in more details the origin of above ambiguities.

Let us write down the eigenvalue of the BFKL kernel near the singularity at $`\gamma =0`$ for the physical conformal spins $`|n|=0,2,\mathrm{}`$ (see Eqs. (42), (43) and (47))

$`\omega =G\left({\displaystyle \frac{1}{\gamma }}\mathrm{\Psi }(1+|n|)+\mathrm{\Psi }(1)\right)+G^2\left({\displaystyle \frac{\mathrm{\Psi }(1)\mathrm{\Psi }(1+|n|)}{\gamma ^2}}+{\displaystyle \frac{c(n)}{2\gamma }}\right).`$ (97)
¿From this equation one can obtain easily the anomalous dimensions (49) at integer positive $`\left|n\right|`$. The analytic continuation of $`\omega `$ to the negative integer points $`|n|r1`$ contains divergencies due to the following formulae

$$\mathrm{\Psi }(1)\mathrm{\Psi }(|n|+1)\frac{1}{|n|+r+1}S_1(r)(|n|+r+1)\widehat{S}_2(r),$$

$$c(n)\frac{1+(1)^r}{\left(|n|+r+1\right)^2}+S_2(r)+(1)^rS_2(r)\frac{1}{2}\left(1+(1)^r\right)\zeta (2),$$

where $`\widehat{S}_2(r)`$ was introduced after Eq. (89).

We interprete the divergencies of $`\omega `$ at $`|n|r1`$ for even $`r=0,2,\mathrm{}`$ as a manifestation of the double-logarithmic contribution $`\mathrm{\Delta }\gamma G^2/\omega ^3`$. Indeed, one can present expression (97) in the kinematical region $`\omega \gamma 1`$ with the same accuracy as a solution of the equation

$$\omega =G\left(\frac{1}{\gamma }\mathrm{\Psi }(j\gamma )+\mathrm{\Psi }(1)\right)+G^2\left(\frac{\mathrm{\Psi }(1)\mathrm{\Psi }(j\gamma )}{\gamma ^2}+\frac{c(j1\gamma )}{2\gamma }\right),$$

where $`j=1+\omega +\left|n\right|`$. It is obvious, that if one neglects the term $`G^2`$ the anomalous dimension $`\gamma `$ contains the double-logarithmic contribution $`\mathrm{\Delta }\gamma G^2/\omega ^3`$. Note, that the above expression corresponds to a small-$`\gamma `$ asymptotics of the BFKL equation written in the separable form (cf. (46))

$`1`$ $`=`$ $`G\left({\displaystyle \frac{2\mathrm{\Psi }(1)\mathrm{\Psi }(\gamma )\mathrm{\Psi }(j\gamma )}{\omega }}+{\displaystyle \frac{c(\gamma 1)+c(j1\gamma )}{2}}+\zeta (2)\right)`$
$`+`$ $`{\displaystyle \frac{G^2}{4}}\left(\varphi (\gamma )+\varphi (j\gamma )\right),`$
where $`\varphi (\gamma )`$ is presented in Eq. (28). On the other hand in the Born approximation at $`jr`$ and $`\gamma 0`$, neglecting the corrections to the relation $`\omega =j+r`$, one can simplify the previous equation as follows

$$1G\left(\frac{1}{\gamma (\omega \gamma )}\frac{S_1(r)+(\omega \gamma )\widehat{S}_2(r)}{\omega }\right).$$

It will be shown below, that for even $`r`$ the double-logarithmic terms $`\mathrm{\Delta }\gamma G^2/\omega ^3`$ obtained from this equation are not cancelled. It is important, that the next-to-leading correction $`G^2`$ to $`\omega `$ was obtained in ref.  by subtracting the first iteration of the Born kernel from the contribution of all one-loop diagrams to avoid the double-counting. At $`\left|n\right|r1`$ we chose another Born approximation which contains the double-logarithmic terms. It means, that now one should subtract from the next-to-leading correction some terms including the most singular pole which appears in the first iteration of this Born contribution to avoid the double-counting. But the subtracted terms should be small in the applicapability region $`\omega \gamma `$ for the BFKL equation. Thus, the coefficient in front of $`G^2`$ after this subtraction can be written for even $`r`$ as follows

$`{\displaystyle \frac{(\omega \gamma )^1S_1(r)(\omega \gamma )\widehat{S}_2(r)}{\omega \gamma ^2}}+{\displaystyle \frac{2(\omega \gamma )^2+S_2(r)+S_2(r)2\zeta (2)}{2\omega \gamma }}`$
$`{\displaystyle \frac{1}{(\omega \gamma )^2\gamma ^2}}+{\displaystyle \frac{a(r)}{(\omega \gamma )\gamma ^2}}+{\displaystyle \frac{b(r)}{(\omega \gamma )^2\gamma }}+{\displaystyle \frac{d(r)}{(\omega \gamma )\gamma }}+{\displaystyle \frac{e(r)}{\gamma ^2}},`$
where the residues $`a(r),b(r),d(r)`$ and $`e(r)`$ are some unknown functions. The second order pole $`(\omega \gamma )^2`$ and less singular contributions are not essential with our accuracy. Thus, we obtain the following modified BFKL equation at even $`r`$ and $`\gamma \omega `$

$`1=G\left({\displaystyle \frac{1}{\gamma \omega }}+{\displaystyle \frac{1}{\omega ^2}}{\displaystyle \frac{S_1(r)+\omega \widehat{S}_2(r)}{\omega }}\right)+`$
$`G^2\left({\displaystyle \frac{a(r)S_1(r)}{\gamma ^2\omega }}+{\displaystyle \frac{a(r)+b(r)}{\gamma \omega ^2}}+{\displaystyle \frac{S_2(r)+3S_2(r)+\zeta (2)+2d(r)}{2\omega \gamma }}+{\displaystyle \frac{e(r)\widehat{S}_2(r)}{\gamma ^2}}\right).`$
If one will take into account, that $`\omega `$ and $`j+r`$ are generally different parameters (see (13)), it will be needed to substitute in the above equation $`\omega `$ by $`j+r`$ and to change the unity at its left hand side by an analytic function of $`j+r`$

$$11+c_1(r)(j+r)+c_2(r)(j+r)^2+\mathrm{},$$

where $`c_1(r)=\stackrel{~}{C}_1(r)`$, $`c_2(r)=\stackrel{~}{C}_1^2(r)\stackrel{~}{C}_2(r)`$ and the parameters $`\stackrel{~}{C}_1(r)`$ and $`\stackrel{~}{C}_2(r)`$ are defined in Eq. (13).

After this modification of the equation we calculate the behavior of the anomalous dimension at $`jr`$ for even $`r`$ in terms of the parameters $`a(r),b(r),d(r),e(r)`$ and $`c_{1,2}(r)`$

$`\gamma =G\left({\displaystyle \frac{1}{j+r}}+K(r)+L(r)(j+r)\right)+G^2\left({\displaystyle \frac{1}{(j+r)^3}}+{\displaystyle \frac{T(r)}{(j+r)^2}}+{\displaystyle \frac{R(r)}{j+r}}\right),`$ (98)
where

$`K(r)`$ $`=`$ $`a(r)S_1(r)c_1(r),T(r)=b(r)+K(r)c_1(r),`$
$`L(r)`$ $`=`$ $`e(r)\widehat{S}_2(r)\left(a(r)S_1(r)\right)^2+c_1^2(r)c_2(r),`$
$`R(r)`$ $`=`$ $`d(r)+{\displaystyle \frac{S_2(r)+S_2(r)\zeta (2)}{2}}\left(a(r)S_1(r)\right)\left(b(r)+S_1(r)\right)`$
$``$ $`c_1(r)(a(r)+b(r)2S_1(r))+3c_1^2(r)2c_2(r).`$
This expression for $`\gamma `$ should be compared with the result (92) obtained in the previous subsection from the DGLAP equation. The leading poles $`G/(j+r)`$ and $`G^2/(j+r)^3`$ coincide in both cases. As for the regular terms $`G`$ and residues of the poles $`G^2/(j+r)^2`$ and $`G^2/(j+r)`$, the agreement can be achieved with an appropriate choice of the parameters $`a(r),b(r),d(r),e(r)`$ and $`c_{1,2}(r)`$. Thus, to verify the possibility, that the anomalous dimension $`\gamma ^{uni}(j)`$ can be derived completely from the BFKL equation one should calculate these parameters from the Feynman diagrams. As it is seen from the above expressions, the parameters $`c_{1,2}(r)`$ enter together with $`a(r),b(r),d(r),e(r)`$ and could be omitted. The vanishing of $`c_{1,2}(r)`$ corresponds to a natural assumption, that to obtain the universal anomalous dimension of the twist-2 operators one should initially put $`\left|n\right|=r1`$ and only after that $`\omega `$ can be pushed to zero. In this case we obtain $`\gamma ^{uni}(j)`$ (92) if the parameters $`a(r),b(r),d(r)`$ are chosen as follows

$`a(r)`$ $`=`$ $`{\displaystyle \frac{1}{r+1}},b(r)=S_1(r+1),e(r)=S_1^2(r+1){\displaystyle \frac{1}{(r+1)^2}},`$
$`d(r)`$ $`=`$ $`{\displaystyle \frac{S_1(r+1)}{r+1}}+{\displaystyle \frac{S_2(r+1)3S_2(r+1)}{2}}.`$
As a result, the equation for the anomalous dimension can be written as follows

$`1`$ $`=`$ $`G\left({\displaystyle \frac{1}{\gamma (\omega \gamma )}}{\displaystyle \frac{S_1(r)+(\omega \gamma )\widehat{S}_2(r)}{\omega }}\right)+`$
$`+G^2({\displaystyle \frac{S_1(r)}{\gamma ^2\omega }}{\displaystyle \frac{(r+1)^1}{\gamma ^2(\omega \gamma )}}{\displaystyle \frac{S_1(r+1)}{\gamma (\omega \gamma )^2}}+{\displaystyle \frac{S_1^2(r+1)\widehat{S}_2(r+1)}{\gamma ^2}}`$
$`+{\displaystyle \frac{S_2(r)+3S_2(r)}{2\omega \gamma }}+{\displaystyle \frac{\frac{2S_1(r+1)}{r+1}+S_2(r+1)3S_2(r+1)}{2\gamma (\omega \gamma )}}+{\displaystyle \frac{f(r)}{\gamma }}),`$
where the function $`f(r)`$ is fixed from the condition, that its solution reproduces correctly the regular terms $`G^2`$ in $`\gamma ^{uni}(j)`$ at $`jr`$ (see (92)). Providing that we shall take into account also non-singular terms at $`\gamma 0`$, the BFKL equation will contain an important information about the singular part of the anomalous dimension in the order $`G^3`$.

For odd $`r`$ the situation is simpler because here the double-logarithmic terms at $`jr`$ are absent. Indeed, in this case one can write the BFKL equation in the form (after the subsequent substitution $`\omega j+r`$)

$`1+\stackrel{~}{c}_1(r)\omega +\stackrel{~}{c}_2(r)\omega ^2=G({\displaystyle \frac{1}{\gamma \omega }}+{\displaystyle \frac{1}{\omega ^2}}{\displaystyle \frac{S_1(r)+\omega \widehat{S}_2(r)}{\omega }})+G^2({\displaystyle \frac{\stackrel{~}{a}(r)S_1(r)}{\gamma ^2\omega }}`$
$`{\displaystyle \frac{1}{\omega ^3\gamma }}+{\displaystyle \frac{\stackrel{~}{a}(r)+\stackrel{~}{b}(r)}{\gamma \omega ^2}}+{\displaystyle \frac{S_2(r)+3S_2(r)+2\zeta (2)+2\stackrel{~}{d}(r)}{2\omega \gamma }}+{\displaystyle \frac{\stackrel{~}{e}(r)\widehat{S}_2(r)}{\gamma ^2}}),`$
where $`\stackrel{~}{a}(r),\stackrel{~}{b}(r),\stackrel{~}{d}(r),\stackrel{~}{e}(r)`$ and $`\stackrel{~}{c}_{1,2}(r)`$ are the corresponding parameters for odd $`r`$. Therefore the behavior of the anomalous dimension at $`jr`$ for odd $`r`$ is

$`\gamma `$ $`=`$ $`G\left({\displaystyle \frac{1}{j+r}}+\stackrel{~}{K}(r)+\stackrel{~}{L}(r)(j+r)\right)+G^2\left({\displaystyle \frac{\stackrel{~}{T}(r)}{(j+r)^2}}+{\displaystyle \frac{\stackrel{~}{R}(r)}{j+r}}\right),`$ (100)
$`\stackrel{~}{K}(r)`$ $`=`$ $`\stackrel{~}{a}(r)S_1(r)\stackrel{~}{c}_1(r),\stackrel{~}{T}(r)=\stackrel{~}{b}(r)+2\stackrel{~}{K}(r)+\stackrel{~}{c}_1(r),`$
$`\stackrel{~}{L}(r)`$ $`=`$ $`\stackrel{~}{e}(r)\widehat{S}_2(r)(\stackrel{~}{a}(r)S_1(r))^2+\stackrel{~}{c}_1^2(r)\stackrel{~}{c}_2(r),`$
$`\stackrel{~}{R}(r)`$ $`=`$ $`\stackrel{~}{d}(r)+\stackrel{~}{e}(r){\displaystyle \frac{S_2(r)+S_2(r)+2\zeta (2)}{2}}\left(\stackrel{~}{a}(r)S_1(r)\right)\left(3\stackrel{~}{a}(r)+\stackrel{~}{b}(r)2S_1(r)\right)`$
$``$ $`\stackrel{~}{c}_1(r)\left(\stackrel{~}{a}(r)+\stackrel{~}{b}(r)2S_1(r)\right)+2\stackrel{~}{c}_1^2(r)\stackrel{~}{c}_2(r).`$
Again we have a qualitative agreement with the result (92) obtained from the DGLAP equation, but in the order $`G^2`$ the residues of the poles $`1/(j+r)^2`$ and $`1/(j+r)`$ depends on parameters $`\stackrel{~}{a}(r),\stackrel{~}{b}(r),\stackrel{~}{d}(r),\stackrel{~}{e}(r)`$ and $`\stackrel{~}{c}_{1,2}(r)`$. If $`\stackrel{~}{c}_{1,2}(r)=0`$, one can have a complete agreement with Eq. (92) provided that

$`\stackrel{~}{a}(r)`$ $`=`$ $`{\displaystyle \frac{1}{r+1}},\stackrel{~}{b}(r)=2S_1(r+1),\stackrel{~}{e}(r)=S_1^2(r+1){\displaystyle \frac{1}{(r+1)^2}},`$
$`\stackrel{~}{d}(r)`$ $`=`$ $`S_1(r)S_1(r+1)+{\displaystyle \frac{S_2(r)S_2(r)}{2}}+\zeta (2)+2S_2(r+1).`$
In this case the equation for the anomalous dimension can be written in the form

$`1=G\left({\displaystyle \frac{1}{\gamma (\omega \gamma )}}{\displaystyle \frac{S_1(r)+(\omega \gamma )\widehat{S}_2(r)}{\omega }}\right)+`$
$`G^2({\displaystyle \frac{1}{(\omega \gamma )^2\omega \gamma }}{\displaystyle \frac{S_1(r)}{\gamma ^2\omega }}{\displaystyle \frac{(r+1)^1}{\gamma ^2(\omega \gamma )}}+{\displaystyle \frac{2S_1(r+1)}{\gamma \omega ^2}}+{\displaystyle \frac{S_1^2(r+1)\widehat{S}_2(r+1)}{\gamma ^2}}`$ (101)
$`+{\displaystyle \frac{3S_2(r)S_2(r)+2\zeta (2)}{2\omega \gamma }}+{\displaystyle \frac{\zeta (2)S_1(r)S_1(r+1)+\frac{S_2(r)S_2(r)}{2}+2S_2(r+1)}{\gamma (\omega \gamma )}}+{\displaystyle \frac{\stackrel{~}{f}(r)}{\gamma }}),`$
where $`\stackrel{~}{f}(r)`$ is obtained from the condition that the regular terms $`\mathrm{\Delta }\gamma G^2`$ at $`jr`$ coincide with the corresponding contributions to $`\mathrm{\Delta }\gamma ^{uni}`$. Taking into account also the non-singular contributions in the BFKL equation at $`\gamma 0`$ one can obtain some information about singularities of $`\gamma `$ in the order $`G^3`$.

Thus, the results for the universal anomalous dimension $`\gamma ^{uni}(j)`$ obtained from the BFKL equation with the use of our hypothesis, that $`\gamma ^{uni}(j)`$ coincides with $`\gamma ^{BFKL}(j)`$ after its analytic continuation to the singular points $`|n|=r1`$, agree in main features with the expressions derived directly from the DGLAP equation, but for a full agreement one should verify by calculating the corresponding Feynman diagrams, that the parameters $`a(r),b(r),d(r),c_{1,2}(r)`$ and $`\stackrel{~}{a}(r),\stackrel{~}{b}(r),\stackrel{~}{d}(r),\stackrel{~}{c}_{1,2}(r)`$ reproduce correctly the non-leading singularities $`G^2/(j+r)^2`$ and $`G^2/(j+r)`$.

## 7  Conclusion

Above we reviewed the LLA results for the anomalous dimensions of twist-2 operators in the N=4 supersymmetric gauge theory and constructed the operators with a multiplicative renormalization (cf. ). These anomalous dimensions can be obtained from the universal anomalous dimension $`\gamma ^{uni}(j)`$ by a shift of its argument $`jj+k`$ because they belong to the same supermultiplet. We calculated the NLO correction to $`\gamma ^{uni}(j)`$ by using a number of plausible arguments. The results have a compact form (see eqs. (70) and (68)) in terms of polylogarithms and associated special functions. They are verified by direct calculations . Note that recently the LLA anomalous dimensions in this theory were constructed in ref.  in the limit $`j\mathrm{}`$ from the superstring model with the use of the AdS/CFT correspondence .

We investigated properties of the next-to-leading corrections to the kernel of the BFKL equation in the N=4 supersymmetric theory. The absence of the coupling constant renormalization in this model leads presumable to the Möbius invariance of the BFKL equation in higher orders of the perturbation theory. But the holomorphic separability of the BFKL kernel is violated in the NLO approximation (see (B10)). Instead we have the hermitial separability of the corresponding Bethe-Salpeter kernel (see (46)). The cancellation of non-analytic contributions proportional to $`\delta _n^0`$ and $`\delta _n^2`$ in $`N=4`$ SUSY is remarkable (such terms contribute to $`\omega `$ in QCD and in N=1,2 supersymmetric models ). Moreover, the hermitian separability of the BFKL kernel is also valid only in N=4 SUSY.

These properties could be a possible manifestation of the integrability of the reggeon dynamics in the Maldacena model  (see also ) which is reduced to the classical supergravity in the limit $`N_c\mathrm{}`$. Indeed, in this N=4 supersymmetric model the eigenvalues of the LLA pair kernels in the evolution equation for the matrix elements of the quasi-partonic operators are proportional to $`\mathrm{\Psi }(j1)\mathrm{\Psi }(1)`$ , which means, that the corresponding Hamiltonian at large $`N_c`$ coincides with the local Hamiltonian for an integrable Heisenberg spin model . The residues of these eigenvalues at the points $`j=r`$ were obtained in Ref.  from the BFKL equation in LLA by an analytic continuation of the anomalous dimensions to negative integer values of the conformal spin $`|n|`$. It was shown above, that the analogous correspondence between the BFKL and DGLAP equations takes place qualitatively in the NLO approximation (cf. Eq. (92) and Eqs. (98, 100)). To verify quantitatively in this approximation the hypothesis, that in N=4 SUSY the DGLAP equation can be obtained from the BFKL equation, one should calculate the next-to-leading corrections to the BFKL equation in the region, where $`|n|+r+10`$ and $`\gamma \omega `$. We hope to return to this problem in our future publications.

Acknowledgments

The authors are supported in part by the INTAS 00-366 grant. A. Kotikov thanks the Alexander von Humboldt Foundation and Landau-Heisenberg program for a possibility to work on this problem in Germany. He was supported also by the RFBR 02-02-17513 grant. L. Lipatov is thankful to the Hamburg, Montpellier-2 and Paris-6 Universities for their hospitality during the period of time when this work was in progress. He was supported also by the RFBR and NATO grants.

We are indebted to V. Fadin, R. Kirschner, E. Kuraev, R. Peschansky, V. Velizhanin and participants of the PNPI Winter School for helpful discussions. We are grateful also to A. Belitsky, A. Gorsky and G. Korchemsky for stimulation of discussion about DRED scheme and to J. Bartels and M. Lyublinsky for their remarks concerning the double-logarithmic behavior of the Bete-Salpeter amplitude $`A(s,k_{}^2)`$ in QED.

## 8 Appendix A

Here we consider the most complicated contribution $`\mathrm{\Phi }(|n|,\gamma )+\mathrm{\Phi }(|n|,1\gamma )`$ to the NLO correction $`\delta (n,\gamma )`$ in the BFKL equation (21) and derive the property similar to the holomorphic separability (cf. ). To begin with, let us split the function $`\mathrm{\Phi }(n,\gamma )`$ in two terms

$$\mathrm{\Phi }(|n|,\gamma )=\mathrm{\Phi }_1(|n|,\gamma )+\mathrm{\Phi }_2(|n|,\gamma ),$$

where

$$\mathrm{\Phi }_1(|n|,\gamma )=\underset{k=0}{\overset{\mathrm{}}{}}\frac{\left(\beta ^{}(k+|n|+1)(1)^k\mathrm{\Psi }^{}(k+|n|+1)\right)}{k+M}$$

$$+\underset{k=0}{\overset{\mathrm{}}{}}\frac{(1)^k\left(\mathrm{\Psi }(|n|+k+1)\mathrm{\Psi }(1)\right)}{(k+M)^2}$$

and

$$\mathrm{\Phi }_2(|n|,\gamma )=\underset{k=0}{\overset{\mathrm{}}{}}\frac{\left(\beta ^{}(k+1)+(1)^k\mathrm{\Psi }^{}(k+1)\right)}{k+M}$$

$$\underset{k=0}{\overset{\mathrm{}}{}}\frac{(1)^k\left(\mathrm{\Psi }(k+1)\mathrm{\Psi }(1)\right)}{(k+M)^2}\mathrm{\Phi }_2(M).$$

It is convenient also to write the functions $`\mathrm{\Phi }_1(|n|,\gamma )`$ and $`\mathrm{\Phi }_2(|n|,\gamma )`$ as

$`\mathrm{\Phi }_1(|n|,\gamma )`$ $`=`$ $`\mathrm{\Phi }^{(1)}(|n|,M)\mathrm{\Phi }^{(2)}(|n|,M),`$
$`\mathrm{\Phi }_2(|n|,\gamma )`$ $`=`$ $`\mathrm{\Phi }^{(1)}(0,M)+\mathrm{\Phi }^{(2)}(0,M),`$ (A1)
where

$`\mathrm{\Phi }^{(1)}(|n|,s)`$ $`=`$ $`{\displaystyle \underset{k=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\beta ^{}(k+|n|+1)}{k+s}},`$ (A2)
$`\mathrm{\Phi }^{(2)}(|n|,s)`$ $`=`$ $`{\displaystyle \underset{k=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{(1)^k}{k+s}}\left(\mathrm{\Psi }^{}(k+|n|+1){\displaystyle \frac{\mathrm{\Psi }(|n|+k+1)\mathrm{\Psi }(1)}{k+s}}\right).`$ (A3)
If we expand $`\beta ^{}(k+|n|+1)`$ in the series

$$\underset{l=0}{\overset{\mathrm{}}{}}(1)^{l+1}/(l+k+|n|+1)^2$$

and insert it in the sum over $`k`$ in the r.h.s of Eq.(A2), we obtain

$$\mathrm{\Phi }^{(1)}(|n|,s)=\mathrm{\Phi }^{(2)}(|n|,|n|s+1)+\beta ^{}(|n|s+1)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(s)\right].$$

In particular,

$`\mathrm{\Phi }^{(1)}(|n|,M)`$ $`=`$ $`\mathrm{\Phi }^{(2)}(|n|,1\stackrel{~}{M})+\beta ^{}(1\stackrel{~}{M})\left[\mathrm{\Psi }(1)\mathrm{\Psi }(M)\right],`$
$`\mathrm{\Phi }^{(1)}(0,M)`$ $`=`$ $`\mathrm{\Phi }^{(2)}(0,1M)+\beta ^{}(1M)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(M)\right].`$
The function $`\mathrm{\Phi }_2(|n|,\gamma )`$ depends only on $`M`$ and therefore the corresponding contribution to $`\omega `$ has the property of the hermitian separability. Further, for $`\mathrm{\Phi }_1(|n|,\gamma )`$ and $`\mathrm{\Phi }_2(|n|,\gamma )`$ we derive the following relation

$`\mathrm{\Phi }_1(|n|,\gamma )`$ $`+`$ $`\mathrm{\Phi }_1(|n|,1\gamma )`$ (A4)
$`=`$ $`\mathrm{\Phi }^{(1)}(|n|,M)\mathrm{\Phi }^{(2)}(|n|,M)+\mathrm{\Phi }^{(1)}(|n|,1\stackrel{~}{M})\mathrm{\Phi }^{(2)}(|n|,1\stackrel{~}{M})`$
$`=`$ $`\beta ^{}(M)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(1\stackrel{~}{M})\right]+\beta ^{}(1\stackrel{~}{M})\left[\mathrm{\Psi }(1)\mathrm{\Psi }(M)\right],`$
$`\mathrm{\Phi }_2(M)`$ $`=`$ $`\mathrm{\Phi }^{(1)}(0,M)+\mathrm{\Phi }^{(1)}(0,1M)\beta ^{}(M)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(1M)\right]`$ (A5)
$`=`$ $`\mathrm{\Phi }^{(2)}(0,M)+\mathrm{\Phi }^{(2)}(0,1M)+\beta ^{}(1M)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(M)\right].`$
Therefore one can obtain

$`\mathrm{\Phi }(|n|,\gamma )+\mathrm{\Phi }(|n|,1\gamma )=\chi (n,\gamma )\left(\beta ^{}(M)+\beta ^{}(1\stackrel{~}{M})\right)`$
$`+\mathrm{\Phi }_2(M)\beta ^{}(M)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(M)\right]+\mathrm{\Phi }_2(1\stackrel{~}{M})\beta ^{}(1\stackrel{~}{M})\left[\mathrm{\Psi }(1)\mathrm{\Psi }(1\stackrel{~}{M})\right],`$
where $`\chi (n,\gamma )`$ is given by Eq.(20).

## 9 Appendix B

In this Appendix we demonstrate a violation of the generalized holomorhical separability in the NLO correction $`\delta (n,\gamma )`$, i.e. we show an impossibility to present Eqs.(25)-(28) in the separable form symmetric to the substitution $`m\stackrel{~}{m}`$ (i.e., respectively, $`M\stackrel{~}{M}`$ for both positive and negative $`|n|`$).

The expressions appearing in $`\delta (n,\gamma )`$ (21) can be written as follows

$`\mathrm{\Psi }^{\prime \prime }(M)+\mathrm{\Psi }^{\prime \prime }(1\stackrel{~}{M})={\displaystyle \frac{1}{2}}\left[\mathrm{\Psi }^{\prime \prime }(M)+\mathrm{\Psi }^{\prime \prime }(1M)+\mathrm{\Psi }^{\prime \prime }(\stackrel{~}{M})+\mathrm{\Psi }^{\prime \prime }(1\stackrel{~}{M})\right],`$ (B1)
$`\mathrm{\Phi }(|n|,\gamma )+\mathrm{\Phi }(|n|,1\gamma )=\mathrm{\Phi }^{(2)}(0,M)+\mathrm{\Phi }^{(2)}(0,1M)+\mathrm{\Phi }^{(2)}(0,\stackrel{~}{M})+\mathrm{\Phi }^{(2)}(0,1\stackrel{~}{M})`$
$`+\left(\beta ^{}(1M)+\beta ^{}(1\stackrel{~}{M})\right)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(M)\right]+\left(\beta ^{}(M)+\beta ^{}(\stackrel{~}{M})\right)\left[\mathrm{\Psi }(1)\mathrm{\Psi }(1\stackrel{~}{M})\right].`$ (B2)
Using the property

$$\beta ^{}(1M)=\beta ^{}(M)+\pi ^2\frac{\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )},$$

we can present the terms containing the functions $`\beta ^{}`$ in the following form

$`\beta ^{}(1M)+\beta ^{}(1\stackrel{~}{M})=`$
$`{\displaystyle \frac{1}{2}}\left[\beta ^{}(1M)+\beta ^{}(1\stackrel{~}{M})+\beta ^{}(M)+\beta ^{}(\stackrel{~}{M})+\pi ^2{\displaystyle \frac{\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )}}\pi ^2{\displaystyle \frac{\mathrm{cos}((1\stackrel{~}{M})\pi )}{\mathrm{sin}^2((1\stackrel{~}{M})\pi )}}\right],`$
$`\beta ^{}(M)+\beta ^{}(\stackrel{~}{M})=`$
$`{\displaystyle \frac{1}{2}}\left[\beta ^{}(M)+\beta ^{}(\stackrel{~}{M})+\beta ^{}(1M)+\beta ^{}(1\stackrel{~}{M})\pi ^2{\displaystyle \frac{\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )}}+\pi ^2{\displaystyle \frac{\mathrm{cos}((1\stackrel{~}{M})\pi )}{\mathrm{sin}^2((1\stackrel{~}{M})\pi )}}\right].`$
Then, the second line in Eq.(B2) can be replaced by

$`{\displaystyle \frac{1}{2}}\left[\beta ^{}(M)+\beta ^{}(\stackrel{~}{M})+\beta ^{}(1M)+\beta ^{}(1\stackrel{~}{M})\right]\chi (n,\gamma )`$
$`+{\displaystyle \frac{1}{2}}\left({\displaystyle \frac{\pi ^2\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )}}{\displaystyle \frac{\pi ^2\mathrm{cos}((1\stackrel{~}{M})\pi )}{\mathrm{sin}^2((1\stackrel{~}{M})\pi )}}\right)\left[\mathrm{\Psi }(1\stackrel{~}{M})\mathrm{\Psi }(M)\right],`$ (B3)
where $`\chi (n,\gamma )`$ from (19) has the symmetric representation:

$`\chi (n,\gamma )`$ $`=`$ $`2\mathrm{\Psi }(1){\displaystyle \frac{1}{2}}\left[\mathrm{\Psi }(M)+\mathrm{\Psi }(1M)+\mathrm{\Psi }(\stackrel{~}{M})+\mathrm{\Psi }(1\stackrel{~}{M})\right].`$ (B4)
Thus, the NLO term $`\delta (n,\gamma )`$ is presented as follows

$`\delta (n,\gamma )=\delta _1(n,\gamma )+\delta _2(n,\gamma ),`$
where

$`\delta _1(n,\gamma )`$ $`=`$ $`\overline{\varphi }(M)+\overline{\varphi }(\stackrel{~}{M}){\displaystyle \frac{\omega _0}{2\widehat{a}}}\left(\overline{\rho }(M)+\overline{\rho }(\stackrel{~}{M})\right)\left(\omega _0=4\widehat{a}\chi (n,\gamma )\right),`$ (B5)
$`\delta _2(n,\gamma )`$ $`=`$ $`\left({\displaystyle \frac{\pi ^2\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )}}{\displaystyle \frac{\pi ^2\mathrm{cos}((1\stackrel{~}{M})\pi )}{\mathrm{sin}^2((1\stackrel{~}{M})\pi )}}\right)\left[\mathrm{\Psi }(M)\mathrm{\Psi }(1\stackrel{~}{M})\right]`$ (B6)
and

$`2\overline{\rho }(M)`$ $`=`$ $`\beta ^{}(M)+\beta ^{}(1M)+\zeta (2),`$ (B7)
$`2\overline{\varphi }(M)`$ $`=`$ $`6\zeta (3)+\mathrm{\Psi }^{\prime \prime }(M)+\mathrm{\Psi }^{\prime \prime }(1M)2\mathrm{\Phi }^{(2)}(0,M)2\mathrm{\Phi }^{(2)}(0,1M).`$ (B8)
All above terms are symmetric to the substitution $`M(1\stackrel{~}{M})`$ in an agreement with previous results. Moreover, the term $`\delta _1(n,\gamma )`$ is symmetric to $`M\stackrel{~}{M}`$ (and, thus, to $`nn`$) due to the representations (B7) and (B8). Therefore for this contribution the property of generalized holomorphic separability is valid.

The term $`\delta _2(n,\gamma )`$ can be presented as

$`\left({\displaystyle \frac{\pi ^2\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )}}+{\displaystyle \frac{\pi ^2\mathrm{cos}(\stackrel{~}{M}\pi )}{\mathrm{sin}^2(\stackrel{~}{M}\pi )}}\right)\left[\mathrm{\Psi }(M)\mathrm{\Psi }(\stackrel{~}{M}){\displaystyle \frac{\pi \mathrm{cos}(\stackrel{~}{M}\pi )}{\mathrm{sin}(\stackrel{~}{M}\pi )}}\right],`$
where

$`{\displaystyle \frac{\pi \mathrm{cos}(\stackrel{~}{M}\pi )}{\mathrm{sin}(\stackrel{~}{M}\pi )}}={\displaystyle \frac{\pi \mathrm{cos}(M\pi )}{\mathrm{sin}(M\pi )}}.`$
It can be splitted in two parts: $`\delta _2(n,\gamma )=\delta _2^{(1)}(n,\gamma )+\delta _2^{(2)}(n,\gamma )`$ where the contribution

$`\delta _2^{(1)}(n,\gamma )`$ $`=`$ $`\left({\displaystyle \frac{\pi ^2\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )}}+{\displaystyle \frac{\pi ^2\mathrm{cos}(\stackrel{~}{M}\pi )}{\mathrm{sin}^2(\stackrel{~}{M}\pi )}}\right){\displaystyle \frac{\pi \mathrm{cos}(\stackrel{~}{M}\pi )}{\mathrm{sin}(\stackrel{~}{M}\pi )}}`$
$`=`$ $`{\displaystyle \frac{\pi ^3\mathrm{cos}^2(M\pi )}{\mathrm{sin}^3(M\pi )}}{\displaystyle \frac{\pi ^3\mathrm{cos}^2(\stackrel{~}{M}\pi )}{\mathrm{sin}^3(\stackrel{~}{M}\pi )}}`$
is symmetric under the substitution $`M\stackrel{~}{M}`$ (and, thus, to $`nn`$). So, the term $`\delta _2^{(1)}(n,\gamma )`$ has also the property of the generalized holomorphic separability.

Thus, the sum of the contributions $`\delta _1(n,\gamma )`$ and $`\delta _2^{(1)}(n,\gamma )`$ has the generalized separable form symmetric to the substitution $`m\stackrel{~}{m}`$:

$`\delta _1(n,\gamma )+\delta _2^{(1)}(n,\gamma )`$ $`=`$ $`\overline{\varphi }(m)+\overline{\varphi }(\stackrel{~}{m}){\displaystyle \frac{\omega _0}{2\widehat{a}}}\left(\overline{\rho }(m)+\overline{\rho }(\stackrel{~}{m})\right)`$ (B9)
$``$ $`\pi ^3\left({\displaystyle \frac{\mathrm{cos}^2(m\pi )}{\mathrm{sin}^3(m\pi )}}+{\displaystyle \frac{\mathrm{cos}^2(\stackrel{~}{m}\pi )}{\mathrm{sin}^3(\stackrel{~}{m}\pi )}}\right)`$
The last term

$`\delta _2^{(2)}(n,\gamma )=\left({\displaystyle \frac{\pi ^2\mathrm{cos}(M\pi )}{\mathrm{sin}^2(M\pi )}}+{\displaystyle \frac{\pi ^2\mathrm{cos}(\stackrel{~}{M}\pi )}{\mathrm{sin}^2(\stackrel{~}{M}\pi )}}\right)\left[\mathrm{\Psi }(M)\mathrm{\Psi }(\stackrel{~}{M})\right],`$
is anti-symmetric under the transformation $`M\stackrel{~}{M}`$, i.e.

when $`n>0`$, and thus, $`M=m,\stackrel{~}{M}=\stackrel{~}{m}`$ it equals to

$`\left({\displaystyle \frac{\pi ^2\mathrm{cos}(m\pi )}{\mathrm{sin}^2(m\pi )}}+{\displaystyle \frac{\pi ^2\mathrm{cos}(\stackrel{~}{m}\pi )}{\mathrm{sin}^2(\stackrel{~}{m}\pi )}}\right)\left[\mathrm{\Psi }(m)\mathrm{\Psi }(\stackrel{~}{m})\right]`$
and when $`n<0`$, and thus, $`M=\stackrel{~}{m},\stackrel{~}{M}=m`$ it is

$`\left({\displaystyle \frac{\pi ^2\mathrm{cos}(m\pi )}{\mathrm{sin}^2(m\pi )}}+{\displaystyle \frac{\pi ^2\mathrm{cos}(\stackrel{~}{m}\pi )}{\mathrm{sin}^2(\stackrel{~}{m}\pi )}}\right)\left[\mathrm{\Psi }(\stackrel{~}{m})\mathrm{\Psi }(m)\right].`$
Therefore the term $`\delta _2^{(2)}(n,\gamma )`$ in the next-to-leading correction

$`\delta (n,\gamma )`$ $`=`$ $`\overline{\varphi }(m)+\overline{\varphi }(\stackrel{~}{m}){\displaystyle \frac{\omega _0}{2\widehat{a}}}\left(\overline{\rho }(m)+\overline{\rho }(\stackrel{~}{m})\right)`$ (B10)
$``$ $`\pi ^3\left({\displaystyle \frac{\mathrm{cos}^2(m\pi )}{\mathrm{sin}^3(m\pi )}}+{\displaystyle \frac{\mathrm{cos}^2(\stackrel{~}{m}\pi )}{\mathrm{sin}^3(\stackrel{~}{m}\pi )}}\right)+\delta _2^{(2)}(n,\gamma )`$
violates the holomorphic separability. Note, that the anomalous term $`\delta _2^{(2)}(n,\gamma )`$ can be written as Eq. (31), i.e. it is zero for odd $`n`$, where $`\delta (n,\gamma )`$ coincides with its analytic continuation to corresponding negative $`|n|`$.

## 10 Appendix C

Here we construct the anomalous dimension matrices $`\gamma _{ab}^{(0)}(j)`$ $`(a,b=g,q,\phi )`$ and $`\stackrel{~}{\gamma }_{ab}^{(0)}(j)`$ $`(a,b=g,q)`$ in terms of their eigenvalues in LLA

$$\gamma _\pm ^{(0)}(j)=4S_1(j2),\gamma _0^{(0)}(j)=4S_1(j),\stackrel{~}{\gamma }_\pm ^{(0)}(j)=4S_1(j1).$$

To illustrate the procedure we consider the polarized and unpolarized cases separately.

Polarized case. From Eqs. (58), (59) and (61) we have

$`\stackrel{~}{\gamma }_{gg}^{(0)}(j)+{\displaystyle \frac{1}{2}}\stackrel{~}{\gamma }_{qg}^{(0)}(j)=\stackrel{~}{\gamma }_{qq}^{(0)}(j)+2\stackrel{~}{\gamma }_{gq}^{(0)}(j)=\stackrel{~}{\gamma }_+^{(0)}(j)4S_1(j1)`$ (C1)
$`\stackrel{~}{\gamma }_{gg}^{(0)}(j){\displaystyle \frac{j+2}{2(j1)}}\stackrel{~}{\gamma }_{qg}^{(0)}(j)=\stackrel{~}{\gamma }_{qq}^{(0)}(j){\displaystyle \frac{2(j1)}{j+2}}\stackrel{~}{\gamma }_{gq}^{(0)}(j)=\stackrel{~}{\gamma }_{}^{(0)}(j)4S_1(j+1)`$ (C2)
It is obvious, that Eqs.(55)-(59) and (C1)-(C2) correspond to a diagonalization of the anomalous dimension matrix. So we can rewrite Eqs. (C1)-(C2) as

$`\widehat{\stackrel{~}{V}}^1\left(\begin{array}{cc}\stackrel{~}{\gamma }_{gg}^{(0)}(j)& \stackrel{~}{\gamma }_{qg}^{(0)}(j)\\ \stackrel{~}{\gamma }_{gq}^{(0)}(j)& \stackrel{~}{\gamma }_{qq}^{(0)}(j)\end{array}\right)\widehat{\stackrel{~}{V}}=\left(\begin{array}{cc}\stackrel{~}{\gamma }_+^{(0)}(j)& 0\\ 0& \stackrel{~}{\gamma }_{}^{(0)}(j)\end{array}\right)\left(\widehat{\stackrel{~}{V}}^1\widehat{\stackrel{~}{V}}=1\right),`$ (C7)
where

$`\widehat{\stackrel{~}{V}}=\left(\begin{array}{cc}\hfill \stackrel{~}{v}_g& \hfill 2\frac{j1}{j+2}\stackrel{~}{v}_q\\ \hfill \frac{1}{2}\stackrel{~}{v}_g& \hfill \stackrel{~}{v}_q\end{array}\right),\widehat{\stackrel{~}{V}}^1={\displaystyle \frac{j+2}{2j+1}}\left(\begin{array}{cc}\hfill \stackrel{~}{v}_g^1& \hfill 2\frac{j1}{j+2}\stackrel{~}{v}_g^1\\ \hfill \frac{1}{2}\stackrel{~}{v}_q^1& \hfill \stackrel{~}{v}_q^1\end{array}\right)`$ (C12)
with arbitrary values of $`\stackrel{~}{v}_g`$ and $`\stackrel{~}{v}_q`$.

Thus, Eq.(C7) leads to the following representation of the anomalous dimension matrix in the polarized case

$`\widehat{\stackrel{~}{\gamma }}^{(0)}(j)\left(\begin{array}{cc}\stackrel{~}{\gamma }_{gg}^{(0)}(j)& \stackrel{~}{\gamma }_{qg}^{(0)}(j)\\ \stackrel{~}{\gamma }_{gq}^{(0)}(j)& \stackrel{~}{\gamma }_{qq}^{(0)}(j)\end{array}\right)=\widehat{\stackrel{~}{V}}\left(\begin{array}{cc}\stackrel{~}{\gamma }_+^{(0)}(j)& 0\\ 0& \stackrel{~}{\gamma }_{}^{(0)}(j)\end{array}\right)\widehat{\stackrel{~}{V}}^1`$ (C17)
$`={\displaystyle \frac{1}{2(2j+1)}}\left[\stackrel{~}{\gamma }_+^{(0)}(j)\left(\begin{array}{cc}2(j+2)& 4(j1)\\ j+2& 2(j1)\end{array}\right)+\stackrel{~}{\gamma }_{}^{(0)}(j)\left(\begin{array}{cc}2(j1)& 4(j1)\\ (j+2)& 2(j+2)\end{array}\right)\right]`$ (C22)
which does not depend on $`\stackrel{~}{v}_g`$ and $`\stackrel{~}{v}_q`$.

Unpolarized case. From Eqs. (58), (59) and (60) we have

$`\gamma _{gg}^{(0)}(j)+\gamma _{qg}^{(0)}(j)+\gamma _{\phi g}^{(0)}(j)=\gamma _{qq}^{(0)}(j)+\gamma _{gq}^{(0)}(j)+\gamma _{\phi q}^{(0)}(j)`$
$`=\gamma _{\phi \phi }^{(0)}(j)+\gamma _{g\phi }^{(0)}(j)+\gamma _{q\phi }^{(0)}(j)=\gamma _+^{(0)}(j)4S_1(j2),`$ (C23)
$`\gamma _{gg}^{(0)}(j){\displaystyle \frac{1}{2(j1)}}\gamma _{qg}^{(0)}(j){\displaystyle \frac{j+1}{3(j1)}}\gamma _{\phi g}^{(0)}(j)`$
$`=\gamma _{qq}^{(0)}(j)2(j1)\gamma _{gq}^{(0)}(j)+{\displaystyle \frac{2(j+1)}{3}}\gamma _{\phi q}^{(0)}(j)`$
$`=\gamma _{\phi \phi }^{(0)}(j){\displaystyle \frac{3(j1)}{j+1}}\gamma _{g\phi }^{(0)}(j)+{\displaystyle \frac{3}{2(j+1)}}\gamma _{q\phi }^{(0)}(j)=\gamma _0^{(0)}(j)4S_1(j),`$ (C24)
$`\gamma _{gg}^{(0)}(j){\displaystyle \frac{j+2}{(j1)}}\gamma _{qg}^{(0)}(j)+{\displaystyle \frac{(j+1)(j+2)}{j(j1)}}\gamma _{\phi g}^{(0)}(j)`$
$`=\gamma _{qq}^{(0)}(j){\displaystyle \frac{j1}{j+2}}\gamma _{gq}^{(0)}(j){\displaystyle \frac{j+1}{j}}\gamma _{\phi q}^{(0)}(j)`$
$`=\gamma _{\phi \phi }^{(0)}(j)+{\displaystyle \frac{j(j1)}{(j+1)(j+2)}}\gamma _{g\phi }^{(0)}(j){\displaystyle \frac{j}{j+1}}\gamma _{q\phi }^{(0)}(j)=\gamma _{}^{(0)}(j)4S_1(j+2).`$ (C25)
The formulae (C23)-(C25) are equivalent to the matrix equation

$`\widehat{V}^1\left(\begin{array}{ccc}\gamma _{gg}^{(0)}(j)& \gamma _{qg}^{(0)}(j)& \gamma _{\phi g}^{(0)}(j)\\ \gamma _{gq}^{(0)}(j)& \gamma _{qq}^{(0)}(j)& \gamma _{\phi q}^{(0)}(j)\\ \gamma _{g\phi }^{(0)}(j)& \gamma _{q\phi }^{(0)}(j)& \gamma _{\phi \phi }^{(0)}(j)\end{array}\right)\widehat{V}=\left(\begin{array}{ccc}\gamma _+^{(0)}(j)& 0& 0\\ 0& \gamma _0^{(0)}(j)& 0\\ 0& 0& \gamma _{}^{(0)}(j)\end{array}\right)\left(\widehat{V}^1\widehat{V}=1\right),`$ (C32)
where

$`\widehat{V}`$ $`=`$ $`\left(\begin{array}{ccc}v_g& 2(j1)v_q& \frac{j(j1)}{(j+1)(j+2)}v_\phi \\ v_g& v_q& \frac{j}{j+1}v_\phi \\ v_g& \frac{2}{3}(j+1)v_q& v_\phi \end{array}\right),`$ (C36)
$`\widehat{V}^1`$ $`=`$ $`{\displaystyle \frac{(j+1)(j+2)}{2(4j^21)(2j+3)}}\left(\begin{array}{ccc}(2j+3)v_g^1& 4\frac{j1}{j+2}(2j+3)v_g^1& 3\frac{j(j1)}{(j+1)(j+2)}(2j+3)v_g^1\\ 3\frac{2j+1}{j+1}v_q^1& 6\frac{2j+1}{(j+1)(j+2)}v_q^1& \frac{j(2j+1)}{(j+1)(j+2)}v_q^1\\ (2j1)v_\phi ^1& 4(2j1)v_\phi ^1& 3(2j1)v_\phi ^1\end{array}\right)`$ (C40)
with arbitrary values of $`v_g`$, $`v_q`$ and $`v_\phi `$.

Analogously to the polarized case one can write the expression for the anomalous dimension matrix in the form

$`\widehat{\gamma }^{(0)}(j)\left(\begin{array}{ccc}\gamma _{gg}^{(0)}(j)& \gamma _{qg}^{(0)}(j)& \gamma _{\phi g}^{(0)}(j)\\ \gamma _{gq}^{(0)}(j)& \gamma _{qq}^{(0)}(j)& \gamma _{\phi q}^{(0)}(j)\\ \gamma _{g\phi }^{(0)}(j)& \gamma _{q\phi }^{(0)}(j)& \gamma _{\phi \phi }^{(0)}(j)\end{array}\right)=\widehat{V}\left(\begin{array}{ccc}\gamma _+^{(0)}(j)& 0& 0\\ 0& \gamma _0^{(0)}(j)& 0\\ 0& 0& \stackrel{~}{\gamma }_{}^{(0)}(j)\end{array}\right)\widehat{V}^1`$ (C47)
$`={\displaystyle \frac{1}{2(4j^21)(2j+3)}}[(2j+3)\gamma _+^{(0)}(j)\left(\begin{array}{ccc}(j+1)(j+2)& 4(j^21)& 3j(j1)\\ (j+1)(j+2)& 4(j^21)& 3j(j1)\\ (j+1)(j+2)& 4(j^21)& 3j(j1)\end{array}\right)`$ (C51)
$`+(2j+1)\gamma _0^{(0)}(j)\left(\begin{array}{ccc}6(j1)(j+2)& 12(j1)& 6j(j1)\\ 3(j+2)& 6& 3j\\ 2(j+1)(j+2)& 4(j+1)& 2j(j+1)\end{array}\right)`$ (C55)
$`+(2j1)\gamma _{}^{(0)}(j)\left(\begin{array}{ccc}j(j1)& 4j(j1)& 3j(j1)\\ j(j+2)& 4j(j+2)& 3j(j+2)\\ (j+1)(j+2)& 4(j+1)(j+2)& 3(j+1)(j+2)\end{array}\right)]`$ (C59)
which does not depend on $`v_g`$, $`v_q`$ and $`v_\phi `$.

## 11 Appendix D

Here we obtain the expression (16) for the $`t`$-channel partial wave $`f_\omega (\kappa )`$ of the process $`e^+e^{}\mu ^+\mu ^{}`$ in the double-logarithmic approximation. Note that our results coincide with those obtained in the original paper .

The Bethe-Salpeter equation (15)

$$A(s,k_{}^2)=1+\lambda _{m^2}^s\frac{dk_{}^2}{k_{}^2}_{k_{}^2}^{\mathrm{min}(s,sk_{}^2/k_{}^2)}\frac{ds^{}}{s^{}}A(s^{},k_{}^2)\left(\lambda =\frac{\alpha _{em}}{2\pi }\right)$$

can be transformed after going in the $`\omega `$-representation

$$A(s,k_{}^2)=_{\sigma i\mathrm{}}^{\sigma +i\mathrm{}}\frac{d\omega }{2\pi i\omega }\left(\frac{s}{k_{}^2}\right)^\omega f_\omega (\kappa )\left(\kappa =\mathrm{ln}\frac{k_{}^2}{m^2}\right)$$

to the form

$$f_\omega (\kappa )=1+\frac{\lambda }{\omega }\left(_\kappa ^\mathrm{\Lambda }𝑑\kappa ^{}e^{\omega (\kappa ^{}\kappa )}f_\omega (\kappa ^{})+_0^\kappa 𝑑\kappa ^{}f_\omega (\kappa ^{})\right).$$

We introduced here the ultraviolet cut-off $`\mathrm{\Lambda }`$ in $`\mathrm{ln}k_{}^2`$. It will be pushed to infinity in the end of calculations. Note, that in QCD such cut-off could be interpreted as an effect of the interaction vanishing at large $`k^2`$ due to the asymptotic freedom. The differentiation of the integral equation gives

$$e^{\omega \kappa }\frac{d}{d\kappa }f_\omega (\kappa )=\lambda _\kappa ^\mathrm{\Lambda }𝑑\kappa ^{}e^{\omega \kappa ^{}}f_\omega (\kappa ^{}).$$

The second differentiation leads to the differential equation

$$\left(e^{\omega \kappa }\frac{d}{d\kappa }e^{\omega \kappa }\frac{d}{d\kappa }+\lambda \right)f_\omega (\kappa )=0.$$

We can find its solution in terms of two arbitrary constants $`a_\omega ^\pm `$

$$f_\omega (\kappa )=a_\omega ^{}e^{\gamma _\omega ^{}\kappa }+a_\omega ^+e^{\gamma _\omega ^+\kappa },$$

where

$`\gamma _\omega ^\pm ={\displaystyle \frac{\omega }{2}}\left(1\pm \sqrt{14{\displaystyle \frac{\lambda }{\omega ^2}}}\right).`$ (D1)
¿From the above expression for $`e^{\omega \kappa }\frac{d}{d\kappa }f_\omega (\kappa )`$ we obtain, that

$$a_\omega ^{}\gamma _\omega ^{}e^{\gamma _\omega ^{}\mathrm{\Lambda }}+a_\omega ^+\gamma _\omega ^+e^{\gamma _\omega ^+\mathrm{\Lambda }}=0,$$

and from the integral equation

$$1=\frac{\lambda }{\omega }\left(\frac{a_\omega ^{}}{\gamma _\omega ^{}}+\frac{a_\omega ^+}{\gamma _\omega ^+}\right).$$

Therefore

$$a_\omega ^+=\frac{\omega }{\lambda }e^{\mathrm{\Lambda }\gamma _\omega ^{}}\left[\frac{1}{\gamma _\omega ^+}e^{\mathrm{\Lambda }\gamma _\omega ^{}}\frac{\gamma _\omega ^+}{\gamma _\omega ^2}e^{\mathrm{\Lambda }\gamma _\omega ^+}\right]^1,a_\omega ^{}=\frac{\omega }{\lambda }e^{\mathrm{\Lambda }\gamma _\omega ^+}\left[\frac{1}{\gamma _\omega ^{}}e^{\mathrm{\Lambda }\gamma _\omega ^+}\frac{\gamma _\omega ^{}}{\gamma _\omega ^{+2}}e^{\mathrm{\Lambda }\gamma _\omega ^{}}\right]^1.$$

In particular, for $`\mathrm{\Lambda }\mathrm{}`$ we obtain

$$f_\omega (\kappa )=f_\omega e^{\gamma _\omega ^{}\kappa },$$

where the $`t`$-channel partial wave on the mass shell

$`f_\omega ={\displaystyle \frac{\omega }{\lambda }}\gamma _\omega ^{}={\displaystyle \frac{\omega ^2}{2\lambda }}\left(1\sqrt{14{\displaystyle \frac{\lambda }{\omega ^2}}}\right)`$ (D2)
satisfies the infrared evolution equation

$`f_\omega =1+{\displaystyle \frac{\lambda }{\omega ^2}}f_\omega ^2.`$ (D3)
In this case the partial wave can be written also as follows

$`f_\omega (\kappa )={\displaystyle _{\frac{\omega }{2}i\mathrm{}}^{\frac{\omega }{2}+i\mathrm{}}}{\displaystyle \frac{d\gamma }{2\pi i}}\left({\displaystyle \frac{k^2}{m^2}}\right)^\gamma {\displaystyle \frac{(\omega 2\gamma )\gamma \omega ^2\lambda ^2}{\omega \frac{\alpha _{em}}{2\pi }\left[\frac{1}{\omega \gamma }+\frac{1}{\gamma }\right]}}`$ (D4)
Note, that providing that $`\mathrm{\Lambda }=\mathrm{}`$ from the beginning the condition $`a_\omega ^+=0`$ follows from the requirement, that the $`t`$-channel partial wave $`f_\omega (\kappa )`$ does not grow at $`\omega \mathrm{}`$ (in an opposite case we can not reduce the initial Bethe-Salpeter equation for $`A(s,k_{}^2)`$ to the equation for $`f_\omega (\kappa )`$ because the lower limit of integration $`s^{}=k_{}^2`$ is not zero).

Thus, we derived the correct expression (16) for $`f_\omega (\gamma )`$. The ’$``$’ component of the value (D1) coincides with one (17) in the section 2.