Print version ISSN 0103-9733
Braz. J. Phys. vol.37 no.1a São Paulo Mar. 2007
J/y production in the saturation regime
F. O. Durães
Departamento de Física, Centro de Ciências e Humanidades, Universidade Presbiteriana Mackenzie, C.P. 01302-907 São Paulo, SP, Brazil
Instituto de Física, Universidade de São Paulo, C.P. 66318, 05389-970 São Paulo, SP, Brazil
In this work we calculate the J/y production in the initial stage of proton-proton, proton-nucleus and nucleus-nucleus collisions at RHIC and LHC energies, taking into account the high parton density regime of QCD, where the physics of parton saturation is expected to be dominant. We perform a quantitative analysis of the xF distributions in these collisions with the Color Glass Condensate (CGC) approach. The ratio between distributions with or without saturated gluons shows that this mechanism produces a suppression on the J/y yield in the forward region and presents a dip in intermediate values of xF, which is visible only at LHC energies.
Keywords: J/y production; Nucleus-nucleus collisions; Color Glass Condensate
Heavy quark production in hard collisions of hadrons has been considered  as a clean test in order to provide some of the most important backgrounds to new physics processes, which have been studied at DESY-HERA, Tevatron and LHC. The main motivation of these studies is the strong dependence of the total cross section on the behavior of the gluon distribution, which determines the QCD dynamics at high energies.
At small Bjorken-x and/or large A (in nuclei collisions) one expects the transition of the regime described by the linear dynamics (DGLAP, BFKL), where only the parton emissions are considered, to a new regime where the physical process of recombination of partons becomes important in the parton cascade and the evolution is given by a nonlinear evolution equation [2, 3].
In a scattering at high energies, the wave function of a hadron (or nucleus), boosted to large rapidity, exhibits a large number of gluons at small Bjorken-x and the density of gluons, per unit of transverse area and of rapidity, provides an intrinsic momentum scale which grows with atomic number (for nuclei) and with rapidity, due to continued gluon radiation as phase space grows. It is denoted by Qs, the so-called saturation momentum scale, which marks the onset of gluon recombination.
The large values of the gluon distribution at saturation (large occupation number of the soft gluon modes) suggests the use of semi-classical methods, which allow to describe the small-x gluons inside a fast moving nucleus by a classical color field. This color field is driven by a classical Yang-Mills equation whose source term is provided by faster partons. In this regime a Color Glass Condensate is expected to be formed [4-7], being characterized by a new type of high density matter, where the transition amplitudes are dominated not by quantum fluctuations, but by the configurations of classical field containing large ( » 1/as) numbers of gluons. Even though the coupling as becomes small due to the high density of partons, the fields interact strongly due to the classical coherence.
The CGC framework is an effective theory describing high energy scattering in QCD and have been applied to study phenomena at a number of existing and upcoming high energy collider facilities, in particular related to the initial conditions and equilibration in heavy ion collisions and to describe heavy ion phenomenology.
Recently, the J/y production in proton (deuteron)-nucleus collisions at high energies was considered in the CGC framework . The main conclusion of that work was that, at high energies (LHC energies), the saturation of gluons in the proton will result in the exact x1 scaling and therefore, because xF = x1 when x2 « 1, an exact xF scaling can be considered as a signature of the saturation in the proton. Apart from this interesting effect, other previous investigations concluded that charm production will not be very much affected by nonlinear effects because, in most of the interesting cases, the saturation scale is never much larger than the charm quark mass . However, even if the overall effect is small, saturation may affect strongly a certain small region of the phase space, or in this particular case the xF space. We found this situation before, studying the xF distribution of J/y's produced in the quark gluon plasma . Motivated by this suspicion and by our previous experience with xF spectra, in this work we calculate the xF distributions of J/y, produced in proton-proton and proton-nucleus collisions at RHIC and LHC energies. We shall use the Color Evaporation Model (CEM)  in two different regimes for the gluon densities: the regime described by the linear dynamics (DGLAP evolution) and the regime where gluon recombination is governed by nonlinear dynamics and leads to the saturation effect. Although our approach is phenomenological, as it will be seen, our results show a clear signal of saturation in the xF spectrum of J/y's produced in proton - nucleus collisions.
II. J/y PRODUCTION IN PROTON-PROTON AND PROTON-NUCLEUS COLLISIONS
In the CEM, charmonium is defined kinematically as a c state with mass below the D threshold. In leading order (LO) the cross section is computed with the use of perturbative QCD for the diagrams of the elementary processes ® c and gg ® c convoluted with the parton densities in the projectile (A) and in the target (B).
Calling xF the fractional momentum of the produced pair (with respect to the momentum of a projectile nucleon in cm frame) and the cm energy of a nucleon - nucleon collision, the cross section for production of a c pair with mass m is just given by:
where x1 and x2 are the nucleon momentum fractions carried respectively by partons in the projectile and target. The function H(x1,x2;m2), which represents the convolution of the elementary cross sections and parton densities is given by:
with the parton densities fi(x,m2) in the nucleon computed at the scale m2 = x1 x2 s.
The LO elementary cross sections in terms of the pair invariant mass (m) are given by [12, 13]:
where mc is the mass of the c quark. The production cross section of the charmed state i ( = J/y, y¢ or ccJ), s i, is then finally obtained by integrating the free pair cross section c over the invariant mass m starting from the production threshold 2 mc ( = 2.4 GeV) up to open charm production threshold 2 mD ( = 3.74 GeV). This model describes well the experimentally measured xF distribution of hidden charm both with LO and NLO cross sections, provided that is defined as multiplied by a theoretical factor k, which is equal to the ratio of the NLO and LO cross sections ( @ 2.54 %) . Therefore the J/y production cross section can be written as:
In the CEM non-perturbative effects are hidden in the factors k, and in the choice of the integration domain in the expression above. There are other ways to deal with these effects, either employing the Color Octect Model  or using QCD sum rules, where cross sections become functions of the vacuum condensates .
In what follows we shall use the CEM to study perturbative J/y production in proton-proton (A = B = 1), proton-nucleus (A = 1) and nucleus-nucleus collisions at RHIC and LHC. As it is well known, in nuclear collisions and in processes involving small values of x, shadowing plays an important role (see, for example, ). However, for simplicity we will, in a first moment, neglect nuclear effects and will make use of parton distribution functions (PDF) in the proton taken from  (GRV98 LO).
We shall also use this approach to study the J/y production in the CGC framework. In the quark/antiquark sectors we will employ the same PDF's above to explicitly investigate the effects of saturation of gluons on the J/y production. For the gluon densities in this regime we will adopt the Ansatz of Ref. [18, 19]:
where the saturation scale is given by:
with = 0.34 GeV2, x0 = 3.10-4 and l = 0.29.
The factor is introduced to account for the fact that the gluon density is small at x ® 1, according to the quark counting rules. We will make use of the dependence of the parameter C2 on Q2, which is typically parameterized in terms of the ln [ln (Q2/L2)/ ln (Q2/µ2)] ( C2 @ 4 at small Q2), taken from the Ref. .
The comparison between the two scenarios for J/y production, namely with and without gluon saturation, only makes sense if both gluon distributions in these regimes obey the same energy-momentum conservation relation in the proton. We will thus fix the normalization constant C1 in our modified gluon density in the saturation regime, at given Q2, requiring that:
Before presenting results a remark is in order. Expressions (1) and (2) imply the validity of collinear factorization, which is violated in many cases in the context of saturation physics. In many cases of interest, expressions analogous to (1) and (2) are valid, in which we have to replace G(x) by the unintegrated (in the gluon transverse momentum) gluon distribution f(x,). This is called kT factorization and was proven to hold in many cases. As shown in , kT factorization is valid for gluon production but is violated in quark production, especially in p-A and A-A collisions. However, from the quantitative point of view, this violation is not very large. Moreover, as shown in , kT and collinear factorization are equivalent at the leading twist level. Given the exploratory nature of this study, we shall assume that (1) and (2) hold also for distributions like (6).
III. RESULTS AND DISCUSSION
In order to investigate the differences between the two J/y production mechanisms we will first construct the ratio between their differential cross sections:
In Figure 1 we show these ratios of J/y momentum distributions in collisions pp,pAu and Au Au at = 0.2 TeV (RHIC) and pp, pPb and PbPb at = 5.5 TeV and = 14 TeV (LHC). It should be noted that the factors k and in eq. (5) are cancelled out in the ratio (9). The same occurs, for proton-nucleus collisions, with the factors A º p ( = 1) and B º Au(Pb) ( = 197 (208)) appearing in eq. (2) because we have neglected nuclear effects and assumed (x,Q2) = A (x,Q2) in both mechanisms for J/y production. In spite of this cancellation there is a nontrivial dependence on the atomic number A coming from the saturation scale (7).
In Figure 1, we observe a suppression on J/y production, which becomes stronger for larger systems and for higher energies. In these situations we access more and more the saturated gluons with smaller x1 and x2 (remember that 4 < m2 = x1x2s < 4 , see eq. (5)) which contribute with larger weight than at RHIC for J/y production. The results in Figure 1 might suggest that gluon saturation produces a significant suppression on J/y production. However this might be only an indication that the KLN and GRV98 PDF's are very different from each other. In order to see how much of this difference comes from the nonlinear regime we calculate the ratios:
where the cross section is always for the process A B ® J/y X and ''linear'' means that the first line of (6) is switched off and we are considering only the case Q2 > , where the evolution is in the linear regime. In this expression ''full'' means that the whole distribution (6) is active. If (10) is equal to one, then saturation plays no role. In Figure 2 we show the ratios (10) for pp (upper pannel) and p A collisions (lower pannel). We see that for pp reactions, even at higher energies the nonlinear effects are negligible. For pA and AA collisions they are negligible at RHIC but become rapidly very strong at LHC energies.
Although our calculation is admitedly crude, the results indicate very clearly that, at RHIC, saturation physics (the Color Glass Condensate) plays no role in J/y suppression. Therefore the J/y suppression at large rapidities measured by PHENIX  must come from some other mechanism. Of course, our conclusion has to be confirmed and the calculations redone with the more rigorous approach proposed in . If this conclusion still holds, then understanding J/y production at RHIC will become a challenge for theorists.
This work has been supported by CNPq and FAPESP. FOD was also supported by Instituto Presbiteriano Mackenzie, through MACKPESQUISA.
 K.J. Eskola, V.J. Kolhine, and R. Vogt, Phys. Lett. B 582, 157 (2004). [ Links ]
 For reviews on the subjects see: E. Iancu and R. Venugopalan, hep-ph/0303204; [ Links ]J. Jalilian-Marian and Y. V. Kovchegov, Prog. Part. Nucl. Phys. 56, 104 (2006) and references therein. [ Links ]
 A.M. Stasto, Acta Phys. Polon. B 35, 3069 (2004). [ Links ]
 L. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994). [ Links ]
 J. Jalilian-Marian, A. Kovner, L. McLerran, and H. Weigert, Phys. Rev. D 55, 5414 (1997). [ Links ]
 Yu. V. Kovchegov, Phys. Rev. D 60, 034008 (1999). [ Links ]
 E. Iancu, A. Leonidov, and L. McLerran, Nucl. Phys. A 692, 583 (2001). [ Links ]
 D. Kharzeev and K. Tuchin, Nucl. Phys. A 770, 40 (2006). [ Links ]
 V. P. Gonçalves and M. V. T. Machado, Eur. Phys. J. C 32, 501 (2004). [ Links ]
 V. Barger, W. Y. Keung, and R. N. Philips, Phys. Lett. B 91, 253 (1980). [ Links ]
 M. Glück and E. Reya, Phys. Lett. B 79, 453 (1978). [ Links ]
 R. V. Gavai et al., Int. J. Mod. Phys. A 10, 3043 (1995). [ Links ]
 C. Y. Wong, Phys. Rev. D 60, 114025 (1999). [ Links ]
 F. O. Durães, S. H. Lee, F. S. Navarra, and M. Nielsen, Phys. Lett. B 564, 97 (2003). [ Links ]
 N. Hammon, L. Gerland, H. Stocker, and W. Greiner, Phys. Rev. C 59, 2744 (1999). [ Links ]
 D. Kharzeev and E. Levin, Phys. Lett. B 523, 79 (2001). [ Links ]
 D. Kharzeev, E. Levin, and M. Nardi, Nucl. Phys. A 730, 448 (2004). [ Links ]
 J. Kwiecinski and A.D. Martin, Phys. Lett. B 353, 123 (1995). [ Links ]
 A. Adare (PHENIX Collaboration), arXiv:nucl-ex/0611020. [ Links ]
Received on 29 September, 2006