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 $${\textrm{N}}^3$$ LO ( $${\textrm{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 $${\textrm{aN}}^3$$ $${\textrm{LO}}_{\textrm{QCD}}\otimes {\textrm{NLO}}_{\textrm{QED}}$$ accuracy, consistent with the latest PDG average and with recent lattice results.

United Kingdom

Richard Ball

Italy

Andrea Barontini

Switzerland

Juan Cruz-Martinez

Stefano Forte

Finland

Felix Hekhorn

Emanuele Nocera

Netherlands

Juan Rojo

Roy Stegeman

A determination of $$\alpha _s(m_Z)$$ at $${{\textrm{aN}}}^3{{\textrm{LO}}}_{{\textrm{QCD}}}\otimes {{\textrm{NLO}}}_{{\textrm{QED}}}$$ accuracy from a global PDF analysis

We present NNPDF 2.0, a new set of collinear helicity parton distribution functions (PDFs) of the proton based on legacy measurements of structure functions in inclusive neutral-current longitudinally polarised deep-inelastic scattering (DIS), and of $\textit{W}$ -boson, single-inclusive, and di-jet production asymmetries in longitudinally polarised proton-proton collisions. The determination is accurate to next-to-next-to-leading order in the strong coupling, and includes heavy quark mass corrections in the analysis of DIS data. Uncertainties due to missing higher-order corrections are systematically incorporated by means of a covariance matrix determined by scale variations. NNPDF 2.0 is based on a machine learning methodology, that makes use of Monte Carlo sampling for the representation of uncertainties into PDFs, of a neural network for the parametrisation of PDFs, of stochastic gradient descent for the optimisation of PDF parameters, and of hyperoptimisation for the selection of the best fitting model. We study the impact on PDFs of higher-order corrections, of the positivity constraint, and of the data. We demonstrate two phenomenological applications of NNPDF 2.0, specifically the determination of the proton spin fraction carried by gluons and quarks, and of theoretical predictions for single-hadron production in longitudinally polarised DIS and proton-proton collisions.

Toon Hasenack

Giacomo Magni

Tanjona Rabemananjara

Tanishq Sharma

Gijs Seeventer

NNPDF

We present a new prescription to account for heavy quark mass effects in the determination of parton distribution functions (PDFs) based on the FONLL scheme. Our prescription makes explicit use of the freedom to choose the number of active flavors at a given scale and, thus, use coexisting PDFs with different active flavor number. This new prescription is perturbatively equivalent to the former but improves the implementation in two ways. First, it can be naturally generalized to account simultaneously for multiple heavy quark effects, such as charm and bottom effects, which can both be relevant at the same scale due to the small mass difference. Second, it can be trivially generalized to use at any fixed-order or collinear resummed accuracy, while previous prescriptions required ad-hoc expansions of the DGLAP evolution kernels for each coefficient. We supplement the paper with codes for the computation of deep inelastic scattering observables in this new prescription.

CERN

Alessandro Candido

An FONLL prescription with coexisting flavor number PDFs

We present yadism, a library for the evaluation of both polarized and unpolarized deep-inelastic scattering (DIS) structure functions and cross sections up to N $$^{3}$$  <math> <mmultiscripts> <mrow></mrow> <mrow></mrow> <mn>3</mn> </mmultiscripts> </math>  LO in perturbative QCD. The package provides computations of observables in fixed-flavor and zero-mass variable flavor number schemes. The implementation of the general mass variable flavor number schemes is supported through the high virtuality limits for the heavy flavor coefficients. In addition, yadism provides a set of tools for the generation of interpolation grids in the PDF-independent PineAPPL format, allowing to test the PDF dependence on any DIS observable without needing to rerun the computation. This work is part of an ongoing effort to standardize the format of theory predictions in high-energy physics within the pineline framework. The code is open source, written in Python and documented to facilitate usage, integrations, and further extensions. Finally, the code has been benchmarked against the widely used APFEL++ and QCDNUM libraries.

Yadism: yet another deep-inelastic scattering module

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}$$  <math> <mrow> <msup> <mtext>aN</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>  ). We construct an approximation to the  $$\hbox {N}^3\hbox {LO}$$  <math> <mrow> <msup> <mtext>N</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   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}$$  <math> <mrow> <msup> <mtext>N</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   corrections to deep-inelastic scattering structure functions and extend the FONLL general-mass scheme to  $$\mathcal {O}\left( \alpha _s^3\right) $$  <math> <mrow> <mi>O</mi> <mfenced> <msubsup> <mi>α</mi> <mi>s</mi> <mn>3</mn> </msubsup> </mfenced> </mrow> </math>   accuracy. We determine a set of  $$\hbox {aN}^3\hbox {LO}$$  <math> <mrow> <msup> <mtext>aN</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   PDFs by accounting both for the uncertainty on splitting functions due to the incomplete knowledge of  $$\hbox {N}^3\hbox {LO}$$  <math> <mrow> <msup> <mtext>N</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   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}$$  <math> <mrow> <msup> <mtext>aN</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   PDFs from the MSHT group. We examine the phenomenological impact of  $$\hbox {aN}^3\hbox {LO}$$  <math> <mrow> <msup> <mtext>aN</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   corrections on parton luminosities at the LHC, and give a first assessment of the impact of  $$\hbox {aN}^3\hbox {LO}$$  <math> <mrow> <msup> <mtext>aN</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   PDFs on the Higgs and Drell–Yan total production cross-sections. We find that the  $$\hbox {aN}^3\hbox {LO}$$  <math> <mrow> <msup> <mtext>aN</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   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.

Stefano Carrazza

Luigi Del Debbio

Tommaso Giani

Zahari Kassabov

Niccolò Laurenti

Germany

Christopher Schwan

Maria Ubiali

The path to  $$\hbox {N}^3\hbox {LO}$$  <math> <mrow> <msup> <mtext>N</mtext> <mn>3</mn> </msup> <mtext>LO</mtext> </mrow> </math>   parton distributions

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  $${\mathcal {O}}\left( \alpha \right) $$  <math> <mrow> <mi>O</mi> <mfenced> <mi>α</mi> </mfenced> </mrow> </math>  ,  $${\mathcal {O}}\left( \alpha \alpha _s\right) $$  <math> <mrow> <mi>O</mi> <mfenced> <mi>α</mi> <msub> <mi>α</mi> <mi>s</mi> </msub> </mfenced> </mrow> </math>  , and  $${\mathcal {O}}\left( \alpha ^2\right) $$  <math> <mrow> <mi>O</mi> <mfenced> <msup> <mi>α</mi> <mn>2</mn> </msup> </mfenced> </mrow> </math>   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.

Photons in the proton: implications for the LHC

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.

Determination of the theory uncertainties from missing higher orders on NNLO parton distributions with percent accuracy

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 of charm quarks and antiquarks in the proton. We propose two novel experimental probes of this intrinsic charm valence component: <math><mi>D</mi></math>-meson asymmetries in <math><mrow><mi>Z</mi><mo>+</mo><mi>c</mi></mrow></math>-jet production at the LHCb experiment and flavor-tagged structure functions at the Electron-Ion Collider.

Richard D. Ball

Emanuele R. Nocera

Intrinsic charm quark valence distribution of the proton

We revisit our argument that shows that parton distribution functions (PDFs) in the  $$\overline{\textrm{MS}}$$  <math> <mover> <mtext>MS</mtext> <mo>¯</mo> </mover> </math>   scheme are non-negative in the perturbative region, with the dual goals of clarifying its theoretical underpinnings and elucidating its domain of validity. We specifically discuss recent results proving that PDFs can turn negative at sufficiently low scale, we clarify quantitatively various aspects of our derivation of positivity in the perturbative region, and we provide an estimate for the scale above which PDF positivity holds.

On the positivity of  $$\overline{\textrm{MS}}$$  <math> <mover> <mtext>MS</mtext> <mo>¯</mo> </mover> </math>   parton distributions

We extend the FONLL general-mass variable-flavour-number scheme to the case of longitudinally polarised DIS structure functions, accounting for perturbative corrections up to  $$\mathcal {O}\left( \alpha _s^2\right) $$  <math> <mrow> <mi>O</mi> <mfenced> <msubsup> <mi>α</mi> <mi>s</mi> <mn>2</mn> </msubsup> </mfenced> </mrow> </math>  . We quantify the impact of charm quark mass and higher-order perturbative corrections on projected measurements of inclusive and charm-tagged longitudinal asymmetries at the Electron-Ion Collider (EIC) and at the Electron-ion collider in China (EicC). We demonstrate how the inclusion of these corrections is essential to compute predictions with an accuracy that matches the projected precision of the measurements. The computation is made publicly available through the open-source EKO and YADISM programs.

Adrianne Schaus

Heavy quarks in polarised deep-inelastic scattering at the electron-ion collider