arxiv: v1 [hep-ph] 14 Apr 2019

Transcription

1 Production of the D 0 (2317) and D 1 (2460) by kaon-induced reaction on a roton target arxiv: v1 [he-h] 14 Ar 2019 HongQiang Zhu 1 and Yin Huang 2, 1 College of Phyic and Electronic Engineering, Chongqing Normal Univerity, Chongqing ,China 2 School of Phyic and Nuclear Energy Engineering, Beihang Univerity, Beijing , China (Dated: Aril 16, 2019) We invetigate the oibility to tudy the charmed-trange meon D 0 (2317) and D 1 (2460) by kaon-induced reaction on a roton target in an effective Lagrangian aroach. The roduction roce i decribed by the t-channel D 0 and D 0 exchange, reectively. Our theoretical aroach i baed on the chiral unitary theory where the D 0 (2317) and D 1 (2460) reonance are dynamically generated. Within the couling contant of the D 0 (2317) to KD and D 1 (2460) to KD channel obtained from the chiral unitary theory, the total and differential cro ection of the Λ c D 0 (2317) and Λ c D 1 (2460) are evaluated. The K initial tate interaction mediated by Pomeron and Reggeon exchange i alo included, which reduce the roduction cro ection of the D 0 (2317) and D 1 (2460). If meaured in future exeriment, the redicted total cro ection and ecific feature of the angular ditribution can be ued to tet the (molecular) nature of the D 0 (2317) and D 1 (2460). PACS number: Le, Mk,13.25.Jx I. INTRODUCTION In recent year, many charmed-trange meon have been oberved [1]. Among them, the D 0 (2317) and D 1 (2460) are two eculiar tate(we abbreviate it a D 0 and D 1 hereafter), ince their mae are about 160 MeV and 70 MeV, reectively, below the quark model redicted value [2]. The charmed-trange meon D 0 wa firt oberved by the BaBar collaboration a a narrow eak in the D π invariant ma ditribution [3]. The tate wa hortly confirmed by CLEO [4] and Belle [5, 6] collaboration. Nowaday it i well etablihed in the PDG [1] with quantum number I(J P )=0(0 + ). The D 1 wa alo oberved in the CLEO exeriment [4] in the D π channel and BaBar [7 9] alo found a ignal in that region. Nowaday it i alo well etablihed in the PDG [1] with quantum number I(J P )=0(1 + ). The ma and width of the D 0 and D 1 tate reorted by the above collaboration [3 9] are conitent with each other, i.e., D 0 (2317) ± : D 1 (2460) ± : M=2317.7±1.3 MeV, Γ <3.8 MeV; M=2459.5±0.6 MeV, Γ <3.5 MeV. The large diagreement between the quark model exectation [2] and the exerimental meaurement [3 9] made it difficult to aign thee two tate a conventional charmedtrange meon. Since the mae of the D 0 and D 1 are about 40 MeV below the DK and D K threhold, reectively, many tudie rooed that the D 0 and D 1 are S -wave DK and D K molecular tate. The tudie in the Bethe-Saleter aroach [10] and otential model [11] howed indeed that the D 0 could be a DK hadronic molecule. In Ref. [12], the correonding author Electronic addre: huangy2017@buaa.edu.cn D 0 and D 1 were conidered a kaonic molecule bound by trong hort range attraction. Auming that the D 0 and D 1 are DK and D K molecular tate, the trong and radiative decay of the D 0 and D 1 were tudied by everal grou [13 16]. The roduction of the D 0 and D 1 from the nonletonic B decay were calculated in Ref. [17], in which the D 0 and D 1 were alo conidered a hadronic molecule of DK and D K, reectively. In the chiral unitary aroach [18 23], the D 0 and D 1 can be dynamically generated from the DK/D K and couled channel interaction. In addition to the interretation of the D 0 and D 1 a DK and D K molecule, the oibility to aign them a a conventional oen charmed meon wa alo dicued in many different aroache, uch a the relativitic quark model [24], the chiral erturbation theory [25], the quark air-creation model [26, 27], and the QCD um rule [28 32]. On the other hand, the large-n c exanion indicated that the D 0 could be a tetraquark meon [33]. The tetraquark interretation wa alo rooed to undertand the ma and decay behavior of the D 0 [34]. We note that the QCD um rule alo uorted the idea that the D 0 doe not eem to be a tandard quarkantiquark meon [35]. The reent knowledge about the D 0 and D 1 wa obtained from the e + e colliion [3 9]. Thu, it will be helful to undertand the nature of the D 0 and D 1 if we can oberve them in other roduction rocee. High-energy kaon beam are available at OKA@U-70 [36] and SPS@CERN [37], which rovide another alternative to tudy D 0 and D 1. The kaon beam at J-PARC can alo be ugraded to the energy region required in charmed-trange meon roduction [38]. Therefore, it i intereting to tudy the D 0 and D 1 roduction in the Λ c D 0 and K Λ c D 1 reaction. Since there exit lenty of exerimental information about the K elatic interaction in the energy region relevant to the D 0 and D 1 roduction [3 9], the effect from the K initial tate interaction (ISI) can be taken into account in order to make a more reliable rediction. Thi aer i organized a follow. In Sec. II, we will reent the theoretical formalim. In Sec. III, the numerical

2 2 reult of the kaon-induced D 0 and D 1 roduction on a roton target will be given, followed by dicuion and concluion in the lat ection. II. THEORETICAL FORMALISM The tree level Feynman diagram for the Λ c D 0 and Λ c D 1 reaction are deicted in Fig. 1, where the t channel D and D exchange are conidered. In thi work, the contribution from and u channel are ignored, becaue the and u channel, which involve the creation of an additional quark air creation in the kaoninduced roduction, are uually trongly ureed, Hence, the Λ c D 0 and K Λ c D 1 reaction hould be dominated by the Born term through the t-channel D and D exchange, which make the background very mall. 1 2 (a) D 3 4 D Λ + c (b) D 4 D 1 Λ + c ioin I= 0, found to be g KDD = in Ref. [18], multilied by the aroriate Clebch-Gordan (CG) coefficient, namely g D 0 D = g KDD / 2. Similarly to Ref. [18], we rely on the chiral unitary aroach [19] to obtain the couling contant g D 0 D = 9.82/ 2. 1 When evaluating the cattering amlitude of the Λ c D and Λ c D 1 reaction, we need to include form factor becaue hadron are not oint-like article. We adot here a common cheme ued in many reviou work [41, 42], F D ( )(q 2 D ( ), M ex )= Λ2 D M 2 D Λ 2 D q 2 D, (5) for the t-channel D ( ) meon exchange. Here the q D ( ) and M D ( ) are the four momentum and the ma of the exchanged D ( ) meon, reectively. In thi model, theλ D ( ) i the hard cutoff and it can be directly related to the hadron ize. Emirically, the cutoff arameterλ D ( ) hould be at leat a few hundred MeV larger than the D ( ) ma. Hence, we choe Λ D ( ) = M D ( )+αλ QCD withλ QCD = 0.22 GeV a ued in reviou work [43 45] for other reaction. The arameter α i uually cloe to unitary, and in thi work a variation of the cutoff will be made to how the enitivity of the reult on the cutoff. FIG. 1: Feynman diagram for the mechanim of the D 0 and D 1 roduction in the D 0 Λ c and D 1 Λ c reaction. We alo how the definition of the kinematic ( 1, 2, 3,and 4) ued in the calculation. A. Lagrangian To comute the diagram hown in Fig. 1, we need the effective Lagrangian denitie for the relevant interaction vertice. For theλ c D andλ c D vertice, we adot the commonly emloyed Lagrangian denitie a follow [39], L Λc D= ig Λc D Λ c γ 5 D 0 + H.c., (1) L Λc D = g Λ c D Λ c γ µ D 0 µ + H.c.. (2) The couling contant g Λc D = and g Λc D =-5.20 are determined from the SU(4) invariant Lagrangian [40] in term of g πnn =13.45 and g ρnn =6.0. In addition to theλ c D andλ c D vertice, we alo need the information on the KDD and KD D 1 vertice. A mentioned in the chiral unitary aroach of Ref. [18 23], the D and D 1 reonance are identified a wave meon-meon molecule that include big KD and KD comonent, reectively. We can write down the KDD and KD D 1 vertice of Fig. 1 a L KDD 0 = g KDD 0 KDD 0, (3) = g KD D KD µ D 1 1,µ, (4) L KD D 1 where the couling of the D 0 to D0, g KDD 0, i obtained from the couling contant of the D 0 to the DK channel in B. Initial tate interaction(isi) Following Ref. [46], the initial tate interaction for the reaction at high energie will be taken into account and the relevant Feynman diagram for the Initial tate interaction(isi) i hown in Fig. 2. The amlitudet K i P,f 2,ρ,ω,a 2 FIG. 2: Feynman diagram for the mechanim of the initial tate interaction of the ytem. written in term of the Pomeron and Reggeon exchange [46] T K = A IP + A f2 + A a2 + A ω + A ρ. (6) When the center-of-ma energy i large, the elatic cattering amlitude i a um of the following term, A i (, t)=η i Ci KN ( ) αi(t) 1 ex( Bi KN 0 2 t), (7) where i = IP for Pomeron and f 2, a 2,ω, andρreggeon. The energy cale 0 = 1 GeV 2. The couling contant Ci KN, the arameter of the Regge linear trajectorie[α i (t)=α i (0)+ α i t], the ignature factor(η i), and the B i KN ued in Ref. [46] rovide a rather good decrition of the exerimental data. The arameter determined in Ref. [46] are lited in Table. I.

3 3 TABLE I: Parameter of Pomeron and Reggeon exchange determined from elatic and total cro ection in Ref. [46]. i η i α i (t) C KN i (mb) B KN i (GeV 2 ) IP i (0.25 GeV 2 )t f i (0.93 GeV 2 )t ρ i (0.93 GeV 2 )t ω i (0.93 GeV 2 )t a i (0.93 GeV 2 )t III. KAON-INDUCED D (2317) AND D 1 (2460) PRODUCTION ON PROTON TARGET Firt, we calculate the total cro ection of the Λ + c D and Λ + c D 1 reaction The correonding unolarized differential cro ection read dσ d coθ = m m Λc 3cm 8π 1cm (1 M 1/2± 2 ), (8) 2 c, 2 where =( ) 2,θi the cattering angle of the outgoing meon relative to the beam direction, 1cm and 3cm are the and D ) three momenta in the center of ma frame, (D 1 1 = λ1/2 (, m 2 K, m 2 ) 2, (9) λ 1/2 (, m 2, m 2 D /1 Λ ) + c 3 = 2, (10) whereλ(x, y, z) i Källen function withλ(x, y, z)=(x y z) 2 4yz. The m, m, and m Λc are the mae of the meon, roton, andλ c, reectively. Here, we take m = MeV, m = MeV, and m Λc = MeV. Taking the ISI of the ytem into account, the full amlitude for the roce Λ + c D 0/1 i a um of the Born and ISI amlitude. With the Lagrangian given in the reviou ection, the Born amlitude of the ( 1 )( 2 ) Λ + c ( 4)D 0/1 ( 3) reaction can be obtained a, 1 B = g Λc D 0ū( 4, c )γ 5 u( 2, 2 ) ( 3 1 ) 2 m 2 D g KD 0 DFD 2 (( 3 1 ) 2, m D ), (11) M D 0 1 B = ig Λc D 1ū( 4, c )γ µ u( 2, 2 ) ( 3 1 ) 2 m 2 D ( g µν + ( 3 1 ) µ ( 3 1 ) ν ) m 2 D M D 1 g KD 1 D F2 D (( 3 1 ) 2, m D )ǫ ν ( 3), (12) where ū( 4, c ) and u( 2, 2 ) are the Dirac inor with c ( 4 ) and 2 ( 2 ) being the in (the four-momenta) of the outgoingλ c and the initial roton, reectively. Theǫ ν( 3 ) i the olarization vector of the D 1. Following the trategy of Ref. [46], the ISI amlitude can be written a M D /1 i IS I = d 2 k 16π 2 t T K K (, kt 2 ) M 1/2± D /1( 2 k t + 4 ), (13) where k t i the momentum tranfer in the reaction. With the formalim and ingredient given above, the total cro ection veru the beam momentum of the ytem for the Λ + c D 0 and K Λ + c D 1 reaction are evaluated. The reult for beam momentum P k from the reaction threhold to 20.0 GeV are hown in Fig. 3. Becaue the cutoff can not be well determined, the reult obtained with everal cutoff are alo reented. It i worthy mentioning that the value of the cro ection i very enitive to the model arameter α when varying the cutoff arameter α from 1.0 to 2.0. Thi i becaue the model arameter we elected are very cloe to the mae of the exchanged article. To make a reliable rediction for the cro ection of the Λ + c D 0 and Λ + c D 1 reaction thu require a good knowledge on the form factor. More and accurate exerimental data can alo be ued to contrain the value of the cutoff arameter. ( b) (a) D (2317) =1.0 =1.5 =2.0 (b) D 1 (2460) P K - (GeV) FIG. 3: (color online). Total cro ectionσfor the Λ + c D 0 (a) and K Λ + c D 1 (b) reaction a a function of the beam momentum. The reult in Fig. 3 alo how that the total cro ection increae harly near the D 0/1 Λ c threhold. At higher energie, the cro ection increae continuouly but relatively lowly comared with that near threhold. However, the total cro ection decreae but very lowly for the D 0 roduction in the Λ + c D 1 reaction when we change the beam energy P from 14.6 GeV to 20.0 GeV. With the increae of the cutoff, the total cro ection increae. Comaring the cro ection of the Λ + c D 0 reaction with that of the Λ + c D 1 reaction, we found that the line hae of the cro ection are very different. A oible exlanation for thi may be that KD interaction to form the D 0 i tronger than that KD interaction to form D 1 due to the D meon decay comletely to the final tate containing the D meon [1]. The reult how that the total cro ection for D 0 roduction i bigger than that for D 1 roduction. At a beam mo-

4 4 mentum about 14.6 GeV and the arameterα=1 the cro ection i of the order of 10 nb, which i quite large for an exerimentally obervation of the D 0 at current and future facilitie. Our reult ugget that it will take high energy, at leat above 14.6 GeV, to oberve the roduction of D 1 in the Λ + c D 1 reaction. ( b) 10 0 (a) D (2317) (b) D 1 (2460) P K - (GeV) FIG. 4: (color online). Total cro ection with or without ISI for the D 0/1 Λ c reaction a a function of the beam momentum for the D 0 cae (left) and the D 1 cae (right). To how the effect from the ISI, we comare the cro ection obtained with and without ISI for the cutoff ofα= 2.0 in Fig. 4, for he Λ + c D 0 (anel A) and K Λ + c D 1 (anel B) reaction, reectively. In Fig. 4, the dahed red line are the ure Born amlitude contribution, while the olid black line are the full reult. It how that the role of the ISI i to reduce the cro ection by aroximately 20%, in agreement with the concluion drawn from Ref. [39, 47, 48] that the ISI for or reaction reduce the cro ection. In addition to the total cro ection, we alo comute the differential cro ection for the Λ + c D 0 and Λ + c D 1 reaction a a function of the cattering angle of the outgoing meon relative to the beam direction at different beam energie,i.e,p =12.0,14.0,16.0, and 18.0 GeV. The theoretical reult are hown in Fig. 5. We note that the differential cro ection i the larget at the extreme forward angle and decreae with the increae of the cattering angle. Thi i becaue we have only conidered the contribution from the t-channel D and D exchange. It hould be ointed out that, if there are contribution from and u channel, there will be a clear bum (or eak) in the total cro ection which can be ditinguihed eaily. IV. SUMMARY In thi work, the roduction of the D 0 and D 1 reonance in the Λ + c D 0 and K Λ + c D 1 reaction wa tudied in an effective Lagrangian aroach. The roduction roce i decribed by the t-channel D 0 and D 0 meon exchange, reectively. The couling contant of the D 0 to KD and D 1 to KD are obtained from the chiral unitary the- d /dco( )( b) 10 0 (a) D (2317) (b) D (2460) co( ) FIG. 5: (color online). The Λ + c D 0 (a) and K Λ + c D 1 (b) differential cro ection at different energie with α = 2.0. The black olid line, red dah line, blue dah dot line, and traight dah dot dot line are obtained at beam energy P =12.0,14.0,16.0, and 18.0 GeV, reectively. ory [18, 19]. where the D 0 and D 1 reonance are dynamically generated. The initial tate interaction(isi) wa included by Pomeron and Reggeon exchange [46], which wa hown to reduce the cro ection by about 20%. the total and differential cro ection comuted can be ued to tet the molecular icture of the D 0 and D 1 meon in facilitie uch a OKA@U-70,SPS@CERN,and future J-PARC. Finally, we would like to tre that, thank to the abence of the channel, u channel, and background contribution in the Λ + c D and Λ + c D 1 reaction, future exerimental data for thee two reaction can be ued to imrove our knowledge of D 0 and D 1 roertie, which are at reent oorly known. Acknowledgment Y.H.thank Li-Sheng Geng for valuable dicuion. Thi work wa uorted by the Science and Technology Reearch Program of Chongqing Municial Education Commiion (Grant No. KJQN ), the Oened Fund of the State Key Laboratory on Integrated Otoelectronic (GrantNo. IOSKL2017KF19), and China Potdoctoral Science Foundation(No.2018M641143). Thi work alo wa uorted in art by the National Natural Science Foundation of China under Grant No and No

5 5 [1] M. Tanabahi et al. [Particle Data Grou], Review of Particle Phyic, Phy. Rev. D 98, (2018). [2] S. Godfrey and N. Igur, Meon in a Relativized Quark Model with Chromodynamic, Phy. Rev. D 32, 189 (1985). [3] B. Aubert et al. [BaBar Collaboration], Obervation of a narrow meon decaying to D + 0 π0 at a ma of 2.32-GeV/c 2, Phy. Rev. Lett. 90, (2003). [4] D. Beon et al. [CLEO Collaboration], Obervation of a narrow reonance of ma 2.46-GeV/c 2 decaying to D + π0 and confirmation of the D J (2317) tate, Phy. Rev. D 68, (2003) Erratum: [Phy. Rev. D 75, (2007)]. [5] Y. Mikami et al. [Belle Collaboration], Meaurement of the D J reonance roertie, Phy. Rev. Lett. 92, (2004). [6] P. Krokovny et al. [Belle Collaboration], Obervation of the D J (2317) and D J (2457) in B decay, Phy. Rev. Lett. 91, (2003). [7] B. Aubert et al. [BaBar Collaboration], Study of B D ( )+ J D ( ) decay, Phy. Rev. Lett. 93, (2004). [8] B. Aubert et al. [BaBar Collaboration], Obervation of a narrow meon decaying to D + 0 π0 γ at a ma of GeV/c 2, Phy. Rev. D 69, (2004). [9] B. Aubert et al. [BaBar Collaboration], A Study of the D J (2317) and D J(2460) Meon in Incluive c anti-c Production Near () (1/2) = 10.6-GeV, Phy. Rev. D 74, (2006). [10] Z. X. Xie, G. Q. Feng and X. H. Guo, Analyzing D 0 (2317)+ in the DK molecule icture in the Beth-Saleter aroach, Phy. Rev. D 81, (2010). [11] Y. J. Zhang, H. C. Chiang, P. N. Shen and B. S. Zou, Poible S -wave bound-tate of two eudocalar meon, Phy. Rev. D 74, (2006). [12] P. Bicudo, The Family of trange multiquark a kaonic molecule bound by hard core attraction, Nucl. Phy. A 748, 537 (2005). [13] A. Faeler, T. Gutche, V. E. Lyubovitkij and Y. L. Ma, Strong and radiative decay of the D 0 (2317) meon in the DKmolecule icture, Phy. Rev. D 76, (2007). [14] A. Faeler, T. Gutche, V. E. Lyubovitkij and Y. L. Ma, D K molecular tructure of the D 1 meon, Phy. Rev. D 76, (2007). [15] M. Cleven, H. W. Griehammer, F. K. Guo, C. Hanhart and U. G. Meißner, Strong and radiative decay of the D 0 (2317) and D 1, Eur. Phy. J. A 50, 149 (2014). [16] c. J. Xiao, D. Y. Chen and Y. L. Ma, Radiative and ionic tranition from the D 1 to the D 0 (2317), Phy. Rev. D 93, (2016). [17] A. Datta and P. J. O donnell, Undertanding the nature of D 0 and D 0 (2460) through nonletonic B decay, Phy. Lett. B 572, 164 (2003). [18] D. Gamermann, E. Oet, D. Strottman and M. J. Vicente Vaca, Dynamically generated oen and hidden charm meon ytem, Phy. Rev. D 76, (2007). [19] D. Gamermann and E. Oet, Axial reonance in the oen and hidden charm ector, Eur. Phy. J. A 33, 119 (2007). [20] M. Altenbuchinger, L.-S. Geng and W. Weie, Scattering length of Nambu-Goldtone boon off D meon and dynamically generated heavy-light meon, Phy. Rev. D 89, (2014). [21] A. Martinez Torre, K. P. Khemchandani and L. S. Geng, Bound tate formation in the DDK ytem, arxiv: [heh]. [22] F. K. Guo, P. N. Shen, H. C. Chiang, R. G. Ping and B. S. Zou, Dynamically generated 0 + heavy meon in a heavy chiral unitary aroach, Phy. Lett. B 641, 278 (2006). [23] F. K. Guo, P. N. Shen and H. C. Chiang, Dynamically generated 1 + heavy meon, Phy. Lett. B 647, 133 (2007). [24] J.B. Liu and M. Z. Yang, Sectrum of the charmed and b- flavored meon in the relativitic otential model, J. High. Energy. Phy 07, (2014) 106. [25] S. Fajfer and A. Praotnik Brdnik, Chiral loo in the ioin violating decay of D 1 and D 0 (2317), Phy. Rev. D 92, (2015). [26] X. Liu, Y. M. Yu, S. M. Zhao and X. Q. Li, Study on decay of D J (2317) and D J(2460) in term of the CQM model, Eur. Phy. J. C 47, 445 (2006). [27] J. Lu, X. L. Chen, W. Z. Deng and S. L. Zhu, Pionic decay of D j (2317), D j (2460) and B j (5718), B j (5765), Phy. Rev. D 73, (2006). [28] Z. G. Wang, Radiative decay of the D 0, D 1 and the related trong couling contant, Phy. Rev. D 75, (2007). [29] Y. B. Dai, C. S. Huang, C. Liu and S. L. Zhu, Undertanding the D + J (2317) and D+ J (2460) with um rule in HQET, Phy. Rev. D 68, (2003). [30] P. Colangelo and F. De Fazio, Undertanding D J (2317), Phy. Lett. B 570, 180 (2003). [31] P. Colangelo, F. De Fazio and A. Ozineci, Radiative tranition of D J (2317) and D J(2460), Phy. Rev. D 72, (2005). [32] P. Colangelo, F. De Fazio, F. Giannuzzi and S. Nicotri, New meon ectrocoy with oen charm and beauty, Phy. Rev. D 86, (2012). [33] Z. H. Guo, U. G. Meißner and D. L. Yao, New inight into the D 0 (2317) and other charm calar meon, Phy. Rev. D 92, no. 9, (2015). [34] M. Nielen, D + J (2317) D+ π0 decay width, Phy. Lett. B 634, 35 (2006). [35] Z. G. Wang and S. L. Wan, D 0 a a tetraquark tate with QCD um rule in heavy quark limit, Nucl. Phy. A 778, 22 (2006). [36] V. Obraztov [OKA Collaboration], High tatitic meaurement of the K + π 0 e + ν(ke3) decay formfactor, Nucl. Part. Phy. Proc , 1330 (2016). [37] B. Velghe [NA62-RK and NA48/2 Collaboration], K ± π ± γγ Studie at NA48/2 and NA62-RK Exeriment at CERN, Nucl. Part. Phy. Proc , 2720 (2016). [38] T. Nagae, The J-PARC roject, Nucl. Phy. A 805, 486 (2008). [39] Y. Dong, A. Faeler, T. Gutche and V. E. Lyubovitkij, Role of the hadron moleculeλ c (2940) in the D 0 Λ c (2286) annihilation reaction, Phy. Rev. D 90, (2014). [40] Y. Dong, A. Faeler, T. Gutche, S. Kumano and V. E. Lyubovitkij, Radiative decay ofλ c (2940) + in a hadronic molecule icture, Phy. Rev. D 82, (2010). [41] J. He, Z. Ouyang, X. Liu and X. Q. Li, Production of charmed baryonλ c (2940) + at PANDA, Phy. Rev. D 84, (2011). [42] J. J. Xie, Y. B. Dong and X. Cao, Role of theλ + c (2940) in the π D D 0 reaction cloe to threhold, Phy. Rev. D 92, (2015). [43] Y. Huang, J. He, H. F. Zhang and X. R. Chen, Dicovery otential of hidden charm baryon reonance via hotoroduction, J. Phy. G 41, (2014). [44] C. J. Xiao and D. Y. Chen, Analyi of the hidden bottom decay of Z b (10610) and Z b (10650) via final tate interaction, Phy. Rev. D 96, (2017). [45] H. Xu, J. J. Xie and X. Liu, Imlication of the oberved e + e π 0 for tudying the ψ(3770)π 0 roce, Eur. Phy. J. C

6 6 76, 192 (2016). [46] P. Lebiedowicz and A. Szczurek, K + reaction at high energie, Phy. Rev. D 85, (2012). [47] C. Hanhart, Meon roduction in nucleon-nucleon colliion cloe to the threhold, Phy. Ret. 397, 155 (2004). [48] V. Baru, A. M. Gaaryan, J. Haidenbauer, C. Hanhart, A. E. Kudryavtev and J. Seth, Production of eta meon in nucleon-nucleon colliion, Phy. Rev. C 67, (2003).