The asymptotic behaviour of parton distributions at small and large x.

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 [Formula: see text], where the values of [Formula: see text] can be deduced from Regge theory, while at asymptotically large values of x the PDFs behave as [Formula: see text], where the values of [Formula: see text] can be deduced from the Brodsky-Farrar quark counting rules. We critically examine these claims by extracting the exponents [Formula: see text] and [Formula: see text] 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 [Formula: see text] at large x, and caution against unrealistic uncertainty estimates due to overconstrained parametrisations.

Introduction

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, , is set at an initial scale , and is then evolved to any other scale 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 vanishes in the elastic limit should also be implicit in the parametrisation. Usually, the following ansatz is adoptedwhere x is the parton momentum fraction and i denotes a given quark flavour (or flavour combination) or the gluon, and is a smooth function which remains finite both when and . The normalisation fractions , the exponents and , and the set of parameters are then determined from the data. Some of the 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] predictswhile the Brodsky–Farrar quark counting rules [16] predictsee also Refs. [17, 18], and references therein. Both Regge theory and the counting rules provide numerical predictions for the values of the exponents and . In Eq. (1), the small- and large-x power-law behaviours are matched at intermediate x values through the function . 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 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 . In principle, the optimal values of should be chosen at the interface between perturbative and nonperturbative hadron dynamics, . It has been shown [25] that GeV by matching the high- and low- behaviour of the strong coupling as predicted respectively by its renormalisation group equation in the 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 , 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].

The effective exponents

In this paper we will compute the effective exponents and , which, when , are asymptotically equal to the exponents and of the input PDF parametrisation Eq. (1). Specifically, we defineso that, at the input parametrisation scale ,andsince 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 tend to zero very quickly at large x, we expect that the definitions Eq. (4) and can be used to accurately determine the asymptotic behaviour of any given PDF .

In order to test this assertion, we have used Eq. (4) to compute the effective asymptotic exponents and for the MSTW08 NLO PDF set [27] (see Appendix 1 for details). Results at GeV, which coincides with the input parametrisation scale , are shown in Fig. 1 for the up valence quark, , the down valence quark, , and the gluon, , PDFs. They are compared to the corresponding fitted exponents and , to which they are expected to approach asymptotically. In Table 1 we show the numerical values computed respectively at and , and we again compare them with the corresponding fitted exponents and .

Fig. 1

The effective exponents (left) and (right) computed with Eq. (4). Results are shown at  GeV for , and g for the MSTW08 NLO PDF set. The shaded horizontal bands indicate the fitted values of the exponents (left) and (right) and their uncertainties. Numerical results at and are collected in Table 1

Table 1

The effective exponents and at  GeV and and computed for the MSTW08 NLO PDF set with Eq. (4), compared to the corresponding fitted exponents and

f_i	α_f_i(x_a,Q^2)	a_f_i	β_f_i(x_b,Q^2)	b_f_i
u_V	+0.29±0.01	+0.291-0.013+0.019	+3.11±0.04	+3.243-0.039+0.062
d_V	+1.02±0.11	+0.968-0.110+0.110	+5.67±0.47	+5.944-0.530+0.510
g	-0.30±0.37	-0.428-0.057+0.066	+2.95±0.39	+3.023-0.360+0.430

From Fig. 1 and Table 1 it is clear that both at and at have converged to the fitted values of and within PDF uncertainties. In addition, by examining the x dependence of and , 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 , we may use them to study the 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, and , the total quark sea, , 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 is a multi-layer feed-forward neural network (also known as perceptron). The power-law term in Eq. (1) is treated as a preprocessing factor that optimises the minimisation process: the exponents and are chosen for each Monte Carlo replica at random in a given range determined iteratively.

MMHT14 The PDFs parametrised are the valence distributions and , the total sea S, the sea asymmetry , the total and valence strange distributions and and the gluon g. The function is taken to be a linear combination of Chebyshev polynomials. The exponents and are fitted, except for .

CT14 The PDFs parametrised are the valence distributions and , the sea quark distributions and , the total strangeness and the gluon g. It is assumed that . The function is a linear combination of Bernstein polynomials. The exponents and are parameters of the fit, but not all of them are free: specifically, it is assumed that , so that as , with k a constant, and that as , which requires .

ABM12 The PDFs parametrised are the valence distributions and , the sea distributions and s, the sea asymmetry and the gluon g. It is assumed that . The function has the form , where is a function of x; for s, . The exponents and are parameters of the fit, except for the condition .

CJ15 The PDFs parametrised are the valence distributions and , the light antiquark sea, , the light antiquark ratio , the total strangeness and the gluon g. It is assumed that . The function is provided by the polynomial for all the distributions except the light antiquark ratio and the total strangeness. Specifically, is parametrised with a simple polynomial which ensures that as , , while it is assumed that ; , and are parameters of the fit. A small admixture of is added to so that as , 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 . 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, 4 we compare both the PDFs and the corresponding effective exponents and for the NNPDF3.0, MMHT14 and CT14 sets at GeV. For NNPDF3.0, PDF uncertainties are computed as 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 for , and S irrespective of the PDF set considered. For the gluon, convergence is achieved at smaller values of x, , at least for MMHT14 for which has an oscillation in the region . Note that at 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 in most cases. The exception is from MMHT14, which exhibits an oscillation in the region .

Fig. 2

The effective exponents (left) and (right), Eq. (4), for the up valence (top) and down valence (bottom) PDFs, as a function of x at  GeV, together with the corresponding PDFs. Results are shown for the NNPDF3.0, CT14 and MMHT14 NNLO PDF sets. The arrows indicate the prediction from Regge theory (Regge) and Brodsky–Farrar quark counting rules (CR)

Fig. 3

Same as Fig. 2 for the sea PDF

Fig. 4

Same as Fig. 2 for the gluon PDF

To the best of our knowledge, this is the first time that the onset of an asymptotic regime in the effective PDF exponents and 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 and 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 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 dependence of the effective exponents. To this end, in Figs. 5, 6 we show the effective exponents and as functions of at fixed values of x in the asymptotic region: and , respectively. We show results for the valence distributions and , the total quark singlet and the gluon. From these plots we can see that as increases the effective exponents become less sensitive to and tend to converge to a finite value asymptotically. This feature is broadly independent of x when x is sufficiently small or large, roughly and . The only exception is again for MMHT14.

Fig. 5

The effective exponents (left) and (right), Eq. (4), for the up (top) and down valence (bottom) PDFs, as a function of at and , respectively, for the NNPDF3.0, CT14 and MMHT14 NNLO sets. At large x, the perturbative QCD prediction Eq. (15) is also displayed for CT14

Fig. 6

Same as Fig. 5 for the quark singlet and the gluon g. For (left plots) we also show the DAS predictions, Eq. (10)

At small x, the dependence of the effective exponents illustrates the transition from a low- region, where PDFs are determined from nonperturbative dynamics, to a high- region, where PDFs are dominated by perturbative QCD evolution. Indeed, as and , 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 and the singlet sector grows aswhere we have definedand the double scaling variablesThe parameters and define the formal boundaries of the asymptotic region, and are normalisation constants, and 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 : at large but fixed one hasNote that both and 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 ,  GeV, and  GeV. The agreement between the expectation from DAS and results from the global fits is excellent at  GeV for both the quark singlet and the gluon.

At large x, the 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 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 formwhere and can be computed perturbatively: for example at NLO

Table 2

The values of the small-x effective exponent computed at  GeV and  GeV at , compared to the values of predicted by Regge theory (and resummation of double logarithms). For the quark sea S and the gluon g we indicate the prediction of the soft Pomeron (and the NLLx perturbative result)

f_i	Q^2	α_f_i(x_a,Q^2)					a_f_i
	(GeV^2)	NNPDF3.0	CT14	MMHT14	ABM12	CJ15
u_V	2.0	+0.48±0.11	+0.72±0.12	+0.65±0.06	+0.76±0.07	+0.61±0.01	+0.5
	10.0	+0.46±0.09	+0.66±0.09	+0.61±0.04	+0.70±0.04	+0.60±0.01	(0.63)
d_V	2.0	+0.41±0.11	+0.73±0.12	+0.79±0.06	+1.39±0.10	+1.11±0.03	+0.5
	10.0	+0.41±0.11	+0.66±0.07	+0.70±0.04	+0.91±0.08	+0.95±0.05	(0.63)
S	2.0	-0.14±0.06	-0.15±0.05	-0.09±0.04	-0.16±0.02	-0.18±0.03	-0.08
	10.0	-0.18±0.04	-0.20±0.05	-0.15±0.04	-0.19±0.01	-0.14±0.02	(-0.2)
g	2.0	-0.16±0.63	+0.06±0.31	-0.79±0.43	+0.18±0.10	+0.08±0.03	-0.08
	10.0	-0.20±0.46	-0.15±0.15	-0.29±0.09	-0.15±0.01	-0.14±0.01	(-0.2)

with coefficients [40]It follows [39] that, if as at a scale , with q either the quark singlet, , or one of the quark valence distributions, or , then this asymptotic behaviour persists at higher scales withGiven our definition Eq. (4) and the asymptotic behaviour Eq. (6) at large x, as one hasThe behaviour predicted by Eq. (15) is displayed for and 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 in Eq. (15) to match the central values obtained from CT14 at  GeV. The agreement between Eq. (15) and the dependence of the large-x effective exponents derived from the PDF fit is excellent. A slight deterioration only appears at small values of due to missing higher order corrections. Similar conclusions can be derived for other PDF sets when the value of in Eq. (15) is assigned consistently.

Comparison with nonperturbative predictions

We now discuss how our findings compare with the expectations from Regge theory and the Brodsky–Farrar quark counting rules. In Tables 2, 3 we show the values of the effective exponents for the NNPDF3.0, CT14, MMHT14, ABM12 and CJ15 PDF sets, computed at and ( for S) at  GeV and  GeV. We also include the values predicted by Regge theory and the Brodsky–Farrar quark counting rules.

Table 3

Same as Table 2 for the large-x effective exponent at (for , and g) and (for S). The values of the exponent predicted by Brodsky–Farrar quark counting rules are also shown

f_i	Q^2	β_f_i(x_b,Q^2)					b_f_i
	(GeV^2)	NNPDF3.0	CT14	MMHT14	ABM12	CJ15
u_V	2.0	+2.94±0.52	+3.11±0.28	+3.37±0.07	+3.38±0.06	+3.50±0.01	∼3
	10.0	+3.30±0.69	+3.38±0.29	+3.62±0.07	+3.61±0.05	+3.78±0.01
d_V	2.0	+3.03±1.96	+3.27±0.37	+2.05±0.59	+4.72±0.43	+3.42±0.06	∼3
	10.0	+3.23±1.88	+3.52±0.36	+2.29±0.59	+4.92±0.42	+3.68±0.05
S	2.0	+6.86±7.25	+6.41±1.22	+8.19±0.68	+8.16±0.38	+7.73±0.18	∼7
	10.0	+6.76±6.71	+6.91±1.14	+6.83±0.88	+8.51±0.38	+8.15±0.18
g	2.0	+2.95±1.25	+5.08±2.18	+1.65±0.23	+4.18±0.06	+6.11±0.33	∼5
	10.0	+3.25±0.98	+5.13±0.51	+2.24±0.23	+4.44±0.06	+4.91±0.33

At small x, Regge theory predicts with a -independent exponent, related to the intercept of the corresponding Regge trajectory. For valence quark distributions, a value of is derived from the non-singlet Regge trajectory intercept . Perturbative calculations which resum the double logarithms of x give a similar value [41, 42]. For the gluon distribution, a value of close to the singlet Pomeron trajectory 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 . Attempts to compute the Pomeron intercept perturbatively by solution of the fixed coupling LLx BFKL equation [45-48] give . 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 . For the total sea distribution, the value of should be similar for large enough to , due to the dominance of the process 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 and 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 and , 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 , where 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 , 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, and thus ; for a gluon, and ; for a sea quark, and ; 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), and it is reasonable to make the comparison at a low scale [25].

Fig. 7

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 and 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, , and of neutron to proton structure functions, . These ratios are particularly interesting because while all PDFs vanish at , their ratio does not necessarily do so, and thus it is a useful discriminator among models of nucleon structure.

In the parametrisation Eq. (1), as , so if , as predicted by the counting rules, then , 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 in their fits. However, while one may expect 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 then will become infinite as , while if , as . These two possibilities result in naive limits on the ratio : if the sea quarks can be ignored at large x, then , , while for , giving for the Nachtmann limits [58]To address these issues empirically, in Fig. 8 we compare the ratios and at  GeV as predicted by the various PDF sets. The neutron and proton structure functions and 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:  [61] describes constituent quarks in the nucleon by wave functions; CQM [62, 63] is the relativistic Constituent Quark Model in which a 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.

Fig. 8

The ratios (left) and (right) at GeV 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. , the predictions for both ratios from all the PDF sets are in reasonable agreement with each other within uncertainties, as might be expected. For , 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 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 at more than two standard deviations (see Table 3), so that as , 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 , 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 . Not all the predictions respect the Nachtmann bound, Eq. (16).

Conclusions and outlook

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 . For the first time, we have unambiguously identified the ranges in x and 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 .

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 . We have also compared the ratios and 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.