Deep inelastic scattering as a probe of entanglement: confronting experimental data

Parton distributions can be defined in terms of the entropy of entanglement between the spatial region probed by deep inelastic scattering (DIS) and the rest of the proton. For very small $x$, the proton becomes a maximally entangled state. This approach leads to a simple relation $S = \ln N $ between the average number $N$ of color-singlet dipoles in the proton wave function and the entropy of the produced hadronic state $S$. At small $x$, the multiplicity of dipoles is given by the gluon structure function, $N = x G(x,Q^2)$. Recently, the H1 Collaboration analyzed the entropy of the produced hadronic state in DIS, and studied its relation to the gluon structure function; poor agreement with the predicted relation was found. In this letter we argue that a more accurate account of the number of color-singlet dipoles in the kinematics of H1 experiment (where hadrons are detected in the current fragmentation region) is given not by $xG(x,Q^2)$ but by the sea quark structure function $x\Sigma(x,Q^2)$. Sea quarks originate from the splitting of gluons, so at small $x$ $x\Sigma(x,Q^2)\,\sim\, xG(x,Q^2)$, but in the current fragmentation region this proportionality is distorted by the contribution of the quark-antiquark pair produced by the virtual photon splitting. In addition, the multiplicity of color-singlet dipoles in the current fragmentation region is quite small, and one needs to include $\sim 1/N$ corrections to $S= \ln N$ asymptotic formula. Taking both of these modifications into account, we find that the data from the H1 Collaboration in fact agree well with the prediction based on entanglement.

Parton distributions can be defined in terms of the entropy of entanglement between the spatial region probed by deep inelastic scattering (DIS) and the rest of the proton. For very small x, the proton becomes a maximally entangled state. This approach leads to a simple relation S = ln N between the average number N of color-singlet dipoles in the proton wave function and the entropy of the produced hadronic state S. At small x, the multiplicity of dipoles is given by the gluon structure function, N = xG(x, Q 2 ). Recently, the H1 Collaboration analyzed the entropy of the produced hadronic state in DIS, and studied its relation to the gluon structure function; poor agreement with the predicted relation was found. In this letter we argue that a more accurate account of the number of color-singlet dipoles in the kinematics of H1 experiment (where hadrons are detected in the current fragmentation region) is given not by xG(x, Q 2 ) but by the sea quark structure function xΣ(x, Q 2 ). Sea quarks originate from the splitting of gluons, so at small x xΣ(x, Q 2 ) ∼ xG(x, Q 2 ), but in the current fragmentation region this proportionality is distorted by the contribution of the quark-antiquark pair produced by the virtual photon splitting. In addition, the multiplicity of color-singlet dipoles in the current fragmentation region is quite small, and one needs to include ∼ 1/N corrections to S = ln N asymptotic formula. Taking both of these modifications into account, we find that the data from the H1 Collaboration in fact agree well with the prediction based on entanglement.
PACS numbers: 13.60.Hb, 12.38.Cy In our paper [1] (see also [2,3]) we computed the von Neumann entropy of the system of partons resolved by deep inelastic scattering (DIS) at a given Bjorken x and momentum transfer q 2 = −Q 2 . We then proposed to interpret it as the entropy of entanglement between the spatial region probed by deep inelastic scattering and the rest of the proton. We found that in the small x, large rapidity Y regime, all partonic micro-states have equal probabilities -the proton is composed by an exponentially large number N of micro-states that occur with equal and small probabilities 1/N . This yields a simple relation between the entanglement entropy and the multiplicity of partons (dominated by gluons at small x): where N (x, Q 2 ) is an average number of color-singlet dipoles. In the region of small x the gluons dominate and N xG(x, Q 2 ) where xG(x, Q 2 ) is the gluon structure function 1 . Assuming that the multiplicity of produced hadrons is proportional to the multiplicity of color-singlet dipoles ("local parton-hadron duality" [4,5]), eq. (1) imposes a relation between the parton structure function (extracted from the inclusive cross section of DIS) and the entropy of produced hadrons; this relation can be directly tested in experiment. The comparison to the experimental data on hadron multiplicity distributions from CMS Collaboration at the LHC provided encouraging results [1,2].
However, recent dedicated experimental analysis performed by the H1 collaboration [6] shows a disagreement with Eq. (1) if one assumes N = xG(x, Q 2 ) (see Fig. 12 in Ref. [6] and the dotted curves in our Fig. 1). In this letter we demonstrate that the H1 data in fact are in a good agreement with our approach, once two important effects are taken into account. Both of them are implied by the kinematics of the H1 measurements that are performed in the current fragmentation region (see also Refs. [6][7][8][9]), and thus at moderate values of Bjorken x. First, since the experimental hadron multiplicities are not large, we need to take into account corrections of the order of 1/N to Eq. (1). Second, because the H1 data are concentrated in the current fragmentation region and not at very small x, we need to reconsider our claim in Ref. [1] that the multiplicity of the color-singlet dipoles is equal to xG(x, Q 2 ). Indeed, the hadrons produced in the current fragmentation region of DIS originate from the hadronization of the struck quark (a constituent of the color-singlet dipole) and the multiplicity of color-singlet dipoles is thus determined by the sea quark structure function, see Fig. 2. Therefore the correct relation between the number of dipoles and the experimentally measured entropy of hadrons (valid for large N ) is The sea quark and gluon distributions are related to each other at small x by (see Fig. 2) where the function C(α s ln Q 2 , α s , α s ln(1/x)) describes the splitting of the virtual photon into the quark-antiquark pair. In the region of very small x, as we will now show, C(ln Q 2 , α s , α s ln(1/x)) → Const, and so the sea quark distribution is proportional to the gluon one. Therefore, with a logarithmic accuracy at very small x we recover our original relation To show that C(ln Q 2 , α s , α s ln(1/x)) → Const at very small x, let us use the leading order DGLAP evolution equation [10] that gives for the sea quark structure function (see Fig. 2) where ξ = ln Q 2 . At small values of x, the gluon structure function takes the form (see Ref. [11] for a review) where N c is the number of colors, C F = N 2 c − 1 /2N c , and the BFKL kernel χ (γ) has the form: g (γ) in should be calculated from the initial conditions. Plugging Eq. (6) in Eq. (5) we can take the integrals over z and ξ . Taking the integral over γ in the diffusion approximation, and using the method of steepest descent , we obtain the following result: We evaluated Eq. (8) for small α s neglecting terms in C(α s ln Q 2 , α s , α s ln(1/x)) proportional to x. Note that Const does not depend on α s and numerically is about 0.3. Actually, the next -to-leading order correction increases the value of this Const and, as we can see from Fig. 1 (panel with 5 ≤ Q 2 ≤ 10 GeV 2 ), Const ∼ 1 in the NNLO fit for the region of smallest x. However, away from the region of very small x and in the kinematics of H1 measurement, the proportionality (8) is violated, and one should use the relation (2) instead of (4).
In addition, away from small x, where the multiplicity of color-singlet dipoles N is not large, one should take into account 1/N corrections to (4). In Refs. [1,3] it is shown that in QCD cascade the multiplicity distribution has the following form:  [6] on the entropy of produced hadrons in DIS [6] with our theoretical predictions, for which we use the sea quark distributions from the NNLO fit [12,13] to the combined H1-ZEUS data.
where N is the average multiplicity of color-singlet dipoles. The distribution (9) leads to the following von Neumann entropy: One can see that at large N we obtain S ln N , but corrections are sizable when N ≤ 10 (see Fig. 3). It should be noted that the distribution of Eq. (9) describes quite well the experimental hadron multiplicity distributions in proton-proton collisions (see Refs. [1][2][3]).
For comparison with the H1 experimental data [6](see Fig. 1), we first assume, following [1], that the hadron multiplicity is equal to the number of color-singlet dipoles. This assumption is based on "parton liberation" picture [5] and on the "local parton-hadron duality" [4]. For sea quark and gluon structure functions in Fig. 1 we use NNLO fit [12,13] to the combined H1 and ZEUS data.
One can see that our approach in fact describes the H1 data quite well -this is the first test of the relation between entanglement and the parton model in DIS enabled by the H1 analysis. We stress that once the data in the target fragmentation region at smaller value of x becomes available at the Electron-Ion Collider, one should be able to use xG(x, Q 2 ) in the relation (1), as it has been done in Refs. [1-3, 14, 15]. However, the general formula is given by Eq. (2) which at small x reduces to S = ln(xG(x, Q 2 )) since xΣ(x, Q 2 ) → xG(x, Q 2 ).