## Is d ∗ a candidate of hexaquark-dominated exotic state?

F. Huang, 1 Z.Y. Zhang, 2, 3 P.N. Shen, 2, 3, 4 and W.L. Wang 5

1

School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 101408, China 2 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China 3 Theoretical Physics Center for Science Facilities, CAS, Beijing 100049, China 4 College of Physics and Technology, Guangxi Normal University, Guilin 541004, China 5 School of Physics and Nuclear Energy Engineering, Beihang University, Beijing 100191, China (Dated: March 16, 2015- dstar˙v7)

We confirm our previous prediction of a d ∗ state with I ( J P ) = 0(3 + ) [Phys. Rev. C 60 , 045203 (1999)] and report for the first time based on a microscopic calculation that d ∗ has about 2/3 hidden color (CC) configurations and thus is a hexaquark-dominated exotic state. By performing a more elaborate dynamical coupled-channels investigation of the ∆∆-CC system within the framework of resonating group method (RGM) in a chiral quark model, we found that the d ∗ state has a mass of about 2 38 . -2 42 GeV, a root-mean-square radius (RMS) of 0 76 . . -0 88 fm, and a CC fraction . of 66% -68%. The last may cause a rather narrow width to d ∗ which, together with the quantum numbers and our calculated mass, is consistent with the newly observed resonance-like structure ( M ≈ 2380 MeV, Γ ≈ 70 MeV) in double-pionic fusion reactions reported by WASA-at-COSY Collaboration.

PACS numbers: 14.20.Pt, 13.75.Cs, 12.39.Jh, 24.10.Eq

The ABC effect has drawn physicists' great attention since its observation in 1961 in the pd reaction [1]. In recent years, much experimental progress in exploring the nature of the ABC effect has been made. In 2009, the CELSIUS/WASA Collaboration measured the most basic double-pionic fusion reaction pn → dπ π 0 0 with an incident proton energy of 1.03 GeV and 1.35 GeV [2], and found significant enhancements in the ππ invariant mass spectrum at ππ invariant mass below 0.32 GeV 2 and also in the dπ invariant mass spectrum at ∆ resonance region. To accommodate these data as well as the energy dependence of the total cross section at √ s &lt; 2 5 . GeV, the conventional t -channel ∆∆ intermediate state is found to be not sufficient, and a new structure, namely an s -channel resonance with mass of about 2.36 GeV and width of about 80 MeV, is expected. In 2011, the WASA-at-COSY Collaboration further measured the pn → dπ π 0 0 reaction with the beam energies of 1 0 . -1 4 GeV which cover the transition region of the con-. ventional t -channel ∆∆ process [3]. They found that neither the t -channel ∆∆ process nor the Roper resonance process can explain the data, and an s -channel resonance with quantum numbers of I ( J P ) = 0(3 + ), mass of about 2.37 GeV and width of about 70 MeV is indeed needed to describe the data. Recently, the WASA-at-COSY Collaboration measured the polarized /vector np scattering through the quasi-free process /vector dp → p spectator np [4, 5]. By incorporating the newly measured A y data into the SAID analysis, they obtained a pole in the 3 D 3 -3 G 3 waves at (2380 ± 10)+ (40 i ± 5) MeV, which again supports the existence of a resonance, called d ∗ , as mentioned in Ref. [3]. Further evidence of this resonance has also been reported in the quasi-free np → npπ π 0 0 reaction [6]. Since its mass is above the threshold of ∆ Nπ channel, while its width is much smaller than the decay width of ∆, this resonance must be a very interesting state involving new physical mechanisms and it is obviously worthwhile investigating.

Theoretically, the possibility of the existence of dibaryon states was first proposed in 1964 by Dyson and Xuong based on SU(6) symmetry [7]. Since then, extensive efforts have been given in exploring the possible existence of a ∆∆ dibaryon on hadronic degrees of freedom. However, no convincing results have ever been released yet. Since the birth of quark model, dramatic progresses on this aspect were pouring in. In 1980, by analyzing the characteristics of the one-gluon-exchange (OGE) interaction between quarks, Oka and Yazaki pointed out that in all non-strange baryon-baryon (BB) systems, the ∆∆ system with I ( J P ) = 0(3 + ) is the only one in which the effective BB interaction induced by OGE shows an attractive feature [8]. In the following years, by including the interaction between the quark field and chiral field into the constituent quark model, one successfully reproduced the data of the nucleon-nucleon (NN) interaction and the binding energy of the deuteron [9], which would provide a much reliable platform to predict the structures of dibaryons in the quark degrees of freedom. In 1999, we carefully performed a dynamical study of the ∆∆ system in the quark degrees of freedom within the framework of Resonating Group Method (RGM), with the hidden color (CC) channel being properly taken into account [10]. Where, by employing the chiral SU(3) quark model with a set of reasonable model parameters which can reproduce the NN scattering phase shifts at relatively lower energies and the binding energy of deuteron, we found that the quark-exchange effect in the ∆∆ system with quantum number I ( J P ) = 0(3 + ) is so important (see also Ref. [11]) that the system should be bound in nature with a binding energy of about 20 -50 MeV in a single ∆∆ channel calculation. The coupling to the CC

channel was also intensively studied and found to play an important role in the binding behavior of the system. It offers an additional binding energy of over 20 MeV to the system, and thus the binding energy of d ∗ , relative to the threshold of ∆∆ channel, would run up to 40 -80 MeV. Unexpectedly, our predicted mass and quantum numbers are quite close to the new observation released recently by WASA-at-COSY Collaboration [4].

In this work, we perform a further elaborate investigation of the ∆∆ -CC system within the framework of RGM in a chiral quark model. Besides the binding energy, we concentrate on a detailed study of the relative wave function of the ∆∆ -CC system, which is crucial to the understanding of the structure and decay property of d ∗ . We find for the first time based on a microscopic calculation that the d ∗ has a CC component of about 2/3, which indicates that d ∗ is a hexaquark-dominated exotic state. The large CC configuration is quite helpful for us to understand that d ∗ , although locating above the thresholds of ∆ Nπ and NNππ channels, has a relatively narrow width as will be discussed later in more detail. In the most recent experimental papers by WASA-at-COSY Collaboration and SAID Data Analysis Center [5, 6], our scenario for d ∗ present here has been cited as a plausible explanation of the resonance observed in the doublepionic fusion reaction, as both our calculated binding energy and the expected width from our picture of d ∗ are in agreement with the experimental findings.

The interaction between i -th quark and j -th quark in our chiral quark model reads

$$V _ { i j } = V _ { i j } ^ { \text{OGE} } + V _ { i j } ^ { \text{conf} } + V _ { i j } ^ { \text{ch} }, \quad \quad ( 1 )$$

where V OGE ij is the OGE interaction which describes the short-range perturbative QCD behavior, and V conf ij the confinement potential describing the long-range nonperturbative QCD effects. V ch ij is the chiral field induced quark-quark interaction which provides the mediumrange non-perturbative QCD effects, and to test the model dependence of our results, we employ two different models for this interaction. In the chiral SU(3) quark model, V ch ij reads

$$V _ { i j } ^ { \text{ch} } = \sum _ { a = 0 } ^ { 8 } \left ( V _ { i j } ^ { \sigma _ { a } } + V _ { i j } ^ { \pi _ { a } } \right ), \quad \ \ ( 2 ) \quad 0 ( 3 )$$

and in the extended chiral SU(3) quark model, V ch ij reads

$$V _ { i j } ^ { \text{ch} } = \sum _ { a = 0 } ^ { 8 } \left ( V _ { i j } ^ { \sigma _ { a } } + V _ { i j } ^ { \pi _ { a } } + V _ { i j } ^ { \rho _ { a } } \right ), \quad \ \ ( 3 ) \quad \text{$\mathrm RM$} \quad \text{$the$} \\ \text{$to 2$} \\ \text{$\mathrm fm$.}$$

with σ a , π a and ρ a ( a = 0 1 , , · · · , 8) being the scalar, pseudo-scalar and vector nonet fields, respectively. Note that the OGE will be largely reduced when vector-meson exchanges are included, i.e. the short-range interaction

TABLE I. Binding energy, root-mean-square radius (RMS) of 6 quarks, and fraction of channel wave function for deuteron in chiral SU(3) quark model and extended chiral SU(3) quark model with ratio of tensor coupling to vector coupling f/g=0 and f/g=2/3 for vector meson fields.

|                        | SU(3)   | Ext. SU(3)   | Ext. SU(3)   |
|------------------------|---------|--------------|--------------|
|                        |         | (f/g=0)      | (f/g=2/3)    |
| Binding energy (MeV)   | 2.09    | 2.24         | 2.20         |
| RMS of 6 q (fm)        | 1.38    | 1.34         | 1.35         |
| Fraction (NN) L =0 (%) | 93.68   | 94.66        | 94.71        |
| Fraction (NN) L =2 (%) | 6.32    | 5.34         | 5.29         |

mechanisms are quite different in these two models. We refer readers to Refs. [9, 12] for further details.

The parameters of both models are fixed by fitting the energies of the octet and decuplet baryon ground states, the NN scattering phase shifts in the low energy region ( √ s /lessorequalslant 2 m N +200 MeV) and the binding energy of deuteron [9, 12]. See Table I for the binding energy, root-mean-square radius (RMS) of 6 q , and the fraction of each partial wave for deuteron obtained from our models. They are all quite reasonable.

With all the parameters being properly fixed, we are eager to a investigation for the properties of the ∆∆-CC system without introducing any new adjustable parameters. Here the hidden color channel CC with isospin I = 0 and spin S = 3 is built as

$$\kappa \, \ln \, \left ( \lambda \right ) \, \lambda \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \cdots \, \\ ( 1 ) \, \left ( \lambda \right ) \, \right ) _ { I S = 0 3 } \equiv - \frac { 1 } { 2 } \left | \Delta \Delta \right \rangle _ { I S = 0 3 } + \frac { \sqrt { 5 } } { 2 } \mathcal { A } ^ { s f c } \left | \Delta \Delta \right \rangle _ { I S = 0 3 }, \ ( 4 )$$

with A sfc being the antisymmetrizer in spin-flavor-color space. We dynamically calculate the binding energy of this system by solving the RGM equation for a bound state problem, and then discuss the structure of this bound state via a systematic analysis of fractions of each channel in the resultant relative wave function. The results obtained in chiral SU(3) quark model and the extended chiral SU(3) quark model with the ratios of the tensor coupling to vector coupling f/g=0 and f/g=2/3 for vector meson fields are listed in Table II, where the values for RMS of 6 quarks are also shown.

One sees from Table II that the ∆∆ state with I ( J P ) = 0(3 + ) has a binding energy of about 30 -60 MeV and a RMS of about 0 80 . -0 96 fm in a ∆∆ ( . L = 0 2) double-, channel calculation. The coupling to the CC channel will further result in an increment of about 20 MeV to the binding energy and a considerable decrement of the RMS, and finally, the mass of this bound state will reach to 2 38 . -2 42 GeV and the RMS will shrink to 0 76 . . -0 88 . fm. This clearly shows that d ∗ is a ∆∆-CC deeply bound and compact state where the coupling to the CC channel plays a significant role.

Apart from the RMS, more information about the 'size' of d ∗ can be acquired from its spacial distribu-

TABLE II. Binding energy, root-mean-square radius (RMS) of 6 quarks, and fraction of channel wave function for d ∗ in chiral SU(3) quark model and extended chiral SU(3) quark model with ratio of tensor coupling to vector coupling f/g=0 and f/g=2/3 for vector meson fields.

|                        | ∆∆ ( L = 0 , 2)   | ∆∆ ( L = 0 , 2)    | ∆∆ ( L = 0 , 2)      | ∆∆ - CC ( L = 0 , 2)   | ∆∆ - CC ( L = 0 , 2)   | ∆∆ - CC ( L = 0 , 2)   |
|------------------------|-------------------|--------------------|----------------------|------------------------|------------------------|------------------------|
|                        | SU(3)             | Ext. SU(3) (f/g=0) | Ext. SU(3) (f/g=2/3) | SU(3)                  | Ext. SU(3) (f/g=0)     | Ext. SU(3) (f/g=2/3)   |
| Binding energy (MeV)   | 28.96             | 62.28              | 47.90                | 47.27                  | 83.95                  | 70.25                  |
| RMS of 6 q (fm)        | 0.96              | 0.80               | 0.84                 | 0.88                   | 0.76                   | 0.78                   |
| Fraction (∆∆) L =0 (%) | 97.18             | 98.01              | 97.71                | 33.11                  | 31.22                  | 32.51                  |
| Fraction (∆∆) L =2 (%) | 2.82              | 1.99               | 2.29                 | 0.62                   | 0.45                   | 0.51                   |
| Fraction (CC) L =0 (%) |                   |                    |                      | 66.25                  | 68.33                  | 66.98                  |
| Fraction (CC) L =2 (%) |                   |                    |                      | 0.02                   | 0.00                   | 0.00                   |

/s32

/s32

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0458v3-figure-3b266540-b7da-424c-9096-6cf2de714caf.jpg)

/s32

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0458v3-figure-7af39187-aa23-4607-9325-08248232c198.jpg)

/s32

FIG. 1. Relative wave functions in the extended chiral SU(3) quark model with f/g=0 for deuteron (left) and d ∗ (right).

tion feature as shown in Fig. 1, where the partial-wave projected relative wave function in physical basis as a function of the distance between two physical states in the extended chiral SU(3) quark model with f/g=0 is plotted. The results in other two models are similar. For comparison, the results for deuteron are also plotted. Here the relative wave function in physical basis, namely the channel wave function in the quark cluster model, is defined in the usual way as [13-15]

$$\psi _ { \text{BB} } ( \boldsymbol r ) \equiv \langle \phi _ { \text{B} } ( \boldsymbol \xi _ { 1 }, \boldsymbol \xi _ { 2 } ) \, \phi _ { \text{B} } ( \boldsymbol \xi _ { 4 }, \boldsymbol \xi _ { 5 } ) \, | \, \Psi _ { 6 q } \rangle \,, \quad ( 5 ) \quad \text{loca} ^ { 1 }$$

with B standing for ∆ or C and Ψ 6 q being the six-quark wave function from RGM calculation,

In Eqs. (5) and (6), φ B is the antisymmetrized internal wave function for cluster B, with ξ i ( i = 1 2 4 5) , , , being the internal Jacobi coordinates for corresponding cluster; A is the antisymmetrizer for quarks from different clusters required by Pauli exclusion principle; and η BB ( r ) is the relative wave function between two clusters of B, which in present investigation is completely determined by the interacting dynamics of the whole sixquark system. From Fig. 1 one sees that d ∗ is rather narrowly distributed and it has a maximal distribution located around 0 7 fm for ∆∆ ( . L = 0) and 0 4 fm for . CC ( L = 0), respectively. While deuteron is widely distributed with a maximal distribution located around 1 4 . fm.

$$\Psi _ { 6 q } = \sum _ { B = \Delta, C } \mathcal { A } \left [ \phi _ { B } ( \xi _ { 1 }, \xi _ { 2 } ) \, \phi _ { B } ( \xi _ { 4 }, \xi _ { 5 } ) \, \eta _ { B B } ( r ) \right ]. \quad ( 6 ) \quad \text{thing} \\ \text{relati}$$

It should be mentioned that an even more interesting thing in this investigation is to study the fractions of relative wave functions for each individual channel, which

will help us get a further understanding of the structure of d ∗ . Our extracted fractions of relative wave functions for ∆∆ and CC channels are tabulated in Table II. Of great interest is that one observes that the fraction of the CC channel in d ∗ is about 66% -68%. Note according to symmetry, a pure hexaquark state of ∆∆-CC system with isospin I = 0 and spin S = 3 reads

$$[ 6 ] _ { o r b } [ 3 3 ] _ { I S = 0 3 } = \sqrt { \frac { 1 } { 5 } } \left | \Delta \Delta \right \rangle _ { I S = 0 3 } + \sqrt { \frac { 4 } { 5 } } \left | \CC \right \rangle _ { I S = 0 3 }, \ ( 7 ) \quad \text{$\text{and} \,m$} \\ \text{partial partial}$$

which indicates that the fraction of CC channel in a pure hexaquark state is 80%. It is thus fair to say that d ∗ is a hexaquark-dominated state as it has a CC configuration of 66% -68%. This finding is of great interest. It helps us to understand that d ∗ , although locating above the thresholds of the ∆ Nπ and NNππ channels, has a relatively narrow width since the CC component cannot subject to a direct break-up decay. Actually, it can only decay to colorless hadrons via the re-combination processes of six color quarks together with quark-antiquark pair creation processes, which are highly suppressed and have much lower probability than those direct break-up decay processes. A conjecture that d ∗ should have an unconventional origin and the CC configuration will suppress its decay width has also been proposed by Bashkanov, Brodsky and Clement in Ref. [16]. Our microscopic calculation presented here supports their argument and our dynamical results show that d ∗ has about 2/3 hiddencolor configurations. However, a detailed calculation of the decay width of d ∗ needs to take into account both the kinematic effects and the effects from the dynamical structure of d ∗ . In the literature the former has been considered by using a simple momentum dependent prescription for d ∗ s decay width [17]. Now an investigation of the later also becomes possible, because the relative wave functions of d ∗ is available in the present paper. Our work along this line is in progress.

It should particularly be stressed that the abovementioned features of d ∗ are due to both the quark exchange effect and the short-range interaction being attractive in the ∆∆-CC system. As we have pointed out in Ref. [11], the quark exchange effect is highly dependent on the quantum numbers of the system. Opportunely, the ∆∆ system with I ( J P ) = 0(3 + ) is one of the few systems which have strong quark exchange effect that drags two baryons together to form a compact state. As for the interaction property, Oka and Yazaki claimed in 1980 that in all the non-strange two-baryon systems, the ∆∆ system with I ( J P ) = 0(3 + ) is the only one in which the OGE provides strong attraction in short range [8]. In our chiral SU(3) quark model, the OGE indeed provides strong short-range attraction. In the extended chiral SU(3) quark model, although the OGE is largely reduced, the short-range interaction is still attractive and the attraction is even much stronger, as the vector meson exchanges (VMEs) are also strongly attrac- tive in short range. Considering both mentioned facts, the ∆∆ system with I ( J P ) = 0(3 + ) is certainly deeply bound and it would couple to the CC channel strongly as two interacting ∆s are dragged closer enough by strong attraction. As a comparison, the NN 3 S 1 partial wave has a rather different feature. In this partial wave, the quark exchange effect is very week (almost negligible), and moreover, the short-range interaction stemming from OGE and VMEs are all repulsive. Therefore, the NN 3 S 1 partial wave can only get attraction in the medium- and long-range through σ and π meson exchanges, and as a result, the deuteron is loosely bound and hardly couples to CC channel. We emphasize that the ∆∆ system with I ( J P ) = 0(3 + ) is a highly special system where both the quark exchange effect and the short-range interaction are attractive, which makes this system deeply bound and promotes a strong coupling to the CC channel.

In summary, it is the first time one finds based on a microscopic calculation with no additional parameters besides those fixed already in the study of NN scattering phase shifts that d ∗ has about 2/3 hidden-color configurations and thus tends to be a hexaquark-dominated exotic state. It has a mass of about 2 38 . -2 42 GeV, a root-. mean-square radius of 0 76 . -0 88 fm, and a CC fraction of . 66% -68% which may cause a relatively narrow width of d ∗ . Our findings are consistent with the newly observed resonance-like structure ( M ≈ 2380 MeV, Γ ≈ 70 MeV) in double-pionic fusion reactions reported by WASA-atCOSY Collaboration. We mention that if its character can be further verified, the d ∗ will be the first hexaquarkdominated exotic state we have ever found, and it may open a door to new physical phenomena.

In literature there are several other interpretations of the recent WASA-at-COSY experiment. Gal et al. carried out a πN ∆ three-body-calculation [18] and found the D 03 state as a dynamically generated pole at the right mass and slightly larger width (see Ref. [5] for a comment of the width). Huang et al. performed a coupled-channel quark model calculation and obtained an energy close to the observed value, but their calculated width is still too large [17].

The πN ∆ three-body resonance scenario proposed by Gal et al. [18] might have a relatively larger RMS than that from our picture. While measuring the 'size' of d ∗ is rather involved, experiments in future might reach d ∗ by electromagnetic transitions, which will give information about d ∗ form factor and further tell us the real structure of the observed resonance [19]. We look forward to experimental progresses along this line.

We are grateful for constructive discussions with Prof. H. Clement and Prof. R. L. Workman.

This work is partly supported by the National Natural Science Foundation of China under grants Nos. 11475181, 11105158, 11035006 and 11165005, the Keyproject by the Chinese Academy of Sciences under project No. KJCX2-EW-N01. F.H. is grateful to the

support of the One Hundred Person Project of the University of Chinese Academy of Sciences.

- [1] N. E. Booth et al. , Phys. Rev. Lett. 7 , 35 (1961).
- [2] M. Bashkanov et al. , Phys. Rev. Lett. 102 , 052301 (2009).
- [3] P. Adlarson et al. (WASA-at-COSY Collaboration), Phys. Rev. Lett. 106 , 242302 (2011).
- [4] P. Adlarson et al. (WASA-at-COSY Collaboration &amp; SAID DAC), Phys. Rev. Lett. 112 , 202301 (2014).
- [5] P. Adlarson et al. (WASA-at-COSY Collaboration &amp; SAID DAC), Phys. Rev. C 90, 035204 (2014).
- [6] P. Adlarson et al. (WASA-at-COSY Collaboration),
- arXiv:1409.2659 [nucl-ex].
- [7] F. J. Dyson and N. H. Xuong, Phys. Rev. Lett. 13 , 815 (1964).
- [8] M. Oka and K. Yazaki, Phys. Lett. B 90 , 41 (1980).
- [9] Z. Y. Zhang et al. , Nucl. Phys. A 625 , 59 (1997).
- [10] X. Q. Yuan et al. , Phys. Rev. C 60 , 045203 (1999).
- [11] Q. B. Li et al. , Nucl. Phys. A 683 , 487 (2001).
- [12] L. R. Dai et al. , Nucl. Phys. A 727 , 321 (2003).
- [13] A. M. Kusainov et al. , Phys. Rev. C 44 , 2343 (1991).
- [14] L. Ya. Glozman, V. G. Neudatchin, and I. T. Obukhovsky, Phys. Rev. C 48 , 389 (1993).
- [15] Fl. Stancu, S. Pepin, and L. Ya. Glozman, Phys. Rev. C 56 , 2779 (1997).
- [16] M. Bashkanov, S.J. Brodsky, and H. Clement, Phys. Lett. B 727 , 438 (2013).
- [17] H. X. Huang et al. , Phys. Rev. C 89 , 034001 (2014).
- [18] A. Gal et al. , Phys. Rev. Lett. 111 , 172301 (2013).
- [19] H. Clement, private communication.
F. Huang Z. Y. Zhang P. N. Shen W. L. Wang
Is $d^*$ a candidate of hexaquark-dominated exotic state?
We confirm our previous prediction of a $d^*$ state with $I(J^P)=0(3^+)$
[Phys. Rev. C 60, 045203 (1999)] and report for the first time based on a
microscopic calculation that $d^*$ is a hexaquark dominated exotic state as it
has about 2/3 hidden color (CC) configurations. By performing a more elaborate
dynamical coupled-channels investigation of the $\Delta\Delta$-CC system within
the framework of resonating group method (RGM) in a chiral quark model, we
found that the $d^*$ state has a mass of about 2.38-2.42 GeV, a
root-mean-square radius (RMS) of 0.76-0.88 fm, and a CC fraction of 66%-68%.
The last ensures that the $d^*$ has a rather narrow width which, together with
the quantum numbers and our calculated mass, is consistent with the newly
observed resonance-like structure ($M\approx 2380$ MeV, $\Gamma\approx 70$ MeV)
in double-pionic fusion reactions reported by WASA-at-COSY Collaboration.