Perturbative evolution: a different approach at small x

We propose an approach to DGLAP evolution at small x that circumvents the usual problem that a perturbation expansion is not valid there. The data for the charm structure function are important to motivate the method, and it describes them much more successfully than the conventional approach.


Introduction
Although fixed-order perturbative QCD is widely applied to data for the proton structure function F 2 (x, Q 2 ), it is well known [1][2][3][4][5][6] that there are theoretical problems at small x. In this paper, we refine our previous approach [7] to overcoming this problem, which combines the DGLAP evolution equation with Regge theory.
An important motivation is provided by the HERA data [8] for charm electroproduction, which have F cc 2 (x, Q 2 ) behaving at small x as the same fixed power x − 0 , with 0 close to 0.4, down to very small Q 2 . This is seen in Fig. 1. We show only data having y < 0.55, so that although strictly the data are for the reduced cross section, this differs little from F cc 2 (x, Q 2 ). The lines in Fig. 1 are the result of the calculations we describe in this paper of F cc 2 (x, Q 2 ) to nonleading order (NLO) in perturbative QCD. They are to be contrasted with calculations [9] in the conventional approach to perturbative evolution, two of which are shown in Fig. 2. The conventional approach ignores the need for resummation at small x, which we do not know how properly to do and which we circumvent. If anything, the photoproduction data in Fig.  1b call for a rise even steeper than (W 2 ) 0.4 , though the line in the figure does not include threshold effects, which will pull it down a little at the lower energies.
Partially supported by STFC grant ST/L000385/1. a e-mail: Sandy.Donnachie@manchester.ac.uk b e-mail: pvl@damtp.cam.ac.uk Data for high energy pp andp p elastic scattering and total cross sections have long been known to be dominated by socalled soft pomeron exchange, giving a power behaviour s 1 to the amplitude. First fits [10,11] gave a value for 1 close to 0.08, but more recent data increase this [12,13] to 0.096. The emergence from HERA of data for F 2 (x, Q 2 ) revealed [14] that at small x the power x − 1 alone did not fit and that a further term behaving as x − 0 needed to be added. The value of 0 is about 0.4, but with quite a large uncertainty, perhaps 10%. For the complete F 2 (x, Q 2 ) both terms are needed to describe the data, but we have seen in Fig. 1 that for the charm component of F 2 (x, Q 2 ) the "hard-pomeron" term x − 0 is dominant, with extremely little room for the addition of the "soft-pomeron" term x − 1 .
In leading order, perturbative QCD charm electroproduction is calculated from the gluon density xg(x, Q 2 ), and to a very good approximation this is true also in nextto-leading order. Hence xg(x, Q 2 ) is dominated by hardpomeron exchange at small x -it behaves as x −0.4 even at quite small Q 2 : If we Mellin transform the quark and gluon densities by applying 1 0 dx x N −1 , a power u(Q 2 )x − in the u-quark density, for example, becomes a pole u( then At nonleading order (NLO) the evolution matrix is The problem with applying this to perturbative evolution at small x lies in that, while the qq and qg elements of P 0 (N ) are finite at N = 0, its other two elements, and those of P 1 (N ), have poles there [15]. This means that at N = 0.096 the NLO evolution matrix is very different from its leading order (LO) version.
Of the four elements of the evolution matrix P(N ), it is [15] P qg that differs most between LO and NLO. As is seen in Fig. 3a, for N between 0.096 and 0.4 good fits to P 0 qg (N ) and P 1 qg (N ) are As a model for the exact P qg (N ) consider where C is a constant. Figure 3b shows the exact P qg in this model, together with the NLO and NNLO approximations to it. The conclusion from this model is that, provided the NLO and the NNLO approximations are not very different, then NLO is already a good approximation. Further, this is much more likely to be true for N = 0.4 than for N = 0.096.
But we have inferred that the contribution from soft pomeron exchange to the gluon density is negligibly small. So we apply NLO perturbative QCD evolution (2b) with N = 0 to the coefficient of the hard-pomeron contribution to F 2 (x, Q 2 ). We have seen in Fig. 1 that the resulting gluon density leads to a highly successful calculation of F cc 2 (x, Q 2 ) to NLO. Figure 4 shows data [16,17] that indicate that the same is true of a calculation to leading order (LO) of the longitudinal structure function F L (x, Q 2 ). Unlike for F 2 (x, Q 2 ), the NLO correction is known to be small [18].

Fit to data for F 2
The H1 and ZEUS experiments at HERA have agreed joint data for the proton structure function F 2 (x, Q 2 ). Their 2010 publication [19] both included results for the so-called reduced cross section and an extraction of F 2 (x, Q 2 ) from it, though only for events with y less than about 0.5. Their 2015 publication [20] gave results only for σ red (x, y, Q 2 ) and did not attempt to extract We fit the small-x 2010 HERA data [19] for F 2 (x, Q 2 ) to a sum of three powers of x: Here, the first two terms are hard and soft pomeron exchange. The last term is reggeon exchange, that is, the families f 2 and a 2 ; its contribution for the small values of x at which we work is small.
The only theoretical constraint we know on the three coefficient functions A i (Q 2 ) is that, as Q 2 → 0, F 2 (x, Q 2 ) vanishes linearly in Q 2 at fixed 2ν = Q 2 /x. We respect this constraint, even though our fit is for values of Q 2 some way from zero. This is sufficient for our purpose, because perturbative evolution is valid only for "large" values of Q 2 . The theory does not tell how large but we shall find that Q 2 must be at least about 5 GeV 2 .
We have tried various forms of the A i (Q 2 ). Several give much the same value of χ 2 in the fit to the data, but rather different values of 0 , which is why there is still a rather large uncertainty in its value. In this paper we give the results for the choices where the last two are extracted [23-25] from pp andp p elastic scattering and total cross section data. We confine attention to the small-x data, which we take to mean x < 10 −3 . But we introduce into each term a power of (1 − x) given by the dimensional counting rule [21], that is, (1 − x) 7 for each of the two pomeron terms, and (1 − x) 3 for the reggeon term. Although this is not strictly correct, it is better than not doing it, at least if Q 2 is not too large. Its effect is small for small x, though not completely negligible, and it makes the fit to the x < 0.001 data valid some way beyond that range and therefore to larger Q 2 : see Fig. 5. The fits in this figure have X 0 = 0.03527, X 1 = 0.03786, X 2 = 0.06445, The χ 2 per data point is 0.95 for all the data with Q 2 between 3.5 and 45 GeV 2 and x < 0.001.

Perturbative QCD evolution
We work to next-to-leading order in M S perturbative QCD, with [22] We follow the standard procedure [15] of changing from 4 to 5 flavours at Q 2 = m 2 b0 and making α S (Q 2 ) continuous there, resulting in For the range of Q 2 values that will concern us the contribution to F 2 (x, Q 2 ) from b quarks is [23-25] negligibly small. Pomerons couple equally to u and d quarks. In the case of the soft pomeron, we know [10,11] that the coupling to s quarks is rather weaker than to u and d and, as may be deduced from the data in Fig. 2, very much so to c. However, the data for F cc 2 (x, Q 2 ) indicate that, for reasons that are not understood and to a very good approximation, the hard pomeron couples equally to all four quarks. To verify this, note that if this is the case, F cc 2 (x, Q 2 ) ∼ 0.4 A 0 (Q 2 )x − 0 at small x, and the lines in Fig. 1, while calculated by NLO perturbative QCD, correspond to this to a very good approximation. As we have said above, the gluon density must also behave as x − 0 at small x.
In order to apply DGLAP evolution, according to (1) we need 1 the elements of the evolution matrix (2c) at N = 0.4. We find that, for values of N close to this, good fits are We choose some valueQ 2 of Q 2 and calculate (Q 2 ) and its derivative (Q 2 ) using the form of A 0 (Q 2 ) we obtained from the fit (6) to F 2 (Q 2 ). The DGLAP Eq. (2b) with N = 0.4 then determines G(Q 2 ). This provides the input to the evolution in both directions from Q 2 =Q 2 . Figure 6 shows the percentage difference for the choiceQ 2 =9 GeV 2 . The plot in this figure is quite sensitive to the choice ofQ 2 . It is well known that perturbative evolution is valid only at sufficiently large Q 2 ; we conclude that this means greater than about 5 GeV 2 . Figure 7 shows how G(Q 2 ) and (Q 2 ) evolve with Q 2 . The valueQ 2 =9 GeV 2 gives agreement with the gluon density at very large Q 2 obtained from the conventional approach to perturbative evolution [26]: see Fig.  8. At these values of Q 2 the effective power (1/x) for the behaviour of xg(x, Q 2 ) has become so large that the conventional approach is now valid.
We observe that the small values of x that we consider in this paper contribute very little to the momentum sum rule. It is not trivial to achieve momentum conservation [27][28][29] but this is done by suitable choice of the parton densities at larger values of x.
(11b) Figure 1a shows the NLO calculation of F cc 2 (x, Q 2 ) using an adaptation of a program 2 originally written by Laenen, Riemersma, Smith and van Neerven [30]. We use the NLO running mass m c (Q 2 ) with [22], A good fit, for 5 < Q 2 < 10 4 GeV 2 , is provided by 2 We are grateful to Johannes Bluemlein for making this program available to us.
In both α s (Q 2 ) and xg(x, Q 2 ) we replace Q 2 with Our calculation, whose results are shown in Fig. 1, includes only the gluon-induced charm production. The cross section for charm electroproduction is the sum of two terms, e g → charm and e q → charm, each convoluted with the corresponding parton density. Figure 1 shows that the first alone, calculated up to order α 2 s , fits the data well. For the second, the order α 2 s term is the first term in the perturbative expansion. It turns out to be small at small Q 2 , but it is negative because of the inclusion of a counterterm to cancel a collinear divergence, 3 and therefore it is a very poor approximation to the exact expression. So it is appropriate to neglect it.
We perform an exactly similar calculation of F bb 2 (x, Q 2 ), with [22] m b (m b0 ) = m b0 m b0 = 4.18 GeV (14a) and the fit The resulting calculations are shown in Fig. 9, together with the data [23 -25].
We have said that, to calculate the F L (x, Q 2 ), it is sufficient to use the LO formula [31] where the contribution of the charm quark to K (x, Q 2 ) is with μ 2 defined in (13) and The light quarks contribute similarly, with zero mass. The calculations are compared with data [16,32] for F L (x, Q 2 ) in Fig. 4. The circles in Fig. 10 show the calculated points for the reduced cross section (4), together with the 2015 joint HERA data [20]. The lines in this figure are the fits to F 2 (x, Q 2 ) from the 2010 HERA data [19]. The χ 2 per data point is 0.96. Although, according to Fig. 6, the DGLAP evolution is not valid below Q 2 = 5 GeV 2 , the evolved gluon density gives good results for F L (x, Q 2 ) down to 3.5 GeV 2 .