Doubly charmed dibaryon states

In this work, we study the doubly charmed dibaryon states with the $qqqqcc$ ($q =u, d, s$) configuration. The mass spectra of doubly charmed dibaryon states are obtained systematically within the chromomagnetic interaction model. In addition to the mass spectrum analysis, we illustrate their two-body strong decay behaviours. Our results suggest that there may be narrow states or even stable states that cannot decay through the strong interaction. We hope that our results will provide valuable information for further experimental searches for doubly charmed dibaryon states.

In this work, we study the doubly charmed dibaryon states with the qqqqcc (q = u, d, s) configuration. The mass spectra of doubly charmed dibaryon states are obtained systematically within the chromomagnetic interaction model. In addition to the mass spectrum analysis, we illustrate their two-body strong decay behaviours. Our results suggest that there may be narrow states or even stable states that cannot decay through the strong interaction. We hope that our results will provide valuable information for further experimental searches for doubly charmed dibaryon states.

I. INTRODUCTION
With the accumulation of high-precision experimental data, the search for exotic hadronic matter is entering a new phase. In the last two decades, numerous of charmoniumlike XY Z states, P c states, and the double charm T + cc state have been reported experimentally, stimulating extensive studies on exotic hadron spectroscopy (see review articles [1][2][3][4][5] for details). Recently, the LHCb Collaboration announced new observations of the T a cs0 (2900) ++,0 state in the processes B 0 →D 0 D + s π − and B − →D − D + s π + [6, 7], and a P Λ ψs (4338) state in B − → J/ψΛp [8]. In addition, the CMS Collaboration confirmed the observation of the X(6900) state and found more structures in the di-J/ψ invariant mass spectrum [9], and the ATALAS Collaboration confirmed the discovery of the P c states and X(6900) [10]. Clearly, we are entering a new era of "Particle Zoo 2.0," which may deepen our understanding of the non-perturbative behaviour of the strong interaction.
Among the various exotic hadrons, the hexaquark or dibaryon state is a special group. The search for the dibaryon has a long history with many twists and turns. One can find that a review of the long-standing search for the dibaryon system, from the early days until the last few years [11]. The best known dibaryon state is the deuteron, which was discovered in 1932 by Urey, Brickwedde, and Murphy [12]. It is a loosely bound state of a proton and neutron with a quark content of uuuddd, I = 0, J = 1.
Another widely believed non-strange dibaryon candidate is d * (2380). Dyson and Xuong first predicted a wellknown candidate ∆∆ resonance state based on the SU(6) symmetry of the strong interaction [13], after proposing the quark model [14]. In 1989, Goldman et al. proposed For the strange dibaryon, Jaffe first used the MIT bag model to predict the existence of the H dibaryon, which is a bound state containing the uuddss quark component, in 1977 [33]. He showed that this particle should be bound with an energy of about 70 MeV relative to the ΛΛ threshold. Meanwhile, the model study based on the Goldstone boson exchange interaction cast doubt on the existence of the H dibaryon [34]. Similarly, experimental searches have ruled out the existence of such a deeply bound state [35,36]. In addition, the H dibaryon has been the subject of numerous lattice QCD calculations in recent years and the results have ruled out a deeply bound state, suggesting rather an unbound or a very weakly bound state if it exists [37][38][39]. In addition to these studies, this state had been extensively researched using various models [40][41][42][43][44][45][46][47]. Another interesting candidate was introduced later by Goldman et al. [48], which is expected to appear as a N Ω resonance. Its quark content (uudsss, I = 1/2, J = 2) is relatively favorable. The STAR experiment at the Relativistic Heavy-Ion Collider (RHIC) has made significant progress in the search for the N Ω state, supporting its existence [49]. Furthermore, the ALICE Collaboration at the LHC has studied the strong interactions between hadrons, providing further support for the possible formation of the N Ω state [50]. J. Leandri and B. Silvestre-Brac conducted a series of investigations on dibaryons in the framework of a pure chromomagnetic Hamiltonian [51][52][53]. They showed that the inclusion of heavy quarks in the dibaryon sector can lead to configurations that are stable against decay into two baryons. Furthermore, they proposed a number of new heavy states that could be stable under the strong interaction. Accordingly, a more realistic way to search arXiv:2209.10440v3 [hep-ph] 16 May 2023 for the dibaryons is to turn to the heavy flavor states, where the structures are easier to detect experimentally.
First, for fully heavy dibaryons, Junnarkar et al. used a state-of-the-art lattice QCD calculation and suggested that Ω c Ω cc (sscscc), Ω b Ω bb (ssbsbb), and Ω ccb Ω cbb (ccbcbb) were stable under strong and electromagnetic interactions [54]. They also found that the binding of these dibaryons became stronger as their mass increased. Lyu et al. used lattice QCD to study the Ω ccc Ω ccc in the 1 S 0 channel and found that this system is loosely bound at about 5.68 MeV [55]. Immediately afterward, the extended one-boson exchange (OBE) model [56] was used to study the existence of ΩΩ, Ω ccc Ω ccc , and Ω bbb Ω bbb systems. Additionally, the QCD sum rule approach suggested that fully heavy dibaryons should exist as weakly bound ones [57]. Huang et al. systematically studied the fully heavy dibaryons, and their results suggested the existence of Ω ccc Ω bbb and Ω ccc Ω ccc and the absence of Ω bbb Ω bbb [58]. However, many studies within quark models have disfavored the bound fully heavy dibaryons [59][60][61][62][63][64]. In particular, Richard et al. explored the possibility of heavy dibaryons with three charm quarks and three beauty quarks (bbbccc) in a potential model, concluding that there is no evidence for a stable state in such a very heavy flavor six-quark system [62]. In addition, Alcaraz-Pelegrina et al. studied the fully heavy dibaryons with the diffusion Monte Carlo method within the quark model and found that all these states are above the thresholds of the two fully heavy baryons [63]. Similarly, Weng et al. systematically studied the mass spectra of the fully heavy dibaryons in an extended chromomagnetic model and found no stable state below the corresponding baryon-baryon thresholds [64].
In addition to the fully heavy dibaryons, the stability of the dibaryon system consisting of six quarks or antiquarks has been studied in the context of a string model [65], where the existence of the six-quark (q 6 ) states and the three-quark-three-antiquark (q 3q3 ) states was discussed, and the investigations were extended to the (q 3 Q 3 ) and (Q 3q3 ) systems. Later, Huang et al. investigated possible N Ω-like dibaryons N Ω ccc and N Ω bbb with quantum numbers I(J P ) = 1/2(2 + ) within the framework of the quark delocalization color screening model [66]. Here, both of these states are bound, and the binding energy increases as the quarks of the system become heavier. In addition, inspired by the LHCb observation of three pentaquark states [67], theoretical studies predicted expect the light-quark structures of theD ( * ) and the Ξ ( * ) cc to be almost identical. Thus, theoretical work has also focused on the triply heavy dibaryons [68][69][70][71][72].
Given this situation, we are interested in further exploring the doubly heavy dibaryon state. The double heavy dibaryon states (qqqqQQ) have been discussed using various approaches, including the constituent quark model (CQM) [83][84][85], the quark delocalization color screening model [86][87][88], the lattice QCD calculation [89], the chiral constituent quark model [90], the nucleonnucleon interaction model [91,92], the chiral effective field theory (EFT) [94][95][96], the QCD sum rules [97,98], the relativistic six-quark equations [99], the complex scaling method [100], and the one-boson-exchange (OBE) model [101][102][103][104][105][106]. Meanwhile, Richard pointed out that there is no obvious gain with respect to the best threshold made of two baryons in a colour-additive model, but flavour independence dictates that before spin corrections, the threshold (QQq) + (qqq) is lower than 2 (Qqq) [107]. In this work, the mass spectrum of these doubly heavy dibaryon states will be obtained in the framework of the chromomagnetic model, which is usually used to calculate the mass spectrum of ordinary hadron states or multiquark states [3]. The present paper discusses the mass spectrum of these doubly charmed dibaryon states and briefly discusses the stability of various states. It may provide valuable information to further experimental searches for doubly charmed dibaryon states. With the running of the High-Luminosity LHC [108], the hunt for this new type of dibaryon may become a reality.
The remainder of this paper is organized as follows. In Sec. II, we introduce the chromomagentic interaction model. According to different permutation symmetries of identical particles, we construct the F f lavor ⊗ φ color ⊗ χ spin wave functions for the S-wave qqqqcc dibaryon state. The mass spectrum and decay mode of the qqqqcc hexaquark state with different constituent quarks are discussed in Sec. III. Finally, a short summary is given in Sec. IV.

II. HAMILTONIAN IN CMI MODEL
In the chromomagnetic interaction (CMI) model, the mass of the ground hadron state can be described by the effective Hamiltonian [64,[75][76][77][109][110][111][112] where V C ij = λ i · λ j and V CMI ij = λ i · λ j σ i · σ j represent the color and chromomagnetic interaction between the quarks, respectively, with σ i and λ i being the Pauli matrices and Gell-Mann matrices, respectively. The mass parameter of the quark pair m ij = 1/4(m i +m j )+4/3A ij contains the effective quark mass m i (m j ) and the color interaction strength A ij . The v ij is the effective coupling constant between the i-th quark and the j-th quark, which determines the mass gaps.
To estimate the mass spectrum of the doubly charmed dibaryon states, we extract the effective coupling parameters m ij and v ij from the conventional baryon masses, which are presented in Table I. To calculate the mass spectrum of the doubly charmed dibaryon states in the CMI model, we need the information on the total wave function, which is composed of the spatial, flavor, color, and spin wave functions, i.e., According to the Pauli principle, this wave function should be completely antisymmetric when identical quarks are exchanged. Here, we consider only the lowlying S-wave dibaryon states, so that their spatial wave functions are symmetric and become trivial. Thus, the F f lavor ⊗ φ color ⊗ χ spin wave functions of the dibaryon states should be completely antisymmetric when identical quarks are exchanged. For the doubly charmed dibaryon system, all the possible flavor combinations are nnnncc (I = 2, 1, 0), nnnscc (I = 3/2, 1/2), nnsscc (I = 1, 0), sssncc, and sssscc. According to the symmetry properties, we can divide the qqqqcc dibaryon system into three groups, which are presented in Table II. Here, we have to construct the corresponding F f lavor ⊗ φ color ⊗ χ spin wave functions with the symmetries {1234}{56}, {123}4{56}, and {12}{34}{56}. We use the notation {123...}, where quarks 1, 2, 3, and so on have antisymmetric properties. In previous works, there are abundant methods on the wave functions of six-quark systems [29,40,83,[113][114][115]. The light flavor SU(3) symmetry and symmetry breaking have been considered in these works. In addition, the same method we used was applied in Refs. [28,116,117].

A. Wave function for {1234}{56} symmetry
First, we construct the wave function of the nnnncc (I = 2, 1, 0) and sssscc dibaryon states. Because the nnnncc (I = 2) dibaryon states and sssscc dibaryon states have exactly the same symmetry requirements, we only need to concentrate on the nnnncc (I = 2, 1, 0) dibaryon states. Here, we introduce the Young tableau to TABLE II. All possible flavor combinations for the qqqqcc dibaryon system. Here, q = n, s and n = u, d.

Symmetry
Flavor combinations represent the irreducible base of the permutation group S n . It can help us to attach certain symmetry properties to the total wave function of the nnnncc dibaryon states. For the nnnncc dibaryon states with isospin I = 2, 1, and 0, the flavored states of four n-quarks can be represented by partitions [4], [3,1], and [2,2], respectively, and specific Young tableaux are given in Table III. Owing to the requirement of color confinement, the color part must be a singlet in the SU(3) symmetry. Thus, only the partition [2,2,2] is satisfied with the color singlet. The spin part can be represented by the partitions [6], [5,1], [4,2], and [3,3] for total spin of J = 3, 2, 1, and 0, respectively. Accordingly, we present the specific Young tableaux in Table III.
Once we have obtained the wave functions represented by the Young tableaux, the next step is to adopt a suitable coupling method to combine these wave functions in different spaces into fully antisymmetric wave functions. Thus, we need to introduce the corresponding Clebsch-Gordan (CG) coefficients of the permutation group S n , similar to the coupling method used in Refs. [28,[116][117][118]] Next, we construct the wave functions of the nnnscc (I = 3/2, 1/2) and sssncc states. Because the nnnscc (I = 3/2) states and sssncc states have exactly the same symmetry requirements, we only need to concentrate on the nnnscc (I = 3/2, 1/2) states. The s quark is far  (3) flavor symmetry breaking effect. We therefore include this effect by distinguishing the s quark from the u and d quarks, while still keeping the SU(2) isospin symmetry in our calculation. Accordingly, it is convenient to use the |[(nn)n][s(cc)] basis to construct the total wave functions of these states. We list all the possible flavours, color singlets, and spin wave functions in Table IV.
In Table IV, S(A) means fully symmetric (antisymmetric), and MS(MA) means that the first two quarks nn or the last two quarks cc are symmetric (antisymmetric) in nnnscc states. In the color part, we use the notation |[(n 1 n 2 ) color1 n 3 ] color2 [s 4 (c 5 c 6 ) color3 ] color4 to describe the color singlet wave functions, where color1, color2, color3, and color4 stand for the color representations of n 1 n 2 , n 1 n 2 n 3 , c 5 c 6 , and s 4 c 5 c 6 , respectively. Similarly, we use the notation |[(n 1 n 2 ) spin1 n 3 ] spin2 [s 4 (c 5 c 6 ) spin3 ] spin4 spin5 to describe the spin wave functions, where spin1, spin2, spin3, spin4, and spin5 represent the spins of n 1 n 2 , n 1 n 2 n 3 , c 5 c 6 , and s 4 c 5 c 6 and total spin, respectively.

C. Wave function for {12}{34}{56} symmetry
Finally, we construct the wave function of the nnsscc (I = 1, 0) dibaryon states. It is easy to see that there are three pairs of identical particles nn, ss, and cc. For these states, it is convenient to use the |(nn)(ss)(cc) basis to construct the nnsscc (I = 1, 0) dibaryon states. We list the corresponding flavours, color singlets, and spin wave functions in Table IV. In color part, we use the notation |(nn) color1 (ss) color2 (cc) color3 to describe the color-singlet wave functions, where color1, color2, and color3 stand for the color representations of n 1 n 2 , s 3 s 4 , and c 5 c 6 , respectively. Similarly, we use the notation |[(n 1 n 2 ) spin1 (s 3 s 4 ) spin2 ] spin3 (c 5 c 6 ) spin4 spin5 to describe the spin wave functions, where spin1, spin2, spin3, spin4, and spin5 represent the spins of n 1 n 2 , s 3 s 4 , n 1 n 2 s 3 s 4 and c 5 c 6 and the total spin, respectively. In Table IV, S(A) means totally symmetric(antisymmetric) for nn, ss, and cc of the discussed nnsscc states.

III. NUMERICAL RESULTS AND DISCUSSION
By constructing all the possible F f lavor ⊗ φ color ⊗ χ spin bases satisfied for {1234}{56}, {123}4{56}, and {12}{34}{56} symmetries and by introducing the parameters collected in Table I, we show the obtained corresponding mass spectra for the studied doubly charmed dibaryon states in Table V. Here, we use the notation H cc,xs (I, J P , Mass) to label a particular hexaquark state, where x(x = 0, 1, 2, 3, 4) represents the number of s quark.
A. The nnnncc states First, we examine the mass spectrum of the nnnncc (I = 2, 1, 0) states. Because the four n quarks are identical particles, we choose the wave functions in Sec. II A as the eigenvector to diagonalize the Hamiltonian matrix. According to the obtained mass spectrum in Table V, we plot the relative positions for the nnnncc dibaryon states in Fig. 1. In addition, by rearranging the six constituent quarks, we plot the possible thresholds of The masses of the ground conventional baryons are taken as the thresholds for these decay channels [112]. The nnnncc states can decay to the corresponding baryonbaryon final states via the S-wave or the D-wave. According to the conservation law in the strong decay process, the initial states and the final states have the same isospin I and spin-parity J P . For convenience, we label the spin (isospin) of the baryon-baryon states with a subscript (superscript) in Fig. 1. One of these states would be a good nnnncc dibaryon candidate if it were observed in the corresponding decay channels.
The positions of these nnnncc states are roughly shown in Fig. 1, and their properties may change accordingly if the masses of nnnncc states are determined by an observed state. Of course, the mass gaps between the different states are reliable.
Because of the complicated mixing between a large number of color-spin wave functions from the strong symmetric constraint, the total number of possible nnnncc states is limited to 17. There are 4 isotensor states, 8 isovector states, and 5 isoscalar states.
In Fig. 1, we can see that for the nnnncc states, the J P = 0 + particles have the smallest and the largest masses among all the quantum numbers. By labeling the nnnncc states for I =2, 1, and 0 with green, red, and blue lines, respectively, we can easily see that the I = 0 nnnncc hexaquarks generally have smaller masses than the I = 1 states, which have smaller masses than the From Fig. 1, there is no state with I(J P ) = 2(3 + ) owing to the constraint of the Pauli principle. Meanwhile, all states with I = 2 are above and beyond possible baryon-baryon thresholds; thus we think that there is no stable state and these states should have relatively large widths for the isotensor sector. In the isovector case, both the heaviest state and the lightest state occur in J P = 1 + . The mass gap between the two states is about 260 MeV. The H cc (1, 1 + , 5107) state can decay through all decay channels via the S-wave except for Λ c Λ c . Similarly, there is no I(J P ) = 0(3 + ) state, because of the symmetry constraint in the isoscalar sectors.
According to Fig. 1, we find that the lowest isoscalar state, i.e., the H cc (0, 0 + , 4577) state can only decay into Λ c Λ c and N Ξ cc final states. Furthermore, owing to the conservation of angular momentum and isospin, the H cc (0, 2 + , 4683) state decays mainly into N Ξ * cc via an Swave, but it can also decay into Λ c Λ c and N Ξ cc via the D-wave. Furthermore, because of the small phase space, this state is expected to be a narrow state.

B. The nnnscc states
Next, we investigate the nnnscc (I = 3/2, 1/2) dibaryon states. Because three n quarks are identical particles, we choose the wave functions in Sec. II B as the eigenvector to diagonalize the Hamiltonian matrix. In contrast to the wave functions of the nnnncc state, the color and spin wave functions are constructed in the "baryon ⊗ baryon" configuration. According to the mass spectrum of nnnscc states in Table V, we similarly plot the relative positions for nnnscc states and all the baryon-baryon thresholds of the associated rearrangement decay patterns in Fig. 2. The decays can occur via the S-or D-wave interactions, and each nnnscc state can decay into these channels from the parity conservation, isospin conversation, and angular momentum conservation. There are three types of rearrangement decays: (nnn + ccs), (nns + ccn), and (cnn + csn), respectively.
In Fig. 2, we label the nnnscc states for I = 3/2 and 1/2 with red and blue lines, respectively. We can easily find that the I = 1/2 states generally have smaller masses than the I = 3/2 nnnscc states. Thus, our results indicate that the states with smaller isospin quantum numbers generally form compact structures easily and have smaller masses. From Fig. 2, we find that the lowest I(J P ) = 1/2(1 + ) H cc,s (4603, 1/2, 1 + ) state is below all the rearrangement decay channels, and its decay via baryon-baryon channels is kinetically forbidden. We therefore consider it to be a relatively stable state. We label this stable state with " * " in Fig. 2 and Table V. In the nnnscc (I = 3/2) case, most states are expected to have relatively large widths, because they all have many different rearrangement decay channels. For example, consider the lowest I(J P ) = 3/2(2 + ) state: H cc,s (3/2, 2 + , 4940) can easily decay to ΣΞ * cc final states via the S wave. It also decays into Σ c Ξ c , ΣΞ cc final states by via the D wave.
As for the nnnscc states with I = 1/2, most of them are unstable particles against strong decay. Both the heaviest state and the lightest states occur in J P = 1 + . Here, the H cc,s (1/2, 0 + , 4755) and H cc,s (1/2, 1 + , 4771) states have the same three quark-rearrangement decay channels Λ c Ξ c , ΛΞ cc , and N Ω cc . Moreover, the H cc,s (1/2, 3 + , 5093) state has only one baryon-baryon decay channel Σ * Ξ * cc via the S wave. Meanwhile, it has many different decay channels via the D wave.

C. The nnsscc states
For dealing with the permutation symmetry of identical quarks n-n, s-s, and c-c, the wave functions are presented in Sec. II C. The flavor, color, and spin wave functions are constructed in the |(nn) ⊗ (ss) ⊗ (cc) configuration. According to the corresponding mass spectrum in Table V, we plot the relative positions for the nnsscc dibaryon states in Fig. 3 (a). Meanwhile, all the possible baryon-baryon thresholds of the associated rear- rangement decay patterns are shown in Fig. 3 (a). There are four types of rearrangement decays: (nns + scc), (nss + ncc), (nsc + nsc), and (nnc + ssc). By labeling the nnsscc states for I = 1 and 0 with red and blue lines, respectively, It is easy to see that the I = 0 states generally have smaller masses than the I = 1 states. Furthermore, it is easy to see that the two lowest nnccss dibaryon states, i.e., H cc,2s (4827, 0, 0 + ) and H cc,2s (4804, 0, 1 + ), should be stable, which are below all the baryon-baryon thresholds, and these two states are probably narrow.
For the isovector nnsscc states, the heaviest state and the lightest state are both the J P = 0 + states, i.e., H cc,2s (1, 0 + , 5550) and H cc,2s (1, 0 + , 4942). The mass gap between the two states is about 610 MeV. All the states are expected to be broad states with many different re-arrangement decay channels.
For the isoscalar nnsscc states, except for the relatively stable states, most of them have many different decay channels. Meanwhile, the H cc,2s (4892, 0, 1 + ) state is lower than the other thresholds and just above the ΛΩ cc threshold; thus it is probably narrow, and we can use the ΛΩ cc channel to identify this state. Furthermore, the lowest I(J P ) = 0(2 + ) state, i.e., H cc,2s (4904, 0, 2 + ), can only decay into ΛΩ cc via a D wave; thus, it should be narrow.

D. The sssncc and sssscc states
Finally, we discuss the sssncc and sssscc dibaryon states. The number of these allowed sssncc dibaryon states is the same as the number of nnnscc (I = 3/2) dibaryon states, owing to the same symmetry requirements according to Table V. Similarly, the number of allowed sssscc dibaryon states is the same as the number of nnnncc (I = 2) dibaryon states. Furthermore, the relative positions for the sssncc and sssscc dibaryon states and their all possible rearrangement decay channels are shown in Fig. 3 (b) and (c), respectively.
The sssncc dibaryon states have 12 possible rearrangement decay channels: cc , ΩΞ cc , Ω * cc Ξ * , Ω * cc Ξ, Ω cc Ξ * , and Ω cc Ξ. The sssscc dibaryon states, they have 5 possible rearrangement decay channels: Ω * c Ω * c , ΩΩ * cc , Ω c Ω * c , ΩΩ cc , and Ω c Ω c . From Fig. 3 (c), we notice that the ground sssscc dibaryon state with quantum number I(J P ) = 0(3 + ) does not exist, owing to the constraint of the Pauli principle. From Fig. 3 (b), all the states are above some possible baryon-baryon thresholds in the sssncc dibaryon state. Thus, they should have relatively large widths. As expected, all the sssscc dibaryon states are above the allowed baryon-baryon thresholds from Fig. 3 (c). They can easily decay into the rearrangement channels owing to the large phase space. Accordingly, all the sssscc dibaryon states are likely to be broad. It does not seem easy to find them experimentally.

IV. SUMMARY
In 2021, the LHCb Collaboration reported the first doubly charmed tetraquark state T + cc (ccūd). Inspired by this discovery, the existence of doubly heavy tetraquark states or pentaquark states with different quark configurations has been studied via various theoretical calculations. According to the similarity between the ccūd tetraquark state and the Ξ ++ cc baryon (ccu), we naturally ask whether the doubly charmed dibaryon states exist. Simultaneously, it is an interesting question whether the dibaryon state with double heavy quarks is also stable against strong decay. In this work, we first introduce the CMI model and construct the F f lavor ⊗ φ color ⊗ χ spin wave functions for the doubly charmed dibaryon system. The main challenge here is how to apply the Pauli principle to the permutation of the three or four identical particles. Therefore, we introduce the Young tableau to represent the irreducible bases of the permutation group. Then, we systematically calculate the corresponding Hamiltonian. Furthermore, we obtain the mass spectra of the dibaryon states nnnncc, nnnscc, nnsscc, sssncc, and sssscc, which are listed in Table V. Accordingly, the relative positions between dibaryon states and possible baryon-baryon thresholds are plotted in Figs. 1-3. Through comparison with baryon-baryon thresholds, we can analyze the stability of dibaryon states and the decay properties of the possible quark rearrangement decay channels. Furthermore, the mass gaps between the relevant states are reliable, and if a qqqqcc dibaryon state is observed, we can use these mass gaps to predict its corresponding multiples.
Hence, we can explore the presence of the three "sta-
Of necessity, owing to the uncertainty of the CMI model, further dynamical calculations are needed to clarify their nature. We hope that the present study will inspire relevant experiments to search for these states in the future.