Impact of SeaQuest data on PDF fits at large x

We evaluate the impact of recent SeaQuest (FNAL-E906 experiment) data on dimuon production in proton-deuteron and proton-proton collisions on parton distribution functions (PDFs). We find these data in a good agreement with the QCD predictions based on PDFs fitted to the Tevatron and LHC data on forward production of W and Z bosons. As a basis for this study we use the ABMP16 PDF fits and show that they turn out to be compatible with the SeaQuest data, and that these data have constraining power, allowing to reduce the uncertainties on the isospin asymmetry of the light-sea-quark distribution at large longitudinal momentum fraction x. We discuss the nuclear corrections needed to describe the deuteron and show that they affect the theoretical description of the proton-deuteron Drell–Yan cross section at the level of O(0.5-1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}(0.5{-}1)$$\end{document}%. We also comment on the compatibility of the SeaQuest results with other state-of-the-art PDF fits and show that these data are in clear disagreement with models proposing an SU(3)-flavor symmetric quark sea. Finally, we perform a comparison between the second Mellin moments of the light-quark PDFs and recent results from various lattice QCD computations, which demonstrates good compatibility, albeit limited by the uncertainties inherent in current lattice QCD simulations.


I. INTRODUCTION
The knowledge of parton distribution functions (PDFs) at large longitudinal momentum fraction x is one of the most urgent open questions [1,2] concerning proton and nuclear structure to which not only theoretical but even experimental efforts are going to be dedicated in the future.While in the long term data from the Electron Ion Collider (EIC) [3,4] are expected to play a very important constraining role, as also emphasized in the Snowmass 2021 EIC-dedicated whitepaper [5], in the near future further experiments might also offer promising opportunities.Among those at the Large Hadron Collider (LHC), we mention here the fixed-target (FT) configurations exploiting one of the LHC beams [6], a possibility already realized by complementing the LHCb detector with the SMOG and SMOG2 apparata [7,8], and also conceptually studied, although not realized, by the ALICE collaboration with the ALICE-FT experiment [9], as well as perspective projects still under discussion, like the Forward Physics Facility [10,11].In particular, the LHCb + SMOG system has already delivered the first data using proton and Pb beams impinging on gaseous nuclei like 4 He, 20 Ne and 40 Ar, at different nucleon-nucleon center-of-mass energies √ s NN ∼ O(50 − 100) GeV, corresponding to various past LHC runs.The LHCb + SMOG2 system, active during Run 3 and 4, can make use of even lighter gases, like deuterium 2 H as well as hydrogen, with increased statistics 1 .These experiments allow to probe the longitudinal momentum fraction interval 0.1 < x < 1 for target partons, on which the constraints from the sets of HERA data [12] which traditionally form the backbone of PDF fits, are quite loose and mostly indirect 2 .
For the time being, constraints on PDFs at large x are imposed by legacy measurements from inclusive deep-inelastic scattering (DIS) experiments at fixed-targets (SLAC, BCDMS, NMC, etc.), semi-inclusive DIS experiments using ν beams and capable of measuring heavy-quark production in DIS (CCFR, NuTeV, CHORUS, NOMAD, etc.) and fixed-target Drell-Yan (DY) experiments (CERN-NA51, FNAL-E605, FNAL-E866, etc.), complemented by measurements of crosssections for DY (+ jets) production and other specific processes in the main detectors at the Tevatron and the LHC in the standard collider-mode configuration (for an overview, see e.g.Ref. [1] and references therein).
The valence quark distributions are constrained by DIS HERA data, up to x < 0.1, and in fixedtarget experiments, up to x ∼ 1.A large-x and relatively low-Q domain is also probed at JLab [14].The DY data from the Tevatron and the LHC (both inclusive cross-sections and charge asymmetries) as well as from fixed-target experiments have also been used to probe up (u) and down (d) quark distributions and their differences (isospin asymmetries).Single-top quark production data have allowed to probe the u/d ratio at x ∼ 0.1, where u = u val + u sea and d = d val + d sea , notwithstanding the big systematic uncertainties still accompanying the experimental cross-sections for this channel of top-quark production [15].The (anti)strange sea quark distributions (s, s) has been constrained by DY (+ jets) LHC data and older (anti)-neutrino-nuclear DIS data, with large uncertainties [1,16,17], and improving their determination remains one of the pressing issues in PDF analysis.The s(x) − s(x) asymmetry [18] can be constrained by semi-inclusive DIS data on dimuon production distinguishing neutrino and antineutrino beams (as discussed e.g. in Ref. [19] and [20]), by W + + c and W − + c data at the LHC [21] and by future DIS experiments using separate beams of neutrinos and anti-neutrinos (e.g. at the Forward Physics Facility).The up and down sea quark distributions are well constrained by DY data.Finally, the gluon PDF at large-x is mostly constrained by measurements of heavy-quark and jet production at the LHC [22].
Recently, the SeaQuest collaboration (FNAL-E906 experiment) has released fixed-target data on dimuon production on 2 H and proton targets through DY, which allow to constrain the difference between down and up sea quarks, i.e. d(x) − ū(x), and the d(x)/ū(x) ratio [23].This experiment can be considered as a continuation of previous ones, FNAL-E866 [24] (NuSea) and FNAL-E605 [25], lowering the center-of-mass energy √ s and extending the kinematic coverage in x3 .Reviews on the flavour structure of the nucleon sea, triggering further investigations, have been provided by e.g.Ref. [28,29].The new experimental results and the present theoretical scenario motivate the present study, where we focus on the light-quark distributions, with particular emphasis on the sea quark case.In Sec.II, we show the impact of the SeaQuest results on the ABMP16 fits, considering both versions, at next-to-leading order (NLO) and at next-to-NLO (NNLO) in perturbative QCD, published in Refs.[30] and [31], respectively.This leads to new PDF fits, dubbed as ABMP16 + SeaQuest NLO and NNLO, performed ab-initio using the same statistical methodology and inputs as for the ABMP16 fits, plus the most recent SeaQuest data.In Sec.III we comment on the compatibility of other state-of-the-art PDF fits with these data and in Sec.IV we discuss nuclear corrections.In Sec.V we compare the second moments of the light-flavor quark PDFs with recent lattice QCD results.Our conclusions are delivered in Sec.VI.

II. CONSTRAINING POWER OF THE SEAQUEST DATA ON THE ABMP16 NLO AND NNLO PDF FITS
The study extends the ABMP16 PDF fits (NLO and NNLO), which have used the combined data from HERA for inclusive DIS, data from the fixed-target experiments NOMAD and CHORUS for neutrino-induced DIS, as well as data from Tevatron and the LHC for the DY process and the hadro-production of single-top and top-quark pairs.The ABMP16 approach uses a fixed-flavor number scheme for n f = 3, 4, 5 and simultaneously determines the PDFs, the value of the strong coupling α s (M z ) and all masses of heavy quarks, fully preserving the correlations among these quantities.
For illustrative purpose, we summarize in Fig. 1 the (x 1 , x 2 ) coverage of most of the DY data used in constraining the up and down sea quark distributions at large x in these fits 4 , together with the (x 1 , x 2 ) coverage of the recently released SeaQuest data.The variables x 1 and x 2 represent the momentum fractions carried by the incident (anti)quarks in beam and target, respectively, which roughly characterize the region of x probed by a particular experiment.Since x 1 and x 2 are not observables and cannot be measured, we detail here how we reconstruct them, assuming leading order (LO) kinematics.For the SeaQuest experiment x 1,2 are computed as follows: where Q is the four-momentum of the virtual photon from the quark-antiquark annihilation in the non-resonant production process, P 1,2 are the four-momenta of the projectile and target hadron, respectively, and P = P 1 + P 2 .Considering γ * → µ + µ − decays, the average values for x 1,2 in the bins of the muon-pair average Feynman variable ⟨x F ⟩ are reported in Ref. [23].These values are plotted in Fig. 1 in comparison with the kinematics of other DY data included in the ABMP16 fits.In particular, for the E605 Fermilab fixed-target data [25], given in the form of a double differential distribution in √ τ = M/ √ s and y, where M and y are the invariant mass and rapidity of the µ + µ − -pair, respectively, and √ s is the collision center-of-mass energy, the values of x 1,2 are computed according to the relation For the E866 experiment [24] the same relation is employed.However, since the muon-pair rapidity is not tabulated in Ref. [24], it is computed from the muon-pair x F and transverse momentum p T using basic definition as follows: where p L = x F p L, max and E = p 2 L + p 2 T + M 2 are the muon-pair longitudinal momentum and energy, respectively, in the center-of-mass frame of the colliding hadrons, with p L,max the maximum longitudinal momentum of the muon-pair, depending on √ s according to the formula p L,max = √ s (1 − M 2 /s)/2.The approach of Eq. ( 2) is also used for the LHCb data on Z-boson production [34][35][36] released in the form of lepton-pair pseudorapidity distributions 5 .The data on W-boson production evidently probe the same kinematics.However, the use of Eq. ( 2) is impossible in this case due to the neutrino escaping detection.Therefore, for W-boson production in the D0 experiment [37,38], we use the following approximate estimate: where M W is the W-boson pole mass and y l is the lepton rapidity.
Both DY data at the Tevatron and the LHC and DY data in fixed-target experiments play a role in constraining the sea quark PDFs at large x.They allow to reach similarly large-x values, although in the case of fixed-target data, relatively large-x partons from both the projectile and the target participate in the same hard interaction, whereas in the case of the LHC, a large-x parton is typically probed simultaneously with a low-x one, as exemplified in the (x 1 , x 2 ) correlation in Fig. 1.The correlation is quite evident for the LHC data and is related to the exchange of heavy bosons in the DY process.In the case of LHC, the largest x 1 values are probed by the LHCb detector with data at large positive rapidity, which covers the interval 2 < y < 4.5.On the other hand, the fixed-target experiments E866 and E605, which have much lower center-ofmass energies than the LHC, probe larger x 2 values and present a less evident (x 1 , x 2 ) correlation, related to the exchange of a γ * with a broad range of mass values in the DY process.In the case of SeaQuest, the (x 1 , x 2 ) correlation is again evident, considering that the invariant mass of the observed γ * decay products is fixed to approximately M ∼ 5 GeV.SeaQuest covers x 2 values higher than the LHC due to the use of a beam with much lower center-of-mass energy ( √ s = 15.1 GeV).FIG.1: The (x 1 , x 2 ) coverage for the SeaQuest experiment [23] (full circles), with x 1,2 given by Eq. ( 1), in comparison to the coverages from DY data of other experiments used in the ABMP16 PDF fits (downoriented triangles: E866 Fermilab fixed-target experiment [24]; up-oriented triangles: E605 [25] with x 1,2 computed from the lepton-pair rapidity using Eq. ( 2) for both the data sets; squares: LHCb, the LHC experiment [34][35][36] with x 1,2 computed from the Z-boson rapidity using Eq. ( 2); open circles: D0, the Tevatron collider experiment [37,38], with x 1,2 estimated from the charged-lepton rapidity using Eq. ( 4)).
The x 2 region covered by SeaQuest extends up to ≤ 0.45.The E605 experiment has a coverage extending even up to slightly higher x 2 values.However, the E605 experiment used a Copper target, thus requiring an evaluation of nuclear corrections (see the end of Section IV).Having only one target material, they could not provide data on cross-section ratios, unlike SeaQuest that has both a deuteron and a proton target.Also, given that Copper is a heavy nucleus close to isoscalarity, the E605 data are much less sensitive to the isospin asymmetry effects, that we investigate in this work.We explicitly verified the very small impact of E605 data, by removing them from our fits where they are included as default.
The green dots along two parallel lines in Fig. 1 refer to the cases of Z-boson production at the LHC at √ s = 7 and 8 TeV, given that data at these center-of-mass energies were included in the ABMP16 PDF fits.The SeaQuest data [23] on the ratio of pd and pp DY distributions over Feynman variable x F with respect to the NNLO predictions obtained using the code VRAP [39] (solid line) and DYNNLO [40] (dashes).

SeaQuest (
In order to compute predictions for the DY cross-sections, we use the FEWZ2.1 code [41] for the collider cases and the VRAP code [39] for the fixed-target cases.In particular, the present analysis of SeaQuest data is based on the x F -distribution, that was directly measured in the experiment and could also be computed to NNLO QCD accuracy using a Monte-Carlo code, like e.g.FEWZ or DYNNLO [40] 6 .However, in the fit we employ the VRAP code, which is based on 2dimensional integration that allows to greatly improve the code performance.To compute VRAP predictions for the SeaQuest data on the x F -distribution we perform a mapping of x F to the rapidity using the basic relation Eq. ( 2) and taking , where s is the center-of-mass energy squared and ⟨x F ⟩, ⟨M⟩ and ⟨P T ⟩ are the averages of muonpair Feynman variable x F , invariant mass and transverse momentum, respectively, over the bins in x F .These averaged quantities are all given in Ref. [23] for each bin in x F .To validate such an approach, we compare its predictions with those obtained with the methodology used in Ref. [23], where the DYNNLO [40] code is employed, instead of VRAP, and the exact information concerning the transverse momentum and the invariant mass of the µ + µ − -pair on an event-by-event basis is considered to build the x F distributions, instead of the average value of these quantities per x F bin.We find that the difference is mostly well below the data uncertainties, cf.Fig. 2, where we compare predictions obtained with VRAP using the approximations outlined above with the DYNNLO predictions based on the exact values for P T and M as input, as in Ref. [23], and applying their same cut M > 4.5 GeV.The latter suppresses the µ + µ − -background contribution from J/ψ and ψ ′ production and decay.We build the x F distributions using Eq. ( 4) of Ref. [23].From Fig. 2 it is evident that only in the smallest x F bin the difference between VRAP and DYNNLO is comparable to the data uncertainty, an observation that might be related to the width of the bin, which is much larger for this bin, than for the other ones.Obviously, such a difference cannot have relevant impact on fit results.Therefore, considering that the approximated procedure with the use of VRAP allows for NNLO simulations much faster than the exact procedure using DYNNLO, and given that the results turn out to be very well compatible, we use VRAP for the analyses and all other plots presented in the rest of this work.
We observe that the corrections related to spectrometer acceptance as a function of x 1 and x 2 reported in Ref. [43] do not impact distributions depending on measured quantities, like e.g., x F , that we consider in this work.On the other hand, their inclusion is relevant for the extraction of the d − ū asymmetry and the d/ū ratio from the SeaQuest data, as performed by the SeaQuest collaboration in approximated form as described in their papers [23,43].
The constraints from the SeaQuest experiment turn out to be compatible with those already imposed by collider data, as shown by the fact that the χ 2 /NDP of the analyses including also the SeaQuest data does not change significantly with respect to the χ 2 /NDP of the original ABMP16 analyses.Here NDP indicates the number of data points and the differences are well within the χ 2 statistical uncertainties, as shown in Tab.I.The χ 2 /NDP for the NNLO analyses turns out to be 1.18, slightly closer to 1 than the χ 2 /NDP of the NLO analyses, which is equal to 1.20.We also observe that incorporating SeaQuest data in the fits has a negligible impact on the values of α s (M Z ) and heavy-quark masses, extracted simultaneously to PDFs in all the fits considered in Tab.I.
Separate χ 2 values for various data sets included in our NNLO QCD analyses are reported in Tab.II.We have considered four variants: (I) the ABMP16 analysis, (II) the ABMP16 + SeaQuest analysis, as well as (III) an analysis, where we consider all data of (II), except the D0 DY data and (IV) an analysis, where we consider all data of (II), except the LHCb DY data.We include variants (III) and (IV) due to the fact that in the past we have observed some tension between the D0 and LHCb DY data.By comparing (I) and (II), we find that, for each considered data set, the addition of SeaQuest data does not introduce significant modifications of the χ 2 .Thus SeaQuest data are well compatible with both the LHCb and the D0 DY data.On the other hand, by comparing (II) and (III), we find that the elimination of the D0 DY datasets from the fit allows to improve the χ 2 of the analysis of the 7 TeV LHCb DY dataset by several units, beyond the statistical χ 2 uncertainty.Vice versa, the elimination of the LHCb DY datasets allows to improve the description of D0 data, as can be understood by comparing (II) and (IV).
The χ 2 values were computed accounting for statistical and systematic uncertainties of the SeaQuest data, assuming that systematic uncertainties are fully correlated bin-by-bin.Detailed information concerning correlations among the uncertainties characterizing the SeaQuest data are, however, not available.Therefore, we also consider a variant of the fit, where the systematic uncertainties are considered as fully uncorrelated.We have found that the χ 2 values related to   the analysis of the SeaQuest data in both analyses are compatible within statistical fluctuations (χ 2 corr.= 7.3 vs. χ 2 uncorr.= 5.9, for NDP = 7).This implies that more details on the precise degree of bin-by-bin correlations of the systematic uncertainties in the SeaQuest data, when available, will not modify the main conclusions of our study.Fig. 3 shows the constraining power of the SeaQuest data on the d(x) − ū(x) difference, increasing towards large x values.At NLO, the uncertainty band of the analysis with SeaQuest data has a large overlap, but is not completely included within the band of the default ABMP16 analysis (not including these data).Additionally, for 0.1 < x < 0.2, the band of the analysis with SeaQuest data turns out to be of the same size of the band of the one without these data.On the other hand, at NNLO, the uncertainty bands are in general smaller than at NLO and the one of the analysis with SeaQuest data is always included and smaller than the band of the analysis without SeaQuest data.These findings confirm that theory predictions at NNLO accuracy are in general more robust and consistent among each other than NLO ones, i.e., the theory description at NLO is still incomplete and hardly provides a simultaneous excellent description of all DY data, like is instead happening at NNLO.This is also reflected in the comparison of the χ 2 /NDP values presented in Tab.I.In any case, in Fig. 3, the constraining power of SeaQuest data is certainly evident for x > 0.3 for both the NLO and NNLO analyses.However, for large x values the difference between the distributions of d(x) and ū(x) diminishes.Analogous observations can be made when examining Fig. 4, whose left panel illustrates the variation of the ratio d(x)/ū(x) with respect to x for µ = 3 GeV.The ratio is larger than unity within a large x interval, up to at least x < 0.5 -0.6.At these x values, both the xū and x d sea distributions are tiny, of the order of 10 −5 .The analysis incorporating SeaQuest data exhibits a high level of compatibility with the analysis that excludes them, and displays a smaller uncertainty band, especially for x > 0.3.This confirms the constraining role of the SeaQuest data.As shown in the right panel of Fig. 4, the results are also very well compatible with the d(x)/ū(x) ratio extracted by the SeaQuest collaboration at the scale Q 2 = 25.5 GeV 2 , which is characteristic of the kinematics of the experiment, using as a starting point the experimentally measured cross-section ratio σ pd /(2σ pp ) and Eqs. ( 8), (10) and (11) of Ref. [23].Although this extraction depends in principle on the PDFs used (the quoted SeaQuest values are those reported in Tab. 8 of Ref. [23], obtained using cross-sections computed with the CT18 PDF fit as input of their Eq.( 11)), this dependence is quite weak, i.e., it comes from subleading terms in Eq. ( 11) of Ref. [23], generating minor corrections to the leading result corresponding to the case x 1 ∼ x 2 .Therefore the extracted d/ū ratio can be considered as a robust quantity, as also already observed in Ref. [23] 7 .We also note that the Seaquest data cover target x values up to 0.45.The uncertainty band of the ABMP16 + SeaQuest PDFs remains small at even larger x values, which is a consequence of assumptions about the parameterization of these PDFs and their extrapolation to large x, performed under assumption of smoothness of the distributions.The same is true for the ABMP16 PDF fit.Only future experimental data in the large x region will be able to check the correctness of this extrapolated result shown here.Regardless, it is important to emphasize that the ABMP16 + SeaQuest fits rely on the identical PDF parameterization employed in the original ABMP16 fits.Remarkably, this parameterization already yielded a satisfactory fit to the new data, without necessitating any postadjustments through the introduction of additional parameters.During the original ABMP16 fit, we employed a strategy that involved investigating the impact of various functional forms while minimizing the number of parameters used.Our aim was to avoid introducing any additional parameters that did not contribute significantly to an improved description of the data.
We also point out that the effect of SeaQuest data, when comparing the ABMP16 PDFs to the ABMP16 + SeaQuest PDFs, is not dramatic, because the ABMP16 fits already included the E866 data, capable of constraining the d/ū ratio up to slightly lower x values than SeaQuest.The main addition of SeaQuest has been to have provided reliable measurements in the interval x ∼ 0.24 -0.45, which have helped to further constrain PDFs with respect to the past.

III. COMPATIBILITY OF SEAQUEST DATA WITH OTHER PDF FITS
The compatibility of the SeaQuest data with a number of modern PDF fits is shown in Fig. 5.The SeaQuest data align well with the predictions based on the NNPDF4.0fit [46], which is not surprising since the NNPDF collaboration incorporated these data into their fitting process.Nevertheless, the uncertainty range associated with this particular fit remains larger compared to our own uncertainty range, in contrast to the uncertainties accompanying the data.We argue that this behaviour can be ascribed to inefficiencies in the statistical estimators used in their analysis.
This issue seems to be confirmed also by the predictions on the d/ū ratio shown in Fig. 4 (right), where the constraining power of the SeaQuest data seems to be only partially reflected in the NNPDF4.0uncertainties.This is particularly visible in the region x ∼ 0.3 -0.45,where the NNPDF4.0uncertainties become large, although this region is still covered by SeaQuest data.The inclusion of the FNAL-E605 data in the NNPDF4.0fit should have imposed additional constraints on the d/ū ratio at these specific values of x.Consequently, one would expect the size of their 1σ band to be smaller compared to the result from their fit.The use of a large fixed number of parameters in the parameterization of these PDFs might be responsible for a relatively large uncertainty in the x region where SeaQuest data are present.The shape of the spike region around x ∼ 0.5 in Fig. 4 (right) seems to be driven by the step functions used in the parameterization of these PDFs.
The large uncertainties of NNPDF4.0 at larger x values in Fig. 4 (right), on the other hand, can be attributed to the lack of data.The smaller uncertainty of the ABMP16 fits (in comparison to NNPDF4.0) in the very large x region, not covered by the SeaQuest data, is neither related to the use of looser W cuts (W > 2 GeV) on the invariant mass of the hadronic system in DIS data 8 (M P is the proton mass), nor to the inclusion of higher-twist corrections in the fit.Moreover, it is important to note that these uncertainties cannot be considered highly informative.This is due to the fact that the uncertainty arises solely from the extrapolation beyond the region where data are available, relying on assumptions of smoothness, as already mentioned in the previous section.It is however true that the sum rules play a role in constraining the shape of PDFs there.We have checked that a shape with even more spikes and larger uncertainties for the d/ū ratio occurs in the case of the predecessor of NNPDF4.0 fit, i.e. the NNPDF3.1 PDF fit [47] (not shown in the plot), not including SeaQuest data.These concerns about the NNPDF4.0PDFs at large x are also manifest in unusual predictions for the forward-backward asymmetry A * FB in the invariant mass of the dilepton final state at the LHC, quite different from those of many others PDF fits particularly for large invariant masses [48,49].The measurement of this quantity and its comparison with theory predictions might be important for improved fits of large-x quark PDFs within the Standard Model (SM) and/or for discovering new physics associated to new gauge sectors beyond the SM, such as a heavy neutral Z ′ -boson, see, e.g.Refs.[50,51].
The SeaQuest data are also compatible with the CT18 fit [54] (not shown in Fig. 5), although the uncertainty of the latter looks particularly large, even due to the tolerance criterion used in this PDF fit (∆χ 2 = 100 at 90% C.L. roughly corresponding to ∆χ 2 ∼ 30 at 68% C.L., vs ∆χ 2 = 1 used in various other PDF fits adopting the Hessian approach, although not in all 9 ), and this prevents any strong conclusion.The CT18 collaboration has investigated the impact of first SeaQuest data of Ref. [43] on their NNLO PDFs in Ref. [56] and they have compared their predictions even to the BNL STAR data on W-boson production [57].Additionally, the CT18A variant of the fit, together with further variants incorporating lattice QCD data on the strangeness asymmetry distribution s(x) − s(x), have also been compared to first SeaQuest data of Ref. [43] in Ref. [20].An advanced study aiming at separating the so-called connected and disconnected sea components, reflecting the topology of the quark lines in the four-point current-current correlator in the nucleon, under the CT18 parameterization, has led to the CT18CS fit [58], using as a basis the original CT18 data sets.The CT18CS PDFs have also been compared with the distributions extracted from the SeaQuest data of Ref. [43], and older E866 data of Ref. [24].FIG.5: The pulls for SeaQuest data [23] on the ratio of pd and pp DY distributions over x F with respect to predictions obtained using the code VRAP [39] in combination with the NNLO ABMP16 PDFs.The 1σ band for prediction (right-tilted hash) is compared to the NLO ABMP16 [30] (left-tilted hash) and NNLO NNPDF4.0 [46] (shaded area) ones.The central values of predictions with other PDFs are shown for comparison (dots: NLO CJ15 [52], long dashes: NNLO epWZ16 [53], suggesting the SU(3)-symmetric quark sea, dashed dots: NNLO MSHT20 [19]).

SeaQuest (
On the other hand, a comparison of the SeaQuest data with predictions obtained with the MSHT20 [19] and CJ15 [52] fits, both shown in Fig. 5, reveals that the d/ū ratio according to the latter has a trend compatible with the data only in part of the (x 1 , x 2 ) range.The CJ collaboration has also investigated the impact of the first SeaQuest data of Ref. [43] plus the aforementioned STAR data on the CJ15 PDFs in Ref. [56] 10 and very recently proposed the new global PDF fit CJ22 in a follow-up paper [60], incorporating the SeaQuest data plus the aforementioned STAR data, including higher-twist effects and nucleon off-shell corrections.It would be interesting to study as well the modification of the MSHT20 fit, after inclusion of the SeaQuest data 11 .
Finally, we remark that the behaviour of the ( d/ū)(x) ratio predicted by the ATLAS 2016 fit [53] turns out to be incompatible with the SeaQuest data, systematically underestimating the latter, pointing to issues in the parameterization of these PDFs and/or shortcomings during the fit.In particular, the comparison of the x F distribution with the SeaQuest data in Fig. 5 confirms the point raised already in Ref. [16] that the assumptions concerning d-quark suppression with respect to the u-quarks in the ATLAS PDF parameterization adopted in that fit, now outdated, are problematic 12 .
The considerations in Ref. [16] were based on the observation that these PDFs already exhibited disagreement with the E866 data, which were already accessible at that particular time.One should in any case not be surprised that this old PDF fit is not in agreement with SeaQuest data, considering that, by definition, it did not include typical non-ATLAS datasets constraining high-x PDFs.In turns, this lack of data required to make more constraining assumptions on the PDF form.Newer ATLAS PDF fits have added more ATLAS data, partially extending x coverage and allowing for more flexible parameterizations.However, we have verified that even the central PDF from a more recent ATLAS fit, ATLASepWZVjet20-EIG [63] (not shown in our plot), including the W, Z/γ * +jet data that are sensitive to partons at larger x's than the inclusive W, Z/γ * data, turns out to be also incompatible with the SeaQuest data, overestimating the data up to several ten percent in the smallest x F bin (corresponding to the largest x).On the other hand, the central PDF from a most recent ATLAS PDF fit, ATLASpdf21 [64] (not shown in our plot), a fit that has included further data and also considered the role of scale uncertainties, largely overestimates the SeaQuest data in the first x F bin, but is compatible with the latter in the other bins, i.e. for 0.2 < x F < 0.8.The lack of agreement in the smallest x F bin can be probably attributed to the fact that ATLAS does not have data constraining ū(x) and d(x) for x > 0.3.On the other hand, the agreement visible at larger x F , corresponding to x < 0.3, remarks the compatibility between SeaQuest and DY ATLAS and Tevatron data.In Ref. [64] the ATLAS collaboration has provided their own comparison of ATLASpdf21 d/ū(x) ratio with that extracted by the NuSea and SeaQuest collaborations in Ref. [24,43].Considering that tha smallest x F correspond to the largest x values, our results and conclusions on compatibility between fixed-target and collider DY datasets are compatible with their ones.

IV. IMPACT OF NUCLEAR CORRECTIONS
SeaQuest data discussed and used in previous sections have been collected for a deuteron target and for this reason the analysis should address the corresponding nuclear corrections.Here we discuss the effect of nuclear corrections on the DY cross section following Ref.[65].This model addresses a number of mechanisms for nuclear corrections including the effect of nuclear momentum distribution (Fermi motion), nuclear binding, the off-shell modification of bound nucleon PDFs, as well as meson-exchange currents and nuclear shadowing corrections.For the kinematics of SeaQuest data the relevant corrections originate from nuclear momentum distribution, binding and off-shell effects on the PDFs.The deuteron PDFs q i/d of type i = u, d, . . .can be written as follows [65] (see also Appendix B of Ref. [61]) where q i/p(n) is the corresponding proton (neutron) PDFs, the integration is performed over the nucleon momentum k, Ψ d (k) is the deuteron wave function in the momentum space, which is normalized as 1, and M is the nucleon mass.We consider the deuteron in the rest frame and the z axis is chosen to be antiparallel to the momentum transfer.The four-momentum of the bound nucleon is k , where M d is the deuteron mass and k 2 = k 2 0 − k 2 is the invariant mass squared (virtuality), while x ′ = xM/(k 0 + k z ) is the Bjorken variable of the off-shell nucleon.
It is convenient to discuss the virtuality dependence of the nucleon PDFs in terms of the dimensionless variable v = (k 2 − M 2 )/M 2 .Since nuclei are weakly bound systems, the value of |v| is small on average.For this reason the off-shell PDFs can be expanded in a power series in v about v = 0 [66].Keeping the terms linear in v we have [67] where the derivative is taken for k 2 = M 2 .The function δ f (x, Q 2 ) measures the modification of the nucleon PDFs in the off-shell region.In Eq. ( 7), in order to simplify notations, we suppress the subscripts referring to the PDF type i.Also, we implicitly assume an average over the proton (p) and neutron (n), q i = (q i/p + q i/n )/2, since Eq. ( 6) for the deuteron depends only on this isoscalar PDF combination.Detailed studies of nuclear DIS, DY lepton-pair and W/Z boson production indicate that the data are consistent with an universal function δ f (x), independent of the parton type, and without significant scale and nucleon isospin dependencies [61,62,65,[67][68][69][70] 13 .In this work we use the results on the function δ f (x) from the recent analysis of some of us in Ref. [62].In Fig. 6 we illustrate the nuclear effects obtained for the valence quark PDFs, antiquark PDFs and the DY cross sections for the kinematics of the SeaQuest experiment.In particular, we show the ratios R val = u val/d /(u val/p + u val/n ), R sea = u d /(u p + u n ) and R DY = σ pd /(σ pp + σ pn ) computed using Eq. ( 6), the NNLO proton PDFs of Ref. [31] and the values of kinematical variables from Tab. 6 of Ref. [23].Note the different shapes of R val and R sea vs x F .This is caused by different x dependencies of the valence and antiquark nucleon PDFs and the smearing effect in the nuclear convolution, Eq. ( 6).The shape and the magnitude of R DY and R sea are similar corresponding to the fact that the DY cross sections σ pd , σ pp and σ pn for SeaQuest kinematics are dominated by the partonic contribution involving a proton beam valence u quark and a target ū, considering PDF x and flavour dependence.However, this dominance is violated for small values of x F (x F < 0.3), causing the different values of nuclear corrections for DY cross sections and the up-quark sea PDFs in this region.The magnitude of the nuclear corrections on the DY σ pd is extremely modest, typically O(0.5 − 1)%, and has a practically negligible impact on the present analysis.This result is consistent with the claim of Ref. [23] that nuclear corrections can be neglected, on the basis of the results of Refs.[71,72].
Nuclear corrections should also be addressed when dealing with data from FNAL-E605 experiment on proton-copper collisions [25].The corresponding corrections on the DY cross sections have been calculated in Ref. [65] (see Fig. 8 and Table 2 there).The rate of nuclear corrections depends on both the target x 2 and the mass of the muon pair as illustrated in Fig. 8 of Ref. [65].Note, however, that the E605 experiment only provides data for copper target and did not take data for the proton target.Since copper is almost an isoscalar target with about 8% of the neutron excess, the E605 cross section data on copper target alone provides a little sensitivity to measuring the ( d − ū)(x) asymmetry of the sea distributions.We also verified that removing the E605 data  The corresponding second Mellin moments from QCD lattice computations turn out to be significantly more uncertain and not yet able to discriminate between the various PDF fits.Taking into account the range of lattice results and the inherent uncertainty associated with each of them, they are presently exhibiting a high level of compatibility with nearly all the PDF fits.The lattice moments ⟨x⟩ u + by the χQCD collaboration and ⟨x⟩ u + − d + by the RQCD collaborations, both computed in 2018, exhibit a slight tension, deviating from their 1σ range, when compared to the majority of PDF fits.Nonetheless, these results carry substantial uncertainties and align with the findings of PDF fits within a 2σ range.Most of the recent lattice results, in particular those obtained by the ETMC collaboration in 2020, turn out to agree very well with almost all PDFs fits.The 2021 result on ⟨x⟩ u + − d + of the Mainz collaboration agrees with moments of some of the global PDF fits, but is slightly smaller, although compatible within 2σ, with the second moments from the ABMP16 (+ SeaQuest) NLO and NNLO PDFs.
We also observe that the addition of SeaQuest data to the ABMP16 NNLO PDF fit has a tiny effect on the values of the considered moments, slightly decreasing the associated uncertainties, while the central values remain approximately the same.The improvement of the uncertainties turns out to be more pronounced in the case of the ABMP16 NLO PDF fit.Overall, the results from NLO and NNLO ABMP16 (+ SeaQuest) PDFs are consistent among each other and, as mentioned, the order of perturbation theory does not have a significant impact on the second Mellin moments, being rather inclusive quantities.In light of the comparisons discussed here, it will be interesting to observe precise lattice calculations of higher Mellin moments (beyond the second/third ones), maybe exploiting concepts and techniques of Ref. [82,83], so as to enable similar comparisons for the fourth, etc. moments.Another valuable improvement would be the ability to distinguish between valence and sea quark PDFs in lattice results.

VI. CONCLUSIONS
We have studied a variant of the ABMP16 NLO and NNLO fits, including the SeaQuest nonresonant data on σ pd /(2σ pp ) as a function of x F .We find that these data reduce uncertainties on the ( d − ū)(x) difference as well as on the ( d/ū)(x) ratio at large x, while leaving essentially unchanged the values of the other quantities, which are simultaneously constrained in these fits (α s (M Z ) and heavy-quark masses).The χ 2 /NDP for the fits including SeaQuest data are within statistical uncertainty of those previously obtained without these data.The simultaneous description of all DY data turns out to be slightly more consistent at NNLO than at NLO, as expected for the improved precision of the theoretical predictions.In particular, we observe the compatibility of SeaQuest data constraints on the d − ū asymmetry, with the corresponding constraints from collider DY data at both the Tevatron and the LHC.This confirms the presence of an asymmetric sea, ruling out PDF fits based on the assumption of (or leading to) a symmetric sea.
Our present results support using the SeaQuest data together with collider DY data in future updated PDF analyses, that would allow further reducing PDF uncertainties, as well as a crosscheck of the compatibility with the data already included there.The inclusion of SeaQuest data is facilitated by the fact that nuclear corrections for the deuteron target, that we have explicitly computed in this work, turned out to be O(0.5 − 1)% in all SeaQuest x F bins, thus having a practically negligible effect on the final PDFs.The smallness of the observed nuclear effects can be attributed to the kinematics of the SeaQuest experiment itself.The experiment combines partons with relatively small but still significant x 2 values (specifically, x target ) and larger x 1 values (specifically, x beam ), with only the target experiencing nuclear corrections.The most substantial corrections occur in the bin with the smallest x F , which corresponds to the largest x 2 values (≤ 0.45).It is worth noting that larger nuclear corrections would be anticipated at larger x 2 values, which correspond to backward kinematics that fall outside the scope of SeaQuest's current detector capabilities.
The second moments of various combinations of light-quark distributions from NLO and NNLO PDF fits are compatible with current lattice QCD results.Although lattice QCD is not yet competitive for distinguishing between different PDF fits, advancements in techniques and increased attention from the lattice community are expected to improve this limitation in the future.
We strongly encourage the SeaQuest collaboration to continue their efforts in reducing the uncertainties associated with their measurements, aiming for values below the current level of approximately 5%.Achieving this would significantly enhance the constraining power of the data on sea quark distributions.It is worth noting that only around one half of the experiment's data has been utilized for the published studies thus far, suggesting the potential for further improvements.Moreover, it would be highly beneficial if the SeaQuest collaboration released separate data on pp and pd cross-sections.Such separate data sets would enable more precise constraints to be obtained for the ū and d quark distributions, facilitating a deeper understanding of their individual characteristics.
FIG. 2:The SeaQuest data[23] on the ratio of pd and pp DY distributions over Feynman variable x F with respect to the NNLO predictions obtained using the code VRAP[39] (solid line) and DYNNLO[40] (dashes).

FIG. 3 :FIG. 4 :
FIG.3:The 1σ band for the n f = 3-flavour isospin asymmetry of the sea distribution x( d − ū)(x) at the scale µ =3 GeV obtained in the present analysis (left-tilted hash) compared to the one of the ABMP16 fit (right-tilted hash).The left panel shows results of the NLO analysis, whereas the right panel refers to the NNLO one.

25 FIG. 7 :
FIG.7: Second Mellin moments of u + (x), d + (x) and the isovector combination (u + − d + )(x) and their uncertainties computed for a range of PDF fits and from lattice QCD.The corresponding numerical values are tabulated in the columns of Tab.III and reported in the panels of this plot for a more immediate visualization.The vertical band in each panel brackets the values from the ABMP16 NNLO fit.

TABLE I :
The total values of χ 2 obtained for the NLO and NNLO ABMP16 fits in comparison with the ones of the present analyses, including all data already considered in the ABMP16 fits plus SeaQuest data.See text for more detail.
[23]ear effects in the deuteron for the valence quark PDFs (R val ), antiquark PDFs (R sea ) and the DY cross sections (R DY ) (see text for more detail on the definition of these ratios) vs x F , computed using Tab.6 of Ref.[23].The upper horizontal axis indicates the corresponding x values of the deuteron target (x 2 ).

TABLE III :
Comparison of second Mellin moments for various combinations of light-quark distributions from different NLO and NNLO PDF fits, including those proposed in this work, with uncertainties due to PDF variations, to corresponding values extrapolated from n f -flavour lattice QCD computations (Q = 2 GeV).In the case of the CT18 fit, the uncertainties refer to the 90% C.L. interval, instead of the 68% one used by other Hessian PDF fits.
and uncertainties from NLO and NNLO PDF fits have almost comparable values, indicating very good perturbative stability.