Proton structure in high-energy high-multiplicity p-p collisions

A few-body proton image, expected to be derivable from QCD in the renormalization group procedure for effective particles, is used within the Monte Carlo Glauber model to calculate the anisotropy coefficients in the initial collision-state of matter in high-energy high-multiplicity proton-proton interaction events. We estimate the ridge-like correlations in the final hadronic state by assuming their proportionality to the initial collision-state anisotropy. In our estimates, some distinct few-body proton structures appear capable of accounting for the magnitude of p-p ridge effect, with potentially discernible differences in dependence on multiplicity.


Introduction
Protons resist precise theoretical description of their internal dynamics in the Minkowski space-time for a long time by now. The simplest such picture, which is provided by the constituent quark model used to classify hadrons, is not precisely derived from QCD. The theory itself uses the Euclidean-space techniques that do not easily yield any real space-time image. In these circumstances, it is of interest to note that high-energy high-multiplicity proton-proton (pp) collisions may shed new light on the issue of proton structure. Namely, the numerous products in such collisions exhibit collective behavior that appears dependent on the initial state of colliding proton matter and the latter depends on the proton structure. Thus, the correlations among products in high-energy high-multiplicity pp collisions may report on the proton structure.
In particular, the CMS [1] and ATLAS [2] collaborations reported the collective flow in pp collisions that resembles the one observed in heavy-ion collisions [3; 4; 5]. Several authors discussed such flows in pp collisions [6; 7; 8; 9; 10; 11; 12; 13; 14; 15; 16; 17; 18] and we follow the insights of Ref. [6] in order to estimate the extent to which the high-energy high-multiplicity pp events are sensitive to the model of proton structure. Following the approach of Ref. [6] and the parameter choice such as in Ref. [15], used here, means making a strong assumption that the parton medium produced in the overlap region of pp collision at the LHC has similar hydrodynamical properties as that in heavy ion collisions at RHIC. The individual proton structures we consider are motivated by the general features of the renormalization group procedure for effective particles (RGPEP) in quantum field theory [19; 20].
We find that, the effective picture of a quark and diquark with a gluon flux between them produces a different dependence of eccentricity and triangularity on multiplicity than the three-quark picture with a star-like junction made of gluons does. According to this finding, the recent data for highenergy high-multiplicity events suggest a significant star-like gluon junction component in the proton structure. Our analysis also indicates a need for assessing the adequacy of the linear relationship used by us between asymmetries in the initial collision state, such as the eccentricity or triangularity, and the final state correlations in high-multiplicity events, such as the elliptic flow.
2 Proton structure in pp collisions As mentioned in Sec. 1, the ridge-effect in pp scattering can be described using the hydrodynamic evolution of the asymmetric state of matter that results from one proton's quark and gluon distribution suddenly colliding with another's. The asymmetric state is meant to evolve according to laws of hydrodynamics until it eventually turns into the detected particles that emerge through hadronization in the final state, in which they exhibit the angular correlation over a long-range in rapidity, called the ridge. The final state ridge-like correlations, such as the elliptic flow, are thus related to the initial stage of pp collision whose nature depends on the proton structure. One of the key issues is thus how to describe the proton structure using QCD in the Minkowski space-time. Fig. 1 Proton structure described using effective quarks of size s that is very much smaller than, smaller and comparable with the constituent quark size sc, with the single large circle indicating the volume available for effective gluons irrespective of their corresponding size, see Ref. [19].
Conceptually, we approach this issue using the RGPEP [19; 20], which is a candidate for providing the mathematical tools for describing protons as bound states of effective quarks and gluons of specific size s. The size parameter s plays the role of an arbitrary renormalization-group scale that can be adjusted to the physical process one wants to accurately describe in simplest possible terms. This condition means choosing the right variables for grasping the essence of physics most economically from the computational point of view. The scale dependence of the proton structure expected in the RGPEP is illustrated in Fig. 1 and the corresponding examples of the color structure are shown in Fig. 2. The expectation is based on the scale-dependent features of effective Hamiltonians, which imply the possibility that a relativistic bound-state eigenvalue problem can be equivalently written in terms of a few-body problem for sizable effective quarks and gluons instead of an infinite combination of bare point-like quarks and gluons of canonical QCD. Hence, the RGPEP provides the scheme in which the distribution of matter in proton can be imagined in terms of wave functions, or probability distributions for the effective quarks and gluons of size s.
The few-body picture of protons in QCD suggested by the RGPEP allows us to preliminarily model the proton structure in terms of shapes illustrated by two typical examples in Fig. 3 [21], knowing that such models can in future be verified in theory. We ask if the ridge effect can phenomenologically distinguish between the effective configurations.
We consider three types of configurations. The proton quark-diquark configuration, denoted by I and shown on the left-hand side of Fig. 3, is motivated by Refs. [17; 18]. It is a superposition of a few Gaussians that represent a quark, a diquark and gluons forming a tube in between. The three-quark configuration, denoted by Y and shown on the right-hand side of Fig. 3, is motivated by Ref. [19]. It is a superposition of Gaussians that represent three quarks and additional gluons forming the Y-shaped junction. The shape of Y configuration is kept fixed. In addition, we consider the Gaussian fluctuating three-quark configuration, denoted by G-f, which is the same as the Y configuration but with the shape parameters generated according to Gaussian probability distributions. Details of all configurations we consider are available in Refs. [21; 22]. 3 Asymmetries in the initial stage of a pp collision and the final state correlations Following Ref. [6], we adapt a simple Glauber model, widely used for modelling high energy nuclear collisions [23], to describe the density of binary partonic collisions in scattering of two systems A and B, which is a function of the coordinates x and y in the plane transverse to the colliding beams, the impact parameter b and the varying parameters Σ that identify the proton structure and its orientation in space. The coefficient σ gg denotes a parton-parton scattering cross-section, in our estimates on the order of 4 mb, and ρ denotes the three-dimensional parton distribution described in Sec. 2. Eccentricity 2 and triangularity 3 in the initial stage of pp collision are calculated using the formula [24] in which the curly brackets denote the expectation value and coordinates are parameterized as x = s cos φ, y = s sin φ. The number of collisions in an event is and the cross-section density in the impact parameter plane is The total pp cross-section is thus where P (Σ) is the probability density for proton configuration Σ. For any quantity Q, its expectation value in many collisions is We used randomly oriented proton configurations in the Monte Carlo generation of about 3 · 10 5 events for each proton model and estimated the averaged eccentricity 2 and triangularity 3 in the resulting samples. In our estimates, σ gg ∼ 4.3 mb [15] and σ pp ∼ 60 mb [25] required the number of scattering partons N g = 9 ± 2 to obtain agreement with data, assuming that the multiplicity N = αN coll and reproducing the value N = 30 [26] for charged particles by choosing α = 8 ± 3. Our minimum bias results for eccentricity and triangularity are shown in Fig. 4 It is visible in Fig. 4 that the initial stage of pp collision is characterized by different multiplicity dependence of the asymmetries for different proton structures. In collisions of quark-diquark (II) and Gaussian-fluctuating (G-f) structures, the asymmetries decrease with multiplicity above about N = 100, while in the collisions of tripod three-quark configurations (YY) the initial stage asymmetries persist or even increase above N = 100.
In order to relate the eccentricity and triangularity to data, we note that the observable multiplicity distributions in transverse momentum p T and pseudorapidity η are conventionally written as where the reaction plane angle Φ RP for colliding spherically symmetric distributions is illustrated in Fig. 5. In the actual events the angles Φ RP are determined from the particle distributions.
Assuming that the averaged minimal bias elliptic flow parameter v 2 is proportional to the minimal bias eccentricity parameter 2 with coefficient order 0.3 [15], we obtain v 2 ∼ 0.11, 0.14 and 0.09 for the proton models II, YY and G-f, respectively, while data indicates v 2 in the range 0.04-0.08 [5]. G-f II YY Fig. 4 Average eccentricity and triangularity obtained using Eq. (7) for three different types of proton structure in pp collisions. Green curves labeled II correspond to collisions of protons in configuration I in Fig. 3, and the red curves labeled YY correspond to collisions in configurations Y in Fig. 3. The blue lines labeled G-f correspond to the Gaussian fluctuating three-quark configuration in which the Y-type proton configurations appear with different shape parameters according to a Gaussian probability distribution [21]. Note the distinct multiplicity dependence in the case of YY configurations. obtained in Eq. (7) with a coefficient on the order of 0.3 [15]. The magnitudes of v2 resulting from different proton models are displayed with quark-diquark = II, triangular = YY.
More generally, taking into account that the initial-stage asymmetry parameters n in Eq. (2) and the final-state coefficients v n in Eq. (8) are small, and assuming that a set of averaged coefficients v n with different n depends approximately linearly on the set of averaged parameters n , one can infer the dependence of v n on multiplicity N from the dependence of n on N . Accordingly, Fig. 4 suggests that the averaged elliptic flow v 2 and higher correlations in high-energy high-multiplicity pp events may indicate which configurations of effective quarks and gluons in the proton structure are most likely to occur. Namely, only the YY configurations lead to 2 and 3 that do not fall off for N exceeding about 100. Actually, recent data [2] from ATLAS Collaboration for pp collisions with √ s ∼ 13 TeV show that v 2 is stable for large N . According to our analysis, this finding favors the configuration YY.
It should be observed that the long-range, near-side angular correlations in pp collisions at LHC energies can also be studied in terms of multiparton interactions [27]. Interactions of four partons to four partons and beyond allow for linking the ridge effect to the models [28] and theory [29] of double parton distributions and their light-front analysis [30]. A unified approach is hence greatly desired.

Conclusion
Simple model estimates suggest that the correlations among final-state particles in high-energy highmultiplicity pp collisions are sensitive to the proton structure. It is not excluded that future sophis-ticated calculations will identify some spatial features of protons through precise interpretation of experimental data on these correlations. The multidimensional linear relationship between initial-state asymmetries and final-state flow parameters that we assume in our estimates requires verification, e.g., in the hydrodynamic model. If confirmed, it would provide an efficient way for studying the proton structure through correlations in high-multiplicity pp collisions using the universal matrices determined by the nature of assumed underlying collision dynamics. Interpreting the most recent LHC data on multiplicity dependence of the elliptic flow coefficient v 2 using our simple estimates, suggests that protons may often occur in the configuration of three effective quarks connected by a Y -shaped gluon junction.