Evidence for the maximally entangled low $x$ proton in Deep Inelastic Scattering from H1 data

We investigate the proposal by Kharzeev and Levin of a maximally entangled proton wave function in Deep Inelastic Scattering at low $x$ and the proposed relation between parton number and final state hadron multiplicity. Contrary to the original formulation we determine partonic entropy from the sum of gluon and quark distribution functions at low $x$, which we obtain from an unintegrated gluon distribution subject to next-to-leading order Balitsky-Fadin-Kuraev-Lipatov evolution. We find for this framework very good agreement with H1 data. We furthermore provide a comparison based on NNPDF parton distribution functions at both next-to-next-to-leading order and next-to-next-to-leading with small $x$ resummation, where the latter provides an acceptable description of data.


Entanglement entropy
The proton is a coherent quantum state with zero von Neumann entropy. However it has been argued in [1,2] that when the proton wave function is observed in Deep Inelastic Scattering (DIS) of electrons and protons, this is no longer true. In DIS, the virtual photon, with momentum q and q 2 = −Q 2 its virtuality, probes only parts of the proton wave function, which gives rise to entanglement entropy, between observed and unobserved parts of the proton wave function, through tracing out inaccessible degrees of freedom of the density matrix. The resulting entanglement is then a measure of the degree to which the probabilities in the two subsystems are correlated; for other approaches where thermodynamical and momentum space entanglement entropy have been studied see [3][4][5][6][7][8][9][10][11]; for studies on Wehrl entropy [12] and jet entropy see [13]). Based on explicit studies of this entanglement entropy , both within a 1+1 dimensional toy model and leading order (LO) Balitsky-Kovchegov evolution [14][15][16], as well as entanglement entropy in conformal field theory, the authors of [1] conclude that DIS probes in the perturbative low x limit a maximally entangled state. With x = Q 2 /2p · q and p the proton momentum, the low x limit corresponds to the perturbative high energy limit, where Q 2 defines the hard scale of the reaction and sets the scale of the strong running coupling constant α s (Q 2 ) ≪ 1. The perturbative low x limit of [1] corresponds then to the scenario where parton densities are high, but not yet saturated and non-linear terms in the QCD evolution equations are therefore sub-leading. This is precisely the kinematic regime, where perturbative low x evolution of the proton is described through Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution, which resums terms [α s ln(1/x)] n to all order in α s ; it is this kinematic regime to which the results of [1] are supposed to apply at first.
The proposal that DIS probes in the low x limit a maximally entangled state is closely related to the emergence of an exponentially large number of partonic micro-states which occur with equal probabilities P n (Y ) = 1/⟨n⟩, with ⟨n(Y, Q)⟩ the average number of partons at Y = ln 1/x and photon virtuality Q. Entropy is then directly obtained as Assuming that the second law of thermodynamics holds for this entanglement entropy, the above expressions yields a lower bound on the entropy of final states hadrons S h through S h ≥ S(x, Q 2 ) [1]. "Local parton-hadron duality" [17] and the "parton liberation" picture [18] then suggest that partonic entropy coincides with the entropy of final state hadrons in DIS, see also the discussion for hadron-hadron collisions in [3]. The hadronic entropy can be further related to the multiplicity distribution of DIS final state hadrons. The latter has been obtained from HERA data in [3], which allows for a direct comparison of Eq. (1) to experimental data.
Confirmation of Eq. (1) is of high interest, since it links hadron structure to final state multiplicities through entropy. If confirmed, it provides an additional constraint on parton distribution functions (PDFs). Moreover, entropy is defined non-perturbatively and the proposed relation is therefore not necessarily limited to perturbative events, unlike PDFs. Last but not least, entropy is subject to quantum bounds [20][21][22] and through Eq. (1) such bounds translate directly on bounds on the number of partons in the proton [1]. This is of particular interest for the search for a saturated gluon state commonly called the Color Glass Condensate at collider facilities such as the Large Hadron Collider and the future Electron Ion Collider.
The explicit model calculations of [1] were based on solutions of purely gluonic LO low x evolution, where quarks appear only as a next-to-leading order (NLO) correction; it is therefore natural to assume that at first the total numbers of partons agrees with the number of gluons. In the following we find that for the kinematic regime explored at HERA, quarks are indeed sub-leading, but nevertheless numerically relevant for a correct description of data. We therefore propose in this letter that the average number of partons in Eq. (1) should be interpreted as the sum of the number of all partonic degrees of freedom, i.e. of quarks and gluons.
Our description is based on the NLO BFKL fit [23,24] (HSS). Initial conditions of the HSS unintegrated gluon distribution have been fitted to HERA data on the proton structure function F 2 and the HSS fit provides therefore a natural framework to verify the validity of Eq. (1) and its conjectured relation to the final state hadron multiplicity. Moreover, the HSS fit is directly subject to NLO BFKL evolution [25] and therefore provides a direct implementation of linear QCD low x evolution.

Results
To compare the HSS unintegrated gluon distribution to data, we need to determine first PDFs, which will yield the total number of partons through where g(x, µ F ) (Σ(x, Q)) denotes the gluon (seaquark) distribution function at the factorization scale µ F . To this end we use the Catani-Hautmann procedure [26] for the determination of high energy resummed PDFs. At leading order, the prescription is straightforward for the gluon distribution function, which is obtained as where µ F denotes the factorization scale which we identify for the current study with the photon virtuality Q, and F(x, k 2 ) the unintegrated gluon distribution, subject to BFKL evolution. To obtain the seaquark distribution, we require a transverse momentum dependent splitting function [26] , where k denotes the gluon momentum and ∆ = q − zk with q the t-channel quark transverse momentum; T F = 1/2. Note that this splitting function reduces in the collinear limit k → 0 to the conventional leading order DGLAP splitting function P qg (z) = The integrated seaquark distribution is then obtained as [26] xΣ Note that in [27,28] a corresponding off-shell gluon-to-gluon splitting function has been determined. Within the current setup, this would allow in principle for the determination of the gluon distribution at next-to-leading order. The use of this splitting function for the determination of the gluon distribution function at NLO has however not been worked out completely so far. Moreover, the HSS fit is based on a leading order virtual photon impact factors, which suggests the use of the leading order prescription Eq. (3) also for this study. The HSS unintegrated gluon density reads [29] F x, k 2 , Q = 1 k 2 Results are compared to the final state hadron entropy derived from the multiplicity distributions measured at H1 [3] whereĝ is an operator in γ space, withᾱ s = α s N c /π, N c the number of colors and χ(γ, Q, Q) the next-to-leading logarithmic (NLL) BFKL kernel which includes a resummation of both collinear enhanced terms as well as a resummation of large terms proportional to the first coefficient of the QCD beta function, see App. A for details. Eq. (3) and Eq. (12) is now used to calculate through Eq. (10) the partonic entropy Eq. (1); the result is then compared to H1 data [3]. To calculate entropy for the H1 Q 2 bins, we employ the following averaging procedure, The results of our study are shown in Fig. 2, where we evaluate all expressions for n f = 4 flavors. We find that the partonic entropy obtained from the total number of partons gives a very good description of H1 data [3] in case of the HSS fit. As anticipated in [1], the purely gluonic contribution is clearly dominant and amounts to approximately 80% of the total contribution; nevertheless the seaquark contribution is needed for an accurate description of H1 data. Given the approximations taken in the derivation of Eq. (10) as well as the possibility that sub-leading corrections are relevant for the determination of hadronic entropy form the multiplicity distribution, we believe that the above result provides an impressive confirmation of Eq. (10) and the results of [1] in general.
In [3] the data shown in Fig. 2 have been compared to Eqs. (1) and (10). Based on the original proposal of [1], only the gluon PDF has been used, for which the LO gluon distribution of the HERAPDF 2.0 set [30] has been chosen. While the use of a LO gluon PDF is somehow natural, since Eq. (1) does at the moment clearly not address questions related to collinear factorization at NLO and beyond, it is well known that the convergence of the gluon distribution is rather poor in the low x region; differences between the LO and NLO gluon amount up to 100% in the low x region, see e.g. Fig. 26 of [30]. While there are still noticeable differences between NLO and NNLO gluon distribution (of the order of 30% at x = 10 −4 ), one can nevertheless argue that the gluon distribution starts to converge beyond leading order and the values provide by the NLO gluon might be taken as a more realistic reflection of the true gluon distribution. To substantiate this point, we show in Fig. 2 also results based on an evaluation of Eqs. (1) and (10) with NNPDF collinear PDFs at NNLO [31]. We further show results obtained using NNLO NNPDF with next-to-leading logarithmic (NLL) low x resummation [32]. In both cases we assume µ F = Q. While both PDF sets allow for an approximate description of data and may therefore serve as an additional confirmation of the correctness of Eq. (10), a satisfactory description of the x-dependence is only possible using the low x resummed NNPDF PDF set, which provides a very good description of the shape, with a slight off-set in normalization.
A different description of these data has been provided in [33] which uses the sea quark distribution only. The authors use however for their LO BFKL description the quark-to-gluon splitting function instead of the required gluon-to-quark splitting. The former is enhanced in the low x limit and yields an incorrect sea quark distribution which is presumably of the order of the gluon distribution. We also could not reproduce the description which is based on the collinear NNLO sea quark distribution.

Conclusions
In this letter we followed the proposal of [1] to treat the low x proton in Deep Inelastic Scattering as a maximally entangled state with an entanglement entropy given as the logarithm of the average number of partons in the proton. Unlike [1,33] we interpret the total number of partons as the sum of quarks and gluon numbers, determined through the regarding PDFs. While we agree with [1,3,9] that the quark distribution is sub-leading in the low x limit, we find that the seaquark distribution provides a numerically relevant contribution of the order of 20%.
Our description is based on the determination of PDFs from an unintegrated low x gluon distribution, subject to BFKL evolution. For the numerical study, the HSS unintegrated gluon, which follows NLO BFKL evolution, has been used. Comparing our result with the final state hadron entropy extracted by the H1 collaboration [3], we find a very good agreement with data, if the total number of partons is taken as the sum of gluons and sea-quarks. We also provided a comparison based on NNLO PDF sets by the NNPDF collaboration. While purely NNLO DGLAP PDFs provide only an approximate description of data, we find that NNLO DGLAP PDFs with NLL low x resummation provide a reasonable description of the slope of H1 data, which emphasis again the role of low x dynamics for the determination of the proton as a maximally entangled state of partons.
Note that such an agreement is not obtained if the comparison is based on leading order collinear PDFs, as used by the H1 collaboration. This clearly hints at the need to further refine the underlying theoretical framework, in particular to clarify in a systematic way the relation between entropy and PDFs within the framework of collinear and/or high energy factorization. This need is immediately apparent if Eq. (10) is evaluated using PDFs beyond leading order, which immediately implies a scheme dependence of the extracted parton number; strictly speaking the latter can be therefore no longer related to the hadron multiplicity which is a physical observable and therefore scheme independent. The description based on NNLO PDFs and NNLO low x resummed PDFs is therefore an approximation at best. Note that a similar limitation does not apply to the description based on the HSS fit, since the resulting PDFs are leading order, from the point of view of collinear factorization and therefore scheme dependent, while similar issues arise due to the use of high energy factorization beyond leading order in that case. Moreover the relation to other frameworks as studied in [3][4][5][6][7][8][9] needs to be clarified. Furthermore it will be interesting to explore possible deviations from this framework at lower values of Q and x due to the onset of nonlinear low x evolution, in particular effects due to saturated parton densities [34,35]. process, while M sets the scale of the running coupling constant. For the current study we set M = M = Q and n f = 4 with Λ QCD = 0.21 GeV. Q 0 = 0.28 GeV, and δ = 6.5. have been determined from a fit to the F 2 structure function in [23]. In this fit the overall running coupling constant has been evaluated at the renormalization scale µ 2 = QQ 0 , with Q the photon virtuality. For the construction of parton distribution function µ 2 = Q 2 is however a more natural choice. We therefore reevaluated the underlying fit and found that data on the proton structure F 2 [37] are equally well described, if we use µ 2 = Q 2 for the photon impact factor with a normalization C = 4.31. It is then this convention which we use in this study.  Figure 2: Partonic entropy versus Bjorken x, as given by Eq. (11) and Eq. (10). We furter show results based on the gluon distribution only as well as on quarks and gluons together. Results are compared to the final state hadron entropy derived from the multiplicity distributions measured at H1 [3] There was a mistake in the scale choice of the running coupling in the gluon density that was was used in the paper [1]; the mistake has been already corrected in [2]. The mistake was difficult to spot since the formulas that we used did not account for the fact that only the charged hadrons were measured. The numerical factors approximately canceled and the net result is only slightly changed. The number of partons in the corrected formulas is [2]: where g(x, µ F ) (Σ(x, Q)) denotes the gluon (seaquark) distribution function at the factorization scale µ F and as described above the factor 2/3 takes into account the fact that only charged partons were measured, see also the more detailed discussion in [2]. To calculate entropy for the H1 Q 2 bins, we employ the following averaging procedure, The corrected results are shown in Fig. 2.
There was also a typo in Eq. (5) of [1] which yields the formula for our determination of the sea quark distribution. The argument of the unintegrated gluon distriubtion has been given as x instead of x/z on the LHS of this equation. The corrected formula, which was actually used in the calculation reads