cpcChinese Physics C1674-1137Chinese Physical Society and the Institute of High Energy Physics of the Chinese Academy of Sciences and the Institute of Modern Physics of the Chinese Academy of Sciences and IOP Publishing Ltdcpc_43_6_06410310.1088/1674-1137/43/6/06410343/6/064103PaperNuclear physicsPhotoproduction of J/ψ in non-single-diffractive p+p collisions *
Supported by the National Natural Science Foundation of China (11775213, 11505180, 11675168), the U.S. DOE Office of Science (DE-SC0012704) and MOST (2016YFE0104800)
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.
photoproductionjpsiRHICLHC25.20.Lj25.75.Dw
Article funded by SCOAP^{3}
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-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 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-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
where
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 point-like 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]:
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.
(color online) Two-dimensional distribution of the photon flux as function of the distance r and the energy of the photon
for p+p collisions at
= 200 GeV
The photoproduction cross-section
depends on the gluon density in the proton [16]. At mid-rapidity,
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-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. For example, the branching ratios of
or (
) have changed with time. To compare different experimental results on equal footing, all measurements are updated with the latest branching fractions (5.961 ± 0.032% for
, 5.971 ± 0.032% for
) [35]. In Ref. [25-27,29-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]:
(color online) Exclusive
photoproduction cross-section as a function of
from the world-wide experimental measurements. The black solid line with the gray band on top of it represents the parametrization discussed in the text.
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-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.
Summary of references for the world-wide data.
experiment
b
collision system
ALICE
[17,18]
p-Pb/pp
LHCb
[19,20]
[20]
PP
H1(2013)
[21]
[21]
ep
H1(2005)
[22]
[22]
ep
H1(2000)
[23]
[23]
ep
ZEUS
[24]
[24]
ep
EMC
[25]
[25]
Fe
BPF
[26,27]
Fe
E516
[28]
[28]
p
E401
[29]
[29]
p/d
E87
[30,31]
Be
E25
[32]
[32]
d
SLAC
[33]
d
cornell
[34]
Be
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:
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.
Results
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 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.
(color online) The rapidity distribution,
, of
from photoproduction in NSD p+p collisions at
= 200 GeV. The solid line is the total production, while the dashed/dotted lines represent the individual cross-section contributions from the two beams.
(color online) 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). 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].
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]:
where
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 world-wide
differential cross-section measurements. The measured
distributions can be phenomenologically described by:
where
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-wide measurements. The references of the data are summarized in 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:
(color online) The slope parameter (b) of exclusive
photoproduction as a function of
from the world-wide experimental measurements. The black solid line with gray band represents the parametrization discussed in the text.
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:
where
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 form-factor 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
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.
(color online) 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) . 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].
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.
ReferencesBraun-MunzingerP.StachelJ.SongL.-H.YanL.-W.LiuY.LuoX.XuN.WangH.-M.et alMatsuiT.SatzH.A. BertulaniC.R. KleinS.NystrandJ.RebyakovaV.StrikmanM.ZhalovM.EskolaK.PaukkunenH.SalgadoC.Adam (ALICE)J.et alZhaW.Kłusek-GawendaM.SzczurekA.ZhaW.et alShiW.ZhaW.ChenB.KraussF.GreinerM.SoffG.J. MohrP.N. TaylorB.B. NewellD.G. RyskinM.Abelev (ALICE)B.et alS. Acharya et al. (ALICE), (2018), arXiv: 1809.03235.Aaij (LHCb)R.et alR. Aaij et al (LHCb), (2018), arXiv: 1806.04079Alexa (H1)C.et alAktas (H1)A.et alAdloff (H1)C.et alChekanov (ZEUS)S.et alAubertJ.et alR. ClarkA.et alR. ClarkA.et alH. DenbyB.et alBinkleyM.et alKnappB.et alD. HolmesS.LeeW.E. WissJ.NashT.et alCameriniU.et alGittelmanB.et alTanabashi (Particle Data Group)M.et alGuzeyV.KryshenE.StrikmanM.et alW. Zha et al, (2018), arXiv: 1810.02064R. KleinS.NystrandJ.Abelev (ALICE)B.et alZha (ALICE)W.et alP. D. B. Collins, An Introduction to Regge Theory and High-Energy Physics, Cambridge Monographs on Mathematical Physics (Cambridge Univ. Press, Cambridge, UK, 2009)
R. KleinS.NystrandJ.Abelev (ALICE)B.et al