Taming of preasymptotic small x evolution within resummation framework

It is well understood that the leading logarithmic approximation for the amplitudes of high energy processes is insufficient and that the next-to-leading logarithmic effects are very large and lead to instability of the solution. The resummation at low x, which includes kinematical constraints and other corrections leads to stable result. Using previously established resummation procedure we study in detail the preasymptotic effects which occur in the solution to the resummed BFKL equation when the energy is not very large. We find that in addition to the well known reduction of the intercept, which governs the energy dependence of the gluon Green’s function, resummation leads to the delay of the onset of its small x growth. Moreover the gluon Green’s function develops a dip or a plateau in wide range of rapidities, which increases for large scales. The preasymptotic region in the gluon Green’s function extends to about 8 units in rapidity for the transverse scales of the order of 30–100 GeV. To visualize the expected behavior of physical processes with two equal hard scales we calculate the cross section of the process γ∗+γ∗→X\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ^{*}+\gamma ^{*}\rightarrow X$$\end{document} to be probed at future very high-energy electron-positron colliders. We find that at γ∗γ∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma ^*\gamma ^*$$\end{document} energies below 100GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$100 \; \mathrm{{GeV}}$$\end{document} the BFKL Pomeron leads to smaller value of the cross section than the Born approximation, and only starts to dominate at energies about 100GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$100 \; \mathrm{{GeV}}$$\end{document}. This pattern is significantly different from the one which we find using LLx approximation. We also analyze the transverse momentum contributions to the cross section for different virtualities of the photons and find that the dominant contributions to the integral over the transverse momenta comes from lower values than the the external scales in the process under consideration.


Introduction and motivation
One of the outstanding and interesting problems in the domain of strong interactions is the behavior of the cross sections of hard processes in QCD in the high energy limit. High energy Deep Inelastic Scattering experiment at HERA established very strong growth of the structure functions with decreasing values of Bjorken x [1,2]. This growth was subsequently understood as a result of the increase of the gluon density at low x due to the non-abelian gluon splitting. The increase of the parton densities towards low x is also responsible for the growth of the various cross sections of hard processes at the LHC. The standard approach to the description of the evolution of parton densities is provided by the collinear picture and the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations [3][4][5][6][7]. These equations have been derived under the assumption of fixed x and very large scale Q 2 , such that they sum all the powers of (α s ) m (log(Q 2 / 2 )) n , where is some hadronic scale. The validity of the DGLAP equations is well established and the collinear parton densities are at the core of most of the perturbative predictions at the Large Hadron Collider.
On the other hand, the Balitsky-Fadin-Kuraev-Lipatov (BFKL) [8][9][10][11][12][13] equation accounts for the evolution in x at small values of x. Within this approximation, one assumes large center of mass energy s, which is much larger than any of the scales in the process, including the momentum transfer t and hard scale Q 2 . BFKL hard-Pomeron approximation is suitably defined for the processes which are characterized by the presence of two comparable, perturbative, scales, with a large rapidity distance in between. For such processes the evolution to infrared pole, to the non-perturbative regime where pQCD becomes doubtful, is kinematically sup-pressed. 1 To give few examples, such configurations can be realized in the production of Muller-Navelet jets in hadronhadron scattering [14], forward jet production in electronhadron DIS (where jet transverse momentum is comparable to photon virtuality) [15,16] and in the γ * γ * scattering in e + e − collisions [17,18]. In these processes the DGLAP configurations, which are strongly ordered in transverse momenta, are suppressed since virtualities in the projectile and the target fragmentation regions are comparable.
It is thus expected that the BFKL evolution, for which the transverse momenta of gluons exchanged in the t channel are not ordered, would be the dominant mechanism for the growth of the cross section for such process with the energy. The LLx (leading-logarithmic) BFKL equation resums terms (α s ln s) n which in turn leads to the power like increase of the cross sections with the center-of-mass energy √ s. The BFKL Pomeron in the leading logarithmic approximation by definition violates energy-momentum conservation. A related effect was observed in QED in the process of production of e + e − pairs by Lipatov and Kuraev [19] (for further discussion, see [20]). This results in NLLx corrections being very large and negative. Thus LLx approximation is unstable and the urgent question arises on the necessity to guarantee stabilization of BFKL Pomeron approximation. BFKL Pomeron is currently available up to next-to-leading logarithmic accuracy in QCD [21,22] and up to next-to-nextto-logarithmic accuracy in N = 4 supersymmetric Yang-Mills theory (sYM) [23][24][25]. The corrections at the NLLx level turned out to be very large, reducing the value of the intercept significantly and also leading to the instabilities like the oscillating cross sections including negative values.
In this paper we shall focus on the Ciafaloni-Colferai-Salam-Stasto (CCSS) resummation [33][34][35][36] and its application to phenomenology. The main sources of the large NLLx corrections to BFKL were identified as originating from the so-called double transverse logarithms, or in Mellin space representation of the kernel, from the triple collinear poles [29]. Large negative terms are also coming from the parts of the leading order DGLAP splitting function which appear at NLLx BFKL, and they manifest themselves as the double collinear poles. The third source of the large corrections at NLLx BFKL is the running of the strong coupling. The strong coupling is fixed in LLx BFKL and practically only starts to run at NLLx level [21]. This leads to the large correction, and the increased importance of the infrared regime. Even if the BFKL evolution is considered for the one scale process, like Mueller-Navelet jets, with two equal, large, perturbative scales, the diffusion of the transverse momenta along the ladder leads to the increased importance of the infrared regime, where the coupling becomes very large. 2 The CCSS resummation procedure was based on combining the LLx and NLLx BFKL together with kinematical constraint and DGLAP anomalous dimension. Appropriate subtractions needed to be performed in order to avoid the double counting and guarantee that the momentum sum rule was satisfied. The resummation procedure was shown to stabilize the result for the BFKL gluon Green's function and significantly reduced the value of the intercept from the LLx result for wide range of values of α s .
Additional important feature of resummation, is the delay in the onset of power like increase, typical of BFKL resummation to lower values of x or consequently higher values of rapidity. This was first discussed in the context of the 'dip' of the splitting function with resummation [35]. It was demonstrated that the resummed splitting function exhibits a decrease in x followed by the increase towards asymptotia. It was understood as the interplay between different terms in the perturbative expansion, notably between the terms which appear at NLLx BFKL (or NNLO DGLAP splitting function) and the terms which appear at LLx BFKL (and would appear at N 3 LO DGLAP). In the region of moderate values of x the splitting function had a small dip, and thus could well be approximated by a flat function.
This may partly explain the success of DGLAP in the description of the current HERA data on structure functions [60] and the fact that NLO fits seem to perform better than NNLO fits, since the x P gg splitting function at NLO is flat as a function of x (there is also effect coming from the F L structure function). The region in which the BFKL growth starts to dominate is pushed towards smaller values of x, leading to relatively large flat region in x P gg function. The question thus remains, as to what is the domain in which the BFKL resummation becomes dominant and DGLAP dynamics ceases to become sufficient as a description of the high energy scattering.
The goal of this paper is to investigate in more detail the preasymptotic regime in the resummed BFKL approach in the context of the one-scale process, like γ * γ * scattering. In particular we want to establish quantitatively the boundary in plane of rapidity and a hard scale at which the BFKL increase starts to occur and leads to the sharp of the cross section which is the signature of the BFKL regime. We shall base our analysis on the CCSS resummation framework. First, we shall recall the solutions to the gluon Green's functions G(Y ; k 0 , k 0 ) using CCSS resummation and perform the study on the level of the BFKL solution itself. We find that the length of the plateau in rapidity depends almost linearly on the logarithm of the scales. The BFKL growth dominates only for rather large rapidity intervals, starting from about 4 units of rapidity for k 0 = 10 GeV to about 8 units of rapidity for k 0 = 100 GeV. Second, for the more realistic case, we then apply this solution to the process with comparable hard external scales. Such processes were argued to be the gold-plated signatures of the BFKL dynamics, since in such cases the DGLAP, or collinear configurations are suppressed. For the purpose of this analysis, we consider γ * γ * scattering [17,18,27,61,62]. Such process could be realized in the future e + e − collider, such as FCC-ee [63,64], or CLIC [65], with doubly tagged electrons in the final state. We compute this process using the BFKL resummed solution and quantify the length of the preasymptotic regime. We find that the cross section which includes BFKL is in fact lower than the Born cross section (i.e. just a two gluon exchange) for large region of energy (this was also observed in [27] which used BFKL with kinematical constraint) which increases with the increase of the scales of the virtual photons. Still, we find that an enhancement over the Born cross section of the order of factor 2 to 3 can be attained for energies of W 300 GeV provided the scales of the virtual photons are less than 10 GeV 2 . We also investigate the contributions of different transverse momenta to this process. We find that even though the impact factors are peaked at the scales of virtualities of the photons, there is significant spread in the transverse momenta along the exchanged Pomeron. This leads to the feature that the energy dependence of the resulting cross section is only weakly dependent on the virtualities of the photons.
The paper is organized as follows. In the next section, Sect. 2 we introduce the LLx and NLLx BFKL formalism, and recall the sources of large corrections in the NLLx BFKL using Mellin space representation. In Sect. 3 we recall the construction of the CCSS resummation scheme. Following that, in Sect. 4, we perform detailed numerical analysis of CCSS, focusing on the preasymptotic features of this solution. In particular, we quantify the dip or a plateau in the rapidity before the evolution starts to be dominated by the power-like growth typical of BFKL. In Sect. 5 we apply the resummed CCSS solution to the selected process of γ * γ * scattering, and analyze the dependence on the energy, polarization, scales and investigate the dominant transverse momenta which contribute to the BFKL Pomeron in this process. Finally, in the last section we briefly summarize the results.

LLx and NLLx BFKL evolution
We shall consider BFKL evolution in the forward case t = 0. The BFKL [12,13] equation can be formulated as an evolution equation differential in rapidity and it can be cast in the following generic form where G(y; k , k 0 ) is the forward BFKL gluon Green's function and K (k, k ) is the BFKL kernel which possesses the following perturbative expansion Here, K 0 is leading logarithmic (LLx) and K 1 is next-toleading logarithmic(NLLx) BFKL kernel respectively. In Eq. (1) y is the evolution variable, the rapidity, and k, k , k 0 are transverse momenta of the reggeized gluons exchanged in the t channel. In the following, we shall also use the notation |k 2 | = k 2 . The (rescaled) strong coupling constant is defined as with N c the number of colors. The kernels K 0 , K 1 can be also transformed to Mellin space, by use of the following relation where we have used the angular averaged version of the kernel. In the Mellin space the BFKL kernel has the expansion In QCD the kernels are known up to NLLx order [21,22] and the corresponding terms at the LLx and NLLx level have the following form in Mellin space and with with N f the number of (active) flavors and is the polygamma function. In N = 4 sYM theory the kernel χ is known up to NNLLx order [23][24][25]. It is also worth pointing out that in QCD the functions (k 2 ) γ are the eigenfunctions of the equation, and χ 0 (γ ) is the eigenvalue at the leading logarithmic order. Strictly speaking this is not the case in the NLLx order in QCD due to the fact that the coupling constant is running and the conformal invariance is broken. In that case, the eigenfunctions are different, and include the running of the coupling, see [66]. This however does not preclude the possibility of performing Mellin transform of the BFKL equation as in Eq. (4) and analyzing the Mellin transform of the NLLx part of the kernel in the form (7).
The above kernel in NLLx has been written for the socalled symmetric choice of scales. As discussed in [21,22,29] at NLLx level the scale choice does matter. Let us briefly recall the scale changing transformations. In general, the generic cross section for two scale process in the high energy factorization can be written in the following form Here h A (Q, k), h B (Q 0 , k 0 ) are the impact factors, for example for the jet production in the case of the Mueller-Navelet process [67] or as will be studied in this paper for the photon-gluon impact factors. In the first case Q, Q 0 would correspond to the transverse momenta of the hard jets or to the virtualities of the photons in the γ * γ * scattering. The relation for the Green's function in Mellin space and in rapidity space in this case is equal to where evolution variable for the BFKL is defined as y = ln ν kk 0 , and ν is the energy available for the BFKL evolution. In this example the energy scale choice is symmetric, but it might as well be asymmetric ν/k 2 or ν/k 2 0 . This will imply the change of the gluon Green's function as well and in consequence there are similarity transformations [22,34] for the BFKL kernel The expression Eq. (7) is given for the symmetric scale choice. In order to transform to the asymmetric scale choice, for example for s 0 = k 2 in the Mellin space one needs to add to Eq. (7), the term which is the NLLx term due to the scale changing transformation. At NNLLx the scale changing transformation is much more complicated as it involves the NLLx kernel as well [68,69]. Coming back to the NLLx kernel, the well known problem that arises at NLLx order in BFKL is due to the presence of double and triple collinear logarithms. The collinear logarithms manifest themselves as poles in the γ plane. The collinear poles, which correspond to the strong ordering when k 2 k 2 appear as 1/γ and the anti-collinear poles which correspond to k 2 k 2 as 1/(1−γ ). The double poles are arising due to the running coupling and the non-singular (in 1/z) part of the DGLAP splitting function which appear at NLLx order. To be precise, keeping most singular γ → 0 and γ → 1 contributions gives for the corresponding terms for the running coupling, for the DGLAP contribution. There are also triple collinear poles which appear due to the kinematical constraints [33,34]. Such constraints were discussed in the BFKL context as originating from the improved kinematics, and more precisely by the requirement that the exchanged momenta are dominated by the transverse components [26,27], for more recent work on different forms of kinematical constraint see [69]. These contributions generate the double transverse logarithms and in the Mellin space they exhibit most singular behavior, the triple collinear poles As it has been demonstrated in [29,30] this collinear approximation accounts for the major part of the NLLx corrections given by χ 1 .

General setup
Let us briefly remind the reader about the content of CCSS resummation scheme. General setup for the CCSS scheme was based on the analysis of poles in the Mellin space, but the final formulation and the solution to the equation was given in the momentum space. A similar idea for the resummation was formulated previously also in Ref. [28] by combining the DGLAP and BFKL evolution with the kinematical constraint. In the CCSS resummation one subtracts triple and double poles and incorporates the full DGLAP splitting function and the kinematical constraint which both resum double and triple poles respectively. In addition, more subtractions are needed to ensure the conservation of the momentum sum rule.
In the original CCSS scheme [33,34] one starts with the LLx + NLLx BFKL kernel with LO DGLAP splitting function and imposes kinematical constraint. The double and triple poles are then subtracted in order to avoid the double counting. Also, the running coupling poles are subtracted from the NLLx expression since the running coupling is later on incorporated directly in front of the LLx kernel with the scale q, which is the transverse momentum of the emitted gluon. The expression for the resummed kernel NLLx part of the kernel in the Mellin space is then Of course one could change the scale choice in the running coupling in front of the LLx kernel, which would change the term at NLLx level. The choice of q 2 is such that no other b dependent terms are present in the LLx kernel. The choice of the scale for the coupling constant in the NLLx kernel is arbitrary in the sense that any change is of NNLLx level which is beyond the control. In CCSS scheme the scale at NLLx was chosen to be the maximum of the momenta squared of the exchanged gluons, i.e. max(k 2 , k 2 ).

Momentum representation
The Mellin space expression for the CCSS resummation was then transformed back to the momentum space, and the resummed evolution equation was directly solved in the momentum space. This was done so that running coupling could be incorporated more easily in the evolution.
With the DGLAP and kinematical constraint included, the equation becomes an integral equation in both longitudinal and transverse momentum (z, k). The final resummed kernel K (z; k, k ) in the CCSS method [33,34] is the sum of three contributions: The first term in Eq. (20) is the LLx BFKL with running coupling and consistency constraint imposed onto real part of the kernel where q = k − k . The consistency or kinematical constraint appears in the form of the ( k z − k ) (k − kz) functions, and is here symmetric, which corresponds to the symmetric scale choice. To be precise, there are different versions of this constraint which appear in the literature, see [26,27,70]. It turns out that all versions are generating the same leading cubic poles in the Mellin space at NLLx, they generate no double poles, and with the difference starting to appear in the single pole level, see [69]. The scale in the running coupling corresponds to the transverse momentum of the emitted gluon q 2 . The second term in (20) is the non-singular (in z ) DGLAP splitting function terms with consistency constraint Finally, the last term is the NLLx part of the BFKL with subtractions given by Eq. (19), all transformed into momentum space The non-singular (as z → 0) splitting function in the DGLAP terms appearing in Eq. (22) is defined as follows:P where the P (0) gg is the DGLAP gluon-gluon splitting function in LO (in this work we only consider purely gluonic channel, Following arguments presented in [34] we note that the argument of the splitting functionP has to be shifted in Eq. (22) in order to reproduce the correct collinear limit when the kinematic constraint (kz < k < k z ) is included. In addition to the terms presented above, one needs also to include the terms which will cancel spurious DGLAP anomalous dimension at NLO. This was done by adding extra terms. Since there is an ambiguity in this procedure, two different schemes were proposed A,B in [33]. In here, we shall utilize scheme B, which importantly also satisfies the momentum sum rule.

Numerical results for the gluon Green's function
We start the presentation of the numerical results for the gluon Green's function G(y; k 1 , k 2 ) which is the solution to the resummed BFKL evolution equation. In this analysis we shall focus on the preasymptotic features of the resummed solution and try to quantify them. The numerical method for the solution of this equation was identical to the one used in [34]. We start from the fixed coupling case, and compare the solutions of the leading order (LLx) BFKL with the resummed evolution. The solutions are demonstrated in Fig. 1, for two values of the fixed couplingᾱ s = 0.1, 0.2. Two sets of curves are presented, one for the resummed case (blue) and one for the LLx case (magenta). In each case, we show the gluon Green's function G(y; k 1 , k 2 ) and G(y; k 1 , k 1 +δk), where δk is a small value. In our case it is equal to the distance between the points on the logarithmic grid in ln k, i.e. δk = k 1 (exp(δ ln k) − 1) where δ ln k is the spacing in the grid. The reason for that, as explained in [34], is that for exactly equal values of scales the preasymptotic behavior is dominated by the numerical 'delta' function, and thus leads to a minimum that partially depends on the numerical details, like the grid spacing. This is evident when we compare the solutions for G(y; k 1 , k 1 ) and G(y; k 1 , k 1 +δk). We observe that as expected, the resummed evolution exhibits much slower rise than the LLx evolution which is related to the lower value of the resummed intercept. However, this is not the only reason that the resummed solution leads to much smaller numerical value for the Green's function. On top of that there are substantial preasymptotic effects, which manifest themselves as the presence of the 'plateau' or even a minimum in the resummed case. Such 'plateau' is almost completely absent in the LLx case, though the preasymptotic effects are present also at this order and are visible through the slower growth in the first few units of rapidity. We observe that the width of the flat region in rapidity in resummed case strongly depends on the value of the strong coupling, and it can extend for many units of rapidity. For values of the coupling of the order ofᾱ s = 0.1, the 'plateau' In Fig. 2 we plot the gluon Green's function as a function of rapidity for the running coupling case. We observe similar plateau/dip as in the calculation with fixed coupling. In this case the position of the dip or a preasymptotic plateau and the onset of the increase depends on the value of the scales for which the gluon Green's function is evaluated. In general for larger scales, the onset of the increase is delayed to larger rapidities.
The preasymtotic plateau can be better illustrated in the two dimensional plot Fig. 3 where we show G(y; k 1 , k 1 + δk) as a two-dimensional surface, as well as contour plot in (y, log(k 1 )) space. The 'dip' in rapidity is most prominent in the low k 2 region, and it is clear that for several units of rapidity the growth is very slow.
In order to quantify the dependence of the preasymptotic region on the external scales, we take the solution with equal scales and evaluate the minimum of this function. As can be seen from Fig. 1 the position of the minimum in this case gives a good estimate of the onset of the rise in rapidity.
In Fig. 4 we show the y dip , which is the value of the rapidity for which the gluon Green's function has a minimum, as a function of scale k. We performed two calculations, where we measure the position of the minimum in G(y, k, k) and in G(y, k, k + δk). Calculations in Fig. 4 were performed with running strong coupling. We see that, to a very good approximation, the dependence on the transverse momentum is logarithmic. In fact the linear fit shown in this figure describes the extracted points very well. This suggests, that the dependence on the value of the strong coupling is like ∼ 1/ᾱ s (k). The approximately linear dependence on ∼ 1/ᾱ s (k) is demonstrated in the plot shown in Fig. 5. Such simple dependence on the value of the coupling constant suggests that the minimum occurs for particular values of the product y dip α s , The linear fit works slightly worse for the G(y; k, k + δk). In the latter case there is some curvature, visible especially when going to higher values of k. In any case, the position of the dip varies from about 3 units of rapidity for k = 2 GeV to about 7 − 9 for k = 300 GeV. This indicates that the preasymptotic effects are rather large and significantly delay the onset of the BFKL regime with Pomeron-like growth.
Thus from this numerical analysis, we see that the LLx and resummed solutions vary not only with respect to the speed of the evolution, but also with respect to where exactly the BFKL evolution sets in. This can lead to much greater differences numerically than one would naively expect if the difference was only in the intercept.
The preasymptotic features will be however modified when convoluted with the impact factors due to the initial spread of the transverse momenta. Before we analyze it in the context of the physical process, we shall test the dependence on the initial conditions, by simply modifying the initial condition to be a Gaussian in log k 2 with some width given by parameter δ, i.e.   Fig. 6 we show the solution for k j = 30 GeV for four cases: delta input, and three Gaussian inputs with different widths: δ = 0.5, 1.0, 2.0. The inputs have been normalized to unity, when integrated over log k. Interestingly all the solutions converge at high rapidity, and evolve with the same power asymptotically as they should. The dip and the slow down of evolution at initial rapidity is most prominent in the case of the delta function, and becomes less visible for broader initial conditions, i.e. when δ is increased. On the other hand the broader initial conditions lead to solutions that initially start to increase faster with rapidity. Overall, we see that the dependence on the initial conditions is almost washed out only for rapidities of about 8 units.
As a last part in this section we analyze the approach to asymptotia by investigating the logarithmic derivative of the gluon Green's function with respect to the rapidity. To simplify the analysis we study here the fixed coupling case with all other  corrections included. In Fig. 7 (left plot) we show the logarithmic derivative of the gluon Green's function as a function of rapidity y. We see that one needs to evolve to very high rapidities for the derivative to go to a limiting values for the intercept. For example, for value ofᾱ s = 0.15, one is close to asymptotic value of the intercept for rapidities larger than 20. The lower the values of α s , the higher the rapidity at which the asymptotic behavior sets in. The same quantity is plotted in a different way in Fig. 7, as a function of combined variableᾱ s y and it is rescaled by 1/ᾱ S . In the case of the leading logarithmic result this quantity should go to χ 0 ( 1 2 ) = 4 ln 2 which is indicated as a dashed horizontal line. We see first that one needs to go to rather high values ofᾱ s y in order for the intercept to settle at an asymptotic value, and moreover these asymptotic values are still away from the leading logarithmic value even at the smallest values of the strong coupling ∼ 0.025. That indicates, that  Fig. 6 The gluon Green's function G(y; k, k + δk) as a function of rapidity for k = 30 GeV. Three different initial conditions are considered: delta function (blue solid), and Gaussian initial conditions with widths equal to δ = 0.5, 1.0, 2.0 (red dashed, black dotted, magenta dashed-dotted) respectively. Initial conditions have been normalized in practice, for realistic energies one cannot expect the leading logarithmic result to be valid, even in approximate form.

γ * γ * scattering at high energy
In the previous section we have discussed the preasymptotic features of the gluon Green's function stemming from the solution to the resummed BFKL evolution within the CCSS scheme. We have identified plateau and even a dip in the solution as a function of rapidity which depends on the values of the external scales. The basic behavior is such that, as the scales become larger the preasymptotic plateau extends to the larger values of rapidities and it is only at several units of rapidity when the onset of BFKL growth takes place. However, the results of the previous section are rather academic since they do not really tell us anything about the preasymptotic region in the physical processes. One could expect from the considerations of the previous section that the onset of the BFKL regime will be significantly delayed, perhaps by several units of rapidity. However, the exact position of the minimum and the extent of the plateau will crucially depend on the details of the initial conditions for the solutions, see Fig. 6. In fact, such features may vary from process to process, given that in the physical reaction the gluon Green's function is convoluted with the impact factor, that allows the coupling of the Pomerons to hadrons on a given process. Even if the impact factor is peaked at some value of the transverse momenta, there will be a significant spread in these momenta. In order to make more realistic predictions, we shall now investigate a process of the scattering of two virtual photons, producing light and heavy quarks. The gluon Green's function needs to be convoluted with the impact factors, which will provide for the suitable initial conditions for the evolution. The predictions will be made for the energies which potentially could be achieved in the future high energy e + e − machines, for example like planned FCC-ee [63,64] or CLIC [65,71]. The diagram for the process of the γ * γ * scattering is indicated in Fig. 8.
This process was calculated in number of works [17,18,28,61,62]. In the first studies of this process at high energy [17,18] the BFKL equation in LLx was used. In [28] the BFKL evolution was supplemented with kinematical constraint, which partially accounts for the higher order corrections. In [61,62] studies were done at NLLx level, in [61] also including the stabilization of the high energy expansion by employing the principle of minimal sensitivity (PMS) [72,73] and Brodsky-Lepage-Mackenzie (BLM) [74] scale fixing scheme. This approach is however rather different from the approach presented here, which utilizes analysis of the collinear singularities at LLx and NLLx level.  We shall now recall basic formulae for the two-scale process of scattering of two virtual photons at very high energies. The advantage of such a process is the possibility of the control of the scales in the process and study of the cross section behavior as a function of energy for fixed virtualities or given masses of the produced quarks. In the present analysis we shall use LO impact factors, but with exact kinematics as implemented in [28,75]. This includes partially higher order corrections in the impact factors. We note that, the program for the calculation of the impact factor at NLLx was started two decades ago in pioneering works [76][77][78]. The original result was obtained however in a form which was rather difficult for practical use. The more recent calculation of NLLx impact factor was performed [79] in the coordinate representation which then was transformed to Mellin space yielded simpler result. Of course for completeness, it would be desirable not only to perform the analysis using the impact factors at NLLx but also perform appropriate resummation and matching to the BFKL gluon Green's function. The latter task is certainly important but goes beyond the scope of this study. In the following we shall analyze many features using ratios of the calculation which includes BFKL with respect to the calculation based on the two gluon exchange, where effects due to the impact factors should largely cancel out.
Let us denote by q 1 , q 2 the four-momenta of the photons. We define Q 2 1 = −q 2 1 and Q 2 2 = −q 2 2 and W 2 = (q 1 + q 2 ) 2 is the center-of-mass energy squared of the photon-photon system. The cross section for the scattering of two virtual photons due to the exchange of the hard Pomeron has the following form [75] where the limits on the integrals are given by and The indices i j denote the two possible polarizations of two virtual photons: T for transverse and L for longitudinal polarization. The mass of the quark in the quark-box is m q . The Bjorken x variable is defined as Functions jq (k 2 , Q 2 2 , ξ) denote the photon-gluon impact factors, i.e. the contributions of the coupling of the two-gluon system to the virtual photon through the box and crossed-box diagrams (at LO). Their exact expressions (see [28,75] for the transverse polarization and for the longitudinal polarization. In the above equations the following variables were defined and the denominators and the maximum value on the ρ integration is The positivity of the argument under the square root leads to the maximum value of k 2 max given by (32). The expression (31) only contains the contribution due to the exchange of the hard Pomeron, i.e. the perturbative part. In principle the cross section can contain some non-perturbative components which should be modeled. For example, the integrals in (31) have a lower cutoff k 2 0 , which can be interpreted as an non-perturbative cutoff. For the phenomenology, these contributions need to be modeled, and they can be numerically important for low energies W < 100 GeV as discussed in [75]. We shall investigate the sensitivity of the results to the cutoff parameter in detail later on in this section.
The function G i (xξ, k 2 , Q 2 ) in (31) is the solution to the BFKL equation with the initial condition given by for transverse polarization and for the longitudinal polarization. The above expressions can be obtained from functions (35) and (36) by integration and they again correspond of the coupling of two gluons to the quark box.
We have implemented the CCSS resummed solution (with contributions up to NLLx BFKL and NLO DGLAP) as presented in the previous section to calculate the γ * γ * cross section using formulae presented in Eqs. (31)- (43).
In Fig. 9 we show the cross section as a function of energy W and for the selected virtualities Q 2 1 , Q 2 2 of both photons. In all calculations we assumed the virtualities of the photons to be equal, i.e. Q 2 1 = Q 2 2 = Q 2 . We consider two cases, where the two photons have the same polarizations: transverse-transverse or longitudinal-longitudinal. There is also longitudinal-transverse case, whose magnitude is between the two cases. The calculation is done for three massless flavors and a massive charm quark with a mass of m c = 1.27 GeV. The solid lines indicate the cross section based on the resummed BFKL evolution. The dashed lines correspond to the Born calculation, i.e. two gluon exchange between impact factors, which is obtained by using (42) and (43) in (31) for G L , G T . The Born calculation therefore tends to a constant value at highest values of the energy W . The initial growth of the Born cross section at low and intermediate energies, W 100 GeV, is due to the opening of the phase space for the impact factors. We see a clear growth of the BFKL calculation with W , which leads to a much higher cross section than the Born exchange for highest energies, W > 100 GeV. However, the BFKL calculation is lower than the Born calculation for low energies, of the order of W 10 − 100 GeV. The point where the cross section with the BFKL resummation becomes larger than the Born calculation depends strongly on the virtualities of the photons, i.e. for higher virtualities it moves to a higher energies.
This effect is better illustrated in Fig. 10 for the transversetransverse case and Fig. 11 for the longitudinal-longitudinal case where the ratio of the resummed BFKL to Born calculation is shown as a function of energy W on the left plots, and as a function of x variable on the right plots. We see that the BFKL based computation is in fact much lower at low energies than the Born calculation. One could argue that perhaps      Fig. 11 The cross section ratio of resummed BFKL to Born for virtual γ * γ * scattering as a function of energy W (left) and x(right) for different virtualities of the photons Q 2 = Q 2 1 = Q 2 2 = 2, 10, 100 GeV 2 . Polarizations of both photons were taken to be longitudinal the BFKL calculation should not be completely trusted at very large values of x ∼ 0.1 and above, nevertheless the effect is prominent even at values of x which are lower than 0.01. This effect, as well as a delay in the onset of the asymptotic small x growth, was observed also in [75] which used kinematical constraint in the LLx BFKL with running coupling, and thus partially accounting for the resummation effects. In the present calculation, which uses full resummation, the effect is more pronounced. Most likely this is due to the additional non-leading effects incorporated into resummation, like the non-singular parts of the DGLAP splitting function which tend to bring the solution even lower. We see that the preasymptotic region extends to quite large energies, of the order of 50-200 GeV, or small x ∼ 5 ÷ 1 × 10 −3 , depending on the virtuality. We have also verified that the calculation based on the LLx BFKL equation (with just running coupling included) is always larger than the Born calculation, even for low energies. This is illustrated in Fig. 12 (left plot). We thus conclude that the strong suppression of the calculation based on the resummed evolution with respect to the Born calculation is indeed due to the non-leading effects and resummation. It is thus connected to the strong preasymptotic effects of the gluon Green's function analyzed in the previous section. We observe, that even though the resummed BFKL cross section is lower than the Born cross section for lower energies, the ratio still grows with energy, starting from the low energies. Only in the case of Q 2 = 100 GeV 2 , particularly for the longitudinal-longitudinal case, a visible delay in the onset of the growth is present, similar to the plateau for the gluon Green's function analyzed in the previous section. This shows that when the gluon Green's function is convoluted with the impact factor, the dependence on the energy (particularly at lowest energies) is strongly modified with respect to the 'bare' gluon Green's function.
We also note that for the case of longitudinal-longitudinal polarization the BFKL calculation dominates the Born approximation, i.e. starting with lower energies, than in the transversetransverse case (compare Figs. 10 and 11). This is related to the difference in the energy behavior of the two polarization cases. In Fig. 12, we show the ratio The transverse-transverse contribution dominates the cross section, and the ratio is rather flat for the energies considered. We checked that, for the higher values of Q 2 the ratio starts to slightly decrease for higher energies W . Overall, the LLx contribution is about 10% of the TT contribution for the range of parameters explored. We also observe, from Figs. 9 and 10 that the exponent governing the energy growth for the transverse-transverse case is almost identical for all the values of photon virtualities. From the analysis of the previous section, and in general expectation from the solutions to the BFKL equation one would conclude that the exponent should depend strongly on the value of Q 2 , since the Pomeron intercept is a function of α s which should be dominated by the typical values of transverse momenta of the order of virtuality in the impact factor. In such scenario, one would expect that the rate of growth at Q 2 = 2 GeV 2 would be faster than at Q 2 = 100 GeV 2 , i.e. the ratio of the intercepts would be approximately ∼ 1.5 − 1.7, see Fig. 1 in [34]. For the transverse-transverse calculation, this does not seem to be the case. The situation should be slightly different for the longitudinal-longitudinal calculation where the impact factor has narrower distribution in the transverse momenta, see Fig. 11. We checked, by computing the derivative, that the growth with x for Q 2 = 100 GeV 2 is a bit slower than for Q 2 = 2 or Q 2 = 10 GeV 2 , with a difference in the exponent being very minimal, of around few percent. Overall, there does not seem to be a lot of dependence of the rate of growth with W on the scales of Q 2 for the impact factors.
To better understand the behavior of the cross section, and in particular which transverse momenta contribute to the integral, we analyzed the distribution in the transverse momenta in the expression for the cross section Eq. (31). That is we define the function which is the integrand in Eq. (31) (44) and plot the normalized distribution for which the areas under the curves are equal to unity. In Figs. 13 and 14 we show these distributions for different sets of Q 2 values. In Fig. 13 we show the calculation with the fixed value of the strong coupling, equalᾱ s = 0.2, both in impact factors and in the BFKL evolution, indicated by the solid lines. For comparison we show also the calculation for the Born case, i.e. without any evolution, just a two gluon exchange. The right plot shows the case for the longitudinal-longitudinal polarization, and the left plot is the case for the transverse-transverse polarization. First, we note that the distribution in the longitudinal-longitudinal case is narrower than in the transverse-transverse case. This is a well known effect, mainly due to the presence of the aligned jet configurations in the impact factors, i.e. when z = 0 or z = 1 in the case of the transverse polarization (see Eqs. (42), (43)). Second, the peak of the distribution is around the value of Q 2 . The exception seems to be in the case of very low Q 2 = 2 GeV 2 , where the distribution is shifted to the slightly higher values of k 2 . We note however, that the integrals over p 2 in (35), (36) are cut off at p 0 = k 0 = 1 GeV.
The calculations are strongly modified by the presence of the running coupling. We note that the running coupling is present both in the BFKL evolution as well as in the impact factors. Following [75] we considered two scenarios of the running ). The difference between the two can be treated as the scale variation in the impact factor. As evident from Fig. 14 the distribution is strongly shifted to the lower values of k 2 even for large values of Q 2 . The effect is more prominent in the BFKL case where the diffusion in k 2 provides additional shift of the distribution towards infrared regime. As a result, in the case of Q 2 = 100 GeV 2 , the peak of the distribution in k 2 is moved down to about 10 GeV 2 .
The relative importance of the contributions from the different k 2 regions is also illustrated in Figs. 15 and 16 where we show the following quantity as a function of the cutoff k 2 cut . This ratio clearly illustrates the typical transverse momenta which are contributing to the cross section. The arrows indicate the median value of the k cut . Figures 13 and 14 can be interpreted as the logarithmic derivatives over k 2 cut of results shown in Figs. 15 and 16. In the fixed coupling case the median is very closely located to the value of the Q 2 , for moderate to high values of Q 2 = 10, 100, 500 GeV 2 . In the case of the running coupling, for Q 2 = 10, 100, 500 GeV 2 , there is significant shift of the median to much lower values of k 2 , which is more pronounced for the transverse-transverse case. The exception is the case when Q 2 is very small, i.e. Q 2 = 2 GeV 2 . In this case the median seems to be shifted to higher values of k 2 , even in the case of the running coupling. This is of course a very low scale, close to the non-perturbative cutoff, and the dependence on the lower cutoff k 2 0 should be significant. To test the sensitivity of the calculations to the parameter k 2 0 , we have performed the calculations with k 2 0 = 0.5 GeV 2 and k 2 0 = 1.5 GeV 2 in addition to k 2 0 = 1 GeV 2 . We show the results in Fig. 17. Again, left is transverse-transverse and right plot is for the longitudinal-longitudinal case. For clarity, we only show three values of Q 2 = 2, 100, 500 GeV 2 . We see significant cutoff dependence, as expected for low Q 2 = 2 GeV 2 . The cutoff dependence becomes smaller for Q 2 = 100 GeV 2 and is almost negligible for higher values. It is also clear that longitudinal-longitudinal case is somewhat less dependent on the value of the non-perturbative cutoff than the transverse-transverse calculation, which is expected effect, due to the smaller width of the k T distribution. Overall, at very small Q 2 2GeV 2 the perturbative (or 'hard' Pomeron) calculation, based on the leading twist, becomes very sensitive to the cutoff on the transverse momenta and thus becomes unreliable at such low scales. In order to make the predictions for the cross section, one would need to perform phenomenological modeling of the 'soft' contribution, i.e. originating from the region of the transverse momenta below the cutoff. Such analysis was done for example in [75] where the cutoff dependence of the soft and hard contributions partially canceled in the total sum, the latter being less sensitive on the cutoff.

Summary and conclusions
Within QCD we found that the small x resummation approach which accounts for the leading and non-leading logarithmic terms to small x phenomena leads to a new portrait of perturbative Pomeron. The contribution of preasymptotic terms for one hard scale processes is relevant for the diminishing of effective intercept in a wide kinematical range and also for the suppressing of onset of BFKL Pomeron at moderate x, or equivalently at moderate energies. We showed that the onset of the BFKL growth can be delayed by 4 − 8 units of rapidity for transverse scales between 3 and 100 GeV. The analysis indicates approximately logarithmic dependence of the 'dip' (or the width of the plateau) in the gluon Green's function on the transverse momentum k. This, rather complicated, energy dependence of the amplitude indicates that in a wide preasymptotic region of x the resummed Pomeron for one hard scale processes looks as the sum of contributions with different energy behaviors. We also found that the asymptotic value of the intercept is achieved only for very large rapidities, of the order of 20 units of rapidity or larger. On top of that, we found that even for the smallest values of the coupling constant, where the resummation effects should be the smallest and one would naively expect the dominance of the leading logarithmic approximation, the preasymptotic effects are still important for very large y values and in practice the solution is never close to the leading logarithmic approximation even at the highest energies and very small values of the coupling constant. This most likely stems from the fact that the LLx approximation violates energy momentum conservation which cannot be remedied by the perturbative expansion in the strong coupling constant.
As an example of the application of the resummed calculation to the physical process we calculated cross section of the process γ * + γ * → X which is very well suited for the hunt for perturbative Pomeron. We found that the cross section at lowest γ * γ * energies is suppressed with respect to the Born calculation which is just the two gluon exchange. The suppression extends to energies of the order of W ∼ 100 GeV. Since we observe no such effect in the case of the calculation based on the LLx BFKL evolution we conclude that this rather striking effect is most likely related to the strong preasymptotic features found in the analysis of the resummed gluon Green's function. The BFKL growth can still be observed, and leads to the enhancement over the Born cross section by a factor 2 − 3 for energies above 200 GeV. This onset of the asymptotic behavior strongly depends on the values of Q 2 , and is delayed to higher energies for higher values of Q 2 . Similarly to the analysis of the gluon Green's function the cross section with the resummed BFKL solution is much lower than the calculation with LLx solution, factor of about 5 for energies of 1000 GeV, for high scales Q 2 100 GeV 2 .
We have also studied the contributions to the cross section stemming from the different regions of the transverse momenta of the gluon density. The integrand in the transverse momentum exhibits a maximum, however its position is shifted towards the infrared with respect to the external scale Q 2 . This is simply related to the running of the coupling and to the fact that the distributions in the transverse momenta are relatively broad, especially for the transverse-transverse configurations, due to the nature of the impact factors. The distribution in the transverse momentum is additionally broadened by the BFKL diffusion. For example, for scales Q 2 = 500 and 100 GeV 2 , the median transverse momentum is shifted to about 80 and 15 GeV 2 for the transverse-transverse case and to 150 and 20 GeV 2 for the longitudinal-longitudinal configuration.
We expect that the phenomena that we have discussed here are quite general and hence would manifest themselves in other processes which are described in terms of single as well as multiple vacuum exchange.