arxiv: v2 [hep-th] 26 Feb 2008

Transcription

1 CERN-PH-TH/ DESY BFKL Pomeron, Reggeized gluons and Bern-Dixon-Smirnov amplitudes J. Bartels, L. N. Lipatov,, A. Sabio Vera 3 arxiv: v [hep-th] 6 Feb 008 II. Institut Theoretical Physics, Hamburg University, Germany St. Petersburg Nuclear Physics Institute, Russia 3 CERN, Geneva, Switzerland March 6, 009 Abstract After a brief review of the BFKL approach to Regge processes in QCD and in supersymmetric SUSY gauge theories we propose a strategy for calculating the nextto-next-to-leading order corrections to the BFKL kernel. They can be obtained in terms of various cross-sections for Reggeized gluon interactions. The corresponding amplitudes can be calculated in the framework of the effective action for high energy scattering. In the case of N = 4 SUSY it is also possible to use the Bern-Dixon-Smirnov BDS ansatz. For this purpose the analytic properties of the BDS amplitudes at high energies are investigated, in order to verify their self-consistency. It is found that, for the number of external particles being larger than five, these amplitudes, beyond one loop, are not in agreement with the BFKL approach which predicts the existence of Regge cuts in some physical channels. Introduction The elastic scattering amplitude in QCD at high energies for particles with color indices A,B and helicities λ A,λ B in the leading logarithmic approximation LLA has the Regge form [] A = gδ λa λ A TAA c s +ωt g TBB c t δ λ B λ B, t = q. The gluon Regge trajectory, jt = + ωt, reads ω q = α sn c π πµǫ d ǫ k q k q k a ln q µ, ǫ where we have introduced dimensional regularization with D = 4 ǫ and the renormalization point µ for the t Hooft coupling constant a = α s N c π 4πe γ ǫ. 3

2 Figure : Multi-Regge kinematics The gluon trajectory is also known in the next-to-leading approximation in QCD [] and in SUSY gauge models [3]. In LLA gluons are produced in the multi-regge kinematics see Fig.. In this kinematics the gluon production amplitude in LLA has the factorized form A +n = sg δ λa λ A T c where s ω q AA q gc µ q,q e µk T d s ω q c c q... sω q n+ n+ q n+ g δ λb λ B T c n+ BB,4 s = p A + p B s r = k r + k r q r, k r = q r+ q r. 5 The matrices T a bc are the generators of the SUN c gauge group in the adjoint representation and C µ q,q are the effective Reggeon-Reggeon-gluon vertices. In the case when the polarization vector e µ k corresponds to a produced gluon with a definite helicity one can obtain [4] C C µ q,q e µk = q q k, 6 where the complex notation q = q x + iq y for the two-dimensional transverse vectors has been used. The elastic scattering amplitude with the vacuum quantum numbers in the t- channel can be calculated with the use of s-channel unitarity []. In this approach the Pomeron appears as a composite state of two Reggeized gluons. It is also convenient to use transverse coordinates in a complex form together with their canonically conjugated momenta as ρ k = x k + iy k, ρ k = x k iy k, p k = i, p k = i ρ k ρ. 7 k In the coordinate representation the Balitsky-Fadin-Kuraev-Lipatov BFKL equation for the Pomeron wave function can then be written as follows [] E Ψ ρ, ρ = H Ψ ρ, ρ, = α sn c π min E, 8

3 where is the Pomeron intercept. The BFKL Hamiltonian has the rather simple representation [5] H = ln p p + p p ln ρ p p + p p ln ρ p p 4ψ, 9 with ρ = ρ ρ. The kinetic energy is proportional to the sum of two gluon Regge trajectories ω p i i =,. The potential energy ln ρ is related to the product of two gluon production vertices C µ. This Hamiltonian is invariant under the Möbius transformation [6] ρ k aρ k + b cρ k + d, 0 where a, b, c and d are complex numbers. The eigenvalues of the corresponding Casimir operators are expressed in terms of the conformal weights m = + iν + n, m = + iν n for the unitary principal series representation of SL,C, with ν being real and n integer. The BFKL Hamiltonian can be iterated in the s channel to account for the exchange of an arbitrary number of Reggeized gluons. In the N c limit this iteration is described by the Bartels-Kwiecinski-Praszalowicz BKP equation [7] for the n-gluon colorless composite state. In this case the Hamiltonian has the property of holomorphic separability [8] in the form H = H k,k+ = h + h, [h,h ] = 0. k The holomorphic Hamiltonian can be written as h = k h k,k+, h = lnp p + p ln ρ p + p ln ρ p ψ, 3 where ψx = ln Γx. Consequently, the wave function Ψ fulfills holomorphic factorization [8] and there exists the remarkable duality symmetry [9] p i ρ i,i+ p i+. 4 Moreover, in the holomorphic and anti-holomorphic sectors, there are integrals of motion commuting among themselves and with h [5, 0]: q r = k <k <...<k r ρ k k ρ k k 3...ρ krk p k p k...p kr, [q r,h] = 0. 5 The integrability of BFKL dynamics was proven in [0] and is related to the fact that h coincides with the local Hamiltonian of the Heisenberg spin model []. In the LLA the Pomeron intercept is = 4 αs π N c ln > 0 and the Froissart bound σ t < cln s for the total cross-section σ t s is violated []. In the next-to-leading logarithmic approximation the integral kernel for the BFKL equation was constructed in Refs. [3, ]. Due to its Möbius invariance, the solution of the BFKL equation can be classified by the anomalous dimension γ = +iν of twist- operators and the conformal spin n, which coincides with the number of transverse indices of the local operators O µ...µ j. 3

4 The eigenvalue of the BFKL kernel in the next-to-leading approximation has the form [] ω = 4â [ ψ ψ γ + n ψ γ + n ] + 4 â n,γ, â = g N c 6π. 6 In QCD, the next-to-leading contribution n, γ is a non-analytic function of the conformal spin n because it contains some terms dependent on the Kronecker symbols δ n,0 and δ n,. However, in N = 4 SUSY this dependence is canceled and we obtain the following hermitially separable expression [3] where and ΦM = n,γ = φm + φm ρm + ρm â/ω ρm = β M + ζ, β z = [ Ψ z + 4, M = γ + n, 7 Ψ z ]. 8 φm = 3ζ3 + ψ M ΦM + β M ψ ψm, 9 k=0 β k + k + M + k ψ k + k + M k=0 ψk + ψ. 0 k + M Very importantly, all these contributions have the property of maximal transcendentality [3]. The behaviour of the gluon Green s function corresponding to this NLO kernel in N = 4 SUSY was investigated in [4]. The NLO conformal spins affect azimuthal angle decorrelations in jet physics as it was originally suggested in [5]. In a different context, the one-loop anomalous dimension matrix for twist- operators in N = 4 SUSY can be easily calculated since it is completely fixed by superconformal invariance. Its eigenvalue is proportional to ψ ψj, which is related to the integrability of the evolution equation for the quasi-partonic operators in this model [6]. The integrability of N = 4 SUSY has also been established for other operators and in higher loops [7, 8]. The maximal transcendentality principle suggested in Ref. [3] made it possible to extract the universal anomalous dimension up to three loops in N = 4 SUSY [9, 0] from the QCD results []. This principle was also helpful to find closed integral equations for the cusp anomalous dimension in this model [, 3] based on the AdS/CFT correspondence [4, 5, 6]. In the framework of the asymptotic Bethe ansatz the maximal transcendentality principle helped to fix the anomalous dimension at four loops [7]. However, the obtained results contradict the predictions stemming from the BFKL equation [3, 3]. The origin of this discrepancy is related to the onset of wrapping effects [7]. In this framework it is, therefore, crucial to obtain more information from the BFKL side through the calculation of its higher order corrections to the integral kernel. We would like to point out that the intercept of the BFKL Pomeron at large t Hooft coupling constant in N = 4 SUSY was found in Refs. [0, 8]. In the present paper we want to formulate a program to calculate the three loop corrections to the BFKL kernel in the t Hooft coupling. Our calculational approach is based on the use of the high energy effective action developed in [9, 30] for the construction of the different Reggeized gluon couplings, and on the BDS ansatz [3] for the scattering amplitudes in the N = 4 super Yang-Mills theory. We begin with a short review of the effective action section and then turn to an analysis of the BDS formula in the Regge limit for amplitudes up to six external gluons section 3. 4

5 An interpretation based upon known results of the high energy limit of scattering amplitudes in QCD is given in section 4. An outlook is presented in the concluding section. Some details of the calculations are presented in several appendices. Effective action for the interactions of Reggeized gluons Initially calculations of scattering amplitudes in Regge kinematics were performed by an iterative method based on analyticity, unitarity and renormalizability of the theory []. The s-channel unitarity was incorporated partly in the form of bootstrap equations for the amplitudes generated by Reggeized gluons exchange. But later it turned out that for this purpose one can also use an effective field theory for Reggeized gluons [9, 30]. We shall write below the effective action valid at high energies for interactions of particles inside each cluster having their rapidities y in a certain interval y = ln ǫ k + k ǫ k k, y y 0 < η, η << ln s. The corresponding gluon and quark fields are v µ x = it a v a µ x, ψx, ψx, [T a,t b ] = if abc T c. In the case of the supersymmetric models one can take into account also the fermion and scalar fields with known Yang-Mills and Yukawa interactions. Let us introduce now the fields describing the production and annihilation of Reggeized gluons [9] A ± x = it a A a ± x. 3 Under the global color group rotations the fields are transformed in the standard way δv µ x = [v µ x,χ], δψx = χψx, δax = [Ax,χ], 4 but under the local gauge transformations with χx 0 at x we have δv µ x = g [D µ,χx], δψx = χxψx, δa ± x = 0. 5 In quasi-multi-regge kinematics particles are produced in groups clusters with fixed masses. These groups have significantly different rapidities corresponding to the multi-regge asymptotics. In this case one obtains the following kinematical constraint on the reggeon fields A ± x = 0, ± = n µ ± µ, 6 n µ ± = δ µ 0 ± δµ 3. For QCD the corresponding effective action local in the rapidity y has the form [9] S = d 4 xl 0 + L ind, 7 where L 0 is the usual Yang-Mills Lagrangian L 0 = i ψ ˆDψ + Tr G µν, D µ = µ + gv µ, G µν = g [D µ,d ν ] 8 and the induced contribution is given below L ind = Tr L k ind + L GR ind, L k ind = µ A a + µ A a. 9 5

6 Here the Reggeon-gluon interaction can be presented in terms of Wilson P-exponents L GR ind = g + P exp g x + v + x dx + σa g P exp g x v x dx + σa + = v + gv + v + + g v + v + v v gv v + g v v v... σ A +A σ. 30 One can formulate the Feynman rules directly in momentum space [30]. For this purpose it is needed to take into account the gluon momentum conservation for induced vertices k 0 ± + k± k± r = 0. 3 Some simple examples of induced Reggeon-gluon vertices are ν 0+ a 0 c = q δ a0 c n + ν 0, ν 0ν + a 0 a c = q T c a a 0 n + ν ν 0ν ν + a 0 a a c = q n+ ν 0 n + ν n + ν T a a a 0 T c a a k + k+ In the general case these vertices factorize in the form k + n + ν 0, 3 + T a a a Ta c 0 a k k+ r ν 0ν...ν r+ a 0 a...a rc = r q n + νs Tr T c G a0 a...a r, 34 s=0 where T c are the color generators in the fundamental representation. In more detail, G a0 a...a r can be written as [30] G a0 a...a r = {i 0,i,...,i r} T a i 0 T a i T a i...t a ir k + i 0 k + i 0 + k + i...k + i 0 + k + i k + i r. 35 These vertices satisfy the following recurrent relations Ward identities [9] k + r ν 0ν...ν r+ a 0 a...a rc k + 0,...,k+ r r = n + νr if aara i ν 0...ν r + a 0...a i aa i+...a r c k+ 0,...,k+ i,k+ i + k r +,k+ i+,...,k r.36 i=0 With the use of this effective theory one can calculate the tree amplitude for the production of a cluster of gluons, or a gluon and a pair of fermions or scalar particles in the case of an extended supersymmetric model in the collision of two Reggeized gluons [30]. The square of the amplitude for three particle production integrated over the momenta of these particles is the main ingredient to construct the corresponding contribution to the BFKL kernel in the next-to-next-to-leading approximation using the methods of [3]. One can go to the helicity basis of produced gluons or fermions [33]. In principle it is also possible to calculate the loop corrections to the above Reggeonparticle vertices with the use of the effective action, however, in the present paper, we will use for this purpose the results for N = 4 SUSY amplitudes presented by Bern, Dixon and Smirnov in [3]. 6

7 3 BDS amplitudes in multi-regge kinematics As we have already remarked in the previous section, to find the next-to-next-toleading corrections to the BFKL kernel in N = 4 SUSY we need to calculate, apart from the amplitude for the transition of two Reggeized gluons to three particles, also the three loop correction to the gluon Regge trajectory, the two loop correction to the Reggeon-Reggeon-gluon vertex, and the one loop correction to the amplitude for the transition of two Reggeized gluons to two gluons or their superpartners. In this section we consider the corrections to the Regge trajectory and to the Reggeon-Reggeon-gluon vertex which can be obtained from the multi-regge asymptotics of the amplitude with the maximal helicity violation, calculated by Bern, Dixon and Smirnov BDS [3]. We also investigate the six point amplitudes 4 and 3 3 in multi-regge kinematics, thus preparing the comparison with QCD calculations to be carried out in the following section. Figure : Elastic scattering in the Regge asymptotics According to Ref. [3], in the case of maximal helicity violation the amplitude A n with n legs in the large-n c limit is factorized in the product of the tree result including the corresponding color structure and the simple scalar quantity M n. In the Regge limit s t the expression for M 4 having the singularities in s and t-channels can be simplified as follows see Fig. and Appendix A M = Γt s ωt µ Γt, 37 where µ is the renormalization point, ωt = γa 4 0 ln t µ + a 0 da a γa 4ǫ + βa, 38 is the all-order gluon Regge trajectory, as obtained from the BDS formula [35] for a verification by comparison with explicit calculations see further below, and ln Γt = ln t a da γa µ 0 a + βa + Ca + γa 8ǫ ζ a da a ln a γa a 4ǫ + βa + δa, 39 ǫ is the vertex for the coupling of the Reggeized gluon to the external particles. The t Hooft coupling is defined as follows a = αn c π 4πe γ ǫ 40 7

8 and the small parameter ǫ is related to the dimensional regularization 4 4 ǫ. The cusp anomalous dimension γa is known to all loops [0, 34, 3] and the functions βa, δa and Ca read [3] γa = 4a 4ζ a + ζ 4 a , 4 βa = ζ 3 a + 6ζ 5 + 5ζ ζ 3 a , δa = ζ 4 a +..., where ζn is the Riemann ζ-function Ca = ζ a +..., 4 ζn = k n. 43 k= Written as in Eq. 37 we can see that the asymptotic behavior of the M BDS amplitude corresponds to the Regge ansatz with the gluon trajectory j = + ωt given by the perturbative expansion ωt = + ln t µ + ǫ [ ζ 4 a + [ζ ln t µ 3ǫ ln t µ ǫ + 6ζ 5 + 5ζ ζ 3 3 ζ ] 3 a ] a The first two terms in this expansion are in agreement with the predictions in Refs. [, 3]. Note that in Ref. [3], where the BFKL kernel at NLO was calculated in N = 4 SUSY, initially the MS-scheme was used, and only later, in Ref. [3] the final result was also presented in the dimensional reduction scheme DRED. The NLO terms in the BDS expression for ωt can be obtained from Ref. [3] by converting it to the DRED scheme, where, apart from the finite renormalization of the coupling constant, one should also take into account in the loop the additional number of ǫ scalar particles for details see Appendix A and the recent paper [36]. The Oa 3 term in ωt, extracted from the BDS amplitude, corresponds to the three-loop correction to the gluon Regge trajectory needed when calculating the next-to-next-to-leading corrections to the BFKL kernel in this model. It is noteworthy to point out that the expression for M in the Regge kinematics coincides in fact with the exact amplitude and therefore it can also be written in the dual form t ωs M = Γs µ Γs. 45 After having described the elastic scattering amplitude we now focus on the BDS production amplitude see Fig. 3 and Appendix B with M 3 = Γt ωt s µ Γt,t,ln κ κ = s s s ωt s µ Γt, 46, 47 We thank A. V. Kotikov and E. M. Levin for helpful discussions regarding these redefinitions. 8

9 Figure 3: Production amplitude in the multi-regge regime in the multi-regge kinematics s s, s t t, 48 where all invariants s,s,s,t and t are negative. Due to the correct factorization properties of this amplitude the Reggeon particle particle vertex Γt in Eq. 46 is exactly the same as in the elastic amplitude in Eq. 39. The gluon Regge trajectory ωt in Eq. 46 also coincides with the one discussed above. The new component is the Reggeon-Reggeon-particle vertex given by ln Γt,t,ln κ = γa 6 γa t 6 ln γa t 6 ζ ln κ µ a 0 ωt + ωt da a ln a γa a 4ǫ + βa ǫ da γa a + βa 4ǫ a 0 + δa ln κ µ. 49 It is now possible to analytically continue this 3 production amplitude to the physical region where the invariants s, s and s are positive see Fig. 4a s iǫ M 3 = Γt µ ωt s iǫ Γt,t,ln κ iπ µ ωt Γt. 50 Figure 4: Physical channels for the one particle production amplitude A similar continuation to another physical region can be performed in the case when s is positive but s and s are negative see Fig. 4b: M 3 = Γt ωt s µ Γt,t,ln κ + iπ 9 ωt s µ Γt. 5

10 Figure 5: Analytic representation for M 3, exhibiting the analytic structure in the energy variables Using a dispersive representation illustrated Fig. 5 in the following, we will refer to this type of representation as analytic representation, it can be easily verified that in all physical regions the amplitude can be written as follows M 3 Γt Γt = s µ ωt ωt sκ ωt µ 4 c + s µ ωt ωt sκ ωt µ 4 c,5 where c c sin πωt φ Γ = Γt,t,ln κ sin πωt ωt, 53 = Γt,t, κ sin πωt φ Γ sin πωt ωt, 54 and φ Γ is the phase of Γ, i.e. Γt,t,ln κ iπ = Γt,t,ln κ e iφ Γ. 55 In this dispersion-type representation we have used in all physical channels the reality condition for the produced particle κ s s s = k, 56 where k is the transverse component of the particle momentum k p A = k p B = 0. The fact that there exists a solution for the coefficients c and c proves that the BDS amplitude for the production of one particle in multi-regge kinematics has the correct multi-regge form. This is encouraging, and we proceed now to study the production of two particles in multi-regge kinematics, for which we use the M 4 BDS scattering amplitudes. We have first checked that the planar BDS amplitude for two particle production having singularities only at positive values of the invariants s,s,s,s 3,t,t,t 3 has the correct multi-regge form in the multi-regge kinematics in the region where all invariants s,s,s,s 3,t,t,t 3 are negative see Fig. 6 and Appendix C M ωt 4 Γt Γt 3 = s µ Γt,t,ln κ where we have used the notation ωt s µ Γt 3,t,ln κ 3 ωt3 s3 µ 57 κ = s s s 0, κ 3 = s s 3 s

11 Figure 6: Production of two particles The s 0 and s 3 invariants are the squared masses of the corresponding three particles in their center-of-mass system. In a similar way to what we did in the 3 case, the BDS 4 amplitude can be analytically continued to several physical channels, each of them corresponding to the different signs of the invariants s,s 0,s 3,s,s,s 3. To begin with, in the region see Fig. 7a the amplitude has the form s,s 0,s 3,s,s,s 3 > 0 59 M 4 Γt Γt 3 = ωt s µ Γt,t,ln κ iπ ωt s µ Γt 3,t,ln κ 3 iπ where we can replace κ and κ 3 by their values on the mass shell ωt3 s3 µ, 60 κ s s s 0 = k, κ 3 s s 3 s 3 = k. 6 Here k = q q and k = q q are the momenta of the produced particles, see Fig. 7. In the physical region, represented in Fig. 7b, where s,s 0,s 3 > 0, s,s,s 3 < 0, 6 we obtain, for the multi-regge asymptotics of the BDS amplitude, the following structure M 4 Γt Γt 3 = ωt s µ Γt,t,ln κ + iπ In a similar way, in the region see Fig. 7c ωt s µ Γt 3,t,ln κ 3 iπ ωt3 s3 µ. 63 s,s 3,s > 0, s 3,s,s 0 < 0 64

12 Figure 7: Physical channels for two particle production we have M 4 Γt Γt 3 = ωt s µ Γt,t,ln κ iπ ωt s µ Γt 3,t,ln κ 3 + iπ ωt3 s3 µ. 65 We can now attempt to present M 4 in terms of the following dispersion ansatz represented in Fig. 8 M 4 Γt Γt 3 = ωt ωt s s 0 κ ωt ωt 3 sκ κ 3 µ µ 4 ωt3 ωt s3 s 3 κ ωt ωt 3 sκ κ 3 µ µ 4 ωt ωt s s 0 κ ωt ωt 3 sκ κ 3 µ µ 4 ωt ωt s 3 s 3 κ ωt3 ωt 3 sκ κ 3 µ s3 µ ωt3 ωt s µ µ 4 ωt ωt sκ κ 3 µ 6 ωt3 µ 6 d ωt µ 6 d ωt3 µ 6 d 3 ωt µ 6 d 4 ωt d 5, 66 where the coefficients d i=,,3,4,5 should be real. By comparing this dispersive representation with the previous expressions for the BDS amplitude in the three physical regions it is possible to extract the coefficients. They read d = c κ c κ 3, d = c κ c κ 3, d 3 + d 4 = c κ 3 c κ, d 5 = c κ c κ 3, 67

13 Figure 8: Analytic representation of the amplitude M 4 where, in fact, c κ and c κ were defined in Eqs. 53, 54: with cf.49 c κ = Γ sin πωt φ Γ sin πωt ωt, 68 c κ = Γ sin πωt φ Γ sin πωt ωt, 69 c κ 3 = Γ 3 sin πωt φ Γ3 sin πωt ωt, 70 c κ 3 = Γ 3 sin πωt φ Γ3 sin πωt ωt, 7 Γ = Γt,t,ln κ, Γ 3 = Γt,t 3,ln κ 3. 7 We note that for the coefficients d 3 and d 4 only their sum can be determined from the three physical regions previously discussed. However, an attempt to fix separately these two coefficients from the multi-regge asymptotics in the physical region see Fig. 7d s,s > 0, s,s 3,s 0,s 3 < 0 leads to a disaster: the corresponding equations do not have any solution. The reason for this is that the BDS amplitude in this region does not have the correct Regge factorization see the discussion in section 4. According to Appendix C its asymptotics here is M 4 Γt Γt 3 = C s µ ωt Γt,t,ln κ iπ ωt s µ Γt 3,t,ln κ 3 iπ where the coefficient C is given by [ γk a C = exp π + iπ ln t ] t 3 4 k + k µ 3 ωt3 s3 µ,73. 74

14 As it will be explained in the next section, this expression does not satisfy the dispersion relation written diagrammatically in Fig. 8 and it is in agreement with the BFKL approach only in the one-loop approximation. In Appendix C we write down the amplitude M 4 also in the quasi-multi-regge kinematics, where the variable s is fixed. Figure 9: Three particle transition To continue our analysis of the BDS six point amplitude we now discuss the asymptotics of M 3 3 see Fig. 9. According to Appendix D in the multi-regge region where all invariants s,s,s 3,s 3,s 0,s t are large and negative, its asymptotics is similar to the corresponding asymptotics of the M 4 amplitude, i.e. M 3 3 Γt Γt 3 = s µ ωt Γt,t,ln κ ωt s µ Γt 3,t,ln κ 3 ωt3 s3 µ. 75 This BDS amplitude can be now analytically continued to the physical region where the invariants s,s,s 3,s,s 0,t are positive see Fig. 0a. The resulting amplitude can be written as M 3 3 Γt Γt 3 = s µ ωt Γt,t,ln κ iπ ωt s µ Γt 3,t,ln κ 3 iπ ωt3 s3 µ, 76 Similarly, the analytic continuation to the region where s,s,t < 0 and s,s 3,s 0 > 0 see Fig. 0b is of the form M 3 3 Γt Γt 3 = s µ ωt Γt,t,ln κ + iπ ωt s µ Γt 3,t,ln κ 3 iπ ωt3 s3 µ. 77 Finally, the continuation to the region where s 3,s 0,t < 0 and s,s,s > 0 see Fig. 0c reads M 3 3 Γt Γt 3 = s µ ωt Γt,t,ln κ iπ ωt s µ Γt 3,t,ln κ 3 + iπ 4 ωt3 s3 µ. 78

15 Figure 0: Physical regions for the amplitude M 3 3 As we did in the M 4 case it is possible to write the dispersion relation for M 3 3, which includes five contributions, shown in Fig., and try to calculate the real coefficients e,e,e 3 + e 4,e 5. Once again, we find that it is impossible to find separately the coefficients e 3 and e 4 from the BDS amplitude M 3 3 Γt Γt 3 = s C µ ωt Γt,t,ln κ + iπ ωt s µ Γt,t,ln κ 3 + iπ ωt3 s3 µ,79 calculated in the physical region where s,s 3,s 3,s 0 < 0 and s,t > 0 see Fig. 0d. The reason for this drawback is the same as before: the absence of a correct Regge factorization for the BDS amplitude. In the next section, using the BFKL approach, we shall discuss the reason for this problem, related to the absence of the BFKL-like contributions plotted in Fig.. 4 Regge cuts and the breakdown of factorization The results of the previous section can best be understood if we confront them with the known high energy behavior of QCD scattering amplitudes in the Regge limit. In the LLA, the high energy behavior of the QCD scattering amplitudes is the same as in the supersymmetric case. A key element in analyzing the high energy limit are the Steinman relations [37] which forbid the existence of simultaneous energy discontinuities in overlapping channels. As an illustrative example of the Steinmann relations, consider the 3 amplitude shown in Fig. 5: obviously, the produced particle in the central region can form resonance states with particle A or with particle B, but not simultaneously with both of them. As a result, the scattering amplitude cannot have 5

16 Figure : Analytic representation of the amplitude M 3 3 simultaneous discontinuities in the energy variables s and s. The way in which this restriction is implemented into physical scattering amplitudes is that, in the double Regge limit, the signatured amplitude can be written as a sum of two pieces, one of them with singularities in the s and in the s channels, the other one in the s and in the s channels. In general, there are cut singularities both in the right and left half energy planes, and one has to form signatured combinations. The general structure is cf. 5: A 3 β A t β B t [ s + µ [ s µ = αt αt αt αt + τ τ s µ + τ τ s µ αt αt ] [ s µ αt + τ s µ αt ] Ṽ t,t,κ αt αt ] [ s αt s αt ] µ + τ µ Ṽ t,t,κ.80 Here α i t i = + ωt i, τ τ are the signatures of the t t channels, and as usual, s α = s α e iπα. In this representation, the energy singularities are explicit, i.e. all phase factors are contained in the energy factors, and the functions βt, V i are real valued functions. With the signature factors ξ i = e iπα it i + τ i, ξ ij = e iπα it i α j t j + τ i τ j sinπα i t i sin πα i t i α j t j 8 we can rewrite the expression for A 3 : A 3 β A t β B t = + s αt αt s αt ξ ξ V t,t,κ µ s µ µ αt αt s µ αt ξ ξ V t,t,κ. 8 6

17 The generalization to the 4 amplitude consisting of five different pieces is illustrated in Fig. 8, and from Eq. 66 one easily obtains the analogue of 8: A 4 = β A t β B t 3 s αt αt s0 αt αt 3 s µ µ s3 αt3 αt s3 µ µ s αt αt s0 µ µ s αt αt 3 s3 µ µ s3 µ αt3 αt s µ µ αt αt s αt3 ξ ξ 3 ξ 3 V t,t,κ V t,t 3,κ 3 µ αt αt 3 s µ αt3 αt s αt αt s µ µ αt ξ 3 ξ ξ V t,t,κ V t,t 3,κ 3 αt3 ξ ξ 3 ξ 3 V t,t,κ V t,t 3,κ 3 αt ξ 3 ξ 3 ξ V t,t,κ V t,t 3,κ 3 αt ξ 3 ξ ξ V t,t,κ V t,t 3,κ 3.83 The analogue for the 3 3 process again consists of five pieces which are shown in Fig.. Decompositions of this kind have first been derived from simple models which contain only Regge poles hybrid Feynman diagrams in scalar ϕ 3 theory [38] or the dual Veneziano 6 point amplitude, B 6 [39], and from studies of dispersion relations and generalized Froissart-Gribov partial wave representations [40]. From these discussions it has also been recognized that the analytic decomposition into a sum of terms as in 8 and 83 is consistent with a factorizing form. For the 3 case we can write: A 3 β A t β B t = s αt µ ξ V t,t,κ s αt µ ξ. 84 Here the important point to be stressed is that the production vertex function V contains phases in contrast to the real-valued functions V i in 8, and it has cut singularities in the κ-plane. Similarly for the 4 case: A 4 β A t β B t 3 = ξ s αt µ s αt V t,t,κξ µ s3 αt3 V t,t 3,κ 3 ξ 3 µ,85 with the production vertex function from 84. As a result, for this class of Regge-pole models the production amplitudes, in the multi-regge limit, can be written either in the analytic form sum of terms with simple analytic properties and real-valued vertex functions V i or in the factorized form with the production vertices V containing phases and singularities in κ. There is an important consequence for planar amplitudes. Planar amplitudes have right hand cut singularities only, and their phases follow from 8 or 84 by simply dropping all terms containing factors τ i. Starting from the factorized representation 84 in the physical region all energies positive, one can analytically continue into the unphysical region, where all energies variables are negative and well separated from their threshold singularities: all phases inside the production vertex then disappear, and the vertex function turns into a real-valued function. This clearly supports the hope that in the unphysical region planar scattering amplitudes have a simple form with real-valued factors. When turning to QCD the situation is slightly more complicated. In the LLA the real part of the n scattering amplitude is well known to have the factorized form of 7

18 Eq. 4, and is in agreement with our previous result in 84 and 85. However, when turning to the imaginary parts i.e. to the energy discontinuities of the production amplitudes, a new piece appears which destroys the simple factorization property. The best way of understanding the appearance of this new piece is the use of s-channel unitarity in the physical region where all energies are positive. Starting from the analytic representation of the scattering amplitude A n, it is possible to determine, in QCD, the partial waves from single energy discontinuities and unitarity equations [4, 4]. As the simplest example, let us consider, in the LLA, the 3 amplitude, consisting of the two terms illustrated in Fig. 5 and written in Eq. 8. Taking the discontinuity in s, only the first term contributes. Making use of the unitarity equation and invoking, for the ladder diagrams in the t channel, the BFKL bootstrap condition we find the partial wave V in the LLA: V = ωt ωt ωt a ln κ µ ǫ. 86 Similarly, the discontinuity in s leads to V = ωt ωt ωt a ln κ µ ǫ. 87 As a further test, one could also compute, from the corresponding unitarity equation, the single discontinuity in s: here both partial waves V and V contribute, and the result is in agreement with 86 and 87. For the 4 amplitude we start from the ansatz 83 with unknown partial wave functions denoted by W i=,,3,4,5. For reasons which will become clear soon we introduce, for the t channel, the Sommerfeld Watson integral dω /πi, replacing everywhere αt + ω. We have the five single discontinuities in s, s, s 3, s 0, and s 3 which allow to find the partial waves W i=,,3,4,5. The results, apart from an overall color tensor and from the common factor are the following: s t t t 3 g δ λa λ A δ λb λ B, 88 W = V t,t,κ ω ωt V t,t 3,κ 3, 89 W = V t,t,κ ω ωt V t,t 3,κ 3, 90 W 3 = V t,t,κ ω ωt V t,t 3,κ 3 + sin πωt 3 sin πω ωt V p + V cut, 9 W 4 = V t,t,κ ω ωt V t,t 3,κ 3 sin πωt sin πω ωt 3 V p + V cut, 9 W 5 = V t,t,κ ω ωt V t,t 3,κ The three amplitudes W, W and W 5 are products of the production vertices V i in 86 and 87, found in the 3 case, whereas the amplitudes W 3 and W 4 contain, We emphasize that the same results are obtained if one starts from the double discontinuities in s and s : using unitarity conditions and making use of generalized bootstrap conditions, one again arrives at 86. This is a crucial test of the selfconsistency of this unitarity-based approach. 8

19 in addition to the products of production vertices V i, the extra pieces, V p + V cut which cancel in the real part to leading order of the sum W 3 + W 4, and they also cancel in the sum over all single discontinuities in the physical region. The term V cut contains Regge cuts and cannot be written as a simple product of vertices for the produced gluons d and d. It takes the form of BFKL-like ladder diagrams in the color octet channel, and it is illustrated in Fig. left figure: V cut = t g 4 d kd k π 6 q k + q C µk + q,k + q G 8A k + q, k;k + q, k ;ω C µ k + q 3,k + q k + q 3, 94 where C µ denotes the effective Reggeon-Reggeon-gluon vertex, and G 8 A is the BFKL Green s function in the color octet channel, satisfying the integral equation ω G 8A k + q, k;k + q, k;ω = π 3 δ k k k k + q + k k + q K G 8 A k + q, k;k + q, k, 95 where K denotes the BFKL kernel in the color octet channel, containing both real emission and the gluon trajectory. In lowest order in the coupling, and for equal helicities of the two produced gluons, V cut equals: with V 0 cut = π t [ ωt + ωt 3 + a ln κ 3 µ ] ǫ q 3 96 κ 3 = k + k. 97 This term V cut not only violates the simple factorization, but also, when computed beyond the one loop approximation, will be shown to disagree with the BDS formula. Finally, the subtraction term V p removes the Regge pole piece inside V cut, and it is of the form: 4sin πω [ ωt + ω + a ln κ µ ] ǫ [ ω ωt ω + ωt 3 + a ln κ 3 µ ].98 ǫ Before we compare with the BDS formula, let us remark on a few consequences of these leading order QCD results. First, by inserting the results of into the full amplitude 83 and approximating the phase factors by, e.g., e iπωt + the leading-order real part, one finds substantial cancellations, after which the final result is in full agreement with the simple factorizing expression 4. Note that in this result not only the extra pieces in V 3 and V 4 cancel, but also the terms proportional to lnκ/µ in the production vertex. Next, if we take the sum over the single discontinuities in s, s 0, s 3, and s, again the extra piece V cut in V 3 and V 4 cancels. The same is true for the full sum of all single discontinuities i.e. after adding the discontinuities in s and s 3 : this means that in the NLO calculation of A 4 which includes, apart from the corrections to the real part, also the leading order imaginary part these pieces, once more, are not visible. This cancellation in the imaginary part, however, will not work if we restrict ourselves to the planar amplitude i.e. neglect all terms proportional to τ and consider the analytic continuation into the physical region: s > 0, s > 0, s,s 0,s 3,s 3 < 0. When forming, in this region, the sum over all nonvanishing discontinuities, this sum includes only s and s and not s, s 3, 9

20 Figure : BFKL contributions to the amplitudes M 4 and M 3 3 s 0, s 3, and thus leads to a nonvanishing coefficient in front of the Regge-cut piece V cut. A completely analogous discussion applies to the case 3 3 Figs. 9, in the multi-regge region for a detailed discussion of the analytic representation see [4]: here, again, the pieces labelled by,,5 are normal and fit into the factorization pattern, whereas the terms 3 and 4 have the extra piece shown in Fig. right figure, which is described in terms of the color octet BFKL equation, contains Regge cut singularitites, and breaks the factorization. In analogy with the 4 case, these cut pieces do not show up in the real part leading order of the full 3 3 amplitude or, in the physical region where all energies are positive, in the sum of all single energy discontinuities. It is, again, only in the planar approximation of the physical region s,s > 0, s,s 3,s 3,s 0 < 0, where the sum of the single energy discontinuities does not include the energies s 3, s 0 and hence leads to a nonzero coefficient in front of the piece V cut. Finally, let us return to the BDS amplitude discussed in the previous section. Beginning with the case 3 in the physical region, we find agreement between the leading order terms of c and c in 53 and 54 and V and V in 86 and 87. The same applies to all other physical regions. In the case of 4 we immediately turn to the physical region s > 0, s > 0, s,s 0,s 3,s 3 < 0 where, in the QCD calculations, the imaginary part contains the factorization breaking term V cut a BFKL ladder with the octet quantum numbers in the t-channel. In the BDS amplitude, we first consider the one loop approximation, i.e. in 74 we expand the phase e iδ +iδ of the coefficient C. Comparison with the QCD in 96, shows agreement. To be more precise, in the discussion after 65 we had found that the factorizing ansatz had no solution in all physical regions. If, on the other hand, we allow for extra pieces in d 3 and d 4 analogous to 9 and 9, we find solutions which, in the one loop approximation, agree with 96 and 98. However, beyond the one loop approximation, the BDS formula does not reproduce the BFKL ladder contained in V cut. The same is found for the 3 3 case: in the one loop approximation, the BDS formula reproduces the factorization calculations, in particular with V 0 cut breaking piece V 0 cut, but in higher loops it fails. We therefore conclude that beyond the five point amplitude and beyond the one loop approximation the BDS formula does 0

21 not agree with the leading logarithmic high energy behavior in all physical regions. Discrepancies in the BDS finite pieces for six gluon amplitudes, starting at two loops, were also found in [44]. In a particular kinematic configuration, and for a very large number of external gluons at strong t Hooft coupling, the finite pieces of the BDS ansatz also failed when compared to the results of [45]. 5 Conclusions In this paper we have assembled the ingredients needed for the three loop corrections NNLO to the BFKL kernel in N = 4 SYM theory at large N c. Following earlier calculations we can obtain the kernel from unitarity sums, i.e. by computing squares of production amplitudes, keeping in mind that at large N c the contributing diagrams belong to the cylinder topology. Figure 3 illustrates the production vertices which enter the three loop calculation. Figure 3: Diagrams contributing to the BFKL kernel in NNLLA

22 Elements in the first two lines are known, whereas the building blocks in the third line are new: they can be and partly have been computed from the effective action summarized in section. For most of the cases, also the BDS formula can be used. The NNLO gluon trajectory function follows from the scattering amplitude first column; details are described in Appendix A. For the 3 case second column, the two loop approximation to the production vertex can be read off from the analysis presented in section 3 c.f. Eq. 49. In column 3, we should take into account the imaginary part of the 4 amplitude in the multi-regge region, and we have to include the vertex: Reggeon + Reggeon gluons in the one loop approximation. Based upon the analysis carried out in sections 3 and 4 we trust that the BDS formula for the maximal helicity violating case can be used. For the non-maximal violating cases in column 3 one can use the results of [46]. Finally, in column 4 we encounter the Born vertex: Reggeon + Reggeon 3 gluons in the quasi-multiregge kinematics. This vertex has been obtained in [30] by means of the effective action. A remark is in place on the Regge-cut contribution illustrated in Fig. and discussed in section 4. In addition to the corrections to the production vertex functions which are illustrated in Fig. 3, we still have to take into account those corrections to the production amplitude in the multi-regge limit which do not fall into the class of loop corrections to the production vertices: in NLO these are just the Regge-cut contributions to the imaginary part in the 4 and the 3 3 cases which we have discussed in the previous section. The diagrams contributing to the BFKL Pomeron in the large N c limit belong to the cylinder topology: two examples are illustrated in Fig. 4, and, to begin with, we consider the discontinuity due to the 4-particle intermediate state. In the left figure, on both sides of the discontinuity cut, we have the 3 3 production amplitudes continued into the physical region of a 4 process c.f. rhs of Fig., and in the right hand figure we recognize a configuration where the 4 amplitude has to be evaluated in a region with negative energies. As discussed before, in the latter case the non-factorizing pieces of the 4 production amplitude do not cancel. If these contributions would survive in the total cross section, the NNLO BFKL Pomeron would receive a four-reggeon cut contribution, and the simple ladder structure would be lost. There are, however, reasons to expect that, in the large N c limit, the sum of these contributions might cancel in the total cross section. Namely, in addition to the contribution of the 4-particle intermediate state, we also need other cuts, which, for example, run across one of the ladders or along one of the reggeized gluons above or below the cylinder. These different cuts provide similar contributions, but they come with different signs. It is likely that, similar to the AFS cancellation of Regge cuts in planar amplitudes, the four reggeon cut contributions cancel in the sum. We will study this in the subsequent part of our investigation. Acknowledgements: Part of this work has been done while one of us ASV has been visiting the II. Institut f. Theoretische Physik, University Hamburg. The hospitality is gratefully acknowledged. One of us LNL wishes to thank the Issac Newton Institute for the invitation to participate in the program String Fields, Integrability and Strings. LNL is supported by the RFBR grants MNTI-a, a, andrsgss Note to be added: After our paper had been submitted, another study appeared which, in some parts, parallels our investigation [39]. Like ours, it studies several different Regge limits of the BDS amplitudes. In the unphysical region negative energies, the results on the energy dependence are fully consistent with ours. In contrast

23 Figure 4: Cylinder-type topologies in the unitarity sums for the total cross section: the intermediate states discontinuity cuts are obtained by slicing the cylinders in all possible ways across the intermediate momenta p A and p B to our paper, however, in [48] the continuation into physical regions has not been investigated, and the conflict with QCD calculations was not found. A The amplitude Let us write the BDS amplitude for the general case of n legs see [3]: ln M n = a l f l ǫ l= Îl n lǫ + F n 0 + C l + E n l ǫ, A. where E n ǫ can be neglected for ǫ 0, the values of the constants are C = 0, A. C = ζ /, A.3 f l ǫ = f l 0 + ǫf + ǫ f l, A.4 f l 0 = 4 γl K, γ K is the cusp anomalous dimension [47], A.5 f = aζ 3 / + a ζ 5 + 5ζ ζ 3 /3, A.6 Î n ǫ = ǫ n µ i= s i,i+ ǫ, A.7 and the finite remainders F n are expressed in terms of logarithms and dilogarithms. For the elastic scattering amplitude case we have Î 4 ǫ = ǫ + s t ln ǫ µ 4 3 ln s t + ln µ µ, A.8

24 Therefore F 4 = s ln t + 4ζ. Î 4 ǫ + F 4 = ǫ + ln t ǫ + ln s ǫ A.9 t ln µ + 4ζ. A.0 As a result we obtain for M 4 Regge-type behaviour, as already discussed in the main part of our paper, with the gluon Regge trajectory given by ωt = a ǫ t ln µ + a ζ ǫ t ln µ ζ A. Note that this result at two loops is in agreement with the direct calculations [, 3] based on the BFKL approach []. Indeed, in Ref. [3] the following expression for the gluon Regge trajectory was obtained in the MS-scheme using the same notations: ω MS t = a ǫ [ t ln µ + a 6 ζ ǫ t ln µ + 9 ζ ] 3. A. The contribution of the scalar loop to this trajectory is proportional to the contribution of the fermion loop [3] ω s MS t = n s 4n q ǫ ωq t = n a [ t MS s ln 4 ǫ µ 8 3 ǫ t + ln µ 5 ], A.3 9 where n s is the number of scalar fields transforming according to the adjoint representation of the gauge group. For the transition from the MS-scheme to the dimensional reduction DRED scheme, which respects N = 4 supersymmetry, one should first increase the number of scalar fields n s 6 + ǫ, A.4 because, in the pure gluonic contribution, n = ǫ for the gluon fields was taken into account after performing the dimensional regularization 4 4 ǫ. This gives the additional contribution to ω MS t ω MS t = a ǫ 8. A.5 3 After that the subsequent finite renormalization of the coupling constant needed for the transition between the MS and DRED schemes a a 6 a A.6 leads to the above result for the trajectory ω MS t ωt = a ǫ [ t ln µ + a ζ ǫ t ln µ ζ ] 3. A.7 Concerning the residues Γt of the Regge pole, they have been calculated in the one-loop approximation in QCD [3]. In supersymmetric models the helicity nonconserving contribution of each of the colliding gluons is cancelled, in accordance with the BDS ansatz. 4

25 B The 3 amplitude For the 3 production amplitude we have Î 5 ǫ = 5 ǫ + ǫ ln s s s t t µ 0 ln s 4 µ + s ln µ + s ln µ + t ln µ + t ln µ F 5 = s ln ln t 4 s s 4 ln t ln s t s 4 ln s s ln s t 4 ln s t ln t s 4 ln t s ln s t ζ. Thus the total contribution in multi-regge kinematics is I 5 ǫ + F 5 = 5 ǫ + ln s µ ǫ ln t µ 4 κ ln µ + κ ln µ ln t t µ 4 ǫ + ln s µ ǫ ln t µ 4 ln t t ζ., B. B. B.3 In this way we obtain the Regge factorization of the production amplitudes, discussed in the main text. In the one-loop approximation in QCD the Reggeon-Reggeon-gluon vertex contains, apart from the Born structure proportional to the vector Cq,q, also the contribution proportional to the gauge-invariant vector p A s p B s [3]. In the supersymmetric theories this contribution is cancelled, in agreement with the BDS ansatz. One can calculate also the production amplitude in the quasi-elastic kinematics, where s s s t,t,k. The amplitude here has the usual Regge factorization. C The 4 amplitude In the case of the 4 transition we have Î 6 ǫ = 3 ǫ + ǫ ln s s s s 3 t t 3 µ ln s 4 µ + s ln µ + s ln µ + s ln 3 µ + t ln µ + t ln 3 µ,c. F 6 = s ln ln t 3 s 0 s 0 ln t 3 ln s 3 t t ln s 3 ln s s 3 s 3 ln s ln s s 0 s 0 ln s ln t t t ln t ln s s 3 s 3 Li ss s 0 s 3 Li t 3s t s 0 Li t s 3 t s ln t s where the dilogarithm function is defined as Li z = ln t s z 0 dt t ln t. C.3 ln s 3 s ζ, C. In multi-regge kinematics it is natural to introduce the independent variables s,s,s 3, κ = s s s 0, κ 3 = s s 3 s 3 5, Φ = s s 3 s 0 s 3. C.4

26 Note that the variable Φ is unity in the region where all above invariants are negative, but Φ = exp πi in the physical region where s,s > 0, s 0,s 3 < 0. In the multi-regge kinematics we obtain the following general result: I 6 ǫ + F 6 = 3 ǫ 4 ln Φ ln Φ ln t t 3 s Li Φ + ln s µ ǫ ln t µ + ln s µ ǫ ln t µ + ln s 3 µ ǫ ln t 3 µ 4 ln κ + ln κ 3 µ µ + ln κ 3 µ + ln κ µ ln t t 3 µ 4 ǫ 4 ln t t µ 4 ǫ ln t + ln t 3 t t + 7 ζ. C.5 At first sight the arguments of the dilogarithm functions in the multi-regge kinematics are either 0 or, and we can use the relations Li 0 = 0, Li = ζ. C.6 However, in the physical region s,s > 0, s,s 3,s 0,s 3 < 0 it is needed to be cautious: we should analytically continue the expression fφ = Li Φ, Φ = ss s 0 s 3 C.7 from the region Φ to the region Φ e πi along a unit circle. After this continuation we obtain fφ = Φ Taking into account also the relation 0 dt t ln t + πi Φ which is valid in the multi-regge kinematics, where dt t πiln Φ. C.8 s Φ = k + k, C.9 s s 0s 3 s k + k, C.0 we obtain the following result in the physical region s,s > 0, s,s 3,s 0,s 3 < 0: I 6 ǫ + F 6 = 3 ǫ + πiln t t 3 k + + ln s k µ µ + ln s µ ǫ ln t µ + ln s 3 µ ǫ ln t 3 µ 4 + ln κ µ ln t µ + ln t µ + ǫ ln κ 3 µ 4 ln t t + ln t 3 t ǫ ln t µ ln κ µ + ln κ 3 µ ln t µ + ln t 3 µ ǫ + 7 ζ + π. C. It is possible to derive, from the BDS amplitude, an expression for M 4 in the one-loop approximation for the quasi-multi-regge kinematics, where s >> s,s 3 >> s t,t,t 3. In this case it is convenient to introduce Sudakov variables for the momenta of the two produced particles k r = β r p A + α r p B + k r, k r = k r, C. 6

27 where β β µ s, α α µ s, sα rβ r = k r q q q 3 µ. C.3 We can express various invariants in terms of these variables s sβ + β α + β k + k, C.4 s α s, s 3 β s, s 0 α + α s, s 3 β + β s. C.5 The expression for the function fǫ = I 6 ǫ + F 6 for Φ = in the quasi-multi- Regge kinematics can be obtained by adding an aditional term fǫ fǫ + f, C.6 where f = ln sκ κ 3 st t 3 ln s s s 3 s 0 s 3 µ + ln s 0κ ln t 3s s s s s t s 0 s 3 + ln s 3κ 3 ln t s 3 s s s 3 s t s 0 s sκ κ ln s 0κ s s s 3 ln s 3κ 3 s s ln s 3 s ln sκ κ 3 ln s 0s 3 κ κ 3 s s s 3 s s 3 s ln s 0κ ln s 3κ 3 + ζ s s s 3 s Li ss s 0 s 3 Li t 3s t s 0 Li t s 3 t s 3. C.7 Note that the expression for f in the quasi-multi-regge kinematics does not contain large logarithms, because the arguments of all logarithms and dilogarithm functions are of the order of unity. It is proportional to the logarithm of the amplitude for the transition of two Reggeized gluons into two particles with the same helicity. The transition to particles with opposite helicity in the one-loop approximation can be found in Ref. [46]. These transition amplitudes are needed for the calculation of the next-to-next-to leading corrections to the BFKL equation. D The 3 3 amplitude Here we consider the BDS amplitude M 6 in the channel corresponding to the transition 3 3 with the following invariants see Fig. 9: s = p A + k + p B, s = p A + k, s 3 = p B + k, D. s 3 = k + p B, s 0 = p A + k, t = p A p A, t 3 = p B p B, t = p A p A k, t = p A + k p A, u = k k. The functions Î 6 ǫ and F 6 in this case are given by [3]: Î D. D.3 D.4 6 ǫ = 3 ǫ + ǫ ln s s 3 s 3 s 0 t t 3 µ ln s 4 µ + s ln 3 µ + ln s 3 µ + s ln 0 µ + ln t µ + t ln 3 µ, D.5 F 6 = ln s s ln s 3 s ln s 3 t ln t 3 t ln t 3 ln s 3 t t 7

28 ln s 3 s ln s 0 s ln s 0 t ln t t ln t ln s t t Li s s 3 st Li s 3s 0 t s Li + ln t + ln t + 4 s 4 t ln t 4 s In multi-regge kinematics t t 3 t t 9 ζ. D.6 s s, s 3, t t, t, t 3 > 0 D.7 we have the relations s = s s t κ κ 3 which allow us to simplify the above expressions I 6 ǫ + F 6 = 3 ǫ + ln s µ + ln s 3 µ + ln κ 4 µ ǫ ln t 3 µ 4 ln t t ln t t + ln t 3 t References µ 4 ǫ, s 0 = s t κ, s 3 = s 3 t κ 3, D.8 ǫ ln t µ ln κ + ln t µ + ln κ 3 µ + ln κ 3 µ µ ǫ ln t µ ln t t 3 µ 4 ǫ + 7 ζ. D.9 [] L. N. Lipatov, Sov. J. Nucl. Phys ; V. S. Fadin, E. A. Kuraev, L. N. Lipatov, Phys. Lett. B ; E. A. Kuraev, L. N. Lipatov, V. S. Fadin, Sov. Phys. JETP ; ; I. I. Balitsky, L. N. Lipatov, Sov. J. Nucl. Phys [] V. S. Fadin, R. Fiore, M. I. Kotsky, Phys. Lett. B [3] A. V. Kotikov, L. N. Lipatov, Nucl. Phys. B [4] L. N. Lipatov, Nucl. Phys. B [5] L. N. Lipatov, Phys. Lett. B [6] L. N. Lipatov, Sov. Phys. JETP [7] J. Bartels, Nucl. Phys. B ; J. Kwiecinski, M. Praszalowicz, Phys. Lett. B [8] L. N. Lipatov, Phys. Lett. B [9] L. N. Lipatov, Nucl. Phys. B [0] L. N. Lipatov High energy asymptotics of multi-colour QCD and exactly solvable lattice models, Padova preprint DFPD/93/TH/70, hep-th/93037, unpublished. [] L. N. Lipatov, JETP Lett ; L. D. Faddeev, G. P. Korchemsky, Phys. Lett. B [] V. S. Fadin, L. N. Lipatov, Phys. Lett. B ; M. Ciafaloni and G. Camici, Phys. Lett. B