Masses and strong decays of open charm hexaquark states

Inspired by the recent discovery of the doubly charmed tetraquark state T + cc by the LHCb Collaboration, we perform a systematic study of masses and strong decays of open charm hexaquark states (cid:2) ( ∗ ) c (cid:2) ( ∗ ) c . Taking into account heavy quark spin symmetry breaking, we predict several bound states of isospin I = 0, I = 1, and I = 2 in the one boson exchange model. Moreover, we adopt the effective Lagrangian approach to estimate the decay widths of (cid:2) ( ∗ ) c (cid:2) ( ∗ ) c → (cid:3) c (cid:3) c and their relevant ratios via the triangle diagram mechanism, which range from a few MeV to a few tens of MeV. We strongly recommend future experimental searches for the (cid:2) ( ∗ ) c (cid:2) ( ∗ ) c hexaquark states in the (cid:3) c (cid:3) c invariant mass distributions.


Introduction
The quenched (or conventional) quark model can well describe the properties of traditional hadrons, especially the ground-state ones [1,2]. However, a large number of states that can not be easily explained by the quenched quark model have been accumulating since D * s0 (2317) and X (3872) were discovered by the BaBar [3] and Belle [4] collaborations in 2003. To unveil the nature of these exotic states, many pictures, such as multiquark states, hadronic molecules, and kinetic effects, have been proposed to understand them from different perspectives, including production mechanisms, mass spectra, and decay widths [5][6][7][8][9][10][11][12]. Among them, the hadronic molecular picture is probably the most popular one since most of these exotic states are located close to the mass thresholds of conventional hadron pairs. Two-body hadronic molecules are expected to exist in three configuraa e-mail: zhengmz11@buaa.edu.cn (corresponding author) b e-mail: lisheng.geng@buaa.edu.cn (corresponding author) tions, i.e., meson-meson, meson-baryon, and baryon-baryon. With more and more XY Z states, P c states, and T cc discovered in the past two decades [11], heavy hadronic molecules have attracted a lot attention.
Heavy quark spin symmetry (HQSS) plays an important role in describing heavy hadronic molecules [13][14][15][16]. HQSS dictates that the interaction between a light (anti)quark and a heavy (anti)quark is independent of the heavy (anti)quark spin in the limit of heavy quark masses. As the heavy mesons (D, D * ) and baryons ( c , * c ) belong to the same HQSS doublets, using them as building blocks, one expects three hidden charm multiplets of hadronic molecules, i.e.,D ( * ) D ( * ) , D ( * ) ( * ) c , and¯ c . With the assumption that X (3872) is aD * D bound state, in Ref. [14], Guo et al. adopted an effective field theory approach to study theD ( * ) D ( * ) system and predicted aD * D * bound state as the HQSS partner of X (3872). In Ref. [16], using the same approach we obtained a complete HQSS multiplet of hadronic molecules, D ( * ) ( * ) c , where P c (4440) and P c (4457) are regarded as theD * c bound states. Although the likely existence of ( * ) c ( * ) c bound states have been explored theoretically [17][18][19][20], they have not been observed experimentally.
Very recently, the LHCb Collaboration observed one open charm tetraquark state T + cc below the D 0 D * + mass threshold [21], which has already been predicted by a lot of theoretical works in the diquark and anti-diquark picture [22][23][24][25][26][27][28][29][30][31][32][33][34][35] as well as in the molecular picture [36][37][38][39][40][41][42][43][44][45]. In our previous work [45], using the one boson exchange (OBE) model we have systematically studied the D ( * ) D ( * ) system and predicted one isocalar D D * molecule with a binding energy of about 3 MeV, which is an open charm tetraquark state and may correspond to the T + cc state. Therefore, it is interesting to investigate the likely existence of open charm pentaquark and hexaquark states. In Ref. [46], we obtained a complete HQSS multiplet of hadronic molecules D ( * ) ( * ) c , consistent with Refs. [47][48][49]. One should note that the minimum quark content of the D ( * ) ( * ) c system is ccq, which may couple to the excited states of cc . Consequently, the D ( * ) ( * ) c molecules can mix with the excited doubly charmed baryons. In the present work, we investigate the mass spectrum and strong decays of c , which cannot couple to conventional baryons of ccq.
Experimental searches for dibaryons containing charm quarks is challenging because of the high energy required and the low production of charmed baryons. On the other hand, many theoretical studies have been performed. In Ref. [73], Liu et al. predicted one bound state generated by c N and c N couple channels in the OBE model, which was later explored in the chiral constituent quark model [74,75], lattice QCD simulations [76], and effective field theories [77]. Later, using the same approach, Li et al. systematically studied possible bound states made of two singly charmed baryons [17,78,79]. Some of these states were investigated in the quark delocalization color screening model [80,81]. After the discovery of the doubly charmed baryon cc by the LHCb Collaboration [82], Meng et al. investigated bound states made of two doubly charmed baryons [83]. Inspired by the discovery of the pentaquark states by the LHCb Collaboration [84], we studied the ( * ) cc ( * ) c system utilizing heavy antiquark diquark symmetry in the OBE model and effective field theory [85,86] and obtained a complete HQSS multiplet of hadronic molecules. Recently, lattice QCD simulations obtained a series of dibaryon states with open charm, such as cc c [87] and ccc ccc [88].
In this work, we investigate the mass spectrum and decay patterns of the c system dictated by HQSS. In Ref. [19], we parameterized the contact-range potential of the c system with three parameters. However, we can not fully determine the mass spectrum of the c system because no experimental data is available. In this work, we turn to the OBE model for help and calculate the mass spectrum. Then based on the so-obtained mass spectrum, we estimate the strong decay widths of c → c c through the effective Lagrangian approach, which has been widely used to explore the strong decays of hadronic molecules [89][90][91][92].
The manuscript is structured as follows. In Sect. 2 we present the details of the OBE model as applied to the

OBE potential
In this work, we employ the OBE model to derive the c interaction, which has been successfully used to investigate hadronic molecules composed of heavy hadrons [93][94][95][96]. The OBE interaction of two heavy hadrons is generated by the exchange of light mesons, π , σ , ρ, and ω. The vector mesons, ρ and ω, provide the short-range interaction, the scalar meson σ provides the medium-range interaction, and the π meson provides the long-range interaction. Given the exploratory nature of the present work, the contributions of other mesons are neglected [97]. To derive the OBE potential, the effective Lagrangians describing the interaction between a charmed baryon and a light meson are necessary. The relevant Lagrangians are the same as those of theD ( * ) ( * ) c interaction [97], which read where S c = ( 1 √ 3 c σ + * c ) denotes the superfield of c and * c dictated by HQSS. The coupling of the π meson to c and * c , g, was determined to be 0.84 in lattice QCD [98], which is smaller than the prediction of the quark model [73]. In this work we take g = 0.84. For the coupling to the σ meson, we estimate it using the quark model. From the nucleon and σ meson coupling of the liner sigma model, g σ N N = 10.2, the corresponding coupling is determined to be g σ = 2 3 g σ N N = 6.8 [99]. The couplings of the light vector mesons contain both electric-type(g v ) and magnetictype( f v ) ones, which are related via f v = κ v g v . The mass scale M represents the nucleon mass to make the coupling f v dimensionless. For the ρ and ω couplings, from SU(3)symmetry and the OZI rule, we obtain g ω = g ρ and f ω = f ρ . The electric-type coupling of the singly charmed baryons to the vector mesons is determined to be g v = 5.8 [73] as well Table 1 Couplings of charmed baryons to light mesons. The magnetic coupling of the ρ and ω mesons are defined as f v = κ v g v , and M refers to the mass scale (in units of MeV) involved in the magnetic-type couplings as the corresponding magnetic-type coupling is estimated to be κ v = 1.7 [100]. All the couplings are given in Table 1 for easy reference.
With the above Lagrangians, the OBE potentials for the c system in coordinate space can be obtained, which read where T 1 · T 2 is the isospin factor for the c system, and a 1 · a 2 and S 12 (r ) denote the spin-spin and tensor terms, respectively, and a denotes the spin operator of the Since the charmed baryons and light mesons involved in our study are not point-like particles, we introduce a form factor in the interaction vertices. Here we use a monopolar form factor (for more details we refer to Refs. [45,97]) where m E and q are the mass and 4-momentum of the exchanged meson. is an unknown cutoff, which is often fixed by reproducing some hadronic molecular candidates that can be related to hadronic molecules of interest via symmetries. Assuming X (3872) and P c (4312) as bound states of D D * andD c , we can determine the corresponding cutoff of the OBE model to be 1.04 GeV [45] and 1.12 GeV [97], respectively. In addition, the cutoff needed to reproduce the binding energy of the deuteron is 0.86 GeV [101]. All these results indicate that a reasonable value for the cutoff of the OBE model is about 1 GeV. In this work, assuming that the c system is similar to the nucleon-nucleon system, we take a cutoff of =0.86 GeV to calculate the binding energies of the c dibaryon system, With the above form factor, the functions δ, W Y , and W T in Eqs. (5)(6)(7)(8) need to be changed accordingly with λ = /m. The corresponding functions d, W Y , and W T read The generic wave function of a baryon-baryon system reads where |I M I denotes the isospin wave function and J M ( r ) are the spin and spatial wave functions. The dynamics in isospin space is embodied in the isospin factor T 1 · T 2 . In this work the total isospin is either 0, 1, or 2 for the ( * ) where L M L SM S |J M are the Clebsch-Gordan coefficients, |SM S the spin wavefunction, and Y L M L (r ) the spherical harmonics.
In the partial wave decomposition of the OBE potential, we encounter both spin-spin and tensor components In the present study, we consider both S and D waves. The relevant matrix elements are listed in Table 4.

Numerical results and discussions
For a pair of identical fermions, the total spin and isospin of  Table 2. Additionally, when deriving the OBE potential, we have assumed HQSS. As is well known, the charm quark mass is not large enough to strictly satisfy HQSS. In the present work, we take into account a HQSS breaking of the order of = 15% in the following way [85]: With the above OBE potentials, the binding energies of the

Effective Lagrangian approach
In this section, we further investigate the strong decays of ( * ) c ( * ) c molecules by the effective Lagrangian approach, which has been widely applied to study the strong decays of two-body and three-body bound states [89][90][91][92][104][105][106]. The transition from c → c c is mediated by the exchange of a π meson, 1 depicted by the triangle diagrams shown in Fig. 1. The isospin and spin of the c c pair are 0 and 0/1, respectively, which can help us select the quantum numbers of initial states because for strong decays isospin and angular momentum are conserved. As a result, in this work, we only study the strong decay of three isoscalar states, J P = 0 + c c , J P = 1 + * c c , and J P = 0 + * c * c , denoted as H cc , H * cc and H * * cc in the following, respectively. 1 In principle, the exchange of ρ is also allowed. However, c can only transit into c ρ via the magnetic term, which is regarded as subleading for interactions involving vector mesons. As a result, we can safely neglect the ρ meson exchange in comparison with the π meson exchange.
To describe the interaction vertices of Fig. 1, we need the following Lagrangians for the c bound states and their constituents where p E is the Euclidean momentum and is the size parameter.
To estimate the couplings between bound states and their constituents, we employ the compositeness condition [107][108][109]. For dibaryon bound states with total angular momentum J = 0 the compositeness condition reads where (m 2 c c ) is the self energy of a c bound states as shown in Fig. 2.
Substituting the transverse term into Eq. (25), we can also determine the couplings of a J = 1 dibaryon state to their constituents.  The mass operators for the c dibaryons read as follows: where In Refs. [110,111], the above cutoff is found to be around 1 GeV. One should note, however, that the couplings of a molecular state to its components are related to its binding energy [112], and hence we take a cutoff of = 0.86 GeV, the same as that used to study the c dibaryons, to determine the couplings in this work. In Table 3, we present the H ( * * ) cc couplings to their components. Because the binding energy of H * * cc is larger than those of H * cc and H cc , g H * * cc * c * c is larger than g H * cc * c c and g H cc c c as well.
The other vertices of the triangle diagrams of Fig. 1 can be classified into two categories, c → c π and * c → c π . The Lagrangians describing these interactions are given by where f π = 132 MeV and the couplings g π c c and g π c * c can be determined by fitting to experimental data. From the decay widths of ( c → c π) = 1.89 MeV and ( * c → c π) = 15.0 MeV [113], we obtain the couplings g π c c = 0.55 and g π c * c = 0.97, consistent with other works [73,114]. In addition, we find that the two couplings approximately satisfy the relationship g π c * c = √ 3g π c c , given by the quark model [73].
With the above Lagrangians, the amplitudes of H ( * * ) cc → c c can be easily written down whereū c andū T c represent the spinors of the final-state c c pair, and F(q, m E , ) is the form factor which not only removes ultraviolet divergence of the loop diagram, but also takes into account the off-shell effects.
The m E is the mass of the exchanged particle. The cutoff is expressed by = m E + α QC D , where QC D is around 200-300 MeV and the dimensionless parameter α is around 1 [115]. Thus we vary the cutoff from 0.4 to 0.6 GeV to estimate the uncertainties induced . With the amplitudes of H ( * * ) cc → c c determined, one can obtain the corresponding partial decay widths as is the mass of the corresponding molecule. At last, we have to multiply a factor 1/2 to the above result to take into account the statistics of identical particles.

Numerical results and discussions
In Fig. 3 Fig. 4, where the ranges of masses are obtained from the upper and lower limits of binding energies of Table 2, and the cutoff is taken to be 0.5 GeV. One can see that the decay widths of H cc and H * cc only vary by several MeV, while the decay width of H * * cc varies by tens of MeV, which shows that the decay width of H ( * * ) cc molecules are not very sensitive to their masses.
All of the H ( * * ) cc molecules can decay into c c , which indicates that all of them could be detected in the c c mass distributions. In Fig. 5 we show the ratios of the decay widths of the H * * cc and H * cc molecules to that of the H cc molecule, and

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: All data generated or analysed during this study are included in this published article.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .