Singlet structure function $$F_1$$F1 in double-logarithmic approximation

The conventional ways to calculate the perturbative component of the DIS singlet structure function $$F_1$$F1 involve approaches based on BFKL which account for the single-logarithmic contributions accompanying the Born factor 1 / x. In contrast, we account for the double-logarithmic (DL) contributions unrelated to 1 / x and because of that they were disregarded as negligibly small. We calculate the singlet $$F_1$$F1 in the double-logarithmic approximation (DLA) and account at the same time for the running $$\alpha _s$$αs effects. We start with a total resummation of both quark and gluon DL contributions and obtain the explicit expression for $$F_1$$F1 in DLA. Then, applying the saddle-point method, we calculate the small-x asymptotics of $$F_1$$F1, which proves to be of the Regge form with the leading singularity $$\omega _0 = 1.066$$ω0=1.066. Its large value compensates for the lack of the factor 1 / x in the DLA contributions. Therefore, this Reggeon can be identified as a new Pomeron, which can be quite important for the description of all QCD processes involving the vacuum (Pomeron) exchanges at very high energies. We prove that the expression for the small-x asymptotics of $$F_1$$F1 scales: it depends on a single variable $$Q^2/x^2$$Q2/x2 only instead of x and $$Q^2$$Q2 separately. Finally, we show that the small-x asymptotics reliably represent $$F_1$$F1 at $$x \le 10^{-6}$$x≤10-6.


Introduction
A description of the singlet structure function F 1 in the framework of Collinear Factorization usually involves DGLAP [1][2][3][4] to calculate the perturbative contributions. In this case F 1 is represented in the form of two convolutions: F 1 = C q (x/y) ⊗ q(y, Q 2 ) + C g (x/y) ⊗ g(y, Q 2 ), (1) where C q and C g are the coefficient functions and q and g denote the evolved (with respect to Q 2 ) quark and gluon disa e-mail: boris.ermolaev@cern.ch tributions, respectively. These distributions are solutions to the DGLAP equations which govern the Q 2 -evolution of the initial quark and gluon distributions δq(x, μ 2 ) and δg(x, μ 2 ), evolving them from the scale μ 2 to Q 2 . Both δq and δg are defined at x ∼ 1 and Q 2 = μ 2 ∼ 1 GeV 2 . The parameter μ is also called the factorization scale. The x-dependence of F 1 is described by the coefficient functions C q,g as well as by the phenomenological factors in δq, δg. In the framework of DGLAP the evolution in the k ⊥ -space is separated from the evolution with respect to x. Such a separation takes place at x ∼ 1 only and breaks at small x as was shown in Refs. [5][6][7]. It is the theoretical reason not to use DGLAP at small x. A practical reason is that DGLAP, by its design, accounts for the total resummation of ln n Q 2 , while contributions ∼ ln n x are present in the DGLAP expressions in few first orders in α s only (through the coefficient functions in NLO, NNLO, etc.).
On the other hand, such contributions are very important at small x, so it would be appropriate to substitute the DGLAP expressions for the DIS structure functions by new ones which include the total resummation of all double-logarithmic (DL) contributions. In the first place there are DL terms ∼ (α s ln 2 (1/x)) n , then the terms ∼ (α s ln(1/x) ln Q 2 ) n , etc. Expressions accounting for resummation of DL contributions and for the running α s effects were obtained for several structure functions with nonvacuum exchanges in the t-channel: the spin structure function g 1 (the singlet and non-singlet components) and the nonsinglet component of F 1 (see the overview of [8] and the references therein). Besides, the expressions were obtained for g 1 and non-singlet F 1 combining the DGLAP results and resummation of the DL contributions, which made it possible to apply these expressions at arbitrary x and Q 2 .
However, a similar generalization of DGLAP was not obtained for the singlet F 1 . The point is that by that time F 1 in the small-x region has been intensively investigated in terms of approaches based on BFKL [9][10][11][12][13] and this looked as the only way to study F 1 at small x. Indeed, the lead-ing x-dependent contributions to F 1 proved to be the singlelogarithmic (SL) terms accompanying the "Born" factor 1/x: while the DL contributions proportional to 1/x, i.e. the terms cancel each other (i.e. c DL k = 0 for k = 1, 2, . . .) as was found first in Ref. [14]. As a result, the common strategy for investigating the QCD processes with vacuum exchanges in the t-channel was based on the use of the BFKL results. In particular, SL contributions to the structure functions F 1,2 was presented in Refs. [15][16][17]; SL contributions to F 2 in Refs. [18,19] were calculated with inclusion of the resummed anomalous dimensions in the renormalization group equation, while F 2 in Ref. [20] was calculated with direct unification of DGLAP and FFKL.
The solution to the BFKL equation is expressed through the series of the high-energy asymptotics of the Regge form, with the leading asymptotics commonly addressed as the BFKL Pomeron, so at x → 0 where P is the Pomeron intercept. As P > 0 for both the LO and the NLO BFKL Pomerons, they are called supercritical. As we are not going to use BFKL or its modifications like [21,22] in the present paper, we just mention that extensive literature on this issue can be found in Ref. [23]. Instead of using the BFKL results or trying to increase the accuracy of the method of Refs. [16,17], in the present paper we account for total resummation of the double-logarithmic contributions to F 1 . In the first place we account for the xdependent contributions and then for DL terms combining logs of x and Q 2 . These DL contributions do not involve the large factor 1/x and for this reason they have been neglected in the BFKL approach. We calculate the singlet structure function F 1 in DLA, summing DL contributions coming from virtual gluon and quark exchanges. As a result, our expressions for coefficient functions and anomalous dimensions contain total resummations of appropriate DL terms. To calculate F 1 we compose and solve Infra-Red Evolution Equations (IREEs) in the same way as we did for calculating the DIS singlet structure function g 1 (see Ref. [8]), investigating the cases of fixed and running α s . We recall that the IREE method was suggested by Lipatov in Refs. [24][25][26][27]. It is based on factorization of DL contributions of the partons with minimal transverse momenta first noticed by Gribov in Ref. [28] in the context of QED of hadrons. The technology of implementation of this method to DIS is described in detail in Ref. [8]. In contrast to DGLAP and BFKL equations, we compose the two-dimensional evolution equations: They control the evolutions in both x and Q 2 . We obtain the explicit expression for F 1 and then, applying the saddle-point method, we calculate the small-x asymptotics of F 1 automatically complemented by the asymptotic Q 2 -dependence. The asymptotics proves to be of the Regge form. The large value of the intercept compensates for the lack of the factor 1/x in the DL contributions and thereby renders the DLA asymptotics of the same order as the BFKL one. This proves that the DL contributions to F 1 at small x are, at least, no less important than the contributions coming from the BFKL Pomeron. Our paper is outlined as follows: in Sect. 2 we compose and solve IREE for the Compton amplitudes A q,g related to F 1 by the optical theorem. In this section we express A q,g through the amplitudes of the 2 → 2 scattering of partons. Those amplitudes are calculated in Sect. 3. In Sect. 4 we apply the saddle-point method to obtain an explicit expression for the small-x asymptotics of F 1 and prove that this asymptotics depends on the single variable Q 2 /x 2 instead of depending separately on Q 2 and x. In Sect. 5 we consider in detail the intercept of the Pomeron in DLA, embracing the cases of fixed and running α s . We also fix the region where the smallx asymptotics can reliably represent F 1 . Finally, Sect. 6 is for our concluding remarks.

IREE for the amplitudes of Compton scattering off partons
Following the DGLAP pattern, we consider F 1 in the framework of Collinear factorization and represent F 1 through the convolutions of the perturbative components T q and T g with non-perturbative initial quark and gluon distributions φ q,g , respectively: Throughout the paper we will consider the perturbative objects F q,g 1 only. It is convenient to consider the Compton amplitudes A q and A g related to T q,g 1 by the optical theorem: where we have introduced the factorization scale μ and used the standard notation x = Q 2 /w, with w = 2 pq and Q 2 = −q 2 . The next step is to represent A q,g in terms of the Mellin transform: where we have introduced the signature factor ξ (+) (ω) = (1 + e −ıω )/2 ≈ 1 and the logarithmic variables ρ, y (using the standard notation w = 2 pq): In what follows we will call F q , F g Mellin amplitudes and will use the same form of the Mellin transform for the other amplitudes as well. For instance, the Mellin transform for the color singlet amplitude A gg of the elastic gluon-gluon scattering in the forward kinematics is We have presumed in Eq. (10) that the virtualities of all external gluons are ∼ μ 2 . Let us notice that the only difference between the Mellin representation for the Compton amplitudes A q,g and similar amplitudes related to the singlet g 1 is in the signature factors only: the signature factor for g 1 is ξ (−) (ω) = (−1 + e −ıω )/2. Otherwise, the technology of composing and solving IREE for A gg and singlet g 1 is the same. Because of that, we present IREE for F q , F g (and for auxiliary amplitudes as well) with short comments only. The full-length derivation of all involved IREE can be found in Ref. [8]. Now all is set up to construct IREEs for F q,g . In the kinematics where the amplitudes F q , F g obey the partial differential equations: where we have used the following convenient notation: with r, r = q, g and f rr being the parton-parton amplitudes. We will calculate h rr in the next section. Actually, the equations in (12) manifest a strong resemblance with the DGLAP equations. Indeed, the first factor in brackets in the l.h.s. of (12) exists in DGLAP too. The second term vanishes when the Mellin factor (s/μ 2 ) ω is replaced by the factor x −ω , which is used in the DGLAP equations. When the parton amplitudes f rr are in the Born approximation, Eq. (12) coincides with the DGLAP equations. A general solution to Eq. (12) is where C (±) (ω) are arbitrary factors whereas and We specify the factors C (±) (ω) by the matching with the Compton amplitudes f q , f g calculated in the kinematics Q 2 ≈ μ 2 , i.e. at y = 0. The matching condition is which leads to the following expressions: Now let us express f q , f g through the parton-parton amplitudes h rr . To this end, we construct IREE for them. As f q , f g do not depend on Q 2 , the IREE for them are algebraic: where a γ q = e 2 , with e 2 being the total electric charge of the involved quacks, so that a γ q /ω is the Born value of amplitude f q (ω). There is no similar term in the equation for f g (ω). The only difference between the r.h.s. of (19) and (12) is the factor a γ q in Eq. (19). The solution to Eq. (19) is with G(ω) being the determinant of the system (19): Combining Eqs. (18) and (20), we express C ± through the parton-parton amplitudes: Combining Eqs. (22,15) and (14), we can easily express F q,g in terms of the parton-parton amplitudes h rr .

Parton-parton amplitudes
In this section we obtain explicit expressions for the parton amplitudes h rr . The IREE for h rr are quite similar to Eq. (19): where the terms b rr include the Born factors a rr and contributions of non-ladder graphs V rr : The Born factors are (see Ref. [8] for details): where A and A stand for the running QCD couplings: with η = ln(μ 2 / 2 QC D ) 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. [8] for details): with and Let us note that D = 0 when the running coupling effects are neglected. It corresponds to the total compensation of DL contributions of non-ladder Feynman graphs to scattering amplitudes with the positive signature as was first noticed in Ref. [29]. When α s is running, such a compensation is only partial. A solution to Eq. (23) is with and The algebraic equations (23) Substituting the expressions of Eq. (30) in (20), we obtain explicit expressions for the amplitudes f q , f g . Combining them with Eqs. (22), (15) and (14), we obtain explicit expressions for F q and F g . Substituting them in Eq. (8), we arrive at the explicit expressions for the Compton amplitudes A q and A g . Finally, applying the optical theorem (7) to A q and A g , we arrive at the singlet structure function F 1 .

Small-x asymptotics of the structure function F 1
The regular way to obtain the small-x asymptotics of A q and A g is to write explicit expressions for F q and F g in Eq. (8), then push x → 0 and apply the saddle-point method. However, before doing this, let us consider in detail how to calculate the asymptotics of the gluon-gluon scattering amplitude A gg , presuming virtualities of all external gluons ∼ μ 2 .

Asymptotics of F 1
The small-x asymptotics of A q and A g can be obtained with applying the saddle-point method to Eq. (8). As (+) > (−) , we neglect the terms C (−) in (14) and represent Eq. (8) in the following form: with ξ = ln(1/x) and and The stationary point at x → 0 of q is given by the rightmost root ω 0 of the following equation: When ξ → ∞, it must be equated by some negative singular contribution in the second term of Eq. (38). Using the explicit formulas for F q,g , one can conclude that such a contribution comes from the factor 1/W . So, the stationary point ω 0 is the rightmost root of the equation We consider in detail solutions to Eq. (39) at fixed and running α s in the next section. In the vicinity of ω 0 we can represent Eq. (38) as with and so in vicinity of the singularity ω 0 Expanding q in a series, we obtain In order to calculate q we notice that the most singular contributions come from differentiation of the numerator in Eq. (40), so and therefore the asymptotics of A q at x → 0 is Repeating the reasoning above for A g and applying to them the optical theorem, we conclude that the small-x asymptotics of F 1 is where the factor (ω 0 , ξ) is with δq and δq being the initial quark and gluon densities. They do not include singular factors ∼ x −a , with positive a. The Regge form of the asymptotics is brought about entirely by the perturbative contributions. Let us notice that ∼ ln −3/2 (1/x). Equation (47) exhibits that the total resummation of DL contributions leads to the Regge behavior of F 1 at small x.

Asymptotic scaling
Substituting the explicit expressions for h rr of Eq. (30) in Eq. (15) and using Eq. (39), we obtain (+) (ω 0 ) = ω 0 /2. This allows us to write the asymptotics of F 1 of Eq. (47) in the following way: Equation (49) manifests that F 1 (x, Q 2 ) at asymptotically high energies depends on the single variable Q 2 /x 2 only. We name such a confluence of the x and Q 2 dependence an asymptotic scaling. The same form of asymptotic scaling was obtained earlier for the structure function g 1 and the non-singlet component of F 1 (see Ref. [8] for details). We stress that the asymptotic scaling for F 1 can be checked by an analysis of available experimental data. Moreover, F 2 = 2x F 1 at very small x, which proves the asymptotic scaling for F 2 . Finally, let us notice that the leading singularity ω 0 in Eq. (49) does not depend on Q 2 .

Anatomy of the leading singularity ω 0
In this section we consider in detail the leading singularity ω 0 , which is the rightmost root of Eq. (39). In order to make the asymptotics of F 1 look similar to Eq. (4), we denote so that ω 0 could look similar to the BFKL leading singularity; see Eq. (4). Now let us discuss different scenarios for calculating ω 0 . In what follows we will call the DL Pomeron intercept. We recall that in the straightforward Reggeology concept = 0 and the Pomeron with > 0 is called the supercritical Pomeron.

Intercept under approximation of fixed QCD coupling
In the first place let us estimate ω 0 for the case of fixed α s . In this case DL contributions of non-ladder graphs totally cancel each other, so that D = 0 and b rr = a rr , with a rr defined in Eq. (25), where A(ω) and A (ω) should be replaced by α s . Then the solution to Eq. (39) is becomes somewhat greater: which again bears a strong resemblance to the LO BFKL intercept. However, we are positive that the approximation of fixed α s can used for rough estimating only, so we will not pursue this approximation any longer.

Intercept for the case of running coupling
Now we account for the running coupling effects in Eq. (39). Because of that, Eq. (39) can be solved only numerically. As the couplings A and A included in the factors b rr depend on μ through η = ln(μ 2 / 2 ), the solution, ω 0 , is also μdependent. Numerical calculations yield the plot of the ηdependence of ω 0 presented in Fig. 1.
The curve in Fig. 1 has the maximum ω DL A 0 = 1.066 at μ/ = 13.8. We call the optimal mass scale and call the intercept of the Pomeron in DLA. It is interesting to notice that is close to the NLO BFKL intercept. In contrast, when the quark contributions are neglected, the purely gluonic intercept DL A g is much greater:

Applicability region of the small-x asymptotics
It is obvious that the small-x asymptotic expressions, like Eq. (47), are always much simpler than non-asymptotic expressions. However, it is important to find at which values of x the asymptotics can reliably be used. To answer this question we numerically investigate R as , defined as follows: The x-dependence of R as at fixed Q 2 is shown in Fig. 2 for the case when Q 2 ≈ μ 2 .  Fig. 3 Approach of A q to its asymptotics A as q at fixed Q 2 = 100μ 2 Figure 2 demonstrates that R as = 0.9 at x ≈ 8.10 −5 , while the curve in Fig. 3, where Q 2 = 100μ 2 , grows slower and reaches the value R as = 0.9 much later, at x ≈ 3.10 −7 .
Therefore, the applicability region of the small-x asymptotics essentially depends on the Q 2 value. The plots in Figs. 2 and 3 lead us to conclude that the small-x asymptotics reliably represent F 1 in the wide range of Q 2 when x < x max , with x max ≈ 10 −6 .

Summary and outlook
In this paper we have calculated the perturbative contributions F q 1 and F g 1 to the structure function F 1 in the Double-Logarithmic Approximation, by collecting the DL contributions and at the same time accounting for the running α s effects. We obtained the explicit expressions for F q,g 1 and then, applying the saddle-point method, calculated the smallx asymptotics of F 1 , arriving at a new DL contribution to the QCD Pomeron. We demonstrated that despite the lack of the factor 1/x in the DL contributions, the impact of their total resummation makes this Pomeron supercritical, albeit the value of the intercept strongly depends on the accuracy of the calculations. The maximal value of the intercept corresponds to the roughest approximation where the quark contributions are neglected and α s is fixed. Then the value of the intercept decreases when the accuracy of the calculations increases: first, when the quark contributions are accounted for and then, notably, when the running α s effects are taken into account. Nevertheless, the Pomeron remains supercritical as ω 0 = 1.066. Such a monotonic decrease allows us to suggest that further accounting for sub-leading contributions can decrease the value of the intercept down to zero, so that eventually the intercept will satisfy the Froissart bound. We proved that the x-and Q 2 -dependencies of F 1 converge at small x in dependence on the single variable Q 2 /x 2 . We call this convergence asymptotic scaling. We stress that this prediction of asymptotic scaling can be confirmed by an analysis of the available experimental data. As asymptotically F 2 ∼ 2x F 1 , the asymptotic scaling should also take place for F 2 . Investigating the applicability region for the asymptotics, we found that F 1 can reliably be represented by its asymptotics at x ≤ x max , with x max ≈ 10 −6 .
Although we have discussed the structure function F 1 , we would like to notice that the experimental date available in the literature are mostly on the structure functions F 2 and F L , so it would be interesting to apply our approach to calculating F 2 and F L as well. Calculating F 2 in DLA can be done in a way quite similar to that we have used for F 1 . As a result, we see that F 2 in DLA can be represented through F 1 : which coincides with the well-known Born relation between F 1 and F 2 . Equation (58) entails that in DLA F L = 0. In order to estimate the deviation of F L from zero, one should account for sub-leading contributions to both F 1 and F 2 . In the first place, such contributions are the single-logarithmic (SL) ones. In this regard we recall that the SL contributions to F 2 following from emission of gluons with momenta widely separated in rapidity and not ordered in transverse momenta were accounted for in Refs. [15][16][17][18][19], which involved dealing with the BFKL characteristic function. However, there are the SL contributions unrelated to BFKL, i.e., for the case of F 1 , the SL terms unaccompanied by the factor 1/x, similar to the DL terms in Eq. (5). In contrast to the DL contributions (5), there is not a general technology in the literature for resummations of such non-BFKL SL terms. On the other hand, we were able to modify the IREE method for the spin structure function g 1 (see Ref. [8] and the references therein) to account for the SL contributions which are complementary to the ones calculated in Refs. [15][16][17][18][19], namely, the SL following from emission of the partons with momenta ordered in the k ⊥ -space and disordered in the longitudinal space. We plan to adapt this approach to calculating the SL contributions to F 1,2 .
Finally, we stress that in contrast to DGLAP we do not need singular factors ∼ x −a in fits for the initial parton distributions for F 1 . Such factors cause a steep rise of the structure functions at small x and lead to the Regge asymptotics of F 1 . However, we have shown in Sect. 4 that the resummation of the DL contributions to F 1 automatically leads to the Regge asymptotics, which makes an inclusion unnecessary of the singular terms into the fits. This result agrees with our earlier results (see Ref. [8]) for the non-singlet structure functions g 1 and F 1 and also with the results of Refs. [18,19] obtained for the small-x behavior of the structure function F 2 . The latter agreement is especially interesting because the approaches used in Refs. [18,19] and in the present paper are totally different.