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.

Italy

Andrea Barontini

CERN

Switzerland

Alessandro Candido

Finland

Felix Hekhorn

Netherlands

Giacomo Magni

United Kingdom

Roy Stegeman

An FONLL prescription with coexisting flavor number PDFs

We discuss a Bayesian methodology for the solution of the inverse problem underlying the determination of parton distribution functions (PDFs). In our approach, Gaussian processes (GPs) are used to model the PDF prior, while Bayes’ theorem is used in order to determine the posterior distribution of the PDFs given a set of data. We discuss the general formalism, the Bayesian inference at the level of both parameters and hyperparameters, and the simplifications which occur when the observable entering the analysis is linear in the PDF. We benchmark the new methodology in two simple examples for the determination of a single PDF flavor from a set of deep inelastic scattering (DIS) data and from a set of equal-time correlators computed using lattice QCD. We discuss our results, showing how the proposed methodology allows for a well-defined statistical interpretation of the different sources of errors entering the PDF uncertainty, and how results can be validated a posteriori.

Luigi Debbio

Tommaso Giani

Giacomo Petrillo

Bayesian inference with Gaussian processes for the determination of parton distribution functions

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.

Tanjona Rabemananjara

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.

Richard Ball

Stefano Carrazza

Juan Cruz-Martinez

Luigi Del Debbio

Stefano Forte

Zahari Kassabov

Niccolò Laurenti

Emanuele Nocera

Juan Rojo

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

The interpretation of present and future neutrino experiments requires accurate theoretical predictions for neutrino-nucleus scattering rates. Neutrino structure functions can be reliably evaluated in the deep-inelastic scattering regime within the perturbative QCD (pQCD) framework. At low momentum transfers (Q 2 ≲ few GeV2), inelastic structure functions are however affected by large uncertainties which distort event rate predictions for neutrino energies E  ν  up to the TeV scale. Here we present a determination of neutrino inelastic structure functions valid for the complete range of energies relevant for phenomenology, from the GeV region entering oscillation analyses to the multi-EeV region accessible at neutrino telescopes. Our NNSFν approach combines a machine-learning parametrisation of experimental data with pQCD calculations based on state-of-the-art analyses of proton and nuclear parton distributions (PDFs). We compare our determination to other calculations, in particular to the popular Bodek-Yang model. We provide updated predictions for inclusive cross sections for a range of energies and target nuclei, including those relevant for LHC far-forward neutrino experiments such as FASERν, SND@LHC, and the Forward Physics Facility. The NNSFν determination is made available as fast interpolation LHAPDF grids, and it can be accessed both through an independent driver code and directly interfaced to neutrino event generators such as GENIE.

United States

Spain

Alfonso Garcia

Neutrino structure functions from GeV to EeV energies

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.

Parton distributions and new physics searches: the Drell–Yan forward–backward asymmetry as a case study