We discuss how the two existing approximate N^3LO (aN^3LO) sets of parton distributions (PDFs) from the MSHT20 and NNPDF4.0 series can be combined for LHC phenomenology, both in the pure QCD case and for the QCD⊗QED sets that include the photon PDF. Using the resulting combinations, we present predictions for the total inclusive cross-section for Higgs production in gluon fusion, vector boson fusion, and associated production at the LHC Run-3. For the gluon fusion and vector boson fusion channels, the corrections that arise when using correctly matched aN^3LO PDFs with N^3LO cross section calculations, compared to using NNLO PDFs, are significant, in many cases larger than the PDF uncertainty, and generally larger than the differences between the two aN^3LO PDF sets entering the combination. The combined aN^3LO PDF sets, MSHT20xNNPDF40_an3lo and MSHT20xNNPDF40_an3lo_qed, are made publicly available in the LHAPDF format and can be readily used for LHC phenomenology.

We present a determination of the strong coupling $\alpha_s(m_Z)$ from a global dataset including both fixed-target and collider data from deep-inelastic scattering and a variety of hadronic processes, with a simultaneous determination of parton distribution functions (PDFs) based on the NNPDF4.0 methodology. This determination is performed at NNLO and approximate N$^3$LO (aN$^3$LO) perturbative QCD accuracy, including QED corrections and a photon PDF up to NLO accuracy. We extract $\alpha_s$ using two independent methodologies, both of which take into account the cross-correlation between $\alpha_s$ and the PDFs. The two methodologies are validated by closure tests that allow us to detect and remove or correct for several sources of bias, and lead to mutually consistent results. We account for all correlated experimental uncertainties, as well as correlated theoretical uncertainties related to missing higher order perturbative corrections (MHOUs). We study the perturbative convergence of our results and the impact of QED corrections. We assess individual sources of uncertainty, specifically MHOUs and the value of the top quark mass. We provide a detailed appraisal of methodological choices, including the choice of input dataset, the form of solution of evolution equation, the treatment of the experimental covariance matrix, and the details of Monte Carlo data generation. We find $\alpha_s(m_Z)=0.1194^{+0.0007}_{-0.0014}$ at aN$^3$LO$_{\rm QCD}\otimes {\rm NLO}_{\rm QED}$ accuracy, consistent with the latest PDG average and with recent lattice results.

We present a determination of the strong coupling α_s(m_Z) from a global dataset including both fixed-target and collider data from deep-inelastic scattering and a variety of hadronic processes, with a simultaneous determination of parton distribution functions (PDFs) based on the NNPDF4.0 methodology. This determination is performed at NNLO and approximate N^3LO (aN^3LO) perturbative QCD accuracy, including QED corrections and a photon PDF up to NLO accuracy. We extract α_s using two independent methodologies, both of which take into account the cross-correlation between α_s and the PDFs. The two methodologies are validated by closure tests that allow us to detect and remove or correct for several sources of bias, and lead to mutually consistent results. We account for all correlated experimental uncertainties, as well as correlated theoretical uncertainties related to missing higher order perturbative corrections (MHOUs). We study the perturbative convergence of our results and the impact of QED corrections. We assess individual sources of uncertainty, specifically MHOUs and the value of the top quark mass. We provide a detailed appraisal of methodological choices, including the choice of input dataset, the form of solution of evolution equation, the treatment of the experimental covariance matrix, and the details of Monte Carlo data generation. We find α_s(m_Z)=0.1194^+0.0007_-0.0014 at aN^3LO_ QCD⊗ NLO_ QED accuracy, consistent with the latest PDG average and with recent lattice results.

Abstract We extend the existing leading (LO), next-to-leading (NLO), and next-to-next-to-leading order (NNLO) NNPDF4.0 sets of parton distribution functions (PDFs) to approximate next-to-next-to-next-to-leading order ( $$\hbox {aN}^3\hbox {LO}$$ aN 3 LO ). We construct an approximation to the $$\hbox {N}^3\hbox {LO}$$ N 3 LO splitting functions that includes all available partial information from both fixed-order computations and from small and large x resummation, and estimate the uncertainty on this approximation by varying the set of basis functions used to construct the approximation. We include known $$\hbox {N}^3\hbox {LO}$$ N 3 LO corrections to deep-inelastic scattering structure functions and extend the FONLL general-mass scheme to $$\mathcal {O}\left( \alpha _s^3\right) $$ O α s 3 accuracy. We determine a set of $$\hbox {aN}^3\hbox {LO}$$ aN 3 LO PDFs by accounting both for the uncertainty on splitting functions due to the incomplete knowledge of $$\hbox {N}^3\hbox {LO}$$ N 3 LO terms, and to the uncertainty related to missing higher corrections (MHOU), estimated by scale variation, through a theory covariance matrix formalism. We assess the perturbative stability of the resulting PDFs, we study the impact of MHOUs on them, and we compare our results to the $$\hbox {aN}^3\hbox {LO}$$ aN 3 LO PDFs from the MSHT group. We examine the phenomenological impact of $$\hbox {aN}^3\hbox {LO}$$ aN 3 LO corrections on parton luminosities at the LHC, and give a first assessment of the impact of $$\hbox {aN}^3\hbox {LO}$$ aN 3 LO PDFs on the Higgs and Drell–Yan total production cross-sections. We find that the $$\hbox {aN}^3\hbox {LO}$$ aN 3 LO NNPDF4.0 PDFs are consistent within uncertainties with their NNLO counterparts, that they improve the description of the global dataset and the perturbative convergence of Higgs and Drell–Yan cross-sections, and that MHOUs on PDFs decrease substantially with the increase of perturbative order.

AbstractWe extend the existing leading (LO), next-to-leading (NLO), and next-to-next-to-leading order (NNLO) NNPDF4.0 sets of parton distribution functions (PDFs) to approximate next-to-next-to-next-to-leading order ($$\hbox {aN}^3\hbox {LO}$$aN3LO). We construct an approximation to the$$\hbox {N}^3\hbox {LO}$$N3LOsplitting functions that includes all available partial information from both fixed-order computations and from small and largexresummation, and estimate the uncertainty on this approximation by varying the set of basis functions used to construct the approximation. We include known$$\hbox {N}^3\hbox {LO}$$N3LOcorrections to deep-inelastic scattering structure functions and extend the FONLL general-mass scheme to$$\mathcal {O}\left( \alpha _s^3\right) $$Oαs3accuracy. We determine a set of$$\hbox {aN}^3\hbox {LO}$$aN3LOPDFs by accounting both for the uncertainty on splitting functions due to the incomplete knowledge of$$\hbox {N}^3\hbox {LO}$$N3LOterms, and to the uncertainty related to missing higher corrections (MHOU), estimated by scale variation, through a theory covariance matrix formalism. We assess the perturbative stability of the resulting PDFs, we study the impact of MHOUs on them, and we compare our results to the$$\hbox {aN}^3\hbox {LO}$$aN3LOPDFs from the MSHT group. We examine the phenomenological impact of$$\hbox {aN}^3\hbox {LO}$$aN3LOcorrections on parton luminosities at the LHC, and give a first assessment of the impact of$$\hbox {aN}^3\hbox {LO}$$aN3LOPDFs on the Higgs and Drell–Yan total production cross-sections. We find that the$$\hbox {aN}^3\hbox {LO}$$aN3LONNPDF4.0 PDFs are consistent within uncertainties with their NNLO counterparts, that they improve the description of the global dataset and the perturbative convergence of Higgs and Drell–Yan cross-sections, and that MHOUs on PDFs decrease substantially with the increase of perturbative order.

We provide a first quantitative indication that the wave function of the proton contains unequal distributions of charm quarks and antiquarks, i.e. a nonvanishing intrinsic valence charm distribution. A significant nonvanishing valence component cannot be perturbatively generated, hence our results reinforce previous evidence that the proton contains an intrinsic (i.e., not radiatively generated) charm quark component. We establish our result through a determination of the parton distribution functions (PDFs) of charm quarks and antiquarks in the proton. We propose two novel experimental probes of this intrinsic charm valence component: D-meson asymmetries in Z+c-jet production at the LHCb experiment, and flavor-tagged structure functions at the Electron-Ion Collider.

We extend the existing leading (LO), next-to-leading (NLO), and next-to-next-to-leading order (NNLO) NNPDF4.0 sets of parton distribution functions (PDFs) to approximate next-to-next-to-next-to-leading order ( aN^3LO ). We construct an approximation to the N^3LO splitting functions that includes all available partial information from both fixed-order computations and from small and large x resummation, and estimate the uncertainty on this approximation by varying the set of basis functions used to construct the approximation. We include known N^3LO corrections to deep-inelastic scattering structure functions and extend the FONLL general-mass scheme to 𝒪( α _s^3) accuracy. We determine a set of aN^3LO PDFs by accounting both for the uncertainty on splitting functions due to the incomplete knowledge of N^3LO terms, and to the uncertainty related to missing higher corrections (MHOU), estimated by scale variation, through a theory covariance matrix formalism. We assess the perturbative stability of the resulting PDFs, we study the impact of MHOUs on them, and we compare our results to the aN^3LO PDFs from the MSHT group. We examine the phenomenological impact of aN^3LO corrections on parton luminosities at the LHC, and give a first assessment of the impact of aN^3LO PDFs on the Higgs and Drell–Yan total production cross-sections. We find that the aN^3LO NNPDF4.0 PDFs are consistent within uncertainties with their NNLO counterparts, that they improve the description of the global dataset and the perturbative convergence of Higgs and Drell–Yan cross-sections, and that MHOUs on PDFs decrease substantially with the increase of perturbative order.

We include uncertainties due to missing higher order corrections to QCD computations (MHOU) used in the determination of parton distributions (PDFs) in the recent NNPDF4.0 set of PDFs. We use our previously published methodology, based on the treatment of MHOUs and their full correlations through a theory covariance matrix determined by scale variation, now fully incorporated in the new NNPDF theory pipeline. We assess the impact of the inclusion of MHOUs on the NNPDF4.0 central values and uncertainties, and specifically show that they lead to improved consistency of the PDF determination. PDF uncertainties on physical predictions in the data region are consequently either unchanged or moderately reduced by the inclusion of MHOUs.

We construct a set of parton distribution functions (PDFs), based on the recent NNPDF4.0 PDF set, that also include a photon PDF. The photon PDF is constructed using the LuxQED formalism, while QED evolution accounting for O(alpha), O(alpha alphas) and O(alpha^2) corrections is implemented and benchmarked by means of the EKO code. We investigate the impact of QED effects on NNPDF4.0, and compare our results both to our previous NNPDF3.1QED PDF set and to other recent PDF sets that include the photon. We assess the impact of photon-initiated processes and electroweak corrections on a variety of representative LHC processes, and find that they can reach the 5% level in vector boson pair production at large invariant mass.

We present a new set of parton distribution functions (PDFs) based on a fully global dataset and machine learning techniques: NNPDF4.0. We expand the NNPDF3.1 determination with 44 new datasets, mostly from the LHC. We derive a novel methodology through hyperparameter optimization, leading to an efficient fitting algorithm built upon stochastic gradient descent. We use NNLO QCD calculations and account for NLO electroweak corrections and nuclear uncertainties. Theoretical improvements in the PDF description include a systematic implementation of positivity constraints and integrability of sum rules. We validate our methodology by means of closure tests and “future tests” (i.e. tests of backward and forward data compatibility), and assess its stability, specifically upon changes of PDF parametrization basis. We study the internal compatibility of our dataset, and investigate the dependence of results both upon the choice of input dataset and of fitting methodology. We perform a first study of the phenomenological implications of NNPDF4.0 on representative LHC processes. The software framework used to produce NNPDF4.0 is made available as an open-source package together with documentation and examples.

An overwhelming number of theoretical predictions for hadron colliders require parton distribution functions (PDFs), which are an important ingredient of theory infrastructure for the next generation of high-energy experiments. This whitepaper summarizes the status and future prospects for determination of high-precision PDFs applicable in a wide range of energies and experiments, in particular in precision tests of the Standard Model and in new physics searches at the high-luminosity Large Hadron Collider and Electron-Ion Collider. We discuss the envisioned advancements in experimental measurements, QCD theory, global analysis methodology, and computing that are necessary to bring unpolarized PDFs in the nucleon to the N2LO and N3LO accuracy in the QCD coupling strength. Special attention is given to the new tasks that emerge in the era of the precision PDF analysis, such as those focusing on the robust control of systematic factors both in experimental measurements and theoretical computations. Various synergies between experimental and theoretical studies of the hadron structure are explored, including opportunities for studying PDFs for nuclear and meson targets, PDFs with electroweak contributions or dependence on the transverse momentum, for incisive comparisons between phenomenological models for the PDFs and computations on discrete lattice, and for cross-fertilization with machine learning/AI approaches. [Submitted to the US Community Study on the Future of Particle Physics (Snowmass 2021).]

We discuss the sensitivity of theoretical predictions of observables used in searches for new physics to parton distributions (PDFs) at large momentum fraction $x$. Specifically, we consider the neutral-current Drell-Yan production of gauge bosons with invariant masses in the TeV range, for which the forward-backward asymmetry of charged leptons from the decay of the gauge boson in its rest frame is a traditional probe of new physics. We show that the qualitative behaviour of the asymmetry depends strongly on the assumptions made in determining the underlying PDFs. We discuss and compare the large-$x$ behaviour of various different PDF sets, and find that they differ significantly. Consequently, the shape of the asymmetry observed at lower dilepton invariant masses, where all PDF sets are in reasonable agreement because of the presence of experimental constraints, is not necessarily reproduced at large masses where the PDFs are mostly unconstrained by data. It follows that the shape of the asymmetry at high masses may depend on assumptions made in the PDF parametrization, and thus deviations from the traditionally expected behaviour cannot be taken as a reliable indication of new physics. We demonstrate that forward-backward asymmetry measurements could help in constraining PDFs at large $x$ and discuss the accuracy that would be required to disentangle the effects of new physics from uncertainties in the PDFs in this region.

As data become more precise, estimating theoretical uncertainties in global PDF determinations is likely to become increasingly necessary to obtain correspondingly precise PDFs. Here we present a soon-to-be-released next generation of global proton PDFs (NNPDF4.0) that include theoretical uncertainties due to the use of heavy nuclear anddeuteron data in the fit. We estimate these uncertainties by comparing the values of thenuclear observables computed with the nuclear PDFs against thosecomputed with proton PDFs. For heavy nuclear PDFs we use thenuclear nNNPDF2.0 set, while for deuteron PDFs we develop an iterativeprocedure to determine proton and deuteron PDFs simultaneously, eachincluding the uncertainties in the other. Accounting for nuclearuncertainties resolves some of the tensions in the global fit of the proton PDFs, especially those between the nuclear data and the extended LHC data set used in NNPDF4.0.

We respond to the criticism raised by Courtoy et al., in which the faithfulness of the NNPDF4.0 sampling is questioned and an under-estimate of the NNPDF4.0 PDF uncertainties is implied. We list, correct, and clarify in detail a number of inaccurate or misleading claims that are made in this Reference. Specifically, we explain and explicitly demonstrate why the central value of the PDF distribution does not generally coincide with the absolute minimum of the χ^2 to the data. We examine some PDFs that have been constructed in the above study and claimed to be "good solutions": we show that similar PDFs are found with the NNPDF methodology, but with very low probability.

A precise knowledge of the quark and gluon structure of the proton, encoded by the parton distribution functions (PDFs), is of paramount importance for the interpretation of high-energy processes at present and future lepton-hadron and hadron-hadron colliders. Motivated by recent progress in the PDF determinations carried out by the CT, MSHT, and NNPDF groups, we present an updated combination of global PDF fits: PDF4LHC21. It is based on the Monte Carlo combination of the CT18, MSHT20, and NNPDF3.1 sets followed by either its Hessian reduction or its replica compression. Extensive benchmark studies are carried out in order to disentangle the origin of the differences between the three global PDF sets. In particular, dedicated fits based on almost identical theory settings and input datasets are performed by the three groups, highlighting the role played by the respective fitting methodologies. We compare the new PDF4LHC21 combination with its predecessor, PDF4LHC15, demonstrating their good overall consistency and a modest reduction of PDF uncertainties for key LHC processes such as electroweak gauge boson production and Higgs boson production in gluon fusion. We study the phenomenological implications of PDF4LHC21 for a representative selection of inclusive, fiducial, and differential cross sections at the LHC. The PDF4LHC21 combination is made available via the LHAPDF library and provides a robust, user-friendly, and efficient method to estimate the PDF uncertainties associated to theoretical calculations for the upcoming Run III of the LHC and beyond.

The theory of the strong force, quantum chromodynamics, describes the proton in terms of quarks and gluons. The proton is a state of two up quarks and one down quark bound by gluons, but quantum theory predicts that in addition there is an infinite number of quark-antiquark pairs. Both light and heavy quarks, whose mass is respectively smaller or bigger than the mass of the proton, are revealed inside the proton in high-energy collisions. However, it is unclear whether heavy quarks also exist as a part of the proton wavefunction, which is determined by non-perturbative dynamics and accordingly unknown: so-called intrinsic heavy quarks. It has been argued for a long time that the proton could have a sizable intrinsic component of the lightest heavy quark, the charm quark. Innumerable efforts to establish intrinsic charm in the proton have remained inconclusive. Here we provide evidence for intrinsic charm by exploiting a high-precision determination of the quark-gluon content of the nucleon based on machine learning and a large experimental dataset. We disentangle the intrinsic charm component from charm-anticharm pairs arising from high-energy radiation. We establish the existence of intrinsic charm at the 3-standard-deviation level, with a momentum distribution in remarkable agreement with model predictions. We confirm these findings by comparing to very recent data on Z-boson production with charm jets from the LHCb experiment.

As data become more precise, estimating theoretical uncertainties in global PDF determinations is likely to become increasingly necessary to obtain correspondingly precise PDFs.Here we present a next generation of global proton PDFs (NNPDF4.0) that include theoretical uncertainties due to the use of heavy nuclear and deuteron data in the fit.We estimate these uncertainties by comparing the values of the nuclear observables computed with the nuclear PDFs against those computed with proton PDFs.For heavy nuclear PDFs we use the nuclear nNNPDF2.0set, while for deuteron PDFs we develop an iterative procedure to determine proton and deuteron PDFs simultaneously, each including the uncertainties in the other.Accounting for nuclear uncertainties resolves some of the tensions in the global fit of the proton PDFs, especially those between the nuclear data and the extended LHC data set used in NNPDF4.0.

We evaluate the uncertainties due to nuclear effects in global fits of proton parton distribution functions (PDFs) that utilise deep-inelastic scattering and Drell–Yan data on deuterium targets. To do this we use an iterative procedure to determine proton and deuteron PDFs simultaneously, each including the uncertainties in the other. We apply this procedure to determine the nuclear uncertainties in the SLAC, BCDMS, NMC and DYE866/NuSea fixed target deuteron data included in the NNPDF3.1 global fit. We show that the effect of the nuclear uncertainty on the proton PDFs is small, and that the increase in overall uncertainties is insignificant once we correct for nuclear effects.

We study the correlation between different sets of parton distributions (PDFs). Specifically, viewing different PDF sets as distinct determinations, generally correlated, of the same underlying physical quantity, we examine the extent to which the correlation between them is due to the underlying data. We do this both for pairs of PDF sets determined using a given fixed methodology, and between sets determined using different methodologies. We show that correlations have a sizable component that is not due to the underlying data, because the data do not determine the PDFs uniquely. We show that the data-driven correlations can be used to assess the efficiency of methodologies used for PDF determination. We also show that the use of data-driven correlations for the combination of different PDFs into a joint set can lead to inconsistent results, and thus that the statistical combination used in constructing the widely used PDF4LHC15 PDF set remains the most reliable method.

We present the software framework underlying the NNPDF4.0 global determination of parton distribution functions (PDFs). The code is released under an open source licence and is accompanied by extensive documentation and examples. The code base is composed by a PDF fitting package, tools to handle experimental data and to efficiently compare it to theoretical predictions, and a versatile analysis framework. In addition to ensuring the reproducibility of the NNPDF4.0 (and subsequent) determination, the public release of the NNPDF fitting framework enables a number of phenomenological applications and the production of PDF fits under user-defined data and theory assumptions.