Home | About | Journals | Submit | Contact Us | Français |

**|**Springer Open Choice**|**PMC5321399

Formats

Article sections

- Abstract
- Introduction
- The effective exponents
- Comparison with nonperturbative predictions
- Conclusions and outlook
- References

Authors

Related links

The European Physical Journal. C, Particles and Fields

Eur Phys J C Part Fields. 2016; 76(7): 383.

Published online 2016 July 8. doi: 10.1140/epjc/s10052-016-4240-4

PMCID: PMC5321399

Emanuele R. Nocera, Email: ku.ca.xo.scisyhp@arecon.eleuname.

Received 2016 May 30; Accepted 2016 June 28.

Copyright © The Author(s) 2016

It has been argued from the earliest days of quantum chromodynamics that at asymptotically small values of *x* the parton distribution functions (PDFs) of the proton behave as *x*^{α}, where the values of *α* can be deduced from Regge theory, while at asymptotically large values of *x* the PDFs behave as (1-*x*)^{β}, where the values of *β* can be deduced from the Brodsky–Farrar quark counting rules. We critically examine these claims by extracting the exponents *α* and *β* from various global fits of parton distributions, analysing their scale dependence, and comparing their values to the naive expectations. We find that for valence distributions both Regge theory and counting rules are confirmed, at least within uncertainties, while for sea quarks and gluons the results are less conclusive. We also compare results from various PDF fits for the structure function ratio ${F}_{2}^{n}/{F}_{2}^{p}$ at large *x*, and caution against unrealistic uncertainty estimates due to overconstrained parametrisations.

An accurate determination of parton distribution functions (PDFs) is an essential building block for the precision physics program at the large hadron collider (LHC) [1–5]. Given current limitations in the understanding of nonperturbative quantum chromodynamics (QCD), such a determination is not achievable from first principles. Instead, PDFs are determined in a global fit to hard-scattering experimental data [6–11], using perturbative QCD to combine information from different processes and scales. In such an analysis, the best-fit values of the input PDF parametrisation are obtained by comparing the PDF-dependent prediction of a suitable set of physical observables with their measured values, and then by minimising a figure of merit which quantifies the agreement between the two.

The parametrisation of the PDFs, $x{f}_{i}(x,{Q}_{0}^{2})$, is set at an initial scale ${Q}_{0}^{2}$, and is then evolved to any other scale *Q*^{2} via DGLAP equations [12–14]. The PDF parametrisation should be as general as possible, and in particular sufficiently smooth and flexible enough to accommodate all of the experimental data included in the fit without artificial bias. The kinematic constraint that $x{f}_{i}(x,{Q}_{0}^{2})$ vanishes in the elastic limit *x* → 1 should also be implicit in the parametrisation. Usually, the following *ansatz* is adopted

$$\begin{array}{c}\hfill x{f}_{i}(x,{Q}_{0}^{2})={A}_{{f}_{i}}\phantom{\rule{0.166667em}{0ex}}{x}^{{a}_{{f}_{i}}}\phantom{\rule{0.166667em}{0ex}}{(1-x)}^{{b}_{{f}_{i}}}\phantom{\rule{0.166667em}{0ex}}\mathcal{F}(x,\{{c}_{{f}_{i}}\}),\end{array}$$

1

where *x* is the parton momentum fraction and *i* denotes a given quark flavour (or flavour combination) or the gluon, and ℱ(*x*, {*c*_{fi}}) is a smooth function which remains finite both when *x* → 0 and *x* → 1. The normalisation fractions *A*_{fi}, the exponents *a*_{fi} and *b*_{fi}, and the set of parameters {*c*_{fi}} are then determined from the data. Some of the *A*_{fi} can be fixed in terms of the other fit parameters by means of the momentum and valence sum rules.

The original motivation for Eq. (1) was the theoretical expectation, based on nonperturbative QCD considerations, of a power-law behaviour of the PDFs at sufficiently small and large values of *x*. Specifically, Regge theory [15] predicts

$$\begin{array}{c}\hfill x{f}_{i}(x,{Q}^{2})\stackrel{x\to 0}{\to}{x}^{{a}_{{f}_{i}}}\u037e\end{array}$$

2

while the Brodsky–Farrar quark counting rules [16] predict

$$\begin{array}{c}\hfill x{f}_{i}(x,{Q}^{2})\stackrel{x\to 1}{\to}{(1-x)}^{{b}_{{f}_{i}}}\u037e\end{array}$$

3

see also Refs. [17, 18], and references therein. Both Regge theory and the counting rules provide numerical predictions for the values of the exponents *a*_{fi} and *b*_{fi}. In Eq. (1), the small- and large-*x* power-law behaviours are matched at intermediate *x* values through the function ℱ(*x*, {*c*_{fi}}). A number of different parametrisations have been used for this function so far, ranging from simple polynomials to more sophisticated Chebyshev [7, 19] and Bernstein [8] polynomials and multi-layer neural networks [20, 21].

It should be emphasised that Eqs. (2)–(3) *cannot* be derived using perturbative QCD, but rather require other more general considerations. For instance, counting rules can be derived from Bloom–Gilman duality [22] or using AdS/QCD methods in nonperturbative QCD [23].^{1} The use of Eqs. (2)–(3) in the input PDF parametrisation, Eq. (1), could therefore lead to theoretical bias. For instance, as we will discuss below, perturbative QCD calculations predict a logarithmic, rather than a power-like, growth of the PDFs at small *x*. Even if Eqs. (2)–(3) were a solid prediction from QCD (which they are not), they would not be particularly useful in the context of a global PDF analysis. First, it is unclear how small or large *x* should be in order for the power laws (2)–(3) to provide a reliably enough approximation of the underlying PDFs. Second, it is unclear at which values of *Q*^{2} Regge theory and Brodsky–Farrar quark counting rules should apply exactly. This is a serious limitation, given the non-negligible PDF scale dependence around the input parametrisation scale ${Q}^{2}\simeq {Q}_{0}^{2}$. In principle, the optimal values of *Q*^{2} should be chosen at the interface between perturbative and nonperturbative hadron dynamics, ${Q}^{2}\simeq {Q}_{0}^{2}={Q}_{\mathrm{in}}^{2}$. It has been shown [25] that ${Q}_{\mathrm{in}}^{2}\simeq 0.75$ GeV^{2} by matching the high- and low-*Q*^{2} behaviour of the strong coupling *α*_{s}(*Q*^{2}) as predicted respectively by its renormalisation group equation in the $\overline{\mathrm{MS}}$ scheme and its analytic form in the light-front holographic approach.

The aim of this study is to present a methodology to quantify the effective asymptotic behaviour of PDFs at small and large values of *x*, and then apply it to compare recent global fits with various perturbative and nonperturbative QCD predictions. The paper is organised as follows. In Sect. 2 we introduce a definition of the effective PDF exponents, and we use them to quantify for which ranges of *x* and *Q*^{2}, if any, PDFs exhibit a power-law behaviour of the form Eqs. (2)–(3). Once the asymptotic range has been determined, in Sect. 3 we investigate to which extent these exponents, as obtained from global PDF fits, are in agreement with the theoretical predictions of their values. In addition to Brodsky–Farrar quark counting rules, we will also compare the global fit predictions with other nonperturbative models of nucleon structure at large *x*. In principle, this comparison will allow us to discriminate among models, in the same way as was done for spin-dependent PDFs in Ref. [26].

In this paper we will compute the effective exponents *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}), which, when ${Q}^{2}={Q}_{0}^{2}$, are asymptotically equal to the exponents *a*_{fi} and *b*_{fi} of the input PDF parametrisation Eq. (1). Specifically, we define

$$\begin{array}{cc}& {\mathit{\alpha}}_{{f}_{i}}(x,{Q}^{2})\equiv \frac{\mathit{\partial}ln[x{f}_{i}(x,{Q}^{2})]}{\mathit{\partial}lnx},\hfill \\ \hfill & {\mathit{\beta}}_{{f}_{i}}(x,{Q}^{2})\equiv \frac{\mathit{\partial}ln[x{f}_{i}(x,{Q}^{2})]}{\mathit{\partial}ln(1-x)}\phantom{\rule{0.166667em}{0ex}},\hfill \end{array}$$

4

so that, at the input parametrisation scale ${Q}_{0}^{2}$,

$$\begin{array}{cc}\hfill {\mathit{\alpha}}_{{f}_{i}}(x,{Q}_{0}^{2})=& {a}_{{f}_{i}}+x\left[\frac{\mathrm{d}ln[\mathcal{F}(x,\{{c}_{{f}_{i}}\})]}{\mathrm{d}x}-\frac{{b}_{{f}_{i}}}{1-x}\right]\stackrel{x\to 0}{\to}{a}_{{f}_{i}}\hfill \\ \hfill & +\phantom{\rule{0.166667em}{0ex}}O(x),\hfill \end{array}$$

5

and

$$\begin{array}{cc}\hfill {\mathit{\beta}}_{{f}_{i}}(x,{Q}_{0}^{2})=& {b}_{{f}_{i}}-(1-x)\left[\frac{\mathrm{d}ln[\mathcal{F}(x,\{{c}_{{f}_{i}}\})]}{\mathrm{d}x}+\frac{{a}_{{f}_{i}}}{x}\right]\stackrel{x\to 1}{\to}{b}_{{f}_{i}}\hfill \\ \hfill & +\phantom{\rule{0.166667em}{0ex}}O(1-x),\hfill \end{array}$$

6

since in both Eqs. (5) and (6) the term in square brackets is by construction of order one in the corresponding limit. Because subasymptotic terms of *O*(*x*) tend to zero very quickly at small *x*, and likewise subasymptotic terms of *O*(1 - *x*) tend to zero very quickly at large *x*, we expect that the definitions Eq. (4) *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}) can be used to accurately determine the asymptotic behaviour of any given PDF *x**f*_{i}(*x*, *Q*^{2}).

In order to test this assertion, we have used Eq. (4) to compute the effective asymptotic exponents *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}) for the MSTW08 NLO PDF set [27] (see Appendix 1 for details). Results at *Q*^{2} = 1 GeV^{2}, which coincides with the input parametrisation scale ${Q}_{0}^{2}$, are shown in Fig. 1 for the up valence quark, ${f}_{i}={u}_{V}=u-\overline{u}$, the down valence quark, ${f}_{i}={d}_{V}=d-\overline{d}$, and the gluon, *f*_{i} = *g*, PDFs. They are compared to the corresponding fitted exponents *a*_{fi} and *b*_{fi}, to which they are expected to approach asymptotically. In Table 1 we show the numerical values computed respectively at *x* = 10^{-5} and *x* = 0.9, and we again compare them with the corresponding fitted exponents *a*_{i} and *b*_{i}.

From Fig. 1 and Table 1 it is clear that both *α*_{fi}(*x*, *Q*^{2}) at *x* = 10^{-5} and *β*_{fi}(*x*, *Q*^{2}) at *x* = 0.9 have converged to the fitted values of *a*_{fi} and *b*_{fi} within PDF uncertainties. In addition, by examining the *x* dependence of *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}), it is possible to identify the asymptotic regions in which they become roughly independent of *x*. Furthermore, since the definitions Eq. (4) may be applied at any value of *Q*^{2}, we may use them to study the *Q*^{2} dependence of the effective exponents.

The definition of the PDF effective exponents, Eq. (4), is robust and we can therefore use it to compare the results of global fits among themselves and with different predictions from perturbative and nonperturbative QCD. We will focus on the up and down valence PDFs, ${u}_{V}=u-\overline{u}$ and ${d}_{V}=d-\overline{d}$, the total quark sea, $S=2(\overline{u}+\overline{d})+s+\overline{s}$, and the gluon, *g*, from the NNPDF3.0 [6], MMHT14 [7], and CT14 [8] NNLO fits. We will also present some results from the ABM12 NNLO [9] and CJ15 NLO [11] sets. A detailed discussion of the similarities and differences between these PDF sets can be found in Refs. [2–4]; here we restrict ourselves to the information relevant for their small and large-*x* behaviour.

- NNPDF3.0 PDFs are parametrised in the basis that diagonalises the DGLAP evolution equations [28]. The function ℱ(
*x*, {*c*_{fi}}) is a multi-layer feed-forward neural network (also known as*perceptron*). The power-law term*x*^{afi}(1-*x*)^{bfi}in Eq. (1) is treated as a preprocessing factor that optimises the minimisation process: the exponents*a*_{fi}and*b*_{fi}are chosen for each Monte Carlo replica at random in a given range determined iteratively. - MMHT14 The PDFs parametrised are the valence distributions
*u*_{V}and*d*_{V}, the total sea*S*, the sea asymmetry ${\mathrm{\Delta}}_{S}=\overline{d}-\overline{u}$, the total and valence strange distributions ${s}^{+}=s+\overline{s}$ and ${s}^{-}=s-\overline{s}$ and the gluon*g*. The function ℱ(*x*, {*c*_{fi}}) is taken to be a linear combination of Chebyshev polynomials. The exponents*a*_{fi}and*b*_{fi}are fitted, except for*a*_{s+}=*a*_{S}. - CT14 The PDFs parametrised are the valence distributions
*u*_{V}and*d*_{V}, the sea quark distributions $\overline{u}$ and $\overline{d}$, the total strangeness*s*^{+}and the gluon*g*. It is assumed that $s=\overline{s}$. The function ℱ(*x*, {*c*_{fi}}) is a linear combination of Bernstein polynomials. The exponents*a*_{fi}and*b*_{fi}are parameters of the fit, but not all of them are free: specifically, it is assumed that*b*_{uV}=*b*_{dV}, so that as*x*→ 1 ${u}_{V}(x,{Q}_{0}^{2})/{d}_{V}(x,{Q}_{0}^{2})\to k$, with*k*a constant, and that as*x*→ 0 $\overline{u}(x,{Q}_{0}^{2})/\overline{d}(x,{Q}_{0}^{2})\to 1$, which requires ${a}_{\overline{u}}={a}_{\overline{d}}$. - ABM12 The PDFs parametrised are the valence distributions
*u*_{V}and*d*_{V}, the sea distributions $\overline{u}$ and*s*, the sea asymmetry Δ_{S}and the gluon*g*. It is assumed that $s=\overline{s}$. The function ℱ(*x*, {*c*_{fi}}) has the form*x*^{Pfi(x)}, where*P*_{fi}(*x*) is a function of*x*; for*s*, ℱ(*x*, {*c*_{fi}} = 1. The exponents*a*_{fi}and*b*_{fi}are parameters of the fit, except for the condition*a*_{ΔS}= 0.7. - CJ15 The PDFs parametrised are the valence distributions
*u*_{V}and*d*_{V}, the light antiquark sea, $\overline{u}+\overline{d}$, the light antiquark ratio $\overline{d}/\overline{u}$, the total strangeness*s*^{+}and the gluon*g*. It is assumed that $s=\overline{s}$. The function ℱ(*x*, {*c*_{fi}}) is provided by the polynomial $(1+{c}_{{f}_{i}}^{(1)}\sqrt{x}+{c}_{{f}_{i}}^{(2)}x)$ for all the distributions except the light antiquark ratio and the total strangeness. Specifically, $\overline{d}/\overline{u}$ is parametrised with a simple polynomial which ensures that as*x*→ 1, $\overline{d}/\overline{u}\to 1$, while it is assumed that ${s}^{+}=\mathit{\kappa}(\overline{u}+\overline{d})$; ${c}_{{f}_{i}}^{(1)}$, ${c}_{{f}_{i}}^{(2)}$ and*κ*are parameters of the fit. A small admixture of*u*_{V}is added to*d*_{V}so that as*x*→ 1*d*_{V}/*u*_{V}→*k*, with*k*a constant.

Although the momentum distributions of strange and antistrange quarks are assumed to be identical in some of these PDF sets, it should be noted that a strange/antistrange asymmetry in the nucleon is predicted based on nonperturbative QCD models; see e.g. Ref. [29] and references therein. Strange and antistrange distributions may also be very different from each other in the polarised case, as it was shown in Ref. [29] based on a light-cone model of energetically favoured meson-baryon fluctuations applied to the *K*^{+}Λ. However, a study of a structured asymmetry in the momentum distributions of strange and antistrange quarks in a global QCD analysis is beyond the scope of this work, and has been addressed elsewhere [6, 7].

In Figs. 2, ,3,3, ,44 we compare both the PDFs and the corresponding effective exponents *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}) for the NNPDF3.0, MMHT14 and CT14 sets at *Q*^{2} = 2 GeV^{2}. For NNPDF3.0, PDF uncertainties are computed as 68% confidence level (CL) intervals, while for MMHT14 and CT14 sets we show the symmetric one-sigma Hessian uncertainties. In most cases it is possible to identify an asymptotic region where the effective exponents become approximately independent of *x*. The onset of this asymptotic regime depends on both the PDF flavour and on the PDF set. At small *x*, the asymptotic regime is reached at *x* ≲ 10^{-3} for *u*_{V}, *d*_{V} and *S* irrespective of the PDF set considered. For the gluon, convergence is achieved at smaller values of *x*, *x* ≲ 10^{-5}, at least for MMHT14 for which *α*_{g}(*x*, *Q*^{2}) has an oscillation in the region 10^{-4} ≲ *x* ≲ 10^{-3}. Note that at *x* ≲ 10^{-4} PDFs are extrapolated into a region with very limited experimental information. This very small-*x* region can be probed at the LHC with forward charm [30, 31] and quarkonium production [32]. At large *x*, the asymptotic region is reached at *x* ≳ 0.7 in most cases. The exception is *β*_{dV}(*x*, *Q*^{2}) from MMHT14, which exhibits an oscillation in the region 0.6 ≲ *x* ≲ 0.8.

To the best of our knowledge, this is the first time that the onset of an asymptotic regime in the effective PDF exponents *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}) has been explicitly demonstrated. Remarkably, this onset takes place at *x* values close to the boundary between the data and extrapolation regions. Our results indicate that the three global PDF sets are broadly consistent among one other within uncertainties not only at the level of PDFs, but also at the level of their small- and large-*x* asymptotic behaviour. The main exceptions are *u*_{V} and *d*_{V} at small *x*, where the effective exponent of NNPDF3.0 is incompatible with those of CT14 and MMHT14. However, this is an extrapolation region where the Hessian approximation has some limitations and non-Gaussian effects are large: indeed, if we compute with NNPDF3.0 the one-sigma PDF interval as opposed to the 68 % CL, the three sets become consistent.

Before we compare our results to the expectations of Regge theory and the Brodsky–Farrar quark counting rules, we first examine the *Q*^{2} dependence of the effective exponents. To this end, in Figs. 5, ,66 we show the effective exponents *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}) as functions of *Q*^{2} at fixed values of *x* in the asymptotic region: *x* = 10^{-4} and *x* = 0.9, respectively. We show results for the valence distributions *u*_{V} and *d*_{V}, the total quark singlet $\mathrm{\Sigma}={\sum}_{i=1}^{{n}_{f}}({q}_{i}+{\overline{q}}_{i})$ and the gluon. From these plots we can see that as *Q*^{2} increases the effective exponents become less sensitive to *Q*^{2} and tend to converge to a finite value asymptotically. This feature is broadly independent of *x* when *x* is sufficiently small or large, roughly *x* ≲ 10^{-3} and *x* ≳ 0.9. The only exception is again *β*_{dV}(*x*, *Q*^{2}) for MMHT14.

At small *x*, the *Q*^{2} dependence of the effective exponents illustrates the transition from a low-*Q*^{2} region, where PDFs are determined from nonperturbative dynamics, to a high-*Q*^{2} region, where PDFs are dominated by perturbative QCD evolution. Indeed, as *x* → 0 and *Q*^{2} → *∞*, PDFs can be solely determined by DGLAP equations [14, 33], provided that their behaviour is sufficiently soft at the input scale. In this limit, it is well known that PDFs exhibit a double asymptotic scaling (DAS) [34–37]. Specifically, as *x* → 0 and *Q*^{2} → *∞* the singlet sector grows as

$$\begin{array}{cc}\hfill x\mathrm{\Sigma}(x,{Q}^{2})\to & {\mathcal{N}}_{\mathrm{\Sigma}}\frac{\mathit{\gamma}}{\mathit{\rho}}\frac{1}{\sqrt{4\mathit{\pi}\mathit{\gamma}\mathit{\sigma}}}{e}^{2\mathit{\gamma}\mathit{\sigma}-\mathit{\delta}\mathit{\sigma}/\mathit{\rho}},\hfill \\ \hfill xg(x,{Q}^{2})\to & {\mathcal{N}}_{g}\frac{1}{\sqrt{4\mathit{\pi}\mathit{\gamma}\mathit{\sigma}}}{e}^{2\mathit{\gamma}\mathit{\sigma}-\mathit{\delta}\mathit{\sigma}/\mathit{\rho}},\hfill \end{array}$$

7

where we have defined

$$\begin{array}{cc}\hfill \mathit{\gamma}\equiv & {\left(\frac{12}{{\mathit{\beta}}_{0}}\right)}^{1/2}\phantom{\rule{0.166667em}{0ex}},\phantom{\rule{2em}{0ex}}\mathit{\delta}\equiv \left(11+\frac{2{n}_{f}}{27}\right)/{\mathit{\beta}}_{0},\hfill \\ \hfill {\mathit{\beta}}_{0}=& 11-\frac{2}{3}{n}_{f}\phantom{\rule{0.166667em}{0ex}},\hfill \end{array}$$

8

and the double scaling variables

$$\begin{array}{cc}& \mathit{\sigma}\equiv {\left[ln\frac{{x}_{0}}{x}ln\frac{ln\left({Q}^{2}/{\mathrm{\Lambda}}^{2}\right)}{ln\left({Q}_{0}^{2}/{\mathrm{\Lambda}}^{2}\right)}\right]}^{1/2},\hfill \\ \hfill & \mathit{\rho}\equiv {\left[\frac{ln\left({x}_{0}/x\right)}{ln\left(ln\left({Q}^{2}/{\mathrm{\Lambda}}^{2}\right)/ln\left({Q}_{0}^{2}/{\mathrm{\Lambda}}^{2}\right)\right)}\right]}^{1/2}\phantom{\rule{0.166667em}{0ex}}\text{.}\phantom{\rule{0.333333em}{0ex}}\hfill \end{array}$$

9

The parameters *x*_{0} and ${Q}_{0}^{2}$ define the formal boundaries of the asymptotic region, 𝒩_{Σ} and 𝒩_{g} are normalisation constants, and *n*_{f} is the number of active flavours. Using the asymptotic form Eq. (7) in the definition of the effective exponents Eq. (4) then gives us a perturbative prediction for the small-*x* exponents *α*_{Σ} and *α*_{g}: at large *σ* but fixed *ρ* one has

$$\begin{array}{c}\hfill {\mathit{\alpha}}_{\mathrm{\Sigma}}(x,{Q}^{2})\to -\frac{\mathit{\gamma}}{\mathit{\rho}}+\frac{3}{4\mathit{\sigma}\mathit{\rho}},\phantom{\rule{1em}{0ex}}{\mathit{\alpha}}_{g}(x,{Q}^{2})\to -\frac{\mathit{\gamma}}{\mathit{\rho}}+\frac{1}{4\mathit{\sigma}\mathit{\rho}}.\end{array}$$

10

Note that both *α*_{Σ}(*x*, *Q*^{2}) and *α*_{g}(*x*, *Q*^{2}) converge asymptotically to the same value - *γ*/*ρ*, as expected since the QCD evolution of the gluon distribution seeds the evolution of the quark singlet distribution. The DAS results Eq. (10), which are a generic prediction of perturbative QCD, are displayed in Fig. 6, where we have used *x*_{0} = 0.1, ${Q}_{0}^{2}=1$ GeV^{2}, *n*_{f} = 5 and Λ^{(nf=5)} = 0.220 GeV. The agreement between the expectation from DAS and results from the global fits is excellent at *Q*^{2} ≳ 10 GeV^{2} for both the quark singlet and the gluon.

At large *x*, the *Q*^{2} dependence of the effective exponents can also be determined from general perturbative QCD considerations, following directly from the universality of the cusp quark anomalous dimension in the $\overline{\mathrm{MS}}$ scheme [38, 39]. Specifically, it can be shown, either by analysing Wilson lines [38], or by using standard results for the exponentiation of soft logarithms in the quark-initiated bare cross sections [39], that the quark anomalous dimension at large *N* takes the universal form

11

where *c*(*α*_{s}(*q*^{2})) and *d*(*α*_{s}(*q*^{2})) can be computed perturbatively: for example at NLO

$$\begin{array}{c}\hfill c({\mathit{\alpha}}_{s}({q}^{2}))=\frac{{\mathit{\alpha}}_{s}({q}^{2})}{2\mathit{\pi}}{c}_{1}+{\left(\frac{{\mathit{\alpha}}_{s}({q}^{2})}{2\mathit{\pi}}\right)}^{2}{c}_{2}+O({\mathit{\alpha}}_{s}^{3}),\end{array}$$

12

with coefficients [40]

$$\begin{array}{c}\hfill {c}_{1}=\frac{8}{3},\phantom{\rule{1em}{0ex}}{c}_{2}=4\left(\frac{67}{9}-2{\mathit{\zeta}}_{2}\right)-\frac{40}{27}{n}_{f}.\end{array}$$

13

It follows [39] that, if $x{f}_{q}(x,{Q}_{0}^{2})\sim {(1-x)}^{b({Q}_{0}^{2})}$ as *x* → 1 at a scale ${Q}_{0}^{2}$, with *q* either the quark singlet, Σ, or one of the quark valence distributions, *u*_{V} or *d*_{V}, then this asymptotic behaviour persists at higher scales *Q*^{2} with

$$\begin{array}{c}\hfill b({Q}^{2})=b({Q}_{0}^{2})+{\int}_{{Q}_{0}^{2}}^{{Q}^{2}}\frac{d{q}^{2}}{{q}^{2}}c({\mathit{\alpha}}_{s}({q}^{2})).\end{array}$$

14

Given our definition Eq. (4) and the asymptotic behaviour Eq. (6) at large *x*, as *x* → 1 one has

$$\begin{array}{c}\hfill {\mathit{\beta}}_{{f}_{i}}(x,{Q}^{2})={\mathit{\beta}}_{{f}_{i}}(x,{Q}_{0}^{2})+{\int}_{{Q}_{0}^{2}}^{{Q}^{2}}\frac{d{q}^{2}}{{q}^{2}}c({\mathit{\alpha}}_{s}({q}^{2})).\end{array}$$

15

The behaviour predicted by Eq. (15) is displayed for *u*_{V} and *d*_{V} in Fig. 5, and for Σ in Fig. 6. Note that Eq. (15) only determines the shape of the curve, not its overall normalisation; for definiteness we fix the value of $\mathit{\beta}(x,{Q}_{0}^{2})$ in Eq. (15) to match the central values obtained from CT14 at *Q*^{2} = 10^{6} GeV^{2}. The agreement between Eq. (15) and the *Q*^{2} dependence of the large-*x* effective exponents derived from the PDF fit is excellent. A slight deterioration only appears at small values of *Q*^{2} due to missing higher order corrections. Similar conclusions can be derived for other PDF sets when the value of ${\mathit{\beta}}_{{f}_{i}}(x,{Q}_{0}^{2})$ in Eq. (15) is assigned consistently.

We now discuss how our findings compare with the expectations from Regge theory and the Brodsky–Farrar quark counting rules. In Tables 2, ,33 we show the values of the effective exponents for the NNPDF3.0, CT14, MMHT14, ABM12 and CJ15 PDF sets, computed at *x*_{a} = 10^{-4} and *x*_{b} = 0.9 (*x*_{b} = 0.5 for *S*) at *Q*^{2} = 2 GeV^{2} and *Q*^{2} = 10 GeV^{2}. We also include the values predicted by Regge theory and the Brodsky–Farrar quark counting rules.

The values of the small-*x* effective exponent *α*_{fi}(*x*_{a}, *Q*^{2}) computed at *Q*^{2} = 2 GeV^{2} and *Q*^{2} = 10 GeV^{2} at *x*_{a} = 10^{-4}, compared to the values of *a*_{fi} predicted by Regge theory (and resummation of double logarithms). For the quark sea *S* and the **...**

At small *x*, Regge theory predicts *x**f*_{i} ∼ *x*^{afi} with *a*_{fi} a *Q*^{2}-independent exponent, related to the intercept of the corresponding Regge trajectory. For valence quark distributions, a value of *a*_{uV} = *a*_{dV} ≃ + 0.5 is derived from the non-singlet Regge trajectory intercept 1 - *α*_{R}(0). Perturbative calculations which resum the double logarithms of *x* give a similar value *a*_{uV} = *a*_{dV} ≃ + 0.63 [41, 42]. For the gluon distribution, a value of *a*_{g} close to the singlet Pomeron trajectory 1 - *α*_{P}(0) is expected; the conventional Regge exchange is that of the soft Pomeron [43] (for a formulation of the parton picture without recourse to perturbation theory see also Ref. [44]), leading to *a*_{g} ≃ - 0.08. Attempts to compute the Pomeron intercept perturbatively by solution of the fixed coupling LLx BFKL equation [45–48] give *a*_{g} ≃ - 0.5. However, this result is destabilised by NLLx corrections [49]. When running coupling effects are taken into account, the perturbative expansion is stabilised [50–54], and the NLLx perturbative prediction becomes *a*_{g} ≃ - 0.2. For the total sea distribution, the value of *a*_{S} should be similar for large enough *Q*^{2} to *a*_{g}, due to the dominance of the process $g\to q\overline{q}$ in the evolution of sea quarks.

In comparing these expectations with the results from PDF fits, we need to choose a scale. Regge predictions are expected to hold only at low scales. For *α*_{uV}(*x*, *Q*^{2}) and *α*_{dV}(*x*, *Q*^{2}) this is not too much of a problem, since the scale dependence of non-singlet distributions is quite weak (see Fig. 5). The values extracted from NNPDF3.0 are accordingly in good agreement with Regge expectations; those from the other global PDF fits are generally a little high (see Table 2). On the other hand, for *α*_{S}(*x*, *Q*^{2}) ≃ *α*_{Σ}(*x*, *Q*^{2}) and *α*_{g}(*x*, *Q*^{2}), the scale dependence is rather strong (see Fig. 6), due to the double scaling behaviour. Making the comparison at low scales, we see reasonable agreement for the sea quarks with the Pomeron prediction, and also with the NLLx perturbative prediction. Uncertainties for the gluon intercept are inevitably large, so here the agreement is only qualitative. Note that for ABM12 and CJ15 the uncertainties are often substantially underestimated due to parametrisation constraints in the extrapolation region.

At large *x*, the Brodsky–Farrar quark counting rules predict that *x**f*_{i} ∼ (1-*x*)^{2ns-1}, where *n*_{s} is the minimum number of *spectator* partons. These are defined to be the partons that are not struck in the hard-scattering process, since it is assumed that, in the limit *x* → 1, there can be no momentum left for any of the partons other than the struck parton. In a proton made of three quarks, one has for a valence quark, *n*_{s} = 2 and thus *b*_{uV} = *b*_{dV} = 3; for a gluon, *n*_{s} = 3 and *b*_{g} = 5; for a sea quark, *n*_{s} = 4 and *b*_{S} = 7; see Fig. 7. Note that the values of the exponents predicted by Brodsky–Farrar quark counting rules are different if the polarisation of the quark with respect to the polarisation of the parent hadron is retained [55]. This also affects the difference between up and down distributions. A detailed comparison between PDFs and quark counting rules in the polarised case was presented in Ref. [26]. Again it is unclear from the quark model argument at which scale these predictions are supposed to apply, but again we are fortunate that the scale dependence of large-*x* PDFs is reasonably moderate (see Figs. 5, ,6),6), and it is reasonable to make the comparison at a low scale [25].

The number of *spectator* partons in a proton consisting of three quarks, whether a valence quark (**a**), a gluon (**b**) or a sea quark (**c**) is struck by a virtual photon *γ*^{∗} in deep-inelastic scattering

The predictions *b*_{uV}(*x*, *Q*^{2}) and *b*_{dV}(*x*, *Q*^{2}) for the valence distributions are then in broad agreement with the effective exponents determined from most of the global PDF fits, though some deviations from Brodsky–Farrar quark counting rule expectations are observed for the MMHT14 down valence quarks: this seems to be a result of the oscillation noted already in Fig. 2. For the quark sea and the gluon, the success is again rather mixed, and only CT14 seems to provide results which agree with the prediction; for NNPDF3.0 the uncertainties on the quark sea are too large for the extraction to be meaningful, while the result for the gluon is a little low; for MMHT14 the result for the gluon is far too low, with a substantially underestimated uncertainty.

In addition to the Brodsky–Farrar quark counting rules, the behaviour of PDFs at large *x* has been predicted by several nonperturbative models of nucleon structure (see *e.g.* [56, 57] and references therein). In many cases, these provide expectations for the ratio of *u* to *d* valence distributions in the proton, *d*_{V}/*u*_{V}, and of neutron to proton structure functions, ${F}_{2}^{n}/{F}_{2}^{p}$. These ratios are particularly interesting because while all PDFs vanish at *x* = 1, their ratio does not necessarily do so, and thus it is a useful discriminator among models of nucleon structure.

In the parametrisation Eq. (1), *d*_{V}/*u*_{V} ∼ (1-*x*)^{bdV-buV} as *x* → 1, so if *b*_{uV} = *b*_{dV}, as predicted by the counting rules, then *d*_{V}/*u*_{V} → *k*, with *k* some constant. Indeed it is the constant *k* that many of the models try to predict. Moreover, as noted above, both CT14 and CJ15 assume *b*_{uV} = *b*_{dV} in their fits. However, while one may expect *b*_{uV} ≃ *b*_{dV} because of isospin symmetry, it is also reasonable to expect that exact equality will be broken by isospin breaking or electromagnetic effects. The sign of these effects is crucial: if *b*_{uV} > *b*_{dV} then *d*_{V}/*u*_{V} will become infinite as *x* → 1, while if *b*_{uV} < *b*_{dV}, as *x* → 1
*d*_{V}/*u*_{V} → 0. These two possibilities result in naive limits on the ratio ${F}_{2}^{n}/{F}_{2}^{p}$: if the sea quarks can be ignored at large *x*, then *d*_{V} ≫ *u*_{V}, ${F}_{2}^{n}/{F}_{2}^{p}\to 4$, while for *d*_{V} ≪ *u*_{V}
${F}_{2}^{n}/{F}_{2}^{p}\to 1/4$, giving for *x* → 1 the Nachtmann limits [58]

$$\begin{array}{c}\hfill \frac{1}{4}\le \frac{{F}_{2}^{n}}{{F}_{2}^{p}}\le 4\phantom{\rule{0.166667em}{0ex}}.\end{array}$$

16

To address these issues empirically, in Fig. 8 we compare the ratios *d*_{V}/*u*_{V}(*x*, *Q*^{2}) and ${F}_{2}^{n}(x,{Q}^{2})/{F}_{2}^{p}(x,{Q}^{2})$ at *Q*^{2} = 2 GeV^{2} as predicted by the various PDF sets. The neutron and proton structure functions ${F}_{2}^{n}(x,{Q}^{2})$ and ${F}_{2}^{p}(x,{Q}^{2})$ have been computed at NNLO accuracy with APFEL [59] using the FONLL-C general-mass scheme [60]. The arrows on the right hand side of each panel indicate the expectations from a representative set of nonperturbative models of nucleon structure: *S*U(6) [61] describes constituent quarks in the nucleon by SU(6) wave functions; CQM [62, 63] is the relativistic Constituent Quark Model in which a SU(6) symmetry breaking is assumed via a color hyperfine interaction between quarks; NJL [64] is a modified Nambu–Jona-Lasinio model in which confinement is simulated by eliminating unphysical thresholds for nucleon decay; pQCD [65] stands for a coloured quark and vector gluon model supplemented with leading order perturbative QCD; DSE1 and DSE2 [66] are two scenarios based on Dyson–Schwinger equations.

The ratios *d*_{V}/*u*_{V} (*left*) and ${F}_{2}^{n}/{F}_{2}^{p}$ (*right*) at *Q*^{2} = 2 GeV^{2} among various PDF sets, compared with the predictions of different nonperturbative models of nucleon structure

From Fig. 8, we see that in the region in which the valence quarks are constrained by experimental data, i.e. *x* ≲ 0.5, the predictions for both ratios from all the PDF sets are in reasonable agreement with each other within uncertainties, as might be expected. For *x* ≳ 0.5, the mutual consistency of PDF sets deteriorates rapidly, and a wide range of different behaviours is observed. This is a consequence of the reduced experimental information in this region: different PDF collaborations extrapolate to large *x* using different assumptions. For those sets with very weak assumptions on the PDF behaviour at large *x*, namely NNPDF3.0 and MMHT14, the uncertainties on the ratios expand rapidly, and at very large *x* there is no predictive power at all. For the two sets which assume that *d*_{V}/*u*_{V} → *k* at large *x*, namely CT14 and CJ15, uncertainties are inevitably much reduced and a value of *k* is predicted. ABM12 is different again, in that they find as a result of their fit that *b*_{dV} > *b*_{uV} at more than two standard deviations (see Table 3), so that *d*_{V}/*u*_{V} → 0 as *x* → 1, and an unrealistically small uncertainty band in a region where there are actually no data.

It follows that all the various model predictions displayed in Fig. 8 are compatible with the NNPDF3.0 and MMHT14 predictions, while ABM12 confirms the Chiral Quark Model but appears to rule out all the others. The CT14 and CJ15 sets favour values of *k* in the region 0 ≲ *k* ≲ 0.25, thus disfavouring the *SU*(6) prediction but unable to discriminate between the others. The preference for smaller values of *k* results in effect from a linear extrapolation of the downwards trend in the data region *x* ≲ 0.5. Not all the predictions respect the Nachtmann bound, Eq. (16).

In summary, in this work we have introduced a novel methodology to determine quantitatively the effective asymptotic behaviour of parton distributions, valid for any value of *x* and *Q*^{2}. For the first time, we have unambiguously identified the ranges in *x* and *Q*^{2} where the asymptotic regime sets in, allowing us to compare in detail perturbative and nonperturbative QCD predictions at large and small *x* with the results of modern global PDF fits.

Concerning the small-*x* region, we have found broad agreement between the results from PDF fits and the predictions from Regge theory for the behaviour of the valence quark distributions. For the singlet and gluon distributions, the agreement with Regge predictions is still only qualitative, due in part to the substantial scale dependence, as well as the limited experimental information available in that region. On the other hand, the perturbative QCD double asymptotic scaling predictions are in excellent agreement with the results of PDFs fits over a wide range of *Q*^{2}.

Concerning the large-*x* region, we have found that the predictions of the Brodsky–Farrar counting rules for the behaviour of the valence quark distributions are in broad agreement with the global fit results, within PDF uncertainties. For the sea and gluon distributions uncertainties are much larger, and the agreement is only qualitative. The scale dependence of the effective exponents based on global PDF fits is in excellent agreement with the perturbative QCD expectation from the cusp anomalous dimension in a wide range of *Q*^{2}. We have also compared the ratios *d*_{V}(*x*, *Q*^{2})/*u*_{V}(*x*, *Q*^{2}) and ${F}_{2}^{n}(x,{Q}^{2})/{F}_{2}^{p}(x,{Q}^{2})$ among PDF fits and with nonperturbative models of nucleon structure, but found that the interpretation of this comparison depends significantly on the assumptions built into the PDF parametrisation, to the extent that it is impossible at present to draw any firm conclusions.

We therefore conclude that, while the ancient wisdom of Regge theory and the Brodsky–Farrar counting rules seems to have some degree of truth, particularly in the valence quark sector, they are no substitute for the precise empirical PDF determinations provided by global analysis, and when used as constraints may lead to unrealistically accurate predictions in kinematic regions where there is no experimental data. Global PDF fits will always be hampered to some extent by the lack of data to constrain PDFs in extrapolation regions, and new measurements from the LHC and other facilities, such as JLab, are required to shed more light on the asymptotic behaviour of parton distributions at small and large *x*. The methodology presented in this work should find applications in future comparisons between different global PDF fits, and between PDF fits and nonperturbative models of nucleon structure.

J.R. is supported by an STFC Rutherford Fellowship ST/K005227/1 and by an European Research Council Starting Grant “PDF4BSM”. The work of J.R. and E.R.N. is supported by a STFC Rutherford Grant ST/M003787/1.

The accurate evaluation of the effective exponents *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}) through Eq. (4) is pivotal in our study. An analytic evaluation of Eq. (4) starting from the explicit PDF parametrisation in Eq. (1), though straightforward, has two main limitations. First, Eq. (1) holds only at the initial parametrisation scale ${Q}_{0}^{2}$; the form of Eq. (1) is rapidly washed out by DGLAP evolution, hence it cannot be used for the analytic computation of the effective exponents at ${Q}^{2}>{Q}_{0}^{2}$. Second, even at ${Q}^{2}={Q}_{0}^{2}$, only the best-fit parameters for the central PDF are provided, and, moreover, for some PDF fits not even a simple analytical parametrisation is used.

To overcome these difficulties, in this work we evaluate Eq. (4) numerically. To this purpose, PDFs in a suitable numerical format and an algorithm for the numerical computation of the logarithmic derivative of the PDF in Eq. (4) are necessary. The first requirement is fulfilled by LHAPDF6 [67], while the second is more delicate. The standard methods used to evaluate numerical derivatives, such as those based on a finite difference approximation or on a polynomial approximation of the function to be derived, see e.g. Sects. 5.7–5.9 in Ref. [68], are found to lead to unstable results. The reason is that PDFs available through LHAPDF6 are tabulated on a grid in (*x*, *Q*^{2}); the values of the PDFs off a grid node are then obtained by a cubic spline interpolation. This interpolation induces small fluctuations of the PDFs with respect to their *true* value, in particular in the small and large-*x* regions, where the grid tabulations are less dense. Such fluctuations are enhanced when the numerical derivative of the PDF is computed, especially if the value of the PDF is very small, thus spoiling the evaluation of Eq. (4).

We overcome this problem and perform the numerical derivative in Eq. (4) by means of a Savitzky–Golay smoothing filter [69]. The idea is the following. Assuming that a function *g*(*x*) is tabulated at *n* + 1 equally spaced intervals, *g*_{i} ≡ *g*(*x*_{i}), with *x*_{i} = *x*_{0} + *i*Δ for some constant sample spacing Δ = (*x*_{n} - *x*_{0})/*n* and *i* = - *n*/2, …, - 2, - 1, 0, 1, 2, …, *n*/2, the filter performs a least-squares fit with a polynomial of some degree *m* at each point, using an additional number *n*_{L} of points to the left and some number *n*_{R} of points to the right of each desired *x* value. The estimated derivative is then the derivative of the resulting fitted polynomial. The values of the parameters *x*_{0}, *x*_{n}, *n*, *n*_{L}, *n*_{R} and *m* are optimised for each flavour and PDF set, so that residual numerical instabilities are minimised.

The robustness of our numerical procedure can be validated by comparing it with an analytic evaluation of Eq. (4). For instance, we consider the MSTW08 NLO PDF set at *Q*^{2} = 1 GeV^{2}, and compute the central value of the effective exponents *α*_{fi}(*x*, *Q*^{2}) and *β*_{fi}(*x*, *Q*^{2}) both analytically and numerically. The relative difference between the two computations, defined as

$$\begin{array}{cc}& {R}^{{\mathit{\alpha}}_{{f}_{i}}}(x,{Q}^{2})=\frac{{\mathit{\alpha}}_{{f}_{i}}^{(\mathrm{num})}(x,{Q}^{2})-{\mathit{\alpha}}_{{f}_{i}}^{(\mathrm{ana})}(x,{Q}^{2})}{{\mathit{\alpha}}_{{f}_{i}}^{(\mathrm{ana})}(x,{Q}^{2})}\hfill \\ \hfill & {R}^{{\mathit{\beta}}_{{f}_{i}}}(x,{Q}^{2})=\frac{{\mathit{\beta}}_{{f}_{i}}^{(\mathrm{num})}(x,{Q}^{2})-{\mathit{\beta}}_{{f}_{i}}^{(\mathrm{ana})}(x,{Q}^{2})}{{\mathit{\beta}}_{{f}_{i}}^{(\mathrm{ana})}(x,{Q}^{2})},\hfill \end{array}$$

17

is displayed in Fig. 9 for the *u*_{V}, *d*_{V} and *g* PDFs. The agreement between the analytic (ana) and numeric (num) computation is excellent: relative differences are at the permille level, with the only exception of *α*_{g}(*x*, *Q*^{2}) around *x* ∼ 10^{-3}, where the input gluon PDF has a node.

1. Forte S, Watt G. Progress in the determination of the partonic structure of the proton. Ann. Rev. Nucl. Part. Sci. 2013;63:291–328. doi: 10.1146/annurev-nucl-102212-170607. [Cross Ref]

2. J. Rojo et al., The PDF4LHC report on PDFs and LHC data: results from Run I and preparation for Run II. J. Phys. G **42**, 103103 (2015). arXiv:1507.00556

3. R.D. Ball, Global parton distributions for the LHC Run II. Nuovo Cim. C **38**(4), 127 (2016). arXiv:1507.07891

4. J. Butterworth et al., PDF4LHC recommendations for LHC Run II. J. Phys. G **43**, 023001 (2016). arXiv:1510.03865

5. Jimenez-Delgado P, Reya E. Delineating parton distributions and the strong coupling. Phys. Rev. D. 2014;89(7):074049. doi: 10.1103/PhysRevD.89.074049. [Cross Ref]

6. **NNPDF** Collaboration, R.D. Ball et al., Parton distributions for the LHC Run II. JHEP **04**, 040 (2015). arXiv:1410.8849

7. Harland-Lang LA, Martin AD, Motylinski P, Thorne RS. Parton distributions in the LHC era: MMHT 2014 PDFs. Eur. Phys. J. C. 2015;75(5):204. doi: 10.1140/epjc/s10052-015-3397-6. [PMC free article] [PubMed] [Cross Ref]

8. S. Dulat, T.-J. Hou, J. Gao, M. Guzzi, J. Huston, P. Nadolsky, J. Pumplin, C. Schmidt, D. Stump, C.P. Yuan, New parton distribution functions from a global analysis of quantum chromodynamics. Phys. Rev. D **93**(3), 033006 (2016). arXiv:1506.07443

9. Alekhin S, Blumlein J, Moch S. The ABM parton distributions tuned to LHC data. Phys. Rev. D. 2014;89(5):054028. doi: 10.1103/PhysRevD.89.054028. [Cross Ref]

10. **ZEUS, H1** Collaboration, H. Abramowicz et al., Combination of measurements of inclusive deep inelastic *e*^{±}*p*${e}^{\pm}p$ scattering cross sections and QCD analysis of HERA data. Eur. Phys. J. C **75**(12), 580 (2015). arXiv:1506.06042

11. A. Accardi, L.T. Brady, W. Melnitchouk, J.F. Owens, N. Sato, Constraints on large-*x*$x$ parton distributions from new weak boson production and deep-inelastic scattering data. arXiv:1602.03154

12. V.N. Gribov, L.N. Lipatov, Deep inelastic e p scattering in perturbation theory. Sov. J. Nucl. Phys. **15**, 438–450 (1972). [Yad. Fiz.15,781(1972)]

13. Altarelli G, Parisi G. Asymptotic freedom in parton language. Nucl. Phys. B. 1977;126:298. doi: 10.1016/0550-3213(77)90384-4. [Cross Ref]

14. Y.L. Dokshitzer, Calculation of the structure functions for deep inelastic scattering and e + $+$ e - $-$ annihilation by Perturbation Theory in quantum chromodynamics. Sov. Phys. JETP **46**, 641–653 (1977). [Zh. Eksp. Teor. Fiz.73,1216(1977)]

15. Regge T. Introduction to complex orbital momenta. Nuovo Cim. 1959;14:951. doi: 10.1007/BF02728177. [Cross Ref]

16. Brodsky SJ, Farrar GR. Scaling laws at large transverse momentum. Phys. Rev. Lett. 1973;31:1153–1156. doi: 10.1103/PhysRevLett.31.1153. [Cross Ref]

17. Roberts RG. The Structure of the Proton: Deep Inelastic Scattering. Cambridge: Cambridge University Press; 1994.

18. Devenish R, Cooper-Sarkar A. Deep Inelastic Scattering. Oxford: Oxford University Press; 2004.

19. Glazov A, Moch S, Radescu V. Parton distribution uncertainties using smoothness prior. Phys. Lett. B. 2011;695:238–241. doi: 10.1016/j.physletb.2010.11.025. [Cross Ref]

20. **NNPDF** Collaboration, L. Del Debbio, S. Forte, J.I. Latorre, A. Piccione, J. Rojo, Unbiased determination of the proton structure function F(2)**p with faithful uncertainty estimation. JHEP **03**, 080 (2005). arXiv:hep-ph/0501067

21. **NNPDF** Collaboration, L. Del Debbio, S. Forte, J.I. Latorre, A. Piccione, J. Rojo, Neural network determination of parton distributions: the nonsinglet case. JHEP **03**, 039 (2007). arXiv:hep-ph/0701127

22. Lepage GP, Brodsky SJ. Exclusive processes in perturbative quantum chromodynamics. Phys. Rev. D. 1980;22:2157. doi: 10.1103/PhysRevD.22.2157. [Cross Ref]

23. J. Polchinski, M.J. Strassler, Hard scattering and gauge/string duality. Phys. Rev. Lett. **88**, 031601 (2002). arXiv:hep-th/0109174 [PubMed]

24. V.A. Matveev, R.M. Muradyan, A.N. Tavkhelidze, Automodelity in strong interactions. Lett. Nuovo Cim. **5S2**, 907–912 (1972). [Lett. Nuovo Cim.5,907(1972)]

25. A. Deur, S.J. Brodsky, G.F. de Teramond, On the interface between perturbative and nonperturbative QCD. Phys. Lett. B **757**, 275–281 (2016). arXiv:1601.06568

26. Nocera ER. Small- and large-*x*$x$ nucleon spin structure from a global QCD analysis of polarized parton distribution functions. Phys. Lett. B. 2015;742:117–125. doi: 10.1016/j.physletb.2015.01.021. [Cross Ref]

27. Martin AD, Stirling WJ, Thorne RS, Watt G. Parton distributions for the LHC. Eur. Phys. J. C. 2009;63:189–285. doi: 10.1140/epjc/s10052-009-1072-5. [Cross Ref]

28. **NNPDF** Collaboration, R.D. Ball, L. Del Debbio, S. Forte, A. Guffanti, J.I. Latorre, A. Piccione, J. Rojo, M. Ubiali, A Determination of parton distributions with faithful uncertainty estimation. Nucl. Phys. B **809**, 1–63 (2009). arXiv:0808.1231. [Erratum: Nucl. Phys.B816,293(2009)]

29. Brodsky SJ, Ma B-Q. The quark/anti-quark asymmetry of the nucleon sea. Phys. Lett. B. 1996;381:317–324. doi: 10.1016/0370-2693(96)00597-7. [Cross Ref]

30. R. Gauld, J. Rojo, L. Rottoli, J. Talbert, Charm production in the forward region: constraints on the small-x gluon and backgrounds for neutrino astronomy. JHEP **11**, 009 (2015). arXiv:1506.08025

31. **PROSA** Collaboration, O. Zenaiev et al., Impact of heavy-flavour production cross sections measured by the LHCb experiment on parton distribution functions at low x. Eur. Phys. J. C **75**(8), 396 (2015). arXiv:1503.04581

32. S.P. Jones, A.D. Martin, M.G. Ryskin, T. Teubner, Exclusive *J*/*ψ*$J/\mathit{\psi}$ and Υ$\mathrm{{\rm Y}}$ photoproduction and the low *x*$x$ gluon. J. Phys. G **43**(3), 035002 (2016). arXiv:1507.06942

33. De Rujula A, Glashow SL, Politzer HD, Treiman SB, Wilczek F, Zee A. Possible nonRegge behavior of electroproduction structure functions. Phys. Rev. D. 1974;10:1649. doi: 10.1103/PhysRevD.10.1649. [Cross Ref]

34. Ball RD, Forte S. Double asymptotic scaling at HERA. Phys. Lett. B. 1994;335:77–86. doi: 10.1016/0370-2693(94)91561-X. [Cross Ref]

35. Ball RD, Forte S. A direct test of perturbative QCD at small x. Phys. Lett. B. 1994;336:77–79. doi: 10.1016/0370-2693(94)00956-2. [Cross Ref]

36. Forte S, Ball RD. Universality and scaling in perturbative QCD at small x. Acta Phys. Polon. B. 1995;26:2097–2134.

37. Ralston JP. Pocket partonometer. Phys. Lett. B. 1986;172:430–434. doi: 10.1016/0370-2693(86)90283-2. [Cross Ref]

38. Korchemsky GP. Asymptotics of the Altarelli–Parisi–Lipatov evolution kernels of parton distributions. Mod. Phys. Lett. A. 1989;4:1257–1276. doi: 10.1142/S0217732389001453. [Cross Ref]

39. Albino S, Ball RD. Soft resummation of quark anomalous dimensions and coefficient functions in MS-bar factorization. Phys. Lett. B. 2001;513:93–102. doi: 10.1016/S0370-2693(01)00742-0. [Cross Ref]

40. Kodaira J, Trentadue L. Summing soft emission in QCD. Phys. Lett. B. 1982;112:66. doi: 10.1016/0370-2693(82)90907-8. [Cross Ref]

41. Ermolaev BI, Greco M, Troyan SI. Intercepts of the nonsinglet structure functions. Nucl. Phys. B. 2001;594:71–88. doi: 10.1016/S0550-3213(00)00647-7. [Cross Ref]

42. Ermolaev BI, Greco M, Troyan SI. Overview of the spin structure function g(1) at arbitrary x and Q2. Riv. Nuovo Cim. 2010;33:57–122.

43. Abarbanel HDI, Goldberger ML, Treiman SB. Asymptotic properties of electroproduction structure functions. Phys. Rev. Lett. 1969;22:500–502. doi: 10.1103/PhysRevLett.22.500. [Cross Ref]

44. Landshoff PV, Polkinghorne JC, Short RD. A nonperturbative parton model of current interactions. Nucl. Phys. B. 1971;28:225–239. doi: 10.1016/0550-3213(71)90375-0. [Cross Ref]

45. L.N. Lipatov, Reggeization of the vector meson and the vacuum singularity in nonabelian gauge theories. Sov. J. Nucl. Phys. **23** (1976) 338–345. [Yad. Fiz.23,642(1976)]

46. Fadin VS, Kuraev EA, Lipatov LN. On the Pomeranchuk singularity in asymptotically free theories. Phys. Lett. B. 1975;60:50–52. doi: 10.1016/0370-2693(75)90524-9. [Cross Ref]

47. E.A. Kuraev, L.N. Lipatov, V.S. Fadin, Multi-Reggeon processes in the Yang–Mills theory. Sov. Phys. JETP **44**, 443–450 (1976). [Zh. Eksp. Teor. Fiz.71,840(1976)]

48. E.A. Kuraev, L.N. Lipatov, V.S. Fadin, The Pomeranchuk singularity in nonabelian gauge theories. Sov. Phys. JETP **45**, 199–204 (1977). [Zh. Eksp. Teor. Fiz.72,377(1977)]

49. Fadin VS, Lipatov LN. BFKL pomeron in the next-to-leading approximation. Phys. Lett. B. 1998;429:127–134. doi: 10.1016/S0370-2693(98)00473-0. [Cross Ref]

50. Altarelli G, Ball RD, Forte S. Factorization and resummation of small x scaling violations with running coupling. Nucl. Phys. B. 2002;621:359–387. doi: 10.1016/S0550-3213(01)00563-6. [Cross Ref]

51. Altarelli G, Ball RD, Forte S. An anomalous dimension for small x evolution. Nucl. Phys. B. 2003;674:459–483. doi: 10.1016/j.nuclphysb.2003.09.040. [Cross Ref]

52. Ciafaloni M, Colferai D, Salam GP, Stasto AM. Renormalization group improved small x Green’s function. Phys. Rev. D. 2003;68:114003. doi: 10.1103/PhysRevD.68.114003. [Cross Ref]

53. Altarelli G, Ball RD, Forte S. Perturbatively stable resummed small x evolution kernels. Nucl. Phys. B. 2006;742:1–40. doi: 10.1016/j.nuclphysb.2006.01.046. [Cross Ref]

54. Ciafaloni M, Colferai D, Salam GP, Stasto AM. A Matrix formulation for small-x singlet evolution. JHEP. 2007;08:046. doi: 10.1088/1126-6708/2007/08/046. [Cross Ref]

55. Brodsky SJ, Burkardt M, Schmidt I. Perturbative QCD constraints on the shape of polarized quark and gluon distributions. Nucl. Phys. B. 1995;441:197–214. doi: 10.1016/0550-3213(95)00009-H. [Cross Ref]

56. W. Melnitchouk, A.W. Thomas, Neutron / proton structure function ratio at large x. Phys. Lett. B **377** (1996) 11–17, arXiv:nucl-th/9602038

57. Holt RJ, Roberts CD. Distribution functions of the nucleon and pion in the valence region. Rev. Mod. Phys. 2010;82:2991–3044. doi: 10.1103/RevModPhys.82.2991. [Cross Ref]

58. Nachtmann O. Inequalities for structure functions of deep inelastic lepton-nucleon scattering giving tests of basic algebraic structures. Nucl. Phys. B. 1972;38:397–417. doi: 10.1016/0550-3213(72)90317-3. [Cross Ref]

59. Bertone V, Carrazza S, Rojo J. APFEL: a PDF evolution library with QED corrections. Comput. Phys. Commun. 2014;185:1647–1668. doi: 10.1016/j.cpc.2014.03.007. [Cross Ref]

60. Forte S, Laenen E, Nason P, Rojo J. Heavy quarks in deep-inelastic scattering. Nucl. Phys. B. 2010;834:116–162. doi: 10.1016/j.nuclphysb.2010.03.014. [Cross Ref]

61. F.E. Close, *An Introduction to Quarks and Partons.* (Academic Press, London, 1979)

62. Close FE. Nu w(2) at small omega’ and resonance form-factors in a quark model with broken su(6) Phys. Lett. B. 1973;43:422–426. doi: 10.1016/0370-2693(73)90389-4. [Cross Ref]

63. Carlitz RD. SU(6) symmetry breaking effects in deep inelastic scattering. Phys. Lett. B. 1975;58:345. doi: 10.1016/0370-2693(75)90670-X. [Cross Ref]

64. Cloet IC, Bentz W, Thomas AW. Nucleon quark distributions in a covariant quark-diquark model. Phys. Lett. B. 2005;621:246–252. doi: 10.1016/j.physletb.2005.06.065. [Cross Ref]

65. Farrar GR, Jackson DR. Pion and nucleon structure functions near x=1. Phys. Rev. Lett. 1975;35:1416. doi: 10.1103/PhysRevLett.35.1416. [Cross Ref]

66. Roberts CD, Holt RJ, Schmidt SM. Nucleon spin structure at very high-x. Phys. Lett. B. 2013;727:249–254. doi: 10.1016/j.physletb.2013.09.038. [Cross Ref]

67. Buckley A, Ferrando J, Lloyd S, Nordström K, Page B, Rüfenacht M, Schönherr M, Watt G. LHAPDF6: parton density access in the LHC precision era. Eur. Phys. J. C. 2015;75:132. doi: 10.1140/epjc/s10052-015-3318-8. [Cross Ref]

68. W.H. Press, B.P. Flannery, S.A. Teukolsky, W.T. Vetterling, *Numerical Recipes in Fortran 77: The Art of Scientific Computing,* 2nd edn. (Cambridge University Press, Cambridge, 1992). ISBN-13: 9780521430647, 10: 052143064X

69. A. Savitzky, M. Golay, Smoothing and differentiation of data by simplified least squares procedures. Anal. Chem. **36**(8), 1627–1639 (1964). http://dx.doi.org/10.1021/ac60214a047

Articles from Springer Open Choice are provided here courtesy of **Springer**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |