Light-by-light scattering in Double-Logarithmic Approximation

In the present paper we consider the elastic 2 ->2 -scattering of virtual photons at high energies in the forward kinematics at zero and non-zero values of t. Accounting for both gluon and quark double-logarithmic (DL) contributions to all orders in the QCD coupling, we obtain explicit expressions for amplitudes of this process in Double-Logarithmic Approximation (DLA). First we keep the QCD coupling fixed and then account for running coupling effects. Applying the saddle-point method to the obtained expressions for the scattering amplitude, we calculate the high-energy asymptotics of the amplitude, which proved to be of the Regge form. The Reggeon bears the vacuum quantum numbers and therefore it is a new, DL contribution to Pomeron. Comparison of the DL Pomeron to the BFKL Pomeron shows that contribution of the DL Pomeron to the high-energy asymptotics is of the same order as contribution of the BFKL Pomeron, so the DL Pomeron should be taken into account together with the BFKL Pomeron. We estimate the applicability region for the asymptotics of the light-by-light scattering amplitude, where the the DL Pomeron can reliably represent the parent amplitude.


I. INTRODUCTION
Elastic scattering of virtual unpolarized photons γ * (p) γ * (q) → γ * (p ′ ) γ * (q ′ ) (1) at high energies in the forward kinematics where Q 2 1 = |p 2 | ≈ |p ′2 | and Q 2 2 = |q 2 | ≈ |q ′2 |, has been an interesting object for theoretical investigation. The point is that theoretical investigation of hadronic reactions involves parton distributions. They inevitably involve non-perturbative phenomenological contributions. In contrast, the process (1) can be studied with the means of Perturbative QCD entirely. Because of that remarkable feature, the light-by-light scattering was the first physics process to apply both BFKL [1,2]. and the approaches [3]- [9] based on development of BFKL. Another attraction to study this process is that this process in studied experimentally by ATLAS Collaboration [10].
The external photons in the process (1) can interact through quark loops only. First of all, it can be just a single quark loop without QCD radiative corrections. Throughout the paper we address this case as the "Born" approximation. This approximation was applied to analysis of the ATLAS results in Refs. [11][12][13].
The Born case can be complemented by insertion of gluon propagators into the same single quark loop. An example of the involved Feynman graphs is shown in Fig. 1. The leading contributions in this case are double-logarithmic (DL). In the simplest version, series of DL contributions to a scattering amplitude A looks as follows: where c DL n are numerical factors and A Born denotes amplitude A in the Born approximation. When the resummation of the series (3) is done, amplitude A is calculated in DL approximation (DLA). Total resummation of DL contributions of the graphs involving the single quark loop (see Fig. 1) to process (1) was done in Ref. [14] for the case of the forward kinematics at t = 0. Calculations in Ref. [15] confirmed the results of Ref. [14] and generalized them on the case of t = 0, and also accounted for the running coupling effects. Further progress in studying the process (1) involves accounting for the graphs where the initial t-channel photons γ * (p), γ * (p ′ ) and the final photons γ * (q), γ * (q ′ ) are attached to different quark loops. These loops are related by an arbitrary intermittence of the t-channel quark and gluon ladder rungs complemented by non-ladder contributions. An example of the involved graphs is shown in Fig. 2, where non-ladder contributions are not shown for simplicity. The u-channel ladder graph with replacement p ⇆ p ′ is also not shown. Such graphs yield DL contributions, however they also bring the other kind of contributions called leading logarithms (LL) where the single-logarithmic contributions are multiplied by the power of s: where c LL n denote numerical factors. It was shown in Ref. [16] that DL contributions in Eq. (4) are absent because the DL contributions from involved Feynman graphs cancel each other. Values of c LL n in high orders in α s are mostly determined by the gluon ladder-like constructions connecting the upper and lowest quark blobs (they are called the impact-factors). The LO BFKL equation [1] sums the leading logarithms to all orders in α s . Solution to this equation is found in the form of the infinite sum of the high-energy asymptotics, each ∼ s ωn . The asymptotic with maximal ω n (the leading asymptotics) is called the BFKL Pomeron. Impact of total resummation of sub-leading contributions ∼ sα s (α s ln s) n on the Pomeron intercept was calculated in Ref. [2].
The fact that DL contributions multiplied by s cancel each other means that A Born in the DL series of Eq. (3) cannot contain contributions ∼ s in contrast to the series (4). By this reason, there is the conventional opinion in the literature that the LL contribution (4) surpasses the DL contribution (3) to such an extent that one can neglect DL contributions (3) compared to (4). In the present paper we prove that this opinion is false. We demonstrate that the high-energy asymptotics of the light-by-light amplitudes in DLA are of the same order as the BFKL result. To this end we calculate the amplitude A γγ of the process (1) in DLA at zero and non-zero values of t. By doing so, we account for both quark and gluon DL contributions and sum them to all orders in α s . First we keep α s fixed and then account for the running coupling effects. Once A γγ is known, applying the Saddle-Point method to A γγ allows us to calculate its asymptotics. As this asymptotics bears the vacuum quantum numbers, it is a new contribution to Pomeron. Throughout the paper we address it as the DL Pomeron. We compare it to the BFKL Pomeron. As we know both the amplitude A γγ and its asymptotics, we can estimate the minimal value of s, where the DL Pomeron, can reliably represent its parent amplitude A γγ .
In order to perform calculations in DLA, we compose and solve Infra-Red Evolution Equations 1 (IREE) forÂ γγ . This approach was suggested by L.N. Lipatov in Refs. [18,19] and generalized on the case of inelastic processes in Ref. [20]. It is based on factorization of partons with minimal transverse momenta first proved by V.N. Gribov [17] in QED and then extended to the non-Abelian field theories. The IREE approach has proved to be an efficient instrument for total resummation of DL contributions to various processes in QED,QCD and other theories. More info on this approach can be found in the overviews [21]. In the context of the light-by-light scattering the IREE approach was used in Ref. [14,15].
Essence of this approach is first to neglect masses of involved quarks so as treat quarks and gluons the same way. Then one should introduce a cut-off µ to regulate IR divergences. It is convenient to introduce such cut-off µ in the transverse momentum space so that transverse momenta of all virtual partons were greater than µ. In order to neglect the involved quark masses m q and guarantee applicability of the Perturbative QCD, µ should satisfy the obvious restrictions: Otherwise a value of µ is not fixed, so one can study evolution with respect to µ and construct thereby evolution equations. All IREEs look much simple when written in terms of the Mellin amplitudes F γγ . We use the Mellin transform in the standard form: where we have introduced the µ-dependent logarithmic variables ρ, y 1 , y 2 : Calculation of A γγ in DLA leads to different results, A we provide A γγ with the superscript M: γγ corresponds to opposite case where We call it the case of Deeply-Virtual photons. In the present paper we calculate both A γγ . Throughout the paper we will use the generic notation A γγ in the cases insensitive to the difference between A , expressing them in terms of the auxiliary amplitudes A γq and A γg (the superscripts q, g refer to quarks and gluon respectively). Amplitudes A γq , A γg are expressed through parton-parton amplitudes in Sect. III. In Sect. IV we obtain explicit expressions for the parton-parton amplitudes and we arrive at explicit expressions for A (M,D) γγ . In Sects. II-IV we considered the particular case of the forward kinematics assuming that t = 0. In Sect. V we study the photon-photon scattering amplitudes in the forward kinematics with non-zero t and arbitrary relations between t and Q 2 1,2 . In Sect.VI we calculate the asymptotics of A (M,D) γγ at s → ∞, applying the saddle-point method and arrive at a new Pomeron. Then we study the applicability region for the use of the asymptotics. In Sect. VII we compare the DL Pomeron to the BFKL Pomeron and discuss their similarities and differences. Finally, Sect. VIII is for concluding remarks.

II. EVOLUTION EQUATIONS FOR PHOTON-PHOTON SCATTERING AMPLITUDES
In this Sect. we study the process (1) in the forward kinematics (2) at t = 0. We construct IREEs for the amplitude A γγ (s, Q 2 1 , Q 2 2 ). They represent A γγ (s, Q 2 1 , Q 2 2 ) through the amplitudes of photon-parton scattering. For constructing the IREEs, we use that the t-channel two-parton state factorizes the graph in Fig. 2 into two parts providing the transverse momentum (with the plane formed by p and q) of the pair is minimal compared to other transverse momenta. Such partons are called the softest. Technically, constructing IREEs for A γγ is pretty similar to the derivation of the IREEs considered in Ref. [15] with one important exception: The softest partons in Ref. [15] were quarks while now they can be both quarks and gluons. Combining the both these opportunities, we come to the IREE depicted in Fig. 3. Applying the standard Feynman rules to the graphs in Fig. 3 and using the Mellin transform (6), we write the IREE in the analytical form in the ω-space: We have used in Eq. (10) the logarithmic variables introduced in Eq. (7). The l.h.s. of Eq. (10) corresponds to differentiation of Eq. (6). It includes the derivatives with respect to y 1,2 and the result of differentiation of the factor (s/µ 2 ) ω . The r.h.s. of Eq. (10) involves the auxiliary photon-parton amplitudes so that amplitude A γg corresponds in the t-channel to the two-photon production of a quark-anti-quark pair while amplitude A γg describes production of two gluons A γg depends on µ through variables ρ, y 1 , y 2 , Technology of solving IREEs for objects with several µ-dependent variables requires that all such variables should be ordered. We use the ordering of Eq. (2), complementing it by the restriction Q 2 1 ≫ Q 2 2 and arriving thereby at but we will represent the final expressions for A γγ (s, Q 2 1 , Q 2 2 ) in the form, where Q 2 1 and Q 2 2 are not ordered. Amplitude A γγ in the IREE (10) depends on three variables. The conventional way to solve such many-parametrical IREEs is firstly to solve IREEs with the lesser number of variables and then to use such solutions for specifying the general solution of (10).

A. All photons are nearly on-shell
We start with calculation of A γγ in the simplest kinematics where Q 2 1 ≈ Q 2 2 µ 2 , i.e. when We denote its f γγ (ω) the Mellin amplitude for the kinematics (14). The l.h.s. of the IREE for f γγ (ω) corresponds to the l.h.s. of Eq. (10) where the y 1,2 -dependence is absent while the r.h.s. coincides with the one of Eq. (10). Therefore we obtain where amplitudes f γq (ω) and f qγ (ω) correspond to the processes (11) and the reversed process respectively. Similarly, f γg (ω) and f gγ (ω) correspond to the process (12) and the reversed process respectively.

B. One of the photons is on-shell and the other is off-shell
Let us consider the more complicated case when i.e. Q 2 1 ≫ Q 2 2 ∼ µ 2 . We denote A γγ (s, y 1 ) the amplitude corresponding to that case and keep the notation F γγ (ω, y 1 ) for the Mellin amplitude of this process. Again, the r.h.s. of the new IREE coincides with the r.h.s. of Eq. (10). The l.h.s. corresponds to the l.h.s. of Eq. (10) without the second derivative because y 2 = 0: Amplitudes F γq (ω, y 1 ), F γg (ω, y 1 ) and f gγ (ω), f gγ (ω) are supposed to be calculated independently. Once they are known, the general solution to Eq. (17) is In order to specify an unknown function C 2 in Eq. (18) we use the matching: where f γγ (ω) is defined in Eq. (15). Therefore, F γγ (ω, y 1 ) in kinematics (16) is represented in terms of the photonquark amplitudes: where r = q, g so that F γr (ω, y) stands for the photon-quark and photon-gluon amplitudes at y = 0.
The IREE for the amplitude F (M) γγ (ω, y 1 , y 2 ) in the kinematic region (21) is In order to use the symmetry with respect to y 1 , y 2 of the differential operator in (22) and simplify the IREE, we have introduced new variables ξ, η: Eq. (22) in terms of ξ, η takes a simpler form: A general solution to Eq. (24) is with C(ω, η) being an arbitrary function and the variables y ′ 1 , y ′ 2 are defined as follows: In order to specify C(ω, η), we use the matching of F (M) γγ (ω, y 1 , y 2 ) with an amplitude F γγ (ω, y 1 ) of the same process but in the simpler kinematic regime (16) considered above: Combining Eqs. (27), (25) and (20), we arrive at the following expression for F γγ : Eq. (28) is obtained under assumption that y 1 > y 2 (see Eq. (13)). Replacing in Eq. (28) η by |η|, allows us to describe the cases y 1 > y 2 and y 1 < y 2 at the same time: Substituting Eq. (29) in (6), we arrive at the expression for the amplitude A (M) γγ at moderate virtualities Q 2 1,2 :

D. Deeply-virtual photons
When Q 2 1 >> µ 2 and Q 2 2 >> µ 2 , and their product is also great, Q 2 1 Q 2 2 >> sµ 2 , the inequality in Eq. (21) is replaced by the opposite one ρ < y 2 + y 1 . (31) We address such photons as deeply-virtual ones. The principal difference between this case and the case of moderately-virtual photons is that the scattering amplitude A (D) γγ (ρ, y 1 , y 2 ) in the kinematics (31) does not depend on µ, so the IREE for it is very simple: A general solution to Eq. (32) can be written in different ways. The most convenient way for us is In order to specify A (D) γγ we use the matching with the amplitude A (M) γγ of the same process but in the region (21). It means that Replacing γγ (y 1 , y 2 ) immediately allows us to express A γγ in the whole the region ρ ≤ y 1 + y 2 :

III. AMPLITUDES OF PHOTON-PARTON SCATTERING
A. Photon-parton amplitude with on-shell photon We consider the case when y = 0 and denote f γq (ω) and f γg (ω) the Mellin amplitudes of processes (11) and (12) respectively. In the ω-space the system of IREE is where a γq = e 2 , so that a γq /ω is the Born value of amplitude f γq (ω). There is no a similar term in the equation for f γg (ω). We have used the following convenient notations in the rhs of Eq. (37): with r, r ′ = q, g and f rr ′ being the parton-parton amplitudes. The solution to Eq. (37) is with K(ω) being the determinant of the system (37): Amplitudes f rγ are given by expressions similar to Eq. (39). They can be obtained from (39) with replacement a γq by a qγ = −e 2 , which means that

B. Off-shell photons
When the photons are off-shell, the system is more involved: and a general solution to it is where and We specify the factors C 1,2 (ω) by the matching with the solution (39) of Eq. (37): with f γq and f γg being calculated in Eq. (39). This matching leads us to the following expressions:

IV. AMPLITUDES OF PARTON-PARTON SCATTERING
Amplitudes h rr ′ satisfy the system of the algebraic IREEs: where the inhomogeneous terms b ik consist of two contributions: with a ik being the Born amplitudes divided by the factor (8π 2 ). They were calculated in Refs. [23,24] : where A and A ′ stand for the running QCD couplings: with η = ln µ 2 /Λ 2 QCD and b being the first coefficient of the Gell-Mann-Low function. When the running effects for the QCD coupling are neglected, A(ω) and A ′ (ω) are replaced by α s . The terms V rr ′ are represented in a similar albeit more involved way (see Ref. [24] for detail): with and Let us note that D = 0 when the running coupling effects are neglected. It corresponds the total compensation of DL contributions of non-ladder Feynman graphs to scattering amplitudes. When α s is running, such compensation is only partial. Solution to Eq. (48) is Eq. (55), with Z given by Eq. (56), presents the explicit expressions for the amplitudes h rr ′ . Substituting them in Eqs. (44,47) and then in Eqs. (39,43), we obtain explicit expressions for the photon-parton amplitudes. Finally, substituting F γr and F rγ in Eqs. (20,29), we arrive at explicit expressions for the photon-photon amplitudes A In this Sect. we generalize the expressions for A γγ on the case of the forward kinematics (2) where t = 0. Technically, this generalization is pretty similar to what we have done in Ref. [15]. To show it, let us consider Fig. 3. We start by considering it at t = 0. The r.h.s. of the IREE in Fig. 3 involves two convolutions, where the vertical lines stand for propagators of the t-channel partons. At t = 0 each of the propagators is ∼ 1/k 2 and dealing with polarizations brings the factor k 2 ⊥ . After integrating over the longitudinal components of k has been performed, the propagators yield the factor with k ⊥ being the component of k transverse to the plane formed by the external momenta. After cancelling k 2 ⊥ in Eq. (58), we are left with the factor 1/k 2 ⊥ which yields a logarithmic contribution when integrated over k ⊥ . We remind here that, by definition of µ, the integration region is When t = 0, one of the propagators on each graph in Fig. 3 remains to be ∼ 1/k 2 while the other becomes ∼ 1/(q − k) 2 , where momentum q is related to t: −t = q 2 . After integrating over the longitudinal component of k has been done, the propagators are . (60) After cancellation of the factors k 2 ⊥ , the remaining factor 1/(|t|+k 2 ⊥ ) in Eq. (60) can yield a logarithmic contribution in the integration region Eq. (61) means that |t| acts as a new IR cut-off. This is true for the both quark and gluon convolutions in Fig. 3 and remains unchanged when the impact of the blobs in the convolutions is accounted for because the blobs bring logarithmic contributions only.
When t is so small that |t| ≤ µ 2 , the restriction (61) does not change the region (59), so there is no t-dependence of the scattering amplitudes and the expressions for A (M,D) γγ obtained for t = 0 can also be used at such small t. The case of large t, where −t > µ 2 , is less trivial because there is an interplay between parameters Q 2 1,2 and t. New photon-photon and photon-parton amplitudes in this case explicitly depend on t. We denote them M γγ and M γr (r = q, g) respectively. In this case |t| really plays the role of a new IR cut-off, so the amplitudes M γγ , M γr do not depend on µ at all. All expressions obtained in the previous Sects. remains true after replacement µ 2 by |t|. It leads to the following redefinitions of the variables ρ, y 1 , y 2 : After this redefinition, the amplitudes A γγ (s, y 1 , y 2 ) can easily be used to obtain explicit expressions for M γγ (s, t, y 1 , y 2 ) . Let us do it, considering step by step all possible relations between t and Q 2 1,2 .
A. Case of Q 2 1 ≫ |t| ≫ Q 2 2 or vice versa In this case the photon-photon amplitude M γγ does not depend on Q 2 2 under the DL accuracy, and is given by the following expression: where A γγ (ρ, y 1 ) is the amplitude is given by Eqs. (6,20).
In the opposite case of Deeply-Virtual photons obeying the photon-photon amplitude M In this case the photon-photon amplitude is insensitive to Q 2 1,2 under the DL accuracy, so amplitude M where amplitude A (0) γγ (ρ) is given by Eqs. (6,15). Replacing s by u in the expressions considered in the cases above, we obtain amplitudes M ′ γγ . Adding them to M γγ , we eventually arrive at amplitudesM γγ describing the photon-photon scattering in the forward kinematics at t = 0.

VI. ASYMPTOTICS OF THE LIGHT-BY-LIGHT AMPLITUDES
When s → ∞ and Q 2 1,2 are fixed, the amplitudes A To begin with, we represent A γγ in the form Now one can expand the exponent in Eq. (69) in the power series at ω = ω 0 , retaining the first three terms: and solve the stationary point equation The rightmost root of this equation, ω 0 is the leading singularity. After that A (M) γγ can be represented as follows: Let us notice in advance that it follows from Eqs. (77,44,55) that λ 1 (ω 0 ) > λ 2 (ω 0 ) and therefore λ 2 (ω 0 ) can be dropped in asymptotic expressions. Eq. (72) can be written as follows: In Eq. (74) ρ → ∞, so it must be equated by another large negative contribution. The only option is that such term corresponds to a singularity of F ′ . Writing where we have denoted and using the explicit expression for F , we see that the rightmost singularity is W = 0, so that the leading singularity is the largest root of the equation In vicinity of ω 0 Eq. (74) is The main contribution to Ψ ′′ comes from differentiation of the factor 1/W in Eq. (75), so Substituting it in Eq. (73), we obtain the explicit expression for the asymptotics A as γγ of amplitude A γγ : When s → ∞ while Q 2 1,2 are fixed, A γγ → A as γγ , where Eq. (81) manifests that the asymptotics of A (M,D) γγ are of the Regge form, with the exponential being a Reggeon. As its quantum numbers in the t-channel are vacuum ones, it is a new contribution to Pomeron. we will address it as DL Pomeron. It universally (save the factors F (ω 0 ) and σ) contributes to high-energy asymptotics of all QCD processes where vacuum quantum numbers are allowed. Eq. (77) was first obtained in the context of the small-x asymptotics of the DIS structure function F 1 in Ref. [25]. Let us discuss Eq. (77) in more detail. First of all, we notice that Eq. (77) can be solved analytically only under the approximation of fixed α s . In this case the rightmost solution of Eq. (77) is In order to fix α s in Eq. (82), we use the value α s = 0.24 obtained in Ref. ([26]). When the quark contribution to ω 0 is neglected, we obtain Accounting for quark contribution at n f = 4 diminishes value of ω 0 . In this case we obtain However, the approximation of fixed α s inevitably involves a procedure of setting α s in the final expressions, which is far from being rigorous (see e.g. Refs. [26,27]). By this reason, the fixed α s approximation can be regarded as too rough to be conclusive. In the more realistic case where α s runs in every vertex of each involved Feynman graph, Eq. (77) can be solved only numerically, which was done in Ref. [25]. When the quark contribution are neglected and α s runs we obtain At last, accounting for both gluon and quark contribution together with accounting for the running coupling effects, we arrive at Eqs. (83-86) demonstrate that DL Pomeron is always supercritical though the intercept is decreasing when accuracy of the calculations increases. The Regge asymptotics always look much simpler than the parent amplitudes (e.g. cf. A as γγ to A γγ ) and by this reason they are often used even at the energies pretty far from asymptotic. On the other hand, it is not always possible to outline applicability regions where the asymptotics can be used. Fortunately, we have such possibility. In order to fix the applicability region for the asymptotics (81), we construct R as defined as the ratio of the asymptotics versus the parent amplitude: and study dependence of R as on s at fixed Q 2 1,2 . Obviously, the asymptotics A as γγ reliably represents the parent amplitude A γγ when R as is not far from 1. In Ref. [25] we showed that R as grows with s, so we suggest that the asymptotics can be used at energies s > s min where s min is the energy where R as ≈ 0.9. Numerical estimates show that At s ≤ s min the parent amplitudes (30,36) should be used instead of their asymptotics (81).

VII. COMPARISON TO THE BFKL APPROACH
The graphs contributing to A γγ are represented in Fig. 4 in the same form as used in the BFKL: the upper and lowest blobs, Φ A (p, p ′ , k A , k ′ A ) and Φ B (k B , k ′ B , q, q ′ ), are connected by the middle blob corresponding to the gluongluon scattering amplitude A gg (k A , k ′ A , k B , k ′ B ). The blobs Φ A,B are called impact-factors. By doing so, we bring A γγ to the form of BFKL-like convolution: The representation (89) dates back to the phenomenological Regge theory where impact-factors depended on specific of the external particles (photons in our case) whereas the middle blob (the phenomenological Pomeron) was universal for all processes with vacuum quantum numbers in the t-channel. This representation holds for the both BFKL and DL Pomerons. However, there is a technical difference: the impact-factors in the BFKL concept are calculated independently of the BFKL Pomeron whereas in DLA they and the middle blob are calculated universally by constructing and solving IREEs. In both DLA and BFKL the impact-factors do not have much influence on the high-energy asymptotic behavior of amplitudes A γγ . The major role here is played by the middle blob n Eq. (89): the gluon-gluon amplitude A gg is the object fastest growing with s. It makes the high-energy asymptotics of A γγ be of the Regge form, so its asymptotics is either BFKL or DL Pomeron, depending on the approach.
There is another similarity between the BFKL and DL Pomerons: The NLO BFKL Pomeron intercepts contains α s at an unknown scale. Setting of the scale was done in Ref. [27]. In contrast, the DL Pomeron intercept does not contain α s at all but contains the IR cut-off µ, which involves a procedure of specifying µ. Such specifying is universal for singlet and non-singlet Reggeons calculated in DLA (see for detail Ref. [21]).
Despite that ω f ix g 0 is close the LO BFKL Pomeron intercept, and ω 0 is pretty close to the NLO BFKL Pomeron intercept, we do not see any theoretical reason for such similarity and think that it is just a coincidence. These approaches account for totally different logarithmic contributions: BFKL deals with single-logarithmic contributions and DLA sums DL contributions.
One more difference between BFKL and our approach is that the parameter s 0 in the BFKL factor (s/s 0 ) ω0 also requires a special setting procedure which is even more complicated than setting of the α s scale (see Ref. [28] and references therein) because it includes dealing with the impact-factors involving many-parton t-channel states. This problem is absent in our approach. For instance, the factor (Q 2 1 Q 2 2 ) 1/2 in Eq. (81) is rigorously fixed by the saddle-point method.
Studying the s-behavior of the ratio R of the DL Pomeron to the parent amplitude allowed us to estimate in Eq. (88) the applicability region of the DL Pomeron. A similar estimate for the applicability region of the BFKL Pomeron has not been done.
The LO BFKL equation [1] and the approaches based on it sum LL terms ∼ sc LO n (α s ln s) n (see Eq. 4). The NLO BFKL equation [2] and the approaches based on it in addition to the LL terms account for the sub-leading terms ∼ sc N LO n α s (α s ln s) n : where c LO k and c N LO k are numerical factors for the leading and sub-leading contributions respectively. Judging by the involved parameters, the NLO contribution should be much less than the LO one. Indeed, s(α s ln s) n ≫ sα s (α s ln s) n .
The estimate (91) predicts that the NLO terms can bring a small impact on the intercept of the BFKL Pomeron. However on the contrary, the NLO contribution turned out to be of the same order as the LO one. We explain the failure of the prediction by the fact that the comparison (91) deals with the involved parameters like s, α s but ignores the numerical factors c LO n and c N LO n because it is impossible to estimate them before the calculations have been done. It drives us to conclude that the by-parameters-considerations like (91) are too rough to predict reliable estimates of numerical characteristics like intercepts. A similar situation occur with estimating the impact of the DL contributions (3): if one considers the involved parameters only, the DL terms look small compared to the LL ones but as a matter of fact the DL impact on the Pomeron intercept is quite essential.

VIII. SUMMARY AND OUTLOOK
We have obtained explicit expressions for the elastic γγ-scattering amplitudes A γγ in DLA and at the same time accounted for the running α s effects. Amplitudes A γγ were calculated first in the collinear kinematics with t = 0, see Eqs. (30, 36), then expressions (65, 67) were obtained for the light-by-light scattering amplitudes M γγ in the forward kinematics with non-zero t and for arbitrary relations between t and the photon virtualities Q 2 1,2 . For A γγ and M γγ we considered separately the cases of deeply-virtual and moderately-virtual external photons.
Applying the saddle-point method to A γγ , we obtained its the high-energy asymptotics A as γγ in Eq. (81) which is the new, DL Pomeron. This Pomeron is supercritical like the BFKL Pomeron. The value of the DL Pomeron intercept depends on accuracy of the calculations, decreasing with the increase of accuracy, as shown in Eqs. (83-86). The maximal intercept (83) corresponds to the roughest approximation where α s is fixed and quark contributions are neglected and the minimal intercept is obtained when the both running α s effects and quark contributions are accounted for. It drives us to conclude that the further increase of accuracy may diminish intercept of the DL Pomeron from supercritical values down to unity, which would agree with the Froissart bound and thereby could prevent violation of Unitarity.
Expressions for the Regge asymptotics are much simpler than for their parent amplitudes and because of that the asymptotics have often been applied instead of their parent amplitudes to description of high-energy processes. However it can be done only if the applicability regions for such approximation are known. To this end we estimated in Eq. (88) the minimal value of s, where the asymptotic A as γγ reliably represent its parent amplitude A γγ . At s ≤ s min one should use A γγ instead of A as γγ . Our results do not mean that the DL Pomeron is supposed to substitute BFKL Pomeron. On the contrary, it would be challenging to study interference between contributions of the both Pomerons to high-energy reactions. Finally, we would like to notice that the explicit expressions for A γγ obtained in our paper should be used at available energies rather than their asymptotics.

IX. ACKNOWLEDGEMENT
We are grateful to V.T. Kim, A.V. Kotikov, D.A. Ross, W. Schafer and especially to D.Yu. Ivanov for interesting discussions.