Photoproduction of J/ψ in non-single-diffractive p+p collisions

Significant enhancements of production at very low transverse momenta were recently observed by the ALICE and STAR collaborations in peripheral hadronic A+A collisions. The anomalous excess points to coherent photon-nucleus interactions in violent hadronic heavy-ion collisions, which were conventionally studied only in ultra-peripheral collisions. Assuming that the coherent photoproduction is the underlying mechanism responsible for the excess observed in peripheral A+A collisions, its contribution in p+p collisions with nuclear overlap, i.e. non-single-diffractive collisions, is of particular interest. In this paper, we perform a calculation of exclusive photoproduction in non-single-diffractive p+p collisions at the RHIC and LHC energies based on the pQCD motivated parametrization using the world-wide experimental data, which could be further employed to improve the precision of the phenomenological calculations for photoproduction in A+A collisions. The differential rapidity and transverse momentum distributions of from photoproduction are presented. In comparison with the production from hadronic interactions, we find that the contribution of photoproduction is negligible.


Introduction
In ultra-relativistic heavy-ion collisions, the aim is the search for a new form of matter -the Quark-Gluon Plasma (QGP), which was predicted by the lattice Quantum Chromodynamics (QCD) calculation [1], and the study of its properties in laboratory [2][3][4]. Among the probes of QGP, suppression in hadronic heavy-ion collisions with respect to elementary p+p collisions has been suggested as a "smoking gun " signature of QGP formation [5] due to the color screening effect in the deconfined medium.
can also be generated by the intense electromagnetic fields that accompany the relativistic heavy ions via coherent photoproduction [6]. The coherently produced are expected to probe the nuclear gluon distribution at low Bjorken-x [7], for which there is still a considerable uncertainty [8]. Conventionally, the associated J/ψ physics of the photoproduction and hadronic production belong to different subject fields, and they are studied in ultra-peripheral collisions (UPC) and hadronic collisions independently. In UPC, only photoproduction and related physics is studied, since there is no hadronic interaction; analogously, in hadronic collisions, only hadronic production is expected.
Is coherent photoproduction really prohibited in hadronic collisions, where violent strong interactions occur? Recently, a significant excess of production at very low transverse momentum ( 0.3 GeV/c) was observed by the ALICE collaboration in peripheral hadronic Pb+Pb collisions at forward-rapidity [9], which can not be described by the hadronic production that is modified by the hot and cold medium effects. STAR made the same measurements in Au+Au collisions at = 200 GeV and U+U collisions at = 193 GeV [10], and also observed significant enhancements at very low in peripheral collision. The observed excesses exhibit the characteristics of coherent photoproduction and can be quantitatively explained by the theoretical calculations with coherent photon-nucleus production mechanism [11][12][13], which strongly suggests the existence of coherent photoproduction in hadronic collisions. If coherent photoproduction is the underlying mechanism responsible for the observed excesses in the hadronic A+A collisions, what is its contribution in hadronic p+p collisions? Can we observe the excess that originates from the same production mechanism in hadronic p+p collisions? If the contribution is significant, it would affect the pp baseline used for nuclear modification factor ( ) of , which would further bias our understanding of QGP extracted from suppression measurements. In this paper, we perform a calculation of exclusive photoproduction in non-single-diffractive (NSD) p+p collisions at the RHIC and LHC energies. The differential rapidity and transverse momentum distributions of from photoproduction are presented, and are compared to those from hadronic production.

Methodology
According to the equivalent photon approximation, the photoproduction rate in p+p collisions can be factorized into two parts: the photon flux, and the photon-proton cross-section. The cross-section can be written as: where is the photon energy, is the photon flux at energy , and is the photonuclear interaction cross-section for . For simplicity, to make an estimate of its contribution, we assume that the photoproduction process in non-single-diffractive p+p collisions is exactly the same as in UPC.
The photon flux induced by a proton can be modeled using the Weizsäcker-Williams method [14]. For the point-like charge distribution, the photon flux is given by the simple formula n(ω, r) is the flux of photons with energy at distance r from the center of a proton, is the electromagnetic coupling constant, , and is the Lorentz factor. Here, is a modified Bessel function. The pointlike assumption is appropriate in UPC. However, in NSD collisions, the two colliding protons come very close to each other and the proton internal structure should be taken into account. A generic formula for any charge distribution can be written as [14]:

F(q)
where the form-factor is the Fourier transform of the internal charge distribution in a proton. A dipole form is employed to describe the form-factor of a proton, defined as: where the parameter a is related to the root mean square charge radius of the proton ( : 0.8768 ± 0.0069 fm [15]) by the equation . Figure 1 shows the two-dimensional distribution of the photon flux induced in p+p collisions at = 200 GeV as a function of distance r and energy w with the dipole form-factor for the proton. One can observe that the photon flux drops rapidly as inside the proton.
The photoproduction cross-section depends on the gluon density in the proton [16]. At midrapidity, production is sensitive to gluons with x down to at RHIC and at the LHC. However, there is still a large uncertainty in such a region of x for different PDF sets. In this calculation, we use the world-wide experimental data [17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34] for the exclusive photoproduction to perform a parametrization of the cross-section estimates. The measurements of photoproduction have been performed for more than forty years. In such a long period, different experimental techniques have been utilized and different input information was available at the time of measurements.
γp E γp , 5.971 ± 0.032% for ) [35]. In Ref. [25][26][27][29][30][31][32][33][34], the cross-section is derived from the measurements of using the relation with the information for distribution from the world-wide measurements, where t is the four-momentum transfer in the process. One should be aware that the effects of nuclear breakup, proton excitation and the potential phase factor are neglected in the extrapolation process, which needs to be further investigated in a future effort. The cross-section so obtained as a function of center-of-mass energy ( ) is shown in Fig. 2. The data are fitted using the following pQCD motivated expression [36]: where the second term on the right-hand side of the equation represents the turning on action near the production threshold, and the last term contains the evolution of gluon distribution on the Bjorken-x. The values of the free parameters and are determined from the fit, resulting in nb and . The systematic uncertainties of the measurements from the same experiment must be highly correlated, however the correlation matrix could not be obtained from the corresponding references. Therefore the correlations are not included in the fit, which would underestimate the error bars of and . The parametrization with the most complete experimental data could also be employed to improve the precision of phenomenological calculations for photoproduction in A+A collisions such as in [11][12][13]37,38]. As shown in the figure, the parametrization describes the experimental measurements very well with . The references of the data are summarized in Table 1.
To effectively relate the NSD cross-section to its corresponding region in the impact parameter space, a Glauber like geometrical picture is employed in the calculations: where is the NSD probability as a function of the impact parameter b, is the density distribution for a proton in the transverse plane, and is the cross-section like parameter determined by the NSD cross-section. The density distribution for a proton is given by: ρ 0 √ s where is the normalization factor. The parametrization formula for the density distribution is consistent with the dipole form-factor given in Eq. (4). There are two components in NSD interactions: colored hadronic interactions, and double-diffractive (DD) interactions. The two classes of interactions have different impact parameters. However, for simplicity, we do not make a distinction between these two types of interactions here, which needs to be further investigated in a future work. In this paper, we perform calculations for NSD p+p collisions at = 0.2, 2.76, 5.02 TeV, and 14 TeV. The corresponding NSD cross-sections are 30, 50, 56, and 64 mb [39], respectively.  With the convolution of equivalent photon spectra and elementary cross-section, the probability to produce a with rapidity y in a collision with an impact parameter b can be given by: where is the effective photon flux with an impact parameter b at photon energy . The effective photon flux, , can be expressed by the photon flux induced by one proton and the effective strength for a photon with a second proton: where b is the impact parameter between the two colliding protons, r is the distance from the proton which emits the photon, and the extra ensures a collision between the photon and the proton. The photon energy, , can be determined from the rapidity of , y: A complication is that either beam particle is equally likely to produce the photon; the cross-sections for these two possibilities from two beam directions are added: where can be obtained from Eq. (6). Figure 3 shows the calculated rapidity distribution of produced from photoproduction in NSD p+p collisions at = 200 GeV. The solid line is the total production cross-section, while the dashed/dotted lines represent individual contributions from the two beam protons. The ω γp E γp rapidity distribution is determined by the evolution of the photon flux with photon energy , and by the elementary cross-section at center-of-mass energy at different rapidities. Figure 4 shows the differential cross-section of from hadronic production and photoproduction as a function of rapidity in p+p collisions at = 0.2 (a), 2.76 (b), 5.02 TeV (c), and 14 TeV (d), respectively. The red and blue dashed lines are predictions from photoproduction with and without interference effect, respectively. The effect of interference will be discussed in detail later. The calculations are performed for NSD collisions, in which violent strong interactions exist that produce . The black solid lines with gray bands in the plots represent cross-sections from hadronic production. The hadronic contributions are extracted from the parametrizations using the world-wide experimental data, as described in Ref. [40]. The rapidity distributions from photoproduction are different at different collision energies due to the evolution of the two component structures (shown in Fig. 3), and the interference from the two beams. In comparison with the contribution from hadronic interactions, the yield from photoproduction is several orders of magnitude smaller, which makes the detection of photoproduction in NSD p+p collisions very difficult. Could we observe an excess of at low in NSD p+p collisions, similar to those in peripheral A+A collisions? Although the total cross-section from photoproduction is very small in comparison to the hadronic contribution, from photoproduction are mainly produced at low , which may have certain significance. of from photoproduction in p+p collisions depends on of the photon and acquired when the vector meson is created; the latter is dominant.
of the photon induced by a proton can be given by the equivalent photon approximation [14]: is the proton form-factor used previously, is the dimensionless normalization factor, and is the transverse momentum of the photon.
from the vector meson production can be estimated from the worldwide differential cross-section measurements. The measured distributions can be phenomenologically described by: is the normalization factor, b is the slope parameter depending on the center-of-mass energy ( ). Figure 5 shows the slope parameter (b) of exclusive photoproduction as a function of from the world-   Table 1. The black solid line with gray band represents the parametrization discussed below. Within the framework of Regge phenomenology [41], the slope parameter b should increase logarithmically with . Therefore, the data are fitted using the following expression: where and are free parameters. The corresponding is uniquely determined by the rapidity of and the pp collision energy. As shown in the figure, this expression describes the data reasonably well. We assume that the photon and that from the vector meson production are randomly oriented.
For , it is impossible to distinguish which proton emits the photon and which acts as a target. Due to the negative parity of , the sign of the two amplitudes are opposite, leading to destructive interference. The interference of vector meson production in UPC has been studied in detail by Klein and Nystrand [42]. We follow the same strategy to calculate the effect of interference: is the amplitude for production at rapidity y with transverse momentum . Figure 6 shows the differential invariant cross-section of from hadronic production and photoproduction as a function of transverse momentum in p+p collisions for mid-rapidity ( ) at = 0.2 (a), 2.76 (b), 5.02 TeV (c), and 14 TeV (d), respectively. The interference effect has been incorporated in the calculations of photoproduction. The red and blue dashed lines are predictions from photoproduction with and without interference effect, respectively. The black solid lines with gray bands represent cross-sections from hadronic production. The hadronic contributions are extracted from the parametrizations using the world-wide experimental data, as described in Ref. [40]. As depicted in the figure, the photoproduction contribution is several orders of magnitude smaller than from hadronic interactions, which means that the excess from photoproduction at low is not visible in NSD p+p collisions at the RHIC and LHC energies. Why is this the case? The photoproduction of is proportional to , which means that the photoproduction in p+p is of that in A+A collisions. The formfactor difference between p and A compensates the gap between p+p and A+A up to a certain extent, but it is not enough. For hadronic production, the production in p+p collisions is of that in A+A collisions, where ranges from 1 to 1000 depending on the collision species and centralities. Thus, the contribution from photoproduction is negligible in comparison with that from hadronic production. However, this is good news for the current measurements for very low in A+A collisions at RHIC and LHC, since the pp baseline for such a p T p T region comes from extrapolations from the relatively high measurements, which ignore the possible excess originating from photoproduction. The coherent photoproduction contribution could be increased by selecting events with a low charged particle multiplicity [43]. However, there are difficulties in relating the event multiplicity to the impact parameter distribution in p+p collisions, which needs to be further explored in a future work.

Conclusion
In summary, we performed a calculation of exclusive photoproduction in NSD p+p collisions at the RHIC and LHC energies. The differential rapidity and transverse momentum distributions of from photoproduction were presented. In comparison with the production from hadronic interactions, the contribution of photoproduction is negligible, which suggests that, in contrast with peripheral A+A collisions, the excess of yield from photoproduction is not visible in NSD p+p collisions.
We thank Dr. Spencer Klein and Prof. Pengfei Zhuang for useful discussions. The red and blue dashed lines are predictions from photoproduction with and without interference effect, respectively. The black solid lines with gray bands represent cross-sections from hadronic production. The hadronic contributions are from parametrizations in Ref. [40].