Rapidity Distributions of Hadrons in Proton-Nucleus Collisions

We study proton-lead collisions with a new model for the Fock states of the incoming proton. The number of collisions which the proton experiences selects the appropriate Fock state of the proton which generates a multiple of pp-like rapidity distributions. We take as input the pp maximum entropy distributions shifting the respective cm-rapidities and reducing the available energies. A comparison with existing data at 5 TeV is made and results for 8 TeV are presented. We also explore the high multiplicity data in this model.


I. INTRODUCTION
The main ideas underlying high energy proton-nucleus collisions are well established. It is easiest to consider the system in the reference frame where the nucleus is at rest. Then high energy excitations in the fast incoming proton become degenerate with the ground state. Their life-time at sufficiently high energy is significantly longer than the nuclear dimension, so these excitations can be treated as Fock components of the proton. In the large N c approximation one can reorganise these excitations in a series of color neutral Fock states consisting of quark-antiquark pairs [1]. We consider only those Fock components of the proton, which are brought to mass shell by interactions, otherwise they remain being a virtual fluctuation of the proton. The n-th Fock state is actualized by n collisions with target nucleons. Data can be explained when in n collisions only (n + 1)/2 as many particles are formed as in a single pp collision [2][3][4]. How can one understand the phenomenon that the fragmentation products do not multiply n times, given the two facts that there are n collisions and that the fragmenting two strings formed in each collision overlap strongly in rapidity space?
In this paper we present a model which can explain this phenomenon. We assume that in proton-proton collisions and in qq-proton collisions the strings hadronize in a similar fashion, producing a symmetric maximum entropy distribution [5,6] in its respective rest frame. With the simplification that pions are produced in the same amount as gluons, the maximum entropy method yields a Bose-type distribution depending on rapidity y and transverse momentum p ⊥ . The maximum entropy distribution has three parameters, namely an effective transverse "temperature" λ(s) and a longitudinal softness w(s), both depending on the center-of-mass (cm) energy. It is symmetrical around the cm-rapidity y cm . Using the rapidity variable we can describe both hemispheres of the distribution: For pp-collisions the transverse phase space is homogeneously distributed over the area Empirical values for this area give sizes L ⊥ ≈ 1.3 fm for pion distributions [5,6]. Invoking parton-hadron duality the multiplicity of produced particles is obtained by integrating the light cone distribution over the respective phase space. Note that the relativistic measure dy = dx/x arises from the large spatial extension in longitudinal direction of the small x partons. With a gluon degeneracy factor g = 2(N 2 c − 1) [5] the total multiplicity becomes: We emphasize that a statistical understanding of the final state in heavy-ion collisions necessitates both a correct description of the momentum and configuration space distribution.
The maximum entropy model was originally conceived for symmetric pp or AA collisions.
It can be adapted to the pA configuration by keeping the symmetry of the fragmentation products with respect to the cm-rapidity which depends on the momentum of the qq substate in the fast proton colliding with a target nucleon.

NUCLEUS COLLISIONS
For proton-nucleus collisions we add the individual contributions of the different projectile Fock states, which lose coherence, interacting with the nuclear participants. This means that a 2-particle Fock state is actualized (brought to mass shell) when the incoming proton interacts with two target nucleons, or a 3-body component interacts with 3 target nucleons and so on. The elastic proton-nucleus cross section amplitude, can be calculated within the Glauber theory, as well as the quasielastic and reaction cross sections (e.g. see in [7]). Thus one can calculate the reaction cross section as where the profile function of the nucleus is obtained from the nuclear density integrating over z direction. Here we employ the optical approximation and neglect the interaction range in comparison with the nuclear radius for the sake of simplicity, more accurate, but lengthy expressions are well known [7].
Apparently the exponentials in Eq. (3) can be expanded as σ pA reac = n σ(n), where This expression is frequently misinterpreted as an n-fold collision of the projectile proton with n bound nucleons. However, the proton can interact inelastically only once, the further collisions of the proton debris may occur with cross sections, quite different from σ pN in .
Multiple inelastic interactions, involved in (4) should be interpreted as independent collisions of different constituents of the projectile Fock component inside the nucleus. However, the cross sections of those collisions are unknown, and the Glauber model alone is unable to predict them and the multiplicity distribution.
Inelastic processes are related through the unitarity relation to the forward elastic amplitude. The latter is given by the Glauber model as a sum of multiple N N elastic amplitudes.
However, the one-to-one correspondence between inelastic processes and different terms in the elastic amplitude is not known within the Glauber model. This problem was solved by Abramovsky, Gribov and Kancheli (AGK) [8], who formulated the unitarity cutting rules for the elastic amplitude. The magnitudes and signs of different multiple interaction terms in the expanded Glauber amplitude were found to be related to a sum of unitarity cuts of different number of elastic N N amplitudes (Pomerons) 1 The key point of the AGK cutting rules is the independence of the proton multi-Pomeron coupling on the number of cut Pomerons. This important result was proven in [8] within the old-fashioned parton model with short rapidity-range correlations.
The AGK cutting rules provide the weight factors (4) for the inelastic pA collision corresponding to n cut Pomerons. Notice that the particles produced on mass shell from n cut Pomerons have to share the total energy (see below), however, the proton multi-Pomeron coupling Eq. (4) is the same as in the elastic cross section (no cut Pomerons), therefore the cross sections σ pN in should be taken at the full collision energy. At 5 TeV we use σ in = 70 mb as pp-inelastic cross-section [4]. The sum over all n, i.e. over all possible collisions gives the total inelastic cross section.
To understand the multiplicity distribution we construct a model for the n-th Fock state in such a way that it reproduces the general features of the phenomenology in pA collisions.
The most important empirical result is the increase of the multiplicity ratio r between pA and pp collisions with the mean number n of collisions: 1 The N N elastic amplitude, which can be called effective Pomeron, already includes multi-Pomeron terms and is subject to the AGK cutting rules. Therefore the hadron multiplicity also depends on the number of cut Pomerons, related to the projectile Fock components, interacting with the same bound nucleon. The hadron multiplicity distribution in pA collisions will still be given by the convolution of the multiplicity distribution in each N N collision with the distribution of the number of collisions. Both are controlled by the AGK rules, but with different weight factors. 4 In order to specify the details of the model, however, more and different reaction channels would have to be analyzed. We assume that the n-th Fock state |n is a coherent superposition of n-particle states with different numbers of "valence-like" and "sea-like" constituents: The valence-like constituents v i carry sizable momentum fractions of the proton. The sealike constituents s i are very slow and do not contribute to the hadronization significantly.
The first "valence"-like constituent in the proton has the structure |qqq and all others constituents are |qq color neutral dipoles. In line with the maximum entropy assumption we assume that in the quantum state with m valence particles each particle carries the same Combining this multiplicity with the Glauber weights gives the empirical multiplicity ratio between pA-and pp-collisions, Eq. (5). This simple additivity is broken by energy conservation and rapidity shifts arising from the lower momenta in the substates. There have been more sophisticated proposals [9] about the distribution function of the different Fock states and how they interact with the target nucleons, but the necessary multi-parton distributions are really unknown and can only be guessed.
The eikonal formula still holds because Gribov inelastic shadowing corrections to the total inelastic cross section from excitations of the fast proton are known to be small [10]. It also has been shown in string motivated models [11] that the inelastic pp cross section is mainly given by the distance (R) of q and diquark in the proton, whereas the string thickness (a) gives only a small contribution. Consequently also the transverse string excitations only give a minor modification of the cross section since the transverse size a is much smaller than the extension of the string a R. using the correct energy s(m) and rapidity y cm (m) to get the pseudorapidity distribution dN (m) d 2 p ⊥ dη which is the basic building block of our model: The n-th Fock state |n interacting with n target nucleon can produce up to n times the proton-proton multiplicity with the probability The first component |1, s, s, s, ... gives a rapidity distribution dN (1) d 2 p ⊥ dη . The second component |1/2, 1/2, s, s, ... gives twice the rapidity distribution dN (2) d 2 p ⊥ dη and similarly for the higher components. The total sum of produced fragmentation products comes from the arithmetic sum of 1 + 2 + 3 + ... + n pp distributions which must be averaged, because each component has equal weight. Averaging gives the correct multiplicity ratio between pp and pA collisions observed for high energy collisions as discussed before. For the charged particle multiplicity we multiply with a factor 2/3: The summation is extended to n max =40, since higher contribution are negligible. To compare with the data we integrate over transverse momentum: In Fig. 1  from the ALICE collaboration [12] together with the theory. The theoretical calculation is close to the data, but slightly underestimates the asymmetry of the measured distribution in pseudo-rapidity. The neglected small sea substates in the Fock states would increase the theoretical result for positive rapidities and decrease the proton side.
A comparison of our theoretical prediction with the data of 8 TeV p-Pb collisions may present an additional test of our model. The theoretical multiplicity distribution at 8 TeV has the same shape as at 5 TeV, but one obtains a larger multiplicity due to the higher energy.
Notice that p ⊥ -dependence Eq. (19) does not exhibit the expected effect of broadening. This is not the problem of the Fock-state picture, but of the Glauber approximation. The unitarity cut of the Glauber elastic pA amplitude suggested by the AGK rules, assumes independent hadronization of the cut Pomerons. This assumption, known as Bethe-Heitler approximation, is subject to corrections for the effects of coherence, which is known as the effect of saturation. Namely, a hadron can be produced coherently by several cut Pomerons, leading to rising p 2 ⊥ . Numerically this is a rather weak effect, as was evaluated in [13]. Even weaker is its influence on the mean p ⊥ -integrated multiplicity, so we neglect this correction. Smallness of this correction is confirmed by successful calculations of the multiplicity which produces n times the proton multiplicity and consequently will increase the overall multiplicity by a factor two: In addition the selection of impact parameter will also give a higher multiplicity. In Fig. 3 we show the highest multiplicity we can get in our model for an impact parameter cut 0 < b < 1.76 fm which corresponds to a centrality of 0-5%. It is impossible to get the factor three enhancement only from a cut in impact parameter space. In our model the average Fock state (cf. Fig. 3) does not suffice to describe the data taken in the ATLAS experiment [15], whereas the selected Fock state gives good agreement in the positive rapidity region. There is an even stronger asymmetry in the data than in In a recent reference [16] predictions for the 8 TeV p+Pb run at the LHC are compiled and compared with each other. They include saturation approaches, Monte Carlo event generators and perturbative QCD-based calculations. In contrast to these basic approaches, which aim at an ab initio description of pp-collisions, our paper uses the phenomenological rapdity distribution of pp-collisions to generate the rapidity distribution in pA-collisions.
Using the maximum entropy distribution for pp-collisions, we can with the same ansatz obtain the rapidity distribution in the most and less central regions. We develop a new view on the fluctuations in the proton, which extends previous work [17] to study fluctuations in the proton wave function explaining the observed azimuthal asymmetries of produced particles. A study of color neutral components in the proton is mostly known from the pion cloud picture. One can interpret the valence qq states as the high energy manifestation of the pion cloud, its chiral partner and other mesons. In the fast proton the container and cloud of the quarks carry quanta with sizeable light cone momenta which can interact with the target nucleons. This is the main new ingredient of our presented model.
Several improvements have to be studied: (i) How does the focusing on the highest Fock state takes place as a function of the trigger in order to understand the variations of the rapidity distribution with centrality. (ii) How do the sea states enter the Fock state composition quantitatively? (iii) Can the maximum entropy distribution for string fragmentation fluctuate? In this paper we only use an average value for the transverse temperature λ, but one expects that fluctuations in λ may be important to boost the multiplicity which is proportional to λ 2 . In statistical physics it is well known that fluctuations of the "temperature" can lead to a power behavior of the high p ⊥ part of the transverse momentum distribution which is encoded in Tsallis distributions. (iv) New experimental tests could be pp-and AA-collisions analyzed along the same lines as outlined in the paper. These points will be studied in further work.