arxiv: v3 [hep-ph] 9 May 2017

Transcription

1 The analysis of the charmonium-like states X (386),X(387), X(3915), X(393) and X(394) according to its strong decay behaviors. Guo-Liang Yu 1, Zhi-Gang Wang 1, and Zhen-Yu Li 1 Department of Mathematics and Physics, North China Electric power university, Baoding 713, People s Republic of China School of Physics and Electronic Science, Guizhou Normal College, arxiv: v3 [hep-ph] 9 May 17 Guiyang 5518, People s Republic of China (Dated: May 1, 17) Inspired by the newly observed state X (386), we analyze the strong decay behaviors of some charmonium-like states X (386),X(387), X(3915), X(393) and X(394) by the 3 P model. We carry out our work based on the hypothesis that these states are all being the charmonium systems. Our analysis indicates that ++ charmonium state can be a good candidate for X (386) and 1 ++ state is the possible assignment for X(387). Considering as the 3 1 S state, the decay behavior of X(394) is inconsistent with the experimental data. So, we can not assign X(394) as the 3 1 S charmonium state by present work. Besides, our analysis implies that it is reasonable to assign X(3915) and X(393) to be the same state, ++. However, combining our analysis with that of Zhou [14], we speculate that X(3915)/X(393) might not be a pure cc systems. PACS numbers: 13.5.Ft; 14.4.Lb 1 Introduction Very recently, the Belle Collaboration observed a new charmonium-like state, X (386), by performing a full amplitude analysis of the process e + e J/ψDD [1]. Its mass is ( MeV/c ) and width is ( )MeV. The JPC = ++ hypothesis is favored over the ++ assignment at the level of.5σ. In reference [], this state was explained to be a Cγ 5 γ5 C type scalar tetraquark state by the method of QCD sum rules(qcdsr). Actually, people once had assigned X(3915) as the ++ charmnium state [3, 4], which was observed by Belle, BABAR Collaboration in B J/ψωK decay mode [4 8]. Its mass and width are listed in Table I. After X(3915) was suggested to be the χ c assignment, it encountered several challenges [9 1]. For example, the decayχ c (P) DD, which was expected to be the dominant decaymode, has not been observed experimentally. In contrast, the decay mode X(3915) Jψω, which should be OZI(Okubo- Zweig-Iizuka) [13] suppressed, was observed instead in experiments. In addition, the mass splitting of Electronic address: yuguoliang11@163.com Electronic address: zgwang@aliyun.com

2 χ c (P) χ c (P) is too small. A reanalysis of the data from Ref. [4], presented in Ref. [14], showed that X(3915) could be the same state as X(393), whose quantum number is ++ [17], due to the degeneracies of their masses and widths. Now, the observation of X (386),which was assigned to be the ++ state, is an important verification of the results in Ref. [14]. ThemassofthisnewlyobservedX (386)isclosetothatofanothercharmonium-likestateX(387). However, these two hadrons are impossible to be the same state because of its different decay modes and widths(see Table I). After X(387) was discovered by Belle Collaboration [18] and confirmed by BABAR[19], CDF[], D[1] and Bell[] Collaborations, its nature has still been very controversial. It was mainly explained to be such structures as a molecule state [4 33], a hybrid charmonium [34 36], a tetraquark state [37 41]. Another important explanation is that it was the charmonium state with quantum of 1 ++ [4, 43], which has a dominant decay mode D D. Belle Collaboration reported another charmonium-like state X(394) from the inclusive process e + x J/ψ+cc at a mass of M = (3.943 ±.6 ±.6)GeV/c [44]. Later, its decay width was confirmed to be Γ = ( ±8)MeV [45]. People have also explored the structure of X(394) with different kinds of methods such as the light-cone formalism [46], the NRQCD factorization formula [47, 48] and QCDSR [49, 5]. According to these studies, there seems to be no doubt that the quantum number of X(394) is 3 1 S. However, its structure is still controversial, which have been explained to be different states such as the charmonium state [46], the molecular state [49 51] and a Mixed Charmonium-Molecule State [5, 53]. TABLE I: The experimental information about the X states in this paper. States Mass(MeV/c ) Width(MeV) J PC Decay channels X (386) [1] ( 3 P ) DD X(3915) [4] ±.±1.6 13±6±3 ++ ( 3 P ), ++ ( 3 P ) J/ψω X(393) 399±5± [15] 9±1± ++ ( 3 P ) DD ±.7±1.1 [16] 1.3±6.8± ( 3 P ) DD 387±.6±.5 [18] J/ψπ + π ±.7±.4 [] J/ψπ + π X(387) ±3.1±3. [1] < ( 3 P 1) J/ψπ + π ±1.4 [19] J/ψπ + π ± [] D D π,j/ψω ± [3] D D +h.c.,j/ψρ X(394) [44, 45] ± ±8 + (3 1 S ) DD In summary, these newly discovered charmonium-like states have inspired many interests about their phyisical natures. In order to further study its structures, we perform an analysis of the strong decay behaviors of X (386), X(387), X(3915), X(393) and X(394) with the 3 P decay model.

3 3 The experimental information about these states are listed in Table I. Since these X states can not be completely ruled out from the cc systems at present, we carry out our calculations by assuming them to be the charmoniums. Our analysis will be helpful to confirm or exclude some cc systems and useful to further determine the quantum numbers of the confirmed charmonium states. As for the strong decays of the hadrons, 3 P decay model [54 56] is an effective method. It has been widely used in this field since it gives a good description of the decay behaviors of many hadrons [57 64]. The article is arranged as follows: In section, we give a brief review of the 3 P decay model; in Sec.3 we study the strong decays of X (386), X(387), X(3915), X(393) and X(394); in Sec.4, we present our conclusions. The decay model 1 1 B B A 4 A 4 3 (a) C (b) 3 C FIG. 1: The two possible diagrams contributing to A BC in the 3 P model. The principle of 3 P decay model is illustrated clearly in Fig.1, where a quark-antiquark pair(q 3 q 4 ) is created from the vacuum with ++ quantum numbers. With the q 1 q within the initial meson, this quark systems regroups into two outgoing mesons via quark rearrangement for the meson decay process A BC. Its transition operator in the nonrelativistic limit reads T = 3γ 1m1 m d 3 p 3 d 3 p 4 δ 3 ( p 3 + p 4 )Y1 m ( p 3 p 4 )χ 34 1 mϕ 34 ω 34 b 3 ( p 3)d 4 ( p 4) (1) m where γ is a dimensionless parameter reflecting the creation strength of the quark-antiquarkq 3 q 4 pair. The solid harmonic polynomial Y m 1 ( p) p 1 Y m 1 (θ p,φ p ) reflects the momentum-space distribution of the q 3 q 4. The helicity amplitude of the decay process A BC in the parent meson A center of mass frame

4 4 is M MJ A MJ B MJ C ( P) =γ 8EA E B E C M LA,M SA, M LB,M SB, L A M LA S A M SA J A M JA L B M LB S B M SB J B M JB M LC,M SC,m L C M LC S C M SC J C M JC 1m1 m χ 14 S BM SB χ 3 S CM SC χ 1 S AM SA χ 34 1 m [ φ 14 B φ 3 C φ 1 Aφ 34 I( P,m 1,m,m 3 ) +( 1) 1+SA+SB+SC φ 3 B φ 14 C φ 1 A φ 34 I( P,m,m 1,m 3 )], where I( P,m 1,m,m 3 ) is the spatial integral which is defined as I( P,m 1,m,m 3 ) = d 3 pψ n m 3 BL BM LB ( PB + p)ψ m 3 n m 1 +m CL CM LC ( PB + p) m +m 3 (3) ψ nal AM LA ( P B + p)y 1 m ( p) where P = P B = P C, p = p 3, m 3 is the mass of the created quark q 3. We employ the simple harmonic oscillator (SHO) approximation as the meson space wave functions in Eq.(3). Ψ nlml ( p) =( 1) n ( i) L n! p R L+3 Γ(n+L+ 3 )exp( R )LL+1 n (R p )Y LML ( p) (4) Where R is the scale parameter of the SHO. With the Jacob-Wick formula, the helicity amplitude can be converted into the partial wave amplitude M JL ( P) 4π(L+1) = LJM JA J A M JA J B M JB J C M JC JM JA M MJAMJBMJC ( P) J A +1 (5) M JB M JC where M JA = M JB +M JC, J A = J B +J C and J A +J P = J B +J C +J L. Finally, the decay width in terms of partial wave amplitudes is where P = P = A, B, and C, respectively. Γ = π P M JL (6) 4 M A [M A (M B+M C) ][M A (MB MC) ] M A 3 The results and discussions JL (), M A, M B, and M C are the masses of the meson The decay width based on 3 P model depends on the following input parameters, the light quark pair(qq) creation strength γ, the SHO wave function scale parameter R, and the masses of the mesons and the constituent quarks. The adopted masses of the hadrons are listed in TABLE II, and m u = m d =. GeV, m s =.419 GeV and m c = 1.65 GeV [65]. As for the scale parameter R, there are mainly two kinds of choices which are the common value and the effective value. The effective value can be fixed to reproduce the realistic root mean square radius by solving the Schrodinger equation with the linear potential [66, 67]. For the cc systems, the R value of P states is estimated to be

5 5 TABLE II: The adopted masses of the hadrons used in our calculations. States M X (386) M X(387) M X(3915) M X(393) M X(394) Mass(MeV) 386 [1] 387 [18] 3919 [4] 397 [16] 394 [44, 45] States M D ± M D M D ± M D Mass(MeV) [65] GeV 1 [68]. For the mesons D andd, its value is taken to be R D [D ± ] = 1.5GeV 1, R D [D ± ] = 1.85GeV 1 [66, 68] in this work. Finally, we choose the value of γ to be 6.5 for the creation of the u/d quark following Ref. [57]. We know that X (386) was favored to be the ++ ( 3 P ) charmonium-like state by Bell Collaboration and X(3915) had also once been explained to be this assignment. Lately, the latter one was corrected to be the same state as another charmonium-like state, X(393) which had been determined to be ++ ( 3 P ) assignment. In order to further confirm these conclusions, we study the strong decay behaviors of X (386) by considering it as the 3 P and 3 P charmoniums. And so does for the X(3915) state. Besides, we also perform an analysis of the decay behaviors of X(393), X(387) and X(394) which have been favored to be ++ ( 3 P ), 1 ++ ( 3 P 1 ) and + (3 1 S ) states, respectively. As mentioned in Ref. [9], the mass difference M X(393) M X(3915) = 9.7±3.7 MeV, is smaller than the fine splitting of 1P states M χc M χc = ±.3MeV [17]. This is an important evidence to recognizing X(3915) and X(393) as the same state. In order to determine its mass precisely, we also calculate the decay widths of 3 P (χ c ) state on different masses. All of the results are illustrated in the form of graphs, which can be seen from Figures to D D D + D D D D + D 5 R X*(386) (GeV 1 ) R X*(386) (GeV 1 ) FIG. : The strong decay of X (386) as the ++ ( 3 P ) state on scale parameter R X (386). FIG. 3: The strong decay of X (386) as the ++ ( 3 P ) state on scale parameter R X (386). Whether we consider X (386) as ++ or ++ charmonium state, there is only one strong decay mode, X(386) DD, where D refers to either D or D +. From Figures and 3, we can clearly

6 D D D + D 15 1 D D D + D D* D D* + D R (GeV 1 ) X(3915) R (GeV 1 ) X(3915) FIG. 4: The strong decay of X(3915) as the ++ ( 3 P ) state on scale parameter R X(3915). FIG. 5: The strong decay of X(3915) as the ++ ( 3 P ) state on scale parameter R X(3915) D D D + D D * D D *+ D DD * DD R (GeV 1 ) X(393) M χ (P) (GeV) FIG. 6: The strong decay of X(397) as the ++ ( 3 P ) state on scale parameter R X(397). FIG. 7: The relation of the strong decay width with the mass of χ c( 3 P ). see the deference between the total decay widths of these two states. Taking R =.3.5GeV 1 discussed above, the total strong decay width of ++ state ranges from 11 to 18MeV, which is compatible with the experimental data in Ref. [1]. The total decay width of ++ state, which ranges from.4 1.9MeV, is much smaller than the experimental data. That means, if we assume X (386) as the ++ charmonium state, its dominant decay mode and total decay width is consistent well with the experimental data. Thus, our present work support X (386) as the ++ (χ c ) charmonium state. Considering X(3915) as the ++ and ++ charmonium respectively, we also observe different strong decay behaviorsfrom Figures4 and 5. For the ++ state, the total strong decay width rangesfrom 159 to MeV, which dominantly decays into DD. Not only its total decay width but also the dominant decay channel is inconsistent with the experimental data in Ref. [4](See Table I). This means that the X(3915) was assumed to be the charmonium state ++ is disfavored. If it is treated as the ++ charmonium, its decay behavior is very similar with that of X(393)(See Figures 5 and 6). They both decay into DD and DD with the total decay width ranging from 1. to 3.MeV. In addition, these

7 values of the decay widths fall in the range of the experimental data. Thus, it seems reasonable to assign both X(3915)and X(393)to be the the ++ (χ c ) charmonium state. If this conclusion is true, the mass of the χ c charmonium state has erros. So, we plot the relations of the strong decay widths on the masses of χ c ( ++ ) in Figure 7, which will be helpful to determine its mass in the experimental and theoretical explorations in the future. Since the decay width of χ c DD is larger than 1MeV, it should be observable in experiments for both X(393) and X(3915). However, it was reported by both Bell and Babar Collaborions that the X(393) and X(3915) were observed in two different decay channels, X(393) DD and X(3915) J/ψω. A reanalysis presented in Ref. [14] shows that if helicity- dominance assumption is abandoned and a sizable helicity- component is allowed, the decay process X(3915) DD may be reproduced in the experimental data. But the large helicity- contribution means that X(393)/X(3915) might not be a pure cc charmonium state. Since X(387) was observed, there have accumulated abundant experimental information, which can be seen in TableI. Belle experiment indicated B(X(387) D D π K + ) = B(X(387) J/ψπ + π K + ) []. Based on these experimental data, we can draw a conclusion that D D is the dominant decay of X(387). Although the underlying structure of this state is very controversial, there is no doubt that its quantum number is As a charmonium state χ c1 (1 ++ ), we show the dependence of the strong decay width on the scale parameter R in Figure 8. Taking R =.3.5GeV 1, the decay width of the inclusive decay channel D D ranges from. to 1.MeV, which falls in the range of the experimental data in Table I and is also consistent with the conclusion of D D being the dominant decay mode. Thus, our present work implies that X(387) is assigned to be the χ c1 (1 ++ ) charmonium state is reasonable D* D 5 4 D * D D *+ D R (GeV 1 ) X(387) R (GeV 1 ) X(394) FIG. 8: The strong decay of X(387) as the 1 ++ ( 3 P 1) state on scale parameter R X(387). FIG. 9: The strong decay of X(394) as the + (3 1 S ) state on scale parameter R X(394). Finally, we can see in Figure 9 that, as a 3 1 S charmonium state, X(394) can decay into D D and D + D final states. This result is consistent with the experiments, where X(394) was truly

8 observed from the inclusive process e + e J/ψD D. However, it can also be seen from Figure 9 that the maximum of the total decay width can only reach up to 1MeV if R is changed from. to 3.GeV 1. The predicted decay width in experiments is Γ = ±8MeV which is much larger than this value. This comparison indicates that 3 1 S charmonium state might not be a good candidate for the X(394). 4 Conclusion In summary, by considering both X (386) and X(3915) as ++ and ++ charmonium states, X(387), X(393), X(394) as 1 ++, ++ and + charmonium separately, we study its two-body open charm strong decay behaviors by the 3 P decay model. According to comparing our results with the experimental data, we find that X (386) and X(387) can be explained to be the χ c ( 3 P ) and χ c1 ( 3 P 1 ) charmonium state separately. The decay width of X(394) is inconsistent with the experimental data if it is supposed to be a 3 1 S charmonium state. Thus, 3 1 S charmonium state can be ruled out at prsent as a candidate for X(394). Treated as a ++ charmonium, the decay behavior of X(3915) is contradictory to experimental data. This indicates that X(3915) is unlikely to be a ++ charmonium state. Supposed as a ++ charmonium, the decay behavior of X(3915) is consistent with not only the experimental data by also that of X(393). Thus, we tentatively assign these two states as the same charmonium χ c. According to a reanalysis of the experimental data, Zhou [14] also suggested them to be the same state ++, but with a significant non-cc component. As a result, the structure of X(3915)/X(393) needs to be further studied according to more experimental and theoretical explorations. Acknowledgment This work has been supported by the Fundamental Research Funds for the Central Universities, Grant Number 16M S [1] K. Chilikin(Belle Collaboration), arxiv: v1 [hep-ex](17). [] Zhi Gang Wang, arxiv: [hep-ph](17). [3] X. Liu, Z.-G. Luo, and Z.-F. Sun, Phys. Rev. Lett. 14, 11 (1). [4] J.P. Lees et al. (BABAR Collaboration), Phys. Rev. D 86, 7 (1). [5] K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 94, 18 (5). [6] S.-K. Choi et al. (Belle Collaboration), Phys. Rev. Lett. 94, 18 (5). [7] P.del Amo Sanchez et al. (BABAR Collaboration), Phys. Rev. D 8, 1111R (1); B. Aubert et al. (BABAR Collaboration), Phys. Rev. Lett. 11, 81 (8). [8] S. Uehara et al. (Belle Collaboration), Phys. Rev. Lett. 14, 91 (1). [9] F.-K. Guo and Ulf-G. Meissner, Phys. Rev. D 86, 9151 (1).F.-K. Guo, C. Hanhart, G. Li, Ulf-G. Meiβer, and Q. Zhao, Phys. Rev. D 83, 3413 (11).

9 9 [1] S. L. Olsen, Phys. Rev. D 91, 5751 (15). [11] E. J. Eichten, K. Lane, and C. Quigg, Phys. Rev. D 73, 1414 (6). [1] Z.G. Wang, Eur. Phys. J. C 77, 78(17); Eur. Phys. J. A 53, 19(17). [13] S. Okubo, Phys. Lett. 5, 165 (1963); G. Zweig, CERN Report Th 41 and 41 (1964); J. Iizuka, K. Okada and O. Shito, Prog. Theor. Phys. 35, 161 (1966); J. Iizuka, Prog. Theor. Phys. Suppl. 37, 1 (1966). [14] Z. Y. Zhou, Z. Xiao and H. Q. Zhou, Phys. Rev. Lett. 115, 1 (15). [15] S.Uehara et al.(bell Collaboration),Phys.Rev.Lett.96,83(6). [16] B.Aubert et al.(babar Collaboration),Phys.Rev.D81,93(1). [17] J. Beringer et al. (Particle Data Group), Phys. Rev. D 86, 11 (1). [18] S.K. Choi et al., Belle Collaboraion, Phys. Rev. Lett. 91, 61 (3). [19] B. Aubert et al., Babar Collaboration, Phys. Rev. D 71, 7113 (5). [] D. Acosta et al., CDFII Collaboration, Phys. Rev. Lett. 93, 71 (4). [1] V.M. Abazov et al., D Collaboration, Phys. Rev. Lett. 93, 16 (4). [] G. Gokhroo et al., Belle Collaboration, Phys. Rev. Lett. 97, 16 (6). [3] Babar Collaboration, talk given by P. Grenier in Moriond QCD 7, 17-4 March, [4] F.E. Close, P.R. Page, Phys. Lett. B 578, 119 (4). [5] M.B. Voloshin, Phys. Lett. B 579, 316 (4). [6] C.Y. Wong, Phys. Rev. C 69, 55 (4). [7] E.S. Swanson, Phys. Lett. B 588, 189 (4); ibid B 598, 197 (4). [8] N.A. Tornqvist, Phys. Lett. B 59, 9 (4). [9] Mohammad T. AlFiky, Fabrizio Gabbiani, Alexey A. Petrov, Phys.Lett. B 64, 38-45(6). [3] A.B. Larionov, M. Strikman, M. Bleicher, Phys. Lett. B 749, (15). [31] Z.G.Wang, T. Huang, Eur.Phys.J.C 74,891 (14). [3] A. Esposito, A. Pilloni, and A.D. Polosa, Phys. Rept. 668, 1 (17). [33] Xian-Wei Kang and J. A. Oller, arxiv:161.84[hep-ph](16). [34] B.A. Li, Phys. Lett. B 65, 36 (5). [35] M. Nielsen, C.M. Zanetti, Phys.Rev.D 8, 116 (1). [36] M. Takizawa, S. Takeuchi, arxiv: [hep-ph]. [37] Y. Cui, X.L. Chen, W.Z. Deng, S.L. Zhu, High Energy Phys. Nucl. Phys. 31, 7 (7). [38] R.D. Matheus, S. Narison, M. Nielsen, J.M. Richard, Phys. Rev. D 75, 145 (7). [39] T.W. Chiu, T.H. Hsieh, Phys. Lett. B 646, 95 (7). [4] S.Dubnicka, A.Z.Dubnickova, M.A.Ivanov, and J.G.Korner, Phys.Rev.D 81,1147(1). [41] Z.G.Wang, T. Huang,Phys. Rev. D 89, 5419 (14). [4] N. N. Achasov and E. V. Rogozina, Mod. Phys. Lett. A 3, (15). [43] Wei-Jun Deng, Hui Liu, Long-Cheng Gui, and Xian-Hui Zhong, Phys. Rev. D 95, 346 (17). [44] K.Abe (Bell Collaboration)Phys. Rev. Lett. 98, 81 (7). [45] P.Pakhlov, (Bell Collaboration)Phys. Rev. Lett. 1, 1 (8). [46] V. V. Braguta, A. K. Likhoded, A. V. Luchinsky, Phys. Rev. D 74, 944 (6). [47] R.L.Zhu, Phys. Rev. D 9, 7417 (15).

10 1 [48] Z.G.He, B.Q.Li, Phys.Lett.B 693, 36-43(1). [49] R. M. Albuquerque, M. E. Bracco and M. Nielsen, Phys. Lett. B 678, 186 (9) [arxiv: [hepph]]. [5] Z.G.Wang Eur.Phys.J.C 74, 963 (14). [51] X. Liu and S. -L. Zhu, Phys. Rev. D 8, 175 (9). [5] R.M. Albuquerque, J.M. Dias, M. Nielsen, C.M. Zanetti, Phys. Rev. D 89, 767 (14). [53] F. Fernandez, P. G. Ortega, D. R. Entem, AIP Conf. Proc. 166, 168 (14). [54] L. Micu, Nucl. Phys. B1, 51 (1969). [55] R. Carlitz and M. Kislinger, Phys. Rev. D, 336 (197); E.W. Colglazier and J. L. Rosner, Nucl. Phys. B7, 349 (1971); W. P. Petersen and J. L. Rosner, Phys. Rev. D 6, 8 (197). [56] A. Le Yaouanc, L. Oliver, O. Pene, and J.-C. Raynal, Phys. Rev. D 8, 3 (1973); 9, 1415 (1974); 11, 17 (1975); Phys. Lett. B 71, 397 (1977);A. Le Yaouanc, L. Oliver, O. Pene, and J. C. Raynal, Phys. Lett. B 7, 57 (1977). [57] H. G. Blundell, arxiv:hep-ph/968473; H. G. Blundell and S. Godfrey, Phys. Rev. D 53, 37 (1996); H. G. Blundell, S. Godfrey, and B. Phelps, Phys. Rev. D 53, 371 (1996). [58] H. Q. Zhou, R. G. Ping, and B. S. Zou, Phys. Lett. B 611, 13 (5). [59] D.-M. Li and S. Zhou, Phys. Rev. D 78, 5413 (8); D.-M. Li and E. Wang, Eur. Phys. J. C 63, 97 (9); D.-M. Li, P.-F. Ji, and B. Ma, Eur. Phys. J. C 71, 158 (11). [6] B. Zhang, X. Liu,W. Z. Deng, and S. L. Zhu, Eur. Phys. J. C 5, 617 (7); Y. Sun, Q. T. Song, D. Y. Chen, X. Liu, and S. L. Zhu, Phys. Rev. D 89, 546 (14). [61] E. S. Ackleh, T. Barnes, and E. S. Swanson, Phys. Rev. D 54, 6811 (1996); T. Barnes, N. Black, and P. R. Page, Phys. Rev. D 68, 5414 (3); T. Barnes, S. Godfrey, and E. S. Swanson, Phys. Rev. D 7, 546 (5). [6] F. E. Close and E. S. Swanson, Phys. Rev. D 7, 944 (5); F. E. Close, C. E. Thomas, O. Lakhina, and E. S. Swanson, Phys. Lett. B 647, 159 (7). [63] J. Ferretti, G. Galata, and E. Santopinto, Phys. Rev. C 88, 157 (13); J. Ferretti, G. Galata, and E. Santopinto, Phys. Rev. D 9, 541 (14); J. Ferretti and E. Santopinto, Phys. Rev. D 9, 94 (14). [64] YU Guo-Liang, WANG Zhi-Gang, LI Zhen-Yu, MENG Gao-Qing, Chin. Phys. C 39, 6311 (15); Guo-Liang Yu, Zhi-Gang Wang, Zhen-Yu Li, Phys. Rev. D 94, 744 (16). [65] C. Patrignani et al. (Particle Data Group), Chin. Phys. C 4, 11 (16). [66] S. Godfrey and R. Kokoski, Phys. Rev. D 43, 1679 (1991). [67] B-Q. Li, and K-T. Chao, Phys. Rev. D 79, 944 (9). [68] You-chang Yang, Zurong Xia, and Jialun Ping, Phys. Rev. D 81, 943 (1).