
In ultrarelativistic heavyion collisions, one aims at searching for a new form of matter  the QuarkGluon Plasma (QGP), which was predicted by the lattice Quantum Chromodynamics (QCD) calculation [1], and studying its properties in laboratory [24]. Among the probes of QGP, J/
$ \psi $ suppression in hadronic heavyion collisions with respect to that in 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. J/$ \psi $ can also be generated by the intense electromagnetic fields accompanied with the relativistic heavy ions via coherent photoproduction [6]. The coherently produced J/$ \psi $ are expected to probe the nuclear gluon distribution at low Bjorkenx [7], for which there is still considerably large uncertainty [8]. Conventionally, the associate physics extracted from J/$ \psi $ photoproduction and hadronic production belong to different subject field, and they are studied in ultraperipheral collisions (UPC) and hadronic collisions individually. In UPC, only photoproduction and related physics were studied, since there is no hadronic interaction; analogously, in hadronic collisions, only hadronic production was expected.Is the coherent photoproduction really prohibited in hadronic collisions, where the violent strong interactions occur? Recently, a significant excess of J/
$ \psi $ production at very low transverse momentum ($ p_{T} < $ 0.3 GeV/c) has been observed by the ALICE collaboration in peripheral hadronic Pb+Pb collisions at forwardrapidity [9], which can not be described by the hadronic production modified by the hot and cold medium effects. STAR made the same measurements in Au+Au collisions at$ \sqrt{s_{\rm{NN}}}$ = 200 GeV and U+U collisions at$ \sqrt{s_{\rm{NN}}}$ = 193 GeV [10], and also observed significant enhancements at very low$ p_{T} $ in peripheral collision. The observed excesses reveal characteristics of coherent photoproduction and can be quantitatively explained by the theoretical calculations with coherent photonnucleus production mechanism [1113], which strongly suggests the existence of coherent photoproduction in hadronic collisions. If coherent photoproduction is the underlying mechanism responsible for the observed excesses in hadronic A+A collisions, how about its contribution in hadronic p+p collision? Can we observe the excess originated 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 ($ R_{\rm{AA}} $ ) of J$ /\psi $ , which would further bias our understanding of QGP extracted from J/$ \psi $ suppression measurements. In this paper, we perform a calculation of exclusive J/$ \psi $ photoproduction in nonsinglediffractive (NSD) p+p collisions at RHIC and LHC energies. The differential rapidity and transverse momentum distributions of J/$ \psi $ from photoproduction are presented, which will be compared to those from hadronic production.According to the equivalent photon approximation, the photoproduction rate in p+p collisions can be factorized into two part: the photon flux, and the photonproton cross section. The cross section can be written as:
$ \sigma({p + p} \rightarrow {p + p} + {\rm {J}}/\psi) = \int {\rm d} \omega n(\omega) \sigma ( \gamma p \rightarrow {\rm {J}}/\psi p), $
(1) where
$ \omega $ is the photon energy,$ n(\omega) $ is the photon flux at energy$ \omega $ , and$ \sigma(\gamma p \rightarrow {\rm {J}}/\psi p) $ is the photonuclear interaction crosssection for J$ /\psi $ . For simplicity, we assume that the photoproduction process in nonsinglediffractive p+p collisions is exactly the same as that in UPCs to make estimations on the contribution.The photon flux induced by proton can be modelled using the WeizsäckerWilliams method [14]. For the pointlike charge distribution, the photon flux is given by the simple formula
$ n(\omega,r) = \frac{{\rm d}^{3}N}{{\rm d}\omega {\rm d}^{2}r} = \frac{Z^{2} \alpha}{ \pi^{2} \omega r^{2}}x^{2}{{K_{1}}^{2}}(x), $
(2) where
$ n(\omega,r) $ is the flux of photons with energy$ \omega $ at distance r from the center of proton,$ \alpha $ is the electromagnetic coupling constant,$ x = \omega r/\gamma $ , and$ \gamma $ is lorentz factor. Here,$ K_{1} $ 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, the proton internal structure should be taken into account. A generic formula for any charge distribution can be written as [14]:$ \begin{split} n(\omega,r) =& \frac{4Z^{2}\alpha}{\omega} \bigg  \int \frac{{\rm d}^{2}q_{\bot}}{(2\pi)^{2}}q_{\bot} \frac{F(q)}{q^{2}} {\rm e}^{{\rm i}q_{\bot} \cdot r} \bigg ^{2}, \\ q =& \left(q_{\bot},\frac{\omega}{\gamma}\right), \end{split} $
(3) where the form factor
$ F(q) $ is Fourier transform of the internal charge distribution in proton. A dipole form is employed to describe the form factor of proton, defined as:$ F(q) = \left(1+\frac{q^{2}}{a^{2}}\right)^{2}, $
(4) where the parameter a is related to the root mean square charge radius of the proton (
$ r_{\rm p} $ : 0.8768 ± 0.0069 fm [15]) by the equation$ a = \displaystyle\frac{\sqrt{12}\hbar c}{r_{\rm p}} $ . Figure 1 shows the twodimensional distributions of the photon flux induced in p+p collisions at$ \sqrt{s} $ = 200 GeV as a function of distant r and energy w with the dipole form factor for proton. One can observe that the photon flux drops rapidly toward$ r \rightarrow 0 $ inside the proton.Figure 1. (color online) Twodimensional distributions of the photon flux in the distant r and in the energy of photon
$ \omega $ for p+p collisions at$ \sqrt{s}$ = 200 GeVThe photoproduction cross sections,
$ \sigma (\gamma p \rightarrow {\rm{J}}/\psi p) $ , depend on the gluon density in the proton [35]. At midrapidity, J/$ \psi $ production is sensitive to gluon with x down to$ 1.5 \times 10^{2} $ at RHIC and$ 6 \times 10^{4} $ at the LHC. However, there is still large uncertainty at such x region for different PDF sets. In this calculation, we use the worldwide experimental data on exclusive J/$ \psi $ photoproduction to do the parametrization for cross section estimates. The measurements of J/$ \psi $ 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. For example, the branching ratio of$ {\rm{J}}/\psi \rightarrow e^{+}e^{} $ or ($ \mu^{+}\mu^{} $ ) have changed with time. To compare the different experimental results on an equal footing, all the measurements are updated with the latest branching fractions (5.961 ± 0.032% for J$ /\psi \rightarrow e^{+} + e^{} $ , 5.971 ± 0.032% for J$ /\psi \rightarrow \mu^{+} + \mu^{} $ ) [36]. In Ref. [2426,2833], the cross sections$ \sigma (\gamma p \rightarrow {\rm{J}}/\psi p) $ are derived from the measurements of$ \gamma A $ using the relation$ \displaystyle\frac{{\rm d}\sigma}{{\rm d}t}(\gamma + p)_{t = 0} = \displaystyle\frac{1}{A^{2}}\displaystyle\frac{{\rm d}\sigma}{{\rm d}t}(\gamma + A)_{t = 0} $ with the information of$ \displaystyle\frac{{\rm d}\sigma}{{\rm d}t}(\gamma + p) $ distribution from worldwide measurements, where t is the fourmomentum transfer in the process. It 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 the future efforts. The treated cross section as a function of$ \gamma p $ center of mass energy ($ E_{\gamma p} $ ) is shown in Fig. 2. The data are fitted using the following pQCD motivated expression [37]:Figure 2. (color online) Exclusive J
$ /\psi $ photoproduction cross section as a function of$ {{E}_{\gamma p}}$ from worldwide experimental measurements. The black solid line with gray band on top of it represents the parametrization discussed in the text.$ \sigma \left( E_{{\gamma p}} \right) = \mathrm{C_0} \left( 1 \frac{\left( m_{\rm p}+M_{\rm{J}/\psi} \right)^2}{E^2_{{\gamma p}}} \right)^{1.5} \left( \frac{E^2_{{\gamma p}}}{100^{2} \rm{GeV}^{2}} \right)^{\delta}, $
(5) where the second term on the righthand side of equation represents the turning on action near the threshold of production, and the last term reveals the evolution of gluon distribution on Bjorkenx. The values of the free parameters
$ C_{0} $ and$ \delta $ are determined from the fit, resulting in$ C_{0} = 80.2 \pm 0.9 $ nb and$ \delta = 0.321 \pm 0.005 $ . 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 fitting, which would underestimate the error bars of$ C_{0} $ and$ \delta $ . The parametrization with the most complete experimental data could also be employed to improve the precision of phenomenal calculations for photoproduction in A+A collisions such as [1113,38,39]. As shown in the figure, the parametrization describes the experimental measurements very well with$ \chi^{2}/NDF = 113.6/116 $ . The references of the data are summarized in Table 1.experiment $ \sigma $ b collision system ALICE [16,17] pPb/pp LHCb [18,19] [19] PP H1(2013) [20] [20] ep H1(2005) [21] [21] ep H1(2000) [22] [22] ep ZEUS [23] [23] ep EMC [24] [24] $ \mu $ FeBPF [25,26] $ \mu $ FeE516 [27] [27] $ \gamma $ pE401 [28] [28] $ \gamma $ p/dE87 [29,30] $ \gamma $ BeE25 [31] [31] $ \gamma $ dSLAC [32] $ \gamma $ dcornell [33] $ \gamma $ BeTable 1. Summary of references for worldwide data.
To efficiently relate the NSD cross section to its corresponding region in impact parameter space, a Glauber like geometrical picture is employed in the calculation:
$ \begin{split} \sigma_{\rm{NSD}} =& \int_{0}^{\infty}2 \pi b P_{\rm{NSD}}(b) {\rm d}b, \\ P_{\rm{NSD}}(b) =& \int {\rm d}^{2}s T(\vec{s})(1{\rm e}^{\sigma_{0}T(\vec{s}\vec{b})}) \end{split}, $
(6) where
$ P_{\rm{NSD}}(b) $ is the NSD probability as a function of impact parameter b,$ T(\vec{s}) = \displaystyle\int_{\infty}^{+\infty}{\rm d}z \rho\left(r = \sqrt{s^{2}+z^{2}}\right) $ is the density distribution for proton in transverse plane, and$ \sigma_{0} $ is the crosssection like parameter determined by the NSD cross section. The density distribution for proton in volume is given by:$ \rho(r) = \rho^{0} {\rm e}^{ar}, $
(7) where
$ \rho^{0} $ is the normalization factor. The parametrization formula for density distribution is consistent with the dipole form factor given in Eq. (4). There are two components in NSD interactions: colored hadronic interactions, and doublediffractive (DD) interactions. The two classes of interactions have different impact parameters, however, for simplicity, we do not make a distinction between these two type of interactions here, which needs to be further investigated in the future work. In this paper, we perform calculations in NSD p+p collisions at$ \sqrt{s}$ = 0.2, 2.76, 5.02 TeV, and 14 TeV, the corresponding NSD cross sections are 30, 50, 56, and 64 mb [40], respectively.With the convolution of equivalent photon spectra and elementary
$ \gamma p \rightarrow {\rm{J}}/\psi p $ cross section, the probability to produce a J/$ \psi $ with rapidity y for a collision at impact parameter b can be given by:$ \frac{{\rm d}P(y,b)}{{\rm d}y} = \omega N(\omega, b) \sigma_{\gamma p \rightarrow {\rm{J}}/\psi p}(E_{\gamma p}), $
(8) where
$ N(\omega,b) $ is the effective photon flux with impact parameter b at photon energy$ \omega $ . The effective photon flux,$ N(\omega, b) $ , can be expressed through the photon flux induced by one proton and effective strength for the photon with the second proton:$ N(\omega,b) = \int n(\omega, r) \frac{\theta(r_{p}(\vec{r}  \vec{b}))}{\pi r_{p}^{2}}{\rm d}^{2}r, $
(9) where b is the impact parameter between the two colliding protons, r is the distant from the proton which emits the photon, and the extra
$ \theta(r_{\rm p}(\vec{r}  \vec{b})) $ ensures collision when the photon hits the proton. The photon energy,$ \omega $ , can be determined from the rapidity of J/$ \psi $ , y:$ \omega = \frac{1}{2} M_{\rm{J}/\psi} {\rm e}^{y}. $
(10) One 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:
$ \frac{{\rm d}\sigma}{{\rm d}y} = \int^{\infty}_{0} (\frac{{\rm d}P(y,b)}{{\rm d}y} + \frac{{\rm d}P(y,b)}{{\rm d}y})P_{\rm{NSD}}(b) 2 \pi b {\rm d}b, $
(11) where
$ P_{\rm{NSD}}(b) $ can be obtained from Eq. (6). Figure 3 shows the calculated rapidity distribution,$ {\rm d}\sigma/{\rm d}y $ , of produced J/$ \psi $ from photoproduction in NSD p+p collisions at$ \sqrt{s}$ = 200 GeV. The solid line is the total production, while the dashed/dotted lines represent the individual cross section contributions from the two beam protons. The rapidity distribution is determined by the evolution of photon flux with photon energy$ \omega $ and elementary$ \gamma p $ cross section with center of mass energy$ E_{\gamma p} $ at different rapidities.Figure 3. (color online) The rapidity distribution,
$ {\rm d}\sigma /{\rm d}y$ , of produced J$ /\psi $ from photoproduction in NSD p+p collisions at$ \sqrt{s}$ = 200 GeV. The solid line is the total production, while the dashed/dotted lines represent the individual cross section contributions from the two beams.Figure 4 shows the differential cross section of J/
$ \psi $ from hadronic production and photoproduction as a function of rapidity in p+p collisions at$ \sqrt{s} $ = 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 in the paper. The calculations are performed in NSD collisions, in which there exists violent strong interactions to produce J/$ \psi $ . The black solid lines with gray bands in the plots represent J/$ \psi $ cross sections from hadronic production. The hadronic contributions are extracted from parameterizations using the worldwide experimental data, as described in Ref [34]. The rapidity distributions from photoproduction are different in different collision energies due to evolution of the two component structures (shown in Fig. 3) and interference from the two beam directions. In comparison with the contribution from hadronic interactions, the yield from photoproduction is several orders of magnitude lower, which makes it very difficult to detect the possible J/$ \psi $ photoproduction in NSD p+p collisions.Figure 4. The differential cross section of J/
$ \psi $ from hadronic production and photoproduction as a function of rapidity in p+p collisions at$ \sqrt{s} $ = 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 black solid lines with gray bands represent cross sections from hadronic production. The hadronic contributions are from parameterizations in Ref [34].Could we observe an excess of J/
$ \psi $ at low$ p_{T} $ 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 that from hadronic contribution, the J/$ \psi $ from photoproduction is mainly produced at low$ p_{T} $ , which may gain certain significance. The J/$ \psi $ $ p_{T} $ from photoproduction in p+p collisions depends on the$ p_{T} $ of the photon and the$ p_{T} $ acquired when the vector meson is created; and the latter is dominant. The$ p_{T} $ of the photon induced by proton can be given by the equivalent photon approximation [14]:$ \frac{{\rm d}^{2}N_{\gamma}}{{\rm d}^{2}\vec{k}_{\gamma\bot}} = K_{0}\frac{F_{\gamma}^{2}(\vec{k}_{\gamma})\vec{k}^{2}_{\gamma\bot}}{(\vec{k}_{\gamma\bot}^{2}+\omega^{2}_{\gamma}/\gamma_{c}^{2})^{2}}, $
(12) where
$ F_{\gamma}(\vec{k}_{\gamma}) $ is the form factor of proton used previously,$ K_{0} $ is the dimensionless normalization factor, and$ \vec{k}_{\gamma\bot} $ is the transverse momentum of the photon. The$ p_{T} $ from the vector meson production can be estimated from the worldwide$ p_{T} $ differential cross section measurements. The measured$ p_{T} $ distributions can be phenomenally described by:$ \frac{{\rm d}\sigma}{{\rm d}p_{T}} = N_{0} p_{T} {\rm e}^{bp_{T}^{2}}, $
(13) where
$ N_{0} $ is the normalization factor, b is the slope parameter depending on the$ \gamma p $ center of mass energy ($ E_{\gamma p} $ ). Figure 5 shows the slope parameter (b) of exclusive J/$ \psi $ photoproduction as a function of$ E_{\gamma p} $ from worldwide experimental measurements. The references of the data are summarized in Table 1. The black solid line with gray band represents the parametrization discussed in the following. Within the framework of Regge phenomenology [41], the slope parameter b should increase logarithmically with$ E_{\gamma p} $ . Therefore the data are fitted using the following expression:Figure 5. (color online) The slope parameter (b) of exclusive J
$ /\psi $ photoproduction as a function of$ {{E}_{\gamma p}}$ from worldwide experimental measurements. The black solid line with gray band represents the parametrization discussed in the text.$ b = C_{0} + C_{1}{\rm{ln}}E_{\gamma p}, $
(14) where
$ C_{0} $ and$ C_{1} $ are free parameters. The corresponding$ E_{\gamma p} $ is uniquely determined by the rapidity of J/$ \psi $ and the collision energy of pp system. As demonstrated in the figure, the expression describes the data reasonably well. We assume that the photon$ p_{T} $ and that from vector meson production are randomly oriented.For
$ p_{T} < \hbar /b $ , it is impossible to distinguish which proton emitted the photon, and which acts as target. Due to the negative parity of J/$ \psi $ , the sign of 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:$\begin{split} \sigma(p_{T},y,b) =& A^{2}(p_{T},y,b) +A^{2}(p_{T},y,b) \\ & 2A(p_{T},y,b)A(p_{T},y,b) \times {\rm cos}(\vec{p}_{T} \cdot \vec{b}), \end{split}$
(15) where
$ A(y,p_{T},b) $ is the magnitude of amplitude for J/$ \psi $ production at rapidity y with transverse momentum$ p_{T} $ .Figure 6 shows the differential invariant cross section of J/
$ \psi $ from hadronic production and photoproduction as a function of transverse momentum in p+p collisions for midrapidity ($ y < 1 $ ) at$ \sqrt{s} $ = 0.2 (a), 2.76 (b), 5.02 TeV (c), and 14 TeV (d), respectively. The interference effect has been incorporated in the calculation of photoproduction. The red and blue dashed lines are predictions from photoproduction with and without interference efect, respectively. The black solid lines with gray bands represent cross sections from hadronic production. The hadronic contributions are extracted from parameterizations using the worldwide experimental data, as described in Ref [34]. As depicted in the figure, the photoproduction contribution is several orders of magnitude smaller than that from hadronic interactions, which means that the excess originated from photoproduction at low$ p_{T} $ is not visible in NSD p+p collisions at RHIC and LHC energies. Why is this case? The photoproduction of J$ /\psi $ is proportional to$ Z^{2}A^{2} $ , which means that the photoproduction in p+p is$ 1/Z^{2}A^{2} $ of that in A+A collisions. The form factor difference between p and A would compensate a certain amount of the gap between p+p and A+A, but, not enough to the same level. For hadronic production, the production in p+p collisions is$ 1/N_{\rm{coll}} $ of that in A+A collisions, where$ N_{\rm{coll}} $ range from 1 to 1000 depending on the collision species and centralities. Thus the contributions from J$ /\psi $ photoproduction is negligible in comparison with that from hadronic production. However, this is a good news for the current$ R_{\rm{AA}} $ measurements for very low$ p_{T} $ in A+A collisions at RHIC and LHIC, since the used pp baseline for such$ p_{T} $ region are from extrapolations utilizing the relative high$ p_{T} $ measurements, which ignores the possible excess originated from photoproduction. The coherent photoproduction contribution could be raised by selecting events with a low charged particle multiplicity [43], however, there are difficulties in relating the event multiplicity to impact parameter distribution in p+p collisions, which needs to be further explored in future work.Figure 6. (color online) The differential invariant cross section of J/
$ \psi $ from hadronic production and photoproduction as a function of transverse momentum in p+p collisions for midrapidity ($ y < 1 $ ) at$ \sqrt{s} $ = 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 black solid lines with gray bands represent cross sections from hadronic production. The hadronic contributions are from parameterizations in Ref [34].In summary, we perform a calculation of exclusive J/
$ \psi $ photoproduction in NSD p+p collisions at RHIC and LHC energies. The differential rapidity and transverse momentum distributions of J/$ \psi $ from photoproduction are presented. In comparison with the J/$ \psi $ production from hadronic interactions, the contribution of photoproduction is negligible, which suggests that, in contrast with the case in peripheral A+A collisions, the excess of J/$ \psi $ yield from photoproduction is not visible in NSD p+p collisions.
Photoproduction of J/$ \psi $ in nonsinglediffractive p+p collisions
 Received Date: 20181026
 Available Online: 20190601
Abstract: Recently, significant enhancements of J/