Production of the spin partner of the X(387) in e + e collisions arxiv:40.4674v [hep-ph] 5 Nov 04 Feng-Kun Guo a,, Ulf-G. Meißner a,b,, Zhi Yang a, a Helmholtz-Institut für Strahlen- und Kernphysik and Bethe Center for Theoretical Physics, Universität Bonn, D-535 Bonn, Germany b Institut für Kernphysik, Institute for Advanced Simulation, and Jülich Center for Hadron Physics, Forschungszentrum Jülich, D-545 Jülich, Germany November 6, 04 Abstract We study the production of the spin partner of the X(387), which is a D D bound state with quantum numbers J P C = ++ and named X (40) here, with the associated emission of a photon in electron positron collisions. The results show that the ideal energy region to observe the X (40) in e + e annihilations is from 4.4 GeV to 4.5 GeV, due to the presence of the S- wave D D (40) and D D (460) thresholds, respectively. We also point out that it will be difficult to observe the γx(40) at the e + e center-of-mass energy around 4.6 GeV. E-mail address: fkguo@hiskp.uni-bonn.de E-mail address: meissner@hiskp.uni-bonn.de E-mail address: zhiyang@hiskp.uni-bonn.de
In the heavy quarkonium mass region, the so called XY Z states have been observed, and many of these quarkonium-like states defy a conventional quark model interpretation. They are therefore suggested to be exotic. The X(387), discovered by the Belle Collaboration [], is the one of the most interesting exotic states. As the mass of the X(387) is extremely close to the D 0 D 0 threshold, it is regarded as one especially promising candidate for a hadronic molecule. Effective field theory (EFT) can cope with the interaction between heavy mesons in bound state systems at low energies. For such a kind of systems, heavy quark symmetry is relevant due to the presence of the heavy quark/antiquark in the meson/antimeson. This fact leads to predictions of new states as partners of the observed XY Z states in the hadron spectrum. For example, with an EFT description of the heavy mesonic molecules, the heavy quark symmetry can be used to predict the existence of the spin and bottom partners of the X(387) [, 3]. The spin partner of the X(387), called X (40) hereafter, is predicted to exist as the S-wave bound state of D D with quantum numbers ++ []. Such a state was also expected to exist in other models, see Refs. [4 8]. It is different from the X(387) in several aspects: first, being an isoscalar state it should decay into the J/ψπππ with a branching fraction much larger than that for the J/ψππ because the J/ψρ and J/ψω thresholds are far below the mass of the X (40) (very different to the case of the X(387)); second, it is expected to decay dominantly into open charm mesons, D D, D D and D D, in a D-wave with a width of the order of a few MeV [9]; third, its mass as set by the D D threshold is higher than the quark model prediction for the first radially excited χ c [0]. The significance of the X (40) state is that its mass should be approximately given by the M X (40) M X(387) + M D M D 40 MeV () as dictated by heavy quark spin symmetry for heavy-flavor hadronic molecules [3, ]. Notice that a state with the same quantum numbers ++ was also predicted in the tetraquark model []. However, the fine splitting between the ++ and ++ tetraquarks, which was predicted to be 70 MeV in Ref. [], is not locked to that between the D and D. Similarly, the splitting between the P c c states in the Godfrey Isgur quark model is 30 MeV [0], also much smaller than M D M D. Therefore, if a ++ state will be observed in experiments with a mass around 40 MeV, the mass by itself would already be a strong support for the hadronic molecular nature of both the X(387) and the tensor state. As a result, searching for a ++ state with a mass around 40 MeV is very important even for understanding the nature of the X(387). However, although the X(387) has been observed by many other experiments after its discovery [3 8], no evidence for the existence of its spin partner has been reported. In Ref. [9], it is shown that the prompt production of the X (40) presents a significant discovery potential at hadron colliders. In this paper, we will investigate the production of the X (40) associated with the photon radiation in electron positron collisions. This work presents an extension of the study on the production of the X(387) as a D D molecule in charmonia radiative transitions reported in Ref. [0]. In that paper, it was shown that the favorite energy regions for the X(387)γ production are around the Y (460) mass and 4.45 GeV. Later on, the BESIII Collaboration observed events for the process Y (460) X(387)γ [], which may be regarded as a support of the dominantly molecular nature of the X(387). Since the existence of the D D bound state, the X (40), is the consequence of the heavy quark spin symmetry of the molecular nature of the X(387), the production of the X (40) in e + e collisions in the energy range of the BESIII experiment [] thus provides an opportunity to search for new charmonium-like states on the one hand and can offer useful information towards understanding the X(387) on the other hand. The production of the X(387) through the radiative decay of the ψ(460) charmonium is considered in Ref. [3] using heavy hadron chiral perturbation theory along with the X-EFT [4]. Then,
γ γ D0 D 0 D0 D 0 D 0 X (40) D 0 X (40) (a) (b) Figure : Relevant triangle diagrams for the production of the X (40) in vector charmonia radiative decays. The charge-conjugated diagrams are not shown here. Ref. [0] studied the X(387) production by considering the contribution from intermediate charmed meson loops, and it was argued that the dominant mechanism is as follows: the initial charmonium is coupled to a pair of charmed mesons with one being S-wave with s P l =, where sl is the total angular momentum of the light-flavor cloud in the charmed meson, and the other being P -wave with s P l = 3 +, and the P -wave charmed meson radiatively transits to a D(D ) which coalesce with the other S-wave charmed meson, D ( D), into the X(387). The spin partner of X(387), the X (40), can be produced by a similar mechanism as shown in Fig.. Notice that the X (40) couples to D D instead of D D + c.c., as it is in the case of the X(387). We will only consider the neutral charmed mesons in the loops because the photonic coupling between the P -wave and S- wave charmed mesons for the neutral ones is much larger than that for the charged. This is due to cancellation of contributions from the charm and down quarks in the charged mesons, see, e.g. [5]. In the loops, the X (40) couples to the D 0 D 0 pair in an S-wave. With the quantum numbers being, the initial charmonium can couple to one P -wave and one S-wave charmed meson in either S- or D-wave. Since both the initial charmonium and the X (40) in the final state are close to the corresponding thresholds of the charmed-meson pairs, we are able to use a power counting in velocity of the intermediate mesons. Following the power counting rules as detailed in Ref. [6] and presented in the case of interest in Refs. [0,7], the dominant contribution comes from the case when the coupling of the initial charmonium to the charmed mesons is in an S-wave. In this case, the initial charmonium should be a D-wave state in the heavy quark limit m c as a consequence of heavy quark spin symmetry [8]. The charmed mesons can be classified according to the total angular momentum of the light degrees of freedom s l and collected in doublets with total spin J = s l ± in the heavy quark limit. The s P l = states correspond to charmed mesons in the doublet (0, ), here denoted as (P, V ), whereas the s P l = 3 + states correspond to charmed mesons in the doublet ( +, + ), denoted as (P, P ). To describe these heavy mesons, we choose the two-component notation introduced in Ref. [9]. The notation uses matrix fields, and is convenient for nonrelativistic calculations. The fields for the relevant heavy meson states are H a = V a σ + P a, T i a = P ij a σj + 3 P i a + i 6 ɛ ijkp j a σk, () for the s P l = (S-wave) and s P l = 3 + (P -wave) heavy mesons, respectively, where σ are the Pauli matrices, and a is the flavor index for the light quarks. In Eq. (), P a and V a annihilate the pseudoscalar 3
and vector heavy mesons, respectively, and P a and P a annihilate the excited axial-vector and tensor heavy mesons, respectively. Under the same phase convention for charge conjugation specified in Ref. [0], the fields annihilating the mesons containing an anticharm quark are [30] H a = Va σ + P a, T a i ij = P a σj + 3 P i a i 6 ɛ ijk P j a σk. (3) In nonrelativistic limit, the field for the D-wave charmonium state can be written as [3] J ij = 3 ( ψ i σ j + ψ j σ i) δ ij ψ σ, (4) 5 5 where ψ annihilates the D-wave vector charmonium, and the spin-0 and spin- states irrelevant for our study are not shown. In order to calculate the triangle diagrams in Fig., we need the Lagrangian for coupling the D-wave charmonia to the - 3 + charmed-meson pair as well as that for the E radiative transitions between the charmed mesons [0] L = g [( 4 Tr T j a σ i H a H ) a σ i Ta j J ij] + c [ ] a Tr TaH i a E i + H.c., (5) a where in the first term the Einstein summation convention is used while for the latter we distinguish the coupling constants for different light flavors because there is no isospin symmetry in the electromagnetic interaction. Moreover, we parametrize the coupling of the X (40) to the pair of vector charm and anticharm mesons as L X = x X ij ( D 0 i D 0 j + D + i D j) + H.c.. (6) With the above preparations, we can now proceed to calculate quantitatively the production of the γx (40) in electron positron collisions. Although in the heavy quark limit the production of the D-wave vector heavy quarkonium or the pair of and 3 + heavy mesons are suppressed due to spin symmetry [8], we can expect a large spin symmetry breaking in the charmonium mass region above 4 GeV. This may be seen from similar values of electronic widths of the excited vector charmonia. Thus, we will assume that the production of the γx (40) occurs through the D-wave charmonia or the D-wave components of excited vector charmonia. Without any detailed information about the values of the coupling constants, we can predict the energy regions with the maximal production cross sections. In Fig., we show the dependence of the decay width of a D-wave charmonium into the γx (40), divided by (g 4 x ), on the mass of the D-wave charmonium or the center-ofmass energy of the e + e collisions. The value of the photonic coupling c u does not affect the shape of the dependence either. Nevertheless, we took c u = 0.4 which is a typical value evaluated from various quark model predictions for the decay widths Γ(D 0 γd( )0 ) [3 33]. In the figure, the dashed curve is obtained neglecting the widths of the D and D states, and the solid curve is the result of evaluating the triangle loop integrals with constant widths for the D and D as done in Ref. [0]. The maximum around 4.447 GeV and the local minimum around 4.49 GeV of the dashed curve are due to the presence of Landau singularities [34] of triangle diagrams in the complex plane at (4.447±i 0.003) GeV (for the D loop) and (4.49±i 0.003) GeV (for the D loop), respectively (for a discussion of the Landau singularities in the triangle diagrams of heavy quarkonium transitions, we refer to Ref. [7]). The two cusps on both sides of the shoulders of the peak show up at the thresholds of the D D and D D. 4
Γ(ψ(nD) γ X)/(g4x GeV) kev 40 30 0 0 0 4.5 4.30 4.35 4.40 4.45 4.50 4.55 4.60 M ψ [GeV] Figure : Dependence of the partial decay width of a D-wave charmonium into γx (40) on the mass of the charmonium. The solid and dotted curves are obtained with and without taking into account the widths of the D (40) and D (460), respectively. Here, c u = 0.4 is used. From the figure, it is clear that the ideal energy regions for producing the γx (40) in e + e collisions are around the D D and D D thresholds, i.e. between 4.4 GeV and 4.5 GeV. It is also clear that the mass region of the Y (460) is not good for the production of the γx (40), contrary to the case of the γx(387). In order to quantify the relative production rate of the γx (40) with respect to the γx(387), we require the Y (460) to couple to the - 3 + meson pair as follows L Y = y Y i ( D i D ) y ( a a D a Di a + i ɛ ijk Y i D j D ) a a k Da k D j a + y ( Y i D ij D ) a a j Da j ij D a + H.c., (7) where we have assumed isospin symmetry in the couplings and the flavor index a runs over up and down quarks. Notice that if the Y (460) is a pure D D (here and in the following the charge conjugated channels are dropped for simplicity) molecule [35, 36], it would not couple to the D D and D D as given by the y and y terms, and thus cannot decay into the γx (40). These two terms are included to allow the decay to occur. Because the X (40) is the spin partner of the X(387), for a rough estimate, we can assume that x takes the same value as the coupling constant of the X(387) to the D D. We also assume that the values of y and y are related to y by a spin symmetry relation for D-wave charmoinia. Comparing Eq. (5) with Eq. (7) one obtains y = y/ and y = 6y/0. Then, the ratio of the partial decay widths of the Y (460) to the γx (40) and the γx(387) can be estimated parameter-free, and is Γ(Y (460) γx (40)) Γ(Y (460) γx(387)) 0. (8) In the above ratio, whether or not to take into account the finite widths of the P -wave charmed mesons only results in a minor change of %. It is clear that unless the Y (460) couples to the D D and/or The possibility of the Y (460) to have D D and D D components was discussed in Ref. [37]. 5
D D with a coupling much larger than that for the D D, which is less possible, the branching fraction of the Y (460) γx (40) is much smaller than that of the Y (460) γx(387). Given that the number of events for the latter process as observed at BESIII is the order of 0 [], it is unlikely to make an observation of the γx (40) at an energy 4.6 GeV at BESIII. To summarize, it is generally expected that the X(387) as a hadronic molecule has a spin partner close to the D D threshold. In this paper, we have investigated the production of the γx (40) in e + e collisions. According to our calculation, we strongly suggest to search for the X (40) associated with a photon in the energy region between 4.4 GeV and 4.5 GeV in e + e collisions. Besides, the width ratio of the Y (460) decaying to γx (40) and γx(387) is quite small, at the order of 0. Thus observing the γx (40) at an energy around 4.6 GeV would be unlikely in the BESIII experiment according to the current result of Y (460) γx(387). Acknowledgments This work is supported in part by the DFG and the NSFC through funds provided to the Sino-German CRC 0 Symmetries and the Emergence of Structure in QCD (NSFC Grant No. 6303), by the EU Integrated Infrastructure Initiative HadronPhysics3 (Grant No. 8386), and by NSFC (Grant No. 65005). References [] S. K. Choi et al. [Belle Collaboration], Phys. Rev. Lett. 9 (003) 600 [hep-ex/030903]. [] J. Nieves and M. P. Valderrama, Phys. Rev. D 86 (0) 056004 [arxiv:04.790 [hep-ph]]. [3] F.-K. Guo, C. Hidalgo-Duque, J. Nieves and M. P. Valderrama, Phys. Rev. D 88 (03) 054007 [arxiv:303.6608 [hep-ph]]. [4] N. A. Törnqvist, Z. Phys. C 6 (994) 55 [hep-ph/93047]. [5] C. Y. Wong, Phys. Rev. C 69 (004) 0550 [hep-ph/03088]. [6] E. S. Swanson, J. Phys. Conf. Ser. 9 (005) 79. [7] R. Molina and E. Oset, Phys. Rev. D 80 (009) 403 [arxiv:0907.3043 [hep-ph]]. [8] Z. F. Sun, Z. G. Luo, J. He, X. Liu and S. L. Zhu, Chin. Phys. C 36 (0) 94. [9] F.-K Guo et al., in preparation. [0] S. Godfrey and N. Isgur, Phys. Rev. D 3 (985) 89. [] F.-K. Guo, C. Hanhart and U.-G. Meißner, Phys. Rev. Lett. 0 (009) 4004 [arxiv:0904.3338 [hep-ph]]. [] L. Maiani, F. Piccinini, A. D. Polosa and V. Riquer, Phys. Rev. D 7 (005) 0408 [hepph/04098]. [3] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 7 (005) 0703 [hep-ex/04060]. 6
[4] D. Acosta et al. [CDF Collaboration], Phys. Rev. Lett. 93 (004) 0700 [hep-ex/030]. [5] V. M. Abazov et al. [D0 Collaboration], Phys. Rev. Lett. 93 (004) 600 [hep-ex/0405004]. [6] R. Aaij et al. [LHCb Collaboration], Eur. Phys. J. C 7 (0) 97 [arxiv:.530 [hep-ex]]. [7] R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 0 (03) 00 [arxiv:30.669 [hepex]]. [8] S. Chatrchyan et al. [CMS Collaboration], JHEP 304 (03) 54 [arxiv:30.3968 [hep-ex]]. [9] F.-K. Guo, U.-G. Meißner, W. Wang and Z. Yang, Eur. Phys. J. C 74 (04) 3063 [arxiv:40.636 [hep-ph]]. [0] F.-K. Guo, C. Hanhart, U.-G. Meißner, Q. Wang and Q. Zhao, Phys. Lett. B 75 (03) 7 [arxiv:306.3096 [hep-ph]]. [] M. Ablikim et al. [BESIII Collaboration], Phys. Rev. Lett. (04) 0900 [arxiv:30.40 [hep-ex]]. [] D. M. Asner, T. Barnes, J. M. Bian, I. I. Bigi, N. Brambilla, I. R. Boyko, V. Bytev and K. T. Chao et al., Int. J. Mod. Phys. A 4 (009) S [arxiv:0809.869 [hep-ex]]. [3] A. Margaryan and R. P. Springer, Phys. Rev. D 88 (03) 0407 [arxiv:304.80 [hep-ph]]. [4] S. Fleming, M. Kusunoki, T. Mehen and U. van Kolck, Phys. Rev. D 76 (007) 034006 [hepph/070368]. [5] Fayyazuddin and Riazuddin, Phys. Rev. D 48 (993) 4. [6] F.-K. Guo, C. Hanhart, G. Li, U.-G. Meißner and Q. Zhao, Phys. Rev. D 83 (0) 03403 [arxiv:008.363 [hep-ph]]. [7] F.-K. Guo, U.-G. Meißner and C. P. Shen, Phys. Lett. B 738 (04) 7 [arxiv:406.6543 [hepph]]. [8] X. Li and M. B. Voloshin, Phys. Rev. D 88 (03) 0340 [arxiv:307.07 [hep-ph]]. [9] J. Hu and T. Mehen, Phys. Rev. D 73 (006) 054003 [hep-ph/053]. [30] S. Fleming and T. Mehen, Phys. Rev. D 78 (008) 09409 [arxiv:0807.674 [hep-ph]]. [3] Fayyazuddin and O. H. Mobarek, Phys. Rev. D 50 (994) 39. [3] J. G. Körner, D. Pirjol and K. Schilcher, Phys. Rev. D 47 (993) 3955 [hep-ph/90]. [33] S. Godfrey, Phys. Rev. D 7 (005) 05409 [hep-ph/0508078]. [34] L. D. Landau, Nucl. Phys. 3 (959) 8. [35] G. J. Ding, Phys. Rev. D 79 (009) 0400 [arxiv:0809.488 [hep-ph]]. [36] Q. Wang, C. Hanhart and Q. Zhao, Phys. Rev. Lett. (03) 3003 [arxiv:303.6355 [hepph]]. [37] M. Cleven, Q. Wang, F.-K. Guo, C. Hanhart, U.-G. Meißner and Q. Zhao, Phys. Rev. D, in press [arxiv:30.90 [hep-ph]]. 7