-
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 charmoniumlike
XYZ states,Pc states, and the double charmT+cc state have been reported experimentally, stimulating extensive studies on exotic hadron spectroscopy (see review articles [1–5] for details). Recently, the LHCb Collaboration announced new observations of theTacˉs0(2900)++,0 state in the processesB0→ˉD0D+sπ− andB−→ˉD−D+sπ+ [6, 7] and aPΛψs(4338) state inB−→J/ψΛˉp [8]. In addition, the CMS Collaboration confirmed the observation of theX(6900) state and found more structures in the di-J/ψ invariant mass spectrum [9], and the ATALAS Collaboration confirmed the discovery of thePc states andX(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 isd∗(2380) . Dyson and Xuong first predicted a well-known candidateΔΔ resonance state according to the SU(6) symmetry of the strong interaction [13], after proposing the quark model [14]. In 1989, Goldman et al. proposed “an inevitable nonstrange dibaryon” namedd∗ because of its unique symmetry features. About 20 years later, thed∗(2380) resonance with theM=2.37 GeV,Γ=70 MeV, andI(JP)=0(3+) was reported by the WASA-at-COSY and CELSIUS/WASA Collaborations [15]. It comes from the exclusive reaction channelpn→dπ0π0 , and its existence was supported by later experiments [16–19]. In addition, many theoretical papers properly described the properties of this resonance [20–32].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–39]. In addition to these studies, this state has been extensively researched using various models [40–47]. Another interesting candidate was introduced later by Goldman et al. [48], which is expected to appear as aNΩ 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 theNΩ 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 theNΩ state [50]. J. Leandri and B. Silvestre-Brac conducted a series of investigations on dibaryons in the framework of a pure chromomagnetic Hamiltonian [51–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 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 the1S0 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–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 (
q6 ) states and the three-quark-three-antiquark (q3ˉq3 ) states was discussed, and the investigations were extended to the (q3Q3 ) and (Q3ˉq3 ) systems. Later, Huang et al. investigated possibleNΩ -like dibaryonsNΩccc andNΩbbb with quantum numbersI(JP)=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 the light-quark structures of theˉD(∗) and theΞ(∗)cc to be almost identical. Thus, theoretical work has also focused on the triply heavy dibaryons [68–72].In particular, inspired by the recent discovery of a narrow structure
T+cc (ccˉuˉd ) in the LHCb experiment [73, 74], on the basis of the similarity between theccˉuˉd configuration and theΞ++cc baryon (ccu ) [73], the existence of other open charm multiquarks has been studied via the extended chromomagnetic interaction (CMI) model [75–77], the molecular state picture [78–80], the Born-Oppenheimer approximation [81], and the chiral quark model [82].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–85], the quark delocalization color screening model [86–88], the lattice QCD calculation [89], the chiral constituent quark model [90], the nucleon-nucleon interaction model [91, 92], the chiral effective field theory (EFT) [93–95], the QCD sum rules [96, 97], the relativistic six-quark equations [98], the complex scaling method [99], and the one-boson-exchange (OBE) model [100–105]. 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 ) [106]. 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 [107], 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
Fflavor⊗ϕcolor⊗χspin wave functions for the S-waveqqqqcc dibaryon state. The mass spectrum and decay mode of theqqqqcc hexaquark state with different constituent quarks are discussed in Sec. III. Finally, a short summary is given in Sec. IV. -
In the chromomagnetic interaction (CMI) model, the mass of the ground hadron state can be described by the effective Hamiltonian [64, 75–77, 108–111]
H=∑imi+HCI+HCMI=∑imi−∑i<jAij→λi⋅→λj−∑i<jvij→λi⋅→λj→σi⋅→σj=−34∑i<jmijVCij−∑i<jvijVCMIij,
(1) where
VCij=→λi⋅→λj andVCMIij=→λi⋅→λj→σi⋅→σj represent the color and chromomagnetic interactions between the quarks, respectively, withσi andλi being the Pauli matrices and Gell-Mann matrices, respectively. The mass parameter of the quark pairmij=1/4(mi+mj)+4/3Aij contains the effective quark massmi (mj ) and the color interaction strengthAij . Thevij 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
mij andvij from the conventional baryon masses, which are presented in Table 1.Parameter mnn mns mss mnc msc mcc Value 181.2 226.7 262.3 520.0 545.9 792.9 Parameter vnn vns vss vnc vsc vcc Value 19.1 13.3 12.2 3.9 4.4 3.5 Table 1. Relevant coupling parameters for the doubly charmed dibaryon system [111] (in units of MeV).
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.,
Φtotal=Ψspace⊗Fflavor⊗ϕcolor⊗χspin.
(2) According to the Pauli principle, this wave function should be completely antisymmetric when identical quarks are exchanged. Here, we consider only the low-lying S-wave dibaryon states, so that their spatial wave functions are symmetric and become trivial. Thus, the
Fflavor⊗ϕ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 arennnncc (I=2,1,0 ),nnnscc (I=3/2,1/2 ),nnsscc (I=1,0 ),sssncc , andsssscc . According to the symmetry properties, we can divide theqqqqcc dibaryon system into three groups, which are presented in Table 2. Here, we have to construct the correspondingFflavor⊗ϕ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.Symmetry Flavor combinations {1234}{56} nnnncc (I=2,1,0 )sssscc {123}4{56} nnnscc (I=3/2,1/2 )sssncc {12}{34}{56} nnsscc (I=1,0 )Table 2. All possible flavor combinations for the
qqqqcc dibaryon system. Here,q=n , s andn=u , d.In previous works, there are abundant methods on the wave functions of six-quark systems [29, 40, 83, 112–114]. 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, 115, 116].
-
First, we construct the wave function of the
nnnncc (I=2,1,0 ) andsssscc dibaryon states. Because thennnncc (I=2 ) dibaryon states andsssscc dibaryon states have exactly the same symmetry requirements, we only need to concentrate on thennnncc (I=2,1,0 ) dibaryon states. Here, we introduce the Young tableau to represent the irreducible base of the permutation groupSn . It can help us to attach certain symmetry properties to the total wave function of thennnncc dibaryon states. For thennnncc dibaryon states with isospinI=2,1, and0 , 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 3. 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 ofJ=3,2,1, and0 , respectively. Accordingly, we present the specific Young tableaux in Table 3.Table 3. The flavor, color, and spin wave functions represented by the partition and the corresponding Young tableaux. Here, the
I/J represents the isospin/spin.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
Sn , similar to the coupling method used in Refs. [28, 115–117]|[f]Y⟩=∑Y′,Y″S([[f′]Y′[f″]Y″]|[f]Y)|[f′]Y′⟩|[f″]Y″⟩.
(3) Here,
[f] ,[f′] , or[f″] denotes an irreducible representation ofSn , Y denotes a Young tableau, and|[f]Y⟩ is a basis vector ofSn .S([[f′]Y′[f″]Y″]|[f]Y)|[f′]Y′⟩ are CG coefficients. As mentioned above, the symmetry of the color wave function is fixed as[2,2,2] . Thus, we adopt a coupling method whereby one can first combine the color-SU(3) singlet[2,2,2] with spin-SU(2) representation[f]S into an SU(6)[f]CS and then combine it with isospin-SU(2) representations[f]I to obtain SU(12) representations[f]CSI . Thus, we can calculate the 6-quark wave function with the desired permutation symmetry[f] as a linear combination ofFflavor⊗ϕcolor⊗χspin . -
Next, we construct the wave functions of the
nnnscc (I=3/2,1/2 ) andsssncc states. Because thennnscc (I=3/2 ) states andsssncc states have exactly the same symmetry requirements, we only need to concentrate on thennnscc (I=3/2,1/2 ) states. The s quark is far heavier than the u and d quarks, leading to an SU(3) flavor symmetry breaking effect. We 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 4.nnnscc statesFlavor I=32 FS,S=|[(nn)I=1n]I=3/2[s(cc)]⟩ I=12 FMS,S=|[(nn)I=1n]I=1/2[s(cc)]⟩ FMA,S=|[(nn)I=0n]I=1/2[s(cc)]⟩ Color ϕMA,MA=|[(nn)ˉ3n]8[s(cc)ˉ3]8⟩ ϕMA,MS=|[(nn)ˉ3n]8[s(cc)6]8⟩ ϕA,A=|[(nn)ˉ3n]1[s(cc)ˉ3]1⟩ ϕMS,MA=|[(nn)6n]8[s(cc)ˉ3]8⟩ ϕMS,MS=|[(nn)6n]8[s(cc)6]8⟩ Spin J=3 χS,S1=|[(nn)1n]3/2[s(cc)1]3/2⟩3 J=2 χS,S2=|[(nn)1n]3/2[s(cc)1]3/2⟩2 χS,MS3=|[(nn)1n]3/2[s(cc)1]1/2⟩2 χMS,S4=|[(nn)1n]1/2[s(cc)1]3/2⟩2 χMS,A5=|[(nn)1n]3/2[s(cc)0]3/2⟩2 χMA,S6=|[(nn)0n]1/2[s(cc)1]3/2⟩2 J=1 χS,S7=|[(nn)1n]3/2[s(cc)1]3/2⟩1 χS,MS8=|[(nn)1n]3/2[s(cc)1]1/2⟩1 χMS,S9=|[(nn)1n]1/2[s(cc)1]3/2⟩1 χMS,MS10=|[(nn)1n]1/2[s(cc)1]1/2⟩1 χS,MA11=|[(nn)1n]3/2[s(cc)0]1/2⟩1 χMA,S12=|[(nn)0n]1/2[s(cc)1]3/2⟩1 χMS,MA13=|[(nn)1n]1/2[s(cc)0]1/2⟩1 χMA,MS14=|[(nn)0n]1/2[s(cc)1]1/2⟩1 χMA,MA15=|[(nn)0n]1/2[s(cc)0]1/2⟩1 J=0 χS,S16=|[(nn)1n]3/2[s(cc)1]3/2⟩0 χMS,MS17=|[(nn)1n]1/2[s(cc)1]1/2⟩0 χMS,MA18=|[(nn)1n]1/2[s(cc)0]1/2⟩0 χMA,MS19=|[(nn)0n]1/2[s(cc)1]1/2⟩1 χMA,MA20=|[(nn)0n]1/2[s(cc)0]1/2⟩0 F⊗ϕ⊗χ I=32 FS,S⊗ϕA,A⊗χS,S ;FS,S⊗ϕA,A⊗χS,MS ;1√2[FS,S⊗(ϕMS,MA⊗χMA,MS−ϕMA,MA⊗χMS,S)] 1√2[FS,S⊗(ϕMS,MA⊗χMA,S−ϕMA,MA⊗χMS,S)] ;1√2[FS,S⊗(ϕMS,MS⊗χMA,MA−ϕMA,MS⊗χMS,MA)] 1√2[FS,S⊗(ϕMS,MA⊗χMA,S−ϕMA,MA⊗χMS,MS)] ;1√2[FS,S⊗(ϕMS,MA⊗χMA,MS−ϕMA,MA⊗χMS,MS)] I=12 1√2[(FMS,S⊗ϕMA,MA−FMA,S⊗ϕMS,MA)⊗χS,S)] ;1√2(FMS,S⊗ϕMA,MA⊗χS,S−FMA,S⊗ϕMS,MS⊗χS,MA) 1√2[(FMS,S⊗ϕMA,MA−FMA,S⊗ϕMS,MA)⊗χS,MS)] ;1√2(FMS,S⊗ϕMA,MS⊗χS,MA−FMA,S⊗ϕMS,MA⊗χS,MS) 1√2[(FMS,S⊗ϕMA,MS−FMA,S⊗ϕMA,MS)⊗χS,MA)] ;1√2(FMS,S⊗ϕMA,MA⊗χS,MS−FMA,S⊗ϕMS,MS⊗χS,MA) 1/2[(FMS,S⊗ϕMS,MS−FMA,S⊗ϕMA,MS)⊗χMA,MA)+(FMS,S⊗ϕMA,MS+FMA,S⊗ϕMS,MS)⊗χMS,MA)] 1/2[(FMS,S⊗ϕMS,MA−FMA,S⊗ϕMA,MA)⊗χMA,MS)+(FMS,S⊗ϕMA,MS+FMA,S⊗ϕMS,MS)⊗χMS,MA)] 1/2[(FMS,S⊗ϕMS,MA−FMA,S⊗ϕMA,MA)⊗χMA,MS)+(FMS,S⊗ϕMA,MA+FMA,S⊗ϕMS,MA)⊗χMS,MS)] 1/2[(FMS,S⊗ϕMS,MS−FMA,S⊗ϕMA,MS)⊗χMA,MA)+(FMS,S⊗ϕMA,MA+FMA,S⊗ϕMS,MA)⊗χMS,MS)] 1/2[(FMS,S⊗ϕMS,MA−FMA,S⊗ϕMA,MA)⊗χMA,S)+(FMS,S⊗ϕMA,MS+FMA,S⊗ϕMS,MS)⊗χMS,MA)] 1/2[(FMS,S⊗ϕMS,MA−FMA,S⊗ϕMA,MA)⊗χMA,S)+(FMS,S⊗ϕMA,MA+FMA,S⊗ϕMS,MA)⊗χMS,MS)] 1/2[(FMS,S⊗ϕMS,MA−FMA,S⊗ϕMA,MA)⊗χMA,MS)+(FMS,S⊗ϕMA,MA+FMA,S⊗ϕMS,MA)⊗χMS,S)] 1/2[(FMS,S⊗ϕMS,MS−FMA,S⊗ϕMA,MS)⊗χMA,MA)+(FMS,S⊗ϕMA,MA+FMA,S⊗ϕMS,MA)⊗χMS,S)] 1/2[(FMS,S⊗ϕMS,MA−FMA,S⊗ϕMA,MA)⊗χMA,S)+(FMS,S⊗ϕMA,MA+FMA,S⊗ϕMS,MA)⊗χMS,S)] nnsscc statesFlavor I=1 FS,S,S=|(nn)I=1(ss)(cc)⟩ I=0 FA,S,S=|(nn)I=0(ss)(cc)⟩ Color ϕA,A,A=|(nn)ˉ3(ss)ˉ3(cc)ˉ3⟩ ϕA,S,A=|(nn)ˉ3(ss)6(cc)ˉ3⟩ ϕA,A,S=|(nn)ˉ3(ss)ˉ3(cc)6⟩ ϕS,A,A=|(nn)6(ss)ˉ3(cc)ˉ3⟩ ϕS,S,S=|(nn)6(ss)6(cc)6⟩ Spin J=3 χS,S,S1=|[(nn)1(ss)1]2(cc)1⟩3 J=2 χS,S,S2=|[(nn)1(ss)1]2(cc)1⟩2 χS,S,S3=|[(nn)1(ss)1]1(cc)1⟩2 χS,A,S4=|[(nn)1(ss)0]1(cc)1⟩2 χS,S,A5=|[(nn)1(ss)1]2(cc)0⟩2 χA,S,S6=|[(nn)0(ss)1]1(cc)1⟩2 J=1 χS,S,S7=|[(nn)1(ss)1]2(cc)1⟩1 χS,S,S8=|[(nn)1(ss)1]1(cc)1⟩1 χS,S,S9=|[(nn)1(ss)1]0(cc)1⟩1 χS,A,S11=|[(nn)1(ss)0]1(cc)1⟩1 χA,S,S12=|[(nn)0(ss)1]1(cc)1⟩1 χS,A,A13=|[(nn)1(ss)0]1(cc)0⟩1 χS,A,A13=|[(nn)1(ss)0]1(cc)0⟩1 χA,S,A14=|[(nn)0(ss)1]1(cc)0⟩1 χA,A,S15=|[(nn)0(ss)0]0(cc)1⟩1 J=0 χS,S,S16=|[(nn)1(ss)1]1(cc)1⟩0 χS,A,S17=|[(nn)1(ss)0]1(cc)1⟩0 χS,S,S18=|[(nn)1(ss)1]0(cc)0⟩0 χA,S,S19=|[(nn)0(ss)1]1(cc)1⟩0 χA,A,A20=|[(nn)0(ss)0]0(cc)0⟩0 F⊗ϕ⊗χ I=1 FS,S,S⊗ϕA,A,A⊗χS,S,S ;FS,S,S⊗ϕA,A,S⊗χS,S,A ;FS,S,S⊗ϕA,S,A⊗χS,A,S ;FS,S,S⊗ϕS,A,A⊗χA,S,S ;FS,S,S⊗ϕS,S,S⊗χA,A,A I=0 FA,S,S⊗ϕS,S,S⊗χS,A,A ;FA,S,S⊗ϕA,A,A⊗χA,S,S ;FA,S,S⊗ϕA,A,S⊗χA,S,A ;FA,S,S⊗ϕA,S,A⊗χA,A,S ;FA,S,S⊗ϕS,A,A⊗χS,S,S Table 4. All the possible flavors, colors, and spin wave functions for the
nnnscc andnnsscc states.F⊗ϕ⊗χ represents the combination of the wave functions in different spaces. Here,I/J represents the isospin/spin.In Table 4, S(A) means fully symmetric (antisymmetric), and MS (MA) means that the first two quarks
nn or the last two quarkscc are symmetric (antisymmetric) innnnscc states. In the color part, we use the notation|[(n1n2)color1n3]color2[s4(c5c6)color3]color4⟩ to describe the color singlet wave functions, where color1, color2, color3, and color4 stand for the color representations ofn1n2 ,n1n2n3 ,c5c6 , ands4c5c6 , respectively. Similarly, we use the notation|[(n1n2)spin1n3]spin2[s4(c5c6)spin3]spin4⟩spin5 to describe the spin wave functions, where spin1, spin2, spin3, spin4, and spin5 represent the spins ofn1n2 ,n1n2n3 ,c5c6 , ands4c5c6 and the total spin, respectively. -
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 particlesnn ,ss , andcc . For these states, it is convenient to use the|(nn)(ss)(cc)⟩ basis to construct thennsscc (I=1,0 ) dibaryon states. We list the corresponding flavors, color singlets, and spin wave functions in Table 4. 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 ofn1n2 ,s3s4 , andc5c6 , respectively. Similarly, we use the notation|[(n1n2)spin1(s3s4)spin2]spin3(c5c6)spin4⟩spin5 to describe the spin wave functions, where spin1, spin2, spin3, spin4, and spin5 represent the spins ofn1n2 ,s3s4 ,n1n2s3s4 , andc5c6 and the total spin, respectively. In Table 4,S(A) means totally symmetric (antisymmetric) fornn ,ss , andcc of the discussednnsscc states. -
By constructing all the possible
Fflavor⊗ϕcolor⊗χspin bases satisfied for{1234}{56} ,{123}4{56} , and{12}{34}{56} symmetries and by introducing the parameters collected in Table 1, we show the obtained corresponding mass spectra for the studied doubly charmed dibaryon states in Table 5. Here, we use the notationHcc,xs(I,JP,Mass) to label a particular hexaquark state, wherex(x=0,1,2,3,4) represents the number of s quarks.Table 5. The masses for the dibaryon states
nnnncc ,nnnscc ,nnsscc ,sssncc , andsssscc in units of MeV. I, J, and P represent the isospin, angular momentum, and parity, respectively. -
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.C as the eigenvector to diagonalize the Hamiltonian matrix. According to the obtained mass spectrum in Table 5, we plot the relative positions for thennnncc dibaryon states in Fig. 1. In addition, by rearranging the six constituent quarks, we plot the possible thresholds of the baryon-baryon decay patterns for thennnncc states in Fig. 1. There are two rearrangement decay types: (single charmed baryon + single charmed baryon) and (light baryon + doubly charmed baryon). Here, all the possible decay channels in the figure areΛcΛc ,ΛcΣc ,ΛcΣ∗c ,ΣcΣc ,Σ∗cΣc ,Σ∗cΣ∗c andNΞcc ,NΞ∗cc ,ΔΞcc ,ΔΞ∗cc . The masses of the ground conventional baryons are taken as the thresholds for these decay channels [111]. Thennnncc states can decay to the corresponding baryon-baryon 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-parityJP . 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 goodnnnncc dibaryon candidate if it were observed in the corresponding decay channels.Figure 1. (color online) Relative positions (units: MeV) for the
nnnncc dibaryon states labeled with solid lines. The green, red, and blue lines represent thennnncc dibaryon states withI=2 ,I=1 , andI=0 , respectively. The dotted lines denote different baryon-baryon thresholds, and the superscripts (subscripts) of the labels represent the possible total angular momentum (isospin) of the channels.The positions of these
nnnncc states are roughly shown in Fig. 1, and their properties may change accordingly if the masses ofnnnncc 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, theJP=0+ particles have the smallest and the largest masses among all the quantum numbers. By labeling thennnncc states forI=2 , 1, and 0 with green, red, and blue lines, respectively, we can easily see that theI=0 nnnncc hexaquarks generally have smaller masses than theI=1 states, which have smaller masses than theI=2 ones. Thus, our results indicate that the states with smaller isospin quantum numbers generally form compact structures easily and have smaller masses.From Fig. 1, there is no state with
I(JP)=2(3+) owing to the constraint of the Pauli principle. Meanwhile, all states withI=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 inJP=1+ . The mass gap between the two states is about 260 MeV. TheHcc(1,1+,5107) state can decay through all decay channels via the S-wave except forΛcΛc . Similarly, there is noI(JP)=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
Hcc(0,0+,4577) state, can only decay intoΛcΛc andNΞcc final states. Furthermore, owing to the conservation of angular momentum and isospin, theHcc(0,2+,4683) state decays mainly intoNΞ∗cc via an S-wave, but it can also decay intoΛcΛc andNΞcc via the D-wave. Furthermore, because of the small phase space, this state is expected to be a narrow state. -
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 thennnncc state, the color and spin wave functions are constructed in the “baryon⊗baryon ” configuration. According to the mass spectrum ofnnnscc states in Table 5, we similarly plot the relative positions fornnnscc 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 eachnnnscc 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.Figure 2. (color online) Relative positions (units: MeV) for the
nnnscc dibaryon states. The red lines and the blue lines show thennnscc dibaryon states withI=3/2 andI=1/2 , respectively. The dotted lines denote different baryon-baryon thresholds, and the superscripts (subscripts) of the labels represent the possible total angular momentum (isospin) of the channels.In Fig. 2, we label the
nnnscc states forI=3/2 and1/2 with red and blue lines, respectively. We can easily find that theI=1/2 states generally have smaller masses than theI=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 lowestI(JP)=1/2(1+)Hcc,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 5.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 lowestI(JP)=3/2(2+) state:Hcc,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 via the D wave.As for the
nnnscc states withI=1/2 , most of them are unstable particles against strong decay. Both the heaviest state and the lightest states occur inJP=1+ . Here, theHcc,s(1/2,0+,4755) andHcc,s(1/2,1+,4771) states have the same three quark-rearrangement decay channelsΛcΞc ,ΛΞcc , andNΩcc . Moreover, theHcc,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. -
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 5, we plot the relative positions for thennsscc dibaryon states in Fig. 3 (a). Meanwhile, all the possible baryon-baryon thresholds of the associated rearrangement decay patterns are shown in Fig. 3 (a). There are four types of rearrangement decays: (nns+scc ), (nss+ncc ), (nsc+nsc ), and (nnc+ssc ).Figure 3. (color online) Relative positions (units: MeV) for three types of hexaquark states. (a) The
nnsscc subsystem: the red lines and the blue lines show thennsscc hexaquark states withI=1 andI=0 respectively; (b)sssncc subsystem; (c)sssscc subsystem. The dotted lines denote different baryon-baryon thresholds, and the superscripts (subscripts) of the labels represent the possible total angular momentum (isospin) of the channels.By labeling the
nnsscc states forI=1 and0 with red and blue lines, respectively, It is easy to see that theI=0 states generally have smaller masses than theI=1 states. Furthermore, it is easy to see that the two lowestnnccss dibaryon states, i.e.,Hcc,2s(4827,0,0+) andHcc,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 theJP=0+ states, i.e.,Hcc,2s(1,0+,5550) andHcc,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 rearrangement decay channels.For the isoscalar
nnsscc states, except for the relatively stable states, most of them have many different decay channels. Meanwhile, theHcc,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 lowestI(JP)=0(2+) state, i.e.,Hcc,2s(4904,0,2+) , can only decay intoΛΩcc via a D wave; thus, it should be narrow. -
Finally, we discuss the
sssncc andsssscc dibaryon states. The number of these allowedsssncc dibaryon states is the same as the number ofnnnscc (I=3/2) dibaryon states, owing to the same symmetry requirements according to Table 5. Similarly, the number of allowedsssscc dibaryon states is the same as the number ofnnnncc (I=2) dibaryon states. Furthermore, the relative positions for thesssncc andsssscc dibaryon states and their all possible rearrangement decay channels are shown in Figs. 3(b) and (c), respectively.The
sssncc dibaryon states have 12 possible rearrangement decay channels:Ω∗cΞ∗c ,Ω∗cΞ′c ,Ω∗cΞc ,ΩcΞ∗c ,ΩcΞ′c ,ΩcΞc ,ΩΞ∗cc ,ΩΞcc ,Ω∗ccΞ∗ ,Ω∗ccΞ ,ΩccΞ∗ , andΩccΞ . Thesssscc dibaryon states have five possible rearrangement decay channels:Ω∗cΩ∗c ,ΩΩ∗cc ,ΩcΩ∗c ,ΩΩcc , andΩcΩc . From Fig. 3(c), we notice that the groundsssscc dibaryon state with quantum numberI(JP)=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 thesssscc 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 thesssscc dibaryon states are likely to be broad. It does not seem easy to find them experimentally. -
In 2021, the LHCb Collaboration reported the first doubly charmed tetraquark state
T+cc (ccˉuˉ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 theccˉuˉ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
Fflavor⊗ϕ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 statesnnnncc ,nnnscc ,nnsscc ,sssncc , andsssscc , which are listed in Table 5. 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 aqqqqcc dibaryon state is observed, we can use these mass gaps to predict its corresponding multiples.According to the obtained mass spectrum, there is no state with
I(JP)=2(3+) and0(3+) among thennnncc andsssscc states, owing to the constraint of the Pauli principle. Moreover, our results show that the states with smaller isospin quantum numbers generally form a compact structure easily and have smaller masses in thennnncc (I=2,1,0 ),nnnscc (I=3/2,1/2 ), andnnsscc (I=1,0 ) states. Meanwhile, we find four relatively "stable" states:Hcc,s(1/2,1+,4603) ,Hcc,2s(0,0+,4827) , andHcc,2s(0,1+,4804) , whose decays through baryon-baryon channels are kinetically forbidden.Hence, we can explore the presence of the three "stable" double charm hexaquark states by investigating their weak decay channels. We assume that only one charm quark decays into
dˉdu ,sˉdu , orsˉsu [118–122], and all quarks can form various hadron final states through quark rearrangement. However, the weak decay process must obey the laws of conservation of charge and angular momentum. In Table 6, we list some possible weak decay modes for the three "stable" states, focusing on electric hadron final states, which could aid in experimental reconstruction of the hexaquark state. For instance, innnnscc , theHcc,s(1/2,1+,4603) states carry one or two units of positive charges and may decay intoΣ+cpK−(π−) ,Ξ+cΣ−ρ+(π+) ,Λ+cΣ+K− ,Σ+cΣ+K− ,Ξ+cpK+K− ,Ξ+cΔ++K− , etc. Similarly, the two statesHcc,2s(0,0+,4827) andHcc,2s(0,1+,4804) share common decay modes such asΞ+cΣ+K− ,Ξ+cΞ−π+ ,Ξ+cΞ−ρ+(π+) , andΛ+cΩ−ρ+(π+) . The weak decay modes presented in Table 6 correspond to final states containing two electric baryons and one electric meson. Theoretically, more electric or neutral mesons could appear in the decay final states.Hexaquark states Weak decay modes Hcc,s(1/2,1+,4603) Σ+cpK−(π−) ,Ξ+cΣ−ρ+(π+) ,Λ+cΣ+K− ,Σ+cΣ+K− ,Ξ+cpK+K− ,Ξ+cΔ++K− ,......Hcc,2s(0,0+,4827) Ξ+cΣ+K− ,Ξ+cΞ−π+ ,Ξ+cΞ−ρ+(π+) ,Λ+cΩ−ρ+(π+) ,......Hcc,2s(0,1+,4804) Table 6. Some possible weak decay modes of the three stable hexaquark states.
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.
Doubly charmed dibaryon states
- Received Date: 2023-03-07
- Available Online: 2023-06-15
Abstract: In this work, we study the doubly charmed dibaryon states with the