Doubly-charmed hexaquarks in the diquark picture

Abstract

We investigate doubly-charmed hexaquark states within the diquark picture, by employing the constituent quark model and the quark-interchange model as our theoretical frameworks. Using the Gaussian expansion method, we systematically study these states, with calculating various properties such as mass spectra, internal contributions of each Hamiltonian component, root-mean-square radii, and two-body strong decay widths. Our analysis of the mass spectra reveals no stable state in this system. Furthermore, the root-mean-square radii suggest that the doubly-charmed hexaquark states exhibit a compact configuration. By examining the decay widths, we identify potentially detectable states and their primary decay channels within each subsystem. Despite the large decay phase space, we still find narrow states with total widths of less than 10 MeV. This study provides a theoretical foundation for understanding the structures and interactions of doubly-charmed hexaquark states and offers valuable insights for future experimental searches.

I Introduction

In the field of hadron physics, the study of exotic states beyond conventional mesons ($q\bar{q}$) and baryons ($qqq$) has been a central research focus Chen:2016qju ; Guo:2017jvc ; Liu:2019zoy ; Hosaka:2016pey ; Jaffe:2004ph . A pivotal breakthrough occurred in 2003 with the Belle collaboration’s discovery of the $X(3872)$ Belle:2003nnu , ushering in an era of rapid progress. Subsequent experiments by Collaborations such as LHCb, BaBar, Belle, and BESIII have unveiled a spectrum of unconventional hadronic states, including the $P_{c}/P_{cs}$ pentaquarks LHCb:2015yax ; LHCb:2019kea ; LHCb:2020jpq ; LHCb:2022ogu and charmonium-like $XYZ$ states Belle:2004lle ; LHCb:2022aki ; LHCb:2021uow ; BESIII:2016bnd ; BESIII:2013ris ; BESIII:2013qmu . These discoveries have motivated extensive theoretical efforts to elucidate their internal structures, exploring configurations such as compact tetraquarks, molecular states, hybrids, and glueballs Chen:2022asf ; Dong:2021juy ; Brambilla:2019esw .

Although the hidden-charm sector has dominated this landscape, recent observations of the doubly-charmed baryon $\Xi^{++}_{cc}(3621)$ LHCb:2017iph ; LHCb:2018pcs and the tetraquark candidate $T_{cc}^{+}(3875)$ LHCb:2021vvq ; LHCb:2021auc by LHCb have shifted attention to the double-charm sector. This raises intriguing questions: Could doubly-charmed pentaquarks exist, analogous to their hidden-charm counterparts? Moreover, what about doubly-charmed hexaquarks? The latter, as a six-quark system, brings computational challenges due to the complexity of its inherent six-body problem. Current theoretical studies often adopt approximation methods. For example, Liu et al. investigated the mass spectra of these systems using the chromomagnetic interaction (CMI) model Liu:2022rzu . Complementary to this approach, the one-boson-exchange potential model has been employed to examine hadron-level interactions and predict the properties of doubly-charmed hexaquarks within the framework of two-body molecular states, as demonstrated in Refs. Meguro:2011nr ; Cheng:2022vgy ; Shah:2024thr ; Dong:2021bvy ; Qi:2024dqz ; Vijande:2016nzk ; Andreev:2024orz . Further insights have emerged from exploratory lattice QCD calculations, which address the binding characteristics of these systems Geng:2024dpk . While these efforts provide valuable insights, the reliance on diverse approximations underscores the necessity for novel methodologies to resolve persistent theoretical ambiguities and deepen our understanding of hexaquark dynamics.

A promising framework to simplify such systems is the diquark model  Ida:1966ev ; Anselmino:1992vg ; Barabanov:2020jvn , proven effective in hadron spectroscopy. For example, many observed charm baryons are bound states of a charm quark and a light diquark ($qq$), with $\lambda$-mode excitations explaining their spectra  Chen:2019ywy ; Chen:2021eyk ; Chen:2014nyo ; Chen:2009tm . The model has shown unique value in studying multi-quark states with other configurations—for instance, in the compact diquark-anti-diquark (diquark+diquark+antiquark) configuration, it offers an effective way to investigate hidden-charm and hidden-bottom $XYZ$ exotic states Zhu:2016arf ; Shi:2021jyr ; Bicudo:2015vta (like $P_{c}$, $P_{cs}$ states Ali:2019clg ; Lebed:2015tna ). When extended to the doubly-charmed hexaquark system (a six-body system), it can be approximated as a four-body system of two charm quarks and two light diquarks. By “freezing” the diquarks’ internal degrees of freedom, this simplification reduces computational complexity while preserving key physical features, akin to charm baryon treatments.

Notably, this framework suggests a structural parallel between the $T_{cc}^{+}(3875)$ tetraquark (composed of two charm quarks and two light antiquarks) and a hypothetical doubly-charmed hexaquark. Replacing two light antiquarks $\bar{q}\bar{q}$ in $T_{cc}^{+}(3875)$ with two light diquarks $[qq][qq]$ would yield a hexaquark configuration ($ccqqqq$), potentially mirroring properties of its tetraquark counterpart. Establishing such connections theoretically reveals universal features of multiquark systems, including potential shared properties like binding mechanisms and decay modes. Experimentally, it provides guidance for detecting elusive hexaquarks by leveraging known tetraquark signatures.

In this work, we systematically investigate the doubly-charmed hexaquark system $ccqqqq$ ($q=n,~{}s;~{}n=u,~{}d$) using the Gaussian expansion method within the framework of the constituent quark model. Our calculations include: the mass spectra for all flavor configurations, corresponding internal mass contributions, root-mean-square (RMS) radii, and partial (total) decay widths.

The paper is organized as follows. After the Introduction, Section II details the theoretical framework, including the effective Hamiltonian, hexaquark configuration, and the calculation methods for root-mean-square radii and two-body strong decays. Then in Section III, we present results for mass spectra, internal structure, and decay properties. Finally, the key findings and implications for future studies are summarized in Section IV.

II Theoretical framework for doubly-charmed hexaquark system

II.1 The effective Hamiltonian

Table 1: Parameters of the Hamiltonian determined by fitting the singly-charmed baryon masses. The $M_{\rm the}$, $M_{\rm exp}$, and Error are the theoretical value, the experimental value, and the error between them, respectively.

Parameter	$m_{n}$	$m_{c}$	$\kappa$		$a_{0}$		$\kappa_{0}$
Value	321.0 MeV	1508.0 MeV	$1.5\times 10^{2}$ MeV fm		$2.0\times 10^{-2}$ $\rm(MeV^{-1}fm)^{1/2}$		$2.8\times 10^{2}$ MeV
Parameter	$m_{s}$	$D$	$\alpha$		$\beta$		$\gamma$
Value	642.0 MeV	1033.0 MeV	1.3 $\rm fm^{-1}$		$-4.4\times 10^{-4}$ $\rm(MeVfm)^{-1}$		$-5.6\times 10^{-4}$ $\rm MeV^{-1}$
Baryon	$\Lambda_{c}$	$\Sigma_{c}$	$\Sigma^{*}_{c}$	$\Xi_{c}$	$\Xi^{\prime}_{c}$	$\Xi^{*}_{c}$	$\Omega_{c}$	$\Omega^{*}_{c}$
$M_{\rm the}$ (MeV)	2279.8	2447.2	2534.2	2487.2	2572.4	2648.7	2678.1	2747.2
$M_{\rm exp}$ (MeV)	2286.5	2452.9	2517.5	2467.8	2577.4	2645.9	2695.2	2765.9
Error (MeV)	-6.7	-6.8	15.8	19.3	5.0	3.2	-17.1	-18.7

In the constituent quark model, the nonrelativistic Hamiltonian incorporates three key components to describe quark interactions and calculate the properties of ground-state hadrons: a linear confinement potential (modeling long-range quark confinement), a Coulomb-like potential (accounting for short-range chromoelectric interactions), and a hyperfine interaction potential (arising from spin-dependent forces). The Hamiltonian is thus expressed as:

$$H=\sum_{i=1}^{4}(m_{i}+T_{i})-T_{\rm cm}-\frac{3}{4}\sum_{i<j}^{4}\frac{\lambda^{c}_{i}}{2}.\frac{\lambda^{c}_{j}}{2}(V^{\rm Con}_{ij}+V^{\rm SS}_{ij}).$$

Here, $m_{i}$ represents the mass of the $i$-th constituent quark; $T_{i}=\textbf{p}^{2}_{i}/(2m_{i})$ stands for the kinetic energy of the $i$-th quark; $T_{\rm cm}$ denotes the center-of-mass kinetic energy of the corresponding hexaquark system; and $\lambda^{c}_{i}$ is the $SU(3)$ color operator associated with the $i$-th quark.

The confinement potential $V^{{\rm Con}}_{ij}$ and spin-spin interaction potential $V_{ij}^{{\rm SS}}$ are defined as

$$V^{{\rm Con}}_{ij} = -\frac{\kappa}{r_{ij}}+\frac{r_{ij}}{a^{2}_{0}}-D,$$

$$V^{{\rm SS}}_{ij} = \frac{\kappa^{\prime}}{m_{i}m_{j}}\frac{1}{r_{0ij}r_{ij}}e^{-r^{2}_{ij}/r^{2}_{0ij}}\sigma_{i}\cdot\sigma_{j},$$ (2)

where $r_{ij}=|\textbf{r}_{i}-\textbf{r}_{j}|$ is the interquark distance between the $i$-th and the $j$-th quarks, and $\sigma_{i}$ denotes the $SU(2)$ spin operator of the $i$-th quark. The parameters $r_{0ij}$ and $\kappa^{\prime}$ incorporate explicit mass dependence

$$r_{0ij} = 1/(\alpha+\beta\frac{m_{i}m_{j}}{m_{i}+m_{j}}),$$

$$\kappa^{\prime} = \kappa_{0}(1+\gamma\frac{m_{i}m_{j}}{m_{i}+m_{j}}).$$ (3)

The numerical values of the parameters in Eqs. (II.1)-(II.1) are listed in Table 1. For completeness, the table also includes the calculated mass of singly-charmed baryons within the heavy-quark-diquark configuration, alongside their experimental values for comparative analysis.

II.2 Hexaquark configuration

When two light quarks are treated as a tightly bound diquark, the doubly-charmed hexaquark system reduces to a four-body system in the [diquark-diquark]-[quark-quark] configuration.

The total wave function comprises flavor, spatial, color, and spin components. In the flavor structure, six distinct flavor configurations exist for the doubly-charmed hexaquark system: $[nn][nn]cc$, $[ss][ss]cc$, $[ns][ns]cc$, $[nn][ns]cc$, $[nn][ss]cc$, and $[ns][ss]cc$.

In the color structure, a diquark in the color-antitriplet representation ($(qq)^{\bar{3}_{c}}$) is considered “good” due to its attractive confinement potential, while the color-sextet $(qq)^{6_{c}}$ (“bad” diquark) exhibits repulsive interactions. Only color-singlet configurations are physically admissible. The color decomposition proceeds as:

$$([3\otimes 3]\otimes[3\otimes 3])\otimes(3\otimes 3)$$ (4)

$$= ([\bar{3}\oplus 6]\otimes[\bar{3}\oplus 6])\otimes(\bar{3}\oplus 6)\to(\bar{3}\otimes\bar{3})\otimes(\bar{3}\oplus 6)$$

$$= (3\oplus\bar{6})\otimes(\bar{3}\oplus 6)=(3\otimes\bar{3})\oplus(\bar{6}\otimes 6)\oplus\cdots$$

From Eq. (4), two color-singlet configurations emerge:

$$\Psi_{1} = |[(qq)^{\bar{3}_{c}}(qq)^{\bar{3}_{c}}]^{3_{c}}(cc)^{\bar{3}_{c}}\rangle,$$

$$\Psi_{2} = |[(qq)^{\bar{3}_{c}}(qq)^{\bar{3}_{c}}]^{\bar{6}_{c}}(cc)^{6_{c}}\rangle.$$ (5)

In the notation $|[(q_{1}q_{2})^{\bar{3}_{c}}(q_{3}q_{4})^{\bar{3}_{c}}]^{\rm color_{1}}(cc)^{\rm color_{2}}\rangle$, $\rm color_{1}$ and $\rm color_{2}$ stand for the color representations of diquark pair and charm-quark pair, respectively.

In the spin structure, the system admits 20 spin configurations. The general spin wave function is expressed as $|[(q_{1}q_{2})_{\rm spin_{1}}(q_{3}q_{4})_{\rm spin_{2}}]_{\rm spin_{3}}(cc)_{\rm spin_{4}}\rangle_{\rm spin_{5}}$, where $\rm spin_{1}$ and $\rm spin_{2}$ represent the spins of the light diquarks ($q_{1}q_{2}$) and ($q_{3}q_{4}$), respectively. While, $\rm spin_{3}$ and $\rm spin_{4}$ are the total spin of charm-quark pair ($cc$) and diquark pair $[q_{1}q_{2}][q_{3}q_{4}]$, respectively. Finally, $\rm spin_{5}$ represents the total spin of the doubly-charmed hexaquark state. All possible spin wave functions are tabulated in Table 2.

Based on the symmetrized configurations of flavor, color, and spin spaces, we systematically construct the combined flavor$\otimes$color$\otimes$spin space for the doubly-charmed hexaquark system. For the ground state, the spatial wave function is symmetric under the exchange of any two identical quarks or diquarks. In accordance with the Spin-Statistics Theorem, this necessitates that the flavor$\otimes$color$\otimes$spin component must be antisymmetric (symmetric) under the permutation of identical quarks (diquarks).

This framework parallels earlier theoretical studies of the $\Theta^{+}$ pentaquark LEPS:2003wug , where Jaffe and Wilczek Jaffe:2003sg treated two diquarks as bosonic point-like constituents within a diquark-diquark-antiquark configuration. Their model predicted the $\Theta^{+}$ as a bound state of two identical diquarks and an anti-strange quark, alongside an isospin $I=3/2$ multiplet including $\Xi^{--}_{5}$ ($S=-2$, $J^{P}=1/2^{+}$) near 1750 MeV Jaffe:2003sg . Subsequent works by Liu et al. extended this approach to compute magnetic moments of analogous pentaquark states Li:2003cb ; Liu:2003ab ; Huang:2004tn ; Zhang:2004xt . While experimental confirmation of $\Theta^{+}$ and $\Xi^{--}_{5}$ remains elusive BES:2004kia ; Hicks:2012zz , the symmetry-driven methodology developed for these systems provides a valuable foundation for constructing hexaquark wave functions.

Guided by these principles, we derive total wave functions for distinct flavor configurations that rigorously satisfy the Spin-Statistics Theorem. The resulting configurations are tabulated in Tables 4–6.

A critical distinction arises in the role of interaction terms: while Coulomb and linear confinement potentials do not induce color-spin mixing, the hyperfine interaction term does. The hyperfine interaction strength scales as $1/m_{i}m_{j}$, where $m_{i}$ and $m_{j}$ denote the masses of the charm quark or diquarks ($[nn]$, $[ns]$, $[ss]$). Given that all constituent masses exceed 1 GeV, color-spin mixing is strongly suppressed in the doubly-charmed hexaquark system. This suppression significantly simplifies the spectroscopic analysis compared to light-quark systems.

Table 2: All possible spin wave functions in the [diquark–diquark]-[quark-quark] configuration.

$J=3$	$\chi_{1}=|[(qq)_{1}(qq)_{1}]_{2}[cc]_{1}\rangle_{3}$
$J=2$	$\chi_{2}=|[(qq)_{1}(qq)_{1}]_{2}[cc]_{1}\rangle_{2}$	$\chi_{3}=|[(qq)_{1}(qq)_{1}]_{2}[cc]_{0}\rangle_{2}$	$\chi_{4}=|[(qq)_{1}(qq)_{1}]_{1}[cc]_{1}\rangle_{2}$	$\chi_{5}=|[(qq)_{1}(qq)_{0}]_{1}[cc]_{1}\rangle_{2}$
	$\chi_{6}=|[(qq)_{0}(qq)_{1}]_{1}[cc]_{1}\rangle_{2}$
$J=1$	$\chi_{7}=|[(qq)_{1}(qq)_{1}]_{2}[cc]_{1}\rangle_{1}$	$\chi_{8}=|[(qq)_{1}(qq)_{1}]_{1}[cc]_{1}\rangle_{1}$	$\chi_{9}=|[(qq)_{1}(qq)_{1}]_{0}[cc]_{1}\rangle_{1}$	$\chi_{10}=|[(qq)_{1}(qq)_{1}]_{1}[cc]_{0}\rangle_{1}$
	$\chi_{11}=|[(qq)_{1}(qq)_{0}]_{1}[cc]_{1}\rangle_{1}$	$\chi_{12}=|[(qq)_{0}(qq)_{1}]_{1}[cc]_{1}\rangle_{1}$	$\chi_{13}=|[(qq)_{1}(qq)_{0}]_{1}[cc]_{0}\rangle_{1}$	$\chi_{14}=|[(qq)_{0}(qq)_{1}]_{1}[cc]_{0}\rangle_{1}$
	$\chi_{15}=|[(qq)_{0}(qq)_{0}]_{0}[cc]_{1}\rangle_{1}$
$J=0$	$\chi_{16}=|[(qq)_{1}(qq)_{1}]_{1}[cc]_{1}\rangle_{0}$	$\chi_{17}=|[(qq)_{1}(qq)_{0}]_{1}[cc]_{1}\rangle_{0}$	$\chi_{18}=|[(qq)_{0}(qq)_{1}]_{1}[cc]_{1}\rangle_{0}$	$\chi_{19}=|[(qq)_{1}(qq)_{1}]_{0}[cc]_{0}\rangle_{0}$
	$\chi_{20}=|[(qq)_{0}(qq)_{0}]_{0}[cc]_{0}\rangle_{0}$

II.3 Numerical calculation method

II.3.1 Gaussian expansion method

In the spatial space, the Jacobi coordinates for a four-body system are defined in terms of single-particle coordinates $\textbf{x}_{i}$ $(i=1,2,3,4)$ as follows:

$$\xi_{1}=\sqrt{1/2}(\textbf{x}_{1}-\textbf{x}_{2}),$$

$$\xi_{2}=\sqrt{1/2}(\textbf{x}_{3}-\textbf{x}_{4}),$$

$$\xi_{3}=(\frac{m_{1}\textbf{x}_{1}+m_{2}\textbf{x}_{2}}{m_{1}+m_{2}})-(\frac{m_{3}\textbf{x}_{3}+m_{4}\textbf{x}_{4}}{m_{3}+m_{4}}),$$ (6)

$$\textbf{R}=\frac{m_{1}\textbf{x}_{1}+m_{2}\textbf{x}_{2}+m_{3}\textbf{x}_{3}+m_{4}\textbf{x}_{4}}{m_{1}+m_{2}+m_{3}+m_{4}}.$$

Here, $\xi_{1}$ represents the relative Jacobi coordinate within the diquarks $[q_{1}q_{2}]$ and $[q_{3}q_{4}]$. Meanwhile, $\xi_{2}$ denotes the relative Jacobi coordinate between the charm quarks $c$. While $\xi_{3}$ describes the relative Jacobi coordinate between the centers of mass of the two diquarks $([qq][qq])$ and the two charm quarks $(cc)$. Using above Jacobi coordinates (Eq. (II.3.1)), the spatial wave function with well-defined symmetry for the pairs (12) and (34) can be readily constructed. In the center-of-mass frame of the four-body system ($\textbf{R}=0$) and the number of Jacobi coordinates reduces to three (see Fig. 1).

The kinetic term in the Hamiltonian (Eq. (II.1)) can be simplified for calculations. Denoted as $T_{c}$, it is expressed as:

$$T_{c}=\sum^{4}_{i=1}\frac{\textbf{p}^{2}_{x_{i}}}{2m_{i}}-\frac{\textbf{p}^{2}_{R}}{2M}=\frac{\textbf{p}^{2}_{\xi_{1}}}{2m^{\prime}_{1}}+\frac{\textbf{p}^{2}_{\xi_{2}}}{2m^{\prime}_{2}}+\frac{\textbf{p}^{2}_{\xi_{3}}}{2m^{\prime}_{3}},$$ (7)

where the reduced masses $m^{\prime}_{i}$ are given by $m^{\prime}_{1}=\frac{2\times(m_{1}m_{2})}{m_{1}+m_{2}}$, $m^{\prime}_{2}=\frac{2\times(m_{3}m_{4})}{m_{3}+m_{4}}$, and $m^{\prime}_{3}=\frac{(m_{1}+m_{2})\times(m_{3}+m_{4})}{m_{1}+m_{2}+m_{3}+m_{4}}$.

The Gaussian expansion method (GEM) Hiyama:2003cu ; Hiyama:2012sma ; Brink:1998as is employed to solve the four-body Schrödinger equation. This method has been widely applied to baryons Luo:2023sra , tetraquarks Wu:2021rrc , pentaquarks Yan:2023iie , and few-body molecular states Wu:2021kbu ; Luo:2021ggs ; Luo:2022cun . The spatial wave function is expanded using a set of correlated Gaussian bases constructed from the Jacobi coordinates (Eq. (II.3.1)):

$$\Psi(\xi_{1},\xi_{2},\xi_{3})$$

$$= \sum^{n_{1max}}_{n_{1}=1}\sum^{n_{2max}}_{n_{2}=1}\sum^{n_{3max}}_{n_{3}=1}c_{n_{1}n_{2}n_{3}}\Psi^{n_{1}n_{2}n_{3}}(\xi_{1},\xi_{2},\xi_{3})$$

$$= \sum^{n_{1max}}_{n_{1}=1}\sum^{n_{2max}}_{n_{2}=1}\sum^{n_{3max}}_{n_{3}=1}c_{n_{1}n_{2}n_{3}}\rm Exp[-v_{n_{1}}\xi^{2}_{1}-v_{n_{2}}\xi^{2}_{2}-v_{n_{3}}\xi^{2}_{3}].$$

Here, $c_{n_{1}n_{2}n_{3}}$ are expansion coefficients determined by the Rayleigh-Ritz variational method, and $v_{n_{1}}$, $v_{n_{2}}$, $v_{n_{3}}$ are Gaussian range parameters chosen via a geometric progression:

$$v_{n_{i}}=\frac{1}{r^{2}_{n_{i}}},\quad\quad\quad\quad r_{n_{i}}=r_{min_{i}}a^{n_{i}-1},$$

$$a=(\frac{r_{max_{i}}}{r_{min_{i}}})^{\frac{1}{n_{max_{i}}-1}}\quad\quad(i=1,2,3),$$ (9)

where $n_{max_{i}}$ is the number of Gaussian functions, and $a$ is the ratio coefficient. There are three parameters {$r_{max_{i}}$, $r_{min_{i}}$, $n_{max_{i}}$} to be determined through the variation method. Stable results are achieved with {5 fm, 0.7 fm, 5}.

After the above preparations, the eigenvalues of doubly-charmed hexaquark system are obtained by solving the four-body Schrödinger equation:

$$\hat{H}\Psi_{\rm total}(\xi_{1},\xi_{2},\xi_{3})=E\Psi_{\rm total}(\xi_{1},\xi_{2},\xi_{3}),$$ (10)

where the Hamiltonian $\hat{H}$ (Eq. (II.1)) includes the kinetic term $T_{c}$ (Eq. (7)) and the two-body interaction term between (di)quarks. The total wave function $\Psi_{\rm total}(\xi_{1},\xi_{2},\xi_{3})$ combines the spatial part (Eq. (II.3.1)) and the flavor-color-spin part.

The normalization, kinetic, and potential matrix elements are calculated as

$$N^{nn^{\prime}} = \langle\Psi^{n_{1}n_{2}n_{3}}(\xi_{1},\xi_{2},\xi_{3})\psi_{cs}|\psi^{\prime}_{cs}\Psi^{n^{\prime}_{1}n^{\prime}_{2}n^{\prime}_{3}}(\xi_{1},\xi_{2},\xi_{3})\rangle,$$

$$V^{nn^{\prime}}_{\alpha} = \langle\Psi^{n_{1}n_{2}n_{3}}(\xi_{1},\xi_{2},\xi_{3})\psi_{cs}|V_{\alpha}|\psi^{\prime}_{cs}\Psi^{n^{\prime}_{1}n^{\prime}_{2}n^{\prime}_{3}}(\xi_{1},\xi_{2},\xi_{3})\rangle,$$

$$T_{c}^{nn^{\prime}} = \langle\Psi^{n_{1}n_{2}n_{3}}(\xi_{1},\xi_{2},\xi_{3})\psi_{cs}|T_{c}|\psi^{\prime}_{cs}\Psi^{n^{\prime}_{1}n^{\prime}_{2}n^{\prime}_{3}}(\xi_{1},\xi_{2},\xi_{3})\rangle,$$

where $\psi_{cs}$ is the spin-color wave function, $V$ represents $V^{\rm Con}$ and $V^{\rm SS}$ in Eq. (II.1), and $n$ simply refers to $\{n_{1},n_{2},n_{3}\}$ in Eq. (II.3.1). The $V_{\alpha}$ ($\alpha=$1-6) imply $V(r_{12})$, $V(r_{13})$, $V(r_{14})$, $V(r_{23})$, $V(r_{24})$, and $V(r_{34})$. According to Eq. (II.3.1), the Schrödinger equation (Eq. (10)) is transformed into a generalized matrix eigenvalue problem:

$$[T^{nn^{\prime}}_{c}+\sum^{6}_{\alpha=1}V^{nn^{\prime}}_{\alpha}-EN^{nn^{\prime}}]C_{nn^{\prime}}=0.$$ (12)

Solving this yields the eigenvalue $E$ and the internal mass contributions, as shown in Tables 4-6.

II.3.2 Root-mean-square radii

To further probe the inner structure of the doubly-charmed hexaquark, we calculate the root-mean-square (RMS) radii between all pairs of (di)quarks. This parameter is crucial for distinguishing between compact multiquark states and hadronic molecular states. Specifically, a molecular state typically exhibits minimal spatial overlap between hadronic constituents, whereas a compact multiquark state shows significant overlap Luo:2022cun .

The RMS radius $R_{ij}$ is defined as:

$$R^{2}_{ij}=\int(\textbf{x}_{i}-\textbf{x}_{j})^{2}|\Psi(\xi_{1},\xi_{2},\xi_{3})|^{2}d\xi_{1}d\xi_{2}d\xi_{3}.$$ (13)

The results for the different flavor combinations are listed in Tables 4-6. Specifically: $R_{12}$ and $R_{34}$ describe the average distance between two charm quarks ($cc$) and two diquarks ($[qq][qq]$), respectively. $R_{13}$, $R_{14}$, $R_{23}$, and $R_{24}$ represent the average distance between the charm quark $c$ and the diquark $[qq]$. $R_{12-34}$ stands for the average distance between the mass centers of charm quark pair ($cc$) and the diquark pair ($[qq][qq]$). $R_{13-24}$ and $R_{14-23}$ represent the average distance between the mass centers of the two pairs of charm quark-diquark $c[qq]$. Since the two charm quarks $c$ are identical particles, we have the following relationships: $R_{13-24}=R_{14-23}$, $R_{13}=R_{23}$, and $R_{14}=R_{24}$.

If the average distance between the mass centers of the two $c[qq]$ pairs ($R_{13-24}$, $R_{14-23}$) is comparable to or even smaller than the intra-pair quark distances ($R_{12}$,$R_{13}$,$R_{24}$, etc.), it implies strong spatial overlap between the $c[qq]$ clusters. Such a configuration supports the interpretation of a compact hexaquark state, as opposed to a loosely bound molecular system.

II.3.3 Two-body strong decay

Figure 1: Coordinates defined for the $H_{cc}$ system and its two-body decays into a baryon pair $BC$ via quark rearrangement. Here, the $BC$ final state can form via two quark rearrangement pathways: $B_{1}C_{1}$ ($[c_{1}(q_{1}q_{2})][c_{2}(q_{3}q_{4})]$) and $B_{2}C_{2}$ ($[c_{1}(q_{3}q_{4})][c_{2}(q_{1}q_{2})]$), as illustrated in the figure.

In addition to mass spectra, we employ the quark-interchange model to calculate the two-body strong decay widths of doubly-charmed hexaquarks. When the phase space permits, the dominant two-body strong decay for these states is the rearrangement process: $cc[q_{1}q_{2}][q_{3}q_{4}]\to c[q_{1}q_{2}]+c[q_{3}q_{4}]$, as illustrated in Fig. 1. While decays into a doubly-charmed baryon and a light baryon are kinematically allowed, they are suppressed due to diquark dissociation constraints. Contributions from three-body strong decays and radiative/weak decays are negligible, so we focus exclusively on the above rearrangement channel.

The quark-interchange model Barnes:1991em ; Wong:2001td describes two-body strong decays via quark rearrangement, driven by the (di)quark-(di)quark interaction $V_{ij}$. This approach has successfully described decays of exotic states like the $X(3872)$ Zhou:2019swr , $X(4630)$ Yang:2021sue , $X(2900)$ Wang:2020prk ; Liu:2022hbk , $X(6900)$ liu:2020eha , $Z_{c}$ and $Z_{b}$ states Wang:2018pwi ; Xiao:2019spy , hidden and double charm-strange tetraquark Liu:2024fnh , hidden-charm pentaquark $P_{c}$ states Wang:2019spc , hidden-charm pentaquarks with triple strangeness Wang:2021hql , and all-heavy pentaquark states Liang:2024met .

The decay width is given by

$$\Gamma=\frac{|\vec{P}_{B}|}{(2J_{A}+1)32\pi^{2}M_{A}}\int d\Omega|\mathcal{M}(A\to BC)|^{2},$$ (14)

where $\vec{P}_{B}$ is the final-state three-momentum in the center-of-mass reference frame, and $M_{A}$ is the initial hexaquark mass. The decay amplitude $\mathcal{M}(A\to BC)$ is

$$\mathcal{M}(A\to BC)=-(2\pi)^{3/2}\sqrt{2M_{A}}\sqrt{2E_{B}}\sqrt{2E_{C}}\times T,$$ (15)

with $T$-matrix:

$$T = \sum_{i}\langle\Psi^{B}\Psi^{C}|V_{i}|\Psi^{A}\rangle$$ (16)

$$= \sum_{i}\langle\Psi^{B}\Psi^{C}|V_{i}|\Psi^{A}_{(cc)}\Psi^{A}_{([qq][qq])}\Psi^{A}_{(cc)-([qq][qq])}\rangle.\quad$$

Here, $\Psi^{A}$, $\Psi^{B}$, $\Psi^{C}$ denote the spatial wave functions of the initial hexaquark and final baryons.

The $T$-matrix in momentum space integrates the effective potential:

$$T=\frac{1}{(2\pi)^{3}}\int d\vec{P}_{\alpha}V_{{\rm eff}}(\vec{P}_{\alpha},\vec{P}_{B})\Psi^{A}_{(cc)-([qq][qq])}(\vec{P}_{\alpha}),$$ (17)

where $V_{{\rm eff}}(\vec{P}_{\alpha},\vec{P}_{B})$ combines contributions from diagrams $C_{1}$, $C_{2}$, $T_{1}$, and $T_{2}$ in Fig. 2. It factorizes as

$$V_{{\rm eff}}(\vec{P}_{\alpha},\vec{P}_{B})=I_{{\rm flavor}}I_{{\rm color}}I_{{\rm spin-space}}.$$ (18)

Here, the flavor factor $I_{{\rm flavor}}$ is taken as 1 for all diagrams. The color factor $I_{{\rm color}}$ is:

$$I_{{\rm color}}=\langle\Psi_{B}\Psi_{C}|\frac{\lambda^{c}_{i}}{2}.\frac{\lambda^{c}_{j}}{2}|\Psi^{\bar{3}_{c}(6_{c})}_{A}(cc)\Psi^{3_{c}(\bar{6}_{c})}_{A}([qq][qq])\rangle.\quad$$ (19)

Numerical values for each diagram are listed in Table 3. As for the spin-space factor $I_{{\rm spin-space}}$, it is decoupled for doubly-charmed hexaquark ground state. The spin factor $I_{{\rm spin}}$ is

$$I_{{\rm spin}}=\langle[\chi^{B}_{s_{3}}\chi^{C}_{s_{4}}]_{S^{\prime}}|\hat{\mathcal{O}}_{s}|[\chi(cc)^{A}_{s_{1}}\chi([qq][qq])^{A}_{s_{2}}]_{S}\rangle,$$ (20)

where $s_{1}/s_{2}$ and $s_{3}/s_{4}$ represent the spins of initial and final components, and $S$ and $S^{\prime}$ represent the total spin of initial and final state. $\hat{\mathcal{O}}_{s}$ stands for the spin operator, taking 1 for the Coulomb and linear confinement potential, and $\frac{\sigma_{i}}{2}\cdot\frac{\sigma_{j}}{2}$ for the hyperfine potential. The space factor $I_{{\rm space}}$ is:

$$I_{{\rm space}}=\langle\Psi^{B}\Psi^{C}|\hat{\mathcal{O}}_{q}|\Psi^{A}_{(cc)}\Psi^{A}_{([qq][qq])}\rangle$$ (21)

$$=\int\int d\textbf{k}_{1}d\textbf{k}_{2}\Psi^{B}(k_{B}+K_{B})\Psi^{C}(k_{C}+K_{C})\hat{\mathcal{O}}_{q}(k_{1}-k_{2})$$

$$\quad\Psi^{*A}_{(cc)}(k_{\alpha}+K_{\alpha})\Psi^{*A}_{([qq][qq])}(k_{\beta}+K_{\beta}),$$

where $\hat{\mathcal{O}}_{q}$ represents the spatial operator. Its corresponding specific forms are as follows: $1/q^{2}$, $1/q^{4}$, and $\exp[-q^{2}]$ represent the Coulomb, linear confinement, and hyperfine potential, respectively. The $\vec{k}_{1}$ ($\vec{k}_{2}$) is the initial (final) three-momenta of the scattered constituent. We use the $\vec{P}_{B}$, $\vec{P}_{C}$, $\vec{P}_{\alpha}$, and $\vec{P}_{\beta}$ to present the three-momenta of final singly-charmed baryons $B$, $C$ and initial doubly-charmed hexaquark ($cc$) and ($[qq][qq]$) components, respectively. For simplicity, we deal with the scattering problem in the center-of-mass frame, so that $\vec{P}_{B}=-\vec{P}_{C}$ and $\vec{P}_{\alpha}=-\vec{P}_{\beta}$. By applying the above relationships, for the four quark exchange diagrams in Fig. 2, the relationships of $\vec{k}_{i}$ and $\vec{K}_{i}$ ($i=B,C,\alpha,\beta$) expressed in terms of $\vec{P}_{B}$, $\vec{P}_{\alpha}$, $\vec{k}_{1}$, and $\vec{k}_{2}$ are shown in Table 3. Here, the constituent quark mass-dependent function $f_{i}$ ($i=\alpha,\beta_{1},\beta_{2},B,C$) is:

$$f_{\alpha} = \frac{m_{c}}{m_{c}+m_{c}}=\frac{1}{2},\quad f_{\beta_{1}}=\frac{m_{[q_{1}q_{2}]}}{m_{[q_{1}q_{2}]}+m_{[q_{3}q_{4}]}},$$

$$f_{\beta_{2}} = \frac{m_{[q_{3}q_{4}]}}{m_{[q_{1}q_{2}]}+m_{[q_{3}q_{4}]}},\quad f_{B}=\frac{m_{[q_{1}q_{2}]}}{m_{c}+m_{[q_{1}q_{2}]}},$$

$$f_{C} = \frac{m_{[q_{3}q_{4}]}}{m_{c}+m_{[q_{3}q_{4}]}}.$$ (22)

The specific derivation of integral simplification for Eq. (21) is referred to Refs. Barnes:1991em ; Wong:2001td . Finally, with the calculated $T$-matrix element $T_{fi}$, the two-body decay widths are computed using Eq. (14) and presented in Tables 4-6.

Figure 2: The quark-interchange diagrams for $H_{cc}$ decaying into two baryons at the quark level. The curve line denotes the (di)quark–(di)quark interactions.

Table 3: The color matrix element $I_{\rm color}=\langle\frac{\lambda^{c}_{i}}{2}.\frac{\lambda^{c}_{j}}{2}\rangle$ and the momentum substitutions in $I_{\rm space}$ for different scattering diagrams.

Diagram   $\langle\frac{\lambda^{c}_{i}}{2}.\frac{\lambda^{c}_{j}}{2}\rangle$   momentum substitutions 11    $(cc)^{\bar{3}_{c}}([qq][qq])^{3_{c}}$ $(cc)^{6_{c}}([qq][qq])^{\bar{6}_{c}}$   $\vec{k}_{\alpha}$ $\vec{K}_{\alpha}$   $\vec{k}_{\beta}$ $\vec{K}_{\beta}$   $\vec{k}_{B}$ $\vec{K}_{B}$   $\vec{k}_{C}$ $\vec{K}_{C}$   $C_{1}$   $\frac{2\sqrt{3}}{9}$ $\frac{\sqrt{6}}{9}$   $\vec{k}_{1}$ $-f_{\alpha}\vec{P}_{\alpha}$   $\vec{k}_{2}$ $-f_{\beta 2}\vec{P}_{\alpha}-\vec{P}_{B}$   $\vec{k}_{2}$ $-f_{B}\vec{P}_{B}$   $\vec{k}_{2}$ $-f_{C}\vec{P}_{B}-\vec{P}_{\alpha}$ $C_{2}$   $\frac{2\sqrt{3}}{9}$ $\frac{\sqrt{6}}{9}$   $\vec{k}_{1}$ $-f_{\alpha}\vec{P}_{\alpha}$   $\vec{k}_{2}$ $-f_{\beta 1}\vec{P}_{\alpha}+\vec{P}_{B}$   $\vec{k}_{1}$ $f_{B}\vec{P}_{B}-\vec{P}_{\alpha}$   $\vec{k}_{1}$ $f_{\beta 1}\vec{P}_{\alpha}$ $T_{1}$   $-\frac{2\sqrt{3}}{9}$ $-\frac{\sqrt{6}}{9}$   $\vec{k}_{1}$ $-f_{\alpha}\vec{P}_{\alpha}$   $\vec{k}_{2}$ $-f_{\beta 2}\vec{P}_{\alpha}-\vec{P}_{B}$   $\vec{k}_{2}$ $-f_{\alpha}\vec{P}_{\alpha}$   $\vec{k}_{1}$ $-f_{C}\vec{P}_{B}-\vec{P}_{\alpha}$ $T_{2}$   $-\frac{2\sqrt{3}}{9}$ $-\frac{\sqrt{6}}{9}$   $\vec{k}_{1}$ $-f_{\alpha}\vec{P}_{\alpha}$   $\vec{k}_{2}$ $-f_{\beta 1}\vec{P}_{\alpha}+\vec{P}_{B}$   $\vec{k}_{1}$ $f_{B}\vec{P}_{B}-\vec{P}_{\alpha}$   $\vec{k}_{2}$ $f_{\beta 1}\vec{P}_{\alpha}$

III Numerical results and discussions

By solving the Schrödinger equation Eq. (12) with the Gaussian expansion method, we systematically calculated the mass spectra of the doubly-charmed hexaquark system. The mass range of the doubly-charmed hexaquark system approximately spans 5000–6000 MeV. The mass spectra and relevant configurations are listed in Tables 4-6. The internal contributions from each part of the Hamiltonian: kinetic energy $\langle T\rangle$, confinement potential $\langle V^{\rm Con}\rangle$, and hyperfine interaction potential $\langle V^{\rm SS}\rangle$, are also presented in the same tables. The results show that the kinetic energy $\langle T\rangle$ and confinement potential $\langle V^{\rm Con}\rangle$ are of the same order of magnitude. Additionally, it is found that the contribution of hyperfine interaction potential $\langle V^{\rm SS}\rangle$—proportional to $(1/m_{i}m_{j})$—significantly suppresses its mass contribution. Owing to its smaller mass contribution, the degree of configuration mixing remains relatively low. As a result, configuration mixing does not induce obvious mass shifts in physical states. Thus, each original configuration maintains high purity in its corresponding physical state. Meanwhile, the mass gaps between different color-spin configurations of the same color configuration ($|\bar{3}3\rangle_{c}$ or $|6\bar{6}\rangle_{c}$) are relatively small, leading to the existence of some partner states in the doubly-charmed hexaquark system.

Besides the mass spectra and internal mass contributions, we provide the corresponding RMS radii according to Eq. (13). Based on the results from relevant tables, most of RMS radii are in the range from 1.2 to 1.6 fm, which are roughly the same order of magnitude. If it is a molecular-state configuration, some distances between charm quark ($c$) and diquark ($[qq]$) are much greater than those between the two charm quarks ($cc$) and between the two diquarks ($[qq]-[qq]$). Meanwhile, the RMS radius of the molecular configuration can reach several femtometers. Therefore, our calculation results are consistent with the expectations of the compact hexaquark configuration. Finally, based on Eq. (14), we presented the partial widths of each state decaying into different two-body final states, as well as its total width. It should be noted that the total decay width here ignores the suppressed three-body strong decays, the two-body strong decays with the final states: doubly-charmed baryon+light-flavor baryon ($ccq+qqq$), as well as the radiative decays and weak decays. Therefore, the actual total widths of these states will be slightly larger than our calculated values.

For clarity, according to Tables 4-6, the relative mass positions of each state, the total decay widths, and the corresponding rearrangement decay channels are plotted in Figs. 3-5. For convenience, we also label all possible spin (isospin) quantum numbers of the rearrangement decay channels with subscripts (superscripts). According to the above-mentioned figures, there is no stable state in the doubly-charmed hexaquark system, and they are all unstable states which can decay into two singly-charmed baryons through two-body strong interaction. The reason is that the pairwise attractive interaction provided by $\langle V^{\rm Con}\rangle$ is far smaller than that of the two singly-charmed baryons in the decay final states. Therefore, their masses are higher than the threshold of the decay final states. For simplicity, we use the notation ${\rm H_{content}}({\rm Mass},I,J^{P})$ to label a particular doubly-charmed hexaquark state.

III.1 The $cc[nn][nn]$ and $cc[ss][ss]$ subsystems

Table 4: The numerical results of the mass spectrum, the mass contributions of each Hamiltonian part (in MeV), the root-mean-square radii (in fm), and the partial decay widths and total decay widths of the fall-apart decay processes (in MeV) for the $cc[nn][nn]$ and $cc[ss][ss]$ hexaquark states.

$cc[nn][nn]$ Internal contribution RMS Radius Fall-apart decay properties 19  $I[J^{P}]$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Sigma^{*}_{c}\Sigma^{*}_{c}$ $\Sigma^{*}_{c}\Sigma_{c}$ $\Sigma_{c}\Sigma_{c}$ $\Lambda_{c}\Sigma^{*}_{c}$ $\Lambda_{c}\Sigma_{c}$ $\Lambda_{c}\Lambda_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $2(0)[2^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=2,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=2}_{s=2}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=2,\bar{6}_{c}}_{s=2}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=2}_{s=2}\end{pmatrix}$ $\begin{pmatrix}5598\\ 5504\end{pmatrix}$ 1377.9 -1571.9 3.9 1.29 1.26 1.36 1.36 1.44 1.28 29.0 0.1 29.0 1369.7 -1656.9 3.6 1.49 1.47 1.31 1.31 1.12 1.48 19.3 19.4 38.7 $2(0)[1^{+}]$ $|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=2,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=2}_{s=1}$ 5577 1389.3 -1583.2 -17.1 1.29 1.25 1.35 1.35 1.43 1.27 6.0 0.7 28.1 34.8 $2(0)[0^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=2,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=2}_{s=0}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=2,\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=2}_{s=0}\\ |[(nn)^{I=0,\bar{3}_{c}}_{s=0}(nn)^{I=0,\bar{3}_{c}}_{s=0}]^{I=0,\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=0}_{s=0}\end{pmatrix}$ $\begin{pmatrix}5582\\ 5512\\ 5043\end{pmatrix}$ 1386.4 -1580.4 -11.8 1.29 1.26 1.35 1.35 1.43 1.27 23.3 21.0 44.3 1365.3 -1652.5 11.3 1.49 1.47 1.31 1.31 1.13 1.48 22.9 4.1 27.0 1-19[1pt/1pt] $0[0^{+}]$ 1389.5 -1471.5 5.8 1.63 1.52 1.40 1.40 1.20 1.58 14.7 14.7 $1[3^{+}]$ $|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=3}$ 5590 1372.3 1632.8 -4.1 1.30 1.26 1.36 1.36 1.43 1.27 36.7 36.7 $1[2^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=2}\\ |[(nn)^{I=0,\bar{3}_{c}}_{s=0}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=2}\end{pmatrix}$ $\begin{pmatrix}5600\\ 5375\end{pmatrix}$ 1376.6 -1570.6 6.4 1.30 1.26 1.36 1.36 1.43 1.28 4.8 9.5 15.3 1383.5 -1475.0 10.4 1.37 1.27 1.45 1.36 1.27 1.48 70.1 70.1 $1[1^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=0}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=1}\\ |[(nn)^{I=0,\bar{3}_{c}}_{s=0}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=1}\\ |[(nn)^{I=0,\bar{3}_{c}}_{s=0}(nn)^{I=1,\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=1}\end{pmatrix}$ $\begin{pmatrix}5590\\ 5586\\ 5510\\ 5365\\ 5283\end{pmatrix}$ 1382.3 -1576.3 -4.1 1.30 1.26 1.36 1.36 1.43 1.27 5.7 0.6 26.6 32.9 1384.3 -1578.2 -7.8 1.29 1.26 1.36 1.36 1.44 1.27 19.6 0.4 0.8 20.8 1366.8 -1654.0 8.7 1.49 1.47 1.31 1.31 1.13 1.48 19.4 2.2 2.8 24.4 1389.0 -1480.5 -0.1 1.37 1.26 1.45 1.36 1.48 1.31 7.8 3.5 11.3 1376.3 -1558.8 9.6 1.56 1.49 1.42 1.30 1.16 1.52 15.4 2.7 18.1 $cc[ss][ss]$ Internal contribution RMS Radius Fall-apart decay properties 19  $I[J^{P}]$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Omega^{*}_{c}\Omega^{*}_{c}$ $\Omega^{*}_{c}\Omega_{c}$ $\Omega_{c}\Omega_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $0[2^{+}]$ $\begin{pmatrix}|[(ss)^{I=0,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=2}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=0,\bar{6}_{c}}_{s=2}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=0}_{s=2}\end{pmatrix}$ $\begin{pmatrix}6038\\ 5945\end{pmatrix}$ 1370.4 -1692.1 5.4 1.20 1.25 1.30 1.30 1.38 1.22 40.0 1.0 41.0 1361.9 -1775.8 4.8 1.40 1.43 1.26 1.26 1.08 1.42 15.8 10.4 26.2 $0[1^{+}]$ $|[(ss)^{I=0,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=1}$ 6021 1379.5 -1701.2 -11.3 1.19 1.24 1.30 1.30 1.37 1.22 2.1 9.1 0.7 11.9 $0[0^{+}]$ $\begin{pmatrix}|[(ss)^{I=0,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=0}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=0,\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=0}_{s=0}\end{pmatrix}$ $\begin{pmatrix}6025\\ 5950\end{pmatrix}$ 1377.3 -1698.9 -7.1 1.19 1.25 1.30 1.30 1.37 1.22 27.7 35.5 63.2 1359.3 -1773.2 9.5 1.40 1.43 1.26 1.26 1.08 1.42 31.0 8.7 39.7

Figure 3: Relative positions for the $cc[nn][nn]$ (a) and $cc[ss][ss]$ (b) hexaquark states labeled with horizontal solid lines, the labels, e.g. $5582~{}(\Gamma=43)$ represents the mass and total decay width of the corresponding state (units: MeV) . In the $cc[nn][nn]$ subsystem, the black, red, and purple horizontal lines represent the hexaquark states with $I=2,0$, $I=1$, and $I=0$, respectively. The dotted lines denote various S-wave baryon-baryon thresholds, and the superscripts (subscript) of the labels, e.g. $(\Sigma^{*}_{c}\Sigma^{*}_{c})^{3,2,1,0}_{2,1,0}$, represent the possible total angular momenta (isospin) of the channels. The solid dots of different colors where the vertical dashed lines with arrows intersect the horizontal solid lines represent the allowed rearranged S-wave decay processes. If a vertical dashed line with an arrow intersects a horizontal solid line without a solid dot, it means that corresponding decay process is forbidden for relevant state.

In the following, we first discuss the $cc[nn][nn]$ and $cc[ss][ss]$ subsystems. For the $cc[nn][nn]$ subsystem, the isospin of two $[nn]$ diquarks can couple to $I=2,1,$ and $0$. According to the Spin-Statistics Theorem, the $I=2$ states and the $I=0$ states have the same symmetry constraints for total wave functions, leading to their degeneracy in mass spectra, RMS radii, and decay properties. A similar situation also occurs in doubly-heavy pentaquark system. For the $ccnn\bar{n}$ and $bbnn\bar{n}$ subsystems with $I_{nn}=1$, the same spectra are obtained for the case of the total isospin $I=1/2$ and $3/2$ Zhou:2018bkn . This arises because both the confinement potential and hyperfine interaction in the Hamiltonian are independent of isospin.

Of course, there is also a special configuration: $|[(nn)^{I=0}_{s=0}(nn)^{I=0}_{s=0}]^{I=0,\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=0}_{s=0}$. In this configuration, the two $[nn]$ diquarks only couple to total isospin $I=0$. This feature leads to an additional allowed state: ${\rm H}_{c^{2}[nn]^{2}}(5043,0,0^{+})$ in the $I(J^{P})=0(0^{+})$ state, which is absent in the $I(J^{P})=0(2^{+})$ sector. Notably, due to the significant suppression of color-spin mixing, the mass gaps between $I(J^{P})=0(0^{+})$ and $0(2^{+})$ states are negligible ($\Delta m<1$ MeV ), maintaining approximate degeneracy. Apart from ${\rm H}_{c^{2}[nn]^{2}}(5043,0,0^{+})$ state, the remaining $0(0^{+})$ and $2(0^{+})$ states still have the same mass, RMS radii, and decay behavior, and are represented by the same data set in Table 4.

For the ${\rm H}_{c^{2}[nn]^{2}}(5043,0,0^{+})$ state, because of the low mass and narrow width, it may be the most ideal in the doubly-charmed hexaquark system. In this scheme, it is potential to observe this state. Its unique properties are as follows: first, both $[nn]$ diquarks in the state are scalar diquarks with $I=0$, $S=0$. Compared with other diquark configurations, the structures of this system have an antisymmetric color part and an antisymmetric spin part, which result in the strongest internal confinement potential and chromomagnetic interaction, and lead to higher stability. Consequently, this configuration achieves the lowest mass (5043 MeV) within the doubly-charmed hexaquark spectrum. Second, this state mainly decays into the $\Lambda_{c}\Lambda_{c}$ final state, while other decay channels are significantly suppressed. Since the $\Lambda_{c}$ is well established, we strongly suggest searching for hexaquark state in the $\Lambda_{c}\Lambda_{c}$ channel. Finally, its total width $\Gamma_{\rm total}=14.3$ MeV. Though it is wider than the sub-MeV-scale width of the experimentally discovered $T_{cc}^{+}(3875)$, it still remains characteristic of a narrow hadronic resonance. This is because $T_{cc}^{+}(3875)$ lies slightly below the $D^{*+}D^{0}$ threshold, making the three-body decay $D^{0}D^{0}\pi$ its dominant channel and resulting in an extremely narrow width. In contrast, ${\rm H}_{c^{2}[nn]^{2}}(5043,0,0^{+})$ state is above the $\Lambda_{c}\Lambda_{c}$ threshold, with a larger decay phase space.

The discovery of $T_{cc}^{+}(3875)$ again demonstrates the experimental detectability of the doubly-charmed multi-quark system. Therefore, we suggest that experimental collaborations like LHCb, CMS, ATLAS, etc. produce the $\Lambda_{c}\Lambda_{c}$ final state via high-luminosity $pp$ collisions and analyze the invariant mass spectrum within the $5000$–$5100$ MeV range to search for this narrow peak structure.

For other two $I[J^{P}]=2(0)[0^{+}]$ states: ${\rm H}_{c^{2}[nn]^{2}}(5582,\\ 2(0),0^{+})$ and ${\rm H}_{c^{2}[nn]^{2}}(5512,2(0),0^{+})$, due to large phase space, they both decay into the $\Sigma^{*}_{c}\Sigma^{*}_{c}$ and $\Sigma_{c}\Sigma_{c}$ final states. Their mass gap reaches 70 MeV, because of different color configurations. Further, we find that compared with the $\bar{3}\otimes 3$ color configuration, the $6\otimes\bar{6}$ color configuration exhibits stronger attractive interactions, resulting in lower masses for these states. Meanwhile, their total decay widths are 44 MeV and 27 MeV, respectively. The corresponding relative partial decay width ratios are as follows:

$$\frac{\Gamma[{\rm H}_{c^{2}[nn]^{2}}(5582,2(0),0^{+})\to\Sigma_{c}\Sigma_{c}]}{\Gamma[{\rm H}_{c^{2}[nn]^{2}}(5582,2(0),0^{+})\to\Sigma^{*}_{c}\Sigma^{*}_{c}]}=1,$$ (23)

and

$$\frac{\Gamma[{\rm H}_{c^{2}[nn]^{2}}(5512,2(0),0^{+})\to\Sigma_{c}\Sigma_{c}]}{\Gamma[{\rm H}_{c^{2}[nn]^{2}}(5512,2(0),0^{+})\to\Sigma^{*}_{c}\Sigma^{*}_{c}]}=\frac{2}{11},$$ (24)

respectively. For other $I=2(0)$ $cc[nn][nn]$ states, one can perform similar discussions on the decay behaviors according to Table 4 and Fig. 3.

For the $I=1$ $cc[nn][nn]$ states, there is no ground $I[J^{P}]=1[0^{+}]$ state due to symmetry constraints. Meanwhile, the total decay widths of most of states are larger than 20 MeV, classifying them as relatively broad states. Among these states, the ${\rm H}_{c^{2}[nn]^{2}}(5365,1,1^{+})$ is the narrowest state. Although its decay phase space is larger than that of ${\rm H}_{c^{2}[nn]^{2}}(5283,1,1^{+})$, its total width is only about 10 MeV. For this state, we have

$$\frac{\Gamma[{\rm H}_{c^{2}[nn]^{2}}(5365,1,1^{+})\to\Sigma_{c}\Lambda_{c}]}{\Gamma[{\rm H}_{c^{2}[nn]^{2}}(5365,1,1^{+})\to\Sigma^{*}_{c}\Lambda_{c}]}=\frac{1}{2}.$$ (25)

This indicates that the $\Sigma^{*}_{c}\Lambda_{c}$ decay channel is dominant. The narrow width of ${\rm H}_{c^{2}[nn]^{2}}(5365,1,1^{+})$ makes its peak shape significantly distinguishable from the background. Therefore, we suggest that experiments scan the invariant mass spectrum of $\Sigma^{*}_{c}\Lambda_{c}$ in the mass range of $5300$-$5400$ MeV, with particular attention to the narrow peak structure near 5360 MeV. If the peak position and width observed in the experiment are consistent with theoretical predictions, and the branching ratio conforms to the 1:2 ratio, it can be confirmed as the signal of ${\rm H}_{c^{2}[nn]^{2}}(5365,1,1^{+})$.

Moreover, for the ${\rm H}_{c^{2}[nn]^{2}}(5590,1,1^{+})$ and ${\rm H}_{c^{2}[nn]^{2}}\\ (5586,1,1^{+})$ states, they are degenerate states, which have same quantum numbers and similar masses ($\Delta M=4$ MeV). Although they have similar decay phase space, the total decay width of ${\rm H}_{c^{2}[nn]^{2}}(5590,1,1^{+})$ is 33 MeV, approximately 1.5 times that of ${\rm H}_{c^{2}[nn]^{2}}(5586,1,1^{+})$. And their relative partial decay width ratios are

$$\Gamma_{\Sigma^{*}_{c}\Sigma^{*}_{c}}:\Gamma_{\Sigma_{c}\Sigma^{*}_{c}}:\Gamma_{\Sigma_{c}\Sigma_{c}}=1:0.1:5,$$ (26)

and

$$\Gamma_{\Sigma^{*}_{c}\Sigma^{*}_{c}}:\Gamma_{\Sigma_{c}\Sigma^{*}_{c}}:\Gamma_{\Sigma_{c}\Sigma_{c}}=25:0.5:1,$$ (27)

respectively. Evidently, $\Sigma_{c}\Sigma_{c}$ and $\Sigma^{*}_{c}\Sigma^{*}_{c}$ are their dominant decay channels, respectively. Although theoretically we can distinguish them by their total decay widths and the branching ratios, current experimental detectors still face great difficulties in distinguishing degenerate states with a mass difference ($\Delta M$) of 4 MeV.

For the $cc[ss][ss]$ subsystem, since it has exactly the same symmetry constraints as the $cc[nn][nn]$ subsystem with $I=2$, the number of allowed states is also identical. Due to symmetry constraints, there is no ground $I[J^{P}]=0[3^{+}]$ state. Among them, the ${\rm H}_{c^{2}[ss]^{2}}(6021,0,1^{+})$ is the narrowest state whose total width is around 12 MeV, even though it has more allowed decay channels: $\Omega_{c}\Omega_{c}$, $\Omega^{*}_{c}\Omega_{c}$, and $\Omega^{*}_{c}\Omega^{*}_{c}$. Here, we obtain the following relative ratio of decay widths:

$$\Gamma_{\Omega^{*}_{c}\Omega^{*}_{c}}:\Gamma_{\Omega^{*}_{c}\Omega_{c}}:\Gamma_{\Omega_{c}\Omega_{c}}=3:13:1.$$ (28)

Our results show that the $\Omega^{*}_{c}\Omega_{c}$ channel is its dominant decay channel. Moreover, all the $cc[ss][ss]$ states can decay into $\Omega^{*}_{c}\Omega^{*}_{c}$ final states, and this decay channel is crucial for identifying $cc[ss][ss]$ states. Although theoretical predictions indicate clear signals in the $cc[ss][ss]$ subsystem, experimental discovery of this subsystem still faces significant challenges. Firstly, the production probability of the $[ss]$ quark pair in proton-proton collisions is much lower than that of the $[nn]$ pair. Additionally, the experimental reconstruction of its decay final states $\Omega^{*}_{c}$ and $\Omega_{c}$ is relatively complex.

III.2 The $cc[nn][ns]$ and $cc[ss][sn]$ subsystems

Table 5: The numerical results of the mass spectrum, the mass contributions of each Hamiltonian part (in MeV), the root-mean-square radii (in fm), and the partial decay widths and total decay widths of the fall-apart decay processes (in MeV) for the $cc[nn][ns]$ and $cc[ss][sn]$ hexaquark states.

$cc[nn][ns]$ Internal contribution RMS Radius Fall-apart decay properties 19  $I[J^{P}]$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Sigma^{*}_{c}\Xi^{*}_{c}$ $\Sigma^{*}_{c}\Xi^{\prime}_{c}$ $\Sigma_{c}\Xi^{*}_{c}$ $\Sigma_{c}\Xi^{\prime}_{c}$ $\Sigma^{*}_{c}\Xi_{c}$ $\Sigma_{c}\Xi_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $\frac{3}{2}(\frac{1}{2})[3^{+}]$ $|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=3}$ 5734 1366.6 -1597.1 20.6 1.27 1.26 1.36 1.33 1.42 1.26 42.1 42.1 $\frac{3}{2}(\frac{1}{2})[2^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=2}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=2}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=2}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{3}{2}}_{s=2}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=2}\end{pmatrix}$ $\begin{pmatrix}5719\\ 5718\\ 5621\\ 5586\end{pmatrix}$ 1374.5 -1605.1 6.0 1.27 1.25 1.36 1.33 1.42 1.26 4.7 11.8 38.9 55.4 1375.3 -1605.9 4.4 1.27 1.26 1.36 1.33 1.42 1.26 35.2 0.2 18.8 54.2 1367.2 -1691.1 3.9 1.46 1.46 1.32 1.27 1.11 1.46 14.4 7.7 10.4 32.5 1374.9 -1592.5 10.4 1.30 1.26 1.36 1.37 1.44 1.28 22.1 22.1 $\frac{3}{2}(\frac{1}{2})[1^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=0}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{3}{2}}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{3}{2}}_{s=1}\end{pmatrix}$ $\begin{pmatrix}5710\\ 5708\\ 5707\\ 5629\\ 5575\\ 5488\end{pmatrix}$ 1379.9 -1610.5 -3.9 1.27 1.25 1.36 1.32 1.41 1.26 6.9 0.9 15.6 32.1 55.5 1380.7 -1611.2 -5.4 1.26 1.25 1.36 1.32 1.42 1.26 2.7 4.7 59.7 0.8 67.9 1381.1 -1611.6 -6.2 1.26 1.25 1.36 1.32 1.42 1.26 23.6 0.7 38.1 0.9 63.3 1364.7 -1688.5 8.5 1.47 1.46 1.32 1.28 1.11 1.46 23.4 2.5 3.4 3.5 32.8 1380.5 -1568.1 0.0 1.30 1.26 1.36 1.36 1.44 1.28 26.0 41.6 67.6 1368.7 -1649.4 6.1 1.49 1.47 1.31 1.32 1.13 1.48 21.0 7.7 28.7 $\frac{3}{2}(\frac{1}{2})[0^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=0}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{3}{2}}_{s=0}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{3}{2}}_{s=0}\end{pmatrix}$ $\begin{pmatrix}5703\\ 5631\\ 5570\end{pmatrix}$ 1383.4 -1613.9 -10.3 1.26 1.25 1.35 1.32 1.41 1.26 28.2 20.0 48.2 1363.4 -1687.2 10.7 1.47 1.46 1.32 1.28 1.11 1.46 27.8 5.2 33.0 1383.4 -1570.9 -5.3 1.30 1.26 1.36 1.36 1.43 1.28 97.5 97.5 $I[J^{P}]$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Lambda_{c}\Xi^{*}_{c}$ $\Lambda_{c}\Xi^{\prime}_{c}$ $\Lambda_{c}\Xi_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $\frac{1}{2}(2^{+})$ $|[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{1}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=2}$ 5496 1380.9 -1506.8 9.7 1.35 1.26 1.45 1.33 1.46 1.30 53.9 53.9 $\frac{1}{2}(1^{+})$ $\begin{pmatrix}|[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{1}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=1}\\ |[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{1}{2},\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{1}{2}}_{s=1}\\ |[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{1}{2},3_{c}}_{s=0}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=1}\end{pmatrix}$ $\begin{pmatrix}5487\\ 5401\\ 5351\end{pmatrix}$ 1385.7 -1511.6 0.6 1.35 1.26 1.45 1.33 1.46 1.30 10.3 6.1 16.4 1376.5 -1593.9 6.1 1.54 1.48 1.43 1.26 1.14 1.50 19.0 3.5 22.5 1386.9 -1472.2 5.2 1.38 1.27 1.45 1.37 1.48 1.32 26.0 26.0 $\frac{1}{2}(0^{+})$ $\begin{pmatrix}|[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{1}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=0}\\ |[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{1}{2}}_{s=0}\end{pmatrix}$ $\begin{pmatrix}5482\\ 5261\end{pmatrix}$ 1388.1 -1514.0 -4.0 1.35 1.26 1.45 1.33 1.46 1.30 4.7 4.7 1378.7 -1555.0 6.0 1.56 1.49 1.42 1.30 1.16 1.53 18.5 18.5 $cc[ss][sn]$ Internal contribution RMS Radius Fall-apart decay properties 19  $I(J^{P})$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Omega^{*}_{c}\Xi^{*}_{c}$ $\Omega_{c}\Xi^{*}_{c}$ $\Omega^{*}_{c}\Xi^{\prime}_{c}$ $\Omega_{c}\Xi_{c}$ $\Omega^{*}_{c}\Xi_{c}$ $\Omega_{c}\Xi_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $\frac{1}{2}(3^{+})$ $|[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=3}$ 5951 1364.5 -1659.1 18.3 1.22 1.25 1.31 1.33 1.39 1.24 51.1 51.1 $\frac{1}{2}(2^{+})$ $\begin{pmatrix}|[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=2}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=2}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=2}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{1}{2}}_{s=2}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=2}\end{pmatrix}$ $\begin{pmatrix}5938\\ 5937\\ 5844\\ 5805\end{pmatrix}$ 1371.6 -1666.2 5.3 1.22 1.25 1.30 1.33 1.39 1.23 3.8 43.6 15.2 62.6 1371.7 -1666.3 5.1 1.22 1.25 1.30 1.33 1.39 1.23 39.8 20.7 0.8 61.3 1363.3 -1750.7 4.5 1.42 1.44 1.25 1.29 1.09 1.43 15.8 13.0 9.7 38.5 1371.7 -1621.6 9.4 1.25 1.25 1.31 1.37 1.41 1.25 27.7 27.7 $\frac{1}{2}(1^{+})$ $\begin{pmatrix}|[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=1}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=1}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=0}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=1}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{1}{2}}_{s=1}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=1}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{1}{2}}_{s=1}\end{pmatrix}$ $\begin{pmatrix}5930\\ 5929\\ 5928\\ 5848\\ 5797\\ 5709\end{pmatrix}$ 1376.4 -1671.0 -3.4 1.22 1.25 1.30 1.32 1.38 1.23 6.6 15.4 1.3 39.6 62.9 1376.5 -1671.0 -3.6 1.22 1.25 1.30 1.32 1.39 1.23 2.2 62.7 7.5 0.7 73.1 1376.5 -1671.0 -3.7 1.21 1.25 1.30 1.32 1.39 1.23 25.2 40.4 2.1 1.4 69.1 1361.3 -1748.7 8.1 1.42 1.44 1.25 1.29 1.09 1.43 25.9 4.3 3.2 5.3 38.7 1376.2 -1626.0 1.1 1.25 1.25 1.30 1.36 1.41 1.25 26.1 45.5 71.6 1365.3 -1708.2 6.2 1.45 1.45 1.25 1.33 1.10 1.45 23.6 9.3 32.9 $\frac{1}{2}(0^{+})$ $\begin{pmatrix}|[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=0}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=\frac{3}{2},\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=\frac{1}{2}}_{s=0}\\ |[(ss)^{I=0,\bar{3}_{c}}_{s=1}(sn)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=\frac{3}{2},3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=\frac{1}{2}}_{s=0}\end{pmatrix}$ $\begin{pmatrix}5925\\ 5850\\ 5793\end{pmatrix}$ 1378.8 -1673.4 -8.0 1.21 1.25 1.30 1.32 1.39 1.23 27.7 32.7 60.4 1360.4 -1747.7 9.8 1.42 1.44 1.26 1.29 1.09 1.43 30.8 7.8 38.6 1378.4 -1628.3 -3.1 1.25 1.25 1.30 1.36 1.41 1.25 99.7 99.7

Figure 4: Relative positions for the $cc[nn][ns]$ (a) and $cc[ss][sn]$ (b) hexaquark states labeled with horizontal solid lines, the labels, e.g. $5703~{}(\Gamma=48)$ represents the mass and total decay width of the corresponding state (units: MeV). In the $cc[nn][ns]$ subsystem, the black and red horizontal lines represent the hexaquark states with $I=3/2$ and $I=1/2$, respectively. The notations are same as those of Fig. 3.

Next, we discuss the $cc[nn][ns]$ and $cc[ss][sn]$ subsystems. For the $cc[nn][ns]$ subsystem, $|[(nn)^{I=0}(ns)^{I=\frac{1}{2}}]\\ (cc)\rangle$ configuration only couples to $I=1/2$ states. On the contrary, $|[(nn)^{I=1}(ns)^{I=\frac{1}{2}}](cc)\rangle$ configuration can couple to $I=3/2$ and $I=1/2$ states. Although color-spin mixing occurs in two above configuration for the $I=1/2$ states, as can be seen from Table 5, the $\langle V^{\rm SS}\rangle$ is significantly suppressed, resulting in the mixing negligible. Thus, the $I=3/2$ and $I=1/2$ states from the $|[(nn)^{I=1}(ns)^{I=\frac{1}{2}}](cc)\rangle$ configuration would have the same mass spectra, RMS radii, and decay behaviors. For these states, they are relatively broad states, and have many different decay channels. Therefore, there are obstacles to experimentally discovering these states. Their resonance peaks are prone to be misjudged as continuous background fluctuations. Thus, we do not recommend that relevant experiments attempt to reconstruct resonance peaks from these decay final states.

Here, we take the ${\rm H}_{c^{2}[nn][ns]}(5710,3/2(1/2),2^{+})$ as an example for discussion. Similar situations also apply to other states. Its total decay width reaches 56 MeV and relative partial decay width ratio is:

$$\Gamma_{\Sigma^{*}_{c}\Xi^{*}_{c}}:\Gamma_{\Sigma^{*}_{c}\Xi^{\prime}_{c}}:\Gamma_{\Sigma_{c}\Xi^{*}_{c}}:\Gamma_{\Sigma_{c}\Xi^{\prime}_{c}}=7.7:1:17:36,$$ (29)

i.e. $\Sigma_{c}\Xi^{*}_{c}$ and $\Sigma_{c}\Xi^{\prime}_{c}$ are its dominant decay channels.

According to Fig. 4, for the $cc[nn][ns]$ states with the configuration $|[(nn)^{I=0}(ns)^{I=\frac{1}{2}}](cc)\rangle$, their masses are generally lower than those of the states with the configuration $|[(nn)^{I=1}(ns)^{I=\frac{1}{2}}](cc)\rangle$. This is because the mass of the scalar diquark $[nn]$ is lower than that of the vector diquark $[nn]$ with $I=1$ and $S=1$. Meanwhile, the total decay widths of $|[(nn)^{I=0}(ns)^{I=\frac{1}{2}}](cc)\rangle$ are generally smaller than those of $|[(nn)^{I=1}(ns)^{I=\frac{1}{2}}](cc)\rangle$, and the same situation exists in the $cc[nn][nn]$ and $cc[nn][ss]$ subsystems. This indicates that compared with other states, those containing the scalar diquark configuration have lower masses, stronger internal interactions, narrower total widths, relatively longer lifetimes, and are more likely to be detected experimentally in the doubly-charmed hexaquark system.

Moreover, there exists a relatively narrow state, ${\rm H}_{c^{2}[nn][ns]}(5482,1/2,0^{+})$, whose total width is less than 5 MeV and it only decays to $\Lambda_{c}\Xi^{\prime}_{c}$ final states. Its relatively narrow width implies that its resonance peak is quite distinct, and the specific decay mode provides a clear signature for experimental searches. Therefore, the ${\rm H}_{c^{2}[nn][ns]}(5482,1/2,0^{+})$ has the potential to be discovered in experiments. Further, we suggest the relevant experiment check for the signal of ${\rm H}_{c^{2}[nn][ns]}(5482,1/2,0^{+})$ in the $5400$-$5500$ MeV energy window, and its lineshape should be relatively prominent in $\Lambda_{c}\Xi^{\prime}_{c}$ mass spectrum.

In addition, the ${\rm H}_{c^{2}[nn][ns]}(5487,1/2,1^{+})$ and ${\rm H}_{c^{2}[nn][ns]}(5261,1/2,0^{+})$ are also relatively narrow states, with total widths of 16 and 19 MeV, respectively. The ${\rm H}_{c^{2}[nn][ns]}(5261,1/2,0^{+})$ only decays to $\Xi_{c}\Lambda_{c}$ final states. For the ${\rm H}_{c^{2}[nn][ns]}(5487,1/2,1^{+})$, we obtain its relative partial decay width ratios as:

$$\frac{\Gamma[{\rm H}_{c^{2}[nn][ns]}(5487,1/2,1^{+})\to\Xi^{*}_{c}\Lambda_{c}]}{\Gamma[{\rm H}_{c^{2}[nn][ns]}(5487,1/2,1^{+})\to\Xi^{\prime}_{c}\Lambda_{c}]}=\frac{5}{3}.$$ (30)

For the $cc[ss][sn]$ subsystem, since it has exactly the same symmetry constraints as $cc[nn][ns]$ subsystem with $I=3/2$, the number of allowed states is also exactly identical. Due to the larger decay phase space and multiple different decay channels, most of the states belong to relatively broad states, with their decay widths ranging from 35 to 100 MeV. Among them, there is a relatively special state: ${\rm H}_{c^{2}[ss][sn]}(5806,1/2,2^{+})$, whose total width is 27 MeV. Therefore, it is the narrowest state in the $cc[ss][sn]$ subsystem. Moreover, it only decays into the $\Omega^{*}_{c}\Xi_{c}$ final states, which means there is a high possibility of observing it in the $\Omega^{*}_{c}\Xi_{c}$ decay channel. Based on the above research on the typical states, and with reference to Table 5 and Fig. 4, one can perform similar discussions on the decay behaviors of other $cc[ss][sn]$ states, and further explore their characteristics in-depth.

III.3 The $cc[nn][ss]$ and $cc[ns][ns]$ subsystems

Figure 5: Relative positions for the $cc[nn][ss]$ (a) and $cc[ns][ns]$ (b) hexaquark states labeled with horizontal solid lines, the labels, e.g. $5703~{}(\Gamma=48)$ represents the mass and total decay width of the corresponding state (units: MeV). In the $cc[nn][ss]$ and $cc[ns][ns]$ subsystems, the black, red, and purple horizontal lines represent the hexaquark states with $I=1$, $I=0$, and $I=1,0$, respectively. The notations are same as those of Fig. 3.

Table 6: The numerical results of the mass spectrum, the mass contributions of each Hamiltonian part (in MeV), the root-mean-square radii (in fm), and the partial decay widths and total decay widths of the fall-apart decay processes (in MeV) for the $cc[nn][ss]$ and $cc[ns][ns]$ hexaquark states.

$cc[nn][ss]$ Each part contribution RMS Radius Fall-apart decay properties 17  $I[J^{P}]$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Sigma^{*}_{c}\Omega^{*}_{c}$ $\Sigma^{*}_{c}\Omega_{c}$ $\Sigma_{c}\Omega^{*}_{c}$ $\Sigma_{c}\Omega_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $1[3^{+}]$ $|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=3}$ 5834 1365.7 -1622.2 19.7 1.25 1.25 1.36 1.31 1.25 1.41 42.2 42.2 $1[2^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=2}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=2}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=2}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=2}\end{pmatrix}$ $\begin{pmatrix}5820\\ 5819\\ 5726\end{pmatrix}$ 1373.4 -1629.9 5.7 1.25 1.25 1.36 1.30 1.40 1.25 4.5 13.4 39.7 57.6 1373.8 -1630.3 4.8 1.25 1.25 1.36 1.30 1.41 1.25 35.3 0.4 18.8 54.5 1366.0 -1715.5 4.2 1.44 1.45 1.32 1.25 1.10 1.44 14.3 9.3 10.7 34.3 $1[1^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=0}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=1}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=1}\end{pmatrix}$ $\begin{pmatrix}5811\\ 5810\\ 5809\\ 5730\end{pmatrix}$ 1378.5 -1635.0 -3.7 1.25 1.25 1.35 1.30 1.40 1.25 7.0 1.1 15.7 35.9 59.7 1378.9 -1635.4 -4.6 1.24 1.25 1.36 1.30 1.40 1.25 2.6 5.7 59.8 0.8 68.9 1379.2 -1635.6 -5.1 1.24 1.25 1.36 1.30 1.41 1.25 23.7 1.1 38.1 1.1 64.0 1363.7 -1713.2 8.3 1.45 1.45 1.32 1.25 1.10 1.45 23.4 2.7 3.5 4.0 33.6 $1[0^{+}]$ $\begin{pmatrix}|[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=0}\\ |[(nn)^{I=1,\bar{3}_{c}}_{s=1}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=0}\end{pmatrix}$ $\begin{pmatrix}5805\\ 5724\end{pmatrix}$ 1381.5 -1638.0 -9.3 1.24 1.25 1.35 1.30 1.40 1.25 28.5 28.2 56.7 1367.1 -1714.5 2.1 1.44 1.45 1.32 1.25 1.10 1.44 28.9 6.2 35.1 $I[J^{P}]$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Lambda_{c}\Omega^{*}_{c}$ $\Lambda_{c}\Omega_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $0[2^{+}]$ $|[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=2}$ 5598 1379.3 -1529.5 9.3 1.33 1.26 1.45 1.31 1.45 1.29 53.6 53.6 $0[1^{+}]$ $\begin{pmatrix}|[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=1}\\ |[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=0}_{s=1}\end{pmatrix}$ $\begin{pmatrix}5590\\ 5503\end{pmatrix}$ 1383.6 -1533.9 1.1 1.33 1.26 1.45 1.31 1.44 1.29 10.6 7.5 18.1 1375.7 -1617.5 6.1 1.52 1.47 1.43 1.24 1.13 1.48 18.8 4.0 22.8 $0[0^{+}]$ $|[(nn)^{I=0,\bar{3}_{c}}_{s=0}(ss)^{I=0,\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=0}$ 5586 1385.8 -1536.1 -3.1 1.33 1.26 1.45 1.30 1.44 1.29 6.5 6.5 $cc[ns][ns]$ Each part contribution RMS Radius Fall-apart decay properties 17  $I[J^{P}]$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Xi^{*}_{c}\Xi^{*}_{c}$ $\Xi^{*}_{c}\Xi^{\prime}_{c}$ $\Xi^{\prime}_{c}\Xi^{\prime}_{c}$ $\Xi_{c}\Xi_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $1[2^{+}]$ $\begin{pmatrix}|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=2}\\ |[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=2}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=2}\end{pmatrix}$ $\begin{pmatrix}5837\\ 5743\end{pmatrix}$ 1373.0 -1641.0 4.8 1.24 1.25 1.33 1.33 1.40 1.24 37.6 0.4 38.0 1364.7 -1725.8 4.3 1.44 1.45 1.28 1.28 1.10 1.44 15.3 8.9 24.2 $1(1^{+})$ $|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=1}$ 5828 1378.0 -1646.0 -4.3 1.24 1.25 1.32 1.32 1.40 1.24 2.4 6.1 0.7 9.2 $1[0^{+}]$ $\begin{pmatrix}|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=0}\\ |[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=0}\\ |[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=1,\bar{6}_{c}}_{s=0}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=0}\end{pmatrix}$ $\begin{pmatrix}5823\\ 5749\\ 5470\end{pmatrix}$ 1380.5 -1648.5 -8.9 1.23 1.25 1.32 1.32 1.40 1.24 27.6 29.5 57.1 1361.4 -1722.5 10.2 1.44 1.45 1.28 1.28 1.10 1.44 30.0 6.7 36.7 1369.1 -1643.4 6.1 1.50 1.47 1.32 1.32 1.13 1.48 26.7 26.7 $0[3^{+}]$ $|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=3}$ 5851 1365.2 -1633.2 19.3 1.24 1.25 1.33 1.33 1.41 1.25 47.1 47.1 $0[2^{+}]$ $|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=2}$ 5838 1372.7 -1640.7 5.6 1.24 1.25 1.33 1.33 1.40 1.25 4.2 13.4 17.6 $0[1^{+}]$ $\begin{pmatrix}|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=2}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=1}\\ |[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=0,3_{c}}_{s=0}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=1}\\ |[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=0,\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=0}_{s=1}\\ |[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}]^{I=0,3_{c}}_{s=0}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=0}_{s=1}\end{pmatrix}$ $\begin{pmatrix}5829\\ 5827\\ 5747\\ 5562\end{pmatrix}$ 1377.7 -1645.7 -3.6 1.24 1.25 1.32 1.32 1.40 1.24 6.6 1.1 36.4 44.1 1378.2 -1646.2 -4.7 1.23 1.25 1.32 1.32 1.40 1.24 28.6 1.4 1.3 31.3 1362.5 -1723.6 8.2 1.44 1.45 1.28 1.28 1.10 1.44 25.0 2.9 4.5 32.4 1378.2 -1559.3 5.2 1.30 1.26 1.36 1.36 1.44 1.28 41.5 41.5 $I(J^{P})$ Configuration Mass $\langle T\rangle$ $\langle V^{\rm Con}\rangle$ $\langle V^{\rm SS}\rangle$ $R_{12}$ $R_{34}$ $R_{13}$ $R_{23}$ $R_{12-34}$ $R_{13-24}$ $\Xi_{c}\Xi^{*}_{c}$ $\Xi_{c}\Xi^{\prime}_{c}$ $\Gamma_{sum}$ $R_{14}$ $R_{24}$ $R_{14-23}$ $1(0)[2^{+}]$ $|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=2}$ 5705 1372.9 -1596.9 9.8 1.27 1.26 1.37 1.33 1.42 1.26 68.4 68.4 $1(0)[1^{+}]$ $\begin{pmatrix}|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=1}\\ |[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,\bar{6}_{c}}_{s=1}(cc)^{I=0,6_{c}}_{s=0}\rangle^{I=1}_{s=1}\end{pmatrix}$ $\begin{pmatrix}5696\\ 5608\end{pmatrix}$ 1377.8 -1601.8 0.6 1.27 1.25 1.36 1.33 1.42 1.26 14.0 13.6 27.6 1366.5 -1683.7 6.2 1.47 1.46 1.33 1.27 1.11 1.46 9.6 9.2 18.8 $1(0)[0^{+}]$ $|[(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=0}(ns)^{I=\frac{1}{2},\bar{3}_{c}}_{s=1}]^{I=1,3_{c}}_{s=1}(cc)^{I=0,\bar{3}_{c}}_{s=1}\rangle^{I=1}_{s=0}$ 5691 1380.3 1604.3 -4.0 1.27 1.25 1.36 1.32 1.42 1.26 13.3 13.3

Finally, we discuss the $cc[nn][ss]$ and $cc[ns][ns]$ subsystems. These two subsystems both have the same quark contents, same quantum numbers, and same mass range. However, we can still distinguish them easily, because their decay final states show obvious differences. The $cc[nn][ss]$ states mainly decay to $\Sigma^{(*)}_{c}\Omega^{(*)}_{c}$ and $\Lambda_{c}\Omega^{(*)}_{c}$ final states, with the $\Xi^{(*)}_{c}\Xi^{(*)}_{c}$ decay channel being extremely suppressed. In contrast, the decay behavior of the $cc[ns][ns]$ states is completely opposite to that of the $cc[nn][ss]$ states.

For isovector $cc[nn][ss]$ states, they are all relatively broad states, with their total widths ranging from 30 to 70 MeV. Then for the four $I[J^{P}]=1[1^{+}]$ states, three of them, namely ${\rm H}_{c^{2}[nn][ss]}(5811,1,1^{+})$, ${\rm H}_{c^{2}[nn][ss]}(5810,1,1^{+})$, and ${\rm H}_{c^{2}[nn][ss]}(5809,1,1^{+})$, are partner states with similar masses and widths. Their total widths are 60, 69, and 64 MeV, respectively, and the mass gaps among them are only 1 MeV. We can distinguish them by the ratios of relative partial widths:

$$\Gamma_{\Sigma^{*}_{c}\Omega^{*}_{c}}:\Gamma_{\Sigma^{*}_{c}\Omega_{c}}:\Gamma_{\Sigma_{c}\Omega^{*}_{c}}:\Gamma_{\Sigma_{c}\Omega_{c}}=6:1:14:33,$$

$$\Gamma_{\Sigma^{*}_{c}\Omega^{*}_{c}}:\Gamma_{\Sigma^{*}_{c}\Omega_{c}}:\Gamma_{\Sigma_{c}\Omega^{*}_{c}}:\Gamma_{\Sigma_{c}\Omega_{c}}=1:2:23:0.3,$$

and

$$\Gamma_{\Sigma^{*}_{c}\Omega^{*}_{c}}:\Gamma_{\Sigma^{*}_{c}\Omega_{c}}:\Gamma_{\Sigma_{c}\Omega^{*}_{c}}:\Gamma_{\Sigma_{c}\Omega_{c}}=22:1:35:1,$$ (31)

respectively. From the above ratios, we notice that $\Sigma_{c}\Omega_{c}$ channel is the dominant decay channel for the ${\rm H}_{c^{2}[nn][ss]}(5811,1,1^{+})$. In contrast, $\Sigma_{c}\Omega_{c}$ channel is suppressed in the ${\rm H}_{c^{2}[nn][ss]}(5810,1,1^{+})$ and ${\rm H}_{c^{2}[nn][ss]}(5809,1,1^{+})$. They mainly decay to $\Sigma_{c}\Omega^{*}_{c}$ final states and $\Sigma_{c}\Omega^{*}_{c}$, $\Sigma^{*}_{c}\Omega^{*}_{c}$ final states, respectively.

For isoscalar states, there exists a narrow state, ${\rm H}_{c^{2}[nn][ss]}(5586,0,0^{+})$. It has a total width of about 7 MeV and only decays to $\Lambda_{c}\Omega_{c}$ final states. Although it has a larger decay phase space, the signs of the Feynman amplitudes $\mathcal{M}(A\to BC)$ (Eq. (15)) from the four quark-interchange diagrams (Fig. 2) are different for the $\Lambda_{c}\Omega_{c}$ decay channel. The contributions among them largely cancel out, leading to the suppression of the decay width. These characteristics, namely the narrow width and the unique $\Lambda_{c}\Omega_{c}$-only decay mode, are highly desirable in experimental searches. They make the lineshape of ${\rm H}_{c^{2}[nn][ss]}(5586,0,0^{+})$ more prominent in relevant experiments, thus increasing the likelihood of its discovery. Therefore, we suggest that experiments prioritize the search for possible resonance peaks in the $5500$-$5600$ MeV range of the $\Lambda_{c}\Omega_{c}$ invariant mass spectrum.

For the $cc[ns][ns]$ subsystem, the two $[ns]$ diquarks need to satisfy the Spin-Statistics theorem as identical particles in the $|[(ns)^{I=1/2}_{s=1}(ns)^{I=1/2}_{s=1}](cc)\rangle$ and $|[(ns)^{I=1/2}_{s=0}(ns)^{I=1/2}_{s=0}](cc)\rangle$ configurations. In the $|[(ns)^{I=1/2}_{s=0}(ns)^{I=1/2}_{s=1}](cc)\rangle$ configuration, since the two diquarks are not identical particles, the constraints imposed on the color-spin wave functions for their coupling to the $I=1$ state and the $I=0$ state are identical. As a result, they have exactly the same masses, RMS radii, and decay behaviors. Among them, ${\rm H}_{c^{2}[ns]^{2}}(5828,1,1^{+})$ is the state with the narrowest width. Although it has multiple rearrangement decay channels: $\Xi^{*}_{c}\Xi^{*}_{c}$, $\Xi^{*}_{c}\Xi^{\prime}_{c}$, and $\Xi^{\prime}_{c}\Xi^{\prime}_{c}$, its total width is still less than 10 MeV. The partial width ratio is:

$$\Gamma_{\Xi^{*}_{c}\Xi^{*}_{c}}:\Gamma_{\Xi^{*}_{c}\Xi^{\prime}_{c}}:\Gamma_{\Xi^{\prime}_{c}\Xi^{\prime}_{c}}=1:2.5:0.3.$$ (32)

Our results show that $\Xi^{*}_{c}\Xi^{\prime}_{c}$ is the dominant decay channel. We propose that future experimental investigations could explore the presence of the signal corresponding to ${\rm H}_{c^{2}[ns]^{2}}(5828,1,1^{+})$ within the mass range of 5800-5900 MeV in the $\Xi^{*}_{c}\Xi^{\prime}_{c}$ final states. For other $cc[ns][ns]$ states, one can perform similar discussions on the decay behaviors according to Table 6 and Fig. 5.

IV summary

In the field of hadron physics, the study of novel exotic states is a central research focus. Inspired by the observed $T_{cc}^{+}(3875)$, this study focuses on the doubly-charmed hexaquark system ($cc[qq][qq]$). We replace the two antiquarks $\bar{q}$ in the $T_{cc}^{+}(3875)$ ($cc\bar{q}\bar{q}$) state with two strongly correlated light diquarks $[qq]$, and conduct an in-depth study within the diquark picture.

In this research, we first construct total wave functions that satisfy the Spin-Statistics Theorem, covering the flavor, spatial, color, and spin parts. Then, within the framework of the constituent quark model, we systematically calculated the mass spectra of the doubly-charmed hexaquark system using the Gaussian expansion method. Finally, through detailed calculations, we also obtained corresponding internal mass contributions, root-mean-square radii, two-body strong decay partial widths, and total decay widths.

The mass spectra of the doubly-charmed hexaquark system are in the range of $5000$–$6000$ MeV, and the analysis of the mass spectra shows that there is no stable state in this system. All states can decay into two final-state singly-charmed baryons through two-body rearrangement strong decay. The analysis of the internal mass contributions shows that the contributions of the kinetic energy and the confinement potential are of the same order of magnitude, while the contribution of the hyperfine interaction potential is suppressed and relatively small, leading to low configuration mixing and small mass gaps between certain configurations, resulting in partner states.

Regarding the RMS radii, most results are in the range of $1.2$-$1.6$ fm, roughly in the same order of magnitude. This implies that the spatial distribution between quarks is relatively close and the internal interactions within the system are strong, which is consistent with the expectations of the compact doubly-charmed hexaquark configuration.

After that, we deeply analyzed the decay behavior of the doubly-charmed hexaquark system. The analysis indicates that most observed states possess total decay widths ranging from 30 to 100 MeV. However, there are still some narrow states. Even though these states generally have a large two-body decay phase space, the total widths of some of them are even less than 10 MeV. Upon further investigation of the reason, we found that for some specific decay channels, the signs of the Feynman amplitudes $\mathcal{M}(A\rightarrow BC)$ from the four quark-exchange diagrams are opposite. This causes their contributions to largely cancel each other out, ultimately leading to the suppression of the decay width.

In the doubly-charmed hexaquark system, the ${\rm H}_{c^{2}[nn]^{2}}(5043,0,0^{+})$ state is an ideal candidate for experimental detection. Composed of two scalar $[nn]$ diquarks ($I=0,S=0$), its structure enhances the internal confinement potential and chromomagnetic interaction, resulting in the lowest mass (5043 MeV) in the mass spectrum. It mainly decays to the $\Lambda_{c}\Lambda_{c}$ final state with a total width of 14.3 MeV, with a narrow hadronic resonance characteristic. Given the experimental detectability of doubly-charmed multi-quark systems confirmed by the discovery of $T_{cc}^{+}(3875)$, we suggest the LHCb, CMS, ATLAS collaborations to analyze the $\Lambda_{c}\Lambda_{c}$ invariant mass spectrum in the $5000$–$5100$ MeV range in high-luminosity $pp$ collisions to search for this state.

Furthermore, several narrow states also show great potential for experimental detection. The ${\rm H}_{c^{2}[nn][ns]}(5482,1/2,0^{+})$ state in the $cc[nn][ns]$ subsystem has a total width of less than 5 MeV and decays only to $\Lambda_{c}\Xi^{\prime}_{c}$, with high detectability in the $5400$-$5500$ MeV window. The ${\rm H}_{c^{2}[nn][ss]}(5586,0,0^{+})$ state in the $cc[nn][ss]$ subsystem has a total width of about 7 MeV and decays only to $\Lambda_{c}\Omega_{c}$. It may form a resonance peak in the 5500-5600 MeV range of the $\Lambda_{c}\Omega_{c}$ invariant mass spectrum. The ${\rm H}_{c^{2}[ns]^{2}}(5828,1,1^{+})$ state in the $cc[ns][ns]$ subsystem has a total width of less than 10 MeV and its dominant decay channel is $\Xi_{c}^{*}\Xi_{c}^{\prime}$, showing distinct lineshape compared to the background in the $5800$-$5900$ MeV range.

In summary, these results comprehensively reveal the mass spectra, internal structures, and decay characteristics of doubly-charmed hexaquark states. We hope they can provide some perspectives for further theoretical research. Meanwhile, we also look forward to more experimental collaborations focusing on doubly-charmed hexaquark states in the future. Conducting more experimental measurements can not only test our research results, but also deepen the understanding of the interactions within multi-quark systems.