Prediction of a missing higher charmonium around 4.26 GeV in $J/\psi$ family

Inspired by the similarity of mass gaps of $J/\psi$ and $\Upsilon$ families, the prediction of missing higher charmonium with mass $4263$ MeV and very narrow width is made. In addition, the properties of two charmonium-like states, $X(3940)$ and $X(4160)$, and charmonium $\psi(4415)$ are discussed. Here, $X(3940)$ as $\eta_c(3S)$ is established while the explanation of $X(4160)$ to be $\eta_c(4S)$ is fully excluded and $\eta_c(4S)$ is typically a very narrow state. These predictions can be accessible at BESIII, Belle and BelleII in near future.

Inspired by the similarity between the mass gaps of the J/ψ and Υ families, the prediction of a missing higher charmonium with mass 4263 MeV and very narrow width is made. In addition, the properties of two charmonium-like states, X(3940) and X(4160), and charmonium ψ(4415) are discussed, where our calculation shows that X(3940) as η c (3S ) is established, while the explanation of X(4160) to be η c (4S ) is fully excluded and that η c (4S ) is typically a very narrow state. These predictions might be accessible at BESIII, Belle, and BelleII in near future.
PACS numbers: 13.25.Gv,12.38.Lg Since the observation of J/ψ in 1974 [1,2], the charmonium family has become abundant with more and more such states announced by the experiments [3]. Especially in the past decade, a series of charmonium-like states have been observed, which have further stimulated theorists' extensive interest in revealing their underlying properties (see a recent review in Ref. [4]), since these novel phenomena reflect nonperturbative behavior of quantum chromodynamics (QCD). Among the studies on these states, it is an important research topic for the whole community of particle physics how to identify the exotic states, whose establishment is, of course, tied with our understanding of the charmonium family. When checking the mass spectra of the observed charmonia with spin-parity J PC = 1 −− and comparing them with those of the corresponding bottomonia, we notice that the mass gap between ψ(2S ) and J/ψ is almost the same as that between Υ(2S ) and Υ(1S ). There also exist the similar mass differ- where ψ(2S ) and ψ(3S ) correspond to ψ(3686) and ψ(4040), respectively. However, if ψ(4415) is ψ(4S ), such a law is violated since the mass gap of ψ(4415) and ψ(3S ) is much larger than that of Υ(4S ) and Υ(3S ). In fact, the properties of the charmonia above 4.1 GeV are still not understood well, which is the possible reason to result in the above puzzling mass gap. In general, compared with the J/ψ family, the bottomonia with the radial quantum numbers n = 1, 2, 3, 4 were well established both by experiment and theory. Thus, the study of the J/ψ family can be borrowed from Υ family. If this law of mass gap relation still holds for states with n = 3, 4 in J/ψ and Υ families, we find that the mass of ψ(4S ) should be located at 4263 MeV, where we add the mass gap between Υ(4S ) and Υ(3S ) to the mass of ψ(3S ). In Fig. 1, we show the details of the mass gaps for J/ψ and Υ families and compare the corresponding gaps with each other. After obtaining the above prediction, we notice the results of the mass spectra of charmonium family given in Ref. [5], where the screening potential was adopted in their calculation. They also found that the mass of ψ(4S ) is about 4,273 MeV, which is consistent with our estimate of the mass for the missing charmonium ψ(4S ). Additionally, another theoretical work [6] also supports existence of a missing ψ(4S ). In Ref. [6], the color-screening effect was considered in the calculation of the mass spectrum of charmonium, where the mass of ψ(4S ) was obtained as 4,247 MeV.
If this predicted state exists in the J/ψ family, we must reveal its underlying properties to answer why there does not have any evidence in the present experiment, which will be the main task of this work.
Thus, both of Y(4260) and Y(4360) have J PC = 1 −− . We will come back later to discuss whether the missing charmonium has the relation to these two charmonium-like states.
In the following, we study the decay behavior of the missing state corresponding to ψ(4S ), which is crucial to explain why this state is still missing at the predicted location and how to search for this in future experiment. In the study we also obtain its full width which is the important information to us. Since the mass of the missing charmonium predicted (in the next discussion we adopt ψ(4S ) to denote this missing charmonium and assume its mass to be 4,263 MeV) is above the thresholds of D ( * )D( * ) and D ( * ) Here, we apply the quark pair creation (QPC) model [9] to calculate the decay widths, which is one of the most popular models for studying the Okubo-Zweig-Iizuka (OZI) allowed two-body strong decays of hadrons. For the OZI-allowed decay A → B + C, the transition matrix element in terms of the amplitude reads BC|T |A = δ 3 (P B + P C )M M J A M J B M J C in the center-of-mass frame of the initial state A, where P B and P C are the three-momenta of the final states B and C, respectively. In the above expression, T is the transition operator, i.e., where p 3 and p 4 denote the momenta carried by the quark and antiquark created from the vacuum. By using the Jacob-Wick formula [10], the partial wave amplitude can be expressed in terms of the amplitude as Finally the decay width is obtained as Γ A→BC = π 2 |P| J,L |M JL | 2 /m 2 A . The above description is just a short summary of the QPC model and the detailed introduction of the QPC model can be found in Refs. [11][12][13]. We need to specify that in the concrete calculation, the simple harmonic oscillator (SHO) wave function is adopted, where the parameter R is taken from Ref. [14] for final states. Moreover, the dimensionless parameter γ, describing the strength of the quarkantiquark pair creation, is taken as the same value as that in Ref. [14] by fitting the experimental data.
We first present the result of charmonium ψ(3S ) decaying into DD, DD * + H.c., D * D * and D sDs that are kinematically allowed by the phase space of ψ(3S ), which is shown in Fig.  2. Our calculated result can well reproduce the experimental widths of ψ(4040) in Ref. [15] when taking R = 1.56 − 1.63 GeV −1 . DD * + H.c. and D * D * are the dominant decay modes of ψ(3S ), while DD and D sDs are the subordinate decays. Here, the conclusion of the D ( * )D( * ) decay modes of ψ(3S ) is in agreement with the BaBar data [16] since the measured ratio Γ DD /Γ D * D +H.c. is 0.24±0.05±0.12. In the same reference, however, the ratio Γ D * D * /Γ D * D +H.c. is given by 0.18±0.14±0.03 [16] which contradicts with a larger value of our calculation as shown in Fig. 2. Hence this ratio must be tested further in future experiment. The above study shows that ψ(4040) assigned as ψ(3S ) is reasonable. Of course, the reliability of the QPC model is also tested here, which enables us to apply this model to safely study the decays of ψ(4S ).  Fig. 2, i.e., the total decay width of ψ(4S ) is stable over the corresponding R range adopted, while its partial decay widths strongly depend on the R value. This phenomenon is due to the node effects. Our result also shows that the node effects are important when discussing the higher charmonium decays because due to these effects we find that the predicted charmonium ψ(4S ) has very narrow width around 6 MeV.
In addition, we also discuss the dependence of ψ(4S ) decay behavior on the mass of ψ(4S ) in Fig. 3 [5,6]. The result listed in Fig. 3 indicates that the decay behavior of ψ(4S ) is weakly dependent on the mass of ψ(4S ). Thus, we come to the definite conclusion that ψ(4S ) is a narrow state.
As a typical higher charmonium with a very narrow width, the predicted ψ(4S ) is difficult to identify by the analysis of the open-charm decay channels [19][20][21][22] and the R value scan [23][24][25][26][27][28][29][30] based on the present experimental data, which can naturally answer why this higher charmonium is still missing in experiment. Thus, we expect future experimental results of the open-charm decays and more precise study of the R value scan, especially from BESIII, Belle, and forthcoming BelleII.
We notice a recent analysis of BESIII data in Ref. [31]. BESIII already realized the measurement of the cross sections of e + e − → h c (1P)π + π − at center-of-mass energies 3.90 − 4.42 GeV [32]. Yuan analyzed the data by fiting the line shape with two Breit-Wigner functions, which indicates that there are a narrow structure with mass 4, 216±18 MeV and width 39±22 MeV and another broad structure with mass 4, 293 ± 9 MeV and width 222 ± 67 MeV [31], where these two charmoniumlike structures have J PC = 1 −− . Comparing the resonance parameters of this narrow structure in Ref. [31] with our result of the predicted ψ(4S ), one finds their similarity, which means that this narrow charmonium-like structure can be as a good candidate of the predicted ψ(4S ) in this work. As mentioned above, more experimental efforts will be necessary to clarify this point.
In the following, we apply the above approach to study the η c family by comparing the mass gaps between the η c and J/ψ families, where the charmonia in η c and J/ψ families belong to S-wave cc states. There are two well established η c states listed in PDG, i.e., η c (1S ) and η c (2S ). Their mass gap is about 657 MeV, which is similar to that between ψ(2S ) and ψ(1S ). Thus, assuming the mass gap between η c (3S ) and η c (2S ) is the same as that between ψ(3S ) and ψ(2S ), we can estimate the mass of η c (3S ) as about 3,992 MeV, which is close to the mass of an observed charmonium-like state X(3940) in the D * D invariant mass spectrum of e + e − → D * D J/ψ [33,34].
Here, the C parity of X(3940) favors +1 since X(3940) is from the double charm production. η c also has the same C parity +1. Thus, we may identify X(3940) as η c (3S ). To test it, we also study the decay behavior of X(3940) under the η c (3S ) assignment, which is listed in Fig. 4. We find that the total decay width of η c (3S ) calculated can be well fitted with the experimental value of X(3940) when we take the same R range as that in the study of ψ(3S ), where DD * + H.c. is the only decay mode of η c (3S ), which can explain why X(3940) was observed in the DD * channel [33,34] under this assumption. In addition, another charmonium-like state X(4160) was reported by the Belle Collaboration through the double charm production, where X(4160) appears in the D * D * invariant mass distribution of e + e − → D * + D * − J/ψ [34]. In this work, we also discuss whether X(4160) can be explained as η c (4S ). In Fig. 4, the partial and total decay widths of η c (4S ) are listed. The situation of the decay behavior of η c (4S ) is very similar to that of ψ(4S ) discussed above. Our study shows that the total width of η c (4S ) is also very narrow, which is not consistent with the measured width of X(4160) that has 139 +111 −61 ± 21 MeV [34]. According to the above analysis, we can fully exclude the η c (4S ) assignment to X(4160).
Before summarizing our work, we would like to pay attention to ψ(4415). Here, we discuss the possibility of ψ(4415) as ψ(5S ).
The two-body strong decay behaviors of ψ(4415) as a ψ(5S ) state are presented in Fig. 5. The total width overlaps with the average width in PDG when R is in the range 19 GeV −1 , which well agree with the BaBar data [16] having the measured values of R 1 = 0.14 ± 0.12 ± 0.03 and R 2 = 0.17 ± 0.25 ± 0.03 [16]. The above analysis indicates that it is very likely that ψ(4415) is the ψ(5S ) state. However, in our calculation the branching ratio Γ DD * 2 /Γ total is about two orders of magnitude smaller than the experimental data from Belle [21]. This incompatibility is a challenge to the explanation of ψ(4415) as ψ(5S ) state. We also suggest the experimental study of ψ(4415) in future experiments, which is important to further understand the properties of ψ(4415).
For the convenience of reader, in Table I we further collect the experimental information of these charmonium-like states or charmonium involved in the present work. In summary, in this work we have predicted a missing charmonium ψ(4S ) with quantum number n 2S +1 L J = 4 3 S 1 , which has mass around 4263 MeV and very narrow width. This observation has been obtained through the similarity of the mass gaps existing in J/ψ and Υ families and further study of its OZI-allowed decay behavior. Comparing this state with the reported higher charmonia ψ(4040), ψ(4160) and ψ(4415), the predicted ψ(4S ) state is the first higher charmonium with such a narrow width, which can explain why there is no evidence in the corresponding analysis of the open-charm decay channels and the R value scan until now. The relation of ψ(4S ) and the recent evidence of a narrow charmonium-like structure in the e + e − → h c (1P)π + π − process [31] has also been discussed.
We have also explored two charmonium-like states X(3940) and X(4160) under the η c (3S ) and η c (4S ) assignments, respectively. Our study has indicated that X(3940) can be well explained as η c (3S ) but X(4160) as η c (4S ) is fully excluded. Here, we have predicted η c (4S ) also has very narrow width similar to the situation of the discussed ψ(4S ). In addition, the properties of ψ(4415) as ψ(5S ) have been given in this work.
The study presented in this work can enhance our understanding of charmonium family especially for J/ψ family, which is valuable to reveal the underlying QCD nonperturbative effects. To test our predictions, we expect further experimental progress on higher charmonia, where BE-SIII, Belle, and BelleII will be good platforms to carry out the search for our predictions.