Next Article in Journal
The Casimir Interaction between Spheres Immersed in Electrolytes
Previous Article in Journal
Dynamical Analysis of Charged Dissipative Cylindrical Collapse in Energy-Momentum Squared Gravity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genuine Resonance of Full-Charm Tetraquarks

College of Science, Jinling Institute of Technology, Nanjing 211169, China
Submission received: 20 March 2021 / Revised: 13 May 2021 / Accepted: 14 May 2021 / Published: 18 May 2021
(This article belongs to the Section High Energy Nuclear and Particle Physics)

Abstract

:
In this work, the genuine resonance states of full-charm tetraquark systems with quantum numbers J P C = 0 + + , 1 + , 2 + + are searched in a nonrelativistic chiral quark model with the help of the Gaussian Expansion Method. In this calculation, two structures, meson-meson and diquark–antidiquark, as well as their mixing with all possible color-spin configurations, are considered. The results show that no bound states can be formed. However, resonances are possible because of the color structure. The genuine resonances are identified by the stabilization method (real scaling method). Several resonances for the full-charm system are proposed, and some of them are reasonable candidates for the full-charm states recently reported by LHCb.

1. Introduction

In the past few decades, many charmonium-like or bottomonium-like X Y Z states [1,2,3,4,5,6,7,8,9] have been observed in experiments, which has generated great challenges and opportunities for researchers to study multiquark states.
Recently, the tetraquark of the all-heavy system, such as c c c ¯ c ¯ and b b b ¯ b ¯ , has received considerable attention due to the development of experiments. If the c c c ¯ c ¯ or b b b ¯ b ¯ states steadily exist, they are most likely to be observed at LHC, Belle II and other facilities. Double J / ψ production becomes possible [10,11,12,13,14] and an enhancement in the differential production cross-section for J / ψ pairs between 6 GeV and 8 GeV can be observed [10,11]. The CMS collaboration measured pair production of Y ( 1 S ) [15]. There was also a claim for the existence of a full-bottom tetraquark state b b b ¯ b ¯ [16], with a global significance of 3.6 σ and a mass of around 18.4 GeV, almost 500 MeV below the threshold of Y Y . Very recently, the LHCb collaboration reported their newest finding that the full-charm states have been observed: there is a broad structure in the range 6.2∼6.8 GeV, a narrower structure at 6.9 GeV with a significance of about 5 σ , and a structure around 7.2 GeV [17]. This discovery brought about widespread theoretical attention [18,19,20,21,22,23,24,25].
In fact, whether or not observable states of fully-heavy tetraquarks exist has been debated for more than forty years. Theoretically, various methods are applied to study the full-heavy tetraquark states. Some work has suggested that stable bound states should exist [26,27,28,29,30,31,32,33]. Iwasaki [26] first argued that the bound state of c 2 c ¯ 2 can exist and estimated its mass to be in the neighborhood of 6 GeV or 6.2 GeV based on a string model. Heller et al. Claimed that the dimensions c c c ¯ c ¯ and b b b ¯ b ¯ are bound, and that the binding energy ranges from 0.16∼0.22 GeV based on the potential energy arising from the MIT bag model [27]. Lloyd et al. have used a parametrized Hamiltonian to calculate the spectrum of the all-charm tetraquark state and found several close-lying bound states with two sets of parameters based on large but finite oscillator bases. For example, the lowest state with quantum number J P C = 0 + + had a mass below the threshold of two η c ( 1 S ) , 5967.2 MeV [28]. Berezhnoy et al. showed that the masses of c c c ¯ c ¯ and b b b ¯ b ¯ states are under the thresholds for J = 1 , 2 by taking a diquark and antidiquark as point particles and employing the hyperfine interaction between them [29]. In a moment QCD sum rule approach, Chen et al. studied the fully-heavy tetraquark states, and discovered that the masses of the b b b ¯ b ¯ tetraquarks are below the thresholds of Y ( 1 S ) Y ( 1 S ) and η b ( 1 S ) η b ( 1 S ) , and the masses of the c c c ¯ c ¯ tetraquarks are above the threshold [30]. Yang Bai et al. and Esposito et al. calculated the mass of the b b b ¯ b ¯ state to be around 100 MeV below the threshold of η b η b [31,32]. In recent research, Debastiani et al. used a non-relativistic model to study the spectroscopy for a tetraquark composed of c c c ¯ c ¯ in a diquark–antidiquark configuration and found that the lowest S-wave c c c ¯ c ¯ tetraquarks might exist below their thresholds [33].
On the contrary, in other work, there is no bound full-heavy tetraquark state but resonance states are possible [34,35,36,37,38,39,40,41,42,43]. Chao predicted full-charm diquark–antidiquark states exist with masses in the range of 6.4∼6.8 GeV [34]. Ader et al. investigated the full-heavy tetraquark state in a potential model and found that there no state exists below the corresponding thresholds [35]. Barnea et al. studied a system consisting of quarks and antiquarks of the same flavor within the hyperspherical formalism, and the mass of c c c ¯ c ¯ was approximately 6038 MeV, which is above the corresponding threshold [36]. In the QCD sum-rule approach, no bound state was found and several resonances for c c c ¯ c ¯ and b b b ¯ b ¯ systems were proposed [37]. Karliner et al. have calculated the mass spectrum of c c c ¯ c ¯ state and found it unlikely to be less than twice the mass of the lowest charmonium state η c [38]. Wu et al. also found no full-heavy tetraquark with the same flavor below the thresholds in a color-magnetic interaction model [39]. Recently in Ref. [40], Liu et al. suggested that no bound states can be formed below the thresholds for meson pairs ( c ¯ c ) - ( c ¯ c ) within a potential model by including the linear confining potential, Coulomb potential and spin-spin interactions. The full-bottom tetraquark was also calculated in the chiral quark model with the help of the Gaussian Expansion Method, and no bound states were found and a resonance was proposed [44]. In the framework of a quark delocalization color screening model and chiral quark model with the help of the resonating group method, the reported state X ( 6900 ) can be explained as a compact resonance state with I J P = 00 + [45].
Hadron spectroscopy always plays an important role in revealing the properties of the dynamics of the strong interaction. The fully-heavy tetraquark system provides a good chance to test our understanding of hadron structure. To date, most calculations for fully-heavy tetraquarks assume that the system has a given structure, meson-meson, diquark–antidiquark, etc. A few studies have carried out calculations with structure mixing. Such as Ref. [30], the authors calculated the mass spectra for doubly hidden-charm/bottom tetraquark states using only the compact diquark–antidiquark configuration. Now most calculations are usually carried out in a finite space. The finite space will discretize the energy of the system. The bound state or genuine resonance state should have a stable energy against an increasing space and structure mixing. To obtain the genuine resonances, various techniques are available: the real scaling method (stabilization) [46,47,48], complex scaling method [49], phase shift analysis [50,51] and so on. In the present work, we systematically investigate the masses of full-charm tetraquarks with J P C = 0 + + , 1 + , 2 + + in the chiral quark model, which can well describe the properties of hadrons and hadron-hadron interactions. The method of the Gaussian Expansion Method (GEM) is employed to carry out a high precision four-body calculation. The dynamical mixing of the meson-meson configuration with the diquark–antidiquark configuration is also considered. One thing that needs to be specified is that, here, the diquark–antidiquark and meson-meson configuration is just a nomenclature for the color saturation of the system and does not correspond to the compact tetraquark and the meson molecule configurations. For example, color singlet-singlet 1 × 1 and color octet-octet 8 × 8 denotes meson-meson structure, and color antitriplet-triplet 3 ¯ × 3 and sextet-antisextet 6 × 6 ¯ denotes diquark–antidiquark structure, and their mixing is taken into account. This mixing occurs via both the spin-independent and the spin-dependent parts of the potential. To obtain the genuine resonances, the real scaling method is employed. In previous work [44], only the lowest resonances of the full-bottom tetraquark were investigated. Inspired by the recent experimental results for the full-charm tetraquark, and unlike our previous work for b b b ¯ b ¯ , the possible resonance states with higher energies are explored in the present work.
This paper is organized as follows. In Section 2, we briefly discuss the chiral quark model, the wave functions of the full-heavy tetraquark, and the Gaussian Expansion Method. In Section 3, the numerical results and discussion are presented. Some conclusions and a summary are given in Section 4.

2. Quark Model and Wave Functions

2.1. The Chiral Quark Model

The chiral quark model has been applied both in explaining the hadron spectra and hadron-hadron interactions successfully. We can find the details for the model in [52,53]. Three parts of the Hamiltonian are included for the c c c ¯ c ¯ fully-heavy system: quark rest mass, kinetic energy, and potential energy:
H = i = 1 4 m i + p 12 2 2 μ 12 + p 34 2 2 μ 34 + p 1234 2 2 μ 1234 + i < j = 1 4 V i j C + V i j G ,
m i is the constituent mass of i-th quark/antiquark, and μ i j is the reduced mass of two interacting quarks, with
μ i j = m i m j m i + m j , μ 1234 = ( m 1 + m 2 ) ( m 3 + m 4 ) m 1 + m 2 + m 3 + m 4 , p i j = m j p i m i p j m i + m j , p 1234 = ( m 3 + m 4 ) p 12 ( m 1 + m 2 ) p 34 m 1 + m 2 + m 3 + m 4 .
V C and V G represents the quark confinement and the one-gluon-exchange potential, respectively. The detailed forms for the two potentials are shown below [52]:
V i j C = ( a c r i j 2 Δ ) λ i c · λ j c , V i j G = α s 4 λ i c · λ j c 1 r i j 2 π 3 m i m j σ i · σ j δ ( r i j ) , δ ( r i j ) = e r i j / r 0 ( μ i j ) 4 π r i j r 0 2 ( μ i j ) .
λ c are the S U ( 3 ) color Gell-Mann matrices and σ are the S U ( 2 ) Pauli matrices; r 0 ( μ i j ) = s 0 / μ i j and α s is an effective scale-dependent running coupling [53],
α s ( μ i j ) = α 0 ln ( μ i j 2 + μ 0 2 ) / Λ 0 2 .
All the model parameters are listed in Table 1, which are determined by fitting the light and heavy meson spectra. In Table 2, the theoretical and experimental masses for the low-lying c c ¯ are demonstrated in the chiral quark model. Because the orbital-spin interactions are not taken into account in the present calculation, the P-wave states χ c J ( J = 0 , 1 , 2 ) have the same mass.

2.2. The Wave Functions of c c c ¯ c ¯ System

The wave functions of four-quark states can be constructed in two steps for both the meson-meson and diquark–antidiquark structure. First we obtain the wave functions for two-body sub-clusters, and, then, secondly, couple the wave functions of these two sub-clusters to construct the total wave functions for a four-quark state with exclusive quantum number I J P C . To save space, only the wave functions for each degree of freedom are shown below.
(1) Diquark–antidiquark structure
We denote α and β as the spin-up and spin-down states of the quarks and the spin wave functions for the four-quark states takes the form of
χ 00 σ 1 = 1 2 ( α β α β α β β α β α α β + β α β α ) , χ 00 σ 2 = 1 12 ( 2 α α β β α β α β α β β α β α α β β α β α + 2 β β α α ) , χ 11 σ 3 = 1 2 ( α β α α β α α α ) , χ 11 σ 4 = 1 2 ( α α α β α α β α ) , χ 11 σ 5 = 1 2 ( α α α β + α α β α α β α α β α α α ) , χ 22 σ 6 = α α α α .
There are six spin wave functions for four-quark states in total, which are marked as the superscript σ i (i = 1∼6) of χ . The subscripts of χ represent the total spin and the third projection of total spin of the system, and only one component ( M S = S ) is needed for a given total spin S.
For the c c c ¯ c ¯ system, the flavor wave function reads as follows:
χ d 0 f = ( c c ) ( c ¯ c ¯ ) .
The marks d 0 of χ represent the diquark–antidiquark structure with an isospin that equals zero.
For the color freedom, we need to obtain the color singlet wave functions for the four-quark states. The detailed coupling pathways can be found in our previous work [54].
χ d c 1 = 3 6 ( r g r ¯ g ¯ r g g ¯ r ¯ + g r g ¯ r ¯ g r r ¯ g ¯ + r b r ¯ b ¯ r b b ¯ r ¯ + b r b ¯ r ¯ b r r ¯ b ¯ + g b g ¯ b ¯ g b b ¯ g ¯ + b g b ¯ g ¯ b g g ¯ b ¯ ) . χ d c 2 = 6 12 ( 2 r r r ¯ r ¯ + 2 g g g ¯ g ¯ + 2 b b b ¯ b ¯ + r g r ¯ g ¯ + r g g ¯ r ¯ + g r g ¯ r ¯ + g r r ¯ g ¯ + r b r ¯ b ¯ + r b b ¯ r ¯ + b r b ¯ r ¯ + b r r ¯ b ¯ + g b g ¯ b ¯ + g b b ¯ g ¯ + b g b ¯ g ¯ + b g g ¯ b ¯ ) .
where, χ d c 1 and χ d c 2 is the color antitriplet-triplet ( 3 ¯ × 3 ) and sextet-antisextet ( 6 × 6 ¯ ) color configuration, respectively.
(2) Meson-meson structure
For the spin freedom, the wave functions are independent of the structure of the system and they are the same as those for the diquark–antidiquark structure, Equation (5).
For the flavor part, the wave function takes the following form:
χ m 0 f = ( c ¯ c ) ( c ¯ c ) .
The marks m 0 of χ represent the meson-meson structure and isospin I = 0 .
For the color freedom, the color singlet wave functions for the four-quark states are,
χ m c 1 = 1 3 ( r ¯ r + g ¯ g + b ¯ b ) ( r ¯ r + g ¯ g + b ¯ b ) , χ m c 2 = 2 12 ( 3 b ¯ r r ¯ b + 3 g ¯ r r ¯ g + 3 b ¯ g g ¯ b + 3 g ¯ b b ¯ g + 3 r ¯ g g ¯ r + 3 r ¯ b b ¯ r + 2 r ¯ r r ¯ r + 2 g ¯ g g ¯ g + 2 b ¯ b b ¯ b r ¯ r g ¯ g g ¯ g r ¯ r b ¯ b g ¯ g b ¯ b r ¯ r g ¯ g b ¯ b r ¯ r b ¯ b ) .
where χ m c 1 and χ m c 2 represents the color singlet-singlet ( 1 × 1 ) and color octet-octet ( 8 × 8 ) configuration, respectively.
The orbital wave functions for the four-quark states can be obtained by coupling the orbital wave function for each relative motion of the system:
Ψ L M L = [ Ψ l 1 ( r 12 ) Ψ l 2 ( r 34 ) ] l 12 Ψ L r ( r 1234 ) L M L ,
where Ψ l 1 and Ψ l 2 are the orbital wave functions of the two sub-clusters, with an angular momentum of l 1 and l 2 , respectively. Ψ L r ( r 1234 ) is the wave function of the relative motion between two sub-clusters with orbital angular momentum L r . The total orbital angular momentum of the four-quark state is L. “[ ]” denotes the Clebsh–Gordan coupling. As a preliminary calculation, here, we take an approximation that all angular momenta ( l 1 , l 2 , L r , L ) are set to zero. The Jacobi coordinates used are defined as follows:
r 12 = r 1 r 2 , r 34 = r 3 r 4 , r 1234 = m 1 r 1 + m 2 r 2 m 1 + m 2 m 3 r 3 + m 4 r 4 m 3 + m 4 .
For the diquark–antidiquark structure, the two quarks c c in one cluster are numbered as 1 , 2 , and the two antiquarks c ¯ c ¯ in the other cluster are numbered as 3 , 4 ; for the meson-meson structure, numbers 1 , 2 denote the antiquark and quark c ¯ c in one cluster, and numbers 3 , 4 denote the antiquark and quark c ¯ c in the other cluster. In the mixing calculations for the two structures, the marks of the quarks, antiquarks in the diquark–antidiquark structure are changed to be consistent with the numbering scheme in the meson-meson structure. In GEM, the orbital wave function is expanded by a series of Gaussians [55]:
Ψ l m ( r ) = n = 1 n max c n ψ n l m G ( r ) , ψ n l m G ( r ) = N n l r l e ν n r 2 Y l m ( r ^ ) ,
where N n l are normalization constants,
N n l = 2 l + 2 ( 2 ν n ) l + 3 2 π ( 2 l + 1 ) 1 2 .
c n are the variational parameters, which are determined dynamically. The Gaussian size parameters are chosen according to the following geometric progression:
ν n = 1 r n 2 , r n = r 1 a n 1 , a = r n max r 1 1 n max 1 .
This procedure enables optimized use of Gaussians, with as small a number of Gaussians used as possible.
So, we obtain the final channel wave function for the four-quark system c c c ¯ c ¯ in the diquark–antidiquark structure:
Ψ I J , i , j M I M J = A 1 [ Ψ L M L χ S M S σ i ] J M J χ d 0 f χ d c j , ( i = 1 6 ; j = 1 , 2 ; S = 0 , 1 , 2 ) ,
where A 1 is the antisymmetrization operator, for the c c c ¯ c ¯ system,
A 1 = 1 2 ( 1 P 12 P 34 + P 12 P 34 ) .
In the meson-meson structure, the final channel wave function for the c c c ¯ c ¯ system is written as
Ψ I J , i , j M I M J = A 2 [ Ψ L M L χ S M S σ i ] J M J χ m 0 f χ m c j , ( i = 1 6 ; j = 1 , 2 ; S = 0 , 1 , 2 ) ,
where A 2 is the antisymmetrization operator, for the c c ¯ c c ¯ system,
A 2 = 1 2 ( 1 P 13 P 24 + P 13 P 24 ) .
Finally, we obtain the eigenenergies for the c c c ¯ c ¯ system by solving a Schrödinger equation:
H Ψ I M I J M J = E I J Ψ I M I J M J .
Ψ I M I J M J is the wave function of the four-quark states, which is the linear combinations of the above channel wave functions, Equation (15) in the pure diquark–antidiquark structure or Equation (17) in the pure meson-meson structure, or both wave functions of Equations (15) and (17), respectively.

3. Results and Discussions

In this work, we estimated the masses of the tetraquark states for c c c ¯ c ¯ with the quantum numbers J P C = 0 + + , 1 + , 2 + + in the chiral quark model by adopting GEM. The pure meson-meson and the pure diquark–antidiquark structure, along with the dynamical mixing of these two structures are considered, respectively. In our calculations, all possible color, and spin configurations are included, and the approximation that all orbital angular momenta are set to 0 is used. In Table 3, we give the possible channels and corresponding wave functions.
The single-channel and channel-coupling calculations are performed in the present work. To determine whether or not any bound states exist, Table 4 gives the low-lying energies for c c c ¯ c ¯ tetraquarks with quantum numbers J P C = 0 + + , 1 + , 2 + + , respectively. From the table, we find that all the energies obtained are above the corresponding theoretical thresholds, so no bound states are formed in our calculations. In addition, the energies of the hidden-color states with color configurations 8 × 8 , 3 ¯ × 3 and 6 × 6 ¯ are higher than those with color 1 × 1 .
Although there is no bound state for the fully-heavy tetraquark system, resonance states with energies higher than the corresponding thresholds are possible because of the color structures of the system. In the present work, we employ the dedicated real scaling (stabilization) method to determine the genuine resonances. Because the calculation is carried out in a finite space, all the energies obtained are discrete. The bound state has a stable energy in the channel coupling calculation, while the resonance has a quasi-stable energy with increasing space after coupling to the scattering states. The real scaling method has been often used for analyzing electron–atom and electron–molecule scattering [46]. In the present approach, the real scaling method is realized by scaling the Gaussian size parameters r n in Equation (14) just for the meson-meson structure with the 1 × 1 color configuration, i.e., r n α r n , where α takes values between 0.9 and 1.7.
(a) J P C = 0 + + : There are six channels in this case. The stabilization plots for some channels are given. The first channel corresponds to the ( η c η c ) 0 configuration. Figure 1 shows the behavior of the energy spectrum of the first channel under the scaling of space. In this figure, most energies for the states decrease with increasing α , and they are scattering states. However, there are several horizontal lines, which correspond to the thresholds, η c ( 1 S ) η c ( 1 S ) , η c ( 1 S ) η c ( 2 S ) , η c ( 1 S ) η c ( 3 S ) and η c ( 2 S ) η c ( 2 S ) . The second channel is the hidden-color channel ( η c 8 η c 8 ) 0 . η c 8 is the color octet state of c c ¯ with spin 0. To check the α dependence of the energies of the hidden color channel, we also changed r n to α r n for the relative motion between two sub-clusters in the hidden color channel. Figure 2 gives a stabilization plot of the second channel, and all the energies are stable against the scaling of space. Figure 3 shows the stabilization plot for the third channel, ( J / ψ J / ψ ) 0 configuration. We observe similar behavior as that found for the first channel ( η c η c ) 0 and some thresholds appear. The remaining three channels ( J / ψ 8 J / ψ 8 ) 0 , ( D A ) 1 0 and ( D A ) 2 0 show similar behavior as observed for the second channel ( η c 8 η c 8 ) 0 , which is omitted here to save space. From these three figures, we can see that the energies of the hidden-color channels are independent with the scaling of space, but for the 1 × 1 color configuration.
To show the coupling effect, the following calculations were carried out. We first studied the coupling effects between hidden-color channels. In Table 5, the energies for the single hidden-color channel and hidden-color channels coupling are shown. To save space, only the first ten energies are listed. From the table, we can see that the coupling has important effects on the energy spectrum because all the energies of the hidden-color channels lie in the same energy range. For example, hidden-color meson-meson structure has the lowest energy of 6328 MeV, while the diquark–antidiquark structure has a lowest energy of 6344 MeV, and the coupling between structures reduces this energy to 6253 MeV. The coupling results for all hidden-channels are stable against the scaling of space, which is shown in Figure 4. The energy of the lowest state is 6253 MeV. There is another state close to it with an energy of 6316 MeV. At around 6500 MeV, there are two states, and, in the range of 6700–7000 MeV, 8 resonance states appear. Which resonance states can be observed among these states? A calculation for the coupling of the hidden-color channels to the scattering states is required.
Figure 5 shows the results obtained by considering the coupling between hidden-color channels ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 and the two scattering states ( η c η c ) 0 + ( J / ψ J / ψ ) 0 . One can find that all the resonance states from the hidden color channels ( ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 ) below 7100 MeV are missing. Only one resonance with an energy of 7138 MeV survives the coupling to the scattering states. Because there is no scattering state in the diquark–antidiquark structure, we also perform a channel coupling calculation involving two diquark–antidiquark structures, ( D A ) 1 0 + ( D A ) 2 0 , and the two scattering states, ( η c η c ) 0 + ( J / ψ J / ψ ) 0 . The results are displayed in Figure 6. Unlike the channel coupling calculations between the hidden-color channels ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 and the two scattering states ( η c η c ) 0 + ( J / ψ J / ψ ) 0 , several resonance states appear in the spectral region of 6700–7200 MeV. So, to identify the resonance states, full channel-coupling calculations are needed. From these calculations, we also find that the low-lying resonance states in the hidden-color channels are all missing.
To understand the reason behind the missing low-lying resonance states, we carried out the following calculations: Figure 7 illustrates the results obtained for the coupling between all hidden-color channels and the color-singlet channel ( η c η c ) 0 . The lowest state with energy of 6253 MeV disappears, which means that the state strongly couples to ( η c η c ) 0 . The resonance states around 6500 MeV are also missing in the coupling. Figure 8 gives the results for the coupling between all hidden-color channels and the color-singlet channel ( J / ψ J / ψ ) 0 . The second lowest state with energy of 6316 MeV disappears, which denotes that this state couples strongly to ( J / ψ J / ψ ) 0 . So, the resonance states below 6714 MeV are all missing after coupling to the scattering states.
To identify the genuine resonances, the full channel-coupling calculations are performed and the results are shown in Figure 9. The first resonance state R ( 6763 ) appears at an energy of 6763 MeV. There are two other resonances below 7000 MeV, with energies of 6849 MeV and 6884 MeV, which can be good candidates for the narrow structure around 6.9 GeV reported by the LHCb collaboration. The broad structure in the range 6.2∼6.8 GeV is due to the effect of the mixture of scattering channels, J / ψ J / ψ , η c ( 1 S ) η c ( 2 S ) and J / ψ ψ ( 2 S ) opening and the resonance R ( 6763 ) . As for the structure around 7.2 GeV, there are too many resonances here and our present calculations cannot give a clear picture. More scattering states with non-zero orbital angular momentum are required to be added.
(b) J P C = 1 + : There are five channels, four with a meson-meson structure and one with a diquark–antidiquark structure. The possible resonances, the results for all hidden-color channels coupling are given in Figure 10. After coupling to the scattering channels, the remaining resonances can be read from Figure 11. From Figure 10, one can see that the energy of the lowest state is 6346 MeV in the hidden-color channels, but this state disappears after coupling the open channel, η c ( 1 S ) J / ψ . However, the second lowest state survives the coupling (see Figure 11), so we obtain a resonance state with an energy of 6383 MeV. There are several resonance states around 6.8 GeV, R ( 6749 ) , R ( 6793 ) , R ( 6828 ) and R ( 6835 ) . There are also several resonances above 7.0 GeV, R ( 7100 ) , R ( 7112 ) , R ( 7173 ) , R ( 7173 ) , R ( 7185 ) , R ( 7247 ) , R ( 7249 ) and R ( 7264 ) . In the meson-meson structure, there are also states with J P C = 1 + + . The calculation shows that there are no resonances below 7.3 GeV, and the results are omitted here to save space.
(c) J P C = 2 + + : There are three channels, two with a meson-meson structure and one with a diquark–antidiquark structure. The results for the two hidden-color channels coupling are given in Figure 12. After coupling to the scattering channel ( J / ψ J / ψ ) 2 , the possible resonances can be read from Figure 13. There are three hidden-color states with energies of 6335 MeV, 6590 MeV and 6742 MeV, below 6800 MeV, respectively. However, they all disappear after coupling to the scattering channel ( J / ψ J / ψ ) 2 . The surviving states after the coupling are shown in Figure 13. These states are R ( 6855 ) , R ( 7085 ) , R ( 7139 ) , R ( 7204 ) and R ( 7270 ) .
Table 6 collects all possible resonance states with quantum numbers J P C = 0 + + , 1 + , 2 + + . The lowest observable resonance is the state with energy of 6383 MeV with J P C = 1 + , and, below 7.0 GeV, there are nine states, which are separated into three groups. The first group has one state with an energy of 6383 MeV; the second group has three states with energies in the range of 6.75∼6.8 GeV, and the final group has five states with energies around 6.85∼6.9 GeV. The broad structure in the region 6.2∼6.8 GeV observed by the LHCb collaboration can be explained by the mixed effects of the resonance states in the first and second groups, and the thresholds for double J / ψ and η c ( 1 S ) η c ( 2 S ) . The narrow structure around 6.9 GeV can be well explained by the resonances R ( 6849 ) and R ( 6884 ) with J P C = 0 + + and/or R ( 6855 ) with J P C = 2 + + . As for the resonances above 7.0 GeV, we obtain nine states below 7.3 GeV with double J / ψ as the final states ( J P C = 0 + + , 2 + + ). Here, we have mentioned that not all the states can be observed, because we have not taken into account angular excited scattering states, such as χ c J χ c J .
In order to identify the structures of these possible resonances, we calculated the distance between c and c ¯ quark, denoted as R c c ¯ , as well as the distance between c and c quark, denoted as R c c for the resonance states, respectively, which are shown in Table 7. From this table, we can see that for all the states with J P C = 0 + + and 2 + + , and some states with J P C = 1 + , R c c and R c c ¯ are all in the range of 0.6∼1.0 fm, which means that the states are in the diquark–antidiquark structure. For the other states with J P C = 1 + , R c c and R c c ¯ are larger than 1 fm, which means that the states are very likely to be molecular. The large R c c ¯ is due to the antisymmetrization and it gives the average distance between c and two c ¯ . The distance between c and c ¯ in one sub-cluster can be extracted from R c c and R c c ¯ , which is round 0.6 fm; this is consistent with the results obtained for charmonium.

4. Summary

In this work, we studied the mass spectra of the fully-charm c c c ¯ c ¯ system with quantum numbers J P C = 0 + + , 1 + , 2 + + in the chiral quark model with the help of GEM. The dynamical mixing of the meson-meson structure and the diquark–antidiquark structure, along with all possible color, spin configurations were taken into account. The predicted masses for the lowest-lying c c c ¯ c ¯ states are all above the corresponding two meson decay thresholds, leaving no space for bound states. By adopting the real scaling method, we identify the genuine resonance states for the J P C = 0 + + , 1 + , 2 + + c c c ¯ c ¯ system. The lowest observable resonance state has quantum numbers J P C = 1 + with energy 6.38 GeV, and there are several states with quantum numbers J P C = 0 + + , 2 + + with energies in the region of 6.85∼6.90 GeV, which can be used to explain the narrow structure around 6.9 GeV in the double J / ψ invariant mass spectrum reported by the LHCb collaboration. There are many resonance states above 7.0 GeV,, which are worth further study by including the angular excited scattering channels.
In general, the Q Q Q ¯ Q ¯ ( Q = b , c ) resonance states mainly decay into two Q Q ¯ meson final states by spontaneous dissociation. For the full-charm c c c ¯ c ¯ tetraquarks, they can decay via the spontaneous dissociation mechanism since they lie above the two-charmonium thresholds. The doubly hidden-charm tetraquarks can be clearly differentiated in experiments because of their much heavier energies than conventional charmonium mesons c c ¯ . However, it is more difficult to produce c c c ¯ c ¯ states because two heavy quark pairs need to be created in vacuum. However, the running of the LHC provides a chance for the observation of double heavy quarkonium. For example, the recent observations of the J / ψ J / ψ [10,11,12,13,14], J / ψ Y ( 1 S ) [56] and Y ( 1 S ) Y ( 1 S ) [15] events, shine a light for the production of c c c ¯ c ¯ and b b b ¯ b ¯ tetraquarks. With the accumulation of experimental data and the advancement of theoretical work, the nature of the fully-heavy tetraquark states will become more clear.

Funding

This research was funded by National Natural Science Foundation of China grant number 11847145.

Acknowledgments

We are grateful for insightful comments from Xiang Liu, Jialun Ping. Work supported by: National Natural Science Foundation of China, under Grant No. 11847145.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Choi, S.-K.; Olsen S., L.; Abe, K.; Abe, T.; Adachi, I.; Ahn, B.; Aihara, H.; Akai, K.; Akatsu, H.; Akemoto, M.; et al. Observation of a narrow charmoniumlike state in exclusive B± → K±π+π J/ψ decays. Phys. Rev. Lett. 2003, 91, 262001. [Google Scholar] [CrossRef] [Green Version]
  2. Aubert, B.; Barate1, R.; Boutigny, D.; Couderc, F.; Karyotakis, Y.; Lees, J.P.; Poireau, V.; Tisserand, V.; Zghiche, A.; Grauges, E.; et al. Observation of a broad structure in the π+π J/ψ mass spectrum around 4.26 GeV/c2. Phys. Rev. Lett. 2005, 95, 142001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yuan, C.Z.; Shen, C.P.; Wang, P.; McOnie, S.; Adachi, I.; Aihara, H.; Aulchenko, V.; Aushev, T.; Bahinipati, S.; Balagura, V.; et al. Measurement of the e+ eπ+π J/ψ cross section via initial-state radiation at Belle. Phys. Rev. Lett. 2007, 99, 182004. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Ablikim, M.; Achasov, M.N.; Ai, X.C.; Albayrak, O.; Ambrose, D.J.; An, F.F.; An, Q.; Bai, J.Z.; Ferroli, R.; Baldini; et al. Observation of a Charged Charmoniumlike Structure in e+ eπ+π J/ψ at s = 4.26 GeV. Phys. Rev. Lett. 2013, 110, 252001. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, Z.Q.; Shen, C.P.; Yuan, C.Z.; Adachi, I.; Aihara, H.; Asner, D.M.; Aulchenko, V.; Aushev, T.; Aziz, T.; Bakich, A.M.; et al. Study of e+ eπ+π J/ψ and Observation of a Charged Charmoniumlike State at Belle. Phys. Rev. Lett. 2013, 110, 252002. [Google Scholar] [CrossRef] [Green Version]
  6. Xiao, T.; Dobbs, S.; Seth, K.K.; Tomaradze, A. Observation of the Charged Hadron Zc (3900) at sqrt (s) = 4170 MeV. Phys. Lett. B 2013, 727, 366. [Google Scholar] [CrossRef] [Green Version]
  7. Ablikim, M.; Achasov, M.N.; Ai, X.C.; Albayrak, O.; Albrecht, M.; Ambrose, D.J.; Amoroso, A.; An, F.F.; An, Q.; Bai, J.Z.; et al. Observation of Zc (3900) 0 in e+ eπ0π0 J/ψ. Phys. Rev. Lett. 2015, 115, 112003. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Bondar, A.; Garmash, A.; Mizuk, R.; Santel, D.; Kinoshita, K.; Adachi, I.; Aihara, H.; Arinstein, K.; Asner, D.M.; Aushev, T.; et al. Observation of Two Charged Bottomoniumlike Resonances in Υ (5 S) Decays. Phys. Rev. Lett. 2012, 108, 122001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Aaij, R.; Adeva, B.; Adinolfi, M.; Affolder, A.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; et al. Observation of J/ψ p Resonances Consistent with Pentaquark States in Λ b 0 → J/ψ K- ρ Decays. Phys. Rev. Lett. 2015, 115, 072001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Aaij, R.; Adeva, B.; Adinolfi, M.; Adrover, C.; Affolder, A.; Ajaltouni, Z.; Albrecht, J.; Alessio, F.; Alexander, M.; Alkhazov, G.; et al. Observation of J/ψ-pair production in pp collisions at s= 7 TeV. Phys. Lett. B 2012, 707, 52. [Google Scholar] [CrossRef]
  11. Aaij, R.; Adeva, B.; Adinolfi, M.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; Alkhazov, G.; et al. Measurement of the J/ψ pair production cross-section in pp collisions at s = 13⎷s = 13 TeV. J. High Energy Phys. 2017, 10, 068. [Google Scholar] [CrossRef] [Green Version]
  12. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Bergauer, T.; Dragicevic, M.; Erö, J.; Fabjan, C.; Friedl, M.; Fruehwirth, R.; et al. Meaurement of prompt J/ψ pair production in pp colliion at ⎷s= 7 Te. J. High Energy Phys. 2014, 09, 094. [Google Scholar] [CrossRef] [Green Version]
  13. Aaboud, M.; Aad, G.; Abbott, B.; Abdallah, J.; Abdinov, O.B.; Abeloos, B.; Abidi, S.H.; AbouZeid, H.; Abraham, N.L.; Abramowicz, H.; et al. Search for dark matter at ⎷s = 13 TeV in final states containing an energetic photon and large missing transverse momentum with the ATLAS detector. Eur. Phys. J. C 2017, 77, 76. [Google Scholar] [CrossRef] [Green Version]
  14. Abe, K.; Abe, R.; Abe, T.; Adachi, I.; Ahn, B.S.; Aihara, H.; Akatsu, M.; Asano, Y.; Aso, T.; Aulchenko, V.; et al. Observation of Double c c Production in e+ e Annihilation at ⎷s ≈ 10.6 GeV. Phys. Rev. Lett. 2002, 89, 142001. [Google Scholar] [CrossRef] [Green Version]
  15. Khachatryan, V.; Sirunyan, A.M.; Tumasyan, A.; Adam, W.; Asilar, E.; Bergauer, T.; Brandstetter, J.; Brondolin, E.; Dragicevic, M.; Erö, J.; et al. Observation of Y (1S) pair production in proton-proton collisions at ⎷s= 8 TeV. J. High Energy Phys. 2017, 1705, 013. [Google Scholar] [CrossRef]
  16. Durgut, S. Search for Exotic Mesons at CMS. Available online: http://meetings.aps.org/Meeting/APR18/Session/U09.6 (accessed on 16 April 2018).
  17. Aaij, R.; Beteta, C.A. Observation of structure in the J/ψ -pair mass spectrum. Sci. Bull. 2020, 65, 1983–1993. [Google Scholar]
  18. Zhu, R. Fully-heavy tetraquark spectra and production at hadron colliders. Nucl. Phys. B 2021, 966, 115393. [Google Scholar]
  19. Zhang, J.R. 0+ fully-charmed tetraquark states. Phys. Rev. D 2021, 103, 014018. [Google Scholar] [CrossRef]
  20. Chen, H.X.; Chen, W.; Liu, X.; Zhu, S.L. Strong decays of fully-charm tetraquarks into di-charmonia. Sci. Bull. 2020, 65, 1994–2000. [Google Scholar] [CrossRef]
  21. Becchi, C.; Ferretti, J.; Giachino, A.; Maiani, L.; Santopinto, E. A study of ccc¯c¯ tetraquark decays in 4 muons and in D(*)D¯(*) at LHC. Phys. Lett. B 2020, 811, 135952. [Google Scholar] [CrossRef]
  22. Weng, X.Z.; Chen, X.L.; Deng, W.Z.; Zhu, S.L. Systematics of fully heavy tetraquarks. Phys. Rev. D 2021, 103, 034001. [Google Scholar] [CrossRef]
  23. Zhao, J.; Shi, S.; Zhuang, P. Fully-heavy tetraquarks in a strongly interacting medium. Phys. Rev. D 2020, 102, 114001. [Google Scholar] [CrossRef]
  24. Lü, Q.F.; Chen, D.Y.; Dong, Y.B. Masses of fully heavy tetraquarks QQQ¯ Q¯ in an extended relativized quark model. Eur. Phys. J. C 2020, 80, 871. [Google Scholar] [CrossRef]
  25. Bedolla, M.A.; Ferretti, J.; Roberts, C.D.; Santopinto, E. Spectrum of fully-heavy tetraquarks from a diquark+antidiquark perspective. Eur. Phys. J. C 2020, 80, 1004. [Google Scholar] [CrossRef]
  26. Iwasaki, Y. A possible model for new resonances: Exotics and hidden charm. Prog. Theor. Phys. 1975, 54, 492. [Google Scholar] [CrossRef] [Green Version]
  27. Heller, L.; Tjon, J.A. Bound states of heavy Q2 Q−2 systems. Phys. Rev. D 1985, 32, 755. [Google Scholar] [CrossRef]
  28. Lloyd, R.J.; Vary, J.P. All-charm tetraquarks. Phys. Rev. D 2004, 70, 014009. [Google Scholar] [CrossRef] [Green Version]
  29. Berezhnoy, A.V.; Luchinsky, A.V.; Novoselov, A.A. Heavy tetraquarks production at the LHC. Phys. Rev. D 2012, 86, 034004. [Google Scholar] [CrossRef] [Green Version]
  30. Chen, W.; Chen, H.X.; Liu, X.; Steele, T.G.; Zhu, S.L. Hunting for exotic doubly hidden-charm/bottom tetraquark states. Phys. Lett. B 2017, 773, 247. [Google Scholar] [CrossRef]
  31. Bai, Y.; Lu, S.; Osborne, J. Beauty-full tetraquarks. Phys. Lett. B 2019, 798, 134930. [Google Scholar] [CrossRef]
  32. Esposito, A.; Polosa, A.D. A b b b ¯ b ¯ di-bottomonium at the LHC? Eur. Phys. J. C 2018, 78, 782. [Google Scholar] [CrossRef]
  33. Debastiani, V.R.; Navarra, F.S. A non-relativistic model for the tetraquark. Chin. Phys. C 2019, 43, 013105. [Google Scholar] [CrossRef]
  34. Chao, K.T. The cccc (diquark-antidiquark) states in e+ e annihilation. Z. Phys. C 1981, 7, 317. [Google Scholar] [CrossRef]
  35. Ader, J.P.; Richard, J.M.; Taxil, P. Do narrow heavy multiquark states exist? Phys. Rev. D 1982, 25, 2370. [Google Scholar] [CrossRef] [Green Version]
  36. Barnea, N.; Vijande, J.; Valcarce, A. Four-quark spectroscopy within the hyperspherical formalism. Phys. Rev. D 2006, 73, 054004. [Google Scholar] [CrossRef] [Green Version]
  37. Wang, Z.G. Analysis of the QQQ¯ Q¯ tetraquark states with QCD sum rules. Eur. Phys. J. C 2017, 77, 432. [Google Scholar] [CrossRef]
  38. Karliner, M.; Nussinov, S.; Rosner, J.L. QQQ¯ Q¯ states: Masses, production, and decays. Phys. Rev. D 2017, 95, 034011. [Google Scholar] [CrossRef] [Green Version]
  39. Wu, J.; Liu, Y.R.; Chen, K.; Liu, X.; Zhu, S.L. Heavy-flavored tetraquark states with the QQQ¯ Q¯ configuration. Phys. Rev. D 2018, 97, 094015. [Google Scholar] [CrossRef] [Green Version]
  40. Liu, M.S.; Lü, Q.F.; Zhong, X.H.; Zhao, Q. All-heavy tetraquarks. Phys. Rev. D 2019, 100, 016006. [Google Scholar] [CrossRef] [Green Version]
  41. Gordillo, M.C.; Soto, F.D.; Segovia, J. Diffusion Monte Carlo calculations of fully-heavy multiquark bound states. Phys. Rev. D 2020, 102, 114007. [Google Scholar] [CrossRef]
  42. Yang, G.; Ping, J.L.; He, L.Y.; Wang, Q. A potential model prediction of fully-heavy tetraquarks QQQ¯ Q¯ (Q = c, b). arXiv 2006, arXiv:2006.13756v3. [Google Scholar]
  43. Yang, G.; Ping, J.L.; Segovia, J. Tetra- and Penta-Quark Structures in the Constituent Quark Model. Symmetry 2020, 12, 1869. [Google Scholar] [CrossRef]
  44. Chen, X. Analysis of hidden-bottom b b b ¯ b ¯ states. Eur. Phys. J. A 2019, 55, 106. [Google Scholar] [CrossRef]
  45. Jin, X.; Xue, Y.Y.; Huang, H.X.; Ping, J.L. Full-heavy tetraquarks in constituent quark models. arXiv 2006, arXiv:2006.13745. [Google Scholar]
  46. Simon, J. Resonance state lifetimes from stabilization graphs. J. Chem. Phys. 1981, 75, 2465. [Google Scholar] [CrossRef]
  47. Hiyama, E.; Kamimura, M.; Hosaka, A.; Toki, H.; Yahiro, M. Five-body calculation of resonance and scattering states of pentaquark system. Phys. Lett. B 2006, 633, 237. [Google Scholar] [CrossRef] [Green Version]
  48. Hiyama, E.; Hosaka, A.; Oka, M.; Richard, J. Quark model estimate of hidden-charm pentaquark resonances. Phys. Rev. C 2018, 98, 045208. [Google Scholar] [CrossRef] [Green Version]
  49. Myo, T.; Kikuchi, Y.; Masui, H.; Kato, K. Scattering phase shifts in the complex scaling method. Prog. Part. Nucl. Phys. 2014, 79, 1. [Google Scholar] [CrossRef] [Green Version]
  50. Kurtz, H.A.; Öhrn, Y. Calculation of 2p shape resonances in Be and Mg. Phys. Rev. A 1979, 19, 43. [Google Scholar] [CrossRef]
  51. Huang, H.X.; Ping, J.L. Investigating the hidden-charm and hidden-bottom pentaquark resonances in scattering process. Phys. Rev. D 2019, 99, 014010. [Google Scholar] [CrossRef] [Green Version]
  52. Chen, X.Y.; Ping, J.L.; Roberts, C.D.; Segovia, J. Light-meson masses in an unquenched quark model. Phys. Rev. D 2018, 97, 094016. [Google Scholar] [CrossRef] [Green Version]
  53. Vijande, J.; Fernández, F.; Valcarce, A. Constituent quark model study of the meson spectra. J. Phys. G. 2005, 31, 481. [Google Scholar] [CrossRef] [Green Version]
  54. Chen, X.Y.; Ping, J.L. Looking for a u d ¯ s ¯ b ¯ bound state in the chiral quark model. Phy. Rev. D 2018, 98, 054022. [Google Scholar] [CrossRef] [Green Version]
  55. Hiyama, E.; Kino, Y.; Kamimura, M. Gaussian expansion method for few-body systems. Prog. Part. Nucl. Phys. 2003, 51, 223. [Google Scholar] [CrossRef]
  56. Abazov, V.M.; Abbott, B.; Acharya, B.S.; Adams, M.; Adams, T.; Agnew, J.P.; Alexeev, G.D.; Alkhazov, G.; Alton, A.; Askew, A.; et al. Evidence for simultaneous production of J/ψ and Υ mesons. Phys. Rev. Lett. 2016, 116, 082002. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The stabilization plots for the energies of ( η c η c ) 0 state of J P C = 0 + + with respect to the scaling factor α . The horizontal axis shows the value of α and the vertical axis represents the energy of the states, which is the same for the following figures.
Figure 1. The stabilization plots for the energies of ( η c η c ) 0 state of J P C = 0 + + with respect to the scaling factor α . The horizontal axis shows the value of α and the vertical axis represents the energy of the states, which is the same for the following figures.
Universe 07 00155 g001
Figure 2. The stabilization plots for the energies of the ( η c 8 η c 8 ) 0 states for J P C = 0 + + with respect to the scaling factor α .
Figure 2. The stabilization plots for the energies of the ( η c 8 η c 8 ) 0 states for J P C = 0 + + with respect to the scaling factor α .
Universe 07 00155 g002
Figure 3. The stabilization plots for the energies of the ( J / ψ J / ψ ) 0 states for J P C = 0 + + with respect to the scaling factor α .
Figure 3. The stabilization plots for the energies of the ( J / ψ J / ψ ) 0 states for J P C = 0 + + with respect to the scaling factor α .
Universe 07 00155 g003
Figure 4. The stabilization plots for the energies of the states for J P C = 0 + + in the full hidden-color channels coupling ( ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 + ( D A ) 1 0 + ( D A ) 2 0 ) with respect to the scaling factor α .
Figure 4. The stabilization plots for the energies of the states for J P C = 0 + + in the full hidden-color channels coupling ( ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 + ( D A ) 1 0 + ( D A ) 2 0 ) with respect to the scaling factor α .
Universe 07 00155 g004
Figure 5. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the coupling between hidden-color channels ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 and the two scattering states ( η c η c ) 0 + ( J / ψ J / ψ ) 0 with respect to the scaling factor α . The red horizontal line represents the resonance states and the blue horizontal line represents the threshold, which is the same for the following figures.
Figure 5. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the coupling between hidden-color channels ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 and the two scattering states ( η c η c ) 0 + ( J / ψ J / ψ ) 0 with respect to the scaling factor α . The red horizontal line represents the resonance states and the blue horizontal line represents the threshold, which is the same for the following figures.
Universe 07 00155 g005
Figure 6. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the coupling between two diquark–antidiquark structures ( D A ) 1 0 + ( D A ) 2 0 and the two scattering states ( η c η c ) 0 + ( J / ψ J / ψ ) 0 with respect to the scaling factor α .
Figure 6. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the coupling between two diquark–antidiquark structures ( D A ) 1 0 + ( D A ) 2 0 and the two scattering states ( η c η c ) 0 + ( J / ψ J / ψ ) 0 with respect to the scaling factor α .
Universe 07 00155 g006
Figure 7. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the coupling between all hidden-color channels and the color-singlet channel ( η c η c ) 0 with respect to the scaling factor α .
Figure 7. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the coupling between all hidden-color channels and the color-singlet channel ( η c η c ) 0 with respect to the scaling factor α .
Universe 07 00155 g007
Figure 8. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the coupling between all hidden-color channels and the color-singlet channel ( J / ψ J / ψ ) 0 with respect to the scaling factor α .
Figure 8. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the coupling between all hidden-color channels and the color-singlet channel ( J / ψ J / ψ ) 0 with respect to the scaling factor α .
Universe 07 00155 g008
Figure 9. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the full-channel couplings with respect to the scaling factor α .
Figure 9. The stabilization plots for the energies of the states for J P C = 0 + + obtained by considering the full-channel couplings with respect to the scaling factor α .
Universe 07 00155 g009
Figure 10. The stabilization plots for the energies of the states for J P C = 1 + obtained by considering the full hidden-channels coupling with respect to the scaling factor α .
Figure 10. The stabilization plots for the energies of the states for J P C = 1 + obtained by considering the full hidden-channels coupling with respect to the scaling factor α .
Universe 07 00155 g010
Figure 11. The stabilization plots for the energies of the states for J P C = 1 + obtained by considering the full-channel couplings with respect to the scaling factor α . Upper for energies 6000–7000 MeV and down for energies 7000–7300 MeV.
Figure 11. The stabilization plots for the energies of the states for J P C = 1 + obtained by considering the full-channel couplings with respect to the scaling factor α . Upper for energies 6000–7000 MeV and down for energies 7000–7300 MeV.
Universe 07 00155 g011
Figure 12. The stabilization plots for the energies of the states for J P C = 2 + + obtained by considering the full hidden-channels coupling with respect to the scaling factor α .
Figure 12. The stabilization plots for the energies of the states for J P C = 2 + + obtained by considering the full hidden-channels coupling with respect to the scaling factor α .
Universe 07 00155 g012
Figure 13. The stabilization plots for the energies of the states for J P C = 2 + + obtained by considering the full-channel couplings with respect to the scaling factor α . Upper for energies 6100–7000 MeV and down for energies 7000–7300 MeV.
Figure 13. The stabilization plots for the energies of the states for J P C = 2 + + obtained by considering the full-channel couplings with respect to the scaling factor α . Upper for energies 6100–7000 MeV and down for energies 7000–7300 MeV.
Universe 07 00155 g013
Table 1. Model parameters, determined by fitting the meson spectrum from light to heavy.
Table 1. Model parameters, determined by fitting the meson spectrum from light to heavy.
Quark masses m u = m d 313
(MeV) m s 536
m c 1728
m b 5112
Confinement a c (MeV fm 2 )101
Δ (MeV)−78.3
OGE α 0 3.67
Λ 0 ( fm 1 ) 0.033
μ 0 (MeV)36.98
s 0 (MeV)28.17
Table 2. The masses of some heavy mesons (unit: MeV). M c a l and M e x p represent the theoretical and experimental masses, respectively.
Table 2. The masses of some heavy mesons (unit: MeV). M c a l and M e x p represent the theoretical and experimental masses, respectively.
Meson η c J / ψ h c χ c 0 χ c 1 χ c 2
M c a l 2986.33096.43417.33416.43416.43416.4
M e x p 2983.4 3096.9 3525.4 3414.8 3510.7 3556.2
Table 3. The possible channels and corresponding wave functions for c c c ¯ c ¯ tetraquarks with quantum numbers J P C = 0 + + , 1 + , 2 + + . The channels ( D A ) represent the diquark–antidiquark structures and the superscripts of the channels represent the total angular momentum J.
Table 3. The possible channels and corresponding wave functions for c c c ¯ c ¯ tetraquarks with quantum numbers J P C = 0 + + , 1 + , 2 + + . The channels ( D A ) represent the diquark–antidiquark structures and the superscripts of the channels represent the total angular momentum J.
J PC = 0 + + J PC = 1 + J PC = 2 + +
ChannelsWave FunctionsChannelsWave FunctionsChannelsWave Functions
( η c η c ) 0 χ 00 σ 1 χ m 0 f χ m c 1 ( η c J / ψ ) 1 χ 11 σ 3 χ m 0 f χ m c 1 ( J / ψ J / ψ ) 2 χ 22 σ 6 χ m 0 f χ m c 1
( η c 8 η c 8 ) 0 χ 00 σ 1 χ m 0 f χ m c 2 ( η c 8 J / ψ 8 ) 1 χ 11 σ 3 χ m 0 f χ m c 2 ( J / ψ 8 J / ψ 8 ) 2 χ 22 σ 6 χ m 0 f χ m c 2
( J / ψ J / ψ ) 0 χ 00 σ 2 χ m 0 f χ m c 1 ( J / ψ η c ) 1 χ 11 σ 4 χ m 0 f χ m c 1 ( D A ) 2 χ 22 σ 6 χ d 0 f χ d c 1
( J / ψ 8 J / ψ 8 ) 0 χ 00 σ 2 χ m 0 f χ m c 2 ( J / ψ 8 η c 8 ) 1 χ 11 σ 4 χ m 0 f χ m c 2
( D A ) 1 0 χ 00 σ 2 χ d 0 f χ d c 1 ( D A ) 1 χ 11 σ 5 χ d 0 f χ d c 1
( D A ) 2 0 χ 00 σ 1 χ d 0 f χ d c 2
Table 4. The low-lying energies of the c c c ¯ c ¯ states with J P C = 0 + + , 1 + , 2 + + in single channel and full-channel coupling calculations. E th theo and E th exp represent the theoretical and experimental thresholds, respectively (unit: MeV).
Table 4. The low-lying energies of the c c c ¯ c ¯ states with J P C = 0 + + , 1 + , 2 + + in single channel and full-channel coupling calculations. E th theo and E th exp represent the theoretical and experimental thresholds, respectively (unit: MeV).
ChannelE E th theo E th exp
J P C = 0 + +
( η c η c ) 0   5973.4    5972.6    5966.8   
( η c 8 η c 8 ) 0 6373.2
( J / ψ J / ψ ) 0 6193.76192.86193.8
( J / ψ 8 J / ψ 8 ) 0 6356.9
( D A ) 1 0 6360.2
( D A ) 2 0 6390.9
full-channel coupling5973.45972.65966.8
J P C = 1 +
( η c J / ψ ) 1   6083.6    6082.7    6080.3  
( η c 8 J / ψ 8 ) 1 6349.8
( J / ψ η c ) 1 6083.66082.76080.3
( J / ψ 8 η c 8 ) 1 6349.8
( D A ) 1 6397.6
full-channel coupling6083.66082.76080.3
J P C = 2 + +
( J / ψ J / ψ ) 2   6193.7    6192.8    6193.8  
( J / ψ 8 J / ψ 8 ) 2 6365.3
( D A ) 2 6410.4
full-channel coupling6193.76192.86193.8
Table 5. The energy spectra for the c c c ¯ c ¯ states with J P C = 0 + + in the single hidden-color channel and hidden-color channels coupling (unit: MeV).
Table 5. The energy spectra for the c c c ¯ c ¯ states with J P C = 0 + + in the single hidden-color channel and hidden-color channels coupling (unit: MeV).
Channel ( η c 8 η c 8 ) 0 ( J / ψ 8 J / ψ 8 ) 0 ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 ( DA ) 1 0 ( DA ) 2 0 ( DA ) 1 0 + ( DA ) 2 0 ( η c 8 η c 8 ) 0 + ( J / ψ 8 J / ψ 8 ) 0 + ( DA ) 1 0 + ( DA ) 2 0
Energy6371635563286388635863446253
6708671764026756669064006316
6973695166576841676366836528
7036706567656848697467596577
7140714269447149702267636714
7301725869787196706368406763
7357739369877200710268496805
7448745371107262726869726849
7507749071307266732670216877
7613757771547281732970636886
Table 6. The resonance states with quantum numbers J P C = 0 + + , 1 + , 2 + + for the c c c ¯ c ¯ system (unit: MeV).
Table 6. The resonance states with quantum numbers J P C = 0 + + , 1 + , 2 + + for the c c c ¯ c ¯ system (unit: MeV).
State 0 + + 1 + 2 + +
Energy    6763        6383        6855    
684967497085
688467937139
706368287204
709868357270
71387100
71567112
72007173
7185
7247
7249
7264
Table 7. The distances between c and c ( c ¯ ) quark for the possible resonance states of the c c c ¯ c ¯ system.
Table 7. The distances between c and c ( c ¯ ) quark for the possible resonance states of the c c c ¯ c ¯ system.
StateResonance (MeV) R c c ¯ (fm) R cc (fm)
0 + +   6763    0.68    0.82  
  6849    0.66    0.70  
  6884    0.80    0.85  
  7063    0.80    1.01  
  7098    0.84    1.03  
  7138    0.96    0.98  
  7156    0.88    0.81  
  7200    0.83    0.73  
1 +   6383    0.85    0.88  
67491.201.39
67930.520.48
68280.750.97
68351.351.31
71001.752.42
71121.001.26
71731.131.51
71851.551.66
72471.591.54
72491.601.50
72641.581.49
2 + + 68550.660.71
70850.940.74
71390.880.80
72040.840.72
72700.760.85
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, X. The Genuine Resonance of Full-Charm Tetraquarks. Universe 2021, 7, 155. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7050155

AMA Style

Chen X. The Genuine Resonance of Full-Charm Tetraquarks. Universe. 2021; 7(5):155. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7050155

Chicago/Turabian Style

Chen, Xiaoyun. 2021. "The Genuine Resonance of Full-Charm Tetraquarks" Universe 7, no. 5: 155. https://0-doi-org.brum.beds.ac.uk/10.3390/universe7050155

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop