Comparing Two Different Initial Conditions AAMQS and Quartic Action in Illustrating the Effect of Valence Color Charges in High-Energy pA Collisions at Forward Rapidity

The inclusive particle productions in proton-proton (pp) and deuton-gold (d+Au) collisions at forward rapidity at the Relativistic Heavy Ion Collider (RHIC) energy are studied in the framework of the color glass condensate (CGC) theory by using two different initial conditions: AAMQS (Albacete-Armesto-Milhano-Quiroga-Salgado) and quartic action. Then, the results obtained by the two different initial conditions in illustrating the effect of valence color charges in high-energy proton-nucleus (pA) collisions at forward energy are compared. Meanwhile, the inclusive particle productions in pA collisions at forward rapidity at the Large Hadron Collider (LHC) energies are predicted.The main dynamical input in our calculations is the use of solutions of the running coupling Balitsky-Kovchegov equation tested in electron-proton (ep) collision data. Particle production is computed via the hybrid formalisms to obtain spectra and yields.These baseline predictions are useful for testing the current understanding of the dynamics of very strong color fields against the upcoming LHC data.


Introduction
As a self-consistent effective perturbative quantum chromodynamics (QCD) theory at high energy (or small momentum fraction ), the color glass condensate (CGC) framework has been successfully applied to many QCD processes in highenergy collisions [1,2].Generally, in the initial state, the collisions involve at least one hadron or nucleus.Including the needed unitarity corrections, the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) renormalization group equations [3][4][5][6][7][8] capture the dynamical features of the CGC, which provides the tools to perform such resummation, although the precise prescription for the resummation may vary from process to process or colliding system.
In the framework of the CGC, as an instrumental tool, the dipole model formulation of deep inelastic scattering (DIS) allows for a relatively simple implementation of saturation effects in the description of the scattering process [9,10].In fact, in the small- range, many experimental data with a number of features are successfully described by the dipole model.By requiring higher-order corrections to the Balitsky-Kovchegov-JIMWLK (BK-JIMWLK) equations [3][4][5][6][7][8][11][12][13], an important progress was achieved in the past years.Recently, through the calculation of higher-order corrections to the BK-JIMWLK equations, an important step towards promoting the CGC to a practical phenomenological tool was performed [14][15][16][17].In particular, the calculation of running coupling (  ) corrections to the evolution kernel of the Balitsky-Kovchegov (BK) equation [18,19] renders the possibility to describe various data at high energies.In particular, the BK equation with running coupling (referred to as rcBK equation [18,19]) grasps most of the higher-order effects and incorporates nonlinear high gluon density effects.These processes close the gap between the first principle theory calculations and the experimental data.
In the CGC formalism, by using varying degrees of approximations and models, many researchers investigated single inclusive hadron productions in proton-nucleus () collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) [20][21][22][23][24][25][26][27][28][29][30][31].The CGC framework describes well single-particle distributions in both 2 Advances in High Energy Physics proton-proton () and deuton-gold (+Au) collisions, as well as the depletion of nuclear modification factors  Au [27,29].In deuton-nucleus () collisions at the RHIC, the observed suppression of single inclusive hadron production and the disappearance of the away-side peak in doublehadron production in the forward rapidity region are perhaps the strongest evidence for the importance and possibly dominance of saturation effects.The CGC calculations firstly predicted these effects [32] and later on quantitatively accounted for such suppression [33,34].This will be tested further at the LHC where, due to the larger energy of the collisions, the CGC dynamics in a much larger kinematic region can be probed.In the theoretical framework of the CGC, the single inclusive hadron production cross-section captures the saturation dynamics.The most important ingredient of the cross-section is the fundamental (or adjoint) dipole crosssection, that is, the imaginary part of the quark-antiquark () scattering amplitude on a proton or nucleus target.This dipole cross-section satisfies the BK-JIMWLK evolution equations [3][4][5][6][7][8][11][12][13].It also resums the small- and high gluon density effects.With next-to-leading-order (NLO) accuracy [14][15][16][17]35], one now knows the evolution equation for the dipole cross-section.
Because of the larger valence charge density in a nucleus, the CGC effect in a nucleus is expected to be more enhanced than that in a proton.However, due to a large amount of highquality DIS data on protons at small-, the most exhaustive searches of the saturation phenomenon are performed using data on proton reactions.The global fits performed by the Albacete-Armesto-Milhano-Quiroga-Salgado (AAMQS) group [36][37][38], in  collision at the HERA [39], show that the rcBK equation successfully accounts for the -dependence of inclusive structure functions at small-.We would like to point out that the AAMQS initial condition is not derived from an effective action (Gaussian action) for the large- valence charges.Based on a non-Gaussian action (quartic action), which involves the subleading correction of valence color charges, Dumitru and Petreska derived a new dipole scattering amplitude [40].Thus, we may compare the two different initial conditions (AAMQS and quartic action), which could illustrate the effect of valence color charges in  collisions at high energy.In particular, as the number of nucleons is finite, the process of scattering then involves the correction of valence color charges effects.
Due to the scientific significances and interests, in this paper, we study the inclusive particle productions in  and +Au collisions at the RHIC energy in the framework of the CGC theory by using the previously mentioned two initial conditions.Then, the inclusive particle productions in  collisions at the LHC energies are predicted by us.As we know, the situations at central rapidity and forward rapidity are different.To focus our attention on fragmentation region, where the leading nucleons and cold spectators are affected by each other, we study only the inclusive particle productions at forward rapidity in this paper.In Sections 2 and 3, we introduce and review the high-energy evolution equation and give single inclusive hadron production, respectively.The conclusions are given in Section 4. 2.1.The BK Equation with Running Coupling.In the CGC framework, one needs the BK-JIMWLK equations to describe the quantum evolution of hadron structure towards small-.The BK-JIMWLK equations are equivalent to an infinite set of coupled nonlinear equations for the different Wilson line correlators, where these correlators systematically incorporate small- gluon emission to all orders [3][4][5][6][7][8].In the large number of colors (large-  ) limit, the JIMWLK evolution equations are reduced to the BK equation [11][12][13], while the later one is a closed-form equation for the evolution of the dipole amplitude.

High-Energy Evolution Equation
According to [11][12][13][14]18], the rcBK equation can be given in the following form: where N  ,  0 , and ⃗  (= ⃗  1 + ⃗  2 ) denote the forward dipoletarget scattering amplitude in the fundamental () representation, the value of  where the evolution starts, and the transverse dipole size, respectively.In order to give some numerical solutions and related applications of (1), we can use a kernel validity to running coupling accuracy.According to Balitsky's prescription [14], the kernel  run in (1) can be presented as follows: For the running coupling in the above equation, we may use the scheme proposed in [37] at one-loop level.
In particular, in the large-  limit, one has the following relation between the adjoint () and fundamental dipoles: The amplitude N () incorporates all multiscattering between a projectile color dipole and the target and encodes the small- dynamics.In the CGC framework, it can be obtained from the solution of rcBK evolution with a given initial condition.

Initial Condition: Beyond the McLerran-Venugopalan
Model.The initial conditions required by the high-energy evolution equations are at a rapidity of  = ln( 0 /) = 0. Generally, one takes  0 = 10 −2 .The dynamical equations were derived by McLerran and Venugopalan [41,42] in the limit of an infinitely large nucleus.Let   ( ⊥ ,  − ) denote the classical color charge density per unit transverse area  ⊥ and longitudinal length  − , where kinetic terms  are neglected because the transverse momenta are assumed to be small.In the McLerran-Venugopalan (MV) model [41,42], the large- valence color charges are then described as recoilless sources on the light cone with   ( ⊥ ,  − ).As  → ∞, the variance of the valence color charge density ∼  1/3 grows large corresponding to a very high-density  of charges.
According to the central limit theorem, the fluctuations of the color charges are described by a Gaussian effective action.Generally, the high multiplicity  collision corresponds to an unusually high valence color charge density [43].Then, in this case, we can use the Gaussian effective action to describe it.However, for the real case with the mass number , the number of valence color charges is finite.Particularly in  and peripheral nucleus-nucleus () collisions, the considered number is very limited.Therefore, one needs to consider extensions of the MV model action to involve higher powers of the color charge density.
In fact, Jeon and Venugopalan [44,45] introduced an "odderon" operator contribution −        / 3 to the effective action with  3 ∼  3 (/ 2 ) 2 , where  denotes the radius of the nucleus.Dumitru et al. [46] introduced the quartic order contribution             / 4 with  4 ∼  4 (/ 2 ) 3 .Finally, to go beyond the limit of infinite valence charge density, the effective action is written as follows [46]: Based on the new effective action, including the first subleading correction in the semiclassical approximation, Dumitru and Petreska [40] introduced the dipole scattering amplitude N() on a dense target.This dipole scattering amplitude can be matched with a phenomenological proton fit by Albacete et al. [36,37] over a broad range of dipole size  and can provide a definite prediction for the -dependence for heavy-ion targets.According to [40], this dipole scattering amplitude has the following form: where   and Λ denote the initial saturation momentum scale and the infrared scale, respectively, and For the nucleus with a mass , we have As discussed in [40], the amplitude decreases sharply at very short distances.This phenomenon may not be reasonable.As we note, that the new term of the valence color charge has an extremely larger contribution than the other one at very small dipole size.This term corresponds to strong interaction caused by higher-order correlations of valence color charges.More carefully, we have N() ∼  2 in the limit of  → 0. Based on this behavior, the valence color charge term is assumed to be controlled by a weight function  2 () to go smoothly to the scattering amplitude.The first term is assumed to be controlled by a weight function  1 ().The new scattering amplitude is given by the following: Following this procedure, the unreasonable phenomena can be avoided.It should be noted that taking the weight function involved is just an experiential treatment.We also assume experientially that the weight functions have the following form: The parameters above are obtained by balanceing the two terms in the amplitude.That means the contribution of the new term is vanished while the first term is domainal at very short distance.Then, the scattering amplitude for a dipole in the adjoint representation with a given  can be obtained.We would like to point out that to regularize the behavior at large  in the calculation, we have replaced ln(1/Λ) by ln( + (1/Λ)) in (7) and assumed an exponentiation of the O( 2 ) expression.This does not affect the behavior at   ≤ 1.
The dipole scattering amplitude for a proton target was fitted by the AAMQS model [36][37][38] which gives the amplitude, in the initial condition for a small- evolution, to be where the infrared scale Λ is taken to be 0.241 GeV.We notice that the case with  = 1 corresponds to the MV model which was derived from an effective action (Gaussian action) for the large- valence charges.The free parameters are parameter  and the initial saturation scale   , which is probed by quarks.For the proton and nuclear targets, we use respectively, the notations  0 and  0 for clearness.From a global fit to proton structure functions in DIS in the small- region and single inclusive hadron data in  collisions at the RHIC and the LHC, the proton initial saturation scale gives  2 0 ≃ 0.168 GeV 2 and  ≃ 1.119.Simultaneously, this model provides a good description of charged hadron transverse momentum distributions in  collision at the center-ofmass energy per nucleon pair √  = 7 TeV [47,48].
Solving the rcBK equation related to the impact-parameter is still an open question.This renders that in (1) the impact parameter dependence of the collisions is ignored.In fact, this question may not be important for the minimumbias analysis because the initial saturation scale  0 of a nucleus can be considered as an averaged value and extracted from the minimum-bias data.In the minimum-bias collisions, one may assume that  2 0 =  1/3  2 0 , where the parameter  is fixed from a fit to data.In [49], it was shown that  ≈ 0.5 if the DIS data is used for the description of heavy nuclear targets.This result is consistent with the inclusive hadron production data in minimum-bias +Au collisions at the RHIC energy, which gives an initial saturation scale within  2 0 ≈ 3 ÷ 4 2 0 .The uncertainties associated to the variation of initial saturation scale in the rcBK evolution equation should be taken into account.
The dependence of dipole scattering amplitudes in the adjoint representation (N  = 2N  − N 2  ) on   ( 0 ) are shown in Figure 1 with three different models: MV (the solid curves), AAMQS with  = 1.119 (the dashed curves), and quartic action (the dashed-dotted curves).Figure 1(a) is for a proton assuming that  2  = 0.168 GeV 2 , Λ 2 = 0.058 GeV 2 , and  = 0.01. Figure 1(b) is for a nucleus with  = 200 (which is between Au and Pb),  2   ∼  1/3 , and   ∼  −2/3 (or  = 0.0003).One can see that both the dipole scattering amplitudes resulted from the quartic action and AAMQS models are beyond the MV model.The result of the quartic action is similar to the AAMQS model over a broad   ( 0 ) range.At a small size of dipole, the result of the quartic action is much smaller than that of AAMQS.Dumitru and Petreska [40] argued that the suppression of the bremsstrahlung tails by the subleading term possibly explains the difference between the AAMQS and MV models within relevant region 0.05 <   < 1.At very small , however, the subleading correction starts to dominate the leading one.This implies the breakdown of the simple perturbative expansion.
To understand more easily the above phenomenon, the scattering amplitude N(, ) for a dipole in coordinate space at different rapidities is given in Figure 2. The solid and dashed-dotted curves represent the results of the AAMQS with  = 1.119 and quartic action, respectively.From the right curve to the left one, the corresponding rapidities are  = 0, 5, and 10, respectively.One can see that, at  = 0, there is a very little difference between the results obtained using the two different initial conditions.At  = 5 and 10, the two types of results are almost the same.These phenomena are caused by the gluon distribution of the proton itself, which is independent of any particular reactions and renders nearly the same result.

Structure Function and Comparison with Experimental
Data.To check the validity of these calculations, let us calculate the proton structure function  2 at small-.At  ≪ 1, the inclusive structure function obtained from DIS process is given by [50]  where  em ,  2 , and   *  , denote the electromagnetic coupling, the transverse resolution scale, and the virtual photonproton ( * ) cross-section for transverse () and longitudinal () polarization states of the virtual photon, respectively.In the dipole model, the total  *  cross-section is given by , (,  2 )   (,   ) , (11) where ,   , Ψ, , and   (,   ) denote the fraction of longitudinal momentum of the virtual photon carried by the quark, the transverse separation between the quark and the antiquark, the light-cone wavefunction for the photon splitting into a  pair [51,52], the flavor of the quark, and the dipole-proton cross-section, respectively.By using the optical theorem,   (,   ) can be obtained from the elastic dipoleproton scattering amplitude A(,   , Δ), where N(,   ,   ),   , and  0 denote the imaginary part of the forward elastic dipole-proton scattering amplitude, the impact parameter of the dipole-target collisions, and the normalization constant obtained by the   integration, respectively.
To see the dependence of structure function  2 on , in Figure 3 we show the comparisons between the experimental data (squares) in  collision [53][54][55] and the results by solving rcBK equation with the two different initial conditions, AAMQS (the solid curves) and quartic action (the dasheddotted curves), by using ( 2) and (3).Figures 3(a for the four  2 values in Figure 4.One can see that the two different initial conditions describe well the  2 data, although their behaviors in Figure 1 on N  (,  0 ) at small   ( 0 ) are different.

Single Inclusive Hadron Production
3.1.Hybrid Formalisms.In asymmetric collisions such as  to hadron (ℎ) and others () at high energy, the scattering occurs between a dilute system of parton and a dense one of nucleus [56,57].The cross section for single inclusive hadron production at forward rapidity is given by [58] where the quantities , ,   , , and   denote the normalization constant, the splitting fraction for a given parton splitting into two others, the small transverse momentum which is strongly affected by target saturation effects, the pseudorapidity of the produced hadron, and the transverse momentum of the considered hadron, respectively; the functions   (,  the final hadron ℎ with a momentum fraction , the inelastic weight function, and the Altarelli and Parisi [59] splitting function, respectively.The last two functions can be found in equations (74-77) of [58] and equations (16.26-16.29) of [60] for more details.The longitudinal momentum fractions  1 and  2 are defined as follows: where   is the Feynman variable in another form.In the above calculation, the hadron masses are neglected since we are only interested in the light hadron production at high   , which results the rapidity to equal approximately to the pseudorapidity.As in the two-dimensional Fourier transform, the amplitude   (  ) in ( 13) can be given by The first integral in (13) denotes the elastic term which was derived in [23] and corresponds to the scattering of collinear partons from the projectile on the target.In the process corresponding to the elastic term, a quark from the projectile scatters on the target and radiates a gluon either before or after the scattering.The incoming quark, the outgoing quark, and the radiated gluon can all multiply scatters on the target.Because the incoming parton initially has zero transverse momentum, but picks up transverse momentum of order   after multiply scattering on the target, the elastic term is the most important one.The other part in (13) denotes the inelastic term which was derived in [58] and is important only when the produced hadron has a much higher transverse momentum than   .The inelastic term corresponds to high transverse momentum partons radiated from the incoming partons in the projectile wave function.The partons with high transverse momentum may scatter off the target with only a small momentum transfer, and finally they are fragmentated into high-  hadrons.
We show the dipole scattering amplitude  () (  , ) and the Fourier transform of 1 − N () (, ) (15), for a proton target and a nucleus target with  = 200 in Figures 4(a) and 4(c), as well as Figures 4(b) and 4(d), respectively.The solid and dashed-dotted curves represents the results of the AAMQS with  = 1.119 and quartic action, respectively.In Figures 4(a) and 4(c), we take  2 0 = 0.168 GeV 2 ; in Figures 4(b) and 4(d), we take  2 0 = 0.5 GeV 2 .From the top curve to the lowest one, the corresponding rapidities are  = 2, 4, and 8, respectively.One can see that in most of the cases the differences between the two types of results obtained using the two initial conditions are small.The value of  () (  , ) for a proton target is much larger than that of  () (  , )/ 1/3 for a nucleus with  = 200.

Numerical Results and Predictions.
To give the numerical results of the hybrid formula, we can use the CT10 NLO [61] parton distribution functions and the DSS NLO [56,62] fragmentation functions.The only input for the rcBK equation is the initial conditions for the evolution of the dipole amplitude.Our numerical results will be given by using the AAMQS and the quartic action, respectively.As already mentioned previously, the inelastic term in (13) involves an additional coupling factor   at NLO or named a running coupling scale which is considered as a free parameter in [57].In this work, we take   = 0.1 which is the same as that in [57].
Our aim is not to fit the experimental data but to use the CGC theoretical tools available in low- physics to highlight the uncertainties involved in making robust predictions for the upcoming  collisions at the LHC.In all the calculations, we take  = 1 in (13).However, the STAR data at very forward pseudorapidity ( = 4) shows that a -factor of ∼ 0.4 ÷ 0.6 may be needed (see [27] for more details).The necessity of such a small value of -factor may correspond to a large- phenomena, which is not included in the CGC at very forward (pseudo)rapidity region.
In Figure 5 we show the single inclusive hadron production yields (transverse momentum distributions) in √  = 0.2 TeV  collisions at different pseudorapidities using the rcBK dipole solution with two different initial conditions, AAMQS (the solid curves) and quartic action (the dasheddotted curves).The symbols (triangles, circles, and squares) represent the experimental data of the STAR [63] and BRAHMS collaborations [64] for different hadrons ( 0 [63], ℎ − [64], and ℎ − [64]) at different pseudorapidities (4.0, 3.2, and 2.2) with different magnifying multiples (1, 50, and 200) in the distributions as shown in the figures.As we have already pointed out previously, the value of  0 at  0 = 0.01 for proton target in the range 0.168 ÷ 0.336 GeV 2 gives a consistent description of the experimental data from the RHIC.We then choose  2 0 = 0.168 GeV 2 (Figure 5(a)) and  2 0 = 0.336 GeV 2 (Figure 5(b)) as input parameters for rcBK equation with the two different initial conditions.From Figure 5 we see that the results of the two initial conditions are similar to each other and approximately in agreement with the data.
In Figure 6, we compare our calculated results of   distributions for  0 and ℎ − produced in minimum-bias +Au collisions at √  = 0.2 TeV with the experimental data of the STAR [63] and BRAHMS collaborations [64] at the RHIC.Except the interacting system and the values of  2 0 , other illuminations are the same as those for Figure 5.According to [27,57], the value of saturation scale for gold within  2 0 = 0.5 ÷ 0.67 GeV 2 is consistent with the RHIC data.It is difficult to pin down the exact value of  0 for nuclei since the DIS data for nuclear targets are limited and have large experimental uncertainties.In this case, we choose  2 0 = 0.5 GeV 2 for Figure 6(a) and  2 0 = 0.67 GeV 2 for Figure 6(b) as input parameters for rcBK equation with the two different initial conditions.We see that at  2 0 = 0.5 GeV 2 the quartic action seems better and at  2 0 = 0.67 GeV 2 the AAMQS seems better when we compare them with the experimental data of the STAR [63] and BRAHMS collaborations [64].
As shown in Figures 5 and 6, the single inclusive hadron production yields at different forward pseudorapidities in  and minimum-bias +Au collisions at the RHIC energy are approximately described by the rcBK dipole solution with the two different initial conditions.We notice that the two different initial conditions are matched together at very forward pseudorapidity ( = 4).However, the discrepancy goes large when the pseudorapidity goes small, which would be caused by the dynamical evolution with increasing the nuclear mass number.
We now present our predictions for single inclusive hadron production at the LHC in terms of the nuclear modification factor   and hope that some of the theoretical uncertainties, such as the sensitivity to -factor, will be reduced.The nuclear modification factor   is defined as where  coll is the number of binary nucleon-nucleon collisions for a given centrality class in  collisions, and   → ℎ / 2    and   → ℎ / 2    denote the single inclusive hadron productions in  and  collisions, respectively.We take  coll = 6.5 and 7.4 for minimum-bias +Pb collisions at √  = 4.4 and 8.8 TeV, respectively, according to a standard Monte Carlo Glauber model [65].
In Figure 7 we show the nuclear modification factor  ℎ  for inclusive charged hadrons ℎ + + ℎ − production in minimum-bias +Pb ( = Pb in ) collisions at √  = 4.4 TeV (Figure 7 using the solution of the rcBK dipole evolution equation with two different initial conditions, AAMQS (the solid curves) and quartic action (the dashed-dotted curves).In both panels, we assume the initial saturation scales of proton and nucleus to be  2 0 = 0.168 GeV 2 and  2 0 = 0.5 GeV 2 , respectively.From Figure 7 one can see that  ℎ  decreases with increasing the pseudorapidity.This is explained from nonlinear evolution which suppresses additional gluon emissions in dense wave functions as compared to the dilute limit.The relationship between  ℎ  and √  is weak in the LHC energy region.We note also from Figure 7 that the quartic action as the solution of rcBK equation leads to a faster rise of  ℎ  with transverse momentum, and the corresponding behavior is obviously different from AAMQS.This is because the particle production at high   in  collisions is mainly dominated by a high valence color charge density in the dense target.
Because the DIS data for nuclear targets are limited and have large experimental uncertainties, we calculate the possible ranges of  ℎ  for forward rapidities in minimum-bias +Pb collisions at two energies using two initial saturation momentum scales with two initial conditions, AAMQS (the double solid curves) and quartic action (the double dasheddotted curves), in Figure 8 (√  = 4.4 TeV) and Figure 9 (√  = 8.8 TeV), respectively.The lower curves correspond to  2 0 = 0.5 GeV 2 and the upper curves correspond to  2 0 = 0.67 GeV 2 for nuclear target.We note that there are large uncertainties in  ℎ  at small pseudorapidity at the LHC due to the choice of the initial saturation scale for the rcBK evolution equation.In particular, if we choose the special non-Gaussian action (quartic action) which involves the subleading correction of the valence color charges, the discrepancy will be clearly larger compared to the other one (AAMQS).
It should be noted that the particle production crosssection given by ( 13) is right at forward rapidity.If we calculate the particle production at central rapidity, we need another formula, the   -factorization formalism [66][67][68].This will dilute our comparison on the two initial conditions.To give prominence to the comparison, we have used only the hybrid formalism (13) at forward rapidity in this paper.

Conclusions
In this work we have used quantitatively a new initial condition (the quartic action initial condition) as input for the rcBK equation [18,19] and presented a description of the hadron yields measured at three different pseudorapidities in  and +Au collisions at the RHIC by using the hybrid formalism which is proposed in [58].The results calculated by the quartic action and AAMQS models [36][37][38] have been compared to each other and some interesting results have been obtained.
At medium and large   ( 0 ), the quartic action gives a comparable result of the scattering amplitude for an adjoint dipole with the MV [41,42] and AAMQS models.At small   ( 0 ), the quartic action gives a much lower scattering amplitude which is beyond the MV and AAMQS models.The results on the proton structure function  2 (,  2 ) calculated by the quartic action and AAMQS models are very similar and describe approximately the HERA data [53][54][55].In a few cases, for example,  2 = 8.5 GeV 2 , the two models present almost the same results on  2 (,  2 ).
In the description of hadron transverse momentum distributions in  collisions at the RHIC energy, both quartic action and AAMQS models present almost the same results at forward pseudorapidity.At the closely mid-pseudorapidity, the results calculated by the two models are similar.In +Au collisions at the RHIC energy, both of the two models present almost the same results at forward pseudorapidity on the   distribution.At the closely mid-pseudorapidity, the two modeling results have an obvious difference in the high   region, which is caused by the dynamical evolution with increasing the nuclear mass number.
Based on the comparisons with the experimental data at the RHIC energy, the nuclear modification factor  ℎ  for charged hadrons produced in  collisions at the LHC energies is predicted by using the previously mentioned two models.The calculated results show that  ℎ  decreases with increasing the pseudorapidity, and a large uncertainty in  ℎ at small pseudorapidity is obtained.The quartic action as the solution of rcBK equation leads to a faster rise of  ℎ  with transverse momentum, and the quartic action shows a larger uncertainty than the AAMQS.These behaviors are related to the correlation and density fluctuation of the valence color charges.Our calculated results in the present work are consistent with the results given in recent works [25,69].

Figure 2 :
Figure 2: Dipole scattering amplitude N(, ) in coordinate space at different rapidities with two different initial conditions indicated in the figure.From the right curve to the left one, the corresponding rapidities are  = 0, 5, and 10, respectively.

Figure 3 :
Figure 3: Proton structure function  2 (,  2 ) versus  at different  2 .The squares represent the experimental data of the HERA Collaboration [53-55].The solid and dashed-dotted curves represent the results calculated by the AAMQS and quartic action initial conditions, respectively.

Figure 4 :
Figure 4: Dependence of the dipole scattering amplitude  () (  , ), the Fourier transform of 1 − N () (, ) (15), on   for (a) a proton target and (b) a nucleus target with  = 200 at different rapidities with two different initial conditions indicated in the figure.From the upper curve to the lowest one, the corresponding rapidities are  = 2, 4, and 8, respectively.

Figure 5 :
Figure 5: Neutral pion and negatively charged hadron yields in  collisions at √  = 0.2 TeV.The symbols represent the experimental data at different pseudorapidities measured by the STAR [56] and BRAHMS Collaborations [57].The solid and dashed-dotted curves represent the results calculated by the AAMQS and quartic action, respectively, with  2 0 = 0.168 GeV 2 (a) and 0.336 GeV 2 (b) in the calculations.

Figure 7 :
Figure 7: Nuclear modification factor  ℎ  versus   for inclusive charged hadrons ℎ + +ℎ − production in minimum-bias +Pb ( = Pb in ) collisions at the LHC energies, with the pseudorapidity being 4, 5, 6, and 7 from the upper curve to lower one.The solid and dashed-dotted curves represent the results calculated by the AAMQS and quartic action, respectively, with  2 0 = 0.168 GeV 2 for  and 0.5 GeV 2 for Pb in the calculations.(a) corresponds to √  = 4.4 TeV with  coll = 6.5, and (b) corresponds to √  = 8.8 TeV with  coll = 7.4.

Figure 8 :
Figure 8: Nuclear modification factor  ℎ  versus   for inclusive charged hadrons ℎ + + ℎ − production in minimum-bias +Pb collisions at √  = 4.4 TeV with different .The double solid curves and the double dashed-dotted curves are calculated, respectively, by the AAMQS and quartic action, which represents the uncertainty ranges in  ℎ  when  2 0 = 0.5 to 0.67 GeV 2 .

Figure 9 :
Figure 9: Same as Figure 8, but for the results at √  = 8.8 TeV.