Perturbative evolution: a different approach at small x

TL;DR

This paper introduces a novel perturbative evolution method for small x in deep inelastic scattering, effectively addressing the limitations of traditional approaches and providing better data description, especially for charm structure functions.

Contribution

It presents a new approach to DGLAP evolution at small x that overcomes the invalidity of perturbation expansion in this regime.

Findings

Better fit to charm structure function data

Successful description of small x behavior

Addresses limitations of conventional methods

Abstract

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.

Equations

$$xg(x,Q^2)=G(Q^2)x^{-{ϵ}_0}$$

$$\mathbf{q}(Q^2)=\left[\text{matrix}\Sigma (Q^2\text{cr}G(Q^2)\text{cr}\right]\Sigma (Q^2)=u(Q^2)+d(Q^2)+s(Q^2)+c(Q^2)+\bar{u}(Q^2)+\bar{d}(Q^2)+\bar{s}(Q^2)+\bar{c}(Q^2)$$

$$Q^2\frac{d}{d{Q}^2}\mathbf{q}(Q^2)=\frac{{\alpha }_s(Q^2)}{2\pi }\mathbf{P}(N=ϵ,{\alpha }_s(Q^2))\mathbf{q}(Q^2)$$

$$\mathbf{P}(N,{\alpha }_s)={\mathbf{P}}^0(N)+\frac{{\alpha }_s}{2\pi }{\mathbf{P}}^1(N)$$

$$P_{qg}^0(N)=2.59-2.21N+1.085N^2$$

$$P_{qg}^1(N)=24.6/N-28.2$$

$$P_{qg}=P_{qg}^0(N)+C\sqrt{C^2{N}^2+N{\alpha }_s{P}_{qg}^1(N)/\pi }-C^{2}N$$

$$=P_{qg}^0(N)+{\alpha }_s{P}_{qg}^1(N)/(2\pi )-{\alpha }_s^2[P_{qg}^1(N){]}^2/(8{\pi }^{2}N{C}^2)+\dots$$

$${\sigma }^{\text{sevenrm}red}(x,y,Q^2)=F_2(x,Q^2)-\frac{y^2}{1+(1-y{)}^2}{F}_L(x,Q^2)$$

$$F_2(x,Q^2)=A_0(Q^2)x^{-{ϵ}_0}+A_1(Q^2)x^{-{ϵ}_1}+A_2(Q^2)x^{-{ϵ}_2}$$

$$A_i(Q^2)=X_i(Q^2/(Q_i^2+Q^2){)}^{1+{ϵ}_i}(1+2Q^2/Q_i^2{)}^{0.15}i=0,1,2$$

$${ϵ}_0=0.4{ϵ}_1=0.096{ϵ}_2=-0.343$$

$$X_0=0.03527X_1=0.03786X2=0.06445$$

$$Q_0^2=7.03262Q_1^2=1.68993\times 10^{-5}{Q}_2^2=0.00192$$

$${\alpha }_S(M_Z^2)=0.1183$$

$${\Lambda }_5=230\; MeV{\Lambda }_4=328\; MeV$$

$${\bf P}^{0}(N)=\left[\matrix{-0.0352299-2.41822N+0.743989N^{2}&2.58893-2.2144N+1.08471N^{2}\cr 18.4131-49.7311N+42.476N^{2}&39.0374-115.451N+96.1001N^{2}\cr}\right]$$

$${\bf P}^{1}(N)=\left[\matrix{-1111.55-3187.40N+-2685.19N^{2}&-943.22+2952.70N-2460.29N^{2}\cr-636.76+2070.57N-1707.29N^{2}&-114.61+383.00N-359.98N^{2}\cr}\right]$$

$$G(Q^{2})=X_{g}\big{(}Q^{2}/(Q_{g}^{2}+Q^{2})\big{)}^{\eta_{1}}\,(1+Q^{2}/Q_{g}^{2})^{\eta_{2}}$$

$$X_g=0.4338{\eta }_1=0.437{\eta }_2=0.252Q_g^2=9.128$$

$$m_c(m_{c0})=m_{c0}{m}_{c0}=1.27\; GeV$$

$$\big{(}m_{c}(Q^{2})\big{)}^{2}=m_{c0}^{2}/\big{(}1+0.3438\log(Q^{2}/m_{c0}^{2})\big{)}$$

$$\mu^{2}=Q^{2}+4\big{(}m_{c}(Q^{2})\big{)}^{2}$$

$$m_b(m_{b0})=m_{b0}{m}_{b0}=4.18GeV$$

$$(m_{b}(Q^{2}))^{2}=m_{b0}^{2}/\big{(}1+0.179\log(Q^{2}/m_{b0}^{2}\big{)}$$

$$F_{L}(x,Q^{2})=K(x,Q^{2})+{4\alpha_{s}(Q^{2})\over 3\pi}\int_{x}^{1}{dy\over y}\Big{(}{x\over y}\Big{)}^{2}F_{2}(y,Q^{2})$$

$$K^{c}(x,Q^{2})=2e_{c}^{2}{\alpha_{s}(\mu^{2})\over\pi}\int_{x\mu^{2}/Q^{2}}^{1}dy\,\Big{(}{x\over y}\Big{)}^{2}\Big{[}\Big{(}1-{x\over y}\Big{)}v-{2m_{c}^{2}(Q^{2})x\over Q^{2}y}L\Big{)}\Big{]}\,g(y,\mu^{2})$$

$$v^{2}(Q^{2})=1-{{4m_{c}^{2}(Q^{2})}\over{Q^{2}(y/x-1)}}~{}~{}~{}~{}~{}~{}~{}L=\log\Big{(}{{1+v}\over{1-v}}\Big{)}$$

Full text

DAMTP-2017-10

MAN/HEP/2017/05

Perturbative evolution: a different approach at small x$$x **partially supported by STFC grant ST/L000385/1

A Donnachie

School of Physics and Astromony, University of Manchester

P V Landshoff

DAMTP, Cambridge University ††† email addresses: Sandy.Donnachie@manchester.ac.uk,  pvl@damtp.cam.ac.uk

Abstract

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.

1 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{}^{{\hbox{\rfont[1]}}}$${}^{{\hbox{\rfont[2]}}} that there are theoretical problems at small $x$. In this paper, we refine our previous approach${}^{{\hbox{\rfont[3]}}}$ to overcoming this problem, which combines the DGLAP evolution equation with Regge theory.

An important motivation is provided by the HERA data${}^{{\hbox{\rfont[4]}}}$ for charm electroproduction, which have $F_{2}^{c\bar{c}}(x,Q^{2})$ behaving at small $x$ as the same fixed power $x^{-\epsilon_{0}}$, with $\epsilon_{0}$ close to 0.4, down to very small $Q^{2}$. This is seen in figure 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_{2}^{c\bar{c}}(x,Q^{2})$. The lines in figure 1 are the result of the calculations we describe in this paper of $F_{2}^{c\bar{c}}(x,Q^{2})$ to nonleading order (NLO) in perturbative QCD. They are to be contrasted with calculations${}^{{\hbox{\rfont[5]}}}$ in the conventional approach to perturbative evolution, two of which are shown in figure 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 figure 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.

Data for high energy $pp$ and $\bar{p}p$ elastic scattering and total cross sections have long been known to be dominated by so-called soft pomeron exchange, giving a power behaviour $s^{\epsilon_{1}}$ to the amplitude. First fits${}^{{\hbox{\rfont[6]}}}$ gave a value for $\epsilon_{1}$ close to 0.08, but more recent data increase this{}^{{\hbox{\rfont[7]}}}$${}^{{\hbox{\rfont[8]}}} to 0.096. The emergence from HERA of data for $F_{2}(x,Q^{2})$ revealed${}^{{\hbox{\rfont[9]}}}$ that at small $x$ the power $x^{-\epsilon_{1}}$ alone did not fit and that a further term behaving as $x^{-\epsilon_{0}}$ needed to be added. The value of $\epsilon_{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 figure 1 that for the charm component of $F_{2}(x,Q^{2})$ the “hard-pomeron” term $x^{-\epsilon_{0}}$ is dominant, with extremely little room for the addition of the “soft-pomeron” term $x^{-\epsilon_{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 next-to-leading order. Hence $xg(x,Q^{2})$ is dominated by hard-pomeron exchange at small $x$: it behaves as $x^{-0.4}$ even at quite small $Q^{2}$:

[TABLE]

If we Mellin transform the quark and gluon densities by applying $\int_{0}^{1}dx\,x^{N-1}$, a power $u(Q^{2})x^{-\epsilon}$ in the $u$-quark density, for example, becomes a pole $u(Q^{2})/(N-\epsilon)$. Write

[TABLE]

then

[TABLE]

At nonleading order (NLO) the evolution matrix is

[TABLE]

The problem with applying this to perturbative evolution at small $x$ lies in that, while the $qq$ and $qg$ elements of ${\bf P}^{0}(N)$ are finite at $N=0$, its other two elements, and those of ${\bf P}^{1}(N)$, have poles there${}^{{\hbox{\rfont[10]}}}$. 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 ${\bf P}(N)$, it is${}^{\hbox{\rfont{[10]}}}$ $P_{qg}$ that differs most between LO and NLO. As is seen in Figure 3a, for $N$ between 0.096 and 0.4 good fits to $P^{0}_{qg}(N)$ and $P^{1}_{qg}(N)$ are

[TABLE]

[TABLE]

As a model for the exact $P_{qg}(N)$ consider

[TABLE]

[TABLE]

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$.

The problem with the perturbative expansion at small $x$ is implicit in all conventional applications of the DGLAP equation: for values of $Q^{2}$ for which the PDFs are fairly flat in $x$ at small $x$, that is the effective power $(1/x)^{\epsilon}$ is small, it is not valid to use the DGLAP equation with ${\bf P}(N,\alpha_{s}(Q^{2}))$ expanded in powers of $\alpha_{s}(Q^{2})$.

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=\epsilon_{0}$ to the coefficient of the hard pomeron contribution to $F_{2}(x,Q^{2})$. We have seen in figure 1 that the resulting gluon density leads to a highly successful calculation of $F_{2}^{c\bar{c}}(x,Q^{2})$ to NLO. Figure 4 shows data ${}^{{\hbox{\rfont[11]}}}$ ${}^{{\hbox{\rfont[12]}}}$ 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${}^{{\hbox{\rfont[13]}}}$.

2 Fit to data for F_{2}$$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${}^{{\hbox{\rfont[14]}}}$ both included results for the so-called reduced cross section

[TABLE]

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${}^{{\hbox{\rfont[15]}}}$ gave results only for $\sigma^{\hbox{{\sevenrm red}}}(x,y,Q^{2})$ and did not attempt to extract $F_{2}(x,Q^{2})$.

We fit the small-$x$ 2010 HERA data${}^{\hbox{\rfont{[14]}}}$ for $F_{2}(x,Q^{2})$ to a sum of three powers of $x$:

[TABLE]

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}\to 0$, $F_{2}(x,Q^{2})$ vanishes linearly in $Q^{2}$ at fixed $2\nu=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 GeV2.

We have tried various forms of the $A_{i}(Q^{2})$. Several give much the same value of $\chi^{2}$ in the fit to the data, but rather different values of $\epsilon_{0}$, which is why there is still a rather large uncertainty in its value. In this paper we give the results for the choices

[TABLE]

We use the values

[TABLE]

where the last two are extracted${}^{\hbox{\rfont{[8]}}}$ from $pp$ and $\bar{p}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${}^{{\hbox{\rfont[16]}}}$, 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 figure 5. The fits in this figure have

[TABLE]

[TABLE]

The $\chi^{2}$ per data point is 0.95 for all the data with $Q^{2}$ between 3.5 and 45 GeV2 and $x<0.001$.

3 Perturbative QCD evolution

We work to next-to-leading order in $\overline{MS}$ perturbative QCD, with${}^{{\hbox{\rfont[17]}}}$

[TABLE]

We follow the standard procedure${}^{\hbox{\rfont{[10]}}}$ of changing from 4 to 5 flavours at $Q^{2}=m_{b0}^{2}$ and making $\alpha_{S}(Q^{2})$ continuous there, resulting in

[TABLE]

For the range of $Q^{2}$ values that will concern us the contribution to $F_{2}(x,Q^{2})$ from $b$ quarks is${}^{{\hbox{\rfont[18]}}}$ negligibly small. Pomerons couple equally to $u$ and $d$ quarks. In the case of the soft pomeron, we know${}^{\hbox{\rfont{[6]}}}$ that the coupling to $s$ quarks is rather weaker than to $u$ and $d$ and, as may be deduced from the data in figure 2, very much so to $c$. However, the data for $F_{2}^{c\bar{c}}(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_{2}^{c\bar{c}}(x,Q^{2})\sim 0.4~{}A_{0}(Q^{2})x^{-\epsilon_{0}}$ at small $x$, and the lines in figure 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^{-\epsilon_{0}}$ at small $x$.

In order to apply DGLAP evolution, according to (1) we need∗††∗ We are grateful to the authors of reference 10 for making available to us the numerical values of the elements of ${\bf P}(N,\alpha_{s}(Q^{2}))$ the elements of the evolution matrix (2c) at $N=0.4$. We find that, for values of $N$ close to this, good fits are

[TABLE]

[TABLE]

We choose some value $\bar{Q}^{2}$ of $Q^{2}$ and calculate $\Sigma(\bar{Q}^{2})$ and its derivative $\Sigma^{\prime}(\bar{Q}^{2})$ using the form of $A_{0}(Q^{2})$ we obtained from the fit (6) to $F_{2}(Q^{2})$. The DGLAP equation (2b) with $N=0.4$ then determines $G(\bar{Q}^{2})$. This provides the input to the evolution in both directions from $Q^{2}=\bar{Q}^{2}$. Figure 6 shows the percentage difference for the choice $\bar{Q}^{2}$=9 GeV2. The plot in this figure is quite sensitive to the choice of $\bar{Q}^{2}$. It is known that perturbative evolution is valid only at sufficiently large $Q^{2}$; we conclude that this means greater than about 5 GeV2. Figure 7 shows how $G(Q^{2})$ and $\Sigma(Q^{2})$ evolve with $Q^{2}$. The value $\bar{Q}^{2}$=9 GeV 2 gives agreement with the gluon density at very large $Q^{2}$ obtained from the conventional approach to perturbative evolution${}^{{\hbox{\rfont[19]}}}$ : see figure 8. At these values of $Q^{2}$ the effective power $(1/x)^{\epsilon}$ for the behaviour of $xg(x,Q^{2})$ has become so large that the conventional approach is now valid.

4 Applications

We now apply the gluon density (1) shown in figure 7 from the NLO calculation of $G(Q^{2})$ using the perturbative evolution formula (1b). A good fit in the range $5<Q^{2}<200$ GeV2 is

[TABLE]

with

[TABLE]

Figure 1a shows the NLO calculation of $F_{2}^{c\bar{c}}(x,Q^{2})$ using an adaptation of a program**We are grateful to Johannes Bluemlein for making this program available to us originally written by Laenen, Riemersma, Smith and van Neerven${}^{{\hbox{\rfont[20]}}}$. We use the NLO running mass $m_{c}(Q^{2})$ with${}^{\hbox{\rfont{[17]}}}$

[TABLE]

A good fit, for $5<Q^{2}<10^{4}$ GeV2, is provided by

[TABLE]

In both $\alpha_{s}(Q^{2})$ and $xg(x,Q^{2})$ we replace $Q^{2}$ with

[TABLE]

Our calculation, whose results are shown in figure 1, includes only the gluon-induced charm production. The cross section for charm electroproduction is the sum of two terms, $e\,g\to\hbox{charm}$ and $e\,q\to\hbox{charm}$, each convoluted with the corresponding parton density. Figure 1 shows that the first alone, calculated up to order $\alpha_{s}^{2}$, fits the data well. For the second, the order $\alpha_{s}^{2}$ 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††† We are grateful to Eric Laenen for explaining this to us, 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_{2}^{b\bar{b}}(x,Q^{2})$, with${}^{\hbox{\rfont{[17]}}}$

[TABLE]

and the fit

[TABLE]

The resulting calculations are shown in figure 9, together with the data${}^{\hbox{\rfont{[18]}}}$.

We have said that, to calculate the $F_{L}(x,Q^{2})$, it is sufficient to use the LO formula${}^{{\hbox{\rfont[21]}}}$

[TABLE]

where the contribution of the charm quark to $K(x,Q^{2})$ is

[TABLE]

with $\mu^{2}$ defined in (13) and

[TABLE]

The light quarks contribute similarly, with zero mass. The calculations are compared with data${}^{{\hbox{\rfont[22]}}}$ for $F_{L}(x,Q^{2})$ in figure 4. The circles in figure 10 show the calculated points for the reduced cross section (4), together with the 2015 joint HERA data${}^{\hbox{\rfont{[15]}}}$. The lines in this figure are the fits to $F_{2}(x,Q^{2})$ from the 2010 HERA data${}^{\hbox{\rfont{[14]}}}$. The $\chi^{2}$ per data point is 0.96. Although, according to figure 6, the DGLAP evolution is not valid below $Q^{2}=5$ GeV2, the evolved gluon density gives good results for $F_{L}(x,Q^{2})$ down to 3.5 GeV2.

Acknowledgments

We thank Marco Guzzi for helpful discussions on DGLAP and the CT14 analysis.

References

relax1 S Catani and F Hautmann, Nuclear Physics B427 (1994) 475 – arXiv:hep-ph/9405388 J R Cudell, A Donnachie and P V Landshoff, Physics Letters B448 (1999) 281 and B526 (2002) 413 G Altarelli, R D Ball and S Forte, Nuclear Physics B599 (2001) 383 A A Almasy, N A Lo Presti and A Vogt, JHEP 01 (2016) 028 – arXiv:1511.08612 relax2 A Donnachie, H G Dosch, P V Landshoff and O Nachtmann, Pomeron Physics and QCD, Cambridge University Press (2002) relax3 A Donnachie and P V Landshoff, Physics Letters B533 (2002) 277 – arXiv:hep-ph/0111427 relax4 H1 and ZEUS collaborations: H Abramowicz et al, European Physics Journal C73 (2013) 2311 – arXiv:1211.1182 and https://www.desy.de/h1zeus/combined_results/index.php relax5 https://www.hepforge.org/archive/lhapdf/pdfsets/6.1/ and A Buckley et al, European Physics Journal C75 (2015) 132 – arXiv:1412.7420) relax6 A Donnachie and P V Landshoff, Nuclear Physics B231 (1984) 189; Physics Letters B296 (1992) 227 relax7 J R Cudell, K Kang and S K Kim, Physics Letters B395 (1997) 311 relax8 A Donnachie and P V Landshoff, Physics Letters B727 (2013) 500 – arXiv:1309.1292 relax9 A Donnachie and P V Landshoff, Physics Letters B437 (1998) 408 – arXiv:hep-ph/9806344 relax10 R K Ellis, W J Stirling and B R Webber, QCD and Collider Physics, Cambridge University Press (1996) relax11 ZEUS Collaboration: S Chekanov et al, Physics Letters B682 (2009) 8 – arXiv:0904.1092 relax12 H1 collaboration: F D Aaron et al, Physics Letters B393 (1997) 452 – arXiv:1012.4355 relax13 S Moch, J A M Vermaseren and A Vogt, Physics Letters B606 (2005) 123 relax14 H1 and ZEUS collaborations: F D Aaron et al, JHEP 01 (2010) 109 relax15 H1 and ZEUS collaborations: H Abramovicz et al, arXiv:1506.06042 relax16 R G Roberts, The Structure of the Proton, Cambridge University Press (1990) relax17 Particle Data Group, http://pdg.lbl.gov relax18 ZEUS collaboration: H Abramowicz et al, JHEP 09 (2014) 127 – arXiv:1405.6915 H1 collaboration: A Aktas et al, European Physics Journal C40 (2005) 349 – arXiv:hep-ex/0411046; F D Aaron et al, European Physics Journal C65 2010) 89 – arXiv:0907.2643 relax19 https://www.hepforge.org/archive/lhapdf/pdfsets/6.1/ relax20 E Laenen, S Riemersma, J Smith and W.L. van Neerven, Nuclear Physics B392 (1993) 162 relax21 V M Budnev et al, Physics Reports C15 (1974) 181 relax22 H1 collaboration: F D Aaron et al, C71 (2011) 1579 – arXiv:1012.4355 ZEUS collaboration: S Chekanov et al, Physics Letters B682 (2009) 8 – arXiv:0904.1092

(a)                                                               (b)

Figure 1: (a) Data${}^{\hbox{\rfont{[4]}}}$ for for charm electroproduction, with lines proportional to $x^{-0.4}$ and (b) for charm photoproduction compared with $(W^{2})^{0.4}$.

(a)                                                               (b)

Figure 2: Calculations${}^{\hbox{\rfont{[19]}}}$ of $F_{2}^{c\bar{c}}(x,Q^{2})$ at NLO and NNLO. In (a) the NNLO calculations are below the NLO, while in (b) they are above.

(a)                                                               (b)

Figure 3: (a) The evolution matrix element${}^{\hbox{\rfont{[10]}}}$ $P_{qg}$ for $\alpha_{s}=0.21$ at LO and NLO together with the fits (3a); (b) the exact $P_{qg}$ in the model (3b) (points) with the NLO and NNLO approximations

Figure 4: LO calculation of $F_{L}(x,Q^{2})$ at the indicated values of $Q^{2}$ compared with data from ZEUS${}^{\hbox{\rfont{[11]}}}$ and H1${}^{\hbox{\rfont{[12]}}}$.

Figure 5: fit to 2010 HERA data${}^{\hbox{\rfont{[14]}}}$ – only data with $x<0.001$ are included in the fit

Figure 6: Percentage difference between the evolved $\Sigma(Q^{2})$ and that obtained from the fit to $A_{0}(Q^{2})$ from the small-$x$ data for $F_{2}(x,Q^{2})$

Figure 7: Evolution of $G(Q^{2})$ and $\Sigma(Q^{2})$

Figure 8: Evolved gluon density at $Q^{2}=10^{4}$ and $10^{5}$ GeV2. The points are the results${}^{\hbox{\rfont{[19]}}}$ of conventional evolution according to HERAPDF20.

Figure 9: NLO calculations of $F_{2}^{b\bar{b}}(x,Q^{2})$ with data${}^{\hbox{\rfont{[18]}}}$ from ZEUS and H1 at various values of $Q^{2}$.

Figure 10: Data${}^{\hbox{\rfont{[15]}}}$ for the reduced cross section. The circles denote the calculations and the lines the fit to $F_{2}(x,Q^{2})$ from the 2010 data${}^{\hbox{\rfont{[14]}}}$ – 0.25 is added to the data points for each successive value of $Q^{2}$.