# Coherent and incoherent Upsilon production in ultraperipheral collisions at the Large Hadron Collider

###### Abstract

The exclusive photoproduction of states were calculated in ultra-peripheral collisions for coherent and incoherent process in PbPb at = 5.5 TeV. Different dipole models were compared in the theoretical framework of light-cone color dipole formalism. Moreover, it was calculated the differential cross section for the Upsilon states and their total cross section for two intervals of rapidity: 4 and y . A systematic study is done on the theoretical uncertainties associated to the production and predictions are presented for the first time for the incoherent cross section of the radially excited states.

###### pacs:

12.38.Bx; 13.60.Hb## I Introduction

The study of exclusive meson photoproduction in UPC upcs () is an essential tool to understand the low-x physics and also to investigate the gluon density in this regime. In the UPC case, the exclusive photoproduction dominates the process through the emission of quasi-virtual photons which interact with the target. The photon-target interaction amplitude, when considering the light-cone dipole formalism nik (), can be written as a convolution between the photon-meson wave functions overlap and the elementary dipole-target cross section Nemchik:1996pp (). In addition, the process considered here is quasi-elastic, 0, and in the region of small-x, the gluon density may increase to the point where gluon fusion, gg g, becomes significant. This kind of fusion produces nonlinear effects in the evolution equations. For instance, at mid-rapidity the typical value of Bjorken variable is for PbPb collisions at the LHC. Differently from DGLAP equations which is a linear equation, dipole models incorporate linear and non-linear effectsCGC (). In pQCD the exclusive meson photoproduction has a differential cross-section [, where the gluon distribution functions do not take into account the effects of saturation since they are evolved by DGLAP equations. The inclusion of parton saturation and nuclear effects are crucial in describing even the experimental observations in heavy-ion collisions at RHIC (for instance, see Ref. heavy ()). Within the color dipole approach one can introduce information on dynamics beyond the leading logarithmic QCD approach for meson production and computing predictions for the radially excited states is a reasonably easy task Nemchik:1996pp (). From the experimental point of view, the considered process is quite clear due to the presence of two rapidity gaps. The rapidity gap describes a region between the axis beam and the decay of meson free of particles. In this way, the experimental observation of this kind of quasi-diffractive process is facilitated. In addition, the absorption corrections in this case is not strong as in the corresponding final state in proton-proton collisions.

In this work, we investigate the exclusive (coherent and incoherent) production of and its radially excited states and in PbPb collisions for LHC energy. In a previous work GKMM () the coherent photoproduction of states at various energies in , pPb PbPb collisions at the LHC has been considered. Those calculations were carried out in the theoretical framework of the color light-cone dipole formalism nik () and focused only on the coherent channel where the initial state particles remain intact after interaction. It was shown that the corresponding predictions describe correctly the experimental results from LHCb Collaboration UpsilonLHCb () for photoproduction in collisions. Those data were obtained for typically large rapidities and the -values to be covered are increasingly smaller for forward rapidities. We roughly get at and it is clear that gluon dynamics is being probed at extremely low- and low perturbative scales GeV. This kinematical range is in the limit of application of usual pQCD and saturation approach should be relevant. For nuclear targets, the nuclear saturation scale (with ) reaches GeV in those cases Dusling (). The main novelty in current work is the detailed study of the incoherent cross section for the upsilon states. This is quite important, as it was pointed out that incoherent diffraction probes the fluctuations in the interactions strengths of multiparton Fock states in the nuclear wavefunctions Lappi (). The connection between incoherent diffraction and fluctuations is a quite rich subject and the pioneering works are found in Refs. FlucHis (). Recently, the topic is very active and we call attention to the following works recworks (); gbwnew ().

The paper is organized as follows. In the next section we give the main theoretical information to obtain the rapidity distribution of coherent and incoherent production of states in PbPb collisions for the future LHC run and energies close to the run2. The main motivation is the successful description of experimental results measured by LHCb Collaboration UpsilonLHCb () for in collisions. In the section III we present the phenomenological calculations, discuss the main theoretical uncertainties and a comparison with other approaches is performed. For instance, we compare the present calculation to the predictions available using the STARlight Monte Carlo STARlight (); KN (). Finally, we show the main conclusions.

## Ii Theoretical Framework

The exclusive meson photoproduction in nucleus-nucleus collisions can be factorized in terms of the equivalent flux of photons of the nucleus projectile and photon-target production cross section upcs (). In UPCs there is the absence of strong interactions between the projectile particle and the target. In this case the reaction is characterized by impact parameter 2 and as the interaction is ultra-relativistic and purely electromagnetic, one can use the Weizsäcker-Williams approximation upcs (). The photon energy spectrum, , which depends on the photon energy , is well known upcs (). The rapidity distribution for states photoproduction in collisions can be written down as,

(1) | |||||

where the photon flux in nucleus is denoted by and (coherent case) or (incoherent case). The symbol denotes the large rapidity gap between the produced meson and the final states nucleus.

The produced state with mass has rapidity and the square of the centre-of-mass energy is given by . The photon-Pomeron interaction will be described within the light-cone dipole frame, where the probing projectile fluctuates into a quark-antiquark pair with transverse separation (and momentum fraction ) long after the interaction, which then scatters off the hadron. The cross section for exclusive photoproduction of states off a nucleon target is given by,

(2) |

where and are the light-cone wavefunction of the photon and of the vector meson (), respectively. The dipole-proton cross section is denoted by and the diffractive slope parameter by . In this context, we are implicitly assuming that the proton shape is Gaussian and that the impact parameter dependence factorizes out from the dipole-nucleon scattering amplitude. Here, we consider the energy dependence of the slope using the Regge motivated expression GKMM ().

The exclusive photoproduction off nuclei for coherent and incoherent processes can be simply computed in high energies where the large coherence length is fairly valid. The expressions for both cases are given by Boris (),

where is the nuclear thickness function. In the numerical evaluations, we have considered the boosted Gaussian wavefunction and several phenomenological saturation models, which encode the main properties of the saturation approaches. Accordingly, the cross sections above include both the skwedness and real part of amplitude corrections. Namely, we multiply the result above by , where is the ratio of real to imaginary parts of the scattering amplitude and incorporates the off-forward correction (see GKMM () for details). The effective power on energy, is determined for each case. In order to take into account the threshold correction for the dipole cross section, we have multiplied them by a factor .

Finally, we set the parameters and phenomenological models to be considered in next section. For the slope parameter it was considered the energy dependency from the Regge phenomenology GKMM (),

(3) |

with GeV, GeV, GeV, GeV and GeV. It will be taken into account only for the incoherent cross section. For the meson wavefuntion, we will use the Boosted-Gaussian model wfbg () because it can be applied in a systematic way for excited states. The corresponding function is given by Sanda2 (),

(4) |

with . The operator is defined by

(5) |

and it acts on the following generatrix function

(6) | |||||

The main physical quantity is the dipole scattering cross section. We consider the following phenomenological models in our analysis: GBW PRD59-014017 (), CGC PLB590-199 () and BCGC PRD74-074016 (). The GBW model is defined by the eikonal shape for the dipole cross section,

(7) |

where is a constant and GeV denotes the saturation scale. We also consider the CGC model PLB590-199 (), based in the Color Glass Condensate framework, in which gluon saturation effects are incorporated via an approximate solution of the Balitsky-Kovchegov equation CGC (). The expression for the CGC model is given by,

where is the effective anomalous dimension and one has the constant .

In order to investigate the theoretical uncertainty associated to the models for the dipole cross section, we use the original values (OLD label) of parameters for the fits including the charm contribution. That is, for GBW-OLD we follow Ref. PRD59-014017 (), for CGC-OLD Ref. Soyez () is considered and bCGC-OLD refers to Ref. PRD74-074016 (). The bCGC model uses the same functional form of Eq. (II) and replaces the saturation scale in the following way: . Moreover, we consider the updated version of those models, GBW-NEW gbwnew (), CGC-NEW Amir () and bCGC-NEW Amir (), respectively. A comment is in order here: the GBW-NEW parametrization is very different from other color dipole fits, as it includes energy evolution of the subnucleonic shape of the proton and it can potentially significantly affect the incoherent cross section. In particular, GBW-OLD and GBW-NEW are qualitatively very different and GBW-NEW was not fitted to all small- data (DESY-HERA) as discussed in details in Ref. gbwnew ().

## Iii Results and discussions

Let us start the analysis by computing the theoretical predictions for the coherent process for PbPb collisions at 5.5 TeV. Here, we disregard any absorptive corrections. In Figure 1, it is presented the results for photoproduction of states, including its radial excitations, taking into account the different models presented in the last section. The theoretical uncertainty is relatively large, being of order 15 % for the state (similar for the remaining and states). We could have an additional uncertainty related to the vector meson wave function, however in Ref. GSSM2 () it was shown that this is not the case for states (the overall theoretical uncertainty is within the experimental error bars in case UpsilonLHCb ()). For the main contribution, we have b for . The relative contribution of the excited states compared to the bound states is . We see that the relative normalization and the overall behavior is changed mostly at mid-rapidity when comparing the old and updated versions of the dipole cross sections (the deviation at large rapidities is less evident). Notice that the LHCb data for upsilon production in collisions is reproduced by all the models in the forward region UpsilonLHCb () as shown in Ref. GKMM (). Therefore,the current level of the experimental uncertainties does not allow us to make definitive statements about the precision of the models considered. For sake of completeness, we present the integrated cross sections considering distinct cuts on rapidity. In Table 1, we present the results for the full rapidity coverage, , and forward rapidities, . In both tables 1 and 2, we present only the updated versions of the dipole cross sections. For sake of completeness, we present the ratio of the cross sections, , as a function of the photon-nucleus centre-of-mass energy, . We present in Fig. 2 the result using the bCGC-NEW and GBW-NEW dipole cross sections. It was verified that the CGC-NEW result is quite similar to the bCGC-NEW one. We see a relative energy dependence, following the same trend as for the states armesto ().

We now focus on the incoherent reaction, . This is a new contribution to the literature concerning the upsilon production. The rapidity distribution is shown in Fig. 3 using the same notation as the previous figure. As already known, the incoherent cross section is smaller than the coherent one. The typical ratio is . For instance, we obtain b for . The theoretical uncertainty seems to be larger than in the coherent case. The integrated cross sections are shown in Table 2 in the rapidity ranges and .

process: PbPb | TeV | ||
---|---|---|---|

GBW | CGC | b-CGC | |

163.7 (60.8) | 171.9 (63.8) | 143 (53.1) | |

20.3 (7.8) | 22.0 (8.2) | 20.5 (7.7) | |

10.3 (3.9) | 11.9 (4.3) | 10.9 (4.1) |

process: PbPb | TeV | ||
---|---|---|---|

GBW | CGC | b-CGC | |

61.2 (25.6) | 58.5 (24.4) | 44.5 (18.6) | |

8.2 (3.4) | 8.0 (3.4) | 6.9 (2.9) | |

4.2 (1.8) | 4.2 (1.8) | 3.8 (1.6) |

The calculations performed above can be compared to another theoretical approaches available in the literature. Let start comparing them to the STARlight Monte Carlo STARlight (). For the coherent production, the predictions for the state ratios are lower that STARlight results. As discussed in Ref. GKMM (), the possible origin comes from the extrapolation of HERA-DATA and taking a fixed ratio for the distinct states in the Monte Carlo, whereas in current case the evolution on energy is dynamically generated by the parton saturation approach models and mostly by the meson wavefuntions for the radially excited states. The nuclear effects are also computed in a different way in the two formalisms. In STARLight, the nuclear shadowing is calculated using vector meson dominance (VDM) plus Glauber model for hadronic collisions. In our case, shadowing comes from the multiple scatterings of color dipoles and is described by the Glauber-Gribov approach. We verified that our results are also smaller that in Ref. STARlight (), which can be related to more shadowing in color dipole models compared to the VDM+Glauber approach. In Ref. AN () only the coherent contribution was computed and the theoretical uncertainty we have found in the color dipole approach is comparable to perturbative QCD formalism. Concerning similar dipole calculations, more recent investigations are available in Refs. GSSM2 (); NavarraUps (). In Ref. NavarraUps () only the coherent production has been considered at 5.02 TeV. The results are smaller than ours and the main reason is the wavefunction chosen (Light Cone Gaussian wavefunction which gives smaller overall normalization compared to Boosted Gaussian one). The authors in NavarraUps () did not investigated the theoretical uncertainty associated to the wavefunction and dipole cross sections as well (only the uncertainty coming from one model for dipole cross section was addressed). In Ref. GSSM2 () the theoretical uncertainty for the coherent and incoherent cross section was investigated. However, predictions for higher energies in PbPb collisions were not presented and only the state was considered (the results are consistent with ours in that case). Finally, we did not consider photonuclear breakup in the present study. We will consider them in future analysis as they are important and the distinct channels have been measured for and photoproduction in UPCs RhoALICE (); CMSJPSI (). This sort of analysis was recently done in Ref.Guzey (), where the coherent production was considered using a pQCD model with NLO accuracy. An important point discussed in Guzey () is that the large region gives the dominant contribution for 0nXn and XnXn channel and they probe larger photon-target centre-of-mass energy than the case without neutron tagging.

## Iv Summary

We presented the predictions of rapidity distribution and integrated cross sections for the states for the LHC run 2 energies. The rapidity intervals used in total cross section were selected to match with the rapidity coverage of LHCb and ALICE detectors both covering 2 . The main contribution is the computation of the incoherent cross section within the color dipole approach and Glauber-Gribov treatment of nuclear shadowing. The cross section for the excited states are also calculated in a consistent formalism where the wavefunction of and states are theoretically well constrained. The usual procedure in the literature involves only an extrapolation of DESY-HERA production ratios to the LHC energies. Our calculations are directly comparable to the STARLight calculation, where distinct procedures are involved in the computation of nuclear shadowing (VDM plus Glauber model versus color dipole plus Glauber-Gribov approach) and how to obtain the incoherent cross section.

###### Acknowledgements.

This work was partially financed by the Brazilian funding agency CNPq and Rio Grande do Sul funding agency FAPERGS.## References

- (1) G. Baur, K. Hencken, D. Trautmann, S. Sadovsky, Y. Kharlov, Phys. Rep. 364, 359 (2002); C. A. Bertulani, S. R. Klein and J. Nystrand, Ann. Rev. Nucl. Part. Sci. 55, 271 (2005).
- (2) N. N. Nikolaev, B. G. Zakharov, Phys. Lett. B 332, 184 (1994); Z. Phys. C 64, 631 (1994).
- (3) J. Nemchik, N. N. Nikolaev, E. Predazzi and B. G. Zakharov, Phys. Lett. B 374, 199 (1996).
- (4) F. Gelis, E. Iancu, J. Jalilian-Marian and R. Venugopalan, Ann. Rev. Nucl. Part. Sci. 60, 463 (2010); H. Weigert, Prog. Part. Nucl. Phys. 55, 461 (2005); J. Jalilian-Marian and Y. V. Kovchegov, Prog. Part. Nucl. Phys. 56, 104 (2006); A.L. Ayala, M.B. Gay Ducati and E.M. Levin, Nucl.Phys. B 493, 305 (1997).
- (5) A. Dainese, in Proceedings of the 38th International Symposium of Multiparticle Dynamics (ISMD2008): Hamburg, Germany. September 15-20 2008, p.118-124, arXiv:0902.0377 [hep-ph].
- (6) M. B. Gay Ducati, F. Köpp, M.V. T. Machado, and S. Martins. Phys. Rev. D 94, 094023 (2016).
- (7) R. Aaij et al. [LHCb Collaboration], JHEP 1509, 084 (2015).
- (8) K. Dusling , F. Gelis, T. Lappi and R. Venugopalan, Nucl. Phys. A 836, 159 (2010) .
- (9) T. Lappi, H. Mäntysaari, R. Venugopalan, Phys. Rev. Lett. 114, 082301 (2015). H. Mäntysaari and B. Schenke, Phys. Rev. Lett. 117, 052301 (2016).
- (10) L. Frankfurt, G.A. Miller and M. Strikman, Phys. Rev. Lett. 71, 2859 (1993); A. Caldwell and H. Kowalski, Phys. Rev. C 81, 025203 (2010).
- (11) Kirill Tuchin, Phys. Rev. C 79, 055206 (2009); T. Lappi and H. Mäntysaari, Phys. Rev. C 83 065202 (2011) ; T. Toll and T. Ullrich, Phys. Rev. C 87, 024913 (2013) ; H. Mäntysaari, B. Schenke, C. Shen and P. Tribedy, Phys. Lett. B 772, 681 (2017).
- (12) J. Cepila, J.G. Contreras, J. D. Tapia Takaki, Phys. Lett. B766 , 186 (2017).
- (13) S.R. Klein, J. Nystrand, J. Seger, Y. Gorbunov and J. Butterworth, Comput. Phys. Commun. 212, 258 (2017).
- (14) S. Klein and J. Nystrand, Phys. Rev. C60, 014903 (1999).
- (15) B. Z. Kopeliovich and B. G. Zakharov, Phys. Rev. D 44, 3466 (1991); Yu.P.Ivanov, J.Huefner, B.Z.Kopeliovich and A.V.Tarasov, AIP Conf.Proc. 660, 283 (2003).
- (16) J. Nemchik, N. N. Nikolaev, E. Predazzi and B. G. Zakharov, Z. Phys. C 75, 71 (1997) 71.
- (17) B.E. Cox, J. R. Forshaw and R. Sandapen, JHEP 0906 , 034 (2009).
- (18) K. Golec-Biernat and M. Wüsthoff, Phys. Rev. D59, 014017 (1998).
- (19) E. Iancu, K. Itakura, S. Munier, Phys. Lett, B590, 199 (2004).
- (20) H. Kowalski, L. Motyka and G. Watt, Phys. Rev. D74, 074016 (2006).
- (21) G. Soyez, Phys. Lett. B655 , 32 (2007).
- (22) A.H. Rezaeain and I. Schmidt, Phys. Rev. D88, 074016 (2013).
- (23) G. Sampaio dos Santos and M.V.T. Machado, J. Phys. G42, 105001 (2015).
- (24) N . Armesto,A. H. Rezaeian, Phys.Rev. D90, 054003 (2014).
- (25) A. Adeluyi and A. Nguyen, Phys. Rev. C 87, 027901 (2013).
- (26) V.P. Gonçalves, B.D. Moreira and F.S. Navarra, Phys. Rev. D95, 054011 (2017).
- (27) J. Adam et al. [ALICE Collaboration], JHEP 1509, 095 (2015).
- (28) V. Khachatryan et al. [CMS Collaboration], Phys. Lett. B 396, 772 (2017).
- (29) V. Guzey, E. Kryshen and M. Zhalov, Phys. Rev. C93, 055206 (2016).