NLO PDFs from the ABMP16 fit

We perform a global fit of parton distribution functions (PDFs) together with the strong coupling constant αs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _s$$\end{document} and the quark masses mc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_c$$\end{document}, mb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_b$$\end{document} and mt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_t$$\end{document} at next-to-leading order (NLO) in QCD. The analysis applies the MS¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {MS}}\, $$\end{document} renormalization scheme for αs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _s$$\end{document} and all quark masses. It is performed in the fixed-flavor number scheme for nf=3,4,5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_f=3, 4, 5$$\end{document} and uses the same data as the previous ABMP16 fit at next-to-next-to-leading order (NNLO). The new NLO PDFs complement the set of ABMP16 PDFs and are to be used consistently with NLO QCD predictions for hard scattering processes. At NLO we obtain the value αs(nf=5)(MZ)=0.1191±0.0011\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _s^{(n_f=5)}(M_Z) = 0.1191 \pm 0.0011$$\end{document} compared to αs(nf=5)(MZ)=0.1147±0.0008\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _s^{(n_f=5)}(M_Z) = 0.1147 \pm 0.0008$$\end{document} at NNLO.

Parton distribution functions (PDFs) are an indispensable ingredient in theory predictions for hadronic scattering processes within perturbative QCD. Currently, the state-of-art calculations for many standard-candle processes at the Large Hadron Collider (LHC) and elsewhere are based on the QCD corrections up to the next-to-next-to-leading order (NNLO) in the strong coupling constant α s [1]. In order to match this theoretical accuracy the PDFs and other input parameters such as α s and the quark masses m c , m b and m t also have to be determined at the same order of perturbation theory, that is with account of the NNLO QCD corrections. In many instances, however, the Wilson coefficient functions or hard partonic scattering cross sections are known to the next-to-leading order (NLO) only. This concerns in particular Monte-Carlo studies at the LHC. Then, to meet the consistency requirements, NLO PDFs and the respective NLO values for α s and the heavy-quark masses are to be used. The NLO fit of PDFs is therefore of immediate practical use and also provides a very good consistency check of the perturbative stability of QCD calculations.
In this article we describe the NLO version of the recent ABMP16 PDF fit, i.e., the NLO analysis, which applies the MS scheme for α s and all heavy-quark masses. It uses the a e-mail: sven-olaf.moch@desy.de same data, their uncertainty treatment and the general theoretical framework, e.g., the fixed-flavor number scheme for n f = 3, 4 and 5, as in the previous fit of the ABMP16 PDFs NNLO. The only difference resides in the order of the perturbative corrections to the QCD evolution equations and for the Wilson coefficients, which are now limited to NLO accuracy. Due to the obvious correlations of the various parameters in the PDFs with the value of α s and those of the quark masses m c , m b and m t , all quantities are extracted simultaneously from the global fit following our previous analyses [2][3][4]. The article discusses in detail the differences in their determinations at NLO and NNLO accuracy. Specific attention is paid to the treatment of power corrections in the description of deep-inelastic scattering (DIS) data, i.e. higher-twist effects which are relevant beyond the leading twist collinear factorization approximation. The final fit results are made available as data grids for use with the LHAPDF library (version 6) [5] and the features of the various grids are briefly discussed.
The values of χ 2 obtained in the present analysis for various data sets are listed in Table 1 in comparison with the earlier ones for the NNLO ABMB16 fit. The overall quality of the data description does not change dramatically between the NLO and the NNLO versions, where the former features a somewhat bigger total value of χ 2 . Of course, the theoretical description at NNLO accuracy comes with a significantly reduced theoretical uncertainty due to variations of the factorization and renormalization scales compared to the NLO one. Nevertheless, for specific scattering reactions the NNLO corrections are crucial for the respective data sets. This holds in particular for the c-quark and, to a lesser extent, for b-quark production in DIS and for hadronic t-quark pair-production, which constrain the heavy-quark masses and which are fitted together with α s simultaneously with the PDFs. The theoretical description at NNLO accuracy is also essential for the parameters of the higher (dynamical) twist, which contribute additively to the leading twist. The x-dependent twist-four contributions to the longitudinal and transverse DIS cross Table 1 The values of χ 2 obtained in the present analysis at NLO for the data on inclusive DIS, the DY process, and on heavy-quark production in comparison with the ones of the ABMP16 fit at NNLO [4]. a This value corrects a misprint in Table V [4]). This provides a better consistency between the PDF sets obtained with different theoretical accuracy. At the same time the uncertainties in the normalization and higher twist parameters are computed in the same way as in the NNLO fit [4], by propagation of the ones in experimental data and simultaneosly with other fit parameters in order to take into account their correlations and, therefore, provide a consistent uncertainty treatment in the NLO and NNLO fits. Therefore, the uncertainties obtained for the data normalization and the twist-four contributions at NLO are only marginally different from those reported at NNLO in Tables I and VIII in Ref. [4]. Closer inspection of the χ 2 -values in Table 1 reveals the largest differences between NLO and NNLO for the fixedtarget DIS data, which can be explained by the kinematic coverage of this data sample, which is predominantly in the low-Q 2 region. For the Drell-Yan (DY) data the impact of the NNLO QCD corrections is less pronounced and the corresponding improvement in the value of χ 2 is small. For the heavy-quark production data on the other hand the trend is not uniform, i.e., for some data sets the χ 2 -values at NLO are larger and vice versa for others. In this context it is worth noting, that to a certain extent the impact of missing NNLO terms in the NLO fit is compensated by tuning the values of heavy-quark masses.
The ABMP16 PDF sets at NLO and NNLO are compared in Fig. 1 for the case of n f = 3 flavors at the scale μ = 3 GeV. Both sets are based on the same flexible parametrization used in Ref. [4]. For the gluon PDF, we see in Fig. 1a that the NLO PDFs are larger by about 15% in the small-x and the large-x region, i.e., for x 10 −4 and x 0.3, respectively. In these kinematic regions for example the DIS coefficient functions receive systematically large corrections at higher orders, which need to be compensated by the gluon PDF if the fit is performed at NLO accuracy. The u-quark PDF in Fig. 1b does not show any big changes, except for large x 0.6, while the d-quark PDF Fig. 1c at NLO is smaller in the entire range x 10 −1 and decreasing more than 20% for x 10 −4 . A similar observation holds for the strange sea displayed in Fig. 1d, which is smaller by even 50% for x 10 −4 at NLO, however, the PDF uncertainties for this quantity are correspondingly larger. On the other hand, the non-strange sea in Fig. 1e does not show big relative differences between NLO and NNLO. There is only a slight decrease of the NLO result by 5% to 10% for x 10 −2 . The small-x sea iso-spin asymmetryd −ū at NLO goes lower than the NNLO one, as The value of α s (M Z ) at the scale of the Z -boson mass M Z obtained in the present analysis at NLO is larger than the one obtained in the NNLO variant of the ABMP16 fit [4], the relative difference amounting to about 4%, which is well comparable to the estimated margin due to variations of the factorization and renormalization scales, cf. Ref. [6]. In the scheme with n f = 5 light flavors we find as listed in the first line of Table 2 together with the kinematic cuts imposed on the DIS data. The description of DIS data at those low values of Q 2 and W for the invariant mass of the hadronic system, where W 2 = M 2 P + Q 2 (1 − x)/x with the proton mass M P , requires modeling of the higher-twist terms. This has been discussed extensively in the ABMP16 analyses at NNLO [4]. Following the theoretical framework there, the fitted twist-four contributions to the longitudinal and transverse DIS cross sections have been used to determine the value of α NLO s (M Z ) in Eq. (1). Alternatively, one can impose cuts both on Q 2 and W 2 to eliminate data from the kinematic regions most sensitive to the higher-twist terms. Then, the fit can be performed with all higher-twist terms set to zero and the results are shown in Table 2. These variants of the fit with substantially higher cuts on Q 2 and W 2 and higher-twist terms set to zero display very good stability of the value of α s (M Z ), both at NLO and NNLO, and therefore very good consistency of the chosen approach.
Finally, in one of the variants of the present analysis we impose the low cuts on Q 2 and W 2 from the first line of Table 2, while fitting also the twist-four contributions. This gives an improvement in the value of χ 2 equal to 86 and α NLO s (M Z ) = 0.1227 ± 0.0011, which is a slightly larger value than those quoted in Table 2 for the fits with higher cuts on Q 2 and W 2 . The magnitude of these shifts in α s (M Z )  Nevertheless, the observed difference between α NLO s (μ) and α NNLO s (μ) is quite essential, particularly at small scales μ, where the NLO and NNLO results differ by more than 10%, as illustrated in Fig. 2 for a wide range of scales. This difference is to a great extent responsible for the perturbative stability of QCD calculations at the hard scales currently probed in scattering processes at colliders. Asymptotic freedom in QCD, i.e. stability of theoretical predictions under higher order perturbative corrections requires very large scales. On the other hand, for realistic kinematics including experiments at the LHC a consistent setting of α s (M Z ) is very important to achieve sensible theoretical predictions. 1 Table 3 The In this context it is worth to mention the conventional choice α NLO , which is adopted as a part of the PDF4LHC recommendations [9] and employed in the CT14 [10] and NNPDF [11] PDF fits. Under this assumption the value of α s obtained at NLO is very close to the NNLO one in a wide range of scales, as shown in Fig. 2. As a result, such an approach has significant limitations when studying the convergence of the perturbative expansion, since the NLO predictions obtained with these PDF sets might be very similar to the NNLO ones simply due to the convention used.
The values for the heavy-quark masses obtained in the NLO and NNLO variants of the ABMP16 analysis are given in Table 3. The comparison of the difference between these values indicates again the limitations of the NLO approximation. The shifts are more significant for m c and m t , since the NNLO corrections are absolutely essential in order to achieve a good description of the data on DIS c-quark production and t-quark hadro-production. On the other hand, the data on DIS b-quark production are less precise, therefore the value of m b extracted from the fit suffers from the larger uncertainties and is less sensitive to the impact of the NNLO corrections, cf. Table 3.
In addition, the value of m t demonstrates a strong correlation with the value of α s , since the Born cross section for ttproduction is proportional to α 2 s , so that changes in the value of α s induce shifts in the fitted value of m t [3,4]. This is quantified in Fig. 3, where the values of m t determined in variants of the present analysis with fixed values of α s demonstrate a nearly linear dependence on α s . It is interesting to note, though, that at NLO this dependence is somewhat shallower than for the similar fit at NNLO since due to important missing QCD corrections of O(α 4 s ) in the hadronic tt-production cross section at NLO, σ (tt) is less sensitive to the α s vari- In contrast, the determination of the charm-quark mass m c is not very sensitive to the value of α s because the variation of α s is compensated by a change in the gluon distribution at small x [12] (cf. also Table B in Ref. [4] for correlations between α S and the quark masses). Indeed, the fitted value of m c only changes by ±20 MeV for a variation of α s in a wide range, cf. Fig. 4. However, the NNLO corrections to the Wilson coefficents for DIS heavy-quark production still have significant impact on m c , moving it up by ∼ 100 MeV and reducing its uncertainty. The relatively weak correlation of m c with α s in Fig. 4 is in contrast to the observed behavior in the MMHT14 PDF set [13], where the particular variable flavor number scheme applied causes a linear relation between m c and α s in those fits and the charm-quark mass has been treated as a variable parameter, while the resulting values of a review of the theoretical treatments of DIS c-quark production used in PDF fits and Table 2 in Ref. [13] as well as Table 12 in Ref. [1] for the respective values of m c , α s and χ 2 in the MMHT14 PDF set. The grids for the NLO PDFs obtained in the present analysis are accessible with the LHAPDF library (version 6) [5] and available for download under http://projects. hepforge.org/lhapdf. For a fixed number of flavors, n f = 3, 4 and 5, we provide ABMP16_3_nlo (0+29), ABMP16_4_nlo (0+29), ABMP16_5_nlo (0+29), which consist of the central fit (set 0) and additional 29 sets for the combined symmetric uncertainties in all parameters (PDFs, α s , m c , m b and m t ). In each PDF set, the strong coupling α s is taken in the corresponding scheme, i.e., α and α (n f =5) s which can be related by the standard decoupling relations in QCD. As usual, the PDF set with three light-quarks n f = 3, ABMP16_3_nlo, is valid at all perturbative scales μ 2 1 GeV 2 , while those with n f = 4 and n f = 5, ABMP16_4_nlo and ABMP16_5_nlo, are subject to minimal cuts in μ 2 ≥ 3 GeV 2 and μ 2 ≥ 20 GeV 2 , respectively, cf. Ref. [4] for additional discussions.
For studies of LHC observables and their dependence on α s we also provide NLO PDF grids for n f =   in the on-shell scheme obtained using RunDec [23] are also given PDF set α where the value of α (n f =5) s (M Z ) has been fixed as indicated in the file names. These grids are determined by re-fitting all PDF parameters for the individual choices of α s , which for technical consistency remains a formal parameter in the fit, but with greatly suppressed uncertainty.
As the heavy-quark masses m c (m c ), m b (m b ) and m t (m t ) have been fitted their numerical values vary for each of the 29 PDF sets and cross section computations involving heavy quarks have to account for this. For reference we list in Table 4  The benchmark cross sections for the Higgs-boson and top-quark pair production at the LHC at NLO and NNLO with consistent use of the PDF sets obtained in the present analysis are given in Table 5. The quoted errors denote the PDF and α s uncertainties derived from the uncertainties in the experimental data. Thus, they are of similar size at NLO and NNLO.
The Higgs boson cross section σ (H ) is computed in the effective theory in the limit m t → ∞, but with full m t dependence in the Born cross section, based on the NNLO results of Refs. [15][16][17]. The NLO value of σ (H ) is about 20% smaller than the NNLO one due to missing large perturbative corrections, which are only partially compensated by a larger value of α s . In the effective theory the Born cross section for σ (H ) is proportional to α 2 s , so that the variant of the NLO fit with the larger value for the strong coupling, α NLO s (M Z ) = 0.1227, gives a NLO cross section increased by 5%, i.e. σ (H ) = 35.2 ± 0.58 pb with the PDF set ABMP16free_5_nlo compared to σ (H ) = 33.59 ± 0.58 pb with the PDF set ABMP16_5_nlo in Table 5.
The inclusive cross section σ (tt) for top-quark pair production uses Ref. [18] based on Refs. [19][20][21][22]. In this case, the NLO and NNLO values of σ (tt) for the range of centerof-mass energies explored at the LHC are similar, since those data have been included in both fits and are accommodated by the corresponding changes in the value of α s and the topquark mass m t , cf. Fig. 3 and Table 3.
In summary, we have completed the determination of the ABMP16 PDF sets at those orders of perturbation theory, which are currently of phenomenological relevance, i.e., at NLO and NNLO. Essential input in the ABMP16 analysis has been the final HERA DIS combination data from run I+II, which has consolidated the available world DIS data. In addition, several new data sets from the fixed-target DIS together with recent LHC and Tevatron data for the DY process and for the top-quark hadro-production have been used.
We have discussed the features of the NLO extraction of PDFs, which in general, have a few limitations due to lacking constraints of the higher order Wilson coefficients and we have emphasized the consistent use of PDFs and an order-dependent value of α s (M Z ), which is absolutely crucial because of correlations. The same holds, to a lesser Table 5 Cross sections at NLO and NNLO in QCD for the Higgs boson production from gluon-gluon fusion (σ  extent, in collider processes also for the values of the heavyquark masses used.
The ABMP16 PDFs establish the baseline for high precision analyses of LHC data from run I and run II, and the NLO variant is now available for computing cross sections of scattering processes with multi-particle final states, for which the NNLO QCD corrections will not be available in the foreseeable future, or for Monte Carlo studies. Precision analyses of LHC data, however, will always require analyses to NNLO accuracy in QCD. This will become even more important with the arrival of the data from the high luminosity runs.