Estimation of $\Upsilon$(1S) production in ep process near threshold

The near-threshold photo-productions of heavy quarkonia are important ways to test the QCD-inspired models and to constrain the gluon distribution of nucleon in the large x region. Investigating the various models, we choose a photon-gluon fusion model and a pomeron exchange model for J/$\Psi$ photo-production near threshold, emphasising on the explanation of the recent experimental measurement by GlueX at JLab. We find that these two models are not only valid in a wide range of the center-of-mass energy of $\gamma$ and proton, but also can be generalized to describe the $\Upsilon$(1S) photo-production. Using these two models, we predict the electro-production cross-sections of $\Upsilon$(1S) at EicC to be 48 fb to 85 fb at the center-of-mass energy of 16.75 GeV.


I. INTRODUCTION
The electro-productions of J/Ψ and Υ(1S) are among the key reactions that are going to be measured on electron-ion colliders [1][2][3][4], because they are important not only for understanding the interaction mechanism of the heavy quarkonium, but also for probing the gluonic properties of the nucleon. These measurements will advance our understanding of quantum chromodynamics (QCD) which governs the properties of hadrons and the interactions involving hadrons. The large statistic of the exclusive J/Ψ production data at the hard scale is very helpful in extracting the generalized parton distribution of the gluon. The generated heavy quarkonia interact with the intact nucleon, hence it provide a good opportunity in the study of the bound state between the quarkonium and the nucleon. Last but not least, the scattering between the heavy quarkonia and the nucleon near threshold can be used to extract the QCD trace anomaly. The heavy quarkonium photo-production off the nucleon helps us to probe the nucleon structure, to test QCD dynamics at short distances and to verify QCD inspired models [5].
There are a lot of experiments that measured the J/Ψ photo-production, while the photo-and electroproduction data of Υ(1S) at low energy are much scarce. Electron ion collider in China (EicC) is optimized to have a 3.5 GeV polarized electron beam and a 20 GeV polarized proton beam [3,4], with the centre-of-mass energy reaching 20 GeV. EicC will fill up the the absence of experimental measurement of Υ(1S) production near threshold. Hence the main purpose of this work is to check the models of J/Ψ production and to make a bridge to the prediction of Υ(1S) production at EicC. Recent discovery of the penta-quark P c make the physicists in predicting the penta-quark P b in the bottom quark sector, reasoned from the flavor symmetry of heavy quarks. The expected penta-quark P b decays into the Υ(1S) and the nucleon. So, for a better searching of penta-quark state P b , a well understanding of the continuum process of exclusive Υ(1S) production is indispensable. The differential cross-section of the channel γp → J/Ψp is linked to the QCD trace anomaly b [6,7], and proton mass decomposition. In principle the study on Υ(1S) photoproduction near threshold better help us in extracting the trace anomaly b and in understanding the origins of the proton mass [8,9] compared to the J/Ψ channel.
Thanks to the electron accelerators, there are a lot of data available on the J/Ψ electro-and photo-production off the nucleon. The data are roughly classified into three energy regions as follows: high and intermediate energy ZEUS and H1 data [10][11][12], low energy Fermilab data [13,14], and near threshold SLAC and GlueX data [15,16]. How to explain the combined data of a wide energy range especially the recently published GlueX data [16] is a challenge.
There are some empirical formulas to explain the process of γp → J/Ψp. In Ref. [17], the two gluon and three gluon exchange models are introduced to describe the processes of correlated quarks γqq → ccqq and γqqq → ccqqq respectively. This model takes into account the fact that the three target quarks recombine into the final proton via its form factor after emission of two or three gluons and the coupling of the photon to the cc pair. Since the power term (1 − x) n in this model is valid at x 1 and the form factor of proton is parameterized as F 2 g = exp(1.13t) according to the experimental data measured at low beam energies [15,18], this model is only valid in a limited energy range near threshold. Besides, a parameterized STZ (Strikman-Tverskoy-Zhalov) model [19] is constructed to study the incoherent J/Ψ photo-production off a nucleus. This model shows a good agreement with almost all data due to the application of two heaviside functions. However such a purely empirical parametrization fitted with the J/Ψ data cannot be applied to the Υ(1S) production.
Proton is a complex system of valence quarks, sea quarks and gluons, hence it is difficult to figure out the explicit process of γp → J/Ψp. Nonetheless, a series of models has been proposed over the passed few decades arXiv:2002.05340v1 [hep-ph] 13 Feb 2020 to explain this process. Some of them are deduced based on QCD theory and the relevant diagrams. As gluon is the gauge boson that transmits the strong interaction, in spirit of a factorization theory, some models exploit the gluon component inside the proton to calculate the heavy vector meson productions. One of them is the photon-gluon fusion (PGF) model, which takes into account the photon fusing with a gluon from the proton and transforming into a pair of cc. Thus, the cross section is a convolution of the gluon distribution function (GDF) with the photon-gluon cross section σ γg→cc . On the one side, this model gives a good description of J/Ψ data for almost the whole energy region (see Ref. [5]). On the other side, this model also provides an experimental approach [20] to test the GDF obtained from the global fit of deep inelastic scattering data, such as the dynamical gluon distribution from IMParton16 [21]. In the dipole model, a colorless object is exchanged between the J/Ψ and the proton. Since J/Ψ does not couple to a meson and the photon, the pomeron (bunch of gluons) exchange is supposed to be dominant. Consequently, several pomeron exchange models are proposed to study the photo-production of vector mesons. Usually, the pomeron exchange model works well at high energy [5] due to the fact that the gluon distribution is larger than the quark distribution in this region. However, there are also some empirical pomeron exchange models which coincide well with the experimental data of all energy range. As Υ(1S) production needs high energy, the pomeron exchange model is appropriate in the study of the photo-production of Υ(1S).
In this paper, we use two different models to study the photo-productions of J/Ψ and Υ(1S), which are discussed in Sec. II and III respectively. On this basis, we predict the electro-production of Υ(1S) at EicC in Sec. IV. Lastly some discussions and a summary is given in Sec. V.

II. PHOTO-PRODUCTION OF J/Ψ
Different model on the cross-section of γp → J/Ψp has its own applicable energy region due to its specific physical mechanism. As discussed in the previous section, we notice that the PGF model and the pomeron exchange model are in good agreement with the old J/Ψ production data at both high and low energy. If they are also consistent with the recent published data by GlueX and SLAC [16], then they are suitable to reproduce the heavy quarkonia production from the threshold to a relatively high energy. As a consequence, these models are then thought to be valid to fit the Υ(1S) production data at high energy and to extrapolate the cross-sections to a lower energy, even near the production threshold. The reason is that both J/Ψ and Υ(1S) have the heavy constituents moving slowly inside and are of the same quantum numbers. In the photon-gluon fusion model, the photon fuses with a gluon emitted from the proton and split into a cc pair, which is shown in Fig. 1. The cross section of charmonium photo-production is given by [5], where σ γg→cc denotes the photon-gluon fusion cross section, g(x, m 2 J/Ψ ) is the gluon distribution of the proton at µ 2 = m 2 J/Ψ , and f is the probability of the cc pair color-neutralized into the J/Ψ meson. σ γg→cc takes the form as [20], in which α = 1/137 is the fine structure constant, e c and m c are the charge and the mass of charm quark respectively, α s = 0.5 is the strong coupling constant at µ 2 = m 2 J/Ψ , s is the squared invariant mass of the photongluon system, and β is defined as β 2 = 1 − 4m 2 c s [5]. The integration range over s is from the charmonium production threshold 4m 2 c to the open charm production threshold 4m 2 D . Once above the DD threshold, the energetic cc from the photon-gluon fusion forms into D mesons. In this work, we use the dynamical gluon distribution from IMParton16 [21].
Because one gluon carries a color octet, the cc pair must radiate one or more soft gluons before the hadronization into a color singlet state of J/Ψ wave function. Yet the PGF model does not implement explicitly this color-neutralization process, an adjustable parameter f is used to describe the probability of the corresponding hadronization. The parameter f is energy-dependent and goes down as the energy goes up. Therefore we apply a power function to describe the energy-dependence as f = f 0 W λ . Performing a fit to the experimental data [10,13,14,16], we get f 0 = 7.36 and λ = −3.08.
As showed in Fig. 3, the PGF model with the dynamical gluon distribution is consistent with the experimental data of J/Ψ photo-production. The phenomenological pomeron exchange model being considered here is the one proposed in Ref. [22] to describe the cross section of γp → V p. The schematic Feynman diagram of the model is shown in Fig. 2. The cross section of γp → V p is parameterized as, where the square factor is derived from the phase space and it decreases to zero at the threshold. The free parameters σ p and are determined to be 17.1 nb and 0.109 respectively, from a fit to the experimental data [10,13,14,16].  the PGF model prediction. The total cross-section data of SLAC are deduced from the differential cross-section data [16]. Both models give good explanations for the J/Ψ photo-production near the threshold.

III. PHOTO-PRODUCTION OF Υ(1S)
We use the same theoretical frameworks discussed in the previous section to predict the photo-production of Υ(1S), for they are successful. Parameters related to the charm quark and J/Ψ should be replaced with the bottom quark and Υ(1S). In the PGF model, the charge and the mass of bottom quark are 1/3 and 4.67 GeV respectively. The mass of Υ(1S) is 9.46 GeV, and the strong coupling constant is taken as α s = 0.2 at this mass scale. The GDF for this channel is at a higher scale as g(x, m 2 Υ ). The upper limit of the integral in PGF model is 4m 2 B accordingly. Then, the free parameters used to describe the hadronization into Υ(1S) are obtained to be f 0 = 0.328 and λ = −2.15, from a fit to the experimental data [11,[23][24][25][26]. For the pomeron exchange model, the free parameters are obtained to be σ p = 5.19 pb and = 0.787, which agrees with a previous study [22]. The comparisons of the model calculations with the experimental data on the Υ(1S) photo-production are shown in Fig. 4. Both models are consistent with the existing data amazingly. The qualities of the fits are χ 2 /N = 0.46 and χ 2 /N = 0.51 for the photon-gluon fusion model and the pomeron exchange model respectively. Therefore, we can use these two models to predict the electro-production of Υ(1S) at EicC.

IV. ELECTRO-PRODUCTION OF J/Ψ AND Υ(1S)
The cross section of ep → epV is expressed as, where the photon flux [27] is, which represents the probability of electron emitting the virtual photons. k denotes the energy of the virtual photon, Q 2 denotes the minus four momentum square of it and the minimum Q 2 is evaluated as Q 2 min = m 2 e k 2 /[E e (E e − k)]. Note that the photon flux from the electron beam is defined in the proton rest frame. The vector meson production by the virtual photon is Q 2dependent, and it is connected with the production by the real photon as [28], where f (M V ) is a Breit-Wigner function [29] and n = c 1 + c 2 (Q 2 + M 2 V ), with c 1 = 2.09 and c 2 = 0.0073 determined by the experimental data [30]. The integration over the virtual photon energy k ranges from the threshold of Υ(1S) production to the maximum k max = E e − 10m e . Since the enormous vector meson productions are at the low virtuality Q 2 , the integration over Q 2 is taken as 0 < Q 2 < 1 GeV 2 in the calculation. Using the formulas mentioned above, we finally obtain the electro-production cross section of Υ(1S) at the optimal energy of EicC, which is 48 fb for the PGF model and 85 fb for the pomeron exchange model. As the luminosity of EicC is about 4 × 10 33 cm −2 s −1 , its integral luminosity is around 50 fb −1 with a year of operation [3,4]. Assuming a detection efficiency of 0.2, there are about 480 to 850 Υ(1S) events which can be collected per year at EicC.

V. DISCUSSIONS AND SUMMARY
We investigate two models on the cross-sections of γp → J/ψp and γp → Υp. The photo-gluon fusion model is on the basis of parton degree of freedom, and the validity of the factorization framework. Our result shows that both the dynamical gluon distribution and the factorization assumption work fine for a wide range of energy of the system. This is due the heavy quark that is connected with the gluon is large (> 1 GeV). For the pomeron exchange model, the colorless component is exchanged. Single gluon information is embedded in the coherent soft gluons. Though it is not related to the parton picture of the hadron, the pomeron exchange model explains well the experimental data with a few parameters. According to Ref. [5,22], the two models in this work also work successfully in explaining the J/ψ photo-production at much higher energy. The pomeron exchange model is a little better in explaining the J/Ψ near-threshold production compared to the PGF model, and it is close to the result based on pomeron trajectory in Ref. [31]. Although there are just a little Υ(1S) photoproduction data in the low energy region, both models in this work match the data successfully.
Electro-production cross-sections of Υ(1S) are calculated with the photo-production cross-sections based on the PGF model and the pomeron exchange model. The extrapolations of the near-threshold cross-sections from the two models exhibit some differences, however they are at the same magnitude. Finally we provide a prediction of the number of Υ(1S) collected at EicC near the threshold with the overall detector efficiency of 20%, which is around 500 events per year at the center-ofmass energy of 16.75 GeV. With years of collections of the Υ(1S) events at several energy points near the threshold, the EicC machine will advance our knowledge about the QCD trace anomaly, the proton mass decomposition, and the exotic penta-quark baryon in the bottom quark sector.