Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Study of the ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT dibaryons in constituent quark model

Pablo Martín-Higueras pablo.higueras@alu.uhu.es Departamento de Ciencias Integradas y Centro de Estudios Avanzados en Física, Matemática y Computación, Universidad de Huelva, 21071 Huelva, Spain    David R. Entem entem@usal.es Grupo de Física Nuclear, Universidad de Salamanca, E-37008 Salamanca, Spain Instituto Universitario de Física Fundamental y Matemáticas (IUFFyM), Universidad de Salamanca, E-37008 Salamanca, Spain    Pablo G. Ortega pgortega@usal.es Instituto Universitario de Física Fundamental y Matemáticas (IUFFyM), Universidad de Salamanca, E-37008 Salamanca, Spain    Jorge Segovia jsegovia@upo.es Departamento de Sistemas Físicos, Químicos y Naturales,
Universidad Pablo de Olavide, E-41013 Sevilla, Spain
   Francisco Fernández fdz@usal.es Grupo de Física Nuclear, Universidad de Salamanca, E-37008 Salamanca, Spain Instituto Universitario de Física Fundamental y Matemáticas (IUFFyM), Universidad de Salamanca, E-37008 Salamanca, Spain
(June 27, 2024)
Abstract

Dibaryons are the simplest system in which the baryon-baryon interaction, and hence the underlying quark-quark interaction, can be studied in a clear way. Although the only dibaryon known today is the deuteron (and possibly the dsuperscript𝑑d^{*}italic_d start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT), fully heavy dibaryons are good candidates for bound states because in such systems the kinetic energy is small and the high symmetry of the wave function favours binding. In this study, the possible existence of ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT dibaryons is investigated in the framework of a constituent quark model that satisfactorily describes the deuteron, the d(2380)superscript𝑑2380d^{*}(2380)italic_d start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( 2380 ) and the NN𝑁𝑁NNitalic_N italic_N interaction. JP=0+superscript𝐽𝑃superscript0J^{P}=0^{+}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT candidates are found in both systems with binding energies of the order of MeV.

Potential models, Quark models, Bottom charmed mesons, Exotic mesons
pacs:
12.39.Pn, 14.40.Lb, 14.40.Rt

I Introduction

Understanding the nucleon-nucleon interaction has been one of the priority problems in Nuclear Physics since Yukawa’s one pion exchange theory. The subsequent development of QCD paved the way to describe the strong interactions in terms of quark degrees of freedom and facilitate to enlarge the field to other flavors like charm en bottom.

Dibaryons are the simplest systems in which these studies can be addressed in a transparent way. Until recently, the only well-established bound state of two baryons was the deuteron. Then, in 2011, another unstable light dibaryon, the d(2380)superscript𝑑2380d^{*}(2380)italic_d start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( 2380 ), was reported by the WASA-at-COSY Collaboration Adlarson et al. (2011) from the double pionic fusion reaction pndπ0π0𝑝𝑛𝑑superscript𝜋0superscript𝜋0pn\to d\pi^{0}\pi^{0}italic_p italic_n → italic_d italic_π start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT italic_π start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT. This resonance can be described as a nonstrange ΔΔΔΔ\Delta\Deltaroman_Δ roman_Δ dibaryon with I(JP)=0(3+)𝐼superscript𝐽𝑃0superscript3I(J^{P})=0(3^{+})italic_I ( italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT ) = 0 ( 3 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ). In 1989, Goldman noted that due to the special symmetry of such a state, any model based on confinement and gluon exchange should predict it Goldman et al. (1989). The long history of the search for dibaryons in the light quark sector can be found in Ref Clement (2017).

It is well known that the binding of the deuteron is due to the coupling of the S13superscriptsubscript𝑆13{}^{3}S_{1}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and D13superscriptsubscript𝐷13{}^{3}D_{1}start_FLOATSUPERSCRIPT 3 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT partial waves by one-pion exchange tensor interactions. Similarly, the binding of the d(2380)superscript𝑑2380d^{*}(2380)italic_d start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( 2380 ) can be explained in terms of Goldstone-boson exchanges Huang et al. (2014). These two systems then prove that the interaction binding these dibaryons arises from QCD chiral symmetry breaking in the light quark sector.

Another interesting system is the fully heavy dibaryon. In such a system the relativistic effects are negligible and the kinetic energy is small. As originally pointed out by Bjorken Bjorken (1985), the triply-charmed baryon ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT is stable against strong interactions. This fact opens the possibility to study systems like ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT or ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT. Moreover, in contrast to the deuteron and the dsuperscript𝑑d^{*}italic_d start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT case, the latter systems provide an ideal scenario to explore the baryon-baryon interaction in an environment free of chiral dynamics.

In this work we will focus on the study of the fully heavy dibaryons. Two recent Lattice QCD calculations have explored these systems: Ref. Lyu et al. (2021) showed that ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT is loosely bound by 5.68(0.77)5.680.775.68(0.77)5.68 ( 0.77 ) MeV, while Ref. Mathur et al. (2023) found a very deep ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT state with a binding energy of 8116+14superscriptsubscript81161481_{-16}^{+14}81 start_POSTSUBSCRIPT - 16 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 14 end_POSTSUPERSCRIPT MeV. These conclusions are confirmed by several quark model calculations, but are contradicted by others. For example, Huang et al. Huang et al. (2022), using a constituent quark model based on the one-gluon exchange interaction and the resonating group method, studied the possible bound states of the ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT, among others. They found a JP=0+superscript𝐽𝑃superscript0J^{P}=0^{+}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT bound state for the ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT system with a binding energy of 2.52.52.52.5 MeV and another ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT state bound by 0.90.90.90.9 MeV, contrary to naive espectations. Deng Deng (2023) performed a study of the di-Δ++superscriptΔabsent\Delta^{++}roman_Δ start_POSTSUPERSCRIPT + + end_POSTSUPERSCRIPT, di-ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and di-ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT systems using a naive one-gluon exchange quark model and a chiral quark model including π𝜋\piitalic_π and σ𝜎\sigmaitalic_σ exchanges between quarks. Obviously, in this case these parts of the interaction apply only to the light quarks, but the set of parameters is different in the two models. Both studies predict very shallow di-ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and di-ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT states with binding energies around 1111 MeV. Using a different model, namely QCD sum rules, Wang Wang (2022) found for each di-ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and di-ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT systems two JP=0+superscript𝐽𝑃superscript0J^{P}=0^{+}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and JP=1superscript𝐽𝑃superscript1J^{P}=1^{-}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 1 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT states that are slightly below their respective thresholds.

On the other hand, several studies within the quark model have ruled out the existence of fully heavy dibaryons. In Ref Richard et al. (2020) the authors investigated the existence of bbbccc𝑏𝑏𝑏𝑐𝑐𝑐bbbcccitalic_b italic_b italic_b italic_c italic_c italic_c dibaryons and extrapolated their results to the properties of the bbbbbb𝑏𝑏𝑏𝑏𝑏𝑏bbbbbbitalic_b italic_b italic_b italic_b italic_b italic_b and cccccc𝑐𝑐𝑐𝑐𝑐𝑐ccccccitalic_c italic_c italic_c italic_c italic_c italic_c systems. They found no bound states for ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT or ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT combinations. On the other hand, Alcaraz-Peregrina et al. Alcaraz-Pelegrina and Gordillo (2022) used the Difussion Montecarlo technique to describe fully heavy compact six-quark arrangements. They found that all the hexaquarks have smaller masses than those of their constituents, i.e., all the hexaquarks are bound systems. However, their masses are also larger than those of any pair of baryons into which they can be divided. This means that each hexaquark is unstable with respect to its splitting into two baryons. Finally, two more calculations, in the framework of the constituent quark model Lü et al. (2022) or the extended chromomagnetic model Weng and Zhu (2024), showed that all the fully heavy dibaryons lie above their corresponding baryon-baryon thresholds.

In view of this controversial situation, since different approaches lead to quite different conclusions, we will study the possible existence of ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT dibaryons using the constituent quark model of Ref. Fernandez et al. (1993) and its extension to the heavy quark sector Vijande et al. (2005); Segovia et al. (2008), which has been able to describe a large variety of hadronic phenomenology. In particular, the model reproduces the properties of the deuteron Valcarce et al. (1994); Juliá-Díaz et al. (2002) and predicts the existence of the d(2380)superscript𝑑2380d^{*}(2380)italic_d start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( 2380 ) as a ΔΔΔΔ\Delta\Deltaroman_Δ roman_Δ dibaryon Garcilazo et al. (1997); Valcarce et al. (2001). Although the binding energy of the dsuperscript𝑑d^{*}italic_d start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT predicted in the latter references is smaller than the experimental value, it is laso worth mentioning that the these studies were performed without coupling to the NN𝑁𝑁NNitalic_N italic_N channel.

The paper is structured as follows. In Sec. II we describe the main aspects of our theoretical model, giving details about the wave functions used to describe ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT (ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT) baryons and the way we derive the ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT interaction using the Resonating Group Method (RGM). Section IV is devoted to presenting our results for the possible dibaryons. Finally, we summarize and give some conclusions in Sec. V.

II Theoretical formalism

II.1 The constituent quark model

Our theoretical framework is a QCD-inspired constituent quark model (CQM) proposed in Ref. Vijande et al. (2005) and extended to the heavy quark sector in Ref. Segovia et al. (2008). The main pieces of the model are the constituent light quark masses and Goldstone-boson exchanges, which appears as consequences of spontaneous chiral symmetry breaking of the QCD Lagrangian together with perturbative one-gluon exchange (OGE) and nonperturbative color confining interactions.

Following Diakonov Diakonov (2003), a simple Lagrangian invariant under chiral transformations can be written as

=ψ¯(i/M(q2)Uγ5)ψ,¯𝜓𝑖𝑀superscript𝑞2superscript𝑈subscript𝛾5𝜓{\mathcal{L}}=\bar{\psi}(i\,{/\!\!\!\partial}-M(q^{2})U^{\gamma_{5}})\,\psi\,,caligraphic_L = over¯ start_ARG italic_ψ end_ARG ( italic_i / ∂ - italic_M ( italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_U start_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ) italic_ψ , (1)

where M(q2)𝑀superscript𝑞2M(q^{2})italic_M ( italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) is the dynamical (constituent) quark mass and Uγ5=eiλaϕaγ5/fπsuperscript𝑈subscript𝛾5superscript𝑒𝑖subscript𝜆𝑎superscriptitalic-ϕ𝑎subscript𝛾5subscript𝑓𝜋U^{\gamma_{5}}=e^{i\lambda_{a}\phi^{a}\gamma_{5}/f_{\pi}}italic_U start_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT = italic_e start_POSTSUPERSCRIPT italic_i italic_λ start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT italic_ϕ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT / italic_f start_POSTSUBSCRIPT italic_π end_POSTSUBSCRIPT end_POSTSUPERSCRIPT is the matrix of Goldstone-boson fields that can be expanded as

Uγ5=1+ifπγ5λaπa12fπ2πaπa+superscript𝑈subscript𝛾51𝑖subscript𝑓𝜋superscript𝛾5superscript𝜆𝑎superscript𝜋𝑎12superscriptsubscript𝑓𝜋2superscript𝜋𝑎superscript𝜋𝑎U^{\gamma_{5}}=1+\frac{i}{f_{\pi}}\gamma^{5}\lambda^{a}\pi^{a}-\frac{1}{2f_{% \pi}^{2}}\pi^{a}\pi^{a}+\ldotsitalic_U start_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT 5 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT = 1 + divide start_ARG italic_i end_ARG start_ARG italic_f start_POSTSUBSCRIPT italic_π end_POSTSUBSCRIPT end_ARG italic_γ start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT italic_λ start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_π start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 italic_f start_POSTSUBSCRIPT italic_π end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_π start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT italic_π start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT + … (2)

The first term of the expansion generates the constituent quark mass, while the second term gives rise to a one-boson exchange interaction between quarks. The main contribution of the third term comes from the two-pion exchange which has been simulated by means of a scalar-meson exchange potential.

In the heavy quark sector, chiral symmetry is explicitly broken and Goldstone-boson exchange does not occur. However, the full interaction constrains the model parameters through the light-meson phenomenology Vijande et al. (2005); Segovia et al. (2008). Thus, OGE and confinement are the only remaining interactions between the heavy quarks.

The OGE potential is generated from the vertex Lagrangian

qqg=i4παsψ¯γμGcμλcψ,subscript𝑞𝑞𝑔𝑖4𝜋subscript𝛼𝑠¯𝜓subscript𝛾𝜇subscriptsuperscript𝐺𝜇𝑐superscript𝜆𝑐𝜓{\mathcal{L}}_{qqg}=i\sqrt{4\pi\alpha_{s}}\,\bar{\psi}\gamma_{\mu}G^{\mu}_{c}% \lambda^{c}\psi,caligraphic_L start_POSTSUBSCRIPT italic_q italic_q italic_g end_POSTSUBSCRIPT = italic_i square-root start_ARG 4 italic_π italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG over¯ start_ARG italic_ψ end_ARG italic_γ start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_G start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_λ start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT italic_ψ , (3)

where λcsuperscript𝜆𝑐\lambda^{c}italic_λ start_POSTSUPERSCRIPT italic_c end_POSTSUPERSCRIPT are the SU(3)𝑆𝑈3SU(3)italic_S italic_U ( 3 ) colour matrices, Gcμsubscriptsuperscript𝐺𝜇𝑐G^{\mu}_{c}italic_G start_POSTSUPERSCRIPT italic_μ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT is the gluon field and αssubscript𝛼𝑠\alpha_{s}italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT is the strong coupling constant. The scale dependence of αssubscript𝛼𝑠\alpha_{s}italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT allows a consistent description of light, strange and heavy mesons. Its explicit expression can be found in, e.g., Ref. Vijande et al. (2005),

αs(μ)=α0ln(μ2+μ02Λ02)subscript𝛼𝑠𝜇subscript𝛼0superscript𝜇2superscriptsubscript𝜇02superscriptsubscriptΛ02\alpha_{s}(\mu)=\frac{\alpha_{0}}{\ln\left(\frac{\mu^{2}+\mu_{0}^{2}}{\Lambda_% {0}^{2}}\right)}italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_μ ) = divide start_ARG italic_α start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG roman_ln ( divide start_ARG italic_μ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) end_ARG (4)

Regarding the confinement potential, it is well known that multi-gluon exchanges produce an attractive linearly rising potential proportional to the distance between infinite-heavy quarks Mateu et al. (2019). However, sea quarks are also important components of the strong interaction dynamics that contribute to the screening of the rising potential at low momenta and eventually to the breaking of the quark-antiquark binding string Bali et al. (2005). Our model tries to mimic this behaviour with a screening potential at high distances.

Then, the full interaction between heavy quarks is given by

Vij(r)subscript𝑉𝑖𝑗𝑟\displaystyle V_{ij}(r)italic_V start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ( italic_r ) =[ac(1eμcr)+Δ+αs(μ)41r](λiλj)absentdelimited-[]subscript𝑎𝑐1superscript𝑒subscript𝜇𝑐𝑟Δsubscript𝛼𝑠𝜇41𝑟subscript𝜆𝑖subscript𝜆𝑗\displaystyle=\bigg{[}-a_{c}(1-e^{-\mu_{c}r})+\Delta+\frac{\alpha_{s}(\mu)}{4}% \frac{1}{r}\bigg{]}(\vec{\lambda}_{i}\cdot\vec{\lambda}_{j})= [ - italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ( 1 - italic_e start_POSTSUPERSCRIPT - italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_r end_POSTSUPERSCRIPT ) + roman_Δ + divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_μ ) end_ARG start_ARG 4 end_ARG divide start_ARG 1 end_ARG start_ARG italic_r end_ARG ] ( over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT )
VijS(r)superscriptsubscript𝑉𝑖𝑗𝑆𝑟\displaystyle V_{ij}^{S}(r)italic_V start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_S end_POSTSUPERSCRIPT ( italic_r ) =αs(μ)416mimjer/r0(μ)rr02(μ)(σiσj)(λiλj)absentsubscript𝛼𝑠𝜇416subscript𝑚𝑖subscript𝑚𝑗superscript𝑒𝑟subscript𝑟0𝜇𝑟superscriptsubscript𝑟02𝜇subscript𝜎𝑖subscript𝜎𝑗subscript𝜆𝑖subscript𝜆𝑗\displaystyle=-\frac{\alpha_{s}(\mu)}{4}\frac{1}{6m_{i}m_{j}}\frac{e^{-r/r_{0}% (\mu)}}{rr_{0}^{2}(\mu)}(\vec{\sigma}_{i}\cdot\vec{\sigma}_{j})(\vec{\lambda}_% {i}\cdot\vec{\lambda}_{j})= - divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_μ ) end_ARG start_ARG 4 end_ARG divide start_ARG 1 end_ARG start_ARG 6 italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG divide start_ARG italic_e start_POSTSUPERSCRIPT - italic_r / italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_μ ) end_POSTSUPERSCRIPT end_ARG start_ARG italic_r italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_μ ) end_ARG ( over→ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) ( over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT )
VijT(r)superscriptsubscript𝑉𝑖𝑗𝑇𝑟\displaystyle V_{ij}^{T}(r)italic_V start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_r ) =116αs(μ)mimjSij(λiλj)×\displaystyle=-\frac{1}{16}\frac{\alpha_{s}(\mu)}{m_{i}m_{j}}S_{ij}\,(\vec{% \lambda}_{i}\cdot\vec{\lambda}_{j})\times= - divide start_ARG 1 end_ARG start_ARG 16 end_ARG divide start_ARG italic_α start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_μ ) end_ARG start_ARG italic_m start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_ARG italic_S start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ( over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) ×
×[1r3er/rg(μ)r(1r2+13rg2(μ)+1rrg(μ))]absentdelimited-[]1superscript𝑟3superscript𝑒𝑟subscript𝑟𝑔𝜇𝑟1superscript𝑟213superscriptsubscript𝑟𝑔2𝜇1𝑟subscript𝑟𝑔𝜇\displaystyle\times\bigg{[}\frac{1}{r^{3}}-\frac{e^{-r/r_{g}(\mu)}}{r}\bigg{(}% \frac{1}{r^{2}}+\frac{1}{3r_{g}^{2}(\mu)}+\frac{1}{rr_{g}(\mu)}\bigg{)}\bigg{]}× [ divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG - divide start_ARG italic_e start_POSTSUPERSCRIPT - italic_r / italic_r start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ( italic_μ ) end_POSTSUPERSCRIPT end_ARG start_ARG italic_r end_ARG ( divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + divide start_ARG 1 end_ARG start_ARG 3 italic_r start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_μ ) end_ARG + divide start_ARG 1 end_ARG start_ARG italic_r italic_r start_POSTSUBSCRIPT italic_g end_POSTSUBSCRIPT ( italic_μ ) end_ARG ) ] (5)

where r0(μ)=r^0mn2μsubscript𝑟0𝜇subscript^𝑟0subscript𝑚𝑛2𝜇r_{0}(\mu)=\hat{r}_{0}\frac{m_{n}}{2\mu}italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_μ ) = over^ start_ARG italic_r end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT divide start_ARG italic_m start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_μ end_ARG with μ𝜇\muitalic_μ the reduced mass of the (ij)𝑖𝑗(ij)( italic_i italic_j ) heavy quark pair, λ𝜆\vec{\lambda}over→ start_ARG italic_λ end_ARG are the colour matrices, σ𝜎\vec{\sigma}over→ start_ARG italic_σ end_ARG the spin matrices and Sij=3(σir^)(σjr^)(σiσj)subscript𝑆𝑖𝑗3subscript𝜎𝑖^𝑟subscript𝜎𝑗^𝑟subscript𝜎𝑖subscript𝜎𝑗S_{ij}=3(\vec{\sigma}_{i}\cdot\hat{r})(\vec{\sigma}_{j}\cdot\hat{r})-(\vec{% \sigma}_{i}\cdot\vec{\sigma}_{j})italic_S start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = 3 ( over→ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over^ start_ARG italic_r end_ARG ) ( over→ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ⋅ over^ start_ARG italic_r end_ARG ) - ( over→ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_σ end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) the tensor operator of the (ij)𝑖𝑗(ij)( italic_i italic_j ) pair with r𝑟\vec{r}over→ start_ARG italic_r end_ARG their relative position.

All the parameters of the model are given in Table 1. We have not included the spin-orbit interaction parts coming from the one-gluon exchange and confinement because they should give small contributions in this calculation. For the same reason, the spin-tensor terms are neglected in the calculation of the ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT (ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT) masses, but are included in the ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPTΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT (ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPTΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT) interaction.

Quark masses (MeV) mcsubscript𝑚𝑐m_{c}italic_m start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT 1763
mbsubscript𝑚𝑏m_{b}italic_m start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT 5110
OGE r^0subscript^𝑟0\hat{r}_{0}over^ start_ARG italic_r end_ARG start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT (fm) 0.181
α0subscript𝛼0\alpha_{0}italic_α start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 2.118
Λ0subscriptΛ0\Lambda_{0}roman_Λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT (fm-1) 0.113
μ0subscript𝜇0\mu_{0}italic_μ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT (MeV) 36.976
Confinement acsubscript𝑎𝑐a_{c}italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT (MeV) 507.4
μcsubscript𝜇𝑐\mu_{c}italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT (fm-1) 0.576
ΔΔ\Deltaroman_Δ (MeV) 184.432
Table 1: Parameters for the quark-quark interaction.

II.2 The wave function of the ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT(ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT)

A precise definition of the wave functions of the ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT baryons (henceforth ΩQQQsubscriptΩ𝑄𝑄𝑄\Omega_{QQQ}roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT) is an essential part of the calculation, because it defines the size of the baryon, which is important for the baryon-baryon interaction.

Once we know the quark-quark interaction, the ΩQQQsubscriptΩ𝑄𝑄𝑄\Omega_{QQQ}roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT wave function can be calculated by solving the Schrödinger equation with the Gaussian Expansion Method (GEM) Hiyama et al. (2003). In the GEM framework one makes an expansion in gaussian wave functions but instead of using only one set of Jacobi coordinates, one includes the lowest orbital angular momentum wave functions using the three sets of possible Jacobi coordinates. The reason to use different sets is that lowest angular momentum wave functions in one set generates higher angular momentum wave functions in the other sets, making a very numerically efficient way to include such high angular momentum components.

However the wave function given by GEM would be quite complicate and would make the calculation of the dibaryon interaction slow. Alternatively, for the calculation of the dibaryon interaction (that would be justified later), the following orbital wave function can be used

ϕ(pξ1,pξ2)italic-ϕsubscript𝑝subscript𝜉1subscript𝑝subscript𝜉2\displaystyle\phi(\vec{p}_{\xi_{1}},\vec{p}_{\xi_{2}})italic_ϕ ( over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) =[2b2π]3/4eb2pξ12[3b22π]3/4e3b24pξ12absentsuperscriptdelimited-[]2superscript𝑏2𝜋34superscript𝑒superscript𝑏2superscriptsubscript𝑝subscript𝜉12superscriptdelimited-[]3superscript𝑏22𝜋34superscript𝑒3superscript𝑏24superscriptsubscript𝑝subscript𝜉12\displaystyle=\bigg{[}\frac{2b^{2}}{\pi}\bigg{]}^{3/4}e^{-b^{2}p_{\xi_{1}}^{2}% }\bigg{[}\frac{3b^{2}}{2\pi}\bigg{]}^{3/4}e^{-\frac{3b^{2}}{4}p_{\xi_{1}}^{2}}= [ divide start_ARG 2 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_π end_ARG ] start_POSTSUPERSCRIPT 3 / 4 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT [ divide start_ARG 3 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_π end_ARG ] start_POSTSUPERSCRIPT 3 / 4 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - divide start_ARG 3 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 end_ARG italic_p start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT (6)

where pξisubscript𝑝subscript𝜉𝑖p_{\xi_{i}}italic_p start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT are the Jacobi coordinates defined as

pξ1subscript𝑝subscript𝜉1\displaystyle\vec{p}_{\xi_{1}}over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT =12(p1p2)absent12subscript𝑝1subscript𝑝2\displaystyle=\frac{1}{2}(\vec{p}_{1}-\vec{p}_{2})= divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT )
pξ2subscript𝑝subscript𝜉2\displaystyle\vec{p}_{\xi_{2}}over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT =23p313(p1+p2)absent23subscript𝑝313subscript𝑝1subscript𝑝2\displaystyle=\frac{2}{3}\vec{p}_{3}-\frac{1}{3}(\vec{p}_{1}+\vec{p}_{2})= divide start_ARG 2 end_ARG start_ARG 3 end_ARG over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT - divide start_ARG 1 end_ARG start_ARG 3 end_ARG ( over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) (7)

In the notation we use for the baryon calculation this corresponds to mode 3 of the GEM basis using only one gaussian with angular momentum zero and the parameters ν=14b2𝜈14superscript𝑏2\nu=\frac{1}{4b^{2}}italic_ν = divide start_ARG 1 end_ARG start_ARG 4 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and λ=13b2𝜆13superscript𝑏2\lambda=\frac{1}{3b^{2}}italic_λ = divide start_ARG 1 end_ARG start_ARG 3 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG. Notice that fixing the relation between the parameters of the gaussians ν𝜈\nuitalic_ν and λ𝜆\lambdaitalic_λ to these values (ν=34λ𝜈34𝜆\nu=\frac{3}{4}\lambdaitalic_ν = divide start_ARG 3 end_ARG start_ARG 4 end_ARG italic_λ), the orbital wave functions is totally symetric, necessary to get a totally antisymetric wave function for the baryon of lowest energy. The spin wave function has to be also symmetric and implies S=32𝑆32S=\frac{3}{2}italic_S = divide start_ARG 3 end_ARG start_ARG 2 end_ARG and the color wave function will be a color singlet.

So our wave function for the baryon is

ψBsubscript𝜓𝐵\displaystyle\psi_{B}italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT =ϕ(pξ1,pξ2)χBξc[13]absentitalic-ϕsubscript𝑝subscript𝜉1subscript𝑝subscript𝜉2subscript𝜒𝐵subscript𝜉𝑐delimited-[]superscript13\displaystyle=\phi(\vec{p}_{\xi_{1}},\vec{p}_{\xi_{2}})\chi_{B}\xi_{c}[1^{3}]= italic_ϕ ( over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , over→ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ) italic_χ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT [ 1 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ] (8)

with χB=((1212)112)32)\chi_{B}=((\frac{1}{2}\frac{1}{2})1\frac{1}{2})\frac{3}{2})italic_χ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT = ( ( divide start_ARG 1 end_ARG start_ARG 2 end_ARG divide start_ARG 1 end_ARG start_ARG 2 end_ARG ) ⁤ 1 divide start_ARG 1 end_ARG start_ARG 2 end_ARG ) divide start_ARG 3 end_ARG start_ARG 2 end_ARG ) the spin wave function and ξc[13]subscript𝜉𝑐delimited-[]superscript13\xi_{c}[1^{3}]italic_ξ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT [ 1 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ] a singlet color wave function.

Using the wave function of Eq. (6), the kinetic energy is given by

T𝑇\displaystyle Titalic_T =ψB|pξ12m+3pξ224m|ψB=32mb2absentquantum-operator-productsubscript𝜓𝐵superscriptsubscript𝑝subscript𝜉12𝑚3superscriptsubscript𝑝subscript𝜉224𝑚subscript𝜓𝐵32𝑚superscript𝑏2\displaystyle=\langle\psi_{B}|\frac{p_{\xi_{1}}^{2}}{m}+\frac{3p_{\xi_{2}}^{2}% }{4m}|\psi_{B}\rangle=\frac{3}{2mb^{2}}= ⟨ italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT | divide start_ARG italic_p start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_m end_ARG + divide start_ARG 3 italic_p start_POSTSUBSCRIPT italic_ξ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_m end_ARG | italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ⟩ = divide start_ARG 3 end_ARG start_ARG 2 italic_m italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (9)

For the interaction energy we can evaluate ψB|V12|ψBquantum-operator-productsubscript𝜓𝐵subscript𝑉12subscript𝜓𝐵\langle\psi_{B}|V_{12}|\psi_{B}\rangle⟨ italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT | italic_V start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ⟩ and multiply by 3, since we have 3 interactions between equivalent quarks. It is easier to evaluate it in coordinate space. The wave function in coordinate space is

ϕB(r3,R3)=[12πb2]3/4er324b2[23πb2]3/4eR323b2subscriptitalic-ϕ𝐵subscript𝑟3subscript𝑅3superscriptdelimited-[]12𝜋superscript𝑏234superscript𝑒superscriptsubscript𝑟324superscript𝑏2superscriptdelimited-[]23𝜋superscript𝑏234superscript𝑒superscriptsubscript𝑅323superscript𝑏2\displaystyle\phi_{B}(r_{3},R_{3})=\bigg{[}\frac{1}{2\pi b^{2}}\bigg{]}^{3/4}e% ^{-\frac{r_{3}^{2}}{4b^{2}}}\bigg{[}\frac{2}{3\pi b^{2}}\bigg{]}^{3/4}e^{-% \frac{R_{3}^{2}}{3b^{2}}}italic_ϕ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT , italic_R start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) = [ divide start_ARG 1 end_ARG start_ARG 2 italic_π italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] start_POSTSUPERSCRIPT 3 / 4 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - divide start_ARG italic_r start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 4 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_POSTSUPERSCRIPT [ divide start_ARG 2 end_ARG start_ARG 3 italic_π italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] start_POSTSUPERSCRIPT 3 / 4 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT - divide start_ARG italic_R start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_POSTSUPERSCRIPT (10)

and so

ψB|V12|ψBquantum-operator-productsubscript𝜓𝐵subscript𝑉12subscript𝜓𝐵\displaystyle\langle\psi_{B}|V_{12}|\psi_{B}\rangle⟨ italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT | italic_V start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ⟩ =4π[12πb2]3/20r32𝑑r3er322b2V(r3)absent4𝜋superscriptdelimited-[]12𝜋superscript𝑏232superscriptsubscript0superscriptsubscript𝑟32differential-dsubscript𝑟3superscript𝑒superscriptsubscript𝑟322superscript𝑏2𝑉subscript𝑟3\displaystyle=4\pi\bigg{[}\frac{1}{2\pi b^{2}}\bigg{]}^{3/2}\int_{0}^{\infty}r% _{3}^{2}dr_{3}e^{-\frac{r_{3}^{2}}{2b^{2}}}V(r_{3})= 4 italic_π [ divide start_ARG 1 end_ARG start_ARG 2 italic_π italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ] start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT italic_r start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_r start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT - divide start_ARG italic_r start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_POSTSUPERSCRIPT italic_V ( italic_r start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) (11)

The mean value distance between quarks is given by

rij2delimited-⟨⟩superscriptsubscript𝑟𝑖𝑗2\displaystyle\sqrt{\langle r_{ij}^{2}\rangle}square-root start_ARG ⟨ italic_r start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⟩ end_ARG =3babsent3𝑏\displaystyle=\sqrt{3}b= square-root start_ARG 3 end_ARG italic_b (12)

and the mass is given by

M=3mb+T+3ψB|V12|ψB𝑀3subscript𝑚𝑏𝑇3quantum-operator-productsubscript𝜓𝐵subscript𝑉12subscript𝜓𝐵\displaystyle M=3m_{b}+T+3\langle\psi_{B}|V_{12}|\psi_{B}\rangleitalic_M = 3 italic_m start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + italic_T + 3 ⟨ italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT | italic_V start_POSTSUBSCRIPT 12 end_POSTSUBSCRIPT | italic_ψ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ⟩ (13)

Finally, the value of the b𝑏bitalic_b parameter is obtained by minimizing the mass

Mb=0𝑀𝑏0\displaystyle\frac{\partial M}{\partial b}=0divide start_ARG ∂ italic_M end_ARG start_ARG ∂ italic_b end_ARG = 0 (14)
ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT
M/b𝑀𝑏\partial M/\partial b∂ italic_M / ∂ italic_b GEM M/b𝑀𝑏\partial M/\partial b∂ italic_M / ∂ italic_b GEM
M𝑀Mitalic_M [MeV] 4810.9 4798.6 14413.8 14396.9
rij2delimited-⟨⟩superscriptsubscript𝑟𝑖𝑗2\sqrt{\langle r_{ij}^{2}\rangle}square-root start_ARG ⟨ italic_r start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⟩ end_ARG [fm] 0.4360 0.4432 0.2716 0.2762
Tdelimited-⟨⟩𝑇\langle T\rangle⟨ italic_T ⟩ [MeV] 522.9 522.8 465.0 471.5
Vdelimited-⟨⟩𝑉\langle V\rangle⟨ italic_V ⟩ [MeV] -333.7 -337.7 -460.4 -468.2
Table 2: Parameters for the ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT baryons obtained from the mass minimization procedure (Mb=0𝑀𝑏0\frac{\partial M}{\partial b}=0divide start_ARG ∂ italic_M end_ARG start_ARG ∂ italic_b end_ARG = 0) and the Gaussian Expansion Method.

Although the GEM method provides a more complete description of the wave function as mentioned before, the calculation is simplified if we use the analytical wave function of Eq. (6). In Table 2 we show the results of the ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT and ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT wave functions using the mass minimization procedure and compare with the GEM solution to justify the use of the simple wave function given by Eq. (6). We see that we get a resonable agreement for the sizes and energies in both cases, although the agreement is better in the beauty sector. The minimal values for b𝑏bitalic_b are given by bmin=0.15679subscript𝑏min0.15679b_{\rm min}=0.15679italic_b start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = 0.15679 fm for the ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT and bmin=0.25172subscript𝑏min0.25172b_{\rm min}=0.25172italic_b start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = 0.25172 fm for the ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT.

III The ΩQQQΩQQQsubscriptΩ𝑄𝑄𝑄subscriptΩ𝑄𝑄𝑄\Omega_{QQQ}\Omega_{QQQ}roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT interaction

The system under study has six identical quarks. Then, the baryon-baryon total wave function must be fully antisymmetric. As the wave function of the baryons is already antisymmetric, the antisymmetrizer operator is just given by,

𝒜=19P36𝒜19subscript𝑃36\displaystyle{\mathcal{A}}=1-9P_{36}caligraphic_A = 1 - 9 italic_P start_POSTSUBSCRIPT 36 end_POSTSUBSCRIPT (15)

In order to obtain the effective baryon-baryon interaction from the underlying quark dynamics we use the Resonating Group Method (RGM) Wheeler (1937); Tang et al. (1978). Then, we need to solve the projected Schrödinger equation,

0=(p22μΩΩE)χ(P)+0superscript𝑝22subscript𝜇ΩΩ𝐸𝜒superscript𝑃\displaystyle 0=\bigg{(}\frac{p^{\prime 2}}{2\mu_{\Omega\Omega}}-E\bigg{)}\chi% (\vec{P}^{\prime})+\int0 = ( divide start_ARG italic_p start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 italic_μ start_POSTSUBSCRIPT roman_Ω roman_Ω end_POSTSUBSCRIPT end_ARG - italic_E ) italic_χ ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) + ∫ (RGMVD(P,Pi)+\displaystyle\bigg{(}\,^{\rm RGM}V_{D}(\vec{P}^{\prime},\vec{P}_{i})+( start_POSTSUPERSCRIPT roman_RGM end_POSTSUPERSCRIPT italic_V start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) +
+RGMK(P,Pi))χ(Pi)d3Pi\displaystyle+\,^{\rm RGM}K(\vec{P}^{\prime},\vec{P}_{i})\bigg{)}\chi(\vec{P}_% {i})d^{3}P_{i}+ start_POSTSUPERSCRIPT roman_RGM end_POSTSUPERSCRIPT italic_K ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ) italic_χ ( over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_P start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT (16)

where Psuperscript𝑃\vec{P}^{\prime}over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT (Pisubscript𝑃𝑖\vec{P}_{i}over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT) is the relative ΩQQQΩQQQsubscriptΩ𝑄𝑄𝑄subscriptΩ𝑄𝑄𝑄\Omega_{QQQ}-\Omega_{QQQ}roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT - roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT final (initial) momentum, E=ET2MΩ𝐸subscript𝐸𝑇2subscript𝑀ΩE=E_{T}-2M_{\Omega}italic_E = italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT - 2 italic_M start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT the relative energy of the system with respect to the threshold, VDRGM(P,Pi)superscriptsubscript𝑉𝐷RGMsuperscript𝑃subscript𝑃𝑖{}^{\rm RGM}V_{D}(\vec{P}^{\prime},\vec{P}_{i})start_FLOATSUPERSCRIPT roman_RGM end_FLOATSUPERSCRIPT italic_V start_POSTSUBSCRIPT italic_D end_POSTSUBSCRIPT ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) is the direct kernel and KRGM(P,Pi)superscript𝐾RGMsuperscript𝑃subscript𝑃𝑖{}^{\rm RGM}K(\vec{P}^{\prime},\vec{P}_{i})start_FLOATSUPERSCRIPT roman_RGM end_FLOATSUPERSCRIPT italic_K ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) is the exchange kernel and μΩΩsubscript𝜇ΩΩ\mu_{\Omega\Omega}italic_μ start_POSTSUBSCRIPT roman_Ω roman_Ω end_POSTSUBSCRIPT is the reduced mass of two ΩQQQsubscriptΩ𝑄𝑄𝑄\Omega_{QQQ}roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT baryons.

Here, MΩsubscript𝑀ΩM_{\Omega}italic_M start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT is

MΩsubscript𝑀Ω\displaystyle M_{\Omega}italic_M start_POSTSUBSCRIPT roman_Ω end_POSTSUBSCRIPT =3mb+32mQb2+3Eintabsent3subscript𝑚𝑏32subscript𝑚𝑄superscript𝑏23subscript𝐸𝑖𝑛𝑡\displaystyle=3m_{b}+\frac{3}{2m_{Q}b^{2}}+3E_{int}= 3 italic_m start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + divide start_ARG 3 end_ARG start_ARG 2 italic_m start_POSTSUBSCRIPT italic_Q end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 3 italic_E start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT (17)
Eintsubscript𝐸𝑖𝑛𝑡\displaystyle E_{int}italic_E start_POSTSUBSCRIPT italic_i italic_n italic_t end_POSTSUBSCRIPT =Vij=d3qeq2b22Vij(q)absentdelimited-⟨⟩subscript𝑉𝑖𝑗superscript𝑑3𝑞superscript𝑒superscript𝑞2superscript𝑏22delimited-⟨⟩subscript𝑉𝑖𝑗𝑞\displaystyle=\langle V_{ij}\rangle=\int d^{3}qe^{-\frac{q^{2}b^{2}}{2}}% \langle V_{ij}(q)\rangle= ⟨ italic_V start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ⟩ = ∫ italic_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_q italic_e start_POSTSUPERSCRIPT - divide start_ARG italic_q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_b start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG end_POSTSUPERSCRIPT ⟨ italic_V start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT ( italic_q ) ⟩ (18)

The direct term will be zero in the present model since the color coefficients, (λiλj)subscript𝜆𝑖subscript𝜆𝑗(\vec{\lambda}_{i}\cdot\vec{\lambda}_{j})( over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ), are zero between color singlets.

Then, the full interaction is driven by exchange diagrams, which take into account the quark rearrangement between baryons. The exchange kernel can be written as,

KRGM(P,Pi)=superscript𝐾RGMsuperscript𝑃subscript𝑃𝑖absent\,{}^{\rm RGM}K(\vec{P}^{\prime},\vec{P}_{i})=start_FLOATSUPERSCRIPT roman_RGM end_FLOATSUPERSCRIPT italic_K ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) = TRGM(P,Pi)+RGMVijE(P,Pi)superscriptRGMsuperscript𝑇RGMsuperscript𝑃subscript𝑃𝑖limit-fromsubscript𝑉𝑖𝑗𝐸superscript𝑃subscript𝑃𝑖\,{}^{\rm RGM}T(\vec{P}^{\prime},\vec{P}_{i})+^{\rm RGM}V_{ij\,E}(\vec{P}^{% \prime},\vec{P}_{i})-start_FLOATSUPERSCRIPT roman_RGM end_FLOATSUPERSCRIPT italic_T ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) + start_POSTSUPERSCRIPT roman_RGM end_POSTSUPERSCRIPT italic_V start_POSTSUBSCRIPT italic_i italic_j italic_E end_POSTSUBSCRIPT ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) -
ETRGMN(P,Pi)superscriptsubscript𝐸𝑇RGM𝑁superscript𝑃subscript𝑃𝑖\displaystyle-E_{T}\,^{\rm RGM}N(\vec{P}^{\prime},\vec{P}_{i})- italic_E start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_RGM end_POSTSUPERSCRIPT italic_N ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) (19)

where TRGM(P,Pi)superscript𝑇RGMsuperscript𝑃subscript𝑃𝑖{}^{\rm RGM}T(\vec{P}^{\prime},\vec{P}_{i})start_FLOATSUPERSCRIPT roman_RGM end_FLOATSUPERSCRIPT italic_T ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) is the exchange kinetic term, NRGM(P,Pi)superscript𝑁RGMsuperscript𝑃subscript𝑃𝑖{}^{\rm RGM}N(\vec{P}^{\prime},\vec{P}_{i})start_FLOATSUPERSCRIPT roman_RGM end_FLOATSUPERSCRIPT italic_N ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) is a normalization term and VijERGM(P,Pi)superscriptsubscript𝑉𝑖𝑗𝐸RGMsuperscript𝑃subscript𝑃𝑖{}^{\rm RGM}V_{ij\,E}(\vec{P}^{\prime},\vec{P}_{i})start_FLOATSUPERSCRIPT roman_RGM end_FLOATSUPERSCRIPT italic_V start_POSTSUBSCRIPT italic_i italic_j italic_E end_POSTSUBSCRIPT ( over→ start_ARG italic_P end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over→ start_ARG italic_P end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) is the exchange potential (for explicit expressions see, e.g., Refs. Ortega et al. (2022); Fernández et al. (2020)).

IV Results

Let us first study the ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT system. One of the states in S𝑆Sitalic_S wave is the JP=0+superscript𝐽𝑃superscript0J^{P}=0^{+}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT, which corresponds to the S01superscriptsubscript𝑆01{}^{1}S_{0}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and D05superscriptsubscript𝐷05{}^{5}D_{0}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT partial waves. As in the case of the deuteron, S𝑆Sitalic_S and D𝐷Ditalic_D waves are mixed. We first calculate the binding energy considering the parameter b𝑏bitalic_b and the reduced mass given by the minimization procedure. Without tensor interactions they are decoupled and only a bound state appears in the S01superscriptsubscript𝑆01{}^{1}S_{0}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT partial wave. The binding energy of this state is E=1.9859𝐸1.9859E=-1.9859italic_E = - 1.9859 MeV. The D05superscriptsubscript𝐷05{}^{5}D_{0}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT partial wave is not bound. If we include the tensor interaction of OGE the partial waves are coupled and the binding energy increases very slightly to E=1.9876𝐸1.9876E=-1.9876italic_E = - 1.9876 MeV. The probability of the D𝐷Ditalic_D wave is only 6.6104%6.6percentsuperscript1046.6\cdot 10^{-4}\%6.6 ⋅ 10 start_POSTSUPERSCRIPT - 4 end_POSTSUPERSCRIPT %. As in the deuteron, the binding energy has a sizeable cancellation between the kinetic and interaction parts. The mean value of the kinetic energy is T=11.2delimited-⟨⟩𝑇11.2\langle T\rangle=11.2⟨ italic_T ⟩ = 11.2 MeV, while for the interaction we have V=13.2delimited-⟨⟩𝑉13.2\langle V\rangle=-13.2⟨ italic_V ⟩ = - 13.2 MeV. The confinement interaction dominates and gives the needed attraction to bind the system. If we exclude the OGE we get E=7.3698𝐸7.3698E=-7.3698italic_E = - 7.3698 MeV with T=20.9delimited-⟨⟩𝑇20.9\langle T\rangle=20.9⟨ italic_T ⟩ = 20.9 MeV and V=28.3delimited-⟨⟩𝑉28.3\langle V\rangle=-28.3⟨ italic_V ⟩ = - 28.3 MeV.

The potential for the S01superscriptsubscript𝑆01{}^{1}S_{0}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT partial wave is given in Fig. 1. The relative wave functions are shown in Fig. 2.

Refer to caption
Figure 1: The potential V(p,p)𝑉superscript𝑝𝑝V(p^{\prime},p)italic_V ( italic_p start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_p ) in the S01superscriptsubscript𝑆01{}^{1}S_{0}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT partial wave for the ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT interaction.
Refer to caption
Refer to caption
Figure 2: The relative wave functions in the S01superscriptsubscript𝑆01{}^{1}S_{0}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and D05superscriptsubscript𝐷05{}^{5}D_{0}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT partial waves for the ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT dibaryon.

If we consider b=rij23𝑏superscriptsubscript𝑟𝑖𝑗23b=\frac{\sqrt{r_{ij}^{2}}}{\sqrt{3}}italic_b = divide start_ARG square-root start_ARG italic_r start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG square-root start_ARG 3 end_ARG end_ARG and the reduced mass given by the GEM calculation we get a binding energy in the coupled case of E=1.81𝐸1.81E=-1.81italic_E = - 1.81 MeV. With b𝑏bitalic_b given by the minimization procedure and the reduced mass is given by GEM we get =1.9754absent1.9754=-1.9754= - 1.9754 MeV. The effect of the different reduced mass is very small and dominates the effect of the different b𝑏bitalic_b parameters. In principle with only one gaussian one should use the value given by the minimization procedure, but this gives us a feeling of the uncertainty due to the simplification of the wave function. Although the binding energy varies a little bit, in both cases the system is bounded.

Another possible state is the JP=2+superscript𝐽𝑃superscript2J^{P}=2^{+}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT which includes the S25superscriptsubscript𝑆25{}^{5}S_{2}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, D21superscriptsubscript𝐷21{}^{1}D_{2}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, D25superscriptsubscript𝐷25{}^{5}D_{2}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and G25superscriptsubscript𝐺25{}^{5}G_{2}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_G start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT partial waves. None of them are bound. One could expect a bound state for the S25superscriptsubscript𝑆25{}^{5}S_{2}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT partial wave, but in this case one can see that the potential coming from the λiλjsubscript𝜆𝑖subscript𝜆𝑗\vec{\lambda}_{i}\cdot\vec{\lambda}_{j}over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ⋅ over→ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT have opposite sign for S=2𝑆2S=2italic_S = 2 with respect to S=0𝑆0S=0italic_S = 0. So if we have attraction for the S=0𝑆0S=0italic_S = 0, this implies repulsion for S=2𝑆2S=2italic_S = 2. Higher partial waves are more difficult to bind.

Antisymmetry implies L+S=𝐿𝑆absentL+S=italic_L + italic_S = even and parity is given by P=(1)L𝑃superscript1𝐿P=(-1)^{L}italic_P = ( - 1 ) start_POSTSUPERSCRIPT italic_L end_POSTSUPERSCRIPT. So for P=+𝑃P=+italic_P = +, the spin S𝑆Sitalic_S has to be even. This means that 1+superscript11^{+}1 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and 3+superscript33^{+}3 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT can be only in D𝐷Ditalic_D or G𝐺Gitalic_G waves, which will be difficult to bind as it was seen for the 0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and 2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT states. In more detail:

  • We start with the 1+superscript11^{+}1 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT state and include D15superscriptsubscript𝐷15{}^{5}D_{1}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT which is the only partial wave. It should be the same as the D25superscriptsubscript𝐷25{}^{5}D_{2}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT partial wave with the exception of the contribution of the OGE tensor interaction. It does not bind.

  • For the 3+superscript33^{+}3 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT state we have the D35superscriptsubscript𝐷35{}^{5}D_{3}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT and G35superscriptsubscript𝐺35{}^{5}G_{3}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_G start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT partial waves and they do not bind.

We give in Fig. 3 the Fredholm determinant for the 4 different J+superscript𝐽J^{+}italic_J start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT quantum numbers, where we can see that only the 0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT channel binds.

Refer to caption
Figure 3: The Fredholm determinant of the ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT system for the 0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT, 1+superscript11^{+}1 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT, 2+superscript22^{+}2 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT and 3+superscript33^{+}3 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT channels. Only the 0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT crosses the zero.

Regarding possible P=𝑃P=-italic_P = - states, this would imply odd partial waves and odd total spin. We have analyzed the JP={0,1,2,3}superscript𝐽𝑃superscript0superscript1superscript2superscript3J^{P}=\{0^{-},1^{-},2^{-},3^{-}\}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = { 0 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , 1 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , 2 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , 3 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT }, finding no additional bound states. Again, in Fig. 4 the Fredholm determinant for the 4 different Jsuperscript𝐽J^{-}italic_J start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT quantum numbers is shown, where we can see that no bound state is predicted.

Refer to caption
Figure 4: The Fredholm determinant of the ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT system for the 0superscript00^{-}0 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT, 1superscript11^{-}1 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT, 2superscript22^{-}2 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT and 3superscript33^{-}3 start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT channels. None of them crosses the zero.

Concerning the ΩcccΩcccsubscriptΩ𝑐𝑐𝑐subscriptΩ𝑐𝑐𝑐\Omega_{ccc}\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT system, the situation is similar to the ΩbbbΩbbbsubscriptΩ𝑏𝑏𝑏subscriptΩ𝑏𝑏𝑏\Omega_{bbb}\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT system, and we only find a bound state in the 0+superscript00^{+}0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT channel. The binding energy is E=0.7104𝐸0.7104E=-0.7104italic_E = - 0.7104 MeV with a D𝐷Ditalic_D-state probability of 1.7103%1.7percentsuperscript1031.7\cdot 10^{-3}\%1.7 ⋅ 10 start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT %. The mean values of kinetic and interaction terms are T=7.46delimited-⟨⟩𝑇7.46\langle T\rangle=7.46⟨ italic_T ⟩ = 7.46 MeV and V=8.17delimited-⟨⟩𝑉8.17\langle V\rangle=-8.17⟨ italic_V ⟩ = - 8.17 MeV. In this case we used the b𝑏bitalic_b parameter from the minimization procedure and the reduced mass from the GEM. Using the reduced mass from the minimization parameter the binding energy changes to E=0.7288𝐸0.7288E=-0.7288italic_E = - 0.7288 MeV and both parameters from the GEM to E=0.62𝐸0.62E=-0.62italic_E = - 0.62 MeV.

IV.1 Dependence on the model parameters

We analyze the dependence on the parameters of the model for the JP=0+superscript𝐽𝑃superscript0J^{P}=0^{+}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT state to see in which parameter space region the system will not bind. In all cases we use the minimization procedure to obtain b𝑏bitalic_b and μΩΩsubscript𝜇ΩΩ\mu_{\Omega\Omega}italic_μ start_POSTSUBSCRIPT roman_Ω roman_Ω end_POSTSUBSCRIPT.

The dependence on the quark mass mqsubscript𝑚𝑞m_{q}italic_m start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT is shown in Fig. 5. Notice that some of the parameters of the potential depends on mqsubscript𝑚𝑞m_{q}italic_m start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT since we use scale dependent parameters. We see that the system binds reducing the quark mass up to mq800900similar-tosubscript𝑚𝑞800900m_{q}\sim 800-900italic_m start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT ∼ 800 - 900 MeV.

Refer to caption
Refer to caption
Refer to caption
Figure 5: Quark mass dependence of b𝑏bitalic_b, MQQQsubscript𝑀𝑄𝑄𝑄M_{QQQ}italic_M start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT and E𝐸Eitalic_E. We show by a dot the result of the Chiral Quark Model.

Our model has an effective string tension given by

σ=83acμc=0.1537GeV2𝜎83subscript𝑎𝑐subscript𝜇𝑐0.1537superscriptGeV2\displaystyle\sigma=\frac{8}{3}a_{c}\mu_{c}=0.1537\,{\rm GeV}^{2}italic_σ = divide start_ARG 8 end_ARG start_ARG 3 end_ARG italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = 0.1537 roman_GeV start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (20)

We plot the parameters b𝑏bitalic_b, MQQQsubscript𝑀𝑄𝑄𝑄M_{QQQ}italic_M start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT and E𝐸Eitalic_E as a function of the string tension in Fig. 6. We vary the value of μcsubscript𝜇𝑐\mu_{c}italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT from 0.15 to 0.85 fm-1 and leave acsubscript𝑎𝑐a_{c}italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT unchanged so the saturation energy does not change.

We see that for higher string tension values (our value is lower than some determinations) the binding energy will increase.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 6: Dependence of b𝑏bitalic_b, MQQQsubscript𝑀𝑄𝑄𝑄M_{QQQ}italic_M start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT and E𝐸Eitalic_E on the effective string tension σ=83acμc𝜎83subscript𝑎𝑐subscript𝜇𝑐\sigma=\frac{8}{3}a_{c}\mu_{c}italic_σ = divide start_ARG 8 end_ARG start_ARG 3 end_ARG italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT for mq=mcsubscript𝑚𝑞subscript𝑚𝑐m_{q}=m_{c}italic_m start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT (left) and mq=mbsubscript𝑚𝑞subscript𝑚𝑏m_{q}=m_{b}italic_m start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT (right). We show by a dot the result of the Chiral Quark Model.

Our confinement effective potential is

V(r)σ1eμcrμcsimilar-to𝑉𝑟𝜎1superscript𝑒subscript𝜇𝑐𝑟subscript𝜇𝑐\displaystyle V(r)\sim\sigma\frac{1-e^{-\mu_{c}r}}{\mu_{c}}italic_V ( italic_r ) ∼ italic_σ divide start_ARG 1 - italic_e start_POSTSUPERSCRIPT - italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT italic_r end_POSTSUPERSCRIPT end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG (21)

We vary the value of μcsubscript𝜇𝑐\mu_{c}italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT from 0.15 to 0.85 fm-1 and change acsubscript𝑎𝑐a_{c}italic_a start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT so that σ𝜎\sigmaitalic_σ does not change. This is the same interval we used when we changed the string tension σ𝜎\sigmaitalic_σ. The saturation energy changes as σμc𝜎subscript𝜇𝑐\frac{\sigma}{\mu_{c}}divide start_ARG italic_σ end_ARG start_ARG italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG. Notice that the interaction region is 3bsimilar-toabsent3𝑏\sim\sqrt{3}b∼ square-root start_ARG 3 end_ARG italic_b so if x3bμc1𝑥3𝑏subscript𝜇𝑐much-less-than1x\equiv\sqrt{3}b\mu_{c}\ll 1italic_x ≡ square-root start_ARG 3 end_ARG italic_b italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ≪ 1 the potential in the interacting region is basically linear. In this calculation we got x=0.062𝑥0.062x=0.062italic_x = 0.062 to x=0.39𝑥0.39x=0.39italic_x = 0.39 in the charm sector and x=0.039𝑥0.039x=0.039italic_x = 0.039 to x=0.24𝑥0.24x=0.24italic_x = 0.24 in the bottom sector. For μc0subscript𝜇𝑐0\mu_{c}\to 0italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT → 0 the potential becomes more linear in the interaction region.

The results varying the saturation are shown in Fig. 7. We see that the dependence on the saturation point of the properties of the ΩQQQsubscriptΩ𝑄𝑄𝑄\Omega_{QQQ}roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT, b𝑏bitalic_b and MQQQsubscript𝑀𝑄𝑄𝑄M_{QQQ}italic_M start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT, is smaller than on the string tension σ𝜎\sigmaitalic_σ as one would expect. For the binding energy of the ΩQQQsubscriptΩ𝑄𝑄𝑄\Omega_{QQQ}roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT dibaryon we see also an smaller dependence.

Notice that when μc0subscript𝜇𝑐0\mu_{c}\to 0italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT → 0 the binding energy increases, so a linear confinement potential should give more binding.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 7: Dependence of b𝑏bitalic_b, MQQQsubscript𝑀𝑄𝑄𝑄M_{QQQ}italic_M start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT and E𝐸Eitalic_E on the saturation parameter μcsubscript𝜇𝑐\mu_{c}italic_μ start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT for mq=mcsubscript𝑚𝑞subscript𝑚𝑐m_{q}=m_{c}italic_m start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT (left) and mq=mbsubscript𝑚𝑞subscript𝑚𝑏m_{q}=m_{b}italic_m start_POSTSUBSCRIPT italic_q end_POSTSUBSCRIPT = italic_m start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT (right). We show by a dot the result of the Chiral Quark Model.

Finally we can study the dependence of the binding energy on the size of the baryon. For that we keep all the parameters unchanged and only vary the parameter b𝑏bitalic_b on the RGM calculation. Results are shown in Fig. 8. With bigger sizes we get less binding but it has to be increased much more than the difference between the sizes of the variational and GEM calculation, which shows that using the exact wave function the system will still bind. This argument is more robust for the bottom sector but it should also work in the charm.

The result should be seen as an upper bound of the binding energy, since we are using a variational calculation. Also other channels may be involved, but since we are considering the lower energy channel, including more channels will provide more attraction. We can conclude that the Chiral Quark Model binds the ΩQQQΩQQQsubscriptΩ𝑄𝑄𝑄subscriptΩ𝑄𝑄𝑄\Omega_{QQQ}\Omega_{QQQ}roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT roman_Ω start_POSTSUBSCRIPT italic_Q italic_Q italic_Q end_POSTSUBSCRIPT system in both cases, when Q𝑄Qitalic_Q is a bottom or a charm quark. These molecular states are analogs of two-atom molecules, where the direct interaction is zero for neutral atoms, as in our model for colorless objects.

Refer to caption
Refer to caption
Figure 8: Dependence of E𝐸Eitalic_E on the size of the baryon given by b𝑏bitalic_b. Charm sector on the left and bottom sector on the right. We show by a dot the result of the Chiral Quark Model.

V Summary

In this work we have studied the possible existence of fully-heavy dibaryons in the charm and bottom sectors. The main conclusion we found is that, using a wave function which minimizes the mass of the ΩcccsubscriptΩ𝑐𝑐𝑐\Omega_{ccc}roman_Ω start_POSTSUBSCRIPT italic_c italic_c italic_c end_POSTSUBSCRIPT (ΩbbbsubscriptΩ𝑏𝑏𝑏\Omega_{bbb}roman_Ω start_POSTSUBSCRIPT italic_b italic_b italic_b end_POSTSUBSCRIPT) baryons, the six c𝑐citalic_c quarks or the six b𝑏bitalic_b quarks can form bound states with JP=0+superscript𝐽𝑃superscript0J^{P}=0^{+}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT quantum numbers. The binding energy of the charm dibaryon is Eb=0.71subscript𝐸𝑏0.71E_{b}=-0.71italic_E start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = - 0.71 MeV, while in the bottom case the binding energy is slightly higher, Eb=1.98subscript𝐸𝑏1.98E_{b}=-1.98italic_E start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = - 1.98 MeV, which is reasonable due to the highest mass of the bottom quark. The JP=0+superscript𝐽𝑃superscript0J^{P}=0^{+}italic_J start_POSTSUPERSCRIPT italic_P end_POSTSUPERSCRIPT = 0 start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT state corresponds to the coupling of S01superscriptsubscript𝑆01{}^{1}S_{0}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and D05superscriptsubscript𝐷05{}^{5}D_{0}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT partial waves, but with a very small D05superscriptsubscript𝐷05{}^{5}D_{0}start_FLOATSUPERSCRIPT 5 end_FLOATSUPERSCRIPT italic_D start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT component. No further bound states are found in other partial waves.

Acknowledgements.
This work has been partially funded by EU Horizon 2020 research and innovation program, STRONG-2020 project, under grant agreement no. 824093; Ministerio Español de Ciencia e Innovación under grant nos. PID2022-141910NB-I00 and PID2022-140440NB-C22; and Junta de Andalucía under contract Nos. PAIDI FQM-370 and PCI+D+i under the title: ”Tecnologías avanzadas para la exploración del universo y sus componentes” (Code AST22-0001).

References