Identifying hidden charm pentaquark signal from non-resonant background in electron–proton scattering

In this study, we analyze the electroproduction of the LHCb pentaquark states with the assumption that they are resonant states. Our main concern is to investigate the final state distribution in the phase space to extract a feeble pentaquark signal from a large non-resonant background. The results indicate that the signal to background ratio will increase significantly with a proper kinematic cut, which will be beneficial for future experimental analysis.


Introduction
In the last few decades, numerous possible candidates for exotic Hadrons have been experimentally established. Specifically, in 2015, the LHCb Collaboration announced the observation of two pentaquark states, one narrow , and one broad , in the invariant mass distribution in the decay [1]. In 2019, the LHCb Collaboration updated knowledge of the pentaquarks, with multiple collected data samples of the same decay [2]. A new narrow pentaquark candidate, , was observed, whereas the old peak was found to be resolved into two narrower structures, and , owing to the larger statistics. In a coupled-channel approach, it is argued that the existence of a narrow is required by heavy quark spin symmetry [3]. After their discovery, numerous discussions in association with their properties were triggered, and various interpretations have been proposed for their internal structure. Because their masses are close to the thresholds, many studies assigned them as molecular states [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17]. Alternative explanations include hadro-charmonium states [18], compact diquark-diquark-antiquark states [11,[19][20][21], etc. A datadriven analysis of found it to be a virtual state [22]. In contrast, several recent works found that the spin parity assignments for (4440) and (4457) are sensitive to the details of the one-pion exchange potential [3,15,23,24]. The fits of the measured -invariant mass distributions indeed point to different quantum numbers for and [3]. Notably, a hidden charm pentaquark was predicted [25,26] before it was observed by the LHCb. Moreover, other possible pentaquarks in a strange and bottom sector were suggested [27][28][29][30]; however, they have not yet been observed experimentally. For more details, refer to the comprehensive reviews in [31][32][33][34][35][36][37][38][39].
However, since its discovery, it has been noted that the narrow peaks of the pentaquarks could be caused by triangle singularities [40,41], and it was recently suggested to distinguish them in isospin breaking decays [39,42]. Furthermore, their decay and production properties are extensively studied in various scenarios [43][44][45][46][47][48][49][50]. To discriminate their nature, the production of pentaquarks has been proposed in photo-induced [51][52][53][54][55] and pion-induced reactions [56][57][58][59], because a triangle singularity cannot be present in the two-body final states of the pro-Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI. Article funded by SCOAP 3 duction process. Thus, if they are observed in the or open charm production [52], they should be genuine states other than kinematic effects. Subsequently, an experimental study of pentaquarks through photoproduction was proposed at the JLab [60]. The GlueX Collaboration searched for the pentaquark states through the nearthreshold exclusive photoproduction off proton [61]. No evidence for pentaquark photoproduction was determined, and model-dependent upper limits on their branching fraction were set. The photoproduction rate was investigated in model-dependent calculations [47,62], where the coupling of the radiative decay was evaluated by the vector meson dominance (VMD) model. Although the extracted branching ratio, , is dependent on the details of the VMD, e.g., off-shell form factor, the photoproduction rate tends to not be high compared to the non-resonant contribution. Double polarization observables were proposed to be useful in the search of pentaquark photoproduction [45]. However, the LHCb results indicate a model-independent lower limit of [62]; thus, it is expected to observe pentaquark eletro and photoproduction after enough events are accumulated, if is a real resonant state. Moreover, it is expected that the distributions of from a pentaquark and Pomeron are different in large angles of the differential cross-sections [47,51], which could be helpful for identifying a pentaquark in cross-sections.

− 20
After the update of the JLab accelerator to 12 GeV, the search for pentaquark electroproduction at JLab12 will continue. Recently, an electron-ion collider at China (EicC) was proposed, where Hadron physics is a main concern [63]. Its designed center of mass (c.m.) energy of GeV covers charmonium electroproduction. In this study, we investigate the electroproduction of pentaquark states in these machines. The main concern here is on the final state distributions, from which a kinematic cut would isolate a feeble pentaquark signal from a large non-resonant background. This paper is organized as follows. In Sec. 2, we will briefly describe the analytic formalism in our computation, following which the results and discussions will be introduced in Sec. 3. Finally, we will present a short summary.

Formalism
As shown in Fig. 1, in the final states of can be produced from a pentaquark decay (b) and non-resonant -channel (a). Here, -channel contribution is from a or exchange, but both of them are negligible because of the highly off-shell intermediate state and the significantly small coupling between , respectively. Several phenomenological models were constructed to parameterize the -channel diagram with a gluon or Pomeron exchange. A detailed comparison of these models can be found in Ref. [64] for the photoproduction. Here, we employ the soft dipole Pomeron model, which can describe vector meson photoproduction from low to high energies [65]. We use a covariant orbital-spin (L-S) scheme to construct the Langrangians of the couplings [66], which has been used widely for normal and resonances [67][68][69].

Pomeron exchange
The Pomeron exchange model [70][71][72] accounts for the dominant contribution in the leptoproduction process. The Pomeron mediates the long-range interaction between the nucleon and confined (anti-)quarks within a quarkonium. This is an effective and useful model to parameterize the diffractive process for the production of neutral vector mesons in a high-energy region. By including a double Regge pole with an intercept equal to one, the soft dipole Pomeron model does not violate the unitarity bounds and can describe nearly all the available cross-section data of the photo and electroproduction of vector mesons, from light to heavy and near the threshold to a high-energy region in a consistent manner [65]. We start from the photoproduction of vector meson off proton in the soft dipole Pomeron model with the formula of the -dependent cross-section [65], where the amplitudes are defined as The and terms indicate the dipole Pomeron and Reggeon. and represent the photon virtuality and mass of a vector meson, respectively. Variable , where is the scattering angle of the final states in the c.m. system of . The nonlinear Pomeron trajectory is with as the pion mass, and the Reggeon trajectory is with and . The parameters, , , , and , can be obtained by fitting the vector meson photoproduction; more details regarding this can be found in Ref. [65]. The advantage of this model is that it includes exclusive photoproduction of all the vector mesons for both real and virtual photons, as shown in the above amplitudes. This is convenient for our calculation of electroproduction. The results for the photoproduction are shown in Fig. 2.
The electroproduction amplitude for the Pomeron exchange is evaluated as is a sub-reaction, , , and is a sub-reaction, . By neglecting the polarization correlations between the two sub-reactions, we obtain the amplitude square, where is the polarization vector of an intermediate photon with spin of -direction . The amplitude for sub-reaction, , can be determined from differential cross-section by the dipole Pomeron model mentioned in Eq. (1) with the relation, . An alternative approach to investigate the electroproduction of a vector meson is using a microscopic description of the Pomeron exchange [79,80].

Pentaquark
Here, we only consider with quantum number in line with GlueX [61], where the branching fractions are determined using the JPAC model [53]. A similar conclusion could be obtained for the alternative assignment of . The effective Lagrangian for the coupling of and is written as [47,67] where and denote resonance and a nucleon, respectively. The coupling constant, , can be determined from the corresponding decay widths. Here, we use the total decay widths of as the measured values by the LHCb and the upper limits of branching fractions, , determined by GlueX [61]. The propagator of can be written as where is the momentum of the propagator, and the mass, and the decay width of . The term, , is defined as γp → P c γp → V p → P c We assume that the pentaquark resonances couples to a photon via a vector meson pole using the VMD model. Therefore, the vertex can be considered as , as shown in Fig. 1(b). The coupling of the vector meson and photon is expressed as is the mass of the vector meson, and and are the vector meson and photon field, respectively.  , can be extracted from the partial decay width, , from the formula, p e where the masses of the electrons and positrons have been neglected, and is the three-vector momentum of an electron in the vector meson rest frame.
For an off-shell vector meson in the VMD, we choose the form factor, where is the cut-off parameter. The choice of the vector meson and cut-off parameter will not change the distribution of the final state, which is the main concern. Thus, we choose the vector meson to be and cutoff to be GeV.

20 z
We explore the electroproduction of the pentaquarks observed by the LHCb Collaboration as listed in Table 1, together with the contributions from the Pomeron exchange in the JLab12 and EicC energy configurations. JLab12 is a fixed-target experiment with 12-GeV electrons and rest protons, whereas the EicC is a colliding experiment with -GeV electrons and -GeV protons. The pentaquark and Pomeron contributions were added incoherently. The interference terms may have large contributions to the total cross-section, but distribute smoothly in the phase space. Furthermore, it is too premature to consider the interference at present, because we do not know the relative phase between the different contributions. Notably, these terms can be neglected for searching a pentaquark because we focus on a pentaquark-dominant phase space area, from which we obtain the main conclusion of this study. In our calculation, we choose the laboratory frame with an electron moving in the opposite direction. The cross-sections were evaluated by the VEGAS program [81], which numerically integrates the kinematic events generated by RAMBO [82] with the dynamics described by the formula above. We also obtain the final state distributions simultaneously.
The total production cross-sections for both JLab12 and EicC are summarized in Table 2. We can observe that the cross-sections of the non-resonant background are a few orders of magnitude larger than those of the pentaquarks. For the cross-sections of the pentaquarks, the model-dependent branching fractions determined by GlueX and the cut-off parameter in the form factor ap-pear as the overall factor, which suggests the cross-sections rather than the final state distributions are considerably model-dependent. Thus, the distributions are our main concern here. t t t The three-momentum and polar angle distributions of the final proton from either the Pomeron exchange process or pentaquark production are shown separately in Fig. 3 for JLab12 and in Fig. 4 for EicC. The left panels are for the final proton from the Pomeron contribution, whereas the right ones are for the proton from the pentaquark decay. For a better comparison, we use the same range in the axes of the two panels for each figure.
Because of the completely different energy configurations, the final proton moves in the electron and proton forward angles in JLab12 and EicC, respectively. Notably, the polar angle distributions of the final proton are significantly different for the Pomeron exchange and pentaquark production. This is because the -dependent cross-section in Eq. (1) is suppressed at large for the Pomeron exchange, whereas the shape from the pentaquarks is completely flat across the full range. This fact has been already noted in several papers [47,51,62], that is, the final particles from different contributions have different behaviours at large angles.
and (4440) overlap with each other because of the closeness of their masses; thus, a good energy resolution is needed to distinguish them. This is a challenge for future detector design. P c In particular, for each energy configuration, the protons from the non-resonant background and pentaquarks present significant differences in the phase space, as shown in Fig. 3 and Fig. 4. Consequently, we can take advantage of this feature to enhance the peaks relative to     Fig. 5, which shows the differential cross-section of the electroproduction process in the energy configurations of JLab12 and EicC. The dashed and solid curves in Fig. 5 show the results with and without the cut on the three-momentum and angle of the final proton, respectively. Note that a simple cut can remove more Pomeron contribution than a pentaquark in JLab12. In comparison, for EicC, a cut GeV and also works well to depress the background. Quantitatively, in the case of in JLab12, the signal to background ratio increases from 0.3 to 19 with the kinematic cut. Therefore, the kinematic cut can make the peaks more prominent and present a huge potential in experimental analysis, although the total number of events would decrease after the cuts are used. Thus, more complex cuts would further improve the situation. W Finally, we would like to emphasize the great potential of both EicC and JLab12 to search for a pentaquark. EicC has a higher signal to background ratio, whereas JLab12 has a much higher luminosity. The center of mass energy of EicC is approximately 16.7 GeV, wich is much larger than the 4.8 GeV of JLab12. As listed in Table 2, the larger center of mass energy would make the total cross-section 8 and 80 times larger for the pentaquark signal and non-resonant background, respectively. However, EicC has 15 times larger phase space in the invariant mass, , for the background than JLab. Therefore, the pentaquark signal could be presented more prominently in the differential cross-section in EicC in Fig. 5. Conversely, the differential cross-section is less reduced in EicC than that in JLab12 after the kinematic cut is employed, which shows that the colliding mode could be more useful to study a pentaquark than the fixed-target mode.

Summary
The GlueX Collaboration at the JLab has searched for pentaquark photoproduction and obtained negative res-ults at the present precision [61]. One possibility is that the pentaquark signal has a smaller total cross-section compared to that of the non-resonant contribution. The production rates have been already investigated in many efforts [47,62], and the signal of a pentaquark in hidden charm photoproduction would be considerably small in cross-sections.
In this study, we calculated the process with both a non-resonant -channel contribution and a hidden charm pentaquark in the -channel. After the nonresonant contribution was normalized using the soft dipole Pomeron model and photoproduction data, the distributions of the final particles from both the sources were investigated. In view of the different shapes of the final proton in the phase space, owing to the different underlying mechanisms, we proposed that the three-momentum and angle cuts on the proton could largely suppress the non-resonant contribution, whereas the signal to background ratio would be significantly increased. For both the energy configurations of JLab12 and EicC, we found promising strategies even with simple cuts, which will be beneficial for future experimental analysis of electroproduction in these machines. In addition, it is promising to search for pentaquarks in a higher energy collider, e.g., US EIC. A similar cut like the one used for EicC would also be helpful. Our criterion is also enlightening for the electroproduction of , bottom analog of states [64]. Here, we focus on the method to supress the background rather than the total production cross-section of the pentaquarks, because the total cross-section of the pentaqurks is greatly model-dependent, owing to the unknown coupling constant and cut-off parameter appearing as the overall factors. Finally, we noted that our framework could be used in a full simulation, the selection criterion of final particles, and the optimization of the detector design in the future.
Z. Yang gratefully acknowledges the hospitality at ITP where part of this work was performed. We are grateful to F. K. Guo, Q. Wang, Q. Zhao, and B. S. Zou for their insightful discussions and comments.