Study of the and dibaryons in constituent quark model
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 ), 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 and dibaryons is investigated in the framework of a constituent quark model that satisfactorily describes the deuteron, the and the interaction. candidates are found in both systems with binding energies of the order of MeV.
pacs:
12.39.Pn, 14.40.Lb, 14.40.RtI 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 , was reported by the WASA-at-COSY Collaboration Adlarson et al. (2011) from the double pionic fusion reaction . This resonance can be described as a nonstrange dibaryon with . 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 and partial waves by one-pion exchange tensor interactions. Similarly, the binding of the 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 is stable against strong interactions. This fact opens the possibility to study systems like or . Moreover, in contrast to the deuteron and the 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 is loosely bound by MeV, while Ref. Mathur et al. (2023) found a very deep state with a binding energy of 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 and , among others. They found a bound state for the system with a binding energy of MeV and another state bound by MeV, contrary to naive espectations. Deng Deng (2023) performed a study of the di-, di- and di- systems using a naive one-gluon exchange quark model and a chiral quark model including and 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- and di- states with binding energies around MeV. Using a different model, namely QCD sum rules, Wang Wang (2022) found for each di- and di- systems two and 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 dibaryons and extrapolated their results to the properties of the and systems. They found no bound states for or 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 and 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 as a dibaryon Garcilazo et al. (1997); Valcarce et al. (2001). Although the binding energy of the 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 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 () baryons and the way we derive the 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
(1) |
where is the dynamical (constituent) quark mass and is the matrix of Goldstone-boson fields that can be expanded as
(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
(3) |
where are the colour matrices, is the gluon field and is the strong coupling constant. The scale dependence of allows a consistent description of light, strange and heavy mesons. Its explicit expression can be found in, e.g., Ref. Vijande et al. (2005),
(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
(5) |
where with the reduced mass of the heavy quark pair, are the colour matrices, the spin matrices and the tensor operator of the pair with 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 () masses, but are included in the () interaction.
Quark masses (MeV) | 1763 | |
5110 | ||
OGE | (fm) | 0.181 |
2.118 | ||
(fm-1) | 0.113 | |
(MeV) | 36.976 | |
Confinement | (MeV) | 507.4 |
(fm-1) | 0.576 | |
(MeV) | 184.432 |
II.2 The wave function of the ()
A precise definition of the wave functions of the and baryons (henceforth ) 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 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
(6) |
where are the Jacobi coordinates defined as
(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 and . Notice that fixing the relation between the parameters of the gaussians and to these values (), 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 and the color wave function will be a color singlet.
So our wave function for the baryon is
(8) |
with the spin wave function and a singlet color wave function.
Using the wave function of Eq. (6), the kinetic energy is given by
(9) |
For the interaction energy we can evaluate 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
(10) |
and so
(11) |
The mean value distance between quarks is given by
(12) |
and the mass is given by
(13) |
Finally, the value of the parameter is obtained by minimizing the mass
(14) |
GEM | GEM | |||
---|---|---|---|---|
[MeV] | 4810.9 | 4798.6 | 14413.8 | 14396.9 |
[fm] | 0.4360 | 0.4432 | 0.2716 | 0.2762 |
[MeV] | 522.9 | 522.8 | 465.0 | 471.5 |
[MeV] | -333.7 | -337.7 | -460.4 | -468.2 |
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 and 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 are given by fm for the and fm for the .
III The 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,
(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,
(16) |
where () is the relative final (initial) momentum, the relative energy of the system with respect to the threshold, is the direct kernel and is the exchange kernel and is the reduced mass of two baryons.
Here, is
(17) | ||||
(18) |
The direct term will be zero in the present model since the color coefficients, , 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,
IV Results
Let us first study the system. One of the states in wave is the , which corresponds to the and partial waves. As in the case of the deuteron, and waves are mixed. We first calculate the binding energy considering the parameter and the reduced mass given by the minimization procedure. Without tensor interactions they are decoupled and only a bound state appears in the partial wave. The binding energy of this state is MeV. The 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 MeV. The probability of the wave is only . 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 MeV, while for the interaction we have MeV. The confinement interaction dominates and gives the needed attraction to bind the system. If we exclude the OGE we get MeV with MeV and MeV.
The potential for the partial wave is given in Fig. 1. The relative wave functions are shown in Fig. 2.
If we consider and the reduced mass given by the GEM calculation we get a binding energy in the coupled case of MeV. With given by the minimization procedure and the reduced mass is given by GEM we get MeV. The effect of the different reduced mass is very small and dominates the effect of the different 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 which includes the , , and partial waves. None of them are bound. One could expect a bound state for the partial wave, but in this case one can see that the potential coming from the have opposite sign for with respect to . So if we have attraction for the , this implies repulsion for . Higher partial waves are more difficult to bind.
Antisymmetry implies even and parity is given by . So for , the spin has to be even. This means that and can be only in or waves, which will be difficult to bind as it was seen for the and states. In more detail:
-
•
We start with the state and include which is the only partial wave. It should be the same as the partial wave with the exception of the contribution of the OGE tensor interaction. It does not bind.
-
•
For the state we have the and partial waves and they do not bind.
We give in Fig. 3 the Fredholm determinant for the 4 different quantum numbers, where we can see that only the channel binds.
Regarding possible states, this would imply odd partial waves and odd total spin. We have analyzed the , finding no additional bound states. Again, in Fig. 4 the Fredholm determinant for the 4 different quantum numbers is shown, where we can see that no bound state is predicted.
Concerning the system, the situation is similar to the system, and we only find a bound state in the channel. The binding energy is MeV with a -state probability of . The mean values of kinetic and interaction terms are MeV and MeV. In this case we used the 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 MeV and both parameters from the GEM to MeV.
IV.1 Dependence on the model parameters
We analyze the dependence on the parameters of the model for the state to see in which parameter space region the system will not bind. In all cases we use the minimization procedure to obtain and .
The dependence on the quark mass is shown in Fig. 5. Notice that some of the parameters of the potential depends on since we use scale dependent parameters. We see that the system binds reducing the quark mass up to MeV.
Our model has an effective string tension given by
(20) |
We plot the parameters , and as a function of the string tension in Fig. 6. We vary the value of from 0.15 to 0.85 fm-1 and leave 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.
Our confinement effective potential is
(21) |
We vary the value of from 0.15 to 0.85 fm-1 and change so that does not change. This is the same interval we used when we changed the string tension . The saturation energy changes as . Notice that the interaction region is so if the potential in the interacting region is basically linear. In this calculation we got to in the charm sector and to in the bottom sector. For 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 , and , is smaller than on the string tension as one would expect. For the binding energy of the dibaryon we see also an smaller dependence.
Notice that when the binding energy increases, so a linear confinement potential should give more binding.
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 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 system in both cases, when 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.
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 () baryons, the six quarks or the six quarks can form bound states with quantum numbers. The binding energy of the charm dibaryon is MeV, while in the bottom case the binding energy is slightly higher, MeV, which is reasonable due to the highest mass of the bottom quark. The state corresponds to the coupling of and partial waves, but with a very small 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
- Adlarson et al. (2011) P. Adlarson et al. (WASA-at-COSY), Phys. Rev. Lett. 106, 242302 (2011), arXiv:1104.0123 [nucl-ex] .
- Goldman et al. (1989) J. T. Goldman, K. Maltman, G. J. Stephenson, Jr., K. E. Schmidt, and F. Wang, Phys. Rev. C 39, 1889 (1989).
- Clement (2017) H. Clement, Prog. Part. Nucl. Phys. 93, 195 (2017), arXiv:1610.05591 [nucl-ex] .
- Huang et al. (2014) H. Huang, J. Ping, and F. Wang, Phys. Rev. C 89, 034001 (2014).
- Bjorken (1985) J. D. Bjorken, AIP Conf. Proc. 132, 390 (1985).
- Lyu et al. (2021) Y. Lyu, H. Tong, T. Sugiura, S. Aoki, T. Doi, T. Hatsuda, J. Meng, and T. Miyamoto, Physical Review Letters 127 (2021), 10.1103/physrevlett.127.072003.
- Mathur et al. (2023) N. Mathur, M. Padmanath, and D. Chakraborty, Physical Review Letters 130 (2023), 10.1103/physrevlett.130.111901.
- Huang et al. (2022) H. Huang, J. Ping, X. Zhu, and F. Wang, Eur. Phys. J. C 82, 805 (2022), arXiv:2011.00513 [hep-ph] .
- Deng (2023) C.-R. Deng, Physical Review D 108 (2023), 10.1103/physrevd.108.054037.
- Wang (2022) Z.-G. Wang, International Journal of Modern Physics A 37 (2022), 10.1142/s0217751x22501664.
- Richard et al. (2020) J.-M. Richard, A. Valcarce, and J. Vijande, Phys. Rev. Lett. 124, 212001 (2020), arXiv:2005.06894 [hep-ph] .
- Alcaraz-Pelegrina and Gordillo (2022) J. Alcaraz-Pelegrina and M. Gordillo, Physical Review D 106 (2022), 10.1103/physrevd.106.114028.
- Lü et al. (2022) Q.-F. Lü, D.-Y. Chen, and Y.-B. Dong, (2022), arXiv:2208.03041 [hep-ph] .
- Weng and Zhu (2024) X.-Z. Weng and S.-L. Zhu, The European Physical Journal C 84 (2024), 10.1140/epjc/s10052-024-12395-4.
- Fernandez et al. (1993) F. Fernandez, A. Valcarce, U. Straub, and A. Faessler, J. Phys. G 19, 2013 (1993).
- Vijande et al. (2005) J. Vijande, F. Fernandez, and A. Valcarce, J. Phys. G31, 481 (2005), arXiv:hep-ph/0411299 [hep-ph] .
- Segovia et al. (2008) J. Segovia, A. M. Yasser, D. R. Entem, and F. Fernández, Phys. Rev. D 78, 114033 (2008).
- Valcarce et al. (1994) A. Valcarce, A. Buchmann, F. Fernández, and A. Faessler, Phys. Rev. C 50, 2246 (1994).
- Juliá-Díaz et al. (2002) B. Juliá-Díaz, D. R. Entem, A. Valcarce, and F. Fernández, Phys. Rev. C 66, 047002 (2002).
- Garcilazo et al. (1997) H. Garcilazo, F. Fernández, A. Valcarce, and R. D. Mota, Phys. Rev. C 56, 84 (1997).
- Valcarce et al. (2001) A. Valcarce, H. Garcilazo, R. D. Mota, and F. Fernandez, J. Phys. G 27, L1 (2001).
- Diakonov (2003) D. Diakonov, Prog. Part. Nucl. Phys. 51, 173 (2003), arXiv:hep-ph/0212026 [hep-ph] .
- Mateu et al. (2019) V. Mateu, P. G. Ortega, D. R. Entem, and F. Fernández, Eur. Phys. J. C 79, 323 (2019), arXiv:1811.01982 [hep-ph] .
- Bali et al. (2005) G. S. Bali, H. Neff, T. Duessel, T. Lippert, and K. Schilling (SESAM), Phys. Rev. D71, 114513 (2005), arXiv:hep-lat/0505012 [hep-lat] .
- Hiyama et al. (2003) E. Hiyama, Y. Kino, and M. Kamimura, Prog. Part. Nucl. Phys. 51, 223 (2003).
- Wheeler (1937) J. A. Wheeler, Phys. Rev. 52, 1083 (1937).
- Tang et al. (1978) Y. C. Tang, M. Lemere, and D. R. Thompson, Phys. Rept. 47, 167 (1978).
- Ortega et al. (2022) P. G. Ortega, D. R. Entem, and F. Fernández, “Quark Models for Baryon-Baryon Interactions,” in Handbook of Nuclear Physics, edited by I. Tanihata, H. Toki, and T. Kajino (2022) pp. 1–26.
- Fernández et al. (2020) F. Fernández, P. G. Ortega, and D. R. Entem, Front. in Phys. 7, 233 (2020).