Dibaryons and where to find them

In recent years there has been tremendous progress in the investigation of bound systems of quarks with multiplicities beyond the more usual two- and three-quark systems. Experimental and theoretical progress has been made in the four-, five- and even six-quark sectors. In this paper, we review the possible lightest six-quark states using a simple ansatz based on SU(3) symmetry and evaluate the most promising decay branches. The work will be useful to help focus future experimental searches in this six-quark sector.


Introduction
Theories of the strong interaction indicate there could be more ways to combine multiple quarks into a single object than established experimentally.The recent experimental discoveries of four-and five-quark systems has led to renewed interest in the field of multiquark states -where "multiquark" refers to objects with "non-standard" quark configurations beyond that of baryons and mesons.Unfortunately, for many of these multiquark systems, the internal structure is not currently established.Due to this, CERN has adopted a new scheme in referring to these states.Tetraquarks refer to all objects with four quarks inside regardless if it is a genuine "tetraquark" or a mesonmeson molecule, pentaquarks for five-quark objects (genuine pentaquarks and mesonbaryon molecules) and hexaquarks (genuine hexaquarks and baryon-baryon molecules).We adopt this system in the current paper -but note that some of the earlier literature adopted different systems with hexaquark referring only to the non-molecular states.
The hexaquark family can be further decomposed by quark content since the six-quark system can be arranged in two possible ways with three quarks and three antiquarks (baryonium) and six quarks (dibaryon).In this paper we consider the latter configuration and develop a simple theoretical ansatz so that we can infer the likely decay properties of such "dibaryon" six quark states.
It is interesting to note that hexaquarks were the first multiquark states to be established.The deuteron, a trivial hexaquark composed (predominantly) of a molecular state of a proton and neutron, was discovered in 1931.The earliest paper on non-trivial hexaquarks is attributed to Dyson and Xuong [1], who predicted the existence of six non-strange hexaquarks based on the SU(6) model, just half a year after the discovery of quarks [2].All of the states predicted by Dyson and Xuong have been the subject of experimental searches in the subsequent years.In Table .1 below we summarise the experimental signals for the states of different isospin I and spin J, along with their associated SU(3) multiplet.[12] As shown in the table, all states in this u,d quark sector have given experimental signatures with properties broadly consistent with the predicted properties.Signatures for some of the members have only been obtained very recently.
Dibaryons, like any other particles, appear in SU(3) multiplets.Since SU(3) symmetry works reasonably well, one can use this symmetry to provide expectations of the properties of the family of particles in the multiplet (spin, isospin, mass).
In their pioneering paper, Dyson and Xuong [1] focused on the u,d quark sector and did not consider the properties of strange quark containing members of dibaryon multiplets or their decay branches.However, modern experimental facilities now offer the prospect of experimental study of dibaryons in the strange quark sector.In this paper, we will review what is known about these mysterious strange quark containing particles, point to the most promising experimental channels and, where feasible, estimate their decay branches.
The paper is structured as follows.In Section 2 we review the current experimental evidence and theoretical interpretation of the deuteron anti-decuplet.In Section 3 we cover the complimentary NN-FSI 27-plet.Section 4 addresses the anti-decuplet based on the d * (2380) hexaquark.Following a review we present our new theoretical estimate based on SU(3) symmetry considerations for the decay branches of the strange quark containing partners in this anti-decuplet under assumptions of pure hexaquark and pure molecular natures of the d * multiplet members.We than use the same theoretical approach to investigate decay properties of the ∆∆ 28-plet (Section 5), N ∆ 27-plet (Section 6) and N ∆ 35-plet (Section 7 ).

Deuteron antidecuplet
The deuteron SU(3) multiplet is probably the most well-studied hexaquark SU(3) family.We know that the deuteron, and hence all the other objects in this multiplet, are predominantly baryon-baryon molecules.The estimated genuine hexaquark |6q > component of the deuteron wavefunction is predicted to be only around 0.15% [13].Fig. 1 shows the members of the deuteron multiplet in standard convention (x-axis -3rd projection of isospin, y-axis -strangeness).For each member of this antidecuplet we also show which 8 ⊕ 8 states can contribute.The SU(3) decomposition of this multiplet Figure 1: Deuteron multiplet is also summarised in Table 2.
Out of all multiplet members, only the deuteron appears to be bound.The Λ − N interaction is not attractive enough to make a bound state e.g. the hyperdeuteron appears to be unbound.‡ The ΛN state therefore presents as a virtual state or a Final-State Interaction(FSI) [14].The Σ−N interaction is suggested to be less attractive than Λ−N [14].Due to the large SU(3) mass splitting in this decuplet (M Σ −M Λ ∼ 77 MeV), compared to the MeV-keV range of bound/virtual states, it is difficult to treat these states in an SU(3) limit.Aside from the deuteron, none of the states from the multiplet form a bound state, so all of them can be observed only in elastic/quasi-elastic reactions or via the correlation observations in heavy ion collisions.There were several recent publications from STAR and ALICE about extraction of the ΛN [15], ΣN [15] and ΞN [16] scattering lengths.Due to very unfavourable decay branches in the Σ family, it is extremely difficult to access any correlation distributions where the Σ baryon is involved.However, one may get some information in ΞN → ΣΣ or ΞN → ΛΣ from KLF or JPARC facilities from tertiary Ξ beams.

NN-FSI 27-plet
The 27-plet 8 ⊕ 8 of spin zero states is shown schematiclaly in Fig. 2 and Table 3 details the quark compositions.A strong attraction in the spin-zero pp-system was probably first realised by Migdal [5] and Watson [4] in their description of NN-FSI (Final State Interactions).Indeed it appears the absence of the tensor force in the NN 1 S 0 channel implies the "demon deuteron" would not form, as this system is unbound by 66 keV.Thus this multiplet is expected to be comprised of unbound molecular-like virtual states since the strangeness zero members are established to be unbound.As is the case for the deuteron multiplet(Section 2), this SU(3) family suffers from large SU(3) mass splitting.All the members can be observed only in elastic/quasi-elastic reactions and/or correlations in heavy-ion collisions.
Although this multiplet in isolation would not be expected to produce bound states the presence of overlapping members with the same quark content but different isospin raises the question about mixing.Effects established in the meson sector, such as ρ − ω − ϕ mixing in the case of meson nonet, could in principle have corollaries here.Several members of this 27-plet can be mixed, but the most famous case proposed is the so-called H-dibaryon, a potential Λ − Λ deeply bound state [17].Experimentally this is ‡ With additional particles the nuclear binding would be expected to increase -as a result hypertritium does becomes bound but by only 150keV .
now shown to be unbound [18,19,20,21].Recent advances in Lattice-QCD calculations for this sector indicate it to be located in the vicinity of the NΞ-threshold, [22].In a simplified picture of this multiplet, we have 3 states with isospins 2,1,0 (plus a SU(3) singlet isosingle H-dibaryon state).These states can mix, but only in an isospin violating way, analagous to Λ − Σ mixing in the baryon octet or ρ − ω mixing in meson nonet.The redistribution of wave function components, similar to the ω − ϕ vs ϕ 0 − ϕ 8 in the meson nonet, is not allowed here since all states have different isospin(besides 27 belongs to a singlet state, but it has exactly the same quantum numbers, including isospin, as the central state of the 27-plet.One could therefore expect some mixing between these states.However, since no H-dibaryon was found in the vicinity of the ΛΛ threshold, the I = 0 state from the 27-plet should be even less bound and should have much weaker experimental signatures.
Since, for an experimental observation of this multiplet we are limited to elastic scattering and heavy-ion correlations, there are very few states which we can currently access.An elastic scattering with either direct beams or rescattering of secondary beams within target material [23] can give us access to NN, ΛN, ΣN, ΞN scattering data.
The latter reaction can be performed at the recently proposed KLF facility [24] with K L N → K + Ξ as a first step and ΞN → X as a second step reaction within a large target volume (or with similar methodology at JPARC).
For the heavy-ion correlation searches, all events with neutral final states need to be excluded from realisable experimental study, which essentially removes all channels with Σ baryons.The NN, ΛN, ΞN and ΛΛ channels were already investigated, see Refs [15,16,20,21].The only other channels pending analysis for this multiplet is ΞΞ (the ΞΛ ALICE analysis was recently published [25], but it suffers from low statistics).We do not expect any breakthrough here -with expectation of loosely unbound states.However, these data would help to fix several ChPT coupling constants in the baryonbaryon sector and also improve our understanding of the nuclear equation of state for a system with strangeness, essential for neutron star physics.
The search for strange SU(3) partners of the d * (2380) is a natural next step for hexaquark studies.Characterisation of additional members of the antidecuplet would provide valuable insights into the underlying physics of hexaquark systems.All previous multiplets mentioned above have at best a very low binding suggesting its molecular nature.In contrast the d * (2380) has a very large binding energy, which unavoidably means it should have a large |6q > component even if one considers it from a molecular picture.Moreover, the majority of theoretical papers in different ansatz suggest that the d * (2380) might indeed be a genuine hexaquark-dominated object [40,41,42].Its evidence in photoproduction (e.g.relative contribution of multipoles) also favours this possibility [35,32].
However, as there is an abundance of states and associated decay channels from the many states in the d * (2380) antidecuplet, some prioritisation of initial searches is beneficial.In this paper, we develop a simple model to identify the most appropriate final states for each member of the antidecuplet, including estimates of the partial decay widths.Due to conservation laws, the production of any of the multiplet members is expected to proceed via associated many-body reactions at rather high energies.Deriving accurate predictions of the production cross-sections are therefore challenging.
In their absence, an experimental strategy to search for such states is to look for structure in the cross sections consistent with the predicted locations, widths and branching ratios derived from theoretical models.In the following sections we present the first such calculations to guide future measurement strategies.The properties of the decouplet are inferred under assumption of a hexaquark (Sections 4.1.1)and a molecular state (Sections 4.1.2).

d * (2380) multiplet masses
From the unitary group theory of the strong interaction for the light quark (u,d,s) sector any strongly interacting particle, such as the d * (2380), should be part of a SU(3) multiplet.For the case of the d * (2380) it would be expected to be a member of an antidecuplet, Fig. 3.The spectroscopic study of the other multiplet members would provide important new constraints on the d * internal structure, complementary to that potentially achievable in Form-Factor studies.In a simple molecular picture, the d * SU(3) multiplet would derive from the coupling of two baryon decuplet members bound by long-range pion exchange: corresponding to ∆∆ for the d * (2380) to ∆Ω for the d sss .However, since the pion does not have a coupling to strange quarks, there is an expectation for such molecular systems that the binding energy should decrease with increasing strangeness content of the state.Conversely, in a genuine hexaquark (non-molecular) picture the binding energy should increase with increasing strangeness content as the presence of heavier s-quarks would imply stronger binding.In these two cases, a rough evaluation can be made for the masses of the d * (2380) multiplet.3) multiplet [2].In this approach the masses originate from two factors -the masses of the constituent quarks (M q ) and a contribution from hyperfine splitting ( K Mq 1 Mq 2 ).In this approach the mass difference between the nucleon and the ∆ (which have the same quark content) can be reproduced: Here M q = 363M eV is the masses of light, unflavored constituent quarks and the K is the parameter responsible for hyperfine splitting K M 2 q = 50M eV .One can immediately see that if we substitute the light quark mass in the hyperfine term with the heavier strange quark mass, the splitting gets smaller and the baryon decuplet members get lighter.A hexaquark system has a lot more permutations between quarks, which is why the splitting term is essential.Using this ansatz we have combined the evaluation of the d * (2380) multiplet masses in Table .4. For the entries in the Table.4, note that M n = 363M eV is the mass of a light constituent = 50M eV is the splitting parameter and B h ∼ −550M eV is the hexaquark binding which reproduces the observed d * (2380) mass.This latter parameter is taken to be the same for all members of the multiplet (similar to a bag constant of Ref [43]).
Note that under these assumptions the mass of the d sss is 236 MeV lower than a ∆Ω threshold and only 160 MeV higher than the ΣΞ threshold.This mass difference with the ΣΞ threshold is not large, especially considering the d sss → ΣΞ decay should proceed via a D-wave in the ΣΞ-system.
Although this is a simplified picture of the antidecuplet, its flexibility does enable lower limits for multiplet masses to be evaluated.We note there are a large variety of estimates with different theoretical ansatz presented in Refs.[44,45,46].4.1.2.d * (2380) multiplet as molecules For the case of a molecular state, the dependence of the mass of the states with strangeness content should be very different than for a genuine hexaquark.Due to the larger spatial extent of molecular states, one can assume that the binding is mainly driven by long-range pion exchange.In the SU(6) approximation, the pion couples only to light quarks, so the strength of the interaction ∆ : Σ * : Ξ * : Ω should scale as 3 : 2 : 1 : 0, see Ref. [38].Since the interaction potential (V ) is proportional to the product of coupling constants and the binding energy B ∼ M Red • V 2 , where M Red is the reduced mass [47].One can therefore estimate binding energies for various states under these assumptions.The ∆Ω state is predicted have a binding energy B = 0 M eV .However, since d sss has not only ∆Ω but also a Ξ * Σ * component, we expect some binding also in this case.The results of these simplified calculations are summarised in Table 5.Note that M Red is the reduced mass and f is an effective meson baryon coupling constant fixed to reproduce the d * (2380) mass.
The observed mass of each member of the antidecuplet would be expected to lie between our estimations for a pure "genuine hexaquark" and a pure "molecule", Fig. 4.
We note that group theory actually prohibits a 100% genuine hexaquark state [40], with the maximum purity state expected to contain 20% of molecular admixture.More elaborate recent calculations based on hidden colour quark models predict a 30% molecular contribution to the d * [48].
In the following sections, we evaluate the widths and branching ratios as a function of binding energy.As a benchmark case, we take the same binding for all hexaquark members which, in this simplified ansatz, would correspond to nearly equal molecular/hexaquark contributions.M o l e c u l e H e x a q u a r k 8 ⊕ 8 are independent of the decaying particle's hypercharge/strangeness.This enables the total width of the decaying state to be expressed as with Γ 8 corresponding to the partial width for hexaquark decay into 8 ⊕ 8 and Γ 10 for 10 ⊕ 10.In case of 8 ⊕ 8 decay, the final state particles are stable against strong decay so the Γ 8 decay width can be expressed as with where g 8 is the coupling constant of the hexaquark decay into 8 ⊕ 8, p is the momentum of the ejectile in the hexaquark rest frame, L is the angular momentum in the system and F (p) is a Form-Factor, usually introduced to account for potential barriers, similar to Ref. [49,41].For a J p = 3 + particle decaying into two J p = 1/2 + baryons there are two possible scenarios -a 3 D 3 partial wave with L = 2 and 3 G 3 partial wave with L = 4.
From partial wave analysis of the d * (2380) → pn reaction we know that the majority (∼ 90%) of the 8 ⊕ 8 decay proceeds through the 3 D 3 partial wave [11,30,50].In our calculation, we will assume that all 8 ⊕ 8 decays proceed with L = 2.This assumption leads to minor corrections but allows a significant reduction in the required number of coupling constants.Since all the 3 + -hexaquarks lie far above the 8 ⊕ 8 threshold the Γ 8 width is expected to be nearly constant throughout the resonance.The g 8 constant is the same for all members of the antidecuplet.It is fixed to the value extracted from Γ(d * → pn) = 8 MeV.The 10 ⊕ 10 channel is much more challenging since the baryon decuplet contains resonant states with a rather sizeable width and an associated strong energy dependence.The Γ 10 width is expected to vary significantly between different resonances.We have calculated the Γ 10 as follows here F (q 10 ) is a Form-Factor which depends on the relative momentum (q 10 ) between the two decuplet baryons in the hexaquark decay.We follow the prescription of Ref. [49], and parameterise this dependence in a monopole form with a cut-off parameter Λ ¶.D D 1 /D D 2 are the propagators for the decuplet of baryons, with m 1 /m 2 the baryon masses (or correspondingly the invariant mass of the meson-baryon system M BM from the decuplet baryon decay into the octet of mesons and the octet of baryons) and m D being their nominal Breit-Wigner masses.The energy-dependent width of the baryon decuplet decay, Γ D , is parameterised in a standard form taking a P -wave decay resonance with Blatt-Weisskopf barrier factors of R = 6.3 GeV/c and γ = 0.74, 0.28, 0.13 for the ∆, Σ * and Ξ * respectively [51].Γ tot D is the sum of all partial widths.The width of the Ω is considered to be zero since it is stable with respect to strong decays.The Form-Factor parameters are fixed based on the d * → ∆∆ → dππ invariant mass distributions of Ref. [10].The γ 10 is taken as a normalisation factor and is fixed to reproduce the width at zero binding energy, e.g Γ(d 4.2.2.d * (2380) multiplet decays in a molecular picture If one would consider the d * multiplet to be a purely molecular state one can expect some difference compared to the calculations above.The 10 ⊕ 10 decay would stay unchanged since it is a fall-apart decay and it is only driven by available phase space and the energy dependence of the widths.However, to get a 8 ⊕ 8 decay, one needs to get a quark rearrangement, so this decay would be, besides other things, dependent on the wave function overlap.In the calculations above, the 8 ⊕ 8 decay width stayed essentially constant over the large range of binding energies, while in the molecular picture the wave function overlap at B = 0 should be zero, hence the 8 ⊕ 8 decay width should be also zero.To account for this effect we have modified Eq. 2 with here P (B) is a wave function overlap and g8 is a modified g 8 constant to reproduce the d * (2380) width.
We have evaluated P (B) in the following assumptions: i) all particles assumed to be spherical with a distribution, similar to the proton charge distribution ρ(r) = exp(−a•r), with a standard a = 4.27f m −1 (a coordinate space analogue of the proton dipole Form-Factor with a cut-off parameter Λ = 0.72GeV /c).The distance between molecular components was taken as 10) ¶ the appropriateness of the value for this cutoff parameter is explored in Ref [10], where it was shown to provide agreement with measured invariant mass distributions in the region of the so-called ABC-effect.
with M Red being reduced mass and B is the binding energy.The wave function overlap P (B) derived under these assumptions is shown on the figure below, Fig. 5.The g8 is again fixed to reproduce the value extracted from Γ(d * → pn) = 8 MeV.With these modified 8 ⊕ 8 decay width both the branching ratios and decay widths would change.

Results for the d*(2380) decuplet under assumption of a genuine hexaquark
We first explored the validity of the adopted cut-off parameter Λ = 0.16 GeV/c in the model.The Form-Factor of the form Eq. 5. with a cut-off parameter Λ = 0.16 GeV/c in the d * → ∆∆ → dππ was first introduced in a Ref. [10] to explain the so-called ABC-effect, an enhancement in M ππ close to the threshold.Indeed, for the case where the nucleons in the deuteron have a very small relative momentum, there is a correspondence in the relative momentum between the ∆'s and the relative momentum between the pions (and hence with the pion invariant mass).It was later speculated by A. Gal, that reduction of the ∆∆ system size within the d * can lead to a further reduction of the d * width [51].To clarify the situation we have studied the predicted d * → ∆∆ width dependence as a function of the cut of parameter Λ, Fig. 6.The results indicate that the adopted Λ = 0.16 GeV/c can not only reproduce the ABC effect but also gives agreement with the measured Γ(d * → ∆∆) = 62 MeV.As evident in Fig 6 we also reproduce the trend in which a higher value adopted for the cut-off parameter leads to a smaller width.
We, therefore, adopt Λ = 0.16 GeV/c for subsequent calculations.The predicted width for all decuplet members is shown in Fig. 7.In Fig. 8 we show the predicted branching ratios as a function of binding energy.The results are summarised in Tables.7, 8.
One can clearly see that the 8 ⊕ 8 decays are predicted to be increasingly important for the higher strangeness states, while both partial and total widths reduce   substantially.For all hexaquark members only 10 ⊕ 10 with ∆ in a final state are important.For the d ss non-∆ channels (Σ * Σ * → anything) covers only 6%.For the d sss the non-∆ (Ξ * Σ * ) decay has a 2% branching ratio.It is interesting to note that the semi-electromagnetic decay branch d sss → Ω∆ → ΩNγ is predicted to have a higher probability than the purely hadronic d sss → Ξ * Σ * → ΞΣππ decay.Hexaquarks are expected to be produced copiously in heavy ion collisions, however, our estimations indicate the width for all d * multiplet states is rather large.Unfortunately, this indicates their clean identification in the tough background conditions in typical heavy ion collisions may be challenging.The most feasible channel for such studies, having both a large partial width and a convenient isospin 3/2 For all channels the most prominent final states have high particle multiplicity, indicating hermetic detector apparatus and exclusivity conditions would need to be established.The d sss is the only member with a dominant octet-octet decay, however, the necessity of associated kaon production to conserve strangeness in the production mechanism also provides challenges.The requirements may be more easily achieved using beams containing intrinsic strangeness (e.g.Kaon beams as proposed in Ref. [24]).

Results for d*(238) decuplet in a molecular picture
The strong variation of the 8 ⊕ 8 decay in the molecular picture leads to very different behaviour in total width and branching ratios, see Table.9, Fig. 9, 10.One also should keep in mind that for this case we generally would expect much smaller binding energies for all strange members than B = −84 M eV of a d * (2380).While at B = −84 M eV binding the decay branchings Figure 9: Same as Fig. 7, but for a molecular picture with nominal molecular masses shown by the cyan dotted line as in Table .5.
are nearly identical for the two scenarios, in a purely molecular picture with reduced binding the importance of the 8 ⊕ 8 decays is also reduced.For the d sss case with a tiny (B = 2 MeV) binding the 8 ⊕ 8, ΣΞ decay is Γ(ΣΞ) = 1.5 MeV, or only 2% of the total decay width.It is also interesting to note that in a molecular picture, the state with zero strangeness, d * (2380) should have the smallest width.One can see that even the stability of the Ω baryon for the case of d sss state cannot compensate the massive reduction of the binding energy, which leads to the prediction of a much larger width for d sss compared to the d * (2380).
Figure 10: same as Fig. 8, but for a molecular picture with nominal molecular masses shown by cyan dotted line as in Table .5.

∆∆ 28-plet
The ∆∆ 28-plet(Fig.11) couples only to a 10 ⊕ 10, hence it is extremely difficult to access it even in a low strangeness state, see Table 10.Indeed, the zero strangeness state has isospin equal to 3. Due to isospin conservation, it does not couple to NN, hence The simplest way to produce even the lowest lying ∆∆ state is a nucleon-nucleon collision with two associated pions, pp → (ππ)(∆∆).A six-body final state and rather high energy required for this reaction make it extremely challenging to predict the background from conventional reactions, which can also potentially produce structures in the measured cross sections.As a result bump-hunting techniques become extremely unreliable.If the ∆∆ 28-plet would be strongly bound (∼ 80 MeV or higher, as the d * antidecuplet) it members would have an even smaller width than corresponding d * multiplet particles, since the 8 ⊕ 8 decay is not allowed here.Unfortunately, in the only measurement from Ref. [12] , where the authors tried to extract the maximum isospin projection state, pp → π − π − ∆ ++ ∆ ++ , it was demonstrated that this multiplet is either loosely bound or loosely unbound.That unavoidably means that the widths of all states would be similar to the combined width of constituent particles, e.g. in ∆∆ case it would be Γ ∼ 240 MeV.Under such conditions, the only state which can be experimentally accessed is the strong-decay-free ΩΩ dibaryon.There is a hope to get the ΩΩ scattering length from heavy ion collisions.Lattice calculations also support the expectation of loosely bound 28-plet.For example, Ref. [52] predict ΩΩ state to be either bound by less than an MeV or unbound.As a first approximation, the 28-plet decay widths can be evaluated by taking the 10 ⊕ 10 decay width from the d * antidecuplet.The only difference is the binding energy dependence of the width for a 28-plet compared to an antidecuplet would be the size of a cut-off parameter Λ.As one can see in Fig. 6, for a given binding energy the width gets larger with an increase of the cut-off and for a small binding energy, close to B = 0 M eV the effect of cut-off is minimal due to the normalisation condition.Since we expect the 28-plet Λ to be larger than for the antidecuplet, the use of Λ 10 ∼ 0.16 GeV give us an upper limit for an expected width, Fig. 12.The ΩΩ state is stable against hadronic decays if bound so is therefore not shown in Fig. 12.

N∆ 27-plet
This multiplet (Fig. 13) is very interesting.The zero strangeness N∆ state is known since the '60s and is often referred to as an R. Arndt resonance [6,7].The best way to observe it is a πd → pp or pp → πd reaction, where it reveals itself as a perfect loop on the Argand plot in a 1 D 2 partial wave.There are plenty of data collected on this resonance experimentally -including differential cross-sections, singleand double-polarisation observables [7].There are even some data on deuteron tensor polarisation for this reaction [53].While experimental evidence for this state is settled, the theoretical interpretation of these data is challenging.We have a very long list of proposed explanations, including N∆-FSI, extra attraction due to ρ-meson exchange, the effect of box diagrams and it being a genuine dibaryon state.The latest analysis of Hoshizaki tends to explain it with a dibaryon state [54,55].Also, modern Fadeev calculations from Gal&Garcilazo tend to explain this state with N∆ molecule [56].Very recently Niskanen put a paper explaining this state by scattering dynamics in presence of attractive N − ∆ potential [57].So after 50 years of research the field is not settled.It is also unclear if its various descriptions derive from semantics due to the different theoretical approaches or if we are really discussing different phenomena.Regardless, under the assumption of a bound state with 20 MeV binding energy, it is hard to expect anything other than a loosely bound molecules, inline with the prediction of Ref. [56].We therefore take this as an assumption in our theoretical analysis.There are two important questions to be addressed -what are the expected binding energies for the other states and what are their branching ratios?Flavour configurations for all members of multiplet for both 8⊕10 and 8⊕8 decays together with their nominal masses are shown in Tables 11,12, 13.Making parallels to the deuteron antidecuplet, where the first state (deuteron) is bound by 2.2 MeV and the second(NΛ) is unbound by 166keV only, we can expect for the 27-plet with the first state bound by 20 MeV that higher strangeness bound states might also exist.Since the mass splitting is comparable or even larger than expected binding energies, the use of SU(6) formalism from Section 4.1.2is not applicable.However, we can use similar arguments to make qualitative evaluations.One can say that the binding should be proportional to the number of light quarks squared, so if a state with all 6 light quarks has a binding energy of 20 MeV, the state with 3 light and 3 strange quarks should have a binding energy of ∼ 1.3 MeV (Table 14).The latter state was recently explored both theoretically and experimentally with the observation of pΩ state in heavy ion collisions [58] and in lattice QCD calculations [59], which predicted it to be bound by 2.5 MeV (1.5 MeV strong binding plus 1.0MeV electromagnetic binding).Hence we can expect a single strangeness state to be bound by about 10 MeV and a double strangeness state by 3.2 MeV + .
In all states with a ∆ component, the width is largely dominated by the width of the ∆ and resulting states appeared to have similar widths e.g.Γ( 27 d) = 111MeV and Γ( 27 d ss ) = 96 MeV.All of these states show interesting, and rather counterintuitive behaviour.An increase of binding energy first leads to an increase of the width.However, the increase of the binding leads to a decrease of particle separation in the state and therefore an increase of the wave function overlap.This unavoidably means an increase of the phase Space available for the octet-octet decay branch.That is why all members of this multiplet have wider widths than free states.Strangeness -2 state with an isospin I = 2 looks very unfavourable due to the large Ξ∆ branch as can be seen in Fig. 14, 15.The isospin I = 1 state also has a Ξ∆ component, somewhat suppressed by the CG coefficients.Only the isospin I = 0 state does not have ∆ components.This state also has rather convenient 8⊕8 decays, N Ξ and ΛΛ which can be explored experimentally.As discussed earlier for another 27-plet case one may expect some mixing between states.Due to this uncertainty, we do not show any width/branching calculations for the inner multiplet states, where SU(3) mixing can occur, Table 15.
It is very unlikely that strangeness -3 isospin I = 3/2 state will be observed.A tiny binding of about 1 MeV would be really difficult to see in a 30 MeV wide state.The inner I = 1/2 state is a lot more interesting.We do not know the composition of that state due to unknown isospin mixing and enormous symmetry breaking due to the mass difference, however, we do know that the pΩ component is by far the lightest one among all(M (pΩ) = 2611 MeV, M (ΛΞ * ) = 2648 MeV, M (ΞΣ * ) = 2705 MeV,M (ΣΞ * ) = 2721 MeV).Both Lattice QCD [59] and the early  15: same as Fig. 14, but for branching ratios heavy ion correlation function analysis [58] claimed the existence of this state with a binding energy of about 2.5 MeV, where half of the binding originates from the electromagnetic attraction.So this state, if it exists, is a mixture of an atomic and a hadronic molecule.The latest heavy ion results with increased statistics tend to suggest this state is unbound.If bound, this state is stable against 8 ⊕ 10 decays, since Ω is stable against strong decays, however, it is still allowed to decay into 8 ⊕ 8 via the most prominent ΛΞ channel due to both phase space and the large, 9/10, associated CG coefficient.For a ∼ 2.5 MeV binding energy one can expect a measurable Γ(pΩ → ΣΞ) ∼ 4 MeV width, see Fig. 16.This coupled channel effect also needs to be taken into account in the analysis of heavy ion correlation data.The strangeness -4 state is very interesting.According to CG coefficients, it is mainly made of ΣΩ with ΞΞ * occupying only a quarter of the wavefunction.However, the mass of a ΣΩ configuration is nearly 10 MeV higher than the one of ΞΞ * .So, if Figure 16: pΩ → ΣΞ decay width it exists, this state will be bound relative to the ΞΞ * threshold and with a ΣΩ branch being exactly zero.Also due to the very small binding, and hence the large size of this molecule, the wave function overlap will be tiny.That is why the 8 ⊕ 8 decay is strongly suppressed.All these factors would lead to an extremely small width of this state, provided it is bound.It should also be mentioned that an atomic Σ + Ω − state could be possible.It probably will be bound by about 1 MeV and have a width in the order of MeV with two possible decay branches ΞΞ * and ΞΞ.Both branches require quark rearrangements hence both will be small.The ΞΞ has large phase space, but it requires additional spin-flip and angular momentum L = 2.The ΞΞ * happens in S−wave, does not require spin flips, but has tiny available phase space < 10 MeV.One can search this state in ΞΞ pair correlations in heavy ion collisions.The study of this atomic state can give a lot of information about hadronic Σ + Ω − interactions which would appear as a small addition to the larger electromagnetic binding.

N∆ 35-plet
In general, the analogy between N∆ 27-and 35-plets(Fig.17) should be similar to that between the 10 * − and 27− NN multiplets.Since the 27-plet is less bound compared to 10 * −, the 35-plet should be less bound compare to N∆ 27-plet.Accessing the members of these multiplet is challenging -they all have large isospin, they do not couple to 8 ⊕ 8, so production of any of these states require associated particles and all decays are fall-apart many body decays.By analogy to NN-multiplets, one can even expect that all states of this multiplet should be unbound (deuteron is bound, while nn/pp-states of "Demon deuteron" are not).However, recently it was found experimentaly that a spin S = 1, isospin I = 2 N∆ state appears to bound by nearly 20 MeV -only slightly smaller compared to the N∆ state of a 27-plet [8].If true, it implies that the other states of this multiplet may also be bound.That is why we decided to calculate possible decay width of these states as well, see Fig. 18, in a similar fashion as a 27-plet.The flavour configurations for these states are presented in Table.16.In all cases we expect the widths of these states to be close to the nominal width of their constituents, with not much chance for an experiment to reliably detect any of these states.The two exceptions are ΛΩ and ΞΩ, which would be stable against strong decays, if bound.These states can be potentially accessed via heavy-ion correlation functions studies.However, accumulating sufficient statistics for the strangeness -4 or -5 states may be extremely challenging.

Summary
We have developed a theoretical model, employing experimentally constrained parameters, to predict the possible decay branches and partial widths for all members of the d * hexaquark antidecuplet.For all strange-quark containing members of the antidecuplet, the predicted widths are rather large, with the most promising decay channels including the broad ∆ resonance in the final state.We demonstrated that a d * Form-Factor, which was first introduced to explain peculiarities in the d * dππ decay, can also explain the smallness of the d * width, in agreement with the qualitative dimensional arguments of A. Gal.The results of the paper will be an important guide for the ongoing search for the d * anti-decuplet members employing photon-, pion-, and kaoninduced reactions, as well as in high-energy collider experiments.We have also extended our model to systematically study all other light and strange dibaryons which can appear as a ground state of 8 ⊕ 8, 8 ⊕ 10 and 10 ⊕ 10 baryon multiplets configurations.Several interesting states were identified for experimental searches and strategies in search of these states were proposed.

Figure 5 :
Figure 5: The density overlap for two diffuse balls in a molecule as a function of the distance(left) and binding energy(right).

Figure 6 :
Figure 6: d * (2380) width as a function of the Form Factor cut-off parameter Λ (left) and partial width Γ d * →∆∆ as a function of binding energy (right) for various Λ values, Λ = 0.1, 0.16, 0.2, 0.5, 1 GeV/c in rainbow order from violet to red.

Figure 7 :
Figure 7: d * multiplet total width as a function of binding energy (relative to the lightest member of the Decuplet-Decuplet pole) for the d * , d s , d ss , d sss from left to right split into major decay branches (note the log scale).The vertical line (common to all figures) shows the nominal expected mass, obtained under the assumption of the same binding (84 MeV) for all multiplet members(black dashed) and for specific binding (green dotted) as specified in Table.4.
Figure 7: d * multiplet total width as a function of binding energy (relative to the lightest member of the Decuplet-Decuplet pole) for the d * , d s , d ss , d sss from left to right split into major decay branches (note the log scale).The vertical line (common to all figures) shows the nominal expected mass, obtained under the assumption of the same binding (84 MeV) for all multiplet members(black dashed) and for specific binding (green dotted) as specified in Table.4.

28 Figure 12 :
Figure 12: 28-plet total width as a function of binding energy (relative to the Decuplet-Decuplet pole) for the various states with increased strangeness from top left (S = 0) to bottom right (S = −5) split into major decay branches.

Table 5 :
Expected masses of the d * (2380) SU(3) multiplet in pure molecular picture

Table 6 :
Expected decay branches of the d * (2380) SU(3) multiplet * SU(3) antidecuplet, the decay couplings are limited to either octet+octet, 8 ⊕ 8 or decuplet+decuplet, 10 ⊕ 10 baryons.The possible decay branches, deriving from the established 8 ⊕ 8 and 10 ⊕ 10 baryonic members are summarised in Table.6.4.2.1.FormalismIn our calculations, we assume that the Breit-Wigner width of all the resonances is energy-dependent and that the coupling constants for the decay into

Table 7 :
Figure 8: d * multiplet branching ratio as a function of binding energy (relative to Decuplet-Decuplet pole) for the d * , d s , d ss , d sss from left to right.The vertical dashed line (common to all figures) shows the nominal expected mass, obtained under the assumption of the same binding (84 MeV) for all multiplet members(black dashed) and for specific binding (green dotted) as specified in Table.4.d * multiplet width results for the B = −84 MeV.
separation, is the ΞΣ branch in d sss decay.This could potentially be tested using Ξ − Σ correlation functions.

Table 8 :
d * multiplet width results for the expected binding energies from Table.4.

Table 9 :
d * multiplet width results for the expected binding energies from Table5.(molecule)

Table 12 :
Expected 8 ⊕ 8 decay branches of the spin J = 2 27-plet 27 d s state shows unusually small width of Γ( 27 d s ) ∼ 40 MeV.The lightest single strangeness component, NΣ * has a mass which is nearly 30 MeV lighter than the next lightest Λ∆ branch.An additional 10 MeV 27 d s binding means that this state is bound by 40 MeV relative to Λ∆ pole.Such large binding leads to a sizeable reduction of a ∆ A+ the strangeness 4 state is likely to be unbound

Table 13 :
[60]es of 8 ⊕ 10 and 8 ⊕ 8 decay branches of the spin J = 2 27-plet width.On top of which, the Λ∆ branch is also suppressed by Clebsh-Gordan coefficients (see Tab. 11), making 27 d s sufficiently narrow, but still wider than a free Σ * state.The experimental search of this state might be challenging.One could search for the Λ∆ branch, however the necessity of a partial wave analysis of a 4 body reaction(Λ, N, π + an associated particle) makes this channel challenging from a theoretical point of view.An 8⊕8 decay, where the state can reveal itself as a Flatte[60]shape with sub-threshold (relative to the 8 ⊕ 10 pole) peaking or as a cusp-like behaviour looks more appealing theoretically, but very inconvenient from an experimental point of view due to necessity to measure Σ in the final state.

Table 14 :
Expected masses of the spin J = 2 d-27 SU(3) multiplet in pure molecular picture Figure 14: d-27 multiplet total width as a function of binding energy (relative to the lightest member of the Octet-Decuplet pole) for the states with various strangeness increasing from left to right split into major decay branches.The vertical dashed line shows the expected mass as specified in Table.14

Table 15 :
27d multiplet width results for the maximum isospin states and expected binding energies from Table14.