Constraints on nuclear parton distributions from dijet photoproduction at the LHC

TL;DR

This study demonstrates that dijet photoproduction data from LHC ultraperipheral collisions can significantly reduce uncertainties in nuclear parton distribution functions, especially at small momentum fractions, using Bayesian reweighting techniques.

Contribution

It introduces a method to incorporate LHC dijet photoproduction pseudo-data into nPDFs, showing substantial uncertainty reduction at small x values.

Findings

Uncertainty in nCTEQ15 nPDFs reduces by about 50% at x=10^{-3} with 5% pseudo-data error.

Reweighting has a smaller impact on EPPS16 nPDFs due to different parametrization.

Pseudo-data can constrain nPDFs effectively at small momentum fractions.

Abstract

Using QCD calculations of the cross section of inclusive dijet photoproduction in Pb-Pb ultraperipheral collisions in the LHC kinematics as pseudo-data, we study the effect of including these data using the Bayesian reweighting technique on nCTEQ15, nCTEQ15np, and EPPS16 nuclear parton distribution functions (nPDFs). We find that, depending on the assumed error of the pseudo-data, it leads to a significant reduction of the nPDF uncertainties at small values of the momentum fraction $x_A$. Taking the error to be 5\%, the uncertainty of nCTEQ15 and nCTEQ15np nPDFs reduces approximately by a factor of two at $x_A=10^{-3}$. At the same time, the reweighting effect on EPPS16 nPDFs is much smaller due to the higher value of the tolerance and a more flexible parametrization form.

Tables

Table 1. Table 1: The effective number of replicas N eff subscript 𝑁 eff N_{\rm eff} for different choices of the experimental error and sets of nPDFs.

$ϵ$	$N_{\mathrm{eff}}$(nCTEQ15)	$N_{\mathrm{eff}}$(nCTEQ15np)	$N_{\mathrm{eff}}$(EPPS16)
0.05	4407	3982	5982
0.1	7483	7742	8727
0.15	8870	9107	9555
0.2	9464	9607	9818

Equations

$$f_{j/A}^k(x,Q^2)=f_{j/A}^0(x,Q^2)+\frac{1}{2}\sum _{i=1}^N\left[f_{j/A}^{i+}(x,Q^2)-f_{j/A}^{i-}(x,Q^2)\right]R_{ki},$$

$$\frac{d{\sigma }^k}{d{x}_A}=\sum _{a,b}{\int }_{y_{\min }}^{y_{\max }}dy{\int }_0^{1}d{x}_{\gamma }{f}_{\gamma /A}(y)f_{a/\gamma }(x_{\gamma },{\mu }^2)f_{b/B}^k(x_A,{\mu }^2)d\hat{\sigma }(ab\to \mathrm{jets}).$$

$$\frac{d{\sigma }^0}{d{x}_A}=\sum _{a,b}{\int }_{y_{\min }}^{y_{\max }}dy{\int }_0^{1}d{x}_{\gamma }{f}_{\gamma /A}(y)f_{a/\gamma }(x_{\gamma },{\mu }^2)f_{b/B}^0(x_A,{\mu }^2)d\hat{\sigma }(ab\to \mathrm{jets}),$$

$${\chi }_k^2=\sum _{j=1}^{N_{\mathrm{data}}}\frac{(d{\sigma }^0/d{x}_A-d{\sigma }^k/d{x}_A{)}^2}{{\sigma }_j^2},$$

$$w_k=\frac{e^{-\frac{1}{2}{\chi }_k^2/T}}{\frac{1}{N_{\mathrm{rep}}}{\sum }_i^{N_{\mathrm{rep}}}{e}^{-\frac{1}{2}{\chi }_i^2/T}},$$

$${⟨\frac{d\sigma }{d{x}_A}⟩}_{\mathrm{new}}$$

$$\delta {⟨\frac{d\sigma }{d{x}_A}⟩}_{\mathrm{new}}$$

$$⟨f_{j/A}(x,Q^2){⟩}_{\mathrm{new}}$$

$$\delta ⟨f_{j/A}(x,Q^2){⟩}_{\mathrm{new}}$$

$$N_{\mathrm{eff}}=\exp \left[\frac{1}{N_{\mathrm{rep}}}\sum _k^{N_{\mathrm{rep}}}{w}_k\ln (N_{\mathrm{rep}}/w_k)\right].$$

Full text

††journal: Eur. Phys. J. C

1 Introduction

Collinear nuclear parton distribution functions (nPDFs) are fundamental quantities of Quantum Chromodynamics (QCD) encoding information on the one-dimensional distributions of quarks and gluons in nuclei in terms of the light-cone momentum fraction $x_{A}$ at a given resolution scale $\mu$. Nuclear PDFs are essential ingredients of QCD calculations of cross sections at high energies involving charged lepton–nucleus and neutrino–nucleus deep inelastic scattering (DIS) with fixed targets and – in the future – in the collider mode and proton–nucleus and nucleus–nucleus scattering at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). While nPDFs are non-perturbative quantities, which cannot be calculated from first principles of QCD, the QCD collinear factorization for hard processes and the Dokshitzer–Gribov–Lipatov–Altarelli–Parisi (DGLAP) evolution equations allow one to set up a framework of global QCD fits, which enables one to extract nPDFs from available data deFlorian:2003qf ; Hirai:2007sx ; Eskola:2009uj ; deFlorian:2011fp ; Kovarik:2015cma ; Khanpour:2016pph ; Eskola:2016oht . Different analyses give noticeably different predictions for nPDFs and carry significant uncertainties originating mostly from the limiting kinematic coverage of the available data, indirect determination of the gluon nPDF from the DIS data using the scaling violations, and different assumptions about the shape of nPDFs at the input scale. As a result, quark nPDFs for small $x$ and the gluon nPDFs for essentially all $x$ are rather poorly known.

Further progress in constraining nPDFs relies on studies of high-energy hard processes with nuclei at collider energies, notably, in proton–nucleus ($pA$) scattering at the LHC Salgado:2011wc ; Paukkunen:2018kmm ; Citron:2018lsq and lepton–nucleus ($eA$) scattering at the future EIC Accardi:2012qut ; Aschenauer:2017oxs and LHeC AbelleiraFernandez:2012cc . However, the QCD analyses of the data on various hard processes in $pA$ scattering at the LHC during Run 1 Eskola:2013aya ; Helenius:2014qla ; Armesto:2015lrg ; dEnterria:2015mgr showed that the data provide only modest restrictions on nPDFs at small $x$. At the same time, it was proposed Brandt:2014vva that measurements of low-mass lepton pair production in proton–lead collisions at the LHC has a large potential to reduce the theoretical uncertainties on nPDFs in a wide range of $x$ or even rule out some parameterizations. While the potential of hard $pA$ scattering at the LHC will certainty continue to be explored, see, e.g. Eskola:2018sxu , it is topical to study complementary probes of nPDFs at the LHC.

It has been realized that collisions of ultrarelativistic ions at large impact parameters, when the strong interaction is suppressed and the ions interact electromagnetically via the emission of quasi-real photons in so-called ultraperipheral collisions (UPCs), give an opportunity to study photon–photon, photon–proton, and photon–nucleus scattering at unprecedentedly high energies Baltz:2007kq . This program was realized during Run 1 at the LHC by measuring exclusive photoproduction of charmonia ($J/\psi$ and $\psi(2S)$ vector mesons) in Pb-Pb UPCs at $\sqrt{s_{NN}}=2.76$ TeV by the ALICE Abbas:2013oua ; Abelev:2012ba ; Adam:2015sia and CMS Khachatryan:2016qhq collaborations, which probes the small-$x$ gluon distribution of the target Ryskin:1992ui . The analyses Guzey:2013xba ; Guzey:2013qza of these data in leading-order (LO) QCD gave first direct and weakly model-dependent evidence of large nuclear gluon shadowing down to $x\approx 10^{-3}$, which agrees very well with the predictions of the leading twist nuclear shadowing model Frankfurt:2011cs and the the EPS09 Eskola:2009uj , nCTEQ15 Kovarik:2015cma , and EPPS16 Eskola:2016oht nPDFs. Note that in next-to-leading order (NLO) perturbative QCD, corrections for this process are large Ivanov:2004vd ; Jones:2015nna and the relation between the gluon PDF and the generalized gluon PDF is model-dependent, which makes it challenging to include the data on $J/\psi$ photoproduction on nuclei into global QCD fits of nPDFs.

During Run 2 at the LHC, in addition to light and heavy vector meson photoproduction in UPCs Klein:2017vua ; Kryshen:2017jfz ; Guzey:2016piu , the ATLAS collaboration for the first time measured inclusive dijet photoproduction in Pb-Pb UPCs Atlas . Predictions for rates of this process in $pA$ and nucleus–nucleus ($AA$) UPCs at the LHC in LO QCD with an emphasis of heavy quark production were made in Ref. Strikman:2005yv . It was found that the rates are very large allowing one to probe deeply into the small-$x$ region. At NLO pQCD, the cross section of inclusive dijet photoproduction in Pb-Pb UPCs in the ATLAS kinematics was calculated in Guzey:2018dlm . It was shown that the used theoretical framework provides a good description of various kinematic distributions measured by the ATLAS collaboration and that the calculated dijet photoproduction cross section is sensitive to nuclear modifications of nPDFs at the level of 10 to 20%.

In this work, we explore the potential of inclusive dijet photoproduction in Pb-Pb UPCs in the LHC kinematics to provide new constraints on nPDFs. In particular, using the results of our NLO QCD calculations of the cross section of this process Guzey:2018dlm as pseudo-data, we study the effect of including these data using the Bayesian reweighting technique Armesto:2013kqa ; Paukkunen:2014zia ; Kusina:2016fxy on nCTEQ15, nCTEQ15np, and EPPS16 nPDFs. We find that depending on the assumed error on the pseudo-data, it leads to a significant reduction of uncertainties of nPDFs at small $x_{A}$. For instance, taking the error to be 5%, we find that the uncertainty of quark and gluon nCTEQ15 and nCTEQ15np nPDFs reduces by approximately a factor of two at $x_{A}=10^{-3}$. The reweighting effect on EPPS16 nPDFs is much smaller due to the higher value of the tolerance and a more flexible parametrization form used in that analysis.

The remainder of this paper is structured as follows. In Sec. 2, we summarize key steps of the Bayesian reweighting method and define our reweighting procedure. We present and discuss our results in Sec. 3 and draw conclusions in Sec. 4.

2 Reweighting of the dijet photoproduction cross section

To quantity the power of inclusive dijet photoproduction in Pb-Pb UPCs at the LHC to constrain nPDFs, we use the standard Bayesian reweighting procedure outlined in the literature Armesto:2013kqa ; Paukkunen:2014zia ; Kusina:2016fxy . In detail, starting with $2N$ error sets of nPDFs ($N=16$ for nCTEQ15 Kovarik:2015cma and $N=20$ for EPPS16 Eskola:2016oht ), one generates $N_{\rm rep}=10,000$ replicas labeled by the index $k$ as follows

[TABLE]

where $f_{j/A}^{0}$ and $f_{j/A}^{i\pm}(x,Q^{2})$ are the central fit and the plus and minus error nPDFs corresponding to the eigenvector direction $i$ and $R_{ki}$ is a random number from the normal distribution centered at zero with the standard deviation of unity.

Next, for each PDF replica, one calculates the observable of interest, which in our case is the dijet photoproduction cross section as a function of $x_{A}$ ($x_{A}$ is the hadron-level estimate for the momentum fraction carried by the interacting nuclear parton) Guzey:2018dlm :

[TABLE]

In Eq. (2), $a,b$ are parton flavors; $f_{\gamma/A}(y)$ is the flux of equivalent photons emitted by ion $A$, which depends on the photon light-cone momentum fraction $y$; $f_{a/\gamma}(x_{\gamma},\mu^{2})$ is the PDF of the photon for the resolved-photon contribution, which depends on the momentum fraction $x_{\gamma}$ and the factorization scale $\mu$; $f_{b/B}(x_{A},\mu^{2})$ is the nuclear PDF with $x_{A}$ being the corresponding parton momentum fraction; and $d\hat{\sigma}(ab\to{\rm jets})$ is the elementary cross section for production of two- and three-parton final states emerging as jets in hard scattering of partons $a$ and $b$. The sum over $a$ involves quarks and gluons for the resolved photon contribution and the photon for the direct photon contribution dominating at $x_{\gamma}\approx 1$. The integration limits are determined by the rapidities and transverse momenta of the produced jets, see Guzey:2018dlm for details. Note that since the inclusive dijet cross section is linear in the nPDFs, which in turn are linear in $R_{ki}$, it is sufficient to evaluate it $2N$ times for each error nPDF.

The essence of the reweighting procedure is the calculation of statistical weights for each replica $w_{k}$, which quantify how well the calculation using Eq. (2) reproduces data or pseudo-data. In our case, for the pseudo-data, we use the results of our calculation of the dijet cross section Guzey:2018dlm

[TABLE]

where $f_{b/B}^{0}(x_{A},\mu^{2})$ corresponds to the central value of the nCTEQ15 nPDFs. Then, the corresponding chi-squared $\chi_{k}^{2}$ is

[TABLE]

where the sum runs over the pseudo-data points; $\sigma_{j}$ is the assumed uncertainty of the pseudo-data. In our case, $N_{\rm data}=9$ corresponding to different bins in $x_{A}$ and we take $\sigma_{j}=\epsilon d\sigma^{0}/dx_{A}$ with $\epsilon=0.05$, 0.1, 0.15, and 0.2. We assume that these errors account only for the statistical uncertainty and that bin-by-bin correlations are neglected. The values of $\epsilon$ span the range of typical uncertainties of measurements of high-$E_{T}$ dijet photoproduction at HERA Chekanov:2007ac .

Finally, with the help of $\chi_{k}^{2}$, one can introduce the weights $w_{k}$ using the following relation:

[TABLE]

where $T$ is the tolerance associated with a given set of nPDFs. Note that $\sum_{k}w_{k}=N_{\rm rep}$. In our analysis, we use $T=35$ for nCTEQ15 and nCTEQ15np Kovarik:2015cma and $T=52$ for EPPS16 Eskola:2016oht .

Using the weights $w_{k}$, one can calculate the new, weighted average cross section and its error:

[TABLE]

Similarly, one can evaluate the reweighted nPDFs and their uncertainties:

[TABLE]

Equations (6) and (7) quantify the effect of the pseudo-data on the calculation of the cross section of inclusive dijet photoproduction and the central value and uncertainties of nPDFs, respectively.

The effective number of replicas contributing to Eqs. (6) and (7) can be estimated using the following expression:

[TABLE]

Table 1 summarizes our values of $N_{\rm eff}$ for $\epsilon=0.05$, 0.1, 0.15, and 0.2 and nCTEQ15, nCTEQ15np, and EPPS16 nPDFs.

3 Results of the reweighting

Using the procedure outlined in Sec. 2, we perform the Bayesian reweighting of the pseudo-data on the cross section of inclusive dijet photoproduction on nuclei in Pb-Pb UPCs in the LHC kinematics. Our results are shown in Figs. 1–4 for nCTEQ15 nPDFs, Figs. 5–8 for nCTEQ15np nPDFs, and Figs. 9–11 for EPPS16 nPDFs.

Figures 1, 5, and 9 show the dijet cross section as a function of $x_{A}$: the pseudo-data points labeled “nCTEQ15” and given by black open squares are the results of the calculation using the central nCTEQ15 fit; red crosses with the associated error bands are the results of the calculations using a given set of nPDFs (the crosses coincide with the open squares in Fig. 1, and, hence, are not shown); finally, the blue filled circles and the associated error bands show the reweighted cross section and its uncertainty, see Eq. (6). The four panels correspond to our four choices of the assumed error $\epsilon$. One can see from these figures that while the reweighting does not noticeably change the central values of the cross section, it reduces its theoretical uncertainty: the effect is largest for the smallest $\epsilon$ and the first small-$x_{A}$ bin. 111This bin has larger statistical uncertainties, which can however be reduced by increasing the precision of the Monte Carlo integration Guzey:2018dlm .

The remaining figures (Figs. 2–4, 6–8, and 10–11) demonstrate the effect of the reweighting on uncertainties of nPDFs: different panels show uncertainty bands of nPDFs normalized to their central value, i.e., the bands spanned by $1\pm\delta f_{j/A}(x,Q^{2})/f_{j/A}(x,Q^{2})$, see Eq. (7), for the gluon, $u$-quark, $d$-quark, and $s$-quark nPDFs before (red, outer band) and after (blue, inner band) the reweighting as a function of the momentum fraction $x_{A}$ at $Q^{2}=400$ GeV2. This is a characteristic value of $Q^{2}$ probed in dijet photoproduction in the ATLAS kinematics. While the central values of nPDFs are essentially not affected by the reweighting, the uncertainty bands for nCTEQ15 and nCTEQ15np are noticeably reduced. As expected, the effect is largest at $\epsilon=0.05$ and much smaller at $\epsilon=0.15$ and $\epsilon=0.2$. (Since the reduction of the uncertainty bands is very similar in the $\epsilon=0.15$ and $\epsilon=0.2$ cases, we only show the results for the former.) For instance, the uncertainty in the gluon and quark nPDFs at $x_{A}=0.001$ reduces by approximately a factor of two. It is interesting to note that the uncertainty of the small-$x_{A}$ gluon distribution in the case of nCTEQ15np after the reweighting is similar to that of nCTEQ15 before the reweighting – it is of the order of 15% in both cases. Therefore, dijet photoproduction should have a similar impact on nCTEQ15 nPDFs as the RHIC inclusive pion production data, which was included in nCTEQ15 and excluded in nCTEQ15np. The advantage of dijet photoproduction is that it does not involve the pion fragmentation functions, which necessarily brings an additional uncertainty in analyses of nPDFs.

In the case of EPPS16 nPDFs (see Figs. 10 and 11), the effect of reweighting is much smaller due to several reasons. First, these nPDFs have been obtained with a higher value of the tolerance $T$, which allows for significantly more replicas to contribute to the reweighted quantities (see Table 1) and which reduces the effectiveness of the reweighting. Second, a more flexible form of the EPPS16 nPDF parametrization also significantly reduces the reweighting effect, which is negligibly small in the $\epsilon=0.15$ and $\epsilon=0.2$ cases. We therefore do not show them here, since the blue and red solid lines and bands completely overlap.

Note that in typical fits of nPDFs, one parametrizes the dependence of the fit parameters on the nuclear mass number $A$ Eskola:2009uj ; Kovarik:2015cma ; Eskola:2016oht , which hence correlates these parameters for different nuclei. While, by construction, nuclear modifications of nPDFs and their uncertainties decrease with a decrease of $A$, the reduction of uncertainties of nPDFs for Pb due to the considered reweighting should also reduce uncertainties of nPDFs for lighter nuclei; the magnitude of the effect depends on numerical values of the fit parameters.

4 Conclusions

In this work, we studied the potential of inclusive dijet photoproduction in Pb-Pb UPCs in the LHC kinematics to give new constraints on nPDFs. Using the results of our NLO QCD calculations of the cross section of this process as pseudo-data, we analyzed the effect of including these data using the Bayesian reweighting technique on nCTEQ15, nCTEQ15np, and EPPS16 nPDFs. We found that depending on the assumed error on the pseudo-data, it leads to a significant reduction of the nPDF uncertainties at small $x_{A}$. For instance, taking the error to be 5%, we find that the uncertainty of quark and gluon nCTEQ15 and nCTEQ15np nPDFs reduces by a factor of two at $x_{A}=10^{-3}$. We observed that the uncertainty of the small-$x_{A}$ gluon distribution in the nCTEQ15np case after the reweighting is similar to that of nCTEQ15 before the reweighting, which indicates that dijet photoproduction should have a similar impact on nCTEQ15 nPDFs as the RHIC inclusive pion production data with the advantage that dijet photoproduction is free from the uncertainty associated with the pion fragmentation functions. At the same time, the reweighting effect on EPPS16 nPDFs is much smaller due to the higher value of the tolerance and a more flexible parametrization form used in the EPPS16 analysis.

Acknowledgements.

VG would like to thank K. Kovarik and H. Paukkunen for useful discussions. VG’s research is supported in part by RFBR, research project 17-52-12070. The authors gratefully acknowledge financial support of DFG through the grant KL 1266/9-1 within the framework of the joint German-Russian project “New constraints on nuclear parton distribution functions at small $x$ from dijet production in $\gamma A$ collisions at the LHC".