To bind or not to bind: The dibaryon in light of chiral effective field theory
Abstract
We analyse the quark mass dependence of the binding energy of the dibaryon in the framework of chiral effective field theory. We show that the SU(3) breaking effects induced by the differences of the pertinent twobaryon thresholds (, , ) have a very pronounced impact that need to be incorporated properly in future lattice QCD simulations. We also point out that if the dibaryon is a twobaryon bound state, its dominant component is rather than , which is a consequence of the approximate SU(3) flavor symmetry of the twobaryon interactions.
pacs:
13.75.Ev,12.39.Fe;14.20.PtI Introduction and summary
In 1977 Jaffe predicted a deeply bound 6quark state with , , from the bagmodel, called the dibaryon (1). Subsequently, many experimental searches for the dibaryon were carried out, but so far no convincing signal was found (2). However, recently evidence for a bound dibaryon was claimed based on lattice QCD calculations (3); (4); (5); (6). Extrapolations of the simulations, performed for pion masses MeV, to the physical mass suggest that the dibaryon could be either loosely bound or move into the continuum (7); (8).
In this paper, we analyze various issues related to these lattice results
in the framework of chiral effective field theory (EFT) for the baryonbaryon ()
interaction at leading order (LO) in the Weinberg counting. This scheme
has proven successful for the few data on hyperonnucleon scattering (9)
and also for the bounds that exist on the interactions with
strangeness (10). There is one lowenergy constant (LEC)
in the sector, corresponding to the SU(3) flavorsinglet channel,
that can not be fitted by hyperonnucleon data and can be
finetuned to produce a bound with a given binding energy. This framework
also allows us to study the quark mass dependence
The pertinent results of our study can be summarized as follows:

We have analyzed the effective range expansion for the channel, assuming a loosely bound dibaryon. It shows a very different behaviour to the case of the deuteron in the neutronproton channel. In fact, any attraction supplied by the flavorsinglet channel contributes with a much larger weight to the interaction than to , according to SU(3) flavor symmetry, so that the dibaryon should predominantly be a bound state.

We observe that, for pion masses below 400 MeV, the dependence of the binding energy of the is linearly decreasing with decreasing pion mass, in agreement with the findings in Ref. (7). In particular, a binding energy (BE) adjusted to the value found by NPLQCD (3) at MeV, is reduced by 7 MeV at the physical pion mass. For larger pion masses, this dependence is weakened. Note, however, that for such large pion masses this should only be considered a trend as the chiral EFT is constructed for masses/momenta well below the chiral symmetry breaking scale.

We find a much more drastic effect caused by the SU(3) breaking related to the values of the three thresholds , and . For physical values the BE of the is reduced by as much as 60 MeV as compared to a calculation based on degenerate (i.e. SU(3) symmetric) thresholds. Translating this observation to the situation in the HAL QCD (4) calculation, we see that the bound state has disappeared at the physical point. For the case of the NPLQCD calculation, a resonance in the system might survive.
Our manuscript is organized as follows: In Sec. II, we recall the basic formalism of the interaction in the framework of chiral EFT. Sec. III contains a detailed discussion of the quark mass dependence of the BE of the and the influence of the SU(3) breaking through the various twobaryon thresholds. In Sec. IV we try to make direct contact to the results published by the NPLQCD and HAL QCD collaborations.
Ii The baryonbaryon interaction to leading order
For details on the derivation of the chiral potentials for the strangeness sector at LO using the Weinberg power counting, we refer the reader to Refs. (9); (10), see also Refs. (15); (16); (17). Here, we just briefly summarize the pertinent ingredients of the chiral EFT for interactions.
Channel  Isospin  Isospin  
The LO potential consists of fourbaryon contact terms without derivatives and of onepseudoscalarmeson exchanges. The LO invariant contact terms for the octet interactions that are Hermitian and invariant under Lorentz transformations follow from the Lagrangians
(1) 
Here and denote the Dirac indices of the particles, is the irreducible octet (matrix) representation of , and the are the usual elements of the Clifford algebra (9). As described in Ref. (9), to LO the Lagrangians in Eq. (1) give rise to only six independent lowenergy coefficients (LECs), the s in Eq. (1), due to constraints. They need to be determined by a fit to experimental data. It is convenient to reexpress the potentials in terms of the irreducible representations, see e.g. Refs. (18); (19). Then the contact interaction is given by
(2) 
and the constraints imposed by the assumed symmetry on the interactions in the various channels for the and partial waves can be readily read off from Table 1.
The lowest order invariant pseudoscalarmeson–baryon interaction Lagrangian embodying the appropriate symmetries was also discussed in (9). The invariance under transformations implies specific relations between the various coupling constants, namely
(3) 
Here , where is the nucleon axialvector strength and is the weak pion decay constant. We use the values and MeV. For , the ratio (9), we adopt the SU(6) value: , which is consistent with recent determinations of the axialvector coupling constants (20).
The spinspace part of the LO onepseudoscalarmesonexchange potential is similar to the static onepionexchange potential in chiral EFT for nucleonnucleon interactions, see e.g. (21) (recoil and relativistic corrections give higher order contributions),
(4) 
where is the mass of the exchanged pseudoscalar meson. The transferred and average momentum, and , are defined in terms of the final and initial centerofmass (c.m.) momenta of the baryons, and , as and . In the calculation we use the physical masses of the exchanged pseudoscalar mesons. The explicit breaking reflected in the mass splitting between the pseudoscalar mesons and, in particular, the small mass of the pion relative to the other members of the octet leads to sizeable differences in the range of the interactions in the different channels and, thus, induces an essential dynamical breaking of symmetry in the interactions. The meson was identified with the octet () and its physical mass was used.
The reaction amplitudes are obtained from the solution of a coupledchannels LippmannSchwinger (LS) equation for the interaction potentials:
(5) 
The label indicates the particle channels and the label the partial wave. is the pertinent reduced mass. The onshell momentum in the intermediate state, , is defined by . Relativistic kinematics is used for relating the laboratory energy of the hyperons to the c.m. momentum.
In (9); (10) the LS equation was solved in the particle basis, in order to incorporate the correct physical thresholds. Since here we are only interested in the dibaryon we work in the isospin basis. Then for , , and we have to consider the three coupled channels , and . We use the following (isospin averaged) masses MeV, MeV, MeV, and MeV so that the , , and thresholds are at 2231.2, 2257.7, and 2385.0 MeV, respectively. Furthermore, the potentials in the LS equation are cut off with a regulator function, , in order to remove highenergy components of the baryon and pseudoscalar meson fields (22). We consider cutoff values in the range 550, …, 700 MeV, similar to what was used for chiral potentials (22).
The imposed flavor symmetry implies that only five of the six LECs contribute to the interaction, namely , , , , and , cf. Table 1. These five contact terms were determined in (9) by a fit to the scattering data. Since the data cannot be described with a LO EFT, constraints from the interaction were not implemented explicitly. As shown in Ref. (9), a good description of the 35 lowenergy scattering can be obtained for cutoff values MeV and for natural values of the LECs. The sixth LEC () is only present in the channels with isospin zero, cf. Table 1. There is scarce experimental information on these channels that could be used to fix this LEC, but it turned out that the quality of the existing data do not really allow to constrain its value reliably (10). Even with the value of the sixth LEC chosen so that , agreement with those data can be achieved. In this case a scattering length of fm (10) is obtained. Analyses of the measured binding energy of the doublestrange hypernucleus (23) suggest that the scattering length could be in the range of to fm (24); (25); (26).
Iii Quark mass dependence of the binding energy and SU(3) breaking effects
Chiral effective field theory itself does not allow one to make any predictions with regard to the existence of the dibaryon because, as stated, one of the contact terms () occurs just in the channel in question, cf. Table 1, and has to be determined from there. If we take over the values for and as fixed from the data and assume that then we find no bound state for  neither in LO (9) nor based on the preliminary NLO results (27). However, if we assume that and vary its value then in both cases a nearthreshold bound state can be produced for values of natural size.
At the same time, the framework of chiral effective field theory in which our interaction is derived is very well suited to shed light on the general characteristics of a dibaryon, should it indeed exist. In particular, it allows us study the implications of the imposed (approximate) SU(3) flavor symmetry and to explore the dependence of the properties of an assumed dibaryon on the masses of the relevant mesons and baryons. The latter aspect is a rather crucial issue in view of the fact that the available lattice QCD calculations were not performed at the physical masses of the involved particles.
To start our discussion, let us assume that the dibaryon is a (loosely) bound state and that its binding energy is similar to that of the deuteron . For pedagogical purposes we fix the value of the flavorsinglet LEC in such a way that fm (, where is either or ), because of the wellknown relation between the binding energy and the effective range parameters (28); (29)
This relation is very well fulfilled for the deuteron and the corresponding neutronproton scattering length (fm) and effective range (fm). One would naively expect that the same should happen for the dibaryon. However, it turns out that the corresponding results for in the partial wave are quite different, namely fm and fm. Specifically, the effective range is much larger and, moreover, negative. Clearly, the properties of the dibaryon are not comparable to those of the deuteron, despite the fact that both bound states are close to the elastic threshold. Indeed, if one recalls the expressions for the relevant potentials as given in Table 1,
one can see that the attraction supplied by the SU(3) flavorsinglet state () contributes with a much larger weight to the channel than to . This indicates that the presumed dibaryon could be predominantly a bound state. We have confirmed this conjecture by evaluating explicitly the phase shifts in the and channels, cf. the discussion in the next section. Indeed, one finds that the phase shift for the channel is rather similar to the case. Specifically, the () phase shift fulfills in agreement with the Levinson theorem. The () phase behaves rather differently and satisfies . Note that there have been earlier discussions on this issue in the context of baryonbaryon interactions derived within the quark model (30); (31).
Let us now consider variations of the masses of the involved particles. The dependence of the binding energy on the pion mass is displayed in Fig. 1 (left). For the case considered above, enlarging the pion mass to around 400 MeV (i.e. a value corresponding to the NPLQCD calculation (3)) increases the binding energy to around 8 MeV and a further change of to 700 MeV (corresponding roughly to the HAL QCD calculation (4)) yields then 13 MeV, cf. the solid line. Readjusting so that we predict a binding energy of 13.2 MeV for MeV, corresponding to the latest result published by NPLQCD (6), yields the dashed curve. Fig. 1 includes also results of a calculation where was fixed in order to reproduce their earlier value of around 16 MeV (3) (dashdotted curve), to facilitate a more direct comparison with chiral extrapolations (7); (8) based on that value. It is obvious that the dependence on we get agrees well – at least on a qualitative level – with that presented in Ref. (7). Specifically, our calculation exhibits the same trend (a decrease of the binding energy with decreasing pion mass) and our binding energy of 9 MeV at the physical pion mass is within the error bars of the results given in (7). On the other hand, we clearly observe a nonlinear dependence of the binding energy on the pion mass. As a consequence, scaling our results to the binding energy reported by the HAL QCD collaboration (4) (3040 MeV for MeV) yields binding energies of more than 20 MeV at the physical point, which is certainly outside of the range suggested in Ref. (7). However, we note that for such large pion masses the LO chiral EFT can not be trusted quantitatively. We remark that in our simulations the curves corresponding to different binding energies remained roughly parallel even up to such large values as suggested by the HAL QCD collaboration.
Our finding that any dibaryon is very likely a bound state rather than a state, which follows from the assumed (approximate) SU(3) symmetry of the interaction, suggests that not only the pion mass but also the masses of the baryons play a significant role for the concrete value of binding energy. Indeed, the physical difference between the and thresholds of around 26 MeV implies that the dibaryon considered above is, in reality, bound by roughly 28 MeV with respect to its “proper” threshold. Accordingly, one intuitively expects that in a fully SU(3) symmetric case, where the masses of all octet baryons coincide, the bound state would remain basically fixed to the threshold and then would lie also about 28 MeV below the threshold. In the concrete case of , , we are dealing with three coupled channels, namely , , and . Since we know from our experience with coupledchannel problems (9); (32); (33); (34) that coupling effects are sizeable and the actual separation of the various thresholds plays a crucial role we investigated also the dependence of the binding energy on the thresholds (i.e. on the , and on the and masses). Corresponding results are displayed in Fig. 1 on the right side. The solid curve is again the result based on our reference case with a binding energy of 1.87 MeV for physical masses of the pion and the baryons. When we now decrease the mass so that its nominal threshold of 2385 MeV moves downwards and finally coincides with the one of the channel (2231.2 MeV), we observe a rather drastic change in the binding energy. Note that the direct interaction in the channel is actually repulsive for the lowenergy coefficients and fixed from the data plus the pseudoscalar meson exchange contributions with coupling constants determined from the SU(3) relations Eq. (3). And it remains repulsive even for values that produce a bound dibaryon. But the coupling between the channels generates a sizeable effective attraction which increases when the channel thresholds come closer. The dashed curve is a calculation with the contact term fixed to simulate the binding energy (13.2 MeV) of the NPLQCD collaboration at MeV. As one can see, the dependence on the binding energy on the mass is rather similar. The curve is simply shifted downwards by around 4.5 MeV, i.e. by the difference in the binding energy observed already at the physical masses. The asterisks and crosses represent results where, besides the variation of the threshold, the threshold is shifted to coincide with that of the channel. This produces an additional increase of the binding energy by 20 MeV at the physical threshold and by 9 MeV for that case where all three threshold coincide. Altogether there is an increase in the binding energy of roughly 60 MeV when going from the physical point to the case of baryons with identical masses. This is significantly larger than the variations due to the pion mass considered before. Note that we have kept the pion mass at its physical value while varying the thresholds.
Iv Application to lattice QCD results
After these exemplary studies let us now try to connect with the published binding energies from the lattice QCD calculations (3); (4). The results obtained by the HAL QCD collaboration are obviously for the SU(3) symmetric case and the corresponding masses are given in Table I of Ref. (4). Thus, we can take those masses and then fix the LEC so that we reproduce their binding energy with those masses. To be concrete, we use MeV and MeV, and fix so that MeV. When we now let the masses of the baryons and mesons go to their physical values the bound state moves upwards, crosses the threshold, crosses also the threshold and then disappears. In fact, qualitatively this outcome can be already read off from the curves in Fig. 1 by combining the effects from the variations in the pion and the baryon masses. Based on those results one would expect a shift of the binding energy in the order of 60 to 70 MeV for the mass parameters of the HAL QCD calculation.
In case of the NPLQCD calculation we take the values provided in Ref. (35), i.e. MeV, MeV, MeV, and MeV. Those yield then 17 MeV for the  threshold separation (to be compared with the physical value of roughly 26 MeV) and 77 MeV for the  separation (physical value around 154 MeV). With those baryon masses we fix again the LEC so that we reproduce the binding energy given by the NPLQCD collaboration, namely MeV (6). We use also MeV, but we take the physical masses for the other pseudoscalar mesons ( and ). Again we let the masses of the baryons and of the pion go to their physical values. Also here the bound state moves upwards and crosses the threshold. However, in the NPLQCD case the state survives and remains below the threshold at the physical point. Specifically, we observe a resonance at a kinetic energy of 21 MeV in the system or, more precisely, a quasibound state in the system around 5 MeV below its threshold.
The results reported so far were all obtained with the LECs as fixed from the data for the cutoff value MeV (9). In order to investigate the stability of our results we considered also variations of the cutoff mass in the LS equation Eq. (5) between 550 and 700 MeV, as in (10), and repeated the exercise described above. Those variations led to changes of the predicted resonance energy in the order of 1 MeV only. Even an exploratory calculation that utilizes the lowenergy coefficients determined in a preliminary NLO study of the system (27) yielded practically the same result.
It is interesting to observe that the chiral extrapolation of the lattice QCD results performed by Shanahan et al. (8) yields results that are similar to ours. In that reference the authors conclude that the dibaryon is likely to be unbound by 1314 MeV at the physical point. Let us emphasize, however, that our values are not really comparable with theirs. In our analysis we assume that the dibaryon is actually a bound state – which seems to be the case also in the lattice QCD studies (3); (4). On the other hand, in Ref. (8) it is assumed that the is a compact, multiquark state rather than a loosely bound molecular state. How such a genuine multiquark state would be influenced by variations of the thresholds is completely unclear. It depends, among other things, on whether and how strongly this state couples to the , , and channels. So far there is no information on this issue from lattice QCD calculations. Clearly, in case of a strong and predominant coupling to the alone, variations of the and would not influence the binding energy significantly. However, should it couple primarily to the and/or channels then we expect a sensitivity to their thresholds comparable to what we found in our study for the case of a bound state.
Finally, for illustrative purposes, let us show phase shifts for the partial wave of the and channels. This is done in Fig. 2. The solid line is the result for our reference interaction that produces a loosely bound dibaryon with MeV. The phase shift for the channel (left side) is rather similar to the one for the partial wave where the deuteron resides, see e.g. (22). Specifically, it starts at , decreases smoothly and eventually approaches zero for large energies, fulfilling the Levinson theorem. The result for () (left side) behaves rather differently. This phase commences at zero degrees, is first negative but becomes positive within 20 MeV and finally turns to zero again for large energies. The dashed curve corresponds to the interaction that was fitted to the result of the HAL QCD collaboration and reproduces their bound dibaryon with their meson and baryon masses. The phase shift of the channel, calculated with physical masses, shows no trace of a bound state anymore. Still the phase shift rises up to around near threshold, a behavior quite similar to that of the partial wave where there is a virtual state (also called antibound state (36)). Indeed, such a virtual state seems to be present too in the channel as a remnant of the original bound state. The effect of this virtual state can be seen in the phase shift where it leads to an impressive cusp at the opening of the channel, cf. the dashed line on the left side. In the phase shifts for the NPLQCD case (dashdotted curve) the presence of a bound state is clearly visible. The corresponding phase shift exhibits a resonancelike behavior at the energy where the (quasibound) dibaryon is located.
Given the present uncertainties of the lattice QCD results for the dibaryon binding energy we have focused here primarily on qualitative features of the that can be inferred from chiral effective field theory. Clearly, once more precise lattice data become available one should also perform a careful assessment of the uncertainties involved in the EFT calculation.
Acknowledgements.
We would like to thank S. R. Beane, A. Gal, and M. J. Savage for their constructive comments. This work is supported by the Helmholtz Association by funds provided to the virtual institute “Spin and Strong QCD” (VHVI231), by the EUResearch Infrastructure Integrating Activity “Study of Strongly Interacting Matter” (HadronPhysics2, grant n. 227431) under the Seventh Framework Program of the EU, and by the DFG (SFB/TR 16 “Subnuclear Structure of Matter”).Footnotes
 Because of the GellMannOakesRenner relation, the pion mass squared is proportional to the average light quark mass. Therefore, the notions “quark mass dependence” and “pion mass dependence” can be used synonymously.
References
 R. L. Jaffe, Phys. Rev. Lett. 38, 195 (1977) [Erratumibid. 38, 617 (1977)].
 C.J. Yoon et al., Phys. Rev. C 75, 022201 (2007).
 S. R. Beane et al., Phys. Rev. Lett. 106, 162001 (2011).
 T. Inoue et al.,, Phys. Rev. Lett. 106, 162002 (2011).
 T. Inoue, [arXiv:1109.1620 [heplat]].
 S. R. Beane et al., [arXiv:1109.2889 [heplat]].
 S. R. Beane et al., [arXiv:1103.2821 [heplat]].
 P. E. Shanahan, A. W. Thomas, R. D. Young, [arXiv:1106.2851 [nuclth]].
 H. Polinder, J. Haidenbauer, U.G. Meißner, Nucl. Phys. A 779, 244 (2006).
 H. Polinder, J. Haidenbauer, U.G. Meißner, Phys. Lett. B 653, 29 (2007).
 S. R. Beane, M. J. Savage, Nucl. Phys. A 713, 148 (2003).
 S. R. Beane, M. J. Savage, Nucl. Phys. A 717, 91 (2003).
 E. Epelbaum, U.G. Meißner, W. Glöckle, Nucl. Phys. A 714, 535 (2003).
 E. Epelbaum, U.G. Meißner, W. Glöckle, [arXiv:nuclth/0208040].
 M. J. Savage, M. B. Wise, Phys. Rev. D 53, 349 (1996).
 C. L. Korpa, A. E. L. Dieperink, R. G. E. Timmermans, Phys. Rev. C 65, 015208 (2002).
 S. R. Beane, P. F. Bedaque, A. Parreño, M. J. Savage, Nucl. Phys. A 747, 55 (2005).
 J. J. de Swart, Rev. Mod. Phys. 35, 916 (1963).
 C. B. Dover, H. Feshbach, Annals Phys. 217, 51 (1992).
 P. G. Ratcliffe, Phys. Lett. B 365, 383 (1996).
 E. Epelbaum, W. Glöckle, U.G. Meißner, Nucl. Phys. A 637, 107 (1998).
 E. Epelbaum, W. Glöckle, U.G. Meißner, Nucl. Phys. A 747, 362 (2005).
 H. Takahashi et al., Phys. Rev. Lett. 87, 212502 (2001).
 I. N. Filikhin, A. Gal, V. M. Suslov, Phys. Rev. C 68, 024002 (2003).
 T. A. Rijken, Y. Yamamoto, Phys. Rev. C 73, 044008 (2006).
 Y. Fujiwara, Y. Suzuki, C. Nakamoto, Prog. Part. Nucl. Phys. 58, 439 (2007).
 J. Haidenbauer, EPJ Web of Conferences 3, 01009 (2010).
 J. Schwinger, Phys. Rev. 72, 742 (1947).
 H. A. Bethe, Phys. Rev. 76, 38 (1949).
 M. Oka, K. Shimizu, K. Yazaki, Phys. Lett. 130B, 365 (1983).
 C. Nakamoto, Y. Suzuki, Y. Fujiwara, Prog. Theor. Phys. 97, 761 (1997).
 J. Haidenbauer, U.G. Meißner, Phys. Rev. C 72, 044005 (2005).
 J. Haidenbauer, U.G. Meißner, Phys. Lett. B684, 275 (2010).
 J. Haidenbauer, G. Krein, U.G. Meißner, L. Tolos, Eur. Phys. J. A 47, 18 (2011).
 S. R. Beane et al., Phys. Rev. D 84, 014507 (2011).
 see, e.g., B. C. Pearce, B. F. Gibson, Phys. Rev. C 40, 902 (1989), and references therein for a discussion of poles and bound states in the context of coupledchannel systems.