Print version ISSN 0103-9733
Braz. J. Phys. vol.37 no.1a São Paulo Mar. 2007
Small-x physics at coherent pA/AA interactions at LHC
Victor P. GonçalvesI; Magno V. T. MachadoII
IInstituto de Física e Matemática, Universidade Federal de Pelotas Caixa Postal 354, CEP 96010-090, Pelotas, RS, Brazil
IICentro de Ciências Exatas e Tecnológicas, Universidade Federal do Pampa - UNIPAMPA Campus de Bagé. Rua Carlos Barbosa. CEP 96400-970. Bagé, RS, Brazil
We report on our recent investigations in photonuclear production of heavy quarks and vector mesons in ultraperipheral collisions at nucleus-nucleus and proton-nucleus reactions. They are initiated by quasi-real photons coming from one of the nuclei or hadron in interactions taking place at large impact parameter. We focus on the role played by the high energy theoretical/phenomenological approaches in photon-nuclei scattering, namely the strength of parton saturation phenomenon and high energy nuclear shadowing. In particular, our theoretical predictions are compared with the recent experimental measurements on coherent r (STAR) and J/Y (PHENIX) photoproduction at RHIC and estimates for LHC are given.
Keywords: Heavy Ion Collisions; Quantum Chromodynamics; Parton Saturation Effects
In the last years, there has been a lot of interest in the description of electron-nucleus collisions at high energies. The results of current analysis show that future electron-nucleus colliders at HERA and RHIC, probably could determine whether parton distributions saturate and constrain the behavior of the nuclear gluon distribution in the full kinematical range [1, 2]. However, until these colliders become reality we need to consider alternative searches in the current and/or scheduled accelerators which allow us to constrain the QCD dynamics. In this contribution, we summarize our investigations in Refs. [3-5] about the possibility of using ultraperipheral heavy ion collisions (UPC's) as a photonuclear collider and study the heavy quark and/or vector meson production assuming distinct approaches for the QCD evolution.
A special class of events in nucleus-nucleus (or proton-nucleus) reactions at high energies is their coherent interaction. The scatterings take place at larger impact parameter of the collision, b » RA + RB, where Ri are the radii of the interacting nuclei (or charged hadrons), i = A, B. This is the reason they are called ultraperipheral collisions. The final state is characterized by low particle multiplicity and dominance of low-pT transverse momenta. The interaction between the nuclei is predominantly electromagnetic in this case, where they have an associated flux of equivalent (quasi-real) photons. This flux is computed through the quantity called number of equivalent photons, dN/dw, where w is the center of mass photon energy. Then, real-photons from one of nuclei initiate photonuclar reactions. In heavy ion collisions, the large number of photons coming from colliding nuclei will allow to study photoproduction, with energies WgA reaching about 1 TeV at LHC. Therefore, while studies of photoproduction at HERA are limited to photon-proton center of mass energies of about 200 GeV, photon-proton interactions at LHC can reach one order of magnitude higher on energy. Consequently, studies of gp(A) interactions at LHC could provide valuable information on the QCD dynamics at high energies. The photonuclear cross section for coherent nucleus-nucleus collisions is given by the convolution between the photon flux from one of the nuclei and the cross section for the photon-nucleus scattering. The expression for the production of a given final state Y (Y = Q, V) ultraperipheral heavy ion collisions is then given by,
where wmin = /4gL mp, is the c.m.s energy of the nucleus-nucleus system and gL is the Lorentz boost of a single beam. The Lorentz factor for LHC is gL = 2930, giving the maximum c.m.s. gA energy Wg A < 950 GeV. For the case of coherent proton-nucleus collision, where the real-photons from the nucleus interacts with a proton, we now have,
The requirement that photoproduction is not accompanied by hadronic interaction (ultraperipheral collision) can be done by restricting the impact parameter b to be larger than twice the nuclear radius, RA = 1.2 A1/3 fm, in the AA case and than (Rp + RA) in the pA case. An analytic approximation for AA collisions can be obtained using as integration limit b > 2 RA, producing
where = 2wRA/gL for nucleus-nucleus and = w (Rp + RA)/gL for proton-nucleus. The quantities K0,1(x) are the modified Bessel functions. Therefore, the main ingredient for computing the production of the specific final state Y at UPC's is the information about its cross section in photon-nuclei interactions, sgA ® YX, as defined in Eq. 1 or in photon-proton interactions, sgp ® YX, as defined in Eq. 2. In what follows, we perform phenomenology on heavy quarks and vector mesons production using the available high energy approaches. The emphasis is placed in the energies and experimental cuts for LHC.
II. HEAVY QUARK PRODUCTION
For our further analysis on photonuclear production of heavy quarks we will consider distinct high energy approaches, which are contrasted against the collinear approach : (a) semihard formalism, (b) saturation model (within the color dipole formalism) and (c) Color Glass Condensate (CGC) approach. The main goal is a comparison among those approaches and show whether the production at UPC's allows disentangle the underlying QCD dynamics at high energies.
(a) Semihard formalism: In the k ^ factorization (semihard) approach, the relevant QCD diagrams are considered with the virtualities and polarizations of the initial partons, carrying information on their transverse momenta. The scattering processes are described through the convolution of off-shell matrix elements with the unintegrated parton distribution, F(x,k^). Considering only the direct component of the photon, the cross section reads as (See, e.g. Ref. ),
where D1 º ^ + and D2 º (k^-p1^)2 + and A(z) = [z2+ (1-z)2] . The transverse momenta of the heavy quark (antiquark) are denoted by p1^ and p2^ = (k^-p1^), respectively. The heavy quark longitudinal momentum fraction is labeled by z. For the scale µ in the strong coupling constant we use the prescription µ2 = + . We use also the simple ansatz for the unintegrated gluon distributions, Fnuc = where xGA(x,Q2) is the nuclear gluon distribution (see Ref.  for details in the numerical calculations).
(b) Saturation model: Based on the color dipole formalism, it can be extended to eA through Glauber-Gribov formalism. In this model the cross section for the heavy quark photoproduction on nuclei targets is given by [7, 8],
where the transverse wave function is known [8, 9]. The nuclear dipole cross section is given by [7, 8],
where b is the impact parameter of the center of the dipole relative to the center of the nucleus and the integrand gives the total dipole-nucleus cross section for fixed impact parameter. The nuclear profile function is labeled by TA(b). The parameterization for the dipole cross section takes the eikonal-like form, (, r2) = s0 [ 1- exp(- () r2/4 ) ], where one has used the parameters from saturation model, which include the charm quark with mass mc = 1.5 GeV and the definition = (Q2+ 4 )/. The saturation scale () = (x0/x)l GeV2, gives the onset of the saturation phenomenon to the process. The equation above sums up all the multiple elastic rescattering diagrams of the pair and is justified for large coherence length, where the transverse separation r of partons in the multiparton Fock state of the photon becomes as good a conserved quantity as the angular momentum, namely the size of the pair r becomes eigenvalue of the scattering matrix.
(c) Color Glass Condensate: The regime of a CGC is characterized by the limitation on the maximum phase-space parton density that can be reached in the hadron/nuclear wavefunction (parton saturation) and very high values of the QCD field strength Fµn » 1/as (For a review see, e.g, ). The large values of the gluon distribution at saturation 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. For our phenomenological purpose here, one takes a parameterization for the dipole cross section which takes into account the main features of the CGC theory. It reads as ,
where = rQsat(x) and the expression for > 2 (saturation region) has the correct functional form, as obtained from the theory of the Color Glass Condensate (CGC) . For the color transparency region near saturation border ( < 2), the behavior is driven by the effective anomalous dimension geff (x, r) = gsat + , where gsat = 0.63 is the LO BFKL anomalous dimension at saturation limit.
Having summarized the theoretical models, let us present the numerical calculation of their total cross section at UPC's in nucleus-nucleus collisions. We focus mostly on LHC domain where small values of x would be probed. In the following, one considers the charm and bottom masses mc = 1.5 GeV and mb = 4.5 GeV , respectively. Moreover, for PbPb collisions at LHC, one has the c.m.s. energy of the ion-ion system = 5500 GeV. The results are presented in Table I. The collinear approach gives a larger rate, followed by the semihard approach (using GRV98 gluon pdf's). The saturation model and CGC formalisms give similar results, including a closer ratio for charm to bottom production. Concerning the CGC approach, our phenomenological educated guess for the color field correlator seems to produce quite reliable estimates. Therefore, the photonuclear production of heavy quarks allow us to constrain already in the current nuclear accelerators the QCD dynamics since the main features from photon-nuclei collisions hold in the UPC reactions. It should be noticed that we have computed the integrated cross section by integrating all rapidities. Experimental cuts should diminish their magnitude, but they are still large. The experimental signal should be clear. Since photon emission is coherent over the entire nucleus and the photon is colorless we expect that the events to be characterized by one rapidity gap in heavy quark production. Finally, the CGC and color dipole calculations are more stable concerning changes in the quark mass (a few percent difference). Moreover, the value of central value for quark mass is determined by a fit to accelerator data rather than being a free parameter.
In Table I we also present the results for coherent proton-nucleus collisions at LHC. For comparison, we present the predictions from the linear dynamics (collinear), which is calculated assuming that the collinear factorization is valid and that the gluon distribution can be described by the GRV98 parameterization. At LHC we can observe a large difference between the predictions and having high rates. For instance, for charm quark CGC gives a cross section a factor 3 lower than collinear and a factor 2 for bottom. This deviation holds even in case of experimental cuts on rapidity. This has been verified in the case of vector meson production when considering experimental cuts . Therefore, photoproduction of heavy quarks should provide a feasible and clear measurement of the underlying QCD dynamics at high energies. Since RHIC is obtaining data for d Au interactions and LHC is to be commissioned, these processes could be analyzed in the next years. The advantages are a clear final state (rapidity gap and low momenta particles) and no competing effect of dense nuclear environment if compared with hadroproduction.
III. VECTOR MESON PRODUCTION
Let us consider the scattering process gA ® V A in the QCD dipole approach, where V stands for both light and heavy mesons. The main motivation for using this approach is the easy intuitive interpretation and its successful data description in the proton case. The scattering process can be seen in the target rest frame as a succession in time of three factorizable subprocesses: i) the photon fluctuates in a quark-antiquark pair (the dipole), ii) this color dipole interacts with the target and, iii) the pair converts into vector meson final state. In the color dipole formalism, the corresponding imaginary part of the amplitude at zero momentum transfer reads as [13, 14],
where (z, r) and (z, r) are the light-cone wavefunctions of the photon and vector meson, respectively. The quark and antiquark helicities are labeled by h and and reference to the meson and photon helicities are implicitly understood. The variable r defines the relative transverse separation of the pair (dipole) and z (1-z) is the longitudinal momentum fractions of the quark (antiquark).
In the dipole formalism, the light-cone wavefunctions (z, r) in the mixed representation (z,r) can be completely determined using light cone perturbation theory. On the other hand, for vector mesons, the light-cone wavefunctions are not known in a systematic way and they are thus obtained through models (For a recent detailed discussion see Ref. ). Here we follows the analytically simple DGKP approach , which assumes that the dependencies on r and z of the wavefunction are factorised, with a Gaussian dependence on r. We keep the original parameters of the model. The main shortcoming of this approach is that it breaks the rotational invariance between transverse and longitudinally polarized vector mesons . However, as it describes reasonably the HERA data for vector meson production, as pointed out in Ref. , we will use it here.
The corresponding parameters for the meson wavefunctions are presented in Table I of Ref. . Following Ref. , we have estimated contribution from real part for the photoproduction of vector mesons: it is about 3% for light mesons and it reaches 13% for J/Y . Additionally for heavy mesons we have taken into account the skewness effects, associated to off-forward features of the process, which are increasingly important in this case. Finally, the photonuclear cross section is given by
with tmin = (/2w)2. The quantity b is the ratio between the imaginary and real part of the amplitude. We have used an analytical approximation of the Woods-Saxon distribution as a hard sphere, with radius RA, convoluted with a Yukawa potential with range a = 0.7 fm in order to compute the nuclear form factor, F(t).
Recently, the STAR Collaboration at RHIC published the first experimental measurement of coherent r production in gold-gold UPC's at = 130 GeV . The energy dependence of the cross section is presented in Fig 1-a. Our theoretical prediction in the curve takes into account the experimental cuts, which gives stheory(|y| < 3) = 410 mb, in good agreement with the STAR measurement sSTAR(|y| < 3) = 370±170 (stat) ±80 (syst) mb. Furthermore, the PHENIX Collaboration has measurements of the cross section for coherent J/Y photoproduction in UPC at midrapidity at = 200 GeV , accompanied by nuclear breakup. The rapidity distribution is shown in Fig 1-b in comparison with our prediction. The theoretical estimation gives ds/dy|y = 0 = 58 µb, which is in agreement with the PHENIX measurement dsPHENIX/dy|y = 0 = 48±16 (stat) ±18 (syst) µb. Finally, in Table II one shows the predictions for the integrated cross sections for coherent nucleus-nucleus interactions at LHC. The experimental feasibility and signal separation on the reaction channels presented here are reasonably clear, namely applying a (low) transverse momentum cut pT < 1 GeV and two rapidity gaps in the final state for the meson case. In contrast, for heavy quarks we have only one rapidity gap.
In Table II we also present the result for coherent proton-nucleus reactions at LHC, considering both light (Ar) and heavy (Pb) ions. It should be noticed that the coherent production of mesons is currently measured at RHIC for the AA case. Therefore, the present estimation could be tested in dAu collisions with a good experimental feasibility. Disregarding the parton saturation phenomenon, we have that the predictions for J/Y production are somewhat similar to the linear pQCD approach, which is expected, since the heavy vector meson production is dominated by small pair separations, where the saturation physics does not contribute significantly. On the other hand, the photoproduction of r mesons is dominated by physics below saturation scale. The difference is huge, as shown in Ref. , reaching 3 orders of magnitude in rapidity distribution, which demonstrate the importance of the saturation physics on this process.
 V. P. Goncalves and M. V. T. Machado, Mod. Phys. Lett. 19, 2525 (2004). [ Links ]
 V. P. Goncalves and M. V. T. Machado, J. Phys. G: Nucl. Part. Phys. 32, 295 (2006). [ Links ]
 V.P. Goncalves and M.V.T. Machado, Eur. Phys. J. C 31, 371 (2003). [ Links ]
 V. P. Goncalves and M. V. T. Machado, Eur. Phys. J. C 40, 519 (2005). [ Links ]
 V. P. Goncalves and M. V. T. Machado, Phys. Rev. C 73, 044902 (2006). [ Links ]
 Bo Andersson et al. [Small x Collaboration], Eur. Phys. J. C 25, 77 (2002). [ Links ]
 N. Armesto, Eur. Phys. J. C 26, 35 (2002). [ Links ]
 N. N. Nikolaev and B. G. Zakharov, Z. Phys. C 49, 607 (1991). [ Links ]
 V. Barone and E. Predazzi, High-Energy Particle Diffraction, Springer-Verlag, 2002. [ Links ]
 E. Iancu, K. Itakura, and S. Munier, Phys. Lett. B 590, 199 (2004). [ Links ]
 V. P. Goncalves and M. V. T. Machado, Eur. Phys. J. C 38, 319 (2004). [ Links ]
 H. G. Dosch, T. Gousset, G. Kulzinger, and H. J. Pirner, Phys. Rev. D55, 2602 (1997). [ Links ]
 J. R. Forshaw, R. Sandapen, and G. Shaw, Phys. Rev. D 69, 094013 (2004). [ Links ]
 S. R. Klein, J. Nystrand, Phys. Rev. C 60, 014903 (1999). [ Links ]
 C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 89, 272302 (2002). [ Links ]
 D. d'Enterria [PHENIX Collaboration], arXiv:nucl-ex/0601001. [ Links ]
Received on 29 September, 2006