Fully heavy dibaryons

The existence of fully heavy dibaryons ΩcccΩbbb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varOmega _{ccc}\varOmega _{bbb}$$\end{document}, ΩccbΩbbc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varOmega _{ccb}\varOmega _{bbc}$$\end{document}, ΩcccΩccc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varOmega _{ccc}\varOmega _{ccc}$$\end{document}, and ΩbbbΩbbb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varOmega _{bbb}\varOmega _{bbb}$$\end{document} with J=0,1,2,3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J=0,~1,~2,~3$$\end{document} and P=±1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P=\pm 1$$\end{document} are investigated in the framework of a constituent quark model with the help of the resonating group method. The dibaryon composed of six c or b quarks with JP=0+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^{P}=0^{+}$$\end{document} is able to be bound, because the requirement for antisymmetrization between the same baryon clusters introduces an attractive interaction between two fully heavy baryons. Although it is difficult for the dibaryon with the color-singlet type ΩcccΩbbb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varOmega _{ccc}\varOmega _{bbb}$$\end{document} to form any bound state, it is possible for the ΩccbΩbbc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varOmega _{ccb}\varOmega _{bbc}$$\end{document} state to be bound. The channel coupling of all channels of both ΩccbΩbbc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varOmega _{ccb}\varOmega _{bbc}$$\end{document} and ΩcccΩbbb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varOmega _{ccc}\varOmega _{bbb}$$\end{document} structures leads to the bound conclusion of this fully heavy system composed of three c quarks and three b quarks.


Introduction
Understanding hadron-hadron interactions and searching for exotic hadron states are important topics in hadron physics, among which questing for dibaryons is a long-standing challenge. Among the various theoretical studies of dibaryons, the lattice quantum chromodynamics (QCD) simulations by the HAL QCD (hadrons to atomic nuclei from lattice QCD) Collaboration have investigated six-quark systems containing light or strange quarks, and have confirmed the existence of the N Ω and ΩΩ bound states with nearly physical quark masses (m π 146 MeV and m K 525 MeV) [1,2]. Recently, Junnarkar and Mathur reported the first lattice QCD study of deuteron-like (np-like) dibaryons with heavy quark flavors, using a state-of-the-art lattice QCD calcua e-mail: hxhuang@njnu.edu.cn (corresponding author) b e-mail: jlping@njnu.edu.cn (corresponding author) c e-mail: zxm_yz@126.com d e-mail: fgwang@nju.edu.cn lation [3]. They suggested that dibaryons Ω c Ω cc (sscscc), Ω b Ω bb (ssbsbb), and Ω ccb Ω cbb (ccbcbb) were stable under strong and electromagnetic interactions, and they also found that the binding of these dibaryons became stronger as they became heavier in mass. A dibaryon with the highest charm number Ω ccc Ω ccc was also studied by the (2+1)-flavor lattice QCD with nearly physical light-quark masses and relativistic heavy-quark action with physical charm quark mass, and their result indicated that Ω ccc Ω ccc was located in the unitary regime [4]. The Ω bbb Ω bbb dibaryon was investigated in the extended one-boson exchange model, and the existence of Ω bbb Ω bbb was proposed [5]. Moreover, the possibility of very heavy dibaryons with three charm quarks and three beauty quarks (bbbccc) was explored using a constituent model in Ref. [6], and no bound state was found below the lowest dissociation threshold.
The discovery of the doubly charmed baryon Ξ cc by the Large Hadron Collider beauty (LHCb) Collaboration [7] provided the crucial experimental input, the interaction between two heavy quarks, for the existence of the stable heavy tetraquark Q Qqq [8]. Eichten and Quigg also predicted the existence of the novel narrow doubly heavy tetraquark Q Qqq based on the information [9]. Inspired by the LHCb Collaboration's observation of the hidden-charm P c pentaquarks [10], the work of Ref. [11] predicted the dibaryons 0 + Ξ cc Σ c , 1 + Ξ cc Σ * c , 2 + Ξ * cc Σ c , and 3 + Ξ * cc Σ * c . Recently, the LHCb Collaboration reported their results on the observations of fully charm states (cccc). A narrow structure X (6900), matching the lineshape of a resonance, and a broad structure next to the di − J/ψ mass threshold were obtained [12]. Such a breakthrough offers more information to assist the search for the tetraquark consisting of four charm quarks. It can naturally be expected that a fully heavy dibaryon cccccc will be accessible experimentally in the future.
Thus, the observation of the fully heavy dibaryons becomes more interesting. Based on the study of the strange dibaryon N Ω, we predicted the N Ω-like dibaryons N Ω ccc and N Ω bbb in the framework of the constituent quark model [13]. The study of the H particle was extended to the heavy sector states Λ c Λ c and Λ b Λ b [14]. In this work, we extend our study to the possibility of fully heavy dibaryon systems and investigate the interaction between two fully heavy baryons. The effective potential and binding energy, as well as the low-energy scattering phase shifts, scattering length, and effective range, are also calculated to explore the existence of fully heavy dibaryon systems.
The paper is structured as follows. Section 2 gives a brief introduction of the quark model and calculation method. The numerical results and discussion are provided in Sect. 3. A summary is presented in the last section.

Quark model and calculation method
We study the fully heavy dibaryon systems within the constituent quark model, which is similar to what was used in our previous work on dibaryons d * , N Ω, and ΩΩ [15][16][17], and is extended to the heavy dibaryons N Ω ccc [13]. When applied to the fully heavy systems, neither the SU (3) scalar octet meson exchange nor the Goldstone boson exchange works. Additionally, although the one-gluon-exchange interaction (V OGE i j ) has a tensor part in Ref. [13], it works in the coupling calculation between the S-wave and D-wave channels. We calculated the threshold energies of the S-wave states in both Ref. [13] and this work, and the coupling calculation between the S-wave and D-wave channels are tentatively not included in this work, as we are interested in the ground states of fully heavy dibaryons, and the spin-orbit and tensor interactions have very small contributions. For example, the tensor force reduced the energy of deuterons 1 or 2 MeV. For the fully heavy dibaryon, the mass shifts will be smaller, because the tensor force is inversely proportional to the masses of interacting quarks. In fact, we performed a calculation of the cccccc system including the tensor force, and the mass shift was less than 0.1 MeV, so it is safe to neglect the spin-orbit and tensor forces in the calculation of lowlying fully heavy dibaryons. To reduce the computational burden, the spin-orbit and tensor part of V OGE i j is not taken into account here. Here, we list the interaction Hamiltonian we used in this work: where T CM is the kinetic energy of the center of mass. V CON i j is the phenomenological confinement potential: It is well known that confinement is a nonperturbative QCD effect. There may be important nonperturbative features missed in every two-body confinement potential model. For example, three-gluon exchange between three quarks and three-body instanton interactions do not contribute within a colorless meson or baryon, but do contribute to a multiquark system [18]. There are many other types of quark interactions resulting from multi-gluon exchanges which cannot be included in a two-body confinement potential. The QCD condensates within a nucleon and between two nucleons might be different. These arguments show that a direct extension of the two-body confinement to the multiquark system has not been justified even though it might be a good approximation in the single hadron case. The absence of an experimentally observed color van der Waals force has long been a problem for the two-body confinement potential model. Multiquark systems provide greater variation in the low-energy behavior of QCD and allow for further testing of the phenomenology built into quark models, especially the nature of confinement.
We assume the following recipe to determine the twobody matrix element of confinement: the interaction takes the normal, unscreened form (quadratic in r) when the interacting quarks always remain in the same baryon orbit, both before and after interaction; otherwise the interaction takes the screening form (second form shown in Eq. (2)). Although this has not been demonstrated to be correct, it is more sophisticated than the usual simple two-body confining interaction, and it is a physically reasonable model of nonlocal, nonperturbative effects. Even though we use a two-body interaction form to evaluate the matrix elements, our model is not a potential model. It is an effective matrix element approach extended from bound states to scattering states. It does reduce to the usual two-body confinement interaction within a single hadron. It also has the usual meaning of a two-body interaction in the asymptotic regime, although not in intermediate regions. The main physics introduced is the recognition that the confining interaction between two nucleons might differ from that within a nucleon. In particular, we represent the nonlocal, nonperturbative backflow of color (and pair creation, as seen in lattice studies [19,20]) by screening of the confining potential.
In Eq. (2), μ q i q j is the color screening parameter, which is related to the flavor of the i-th quark q i and the j-th For the light quark systems, μ uu , μ us , and μ ss are determined by fitting the deuteron properties, nucleonnucleon scattering phase shifts, and hyperon-nucleon scattering phase shifts, respectively [21]. We also found that the heavier the quark, the smaller the parameter μ q i q j . When extending to the heavy quark system, the hidden-charm pentaquark system, we took μ cc as an adjustable parameter from 0.01 to 0.001 fm −2 , and found that the results were insensitive to the value of μ cc [22], and the P c states were well described in the work of Refs. [22,23]. In fact, if the value of the μ q i q j is very small, the exponential function can be approximated to Then, the confinement expression in Eq. (2) (when two quarks are in different clusters) is which is the same as the expression for two quarks in the same cluster. Thus, when the value of the μ q i q j is very small, the screened confinement will return to the quadratic form. This is why the results are insensitive to the value of μ cc and μ bb . Therefore, we can use the quadratic form of the confinement expression for the fully heavy dibaryon systems.
is the one-gluon-exchange interaction: where α s (μ) is the quark-gluon coupling constant. In Ref. [24], it is written as where μ is the reduced mass of the interacting quark pair. The other symbols in the above expressions have their usual meanings. All parameters are fixed by fitting to the masses of Table 2 The masses (in MeV) of the light, charmed, and bottomed baryons. Experimental values are taken from the Particle Data Group (PDG) [32] Baryon  28] are used to carry out a dynamical calculation. The main feature of the RGM for two-cluster systems is that it assumes that two clusters are frozen inside, and only considers the relative motion between the two clusters. Thus, the conventional ansatz for the twocluster wavefunctions is where the symbol A is the antisymmetrization operator.
gives the total color symmetry, and all other symbols have their usual meanings. φ B i is the three-quark cluster wavefunction, and χ L (R) is the relative motion wavefunction, which is expanded by a set of Gaussians with different centers [29], where L is the orbital angular momentum between two clusters. n is the number of Gaussians used in the expansion, which is fixed by the requirement that the results are stable against the increasing of n. By including the center-of-mass motion the ansatz Eq. (7) can be rewritten as where φ α (S i ) and φ β (−S i ) are the single-particle orbital wavefunctions with different reference centers: In the present work, to reduce the computational burden, perturbation theory is employed, and the zeroth-order wavefunctions are used to obtain the energy of the system in first order. The relative motion and center-of-mass motion with different masses of quarks are calculated, taking proper care.
We used this method to study the dibaryon N Ω [16], and the result was consistent with the lattice QCD [30]. The experimental result from the STAR Collaboration also supported our prediction of N Ω resonance [31].
From the variational principle, we obtain the generalized eigenvalue equation where H i, j and N i, j are the Hamiltonian matrix elements and overlaps, respectively. By solving the generalized eigen problem, we can obtain the energy and the corresponding wavefunctions of the dibaryon system.

Results and discussion
As the first step, we calculate the energy of the color-singlet states (two clusters with color singlet) with the quantum numbers J = 0, 1, 2, and 3. The results are listed in Table  3. Here, we should mention that in the RGM, the internal motions of clusters are frozen, and the relative motion wavefunction is expanded by a set of Gaussians (see Eq. (8)). So, in this method, the orbital angular momentum between two clusters is also the total orbital angular momentum of the system. For the Ω ccc Ω ccc or Ω bbb Ω bbb system, which includes six identical quarks, totally antisymmetric wavefunctions must be used. The total symmetry is determined by the symmetry of four degrees of freedom (color, isospin, spin, and orbital). Here, the symmetry of both color and isospin is symmetrical; the spin (S) is symmetrical for S = 1 or 3, and antisymmetric for S = 0 or 2. To obtain the totally antisymmetric wavefunction, the orbital angular momentum L should be odd for the S = 1 or 3, and even for S = 0 or 2. If we consider the lower partial wave, we can take L = 1 for S = 1 or 3, and L = 0 for S = 0 or 2. Then, we can obtain the total angular momentum J = 0, 1, 2, 3, 4 with negative parity and J = 0, 2 with positive parity. Also, the spin-orbit interaction is not included in the present calculation. So, for the case of L = 1, S = 1, the results of J P = 0 − are the same with J P = 1 − and 2 − , and for the case of L , = 1, S = 3, the results of J P = 2 − are the same with J P = 3 − and 4 − . Therefore, we only show the results of J P = 1 − and J P = 3 − in Table 3. However, nonidentical quarks in two subclusters of the system Ω ccc Ω bbb leading to both positive and negative parity are possible for this system with J = 0, 1, 2, and 3. From Table 3, one can see that the energies of all states are above the corresponding theoretical threshold, except the states Ω ccc Ω ccc , Ω bbb Ω bbb , and Ω ccb Ω bbc with J P = 0 + , whose energy is −2.5 MeV, −0.9 MeV, and −1.1 MeV lower than the corresponding thresholds, respectively. Therefore, for the J P = 0 + systems, the dibaryon with six c quarks or six b quarks, or with the ccb − bbc structure, can form bound states, while the one with the ccc −bbb structure is unbound. This is due to the different symmetry requirements for the different dibaryon structures. An antisymmetrization operator is required for the systems Ω ccc Ω ccc , Ω bbb Ω bbb , and Ω ccb Ω bbc , while there is no symmetry requirement between the clusters Ω ccc and Ω bbb . We also note that the energies of the states Ω ccc Ω bbb with different quantum numbers are the same. This is mainly because there is no interaction between two color-singlet subclusters because there is no exchange term. In this case, the color-related interaction does not work between (ccc) and (bbb), so the energies with different total angular momenta are the same. Besides, the spin-orbit interaction is not included in present calculation, and the states with J P = 0 + (S = 0, L = 0) and J P = 0 − (S = 1, L = 1) have the same energy, because the two clusters are well separated to minimize the kinetic energy. However, for the hiddencolor channels (two clusters with a color octet), the colorrelated interaction works between (ccc) and (bbb), so the results will be different for these states. The results of the hidden-color channels and the channel-coupling calculation will be shown later.
The effective potentials are also calculated to understand the interaction between two fully heavy baryons, which are shown in Figs. 1 and 2

. The effective potential between two colorless clusters is defined as V (S) = E(S)− E(∞)
, where E(S) is the diagonal matrix element of the Hamiltonian of the system in the generating coordinate. It is clear that the potentials for the Ω ccc Ω ccc , Ω bbb Ω bbb , and Ω ccb Ω bbc with J P = 0 + are attractive, which leads to the possibility for these three channels to form bound states. However, the potentials for other channels are all repulsive, and therefore the other states are unbound.
To investigate the interactions between two fully heavy baryons in detail, we calculate the contributions to the effective potential from the kinetic energy (V vk ), the confinement (V con ), and the one-gluon exchange (V oge ) of the Ω ccc Ω ccc (a) (b) Fig. 1 The effective potential of the Ω ccc Ω ccc and Ω bbb Ω bbb systems (a) (b) Fig. 2 The effective potential of the Ω ccb Ω bbc and Ω ccc Ω bbb systems and Ω ccc Ω bbb systems with J P = 0 + , which are shown in Fig. 3. One can see that there is no contribution of color confinement interaction between two color-singlet clusters. For the Ω ccc Ω bbb system, the contribution of the one-gluon exchange interaction between Ω ccc and Ω bbb is also zero, because there is no quark exchange between these two baryons. Although the one-gluon exchange interaction works for the Ω ccc Ω ccc , it is very small due to the small quarkgluon coupling constant α s (μ) used for the fully heavy system. Besides, the effective potential of the kinetic energy shows that the kinetic energy term provides attractive interactions for the Ω ccc Ω ccc system, while it provides repulsive interactions for the Ω ccc Ω bbb system. This is due to the different symmetry requirement for the different dibaryon structures. At finite separation, the kinetic energy of the relative motion between two clusters generally makes a positive contribution. However, the antisymmetrization operator for (a) (b) Fig. 3 The contributions to the effective potential from the kinetic energy (V vk ), the confinement (V con ), and the one-gluon exchange (V oge ) of the Ω ccc Ω ccc and Ω ccc Ω bbb systems with J P = 0 + the Ω ccc Ω ccc and Ω bbb Ω bbb systems causes the exchange of quarks between two clusters, which increases the motion space of each quark for the six identical quarks systems, and the internal kinetic energy of each cluster is reduced. Thus, the total kinetic energy of the six-quark system may be smaller than the sum of the internal kinetic energies of two free three-quark systems, although there is relative kinetic energy between the two clusters. Therefore, it is possible that the contribution of the kinetic energy term for the Ω ccc Ω ccc system will be negative. For the Ω ccc Ω bbb system, there is no quark exchange between Ω ccc and Ω bbb , so the contribution of the kinetic energy term for the Ω ccc Ω bbb system is positive.
In the above calculation, we only consider the colorsinglet channels of the systems. However, the hidden-color channels and the channel-coupling effect are important in multiquark systems. Therefore, we calculate the hiddencolor channels and the channel coupling of the systems. The results are presented in Tables 4 and 5, in which the binding energy is obtained by subtracting the threshold of the corresponding channel from the energy of the system. E sc1 /B sc1 and E sc2 /B sc2 represent the energies/binding energies of the color-singlet channels and hidden-color channels, respectively; E cc /B cc represents the energies/binding energies of the channel coupling of both the color-singlet and hiddencolor channels; E cc /B cc represents the energies/binding energies of the channel coupling of all channels of the Ω ccb Ω bbc and Ω ccc Ω bbb systems; and ub means that the energy is above the threshold of the corresponding channel.
We can see from Tables 4 and 5 that the hidden-color channels have lower energies than the corresponding color-singlet channels, and the channel coupling of the color-singlet and hidden-color channels will increase the binding energy of the systems considerably. For the cccccc system with J P = 0 + , the binding energy in the single-channel calculation is −2.5 MeV for the color-singlet channel and −14.7 MeV for the hidden-color channel, and it will reach −30.1 MeV in the channel-coupling calculation. For the bbbbbb system, the results are similar to those for the cccccc system, −0.9 MeV for the color-singlet channel, −1.3 MeV for the hidden-color channel, and −7.6 MeV in the channel-coupling calculation. For the cccbbb system, the state with the Ω ccc Ω bbb structure is unbound, but the state with Ω ccb Ω bbc is bound, and the binding energies have similar behavior as that of the cccccc system, −1.1 MeV for the color-singlet channel, −11.6 MeV for the hidden-color channel, and −35.6 MeV in the channel-coupling calculation. These two structures have the same quark content, and they can be coupled. The structure-coupling calculation obtains binding energy of −80.1 MeV for the cccbbb system with J P = 0 + . Thus, the fully heavy dibaryons with J P = 0 + are bound states. In determining the model parameters, the same b, which relates to the size of baryon, is used for all baryons. This is not reasonable, because the light baryons and the fully heavy baryons differ in size. A value of b = 0.518 is good for describing the properties of the proton, while smaller b Fig. 4 The phase shifts of the S-wave Ω ccc Ω ccc and Ω bbb Ω bbb systems with J P = 0 + is expected to describe the fully heavy baryons. To test the dependence of the result on the baryon size parameter b, we calculate the fully heavy six-quark systems with b = 0.3 fm. The results are as follows. The binding energy for the cccccc system with J P = 0 + is −56 MeV; it is −8.0 MeV for the bbbbbb system, and −15.4 MeV for the cccbbb system. Therefore, the existence of fully heavy dibaryons is weakly dependent on the model parameters.
To further check the possibility of the bound state, we can also study the low-energy scattering phase shifts of the state. Here, we calculate the low-energy scattering phase shifts of the color-singlet channel of the S−wave Ω ccc Ω ccc and Ω bbb Ω bbb states, which are shown in Fig. 4. The welldeveloped Kohn-Hulthen-Kato (KHK) variational method is used here, the details of which can be found in Ref. [26]. It is obvious that the scattering phase shifts of the two systems reach 180 degrees at E c.m. ∼ 0 and rapidly decrease as E c.m. increases, suggesting the presence of a bound state.
Then, we can extract the low-energy scattering parameters, scattering length a 0 , and the effective range r 0 for both Ω ccc Ω ccc and Ω bbb Ω bbb systems using the formula where δ is the low-energy phase shifts obtained above, k is the momentum of relative motion with k = √ 2μE c.m. , in which μ is the reduced mass of two baryons, and E c.m. is the incident energy. The binding energy B can be calculated according to the relation where α is the wave number which can be obtained from the relation [33] For the bound state, the wave number α is the inverse of the spatial dimension of the system [34], and the required solution of Eq. (15) can be fixed accordingly. Please note that we use another method to calculate the binding energy here, so we label it B . The results are listed in Table 6, from which we can see that the scattering length is positive for both dibaryons Ω ccc Ω ccc and Ω bbb Ω bbb , which implies the existence of a bound state of these two fully heavy dibaryons. The binding energies obtained here are coincident with those calculated from Table 3.

Summary
The main conclusion of our dynamical investigation of fully heavy dibaryons is that the dibaryons composed of six c quarks, six b quarks, or three c quarks and three b quarks can form bound states, and the quantum number is J P = 0 + . However, the nature of these systems differs. For the Ω ccc Ω ccc or Ω bbb Ω bbb , the requirement of antisymmetrization between the same baryon clusters introduces attractive interaction between two fully heavy baryons, which leads to a super-heavy bound dibaryon Ω ccc Ω ccc or Ω bbb Ω bbb . The channel coupling between the color-singlet and hidden-color channels causes the binding energy to increase.
For the system with three c quarks and three b quarks, there is no evidence for any stable dibaryon with the structure of color-singlet Ω ccc Ω bbb . This outcome is based on the concept of "QCD-inspired" constituent quark model calculation, where the interactions among quarks are color-related. The principal reason for the instability of the Ω ccc Ω bbb system is that there is no symmetry requirement between the subclusters Ω ccc and Ω bbb because quark c and quark b are not identical quarks. Therefore, the interaction between two subclusters is zero. In contrast, for the Ω ccb Ω bbc , the exchange of identical quarks can occur between two subclusters, resulting in nonzero exchange terms of the interaction, and thus the formation of a bound state is possible. The channel coupling of all channels of both Ω ccb Ω bbc and Ω ccc Ω bbb structures leads to the bound state of this fully heavy system, which is composed of three c quarks and three b quarks.
The study of the interaction between two fully heavy baryons supports this conclusion, because the effective potentials are attractive. The behavior of low-energy scattering phase shifts and the scattering length also confirm the existence of stable fully heavy dibaryons Ω ccc Ω ccc and Ω bbb Ω bbb . The lattice QCD was used to study the highest charm dibaryon Ω ccc Ω ccc , and indicated the existence of this state [4]. We suggest that the lattice QCD indicates the existence of the fully heavy dibaryon with six b quarks, and the one with three c quarks and three b quarks.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The data have been illustrated in the figures and tables, so they are not necessary to be deposited. Data may be made available upon request.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.   The wavefunction ψ is the six-quark one in Eq.(10) without the antisymmetrization operator. The antisymmetrization operator A = 1 − A , where A = 9P 36 for Ω ccc Ω ccc and A = P 16 + P 26 + P 34 + P 35 − P 16 P 34 − P 16 P 35 − P 26 P 34 − P 26 P 35 for Ω ccb Ω bbc . For the Ω ccc Ω ccc system, the direct and exchange terms of overlap for the color-singlet channel are N D = 1, N E = −1; for the hidden-color channel they are N D = 1, N E = −7; and the cross ones between color-singlet and hidden-color channels are N D = 0, N E = −4. For the Ω ccb Ω bbc system, the direct and exchange terms of overlap for the color-singlet channel are N D = 1, N E = −8/9; for the hidden-color channel they are N D = 1, N E = −56/9; and the cross ones between the color-singlet and hidden-color channels are N D = 0, N E = −32/9.