arxiv:hep-ph/ v1 28 Apr 1999

Transcription

1 Reconstruction of non-forward evolution kernels. A.V. Belitsky 1, D. Müller Institut für Theoretische Physik, Universität Regensburg D-934 Regensburg, Germany arxiv:hep-ph/994477v1 8 Apr 1999 A. Freund INFN, Sezione di Firenze, Largo E. Fermi 515 Firenze, Italy Abstract We develop a framework for the reconstruction of the non-forward kernels which govern the evolution of twist-two distribution amplitudes and off-forward parton distributions beyond leading order. It is based on the knowledge of the special conformal symmetry breaking part induced by the one-loop anomaly and conformal terms generated by forward next-to-leading order splitting functions, and thus avoids an explicit two-loop calculation. We demonstrate the formalism by applying it to the chiral odd and flavour singlet parity odd sectors. Keywords: evolution equation, two-loop exclusive kernels, conformal constraints PACS numbers: 11.1.Hi, 11.3.Ly, 1.38.Bx 1 Alexander von Humboldt Fellow.

2 1 Introduction The possibility to access new characteristics of hadrons by means of the deeply virtual Compton scattering [1,, 3] and the hard diffractive hadron electroproduction [3, 4] processes has recently initiated a growing phenomenological interest in the underlying non-perturbative elements the so-called off-forward parton distributions OFPD) which parametrize hadronic structure in these reaction making use of the QCD factorization theorems. The main feature of the processes is a non-zero skewedness, i.e. plus component, = η, of the t-channel momentum transfer. One of the central issues which has been addressed in this context is the description of the scaling violation phenomena in the cross section via the evolution of the off-forward parton distributions. Since the OFPD is defined as an expectation value of a non-local string operator, its Q -dependence is governed by the renormalization of this operator Fourier transformed to the momentum fraction space. Inasmuch as the generalized skewed kinematics can be unambiguously restored [5] from the conventional exclusive one, known as Efremov-Radyushkin-Brodsky-Lepage ER-BL) region η = 1, in what follows we deal formally with renormalization of the ordinary distribution amplitudes which obey the ER-BL equation [6, 7] Here we have introduced the exclusive convolution d d lnq φx, Q) = V x, y α sq)) e φy, Q). 1) e 1 dy, to distinguish it from the inclusive one used later. Here φ = Q φ G φ) is the two-dimensional vector and V x, y α s ) is a -matrix of evolution kernels given by a series in the coupling. Several methods have been offered so far to solve the off-forward evolution equation: numerical integration [8], expansion of OFPD w.r.t. an appropriate basis of polynomials [9, 1], mapping to the forward case [1, 13] and solution in the configuration space [14, 15]. The last three methods are based on the well-known fact that operators with definite conformal spin do not mix in the one-loop approximation. Beyond leading order the latter two methods can only be applied in the formal conformal limit of QCD where the β-function is set equal to zero and making use of the conformal subtraction scheme which removes the special conformal symmetry breaking anomaly appearing in the minimal subtraction scheme. Thus, we are only left with the former two methods which allows for a successive improvement of the perturbative approximations involved. Up to now only orthogonal polynomial reconstruction method has allowed the analysis of the scaling violation in the singlet sector in two-loop approximation since only the anomalous This idea has earlier been applied directly to the kernels in [11]. 1

3 dimensions required in the formalism were available so far [11, 16, 17]. This was sufficient to get a first insight into the NLO evolution corrections. However, in order to have an access to the whole kinematical region, especially for small x, η and high precision handling of the x η domain, one should look for a more efficient numerical treatment. This can be achieved with the first method alluded to above. To do the direct numerical integration of the evolution equation one needs the corresponding evolution kernels whose Gegenbauer moments define the anomalous dimensions mentioned earlier. For the time being the former were available at LO order only. The flavour non-singlet ER-BL kernel η = 1) was obtained in NLO by a cumbersome analytical calculation [18, 19, ]. As we have mentioned above the continuation to η [, 1] is a unique procedure [5], so that one can obtain in a simple way the evolution kernels for OFPD. The goal of this paper is to outline a method that allows one to construct the singlet ER-BL kernels by applying conformal and supersymmetric constraints where the latter ones arise from the N = 1 super Yang-Mills theory [1, ]. In this way we can avoid the direct diagrammatical calculation which would be very difficult to handle otherwise since no appropriate technology has been developed yet. The derivation is based on the fairly well established structure of the ER-BL kernel in NLO. Up to two-loop order we have V x, y α s ) = α s π V ) x, y) ) αs V 1) x, y) Oαs 3 ), ) π with the purely diagonal LO kernel V ) in the basis of Gegenbauer polynomials and NLO one separated in two parts: V 1) x, y) = V D1) x, y) V ND1) x, y), with the diagonal part which is entirely determined by the well-known forward DGLAP splitting functions Pz) [11] AB V D x, y) = 1 dz wy ν) j= N j ν) CνA) j x 1)z j AB Pz)C νb) j y 1), 3) where N j ν) = 4ν1 Γ 1 )Γνj) Γ ν)νj)j! and wy ν) = yȳ) ν 1/ are the normalization and weight factors, respectively. The non-diagonal piece is fixed completely by the conformal constraints [17] { V ND1) x, y) = I D) V e V ) β ) 1l [g e ] } V ) x, y), 4), in terms of V ) = C F QQ V ) T F N QG f V ) C GQ F V ) C GG A V ), g = C F QQ g C F GQ g C A GG g, 5) the ER-BL kernels at LO and the special conformal symmetry breaking matrix g. Here β = 4 T 3 FN f 11C 3 A is the first expansion coefficient of the QCD β-function. In the parity odd sector the dotted kernel, V, is simply given by a logarithmic modification of the V ). Due to subtleties,

4 appearing in the parity even case [11], we deal here, for the sake of simplicity, only with the parity odd and transversity sectors. The main problem is thus to restore the diagonal part of the NLO kernels. Since the use of Eq. 3) beyond LO is extremely complicated in practice, we are forced to look for other solutions. It turns out that the bulk of contributions in the ER-BL kernel can be deduced by going to the forward limit making use of the reduction Then the difference 3 1 Pz) = LIMV x, y) lim τ τ τ z QQ V GQ V 1 z 1 QG V τ GG V Pz) P cross ladder z) LIM V ND x, y) ext z τ, 1 τ ). 6) can be represented in terms of inclusive convolutions of simple splitting functions and the back transformation to the exclusive kinematics is trivial. The contributions of the purely diagonal cross-ladder diagrams V cross ladder z) can be found from the known QQ sector [18, 19, ] exploiting the N = 1 supersymmetric constraints. The paper is organized as follows. In the next section we analyze the structure of the known flavour non-singlet ER-BL kernel and state the benchmarks of the formalism. The structure observed will give us a guideline for construction of all other kernels: quark chiral odd sector is considered in Section 3 and parity odd flavour singlet one is discussed in Section 4. Finally, we give our conclusions and an outlook. Structure of ER-BL kernel in non-singlet sector. It is very instructive to demonstrate the machinery in the simplest case of non-singlet sector. Since the explicit two-loop calculation is available [18, 19, ] the direct comparison can be made. The NLO QQ-kernel can be decomposed in colour structures as 4 V x, y α s ) = α s π C FV ) x, y) ) [ αs CF C F V F x, y) β π V βx, y) O α 3 s), C F C A ) ] V G x, y) 7) 3 Here V ND is understood without the I D)-projector. 4 We omit the superscript QQ later in this section. 3

5 with the LO kernel V ) x, y) = [vx, y)], where vx, y) = θy x)fx, y) θx y)fx, y), and fx, y) = x 1 1 ). 8) y y x The shorthand notations x = 1 x and f = f x, ȳ) are used throughout the paper. The - prescription is conventionally defined by 1 [V x, y)] = V x, y) δx y) dz V z, y). Let us now recall a few properties of the kernel that are useful for the following considerations. Due to absence of the conformal symmetry breaking counterterms at leading order for the renormalization of the composite operators with total derivatives, one can use its consequences to fix the eigenfunctions which turn out to be the Gegenbauer polynomials C 3/ j x 1) [6, 7]. Thus, the LO kernel is symmetric with respect to the weight function x x: yȳv ) x, y) = x xv ) y, x). Its eigenvalues are given by the anomalous dimensions appeared in the analysis of deep inelastic scattering. Thus, it is not surprising that a simple limit already mentioned in Eq. 6) gives us the DGLAP kernel [5]: 1 z Pz) = LIMV x, y) lim τ τ V ext τ, 1 ). 9) τ To perform this limit, we have to extend at first the ER-BL kernel, originally defined in the domain x, y 1, to the whole region x, y, ) by a unique procedure which is given in practice by the replacement, e.g. at leading order, of the θ-function by θy x) θ 1 x ) ) x θ signy). 1) y y If a kernel is diagonal in the ER-BL representation, we can restore it from the known DGLAP kernel by the integral transformation 3). Because of branch cuts appearing in the convolutions of the NLO terms with the transformation kernel, it is highly nontrivial to handle the inverse reduction to the exclusive kinematics. At NLO the kernel 7) contains besides a pure diagonal part with respect to the Gegenbauer polynomials also a non-diagonal part located in V F x, y) and V β x, y). These parts are predicted by conformal constraints see QQ-entry of Eq. 4)) and are fixed by the one-loop special conformal anomaly kernels [11, 16, 17]: vx, y) = θy x)fx, y) ln x { } x x y y ȳ, gx, y) = θy x) ln 1 x y y x ) { } x x. 11) y ȳ Let us now analyze in detail the contributions to the NLO kernel from different colour structures. The expressions for C F terms arise from Feynman diagrams containing quark self-energy 4

6 insertions and ladder graphs 5. In order to subtract the ultraviolet UV) divergences in subgraphs it requires the LO renormalization of the composite operator to which these lines are attached to. The explicit calculation gives [18, 19, ] V F x, y) = θy x) f x ȳ { 4 )f 3 ζ) 3 x 3 y f x ) ln x ȳ y f f) ln x y ln 1 x ) y ) } ln x x { } x x y ȳ ln x 1 ln x ln x). 1) y ȳ Making use of the known non-diagonal part 4), V F can be represented up to a pure diagonal term, denoted as D F x, y), by the convolution V F x, y) = v v e g e v v g e ) x, y) D F x, y). 13) To find an appropriate representation of this missing diagonal element we first take the forward limit. Since the forward limit of the convolution is 6 } e LIM {[ v] [v] = {LIM [ v] } i {LIM [v] }, 14) where we have introduced the inclusive convolution P 1 z) i P z) 1 1 dx dyδz xy)p 1 x)p y), and the commutator g V ) V ) g drops out in the forward limit, we obtain LIM V F x, y) = ṗ i p LIM D F x, y), 15) where ṗ = LIM v = pz) ln z 1 z and pz) = LIMvx, y) = 1 z )/1 z). The comparison of LIM V F x, y) with the corresponding part of the DGLAP kernel [3] P F z) = { 4 3 ζ) 3 ln z ln z lnz ln1 z)} pz) 1 z 1 3z ln z 1 z ln z, 16) yields the result in which all double log terms are contained in the convolution ṗ p and, therefore, only single logs survive in D F z) = LIMD F x, y): D F z) = P F ṗ i p = 1 pa z) ln z p a z) 5 For simplicity we imply the diagrams in the light-cone gauge [18, 19]. 6 We remind as well that [A] [B] = [C]. { lnz ln1 z) 1 } 5 pz). 17) 1 5

7 Here we have introduced for convenience the kernel p a z) = 1 z. The next important point is that the remaining log terms can be represented as convolutions of p a and p. Thus, we have finally D F z) = 1 pa i { p p a }z) 1 1 pz) 5 pa z). 18) Since D F x, y) is by definition diagonal, the extension of D F z) towards the ER-BL kinematics is trivial: D F z) D F x, y) = 1 va e v v a ) x, y) 1 1 vx, y) 5 va x, y), 19) where a new diagonal element is v a x, y) = θy x) x θx y) x. Evaluating the convolutions y ȳ one can establish the equivalence of our prediction with Eq. 1). Next, the Feynman diagrams containing vertex and self-energy corrections provide V β proportional to β. Its off-diagonal part is induced by the renormalization of the coupling and is contained in the dotted kernel 11) V β x, y) = vx, y) D β x, y). ) The remaining diagonal piece, D β, is deduced from the known NLO DGLAP kernel [3] P β z) = 5 3 pz) pa z) ṗz) 1) by going to the forward kinematics and restoring then the missed contributions from it. Thus, Indeed, the final result coincides with [18, 19, ]. D β x, y) = 5 3 vx, y) va x, y). ) Finally, we come to the contribution which mainly originates from the crossed ladder diagram proportional to C F C A /): V G x, y) = v a x, y) 4 3 vx, y) Gx, y) { }) x x. 3) y ȳ Since this diagram has no UV divergent subgraph and thus requires no subtraction, its contribution has to be diagonal w.r.t. the Gegenbauer polynomials. This is obvious for the first two terms appearing in Eq. 3). The function 7 Gx, y) contains in the unusual θy x)-structure the mixing between quarks and antiquarks Gx, y) = θy x)hx, y) θy x)hx, y), 4) 7 We have slightly changed the original definition given in [18, 19] by Gx, y) θy x)f lny ln x Gx, y). 6

8 with Hx, y) = [ f Li x) ln y ln x) f Li ȳ) ], 5) [ Hx, y) = f f) Li 1 x ) 1 ) ] y ln y f Li ȳ) Li x) ln y ln x), 6) where Li is the dilogarithm. It can be easily checked that the G-contribution not the terms H and H separately) is symmetrical w.r.t. the weight x x. Performing the limit 9) we obtain the following correspondence with the non-singlet DGLAP kernel [1]: where Gz) LIMGx, y) = θz)θ1 z)hz) θ z)θ1 z)hz), 7) Hz) LIM Hx, y) = pz) ln z ζ) ) Tz), 8) Hz) LIM Hx, y) = pz)s z) T z). 9) Here S z) = 1/1z) dx 1 x z/1z) ln and Tz) = 1 z) ln z 41 z). We should emphasize that x x this structure of G is the most general, especially, it is the only contribution that contains Spence functions. In the forward limit we obtain therefore a typical combinations given in 8) and 9), which can be found in all other channels as well. This observation provides us with a hint for the construction of all singlet G kernels in the ER-BL representation. For completeness, we give the corresponding part of the DGLAP kernel [3] and Gz) defined above in Eqs. 7)-9). P G z) = LIMV G x, y) = p a z) 4 pz) Gz), 3) 3 Recapitulating the results obtained in this section, we have observed a rather simple structure of the non-singlet NLO kernel in the QQ-channel. Up to the diagonal G-function, which is in fact the only new element in the two-loop approximation, we can represent all other terms by a simple convolution of LO kernels already known. It is not accidental but a mere consequence of the topology of contributing Feynman graphs at Oαs ). Thus, we anticipate the same feature to appear in all other channels as well. 3 Quark kernel in chiral odd sector. After we have outlined and tested in the preceding section our formalism, we can apply it to the previously unknown transversity two-loop ER-BL kernel. We decompose the transversity kernel 7

9 analogous to the chiral even case 7). We also use the same decomposition for the DGLAP kernels [4]. The leading order kernel is with V )T x, y) = [ v b x, y) ] 1 δx y), 31) v b x, y) = θy x)f b x, y) θx y)f b x, y), and f b x, y) = x y 1 y x. 3) The non-diagonal part has been analyzed in Ref. [17] and is completely analogous to the chiral even case discussed above. Thus, { VF T x, y) = [ v ] b e V )T e [g] V )T V )T e } [g] x, y) DFx, T y), 33) where v b is obtained from Eq. 11) by replacing f by f b and the g kernel is the same as in Eq. 11). Taking the forward limit of VF T x, y) and comparing it with the known result for the DGLAP kernel in Ref. [4], we find the following trivial representation of the remaining diagonal part D T F x, y) = 3 [ v b x, y) ] 19 δx y). 34) 4 There is essentially no extra work required to find the contribution proportional to the β - function, since it can be easily traced from the DGLAP kernel [4] to be V T β x, y) = 5 3 [ v b x, y) ] [ v b x, y) ] 13 δx y). 35) 1 The case of the G function is easy to handle as well. If we replace f by f b in the definition 4), we obtain the diagonal G T x, y) kernel. Taking the forward limit and comparing it with the DGLAP kernel, we immediately find the remaining δ-function contribution, so that the whole result reads V T G x, y) = [ G T x, y) { }] x x y ȳ This completes the discussion of the quark chiral-odd channel. 19 δx y). 36) 6 4 Flavour singlet parity odd sector. Let us now address the flavour singlet parity odd sector responsible for the evolution of axial-vector distribution amplitudes. For even parity there are few subtleties, which will be discussed elsewhere. Here we would only like to note that in the latter case a direct leading order calculation provides a result that suffers for the mixed channel from off-diagonal matrix elements in the unphysical 8

10 sector. Although the improved result has been found in Ref. [11], it still remains a difficult task to find an appropriate representation for the dotted kernels and the two-loop G functions. Making use of the known non-diagonal part of the ER-BL kernel 4), the whole NLO result in the axial-vector case reads V 1)A = V A e V )A β ) 1l g e V )A V )A e g D A G A, 37) where the kernels D A x, y) and G A x, y) are purely diagonal. Here the matrix of the LO kernels is given in a compact form by [ V )A C QQ F vx, y) ] x, y) = T FN QG f v a x, y) [ C GQ F v a x, y) C GG A v A x, y) ] β δx y), 38) where QQ v QQ v a QQ v b and GG v A GG v a GG v b. The general structure of the functions v i is { } { x x A = B AB v i x, y) = θy x) AB f i x, y) ± for y ȳ A B, 39) with { AB f a } AB f b = xνa) 1/ y νb) 1/ ) y x The index νa) coincides with the index of Gegenbauer polynomials in the corresponding channel, i.e. νq) = 3/ and νg) = 5/. The dotted kernels involved in the definition 37) can simply be obtained by differentiating LO results w.r.t. the index ν which gives rise to the additional lnx/y)-multiplier in front of the former [ V )A C QQ F vx, y) ] x, y) = T QG FN f v a x, y) GQ C F v a [ x, y) C GG A v A x, y) ], 41) with the matrix elements AB vx, y) = θy x) AB fx, y) ln x y ± { x x y ȳ }, for { A = B A B. 4) Note that for A = B the dotted kernels are defined with the -prescription. The g function is given by gx, y) = θy x) [ ] ln1 C x y) F y x [ ] x ln1 C F C x y) y A y x { } x x ±, 43) y ȳ with ) sign corresponding to non-) diagonal elements. Note, that we have used the property I D) ln1 x y ) = I D)x y of Ref. [17]. for the element of GQ-channel to make contact with the results 9

11 Next we construct the diagonal Gx, y) kernel. At first glance one would naively expect that one can obtain these kernels by only inserting appropriate AB f functions in the definition 4), so that the symmetry properties of the f functions w.r.t. the weight induce then the desired symmetry of the AB G functions. Unfortunately, the symmetry is not sufficient for the diagonal form of the Gx, y) kernel. To ensure the diagonality, we have to add terms containing single logs and rational functions. Let us define the matrix ) [ G A x, y) = 1 C F CF C A QQ G A x, y) ] C AT F N QG f G A x, y) C F C GQ A G A [ x, y) CA GG G A x, y) ], 44) with the following general structure AB G A x, y) = θy x) AB H A AB H A) x, y) θy x) AB H A AB H A) x, y). 45) Here analogous to the non-singlet case we set AB H A x, y) = [ ± AB f A Li x) ln y ln x) AB f A Li ȳ) ], 46) [ ) ) ] AB f A AB f A) Li AB H A x, y) = 1 x y 1 ln y AB f A Li ȳ) Li x) ln y ln x) where the upper lower) sign corresponds to the A = B A B) channels. An explicit use of the reduction P V D procedure 3) to restore the H contributions is rather involved due to complexity of the integrand function. Rather we have succeeded to deduce them using different arguments. Since the crossed ladder diagrams have no UV divergent subgraphs the kernels AB G in different channels are related in a scheme independent way by supersymmetry and conformal covariance of N = 1 super Yang-Mills theory []. Employing these symmetries we restore 8 all necessary terms in a straightforward manner to be QQ H A x, y) = QQ H A x, y) =, 48) QG H A x, y) = x yȳ ln x x ln y, yȳ QG H A x, y) = x yȳ ln x x ln y, yȳ 49) GQ H A x, y) = x x y ln x x x ȳ lny, QG H A x, y) = x x y ln x x x ln y, ȳ 5) GG H A x, y) = x 1 x y) x x y y ȳ ȳ ln x y x x x) yȳ ln x GG H A x, y) = x y x x x 1 x x ȳ yȳ y ln x ȳ) x x ln x x yȳ y x x x) yȳ x x 1 yȳ, 47) ln y, 51) ln x 6 x x ln y. yȳ Finally, we have to extract the remaining diagonal piece D A of V A in the forward limit 6) from the known DGLAP kernel P A [5]. We take into account the underlying symmetry of the 8 The details will be presented elsewhere. 1

12 singlet parton distributions to map the antiparticle contribution, i.e. z <, into the region z >. As expected we find from D A z) = P A z) LIM { V e V )A β ) } 1l g V e )A V )A g e G A a simple convolution-type representation for ER-BL kernels which can be immediately deduced from the forward results for singlet QQ-channel QQ D A = CF [D β F] C F [D β] C F C F C A 5) ) [ 4 QQ v QQ v a] 6 C F T F N QQ f v a, 53) 3 where D F, D β are given by Eqs. 19) and ), respectively. The rest of channels is expressed as { QG D A = 3 C F T F N QQ f v a e QG v a 1QG v a} 54) GQ D A GG D A C A T F N f {3 QQ v a e QG v a [1 ζ)] QG v a}, { [GG = CF v A] e GQ v a 3GQ v a} β C F { 3 [ C F C GG A v A] [ e GQ v c [ GG v A] GG 1 v a] e GQ v a {[ [GG = CA v A] GG 1 v a] e GG v a 4 3 { β C A 1GG v a e GG v a 5 [ GG v A] 3 GG v a δx y)} {GG } C F T F N f v a e GG v a GG v a δx y), { GQ v a e [ QQ v ] 1GQ v a}, 55) 6 [ GG v A] 1GG v a δx y) 4 [ 7 3 ζ) ] GQ v a}, } 56) where we have introduced a new kernel GQ v c x, y) = θy x) x y xy ȳ) These results provide us with the explicit parity odd singlet evolution kernels. { } x x. 57) y ȳ 5 Conclusions. In this paper, we have presented a simple method for construction of the exclusive evolution kernels in NLO from the knowledge of the conformal anomalies and the available two-loop splitting functions. The main task was, of course, the reconstruction of the diagonal part of the kernel in the basis of Gegenbauer polynomials. In the course of study we have established convolution-type formulae for the bulk of contributing two-loop graphs with an exception of cross-ladder diagrams. The complications which arise in the restoration of the latter from the known forward kernels has been overcome making use of 11

13 N = 1 supersymmetric constraints []. The former feature suggests that by disentangling the topology of corresponding diagrams, it might allow for an effective and facilitated way of explicit calculation. One may expect that this property persists for a subset of diagrams at higher orders and can be used, e.g. for diagrammatical derivation of NNLO splitting functions. The details of the present formalism together with the flavour singlet parity even case, where new subtleties appear, will be discussed elsewhere. A.B. was supported by the Alexander von Humboldt Foundation. References [1] D. Müller, D. Robaschik, B. Geyer, F.-M. Dittes, J. Hořejši, Fortschr. Phys ) 11. [] X. Ji, Phys. Rev. D ) 7114; J. Phys. G ) 1181; [3] A.V. Radyushkin, Phys. Rev. D ) 554. [4] J.C. Collins, L. Frankfurt, M. Strikman, Phys. Rev. D ) 98. [5] F.-M. Dittes, B. Geyer, D. Müller, D. Robaschik, J. Hořejši, Phys. Lett. B ) 35; [6] A.V. Efremov, A.V. Radyushkin, Theor. Math. Phys ) 97; Phys. Lett. B ) 45. [7] S.J. Brodsky, G.P. Lepage, Phys. Lett. B ) 359; Phys. Rev. D 198) 157. [8] L. Frankfurt, A. Freund, V. Guzey, M. Strikman, Phys. Lett. B ) 345; Phys. Lett. B ) 414 E); K.J. Golec-Biernat, A.D. Martin, Phys. Rev. D ) 149. [9] A.V. Belitsky, B. Geyer, D. Müller, A. Schäfer, Phys. Lett. B ) 31; A.V. Belitsky, D. Müller, L. Niedermeier, A. Schäfer, Phys. Lett. B ) 16; Nucl. Phys. B ) 79. [1] L. Mankiewicz, G. Piller, T. Weigl, Eur. J. Phys. C ) 119. [11] A.V. Belitsky, D. Müller, Nucl. Phys. B ) 7. [1] A.G. Shuvaev, Solution of the off forward leading logarithmic evolution equation based on the Gegenbauer moments inversion, hep-ph/

14 [13] A.G. Shuvaev, K.J. Golec-Biernat, A.D. Martin, M.G. Ryskin, Off diagonal distributions fixed by diagonal partons at small x and ξ, hep-ph/9941. [14] I.I. Balitsky, V.M. Braun, Nucl. Phys. B ) 541. [15] N. Kivel, L. Mankiewicz, Conformal string operators and evolution of skewed parton distributions, hep-ph/ [16] D. Müller, Phys. Rev. D ) 55. [17] A.V. Belitsky, D. Müller, Nucl. Phys. B ) 397. [18] M.H. Sarmadi, Phys. Lett. B ) 471. [19] F.-M. Dittes, A.V. Radyushkin, Phys. Lett. B ) 359. [] S.V. Mikhailov, A.V. Radyushkin, Nucl. Phys. B ) 89. [1] A.P. Bukhvostov, G.V. Frolov, E.A. Kuraev, L.N. Lipatov, Nucl. Phys. B ) 61. [] A.V. Belitsky, D. Müller, A. Schäfer, Phys. Lett. B ) 16. [3] G. Curci, W. Furmanski, R. Petronzio, Nucl. Phys. B ) 7. [4] W. Vogelsang, Phys. Rev. D ) [5] R. Mertig, W.L. van Neerven, Z. Phys. C ) 637; W. Vogelsang, Phys. Rev. D ) 3; Nucl. Phys. B )