#### Abstract

We investigate charmonium production in Pb + Pb collisions at LHC beam energy TeV at fixed-target experiment ( = 72 GeV). In the frame of a transport approach including cold and hot nuclear matter effects on charmonium evolution, we focus on the antishadowing effect on the nuclear modification factors and for the yield and transverse momentum. The yield is more suppressed at less forward rapidity () than that at very forward rapidity () due to the shadowing and antishadowing in different rapidity bins.

#### 1. Introduction

Recently a fixed-target experiment using the LHC beams has been proposed [1], where the study on quarkonia in nuclear collisions becomes specifically important, due to the wide parton distributions in phase space which is helpful to reveal the charmonium production mechanism [2]. Corresponding to the LHC beam energy TeV, where is the nucleon number of the incident nucleus, the center-of-mass energy GeV is in between the SPS and RHIC energies, and a quark-gluon plasma is expected to be created in the early stage of heavy ion collisions. Taking into account the advantage of high luminosity in fixed-target experiments, which is helpful for detailed study of rare particles, the yield in Pb + Pb collisions at TeV per LHC run year is about 100 times larger than the yield in Au + Au collisions at GeV per RHIC run year [1]. With the high statistics, one may precisely distinguish between different cold and hot nuclear matter effects on charmonium production [3]. As is well known, the shadowing effect [4, 5], namely, the difference between the parton distributions in a nucleus and in a free nucleon, depends strongly on the parton momentum fraction . Since runs in a wide region, , in the fixed-target experiments, it provides a chance to see clearly the shadowing effect on the charmonium distributions in different rapidity bins. In this paper, we study the shadowing effect on the nuclear modification factors for yield and transverse momentum in Pb + Pb collisions at LHC beam energy TeV.

#### 2. Evolution of Quark-Gluon Plasma

The medium created in heavy ion collision at GeV is assumed to reach local equilibrium at a proper time fm/c [6]; its consequent space-time evolution can be controlled by the ideal hydrodynamic equations:where , , and , , , and are, respectively, the energy-momentum tensor, baryon current, four-velocity of the fluid cell, energy density, pressure, and baryon density of the system. The solution of the hydrodynamic equations provides the local temperature , baryon chemical potential , and fluid velocity of the medium which will be used in the calculation of the charmonium suppression and regeneration rates [7]. Taking the assumption of Hubble-like expansion and initial boost invariance along the colliding direction for high energy nuclear collisions, we can employ the well tested -dimensional version of the hydrodynamics in describing the evolution of the medium created at GeV. Introducing the proper time and space-time rapidity instead of the time and longitudinal coordinate , the conservation equations can be simplified as [8]with the definitions , **,** and , where and are the three-velocity of the fluid cell and Lorentz factor in the transverse plane.

To close the hydrodynamical equations one needs to know the equation of state of the medium. From recent studies on particle elliptic flow and shear viscosity, the matter created in heavy ion collisions at RHIC and LHC energies is very close to a perfect fluid [9]. Considering that the momentum integrated particle yield, especially for heavy quarkonia, is not sensitive to the equation of state, we follow [10] where the deconfined phase at high temperature is an ideal gas of gluons and massless and quarks plus 150 MeV massed quarks, and the hadron phase at low temperature is an ideal gas of all known hadrons and resonances with mass up to 2 GeV [11]. There is a first-order phase transition between these two phases. In the mixed phase, the Maxwell construction is used. The mean field repulsion parameter and the bag constant are chosen as MeV fm^{3} and MeV to obtain the critical temperature MeV [10] at vanishing baryon number density. Note that when one calculates the rapidity or transverse momentum distribution of quarkonia, the choice of the equation of state may result in sizeable difference.

The initialization of the hot medium is taken as the same treatment in [8]. We use the final charged multiplicity to determine the initial entropy density. For GeV, the charged multiplicity at central rapidity in center-of-mass frame is estimated to be based on the empirical formula [12]: The initial baryon density is obtained by adjusting the entropy per baryon to be 250 [13]. From the empirical relation [14] between the inelastic nucleon-nucleon cross section in unit of mb and the colliding energy in unit of GeV, we have mb at GeV. These initial conditions lead to a maximum medium temperature MeV at the initial time fm/c. The medium maintains local chemical and thermal equilibrium during the evolution. If we do not consider the charmonium interaction with the hadron gas, the charmonium distributions in the final state will be fixed at time corresponding to the critical temperature of the deconfinement phase transition.

#### 3. Charmonium Transport in Quark-Gluon Plasma

Since a charmonium is so heavy, its equilibrium with the medium can hardly be reached; we use a Boltzmann transport equation to describe its phase space distribution function in heavy ion collisions at impact parameter : where the loss and gain terms and come from the charmonium dissociation and regeneration in the created hot medium. We have neglected here the elastic scattering, since the charmonium mass is much larger than the typical medium temperature. Considering that the feed-down from the excited states and to the ground state [15] happens after the medium evolution, we should take transport equations for and when we calculate the distribution in the final state.

Introducing the momentum rapidity and transverse energy to replace the longitudinal momentum and energy , the transport equation can be rewritten as with the dissociation and regeneration rates and , where the third term in the square bracket arises from the free streaming of with transverse velocity which leads to a strong leakage effect at SPS energy [16].

Considering the gluon dissociation in the quark-gluon plasma, the dissociation rate can be expressed as where is the gluon energy, the flux factor, and the gluon thermal distribution as a function of the local temperature and fluid velocity determined by the hydrodynamics. The dissociation cross section in vacuum can be derived through the operator production expansion (OPE) method with a perturbative Coulomb wave function [17–21]. However, the method is no longer valid for loosely bound states at high temperature. To reasonably describe the temperature dependence of the cross section, we take the geometric relation between the averaged charmonium size and the cross section: The averaged radial square is calculated via potential model [22] with lattice simulated heavy quark potential [23] at finite temperature. When approaches the charmonium dissociation temperature , the averaged radius square and in turn the cross section go to infinity, which means a complete charmonium melting induced by color screening [24]. Using the internal energy as the heavy quark potential , the dissociation temperature is calculated to be , , and for , , and , respectively [22].

The regeneration rate is connected to the dissociation rate via the detailed balance between the gluon dissociation process and its inverse process [25, 26]. To obtain the regeneration rate, we also need the charm quark distribution function in medium. Although the initially produced charm quarks would carry high transverse momentum, they lose energy (momentum) when passing through the medium. Considering the experimentally observed large open charm quench factor [27–29] and elliptic flow [30, 31], we take as a first approximation a kinetically thermalized momentum spectrum for the charm quark distribution . Neglecting the creation and annihilation of charm-anticharm pairs inside the medium, the spatial density of charm quark number satisfies the conservation law with the initial density determined by the nuclear geometry , where are the thickness functions and is the charm quark rapidity distribution in p + p collisions.

For the regeneration rate , we also considered the canonical effect which is shown to be important in explaining the suppression of strange mesons [32]. When there are only few pairs or even less than one pair of charm quarks produced in an event, one needs to consider the canonical effect to guarantee the exact charm number conservation. Taking into account the fact that the charm and anticharm quarks inside a pair are produced at the same rapidity, we simply multiply the regeneration rate in a unit rapidity bin by a canonical enhancement factor [33]

To take into account the relativistic effect on the dissociation cross section to avoid the divergence in the regeneration cross section, we should replace the charmonium binding energy by the gluon threshold energy in the calculations of and [34].

In the hadron phase of the fireball with temperature , there are many effective models that can be used to calculate the inelastic cross sections between charmonia and hadrons [35]. For the dissociation cross section is about a few mb which is comparable with the gluon dissociation cross section. However, considering that the hadron phase appears in the later evolution of the fireball, the ingredient density of the system is much more dilute in comparison with the early hot and dense period [7]. Taking, for instance, the regeneration processes in quark matter and in hadron matter, the density ratio between charm quarks at initial temperature MeV and mesons at critical temperature MeV is around . Considering further the life time of the quark matter ~6 fm/c and the life time of the hadron matter ~2 fm/c calculated from the hydrodynamics in Section 2, we neglect the charmonium production and suppression in hadron gas, to simplify the numerical calculations. Note that the suppression and regeneration in hadron gas may become important for excited charmonium states [36].

The transport equation can be solved analytically with the explicit solution [7, 37]withThe first and second terms on the right-hand side of solution (10) indicate the contributions from the initial production and continuous regeneration, respectively, and both suffer from the gluon dissociation in the medium. Since the regeneration happens in the deconfined phase, the regenerated quarkonia would have probability to be dissociated again by the surrounding gluons. The coordinate shifts and in solution (10) reflect the leakage effect in the transverse and longitudinal directions.

For fixed-target nuclear collisions at TeV, the collision time for the two Pb nuclei to pass through each other in the center of mass frame is fm/c, which is compatible with the charmonium formation time but shorter than the QGP formation time fm. Therefore, all the cold nuclear matter effects can be reflected in the initial charmonium distribution at time . We take into account nuclear absorption, nuclear shadowing, and Cronin effect. The initial distribution in solution (10) can be obtained from a superposition of p + p collisions, along with the modifications from these cold nuclear matter effects.

The nuclear absorption is important in explaining the suppression in p + A and A + A collisions at low energies. It is due to the inelastic collision between the initially produced charmonia and the surrounding nucleons, and its effect on the charmonium surviving probability can be described by an effective absorption cross section . The value of is usually measured in p + A collisions and is several mb at SPS energy. Since the nuclear absorption becomes weaker at higher colliding energy due to the shorter collision time [38, 39], we take mb at TeV [39] and the nuclear absorption factor

The Cronin effect broadens the momentum distribution of the initially produced charmonia in heavy ion collisions [7]. In p + A and A + A collisions, the incoming partons (both gluons and quarks) experience multiple scatterings with surrounding nucleons via soft gluon exchanges. The initial scatterings lead to an additional transverse momentum broadening of partons which is then inherited by produced hadrons [40]. Since the Cronin effect is caused by soft interactions, rigorous calculations for the effect are not available. However, the effect is often treated as a random motion. Inspired from a random-walk picture, we take a Gaussian smearing [41, 42] for the modified transverse momentum distribution: whereis the path length of the two initial gluons in nuclei before fusing into a charmonium at , and , is the averaged charmonium transverse momentum square gained from the gluon scattering with a unit of length of nucleons, and is the charmonium momentum distribution in a free p + p collision. The Cronin parameter is usually extracted from corresponding p + A collisions. Considering the absence of p + A collision data at GeV, we take /fm from some empirical estimations [4, 25, 43]. As a comparison, for collisions at SPS (GeV) and RHIC (GeV) we take [8] and 0.1 [44] /fm, respectively.

Assuming that the emitted gluon in the gluon fusion process is soft in comparison with the initial gluons and the produced charmonium and can be neglected in kinematics, the charmonium production becomes a process approximately, and the longitudinal momentum fractions of the two initial gluons are calculated from the momentum conservation: The free distribution can be obtained by integrating the elementary partonic process: where is the gluon distribution in a free proton, the emitted gluon rapidity, the charmonium momentum distribution produced from a gluon fusion process, and the factorization scale of the fusion process.

Now we consider the shadowing effect. The distribution function for parton in a nucleus differs from a superposition of the distribution in a free nucleon. The nuclear shadowing can be described by the modification factor . To account for the spatial dependence of the shadowing in a finite nucleus, one assumes that the inhomogeneous shadowing is proportional to the parton path length through the nucleus [45], which amounts to considering the coherent interaction of the incident parton with all the target partons along its path length. Therefore, we replace the homogeneous modification factor by an inhomogeneous one [46] with the definition . We employ in the following the EKS98 package [5] to evaluate the homogeneous ratio , and the factorization scale is taken as .

Replacing the free distribution in (16) by the modified distribution and then taking into account the Cronin effect (13), we finally get the initial charmonium distribution for solution (10): Now the only thing left is the distribution in a free p + p collision which can be fixed by experimental data or some model simulations.

#### 4. Numerical Results

The beam energy TeV in fixed-target experiments corresponds to a colliding energy GeV, and the rapidity in the center-of-mass frame is boosted in the laboratory frame with a rapidity shift . Let us first focus on the central rapidity region around in the center-of-mass frame, which corresponds to in the laboratory frame. The centrality and momentum dependent antishadowing for initially produced charmonia is reflected in the inhomogeneous modification factor for gluons. The longitudinal momentum fractions are for the two gluons, which are located at the strong antishadowing region [47] by some parametrization of parton distribution shadowing like EKS98 [5], EPS08 [48], and EPS09 [49]. The antishadowing changes not only the gluon distribution but also the charm quark production cross section used in the regeneration. For the process , the antishadowing for gluons leads to an antishadowing factor, for the cross section. Considering that in peripheral collisions the regeneration is weak and its contribution is not remarkably affected by the antishadowing, we take a centrality averaged antishadowing factor for the cross section to simplify the numerical calculation for regeneration. Estimated from the EKS98 evolution [5], we take a enhancement of the charm quark production cross section compared to free p + p collisions. From FONLL calculation [50], the upper limit for is 0.047 mb at GeV. Note that the experimental data for charm quark cross section in free p + p collisions are close to the upper limit of perturbative calculation; we take mb at GeV. After taking into account the antishadowing effect in A + A collisions, it becomes 0.06 mb. For p + p collisions, we assume a constant hidden to open charm ratio at any colliding energy. From the ratio extracted from the RHIC data [51], we have *μ*b at GeV. The transverse momentum distribution for in free p + p collisions can be simulated by PYTHIA [52] and the mean transverse momentum square is .

Figure 1 shows our calculated centrality dependence of nuclear modification factor in Pb + Pb collisions at LHC beam energy TeV in laboratory frame ( GeV in center-of-mass frame) at forward rapidity (central rapidity ), where and are charmonium yields in p + p and A + A collisions, and and are numbers of binary collisions and participants. For comparison, we show also the RHIC data at GeV [53] at central rapidity. Since the shadowing/antishadowing effect is still an open question and its degree depends strongly on the models we used, we show in Figure 1 two calculations for the total in Pb + Pb collisions at GeV: one is with the above discussed antishadowing and the other is without antishadowing. The hatched band is due to this uncertainty in the antishadowing. With increasing collision centrality, the initial contribution drops down, while the regeneration goes up. The canonical effect is important in peripheral collisions where the number of charm quark pairs is less than one and the inclusion of the canonical effect enhances sizeably the charmonium yield. In most central collisions, the regeneration can contribute about to the total charmonium yield. The antishadowing at very forward rapidity in the laboratory frame (central rapidity in the center-of-mass frame) enhances the charm quark cross section and in turn the initial charmonium yield by a factor of 1.2. As a consequence, the enhancement factor for the regenerated charmonium number is which leads to a strong charmonium enhancement! If we do not consider the antishadowing effect on the charmonium regeneration and initial production, the total is significantly reduced.

To see more clearly the charmonium production mechanism, we turn to the transverse momentum information. In Figure 2 we show the nuclear modification factor [54] in Pb + Pb collisions at beam energy TeV, where and are averaged transverse momentum square in Pb + Pb and p + p collisions at very forward rapidity . If we neglect the contribution from the regeneration and consider only the initial production, the ratio goes up monotonously with centrality due to the Cronin effect and leakage effect [54]. The inclusion of regeneration (upper border of the band) remarkably reduces the averaged transverse momentum, because the regenerated charmonia possess a soft momentum distribution induced by the charm quark energy loss. Since the degree of regeneration increases with centrality, the increased soft component leads to a decreasing in most central collisions. The canonical effect can reduce the further, since it enhances the regeneration especially in peripheral collisions. However, we should note that the assumption of charm quark thermalization indicates a full energy loss and it may not be reached in peripheral and semicentral collisions at beam energy TeV. When we switch off the antishadowing (lower border of the band), both the hard component controlled by the initial production and the soft component dominated by the regeneration would be reduced. Considering that the enhancement factor resulted from the antishadowing is for the initial production but for the regeneration, the strong antishadowing in the soft component leads to only a slight difference between with and without considering the antishadowing, shown in Figure 2. It is obvious that, compared to the nuclear modification factor for the yield, the modification factor for the transverse momentum is less sensitive to the shadowing effect [54].

From the simulations of parton distributions in cold nuclear matter [5, 48, 49], the nuclear shadowing region is located at very small . In the following we consider the shadowing and see its difference from the antishadowing in and in fixed-target Pb + Pb collisions. The maximum rapidity in the center-of-mass frame is at GeV. Considering the expected amount of measured events, we focus on the backward rapidity region around which corresponds to the less forward rapidity in laboratory frame. From the kinematics, the momentum fractions for the two gluons involved in the gluon fusion process are and . One is located in the EMC region and the other in the shadowing region [5, 48, 49], leading to a reduction of for the charm quark production cross section from EKS98 evolution [5] ( from EPS09 NLO evolution [49]). Taking the same ratio of charm quark cross section between and calculated from FONLL [50] and including the shadowing reduction, we obtain mb at . For the medium evolution at this backward rapidity region, we initialize the entropy density to be half of that at central rapidity [6, 55] which leads to a maximum temperature of MeV. Figure 3 shows the two double ratios and of ; the upper and lower borders of the two bands correspond to the calculations with and without considering the nuclear shadowing and antishadowing. While the double ratio for the transverse momentum is not sensitive to the shadowing and antishadowing, as we discussed above, the strong antishadowing at and shadowing at lead to a strong enhancement of the double ratio for the yield. Without considering the shadowing and antishadowing, the stronger charmonium suppression in the hotter medium at ( MeV) compared with the weaker suppression in the relatively colder medium at ( MeV) makes the double ratio less than unit. However, the inclusion of the yield enhancement due to the antishadowing at and the yield suppression due to the shadowing at changes significantly the behavior of the double rati; it becomes larger than unit and can reach 1.3 in most central collisions. Note that the rapidity dependent shadowing effect was used to qualitatively interpret the stronger suppression at forward rapidity than that at midrapidity in Au + Au collisions at RHIC [56, 57].

#### 5. Summary

We investigated with a transport approach the charmonium production in fixed-target Pb + Pb collisions at LHC beam energy TeV. We focused on the rapidity dependent shadowing effect on the nuclear modification factors for the charmonium yield and transverse momentum. While the averaged transverse momentum is not sensitive to the shadowing effect, the antishadowing leads to a strong yield enhancement at very forward rapidity , and the shadowing results in a strong yield suppression at less forward rapidity . The double ratio between the nuclear modification factors in the two rapidity regions amplifies the shadowing effect; it is larger than unit and can reach 1.3 in most central collisions.

From the model studies on gluon distribution in nuclei, see, for instance [5, 47–49], there are large uncertainties in the domain of large , which is probably due to the unknown EMC effect. From our calculation here, the double ratio of the nuclear modification factor for yield is very sensitive to the gluon shadowing effect in different region. A precise measurement of the ratio may provide a sensitive probe to the gluon distribution.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The work is supported by the NSFC under Grant no. 11335005 and the MOST under Grant nos. 2013CB922000 and 2014CB845400.